| 9 | 1/1 | 返回列表 |
| 查看: 2117 | 回復(fù): 8 | ||||
| 本帖產(chǎn)生 1 個 仿真EPI ,點擊這里進(jìn)行查看 | ||||
bdwhq新蟲 (著名寫手)
|
[求助]
請各位前輩指點一下,如何用MATLAB求解如下的微分方程,謝謝!
|
|||
|
請各位前輩指點一下,如何用MATLAB求解如下的微分方程,謝謝! 我想用matlab求解下面的微分方程組: x'_1=x_2+x_1*sin(x_1); x'_2=w+x_1*x_2; w'=|v'|*(0.5*v-w)+0.3*v'; 其中V是可以給出的函數(shù):v=-5*(x_2-u)-x_1*x_2,v函數(shù)中的u=-10*x_1-x_1*sin(x_1). 另外,x'_1,x'_2,w',v'表示各個變量對時間t的導(dǎo)數(shù),X_1(0),x_2(0),w(0),v(0),也就是各個變量在t=0時的初值是給定的,比如X_1(0)=0.1,x_2(0)=0.1,w(0)=0,v(0)=0,也可以是其他的值。想畫出t從0到20秒x_1,x_2的圖像,敬請大家給予指點,謝謝! 有的前輩可能說把v及u的導(dǎo)數(shù)具體算出來,然后利用ode45來求解,可是這樣來算是特別麻煩,并且結(jié)果還不理想。煩請各位看看能否好的辦法,小弟在此謝謝各位了。。 若懸賞金幣少,可追加。 |
榮譽版主 (著名寫手)
![]() |
專家經(jīng)驗: +4 |
|
>> x1=sym('x1(t)');x2=sym('x2(t)'); >> vv=-5*(x2+10*x1+x1*sin(x1))-x1*x2 vv = - 50*x1(t) - 5*x2(t) - x1(t)*x2(t) - 5*sin(x1(t))*x1(t) >> syms t >> dv=diff(vv,t) dv = - x1(t)*diff(x2(t), t) - x2(t)*diff(x1(t), t) - 5*sin(x1(t))*diff(x1(t), t) - 50*diff(x1(t), t) - 5*diff(x2(t), t) - 5*cos(x1(t))*x1(t)*diff(x1(t), t) ================================ v‘=-x1*x2'-x2*x1'-5*sin(x1)*x1'-50*x1'-5*x2'-5*cos(x1)*x1*x1' 微分方程組變?yōu)?br /> x1=x(1) x2=x(2) v=x(3) w=x(4) =========== x1’=x2+x1*sin(x1); x2‘=w+x1*x2; v‘=-x1*x2'-x2*x1'-5*sin(x1)*x1'-50*x1'-5*x2'-5*cos(x1)*x1*x1' w'=|v'|*(0.5*v-w)+0.3*v'; matlab 代碼 |

|
本帖內(nèi)容被屏蔽 |
|
本帖內(nèi)容被屏蔽 |
新蟲 (著名寫手)
榮譽版主 (著名寫手)
![]() |
專家經(jīng)驗: +4 |

|
本帖內(nèi)容被屏蔽 |
|
|
本帖內(nèi)容被屏蔽 |
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 294求調(diào)劑材料與化工專碩 +8 | 陌の森林 2026-03-18 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 0817調(diào)劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 311求調(diào)劑 +11 | 冬十三 2026-03-15 | 12/600 |
|
|
[考研] 302求調(diào)劑 +10 | 呼呼呼。。。。 2026-03-17 | 10/500 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 293求調(diào)劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|
|
[考研] 085600材料與化工 309分請求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|