求助:fsolve求解非线性方程组的问题
需要求解的方程组如下:
function y=fc(x)
pr=15e-6
ps=17e-6
ec=100e3
ef=200
tf=150e-7
ti=2e-7
q=1.602e-19
k=1.380e-23
t=300
eo=8.854e-14
esi=11.9
nd=1e16
ni=1.5e10
ei=3.9
dat=ec/log((1+pr/ps)/(1-pr/ps))
bat=q/(k*t)
ld=sqrt(esi*eo/(q*nd*bat))
y(1)=-sqrt(2)*esi*eo/(bat*ld)*sqrt((ni^2)/(nd^2)*(exp(-bat*x(1))+bat*x(1)-1)+exp(bat*x(1))-bat*x(1)-1)+ps*tanh((x(2)-ec)/(2*dat))
y(2)=x(1)+sqrt(2)*esi*eo/(bat*ld)*sqrt((ni^2)/(nd^2)*(exp(-bat*x(1))+bat*x(1)-1)+exp(bat*x(1))-bat*x(1)-1)*ti/(ei*eo)+x(2)*tf-5
y=[y(1) y(2)]
其中x(1)的范围为1左右,x(2)的范围在3e5左右。
我解的时候是这样的:
>>x0=[1 3e5];
>>fsolve('fc',x0)
运行以后,最后显示的是:
Optimization terminated: first-order optimality is less than options.TolFun.
ans =
1.0e+005 *
0.0000 1.0003
结果明显是不对的。方程是没有问题的。大家能够帮我解决一下这个问题吗??我很着急啊~~