| 查看: 2057 | 回复: 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]=... |
回复此楼» 猜你喜欢
材料工程085601，270求调剂
已经有42人回复
-
279求调剂
已经有11人回复
-
271求调剂
已经有24人回复
-
297工科调剂?
已经有6人回复
-
材料工程281还有调剂机会吗
已经有35人回复
-
0856专硕求调剂 希望是a区院校
已经有27人回复
-
考研调剂
已经有10人回复
-
求助调剂，跨调
已经有14人回复
-
085600材料与化工，求调剂
已经有14人回复
-
求调剂
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- matlab如何求解一个非线性微分方程组
已经有8人回复
- matlab 计算非线性方程组求解
已经有4人回复
- MATLAB求解非线性方程组
已经有5人回复
- matlab解非线性方程组解不开
已经有3人回复
- matlab求解非线性方程组
已经有3人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- matlab非线性方程组求解，急！！！
已经有5人回复
- 用不动点迭代求非线性方程组，求得的不是想要的解，怎么办
已经有12人回复
- matlab求解非线性方程组，求助！
已经有6人回复
- 如何使用matlab求解非线性方程组的所有整数解？
已经有9人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab求解非线性方程组
已经有16人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- MATLAB非线性方程组求解~~~反应釜模拟
已经有9人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
- 非线性方程组的迭代法（数值计算高手请进）
已经有7人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
4楼2013-07-03 15:46:50
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68129.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
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68129.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
|
新建一个m文件，把以下程序全部复制进去，运行即可。 至于你说赋值1863的问题，一种解决办法就是定义一个全局变量k，对k赋值即可，以下的程序就才用这个办法；更简洁一些的话，直接把k=1863写到“function m=slag(n)”的下面。 function solnol clear all;clc global k k=1863; 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) global k 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*k)/8.314*1863)*x1^3*x2-x5; m(2)=exp((35288+204.59*k)/8.314*1863)*x1^12*x2^7-x6; m(3)=exp((18946+25.44*k)/8.314*1863)*x1*x2-x7; m(4)=exp((18566+23.68*k)/8.314*1863)*x1*x2^2-x8; m(5)=x9*0; m(6)=exp((107800+13.44*k)/8.314*k)*x1^3*x3^2-x10; m(7)=exp((122731+8.85*k)/8.314*k)*x1^2*x3-x11; m(8)=exp((302587+66.24*k)/8.314*k)*x1^3*x3^2-x12; m(9)=exp((86231+0.995*k)/8.314*k)*x1*x3-x13; m(10)=exp((19512+28.19*k)/8.314*k)*x2^3*x3^2-x14; m(11)=exp((354932-417.63*k)/8.314*k)*x1^2*x2*x3-x15; m(12)=exp((1252543-277.06*k)/8.314*k)*x1*x2*x3^3-x16; m(13)=exp((203172+21.38*k)/8.314*k)*x1^3*x4^2-x17; m(14)=exp((75158+21.38*k)/8.314*k)*x1*x4-x18; m(15)=exp((313417-282.66*k)/8.314*k)*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; |
5楼2013-07-03 16:49:39
5