| 查看: 1925 | 回复: 10 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
matlab 求解非线性方程组未果
|
||||
|
非软件专业，只想求一个非线性方程组的解，就显示错误，求指导~..... 程序如下 function m=slag(n) x1=n(1); x2=n(2); x3=n(3); x4=n(4); x5=n(5); x6=n(6); x7=n(7); x8=n(8); x9=n(9); x10=n(10); x11=n(11); x12=n(12); x13=n(13); x14=n(14); x15=n(15); x16=n(16); x17=n(17); x18=n(18); x19=n(19); m(1)=exp((13836+25.44*1863)/8.314*1863)*x1^3*x2-x5; m(2)=exp((35288+204.59*1863)/8.314*1863)*x1^12*x2^7-x6; m(3)=exp((18946+25.44*1863)/8.314*1863)*x1*x2-x7; m(4)=exp((18566+23.68*1863)/8.314*1863)*x1*x2^2-x8; m(5)=x9*0; m(6)=exp((107800+13.44*1863)/8.314*1863)*x1^3*x3^2-x10; m(7)=exp((122731+8.85*1863)/8.314*1863)*x1^2*x3-x11; m(8)=exp((302587+66.24*1863)/8.314*1863)*x1^3*x3^2-x12; m(9)=exp((86231+0.995*1863)/8.314*1863)*x1*x3-x13; m(10)=exp((19512+28.19*1863)/8.314*1863)*x2^3*x3^2-x14; m(11)=exp((354932-417.63*1863)/8.314*1863)*x1^2*x2*x3-x15; m(12)=exp((1252543-277.06*1863)/8.314*1863)*x1*x2*x3^3-x16; m(13)=exp((203172+21.38*1863)/8.314*1863)*x1^3*x4^2-x17; m(14)=exp((75158+21.38*1863)/8.314*1863)*x1*x4-x18; m(15)=exp((313417-282.66*1863)/8.314*1863)*x1*x2*x3-x19; m(16)=(0.5*x1+3*x5+12*x6+x7+x8+x9+3*x10+2*x11+3*x12+x13+2*x15+x16+3*x17+x18+x19)-0.25; m(17)=(x2+x5+7*x6+x7+2*x8+6*x9+3*x14+x15+x16)-0.25; m(18)=(x3+x10+x11+2*x12+x13+2*x14+x15+2*x16+x19)-0.25; m(19)=(x4+2*x17+x18+x19)-0.25; end x=fsolve(@slag,[0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.0,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5],optimset('display','off')) Error using trustnleqn (line 28) Objective function is returning undefined values at initial point. FSOLVE cannot continue. Error in fsolve (line 376) [x,FVAL,JACOB,EXITFLAG,OUTPUT,msgData]=... |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有8人回复
-
存款400万可以在学校里躺平吗
已经有24人回复
-
最失望的一年
已经有7人回复
-
推荐一本书
已经有16人回复
-
国自然申请面上模板最新2026版出了吗？
已经有20人回复
-
26申博
已经有3人回复
-
请教限项目规定
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
疑惑？
已经有5人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
11楼2013-07-04 15:26:21
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68123.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★
感谢参与，应助指数 +1
xiao__guan: 金币+3, ★★★很有帮助, 非常感谢 2013-07-05 09:33:14
感谢参与，应助指数 +1
xiao__guan: 金币+3, ★★★很有帮助, 非常感谢 2013-07-05 09:33:14
|
程序中，end是多余的。修改后如下，可运行，不过结果不好，也许和初值有关。以下程序复制进一个m文件中运行即可。 function solnol clear all;clc x0=[0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.0,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5]; x=fsolve(@slag,x0,optimset('display','off')) function m=slag(n) x1=n(1); x2=n(2); x3=n(3); x4=n(4); x5=n(5); x6=n(6); x7=n(7); x8=n(8); x9=n(9); x10=n(10); x11=n(11); x12=n(12); x13=n(13); x14=n(14); x15=n(15); x16=n(16); x17=n(17); x18=n(18); x19=n(19); m(1)=exp((13836+25.44*1863)/8.314*1863)*x1^3*x2-x5; m(2)=exp((35288+204.59*1863)/8.314*1863)*x1^12*x2^7-x6; m(3)=exp((18946+25.44*1863)/8.314*1863)*x1*x2-x7; m(4)=exp((18566+23.68*1863)/8.314*1863)*x1*x2^2-x8; m(5)=x9*0; m(6)=exp((107800+13.44*1863)/8.314*1863)*x1^3*x3^2-x10; m(7)=exp((122731+8.85*1863)/8.314*1863)*x1^2*x3-x11; m(8)=exp((302587+66.24*1863)/8.314*1863)*x1^3*x3^2-x12; m(9)=exp((86231+0.995*1863)/8.314*1863)*x1*x3-x13; m(10)=exp((19512+28.19*1863)/8.314*1863)*x2^3*x3^2-x14; m(11)=exp((354932-417.63*1863)/8.314*1863)*x1^2*x2*x3-x15; m(12)=exp((1252543-277.06*1863)/8.314*1863)*x1*x2*x3^3-x16; m(13)=exp((203172+21.38*1863)/8.314*1863)*x1^3*x4^2-x17; m(14)=exp((75158+21.38*1863)/8.314*1863)*x1*x4-x18; m(15)=exp((313417-282.66*1863)/8.314*1863)*x1*x2*x3-x19; m(16)=(0.5*x1+3*x5+12*x6+x7+x8+x9+3*x10+2*x11+3*x12+x13+2*x15+x16+3*x17+x18+x19)-0.25; m(17)=(x2+x5+7*x6+x7+2*x8+6*x9+3*x14+x15+x16)-0.25; m(18)=(x3+x10+x11+2*x12+x13+2*x14+x15+2*x16+x19)-0.25; m(19)=(x4+2*x17+x18+x19)-0.25; |
2楼2013-07-03 11:32:44
3楼2013-07-03 15:44:41
4楼2013-07-03 15:46:50