| 查看: 1717 | 回复: 5 | ||
bobo1314银虫 (小有名气)
|
[求助]
fsolve求解非线性方程组 已有1人参与
|
|
我的程序，总有这样那样的问题，最终还是会出现错误 M文件： function F=myfunfsolve(x) Iph=x(1); I0=x(2); Rs=x(3); Rp=x(4); Vt=x(5); Voc=44.6/72; Vm=35.4/72; Isc=5.43; Im=4.95; Ns=1; Np=1; q=1.602E-19; K=1.381E-23; Alpha=0.0008; beita=-0.145; T0=298.15; T01=323.15; Eg0=1.16-7.02e-4*T0^2/(T0-1108); Eg01=1.16-7.02e-4*T01^2/(T01-1108); I01bi=T01^3/T0^3*exp(q*Eg0/(Vt)-q*T0*Eg01/(Vt)); F=[Np*Iph-Np*I0*(exp(Voc/Ns/Vt)-1)-Np/Ns*Voc/Rp; Np*Iph-Np*I0*(exp(Isc*Rs/Np/Vt)-1)-Isc*Rs/Rp-Isc; Np*Iph-Np*I0*(exp((Vm/Ns+Im/Np*Rs)/Vt)-1)-Np*(Vm/Ns+Im/Np*Rs)/Rp-Im; (Np/Ns/Vt*I0*exp((Vm+Im*Ns/Np*Rs)/Ns/Vt)+1/(Ns/Np*Rp))/(1+Rs/Vt*I0*exp((Vm+Im*Ns/Np*Rs)/Ns/Vt)+Rs/Rp)-Im/Vm; (-Np/Ns/Vt*I0*exp(Isc*Ns*Rs/Np/Ns/Vt)-1/(Ns/Np*Rp))/(1+Rs/Vt*I0*exp(Isc*Ns*Rs/Np/Ns/Vt)+Rs/Rp)+1/Rp; Np*Iph*(1+(T01-T0)*Alpha/Isc)-Np*I0*I01bi*(exp(Voc*(1+(T01-T0)*beita/Voc)/Ns/Vt*T0/T01)-1)-Np/Ns*Voc*(1+(T01-T0)*beita/Voc)/Rp]; 命令行输入： clear; x0 = [5.4 2e-9 0.02 2 0.05]; % Make a starting guess at the solution options=optimset('Display','iter','TolFun',1e-3,'TolX',1e-100,'Algorithm','levenberg-marquardt', 'MaxIter',1e4,'MaxFunEvals',1e5,'Diagnostics','on'); %option to display output [x,fval] = fsolve(@myfunfsolve5,x0,options); % Call solver 出现下列提示： No solution found. fsolve stopped because the problem appears regular as measured by the gradient, but the vector of function values is not near zero as measured by the selected value of the function tolerance. 我就不知道该如何处理了，如何才能得到解，顺便说一声，这个应该有前人的正确解，[5.449 1.2e-9 0.092 2.725 0.028] |
回复此楼» 猜你喜欢
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有284人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
bobo13142楼2017-12-02 18:28:52
bobo1314
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 493.6
- 帖子: 73
- 在线: 46.7小时
- 虫号: 3100332
- 注册: 2014-03-30
- 性别: GG
- 专业: 可再生与替代能源利用中的
3楼2017-12-02 21:02:48
bobo1314
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 493.6
- 帖子: 73
- 在线: 46.7小时
- 虫号: 3100332
- 注册: 2014-03-30
- 性别: GG
- 专业: 可再生与替代能源利用中的
4楼2017-12-02 21:03:16
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
5楼2017-12-03 13:16:50
bobo1314
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 493.6
- 帖子: 73
- 在线: 46.7小时
- 虫号: 3100332
- 注册: 2014-03-30
- 性别: GG
- 专业: 可再生与替代能源利用中的
6楼2017-12-03 16:20:51