| 5 | 1/1 | 返回列表 |
| 查看: 2138 | 回復(fù): 8 | ||||
| 本帖產(chǎn)生 1 個(gè) 仿真EPI ,點(diǎn)擊這里進(jìn)行查看 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
bdwhq新蟲 (著名寫手)
|
[求助]
請(qǐng)各位前輩指點(diǎn)一下,如何用MATLAB求解如下的微分方程,謝謝!
|
|||
|
請(qǐng)各位前輩指點(diǎn)一下,如何用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'表示各個(gè)變量對(duì)時(shí)間t的導(dǎo)數(shù),X_1(0),x_2(0),w(0),v(0),也就是各個(gè)變量在t=0時(shí)的初值是給定的,比如X_1(0)=0.1,x_2(0)=0.1,w(0)=0,v(0)=0,也可以是其他的值。想畫出t從0到20秒x_1,x_2的圖像,敬請(qǐng)大家給予指點(diǎn),謝謝! 有的前輩可能說(shuō)把v及u的導(dǎo)數(shù)具體算出來(lái),然后利用ode45來(lái)求解,可是這樣來(lái)算是特別麻煩,并且結(jié)果還不理想。煩請(qǐng)各位看看能否好的辦法,小弟在此謝謝各位了。! 若懸賞金幣少,可追加。 |
新蟲 (著名寫手)
|
本帖內(nèi)容被屏蔽 |
|
本帖內(nèi)容被屏蔽 |
榮譽(yù)版主 (著名寫手)
![]() |
專家經(jīng)驗(yàn): +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 代碼 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 求調(diào)劑 +7 | 爭(zhēng)取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
[考研]
|
平樂(lè)樂(lè)樂(lè) 2026-03-26 | 4/200 |
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分,求調(diào)劑 +6 | jiajunX 2026-03-22 | 6/300 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[考研] 070300,一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|