| 5 | 1/1 | 返回列表 |
| 查看: 1603 | 回復(fù): 4 | ||
周寧_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è)方程,否則會提示W(wǎng)arning: Trust-region-dogleg algorithm of FSOLVE cannot handle non-square systems; using Levenberg-Marquardt algorithm instead. 用前三個(gè)方程來解方程,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 方法來解這個(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. 我對matlab不懂,請大神幫忙解一下這個(gè)方程組,望附上程序,謝謝! |
鐵桿木蟲 (職業(yè)作家)
新蟲 (初入文壇)
木蟲 (知名作家)
木蟲 (知名作家)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 324分 +5 | llllkkkhh 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[教師之家] 焦慮 +8 | 水冰月月野兔 2026-03-13 | 11/550 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] [0860]321分求調(diào)劑,ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
[考研] 材料301分求調(diào)劑 +5 | Liyouyumairs 2026-03-12 | 5/250 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|
|
[考研] 321求調(diào)劑(食品/專碩) +3 | xc321 2026-03-12 | 6/300 |
|