| 5 | 1/1 | 返回列表 |
| 查看: 1627 | 回復(fù): 4 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
周寧_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ù)求解的時(shí)候3個(gè)未知數(shù)只能編輯三個(gè)方程,否則會(huì)提示W(wǎng)arning: Trust-region-dogleg algorithm of FSOLVE cannot handle non-square systems; using Levenberg-Marquardt algorithm instead. 用前三個(gè)方程來解方程,x的解和所選的初值完全相同。請(qǐng)問這是什么原因,程序如下請(qǐng)問解出來的解對(duì)不對(duì)。 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 方法來解這個(gè)方程組,程序如下 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) 運(yùn)行時(shí)顯示Error: File: Untitled4.m Line: 1 Column: 2 An indexing expression on the left side of an assignment must have at least one subscript. 我對(duì)matlab不懂,請(qǐng)大神幫忙解一下這個(gè)方程組,望附上程序,謝謝! |
新蟲 (初入文壇)
鐵桿木蟲 (職業(yè)作家)
木蟲 (知名作家)
木蟲 (知名作家)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京理工大學(xué)本科211材料工程294求調(diào)劑 +6 | mikasa的圍巾 2026-03-28 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑310 +7 | 溫柔的晚安 2026-03-25 | 8/400 |
|
|
[考研] 332求調(diào)劑 +4 | @MZB382400 2026-03-28 | 4/200 |
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
[考研] 085404求調(diào)劑,總分309,本科經(jīng)歷較為豐富 +4 | 來財(cái)aa 2026-03-25 | 4/200 |
|
|
[考研] 一志愿鄭大085600,310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap- 2026-03-22 | 10/500 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 334分 一志愿武理-080500 材料求調(diào)劑 +4 | 李李不服輸 2026-03-25 | 4/200 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生(材料、化學(xué)、電化學(xué),環(huán)境) +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境,總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|