| 查看: 2144 | 回复: 12 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
Matlab：Newton-Raphson 解线非线性性方程组的问题
|
|||
|
大家好，能否帮我解决一个Matlab用Newton-Raphson 解线非线性性方程组的问题，可以肯定这五个方程组肯定是有解的。程序代码在最后。 但是运行之后出现了问题，提示消息是：Warning: Matrix is singular, close to singular or badly scaled. Results may be inaccurate. RCOND = NaN. 因为刚学matlab，不是很熟悉这些东西和算法，能否帮我看一下问题出在哪？谢谢大家。 maxiter=1000; epsilon=0.0000000001; nsys=5; Voc=33.4; Isc=8.12; Vm=26.2; Im=7.63; Ns=54; Np=1; %equations F1 = @(x) Np*x(1)-Np*x(2)*(exp(Voc/Ns/x(5))-1)-Np*Voc/Ns/x(4); F2 = @(x) Np*x(1)-Np*x(2)*(exp(Isc*x(3)/Np/x(5))-1)-Isc*x(3)/x(4)-Isc; F3 = @(x) Np*x(1)-Np*x(2)*(exp((Vm/Ns+Im*x(3)/Np)/x(5))-1)-Np/x(4)*(Vm/Ns+Im*x(3)/Np)-Im; F4 = @(x) (Np*x(2)*exp((Vm+Im*Ns*x(3)/Np)/Ns/x(5))/Ns/x(5)+Np/Ns/x(4))/(1+x(3)*x(2)*exp((Vm+Im*Ns*x(3)/Np)/Ns/x(5))/x(5)+x(3)/x(4))-Im/Vm; F5 = @(x) (Np*x(2)*exp(Isc*Ns*x(3)/Np/Ns/x(5))/Ns/x(5)+Np/Ns/x(4))/(1+x(3)*x(2)*exp(Isc*Ns*x(3)/Np/Ns/x(5))/x(5)+x(3)/x(4))-1/x(4); F={F1,F2,F3,F4,F5}; %initial values x(1)=8.12; x(2)=2e-9; x(3)=6e-7; x(4)=8700; x(5)=0.0275; nloop=0; nconverg=0; while nloop dFidxj=zeros(nsys,nsys); for i=1:1:nsys for j=1:1:nsys fx=F{i}; xj=x; xj(j)=x(j)+epsilon/2; y1=fx(xj); xj(j)=x(j)-epsilon/2; y0=fx(xj); dFidxj(i,j)=(y1-y0)/epsilon; end end %% for i=1:1:nsys j=F{i}; D(i,1)=j(x); end %% A=inv(dFidxj); error=A*D; dx=error'; %% x=x-dx; %% nloop=nloop+1; nconverg=nconverg+1; %% maxdx=max(abs(error)); if maxdx end end %% %results output x nconverg |
回复此楼» 猜你喜欢
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
今天务委会开完了，明天出结果吗
已经有13人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有8人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有22人回复
-
什么时候开奖？
已经有13人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助 MATLAB解方程组-fslove
已经有7人回复
- 用不动点迭代求非线性方程组，求得的不是想要的解，怎么办
已经有12人回复
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- matlab求解非线性方程组
已经有16人回复
- matlab解微分方程组
已经有15人回复
- 求助有关Matlab有约束非线性最优化问题
已经有9人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab非线性参数拟合问题
已经有7人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 非线性方程组的迭代法（数值计算高手请进）
已经有7人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★
mataohk: 金币+1, ★★★很有帮助 2012-10-16 18:43:58
jjdg: 金币+1, 辛苦了 2012-10-17 00:53:24
mataohk: 金币+1, ★★★很有帮助 2012-10-16 18:43:58
jjdg: 金币+1, 辛苦了 2012-10-17 00:53:24
|
文献答案是针对你一楼还是你修改后三楼的公式？ 如果对一楼： 1stOpt验证（2楼结果）：公式左右相减 -3.375077995E-14 -5.329070518E-15 1.776356839E-15 6.77236045E-15 8.515250924E-9 文献验证：公式左右相减 -0.001784587439 -2.420748046E-10 -0.0001076813772 5.662268229E-5 -0.000112805466 如果对三楼： 1stOpt验证（4楼结果）：公式左右相减 8.679745341E-5 0.0378520978 -0.03793594713 -0.005267011529 0.0180096136 文献验证：公式左右相减 -0.001784587439 -2.420748046E-10 -0.0001076813772 5.662268229E-5 1.72645739 无论那种情况，文献值都要比1stOpt得到的结果差。 |
9楼2012-10-16 17:14:15
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2012-10-16 08:48:36
|
1. 谢谢你的回答和计算。但是我用你的程序到进去之后算出来的结果和你算出来的结果差别很大，我用的是1stOpt15PRO。是不是版本不同的原因？麻烦你告诉我问题在哪？ 2. 对于你的结果，x4的结算结果差别很大，数量级应该在10^3，不会到10^34. 其他结果还比较准确。 3. 现在我把第五个方程修改一下，你再帮我算算这一次的结果，我想看一下哪一种方法更准确。谢谢你。修改后的方程组为： Function np*x1-np*x2*(exp(voc/ns/x5)-1)-np*voc/ns/x4; np*x1-np*x2*(exp(isc*x3/np/x5)-1)-isc*x3/x4-isc; np*x1-np*x2*(exp((vm/ns+im*x3/np)/x5)-1)-np/x4*(vm/ns+im*x3/np)-im; (np*x2*exp((vm+im*ns*x3/np)/ns/x5)/ns/x5+np/ns/x4)/(1+x3*x2*exp((vm+im*ns*x3/np)/ns/x5)/x5+x3/x4)-im/vm; (np*x2*exp(Voc/ns/x5)/ns/x5+np/ns/x4)/(1+x3*x2*exp(Voc/ns/x5)/x5+x3/x4)-1/x4; |
3楼2012-10-16 09:57:46
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2012-10-16 11:58:18
28