| 查看: 2012 | 回復(fù): 10 | |||||
xiao__guan新蟲(chóng) (初入文壇)
|
[求助]
matlab 求解非線性方程組未果
|
|
非軟件專業(yè),只想求一個(gè)非線性方程組的解,就顯示錯(cuò)誤,求指導(dǎo)~..... 程序如下 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]=... |
主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |
|
程序中,end是多余的。修改后如下,可運(yùn)行,不過(guò)結(jié)果不好,也許和初值有關(guān)。以下程序復(fù)制進(jìn)一個(gè)m文件中運(yùn)行即可。 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; |

新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |
|
新建一個(gè)m文件,把以下程序全部復(fù)制進(jìn)去,運(yùn)行即可。 至于你說(shuō)賦值1863的問(wèn)題,一種解決辦法就是定義一個(gè)全局變量k,對(duì)k賦值即可,以下的程序就才用這個(gè)辦法;更簡(jiǎn)潔一些的話,直接把k=1863寫(xiě)到“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; |

主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |

新蟲(chóng) (初入文壇)
新蟲(chóng) (初入文壇)
主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |

新蟲(chóng) (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 086003食品工程求調(diào)劑 +3 | 淼淼111 2026-03-24 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 材料專業(yè)求調(diào)劑 +11 | hanamiko 2026-03-18 | 11/550 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 奶油草莓. 2026-03-22 | 5/250 |
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +7 | 你好明天你好 2026-03-23 | 8/400 |
|
|
[考研] 一志愿中國(guó)石油大學(xué)(華東) 本科齊魯工業(yè)大學(xué) +4 | 石能偉 2026-03-17 | 4/200 |
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 9/450 |
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 085601調(diào)劑 358分 +3 | zzzzggh 2026-03-20 | 4/200 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚(yú) 2026-03-18 | 7/350 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 288求調(diào)劑 +16 | 于海海海海 2026-03-19 | 16/800 |
|
|
[考研]
|
簡(jiǎn)木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3 | 葵梓衛(wèi)隊(duì) 2026-03-20 | 3/150 |
|
|
[考研] 一志愿福大288有機(jī)化學(xué),求調(diào)劑 +3 | 小木蟲(chóng)200408204 2026-03-18 | 3/150 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|