| 5 | 1/1 | 返回列表 |
| 查看: 1629 | 回復(fù): 4 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
周寧_fish新蟲 (初入文壇)
|
[求助]
matlab解方程組求助攻 已有2人參與
|
|
|
方程組如下,求解x y z f1=0.00124*y+0.114094*z-0.1*y+0.1*y*exp(-900*x)=0 f2=0.00549*y+0.648922*z-0.1*y+0.1*y*exp(-1800*x)=0 f3=0.00607*y+0.750408*z-0.1*y+0.1*y*exp(-2700*x)=0 f4=0.00821*y+1.251459471*z-0.1*y+0.1*y*exp(-3600*x)=0 還有f5... 用matlab的fsolve函數(shù)求解的時候3個未知數(shù)只能編輯三個方程,否則會提示W(wǎng)arning: Trust-region-dogleg algorithm of FSOLVE cannot handle non-square systems; using Levenberg-Marquardt algorithm instead. 用前三個方程來解方程,x的解和所選的初值完全相同。請問這是什么原因,程序如下請問解出來的解對不對。 function q=myfun(p) x=p(1); y=p(2); z=p(3); q(1)=0.00124*y+0.114094*z-0.1*y+0.1*y*exp(-900*x); q(2)=0.00549*y+0.648922*z-0.1*y+0.1*y*exp(-1800*x); q(3)=0.00607*y+0.750408*z-0.1*y+0.1*y*exp(-2700*x); x = fsolve('myfun',[10,100000,3000]',optimset('Largescale','off')) 在網(wǎng)上搜的 Levenberg-Marquardt 方法來解這個方程組,程序如下 f()={"0.00124*y+0.114094*z-0.1*y+0.1*y*exp(-900*x)","0.00549*y+0.648922*z-0.1*y+0.1*y*exp(-1800*x)","0.00607*y+0.750408*z-0.1*y+0.1*y*exp(-2700*x)"} x()={"x","y","z"} x0()={10,100000,3000} d()={0,0,0} erro=1e-8 Loopn=3000 nonlin_LM(f,x,x0,d,erro,Loopn) 運行時顯示Error: File: Untitled4.m Line: 1 Column: 2 An indexing expression on the left side of an assignment must have at least one subscript. 我對matlab不懂,請大神幫忙解一下這個方程組,望附上程序,謝謝! |
木蟲 (知名作家)
鐵桿木蟲 (職業(yè)作家)
新蟲 (初入文壇)
木蟲 (知名作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 277跪求調(diào)劑 +6 | 1915668 2026-03-27 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑310 +10 | 溫柔的晚安 2026-03-25 | 11/550 |
|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +8 | mikasa的圍巾 2026-03-28 | 8/400 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +7 | @taotao 2026-03-29 | 7/350 |
|
|
[考研] 348求調(diào)劑 +5 | 小懶蟲不懶了 2026-03-28 | 5/250 |
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 求化學(xué)調(diào)劑 +4 | wulanna 2026-03-28 | 4/200 |
|
|
[考研] 311求調(diào)劑 +9 | lin0039 2026-03-26 | 9/450 |
|
|
[考研] 0856,材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分,前三科223分 +5 | dongfang59 2026-03-27 | 5/250 |
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|