应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » fsolve求解非线性方程组
61/1返回列表
查看: 1691  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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]
回复此楼

» 猜你喜欢
一个矛盾无知的人
1楼 2017-12-01 17:17:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

儒雅圣士

新虫 (正式写手)
得到的解是【-0.076195037566324  -0.000000573386969  -2.251437071748727   2.440695364463143   0.046798441863606】  请查看方程并确认是否写正确了，估计是你写错了方程
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

bobo1314
一起努力
2楼2017-12-02 18:28:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bobo1314

银虫 (小有名气)
方程已经过了无数次的检查，而且将前人迭代的正确解反带入原档程，几乎每个式子都在0附近

发自小木虫Android客户端
赞一下 回复此楼
一个矛盾无知的人
3楼2017-12-02 21:02:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bobo1314

银虫 (小有名气)
送红花一朵
引用回帖:
2楼: Originally posted by 儒雅圣士 at 2017-12-02 18:28:52
得到的解是【-0.076195037566324  -0.000000573386969  -2.251437071748727   2.440695364463143   0.046798441863606】  请查看方程并确认是否写正确了，估计是你写错了方程

方程已经过了无数次的检查，而且将前人迭代的正确解反带入原档程，几乎每个式子都在0附近

发自小木虫Android客户端
赞一下 回复此楼
一个矛盾无知的人
4楼2017-12-02 21:03:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
6个方程5个未知数，本身就不存在精确解吧。
而且，所谓的“这个应该有前人的正确解，[5.449 1.2e-9 0.092 2.725 0.028]”，带回原方程验证，误差是非常之大，下面随便一组近似解都要比这个“正解”好的多：

iph: -2.44512940626882E-8
i0: 234876089.693707
rs: 2.65964117382912E15
rp: 2.01138445074181E24
vt: -6.24686104906122E23
赞一下(1人) 回复此楼
5楼2017-12-03 13:16:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bobo1314

银虫 (小有名气)
引用回帖:
5楼: Originally posted by dingd at 2017-12-03 13:16:50
6个方程5个未知数，本身就不存在精确解吧。
而且，所谓的“这个应该有前人的正确解，”，带回原方程验证，误差是非常之大，下面随便一组近似解都要比这个“正解”好的多：

iph: -2.44512940626882E-8
i0: 2348 ...

非常感谢回复，主要是求其优化解，用的fsolve的l-m算法

发自小木虫Android客户端
赞一下 回复此楼
一个矛盾无知的人
6楼2017-12-03 16:20:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 bobo1314 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +45 医学老男孩 2026-08-13 99/4950 2026-08-18 12:56 by dotnetchen9056
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 KZyU1THERRVo 2026-08-18 3/150 2026-08-18 12:52 by vHWVV9D5Deb5
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 5/250 2026-08-18 11:25 by 無關想念
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
关闭
请填处理意见
关闭