| 查看: 1568 | 回复: 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] |
回复此楼» 猜你喜欢
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有95人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有0人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有0人回复
bobo1314
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 493.6
- 帖子: 73
- 在线: 46.7小时
- 虫号: 3100332
- 注册: 2014-03-30
- 性别: GG
- 专业: 可再生与替代能源利用中的
3楼2017-12-02 21:02:48
bobo13142楼2017-12-02 18:28:52
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