| 查看: 2138 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
人气不行了
已经有10人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助 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
20