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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » matlab求解非线性方程组
51/1返回列表
查看: 4402  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

sunjiahuan

铜虫 (小有名气)

[求助] matlab求解非线性方程组 已有1人参与

I1=(1.3403/0.323)*(354.3091/443.6719-lambertw((1.0885e-004/443.6719)*exp((330.7*x1+354.3091)/443.6719)))-x1/(0.323+330.7)+1.019*10^(-6)*(exp(x1/1.3403)-1);
I2=(1.3403/0.323)*(349.7334/437.9477-lambertw((1.0744e-004/437.9477)*exp((326.4292*x2+349.7334)/437.9477)))-x2/(0.323+326.4292)+1.019*10^(-6)*(exp(x2/1.3403)-1);
I3=(1.3403/0.323)*(337.973/423.2354-lambertw((1.0383e-004/423.2354)*exp((315.4524*x3+337.973)/423.2354)))-x3/(0.323+315.4524)+1.019*10^(-6)*(exp(x3/1.3403)-1);
I4=(1.3403/0.323)*(324.4644/406.3363-lambertw((9.9677e-005/406.3363)*exp((302.844*x4+324.4644)/406.3363)))-x4/(0.323+302.844)+1.019*10^(-6)*(exp(x4/1.3403)-1);
I5=1.019*10^(-6)*(exp(x5/1.3403)-1)；
上面是5个方程，I1=f(x1)，I2=f(x2)，I3=f(x3)，I4=f(x4),I5=f(x5); 关于方程里的lambertw函数，可以在附件（太阳电池I-V方程显示求解原理研究及应用中的第二章中找到）。

而要求解的方程组是由它们5个组成的,如下：
x1+x2+x3+x4+x5-30=0
I1-I2=0;
I1-I3=0;
I1-I4=0;
I1-I5=0;
上面的方程组已经写好，在fun.m里。

而在dfun.m 中是方程组里的每一个方程分别对x1,x2 x3 x4 x5 求一阶导数。其实就是求方程组的雅克比矩阵；

在newton.m里是求解的编程，是我在这个网上（http://hi.baidu.com/aillieo/blog ... c9a59647106493.html）找到的，只是把里面的方程及相关部分改了一下。

我在matlab里输入newton([0.01 0.01 0.01 0.01 0.01],0.00001,200)
最后得到的结果是
In newton at 6
Warning: Matrix is singular, close to singular or badly scaled.
         Results may be inaccurate. RCOND = NaN.
> In newton at 6

ans =

   NaN   NaN   NaN   NaN   NaN。
这种方法是好多学者用过的，得到的结果都挺好。所以方法应该没问题。因为我的matlab只学了个皮毛，不知道这里面哪出错了，还请各位高手给指点一下。
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : newton.m
    • 2012-02-22 16:56:51, 775 bytes
  • 附件 2 : dfun.m
    • 2012-02-22 16:56:51, 120 bytes
  • 附件 3 : fun.m
    • 2012-02-22 16:56:52, 969 bytes
  • 附件 4 : 太阳电池I_V方程显式求解原理研究及应用.nh
    • 2012-02-22 16:57:53, 4.29 M

» 收录本帖的淘帖专辑推荐

材料+新能源+催化+外语

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-02-22 17:39:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunjiahuan

铜虫 (小有名气)
送鲜花一朵
引用回帖:
楼: Originally posted by dbb627 at 2012-02-22 19:55:48:
错误提示信息是矩阵奇异，
可能初值不太合适
我用自带fsolve解了下
结果不是很好
[/code]
function aa
[X,FVAL,EXITFLAG,OUTPUT] =fsolve(@fun,[18.7783   0.01   18.1315   18.0062  -2],optimset('MaxFun ...

您好，非常感谢你的帮助，您能不能把你的matlab代码发到我的邮箱里么sunjiahuan0370@sina.com？我把你贴出来的代码移到我的matlab里居然没成功。谢谢了。
赞一下 回复此楼
3楼2012-02-22 20:29:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
ben_ladeng(金币+2): 谢谢应助 2012-02-22 21:59:15
ben_ladeng(专家考核): 2012-02-22 21:59:24
caemechanics(金币+5): 感谢应助 2012-02-23 16:18:30
错误提示信息是矩阵奇异，
可能初值不太合适
我用自带fsolve解了下
结果不是很好
[/code]
function aa
[X,FVAL,EXITFLAG,OUTPUT] =fsolve(@fun,[18.7783   0.01   18.1315   18.0062  -2],optimset('MaxFunEvals',1000))
function f=fun(x)
x1=x(1);
x2=x(2);
x3=x(3);
x4=x(4);
x5=x(5);
f1=x1+x2+x3+x4+x5-30;
f2=4.1495*(0.7986-lambertw(2.4534e-007*exp((330.7*x1+354.3091)/443.6719)))-x1/331.0230+1.019e-6*(exp(x1/1.3403)-1)-4.1495*(0.7986-lambertw(2.4533e-007*exp((326.4292*x2+349.7334)/437.9477)))+x2/(0.323+326.4292)-1.019*10^(-6)*(exp(x2/1.3403)-1);
f3=4.1495*(0.7986-lambertw(2.4534e-007*exp((330.7*x1+354.3091)/443.6719)))-x1/331.0230+1.019e-6*(exp(x1/1.3403)-1)-4.1495*(0.7985-lambertw(2.4532e-007*exp((315.4524*x3+337.973)/423.2354)))+x3/(0.323+315.4524)-1.019*10^(-6)*(exp(x3/1.3403)-1);
f4=4.1495*(0.7986-lambertw(2.4534e-007*exp((330.7*x1+354.3091)/443.6719)))-x1/331.0230+1.019e-6*(exp(x1/1.3403)-1)-4.1495*(0.7985-lambertw( 2.4531e-007*exp((302.844*x4+324.4644)/406.3363)))+x4/(0.323+302.844)-1.019*10^(-6)*(exp(x4/1.3403)-1);
f5=4.1495*(0.7986-lambertw(2.4534e-007*exp((330.7*x1+354.3091)/443.6719)))-x1/331.0230+1.019e-6*(exp(x1/1.3403)-1)-1.019*10^(-6)*(exp(x5/1.3403)-1);
f=[f1 f2 f3 f4 f5];
[/code]
Solver stopped prematurely.

fsolve stopped because it exceeded the function evaluation limit,
options.MaxFunEvals = 1000 (the selected value).


X =

   18.7895   19.0894   19.0963   19.1094  -46.0868


FVAL =

   -0.0022   -0.0180   -0.0172   -0.0175    2.7252


EXITFLAG =

     0


OUTPUT =

       iterations: 180
        funcCount: 1001
        algorithm: 'trust-region dogleg'
    firstorderopt: 0.0022
          message: [1x143 char]
赞一下(2人) 回复此楼

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

sunjiahuan
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
2楼2012-02-22 19:55:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunjiahuan

铜虫 (小有名气)
sunjiahuan: 回帖置顶 2012-02-22 20:46:39
不好意思，刚看到I1，I2,I3,I4方程中最后一项的里的x1，x2,x3,x4, 应该为-x1,-x2,-x3,-x4,而x5不变。所以我附件里fun.m的方程组里的也要相应改变，因无法编辑，在此说明一下。比如I1中的最后一项1.019*10^(-6)*(exp(x1/1.3403)-1);应改为1.019*10^(-6)*(exp(-x1/1.3403)-1);
赞一下 回复此楼
4楼2012-02-22 20:46:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunjiahuan

铜虫 (小有名气)
引用回帖:
楼: Originally posted by dbb627 at 2012-02-22 19:55:48:
错误提示信息是矩阵奇异，
可能初值不太合适
我用自带fsolve解了下
结果不是很好
[/code]
function aa
[X,FVAL,EXITFLAG,OUTPUT] =fsolve(@fun,[18.7783   0.01   18.1315   18.0062  -2],optimset('MaxFun ...

不好意思，刚看到I1，I2,I3,I4方程中最后一项的里的x1，x2,x3,x4, 应该为-x1,-x2,-x3,-x4,而x5不变。所以我附件里fun.m的方程组里的也要相应改变，因无法编辑，在此说明一下。比如I1中的最后一项1.019*10^(-6)*(exp(x1/1.3403)-1);应改为1.019*10^(-6)*(exp(-x1/1.3403)-1);非常感谢您的帮助。
赞一下 回复此楼
5楼2012-02-22 20:53:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297，工科调剂? +11 河南农业大学-能 2026-04-14 11/550 2026-04-19 20:07 by Equinoxhua
[考博] 申博自荐 +4 Linxia林夏 2026-04-13 4/200 2026-04-19 19:55 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 289 分105500药学专硕求调剂(找B区学校) +5 白云123456789 2026-04-13 5/250 2026-04-19 18:12 by Equinoxhua
[论文投稿] 有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业 20+4 之护着 2026-04-16 6/300 2026-04-19 13:00 by Aaron_zyn
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 320求调剂 +5 深郊akm 2026-04-17 5/250 2026-04-18 19:52 by 王珺璞
[考研] 接受任何调剂 +6 也就是栗子 2026-04-17 7/350 2026-04-18 17:20 by 涵竹刘
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[考研] 一志愿华中农业071010，320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 一志愿中科大材料与化工，353分还有调剂学校吗 +10 否极泰来2026 2026-04-15 12/600 2026-04-17 17:54 by mapenggao
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 一志愿沪9，生物学326求调剂 +9 刘墨墨 2026-04-15 9/450 2026-04-16 17:14 by 崔崔崔cccc
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
[考研] 302求调剂 +10 易！? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
信息提示
关闭
请填处理意见
关闭