解决一系列非线性方程遇到麻烦 . 所以我试图用scipy.optimize.fsolve在python中完成它 . 我的等式如下
$$
-x^{(\alpha+1)}+a*\rho^2*x^\alpha+(1-\alpha)a(1-\rho^2)*(xm/(a(1-\rho)^2))^alpha
$$
其中$ \ alpha $，$ \ rho $和$ xm $是参数 . 我想从lognormal中得到一些$ a $来解决$ x $ . 为此，我编写python代码如下 .
from scipy.optimize import *
from numpy.random import *
import matplotlib.pyplot as plt
draws=[np.random.lognormal(1,0.5) for i in range(n)]
x_1=[fsolve(lambda x:(-x**(alpha+1)+a*rho**2*x**alpha+(1-alpha)*a*(1-rho**2)*(xm/(a*(1-rho)**2))**alpha),0) for a in draws]
没运气！ . 有人可以帮帮我吗？