| 5 | 1/1 | 返回列表 |
| 查看: 2145 | 回復(fù): 8 | |||||
| 本帖產(chǎn)生 1 個 仿真EPI ,點擊這里進行查看 | |||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||||
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)容被屏蔽 |
新蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 309求調(diào)劑 +4 | 呆菇不是戴夫 2026-04-02 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 337求調(diào)劑 +10 | 《樹》 2026-03-29 | 10/500 |
|
|
[考研] 材料求調(diào)劑 +8 | 呢呢妮妮 2026-04-01 | 8/400 |
|
|
[考研] 366求調(diào)劑一志愿東北大學 +3 | 運氣來得若有似?/a> 2026-04-02 | 3/150 |
|
|
[考研] 材料調(diào)劑 +12 | 一樣YWY 2026-04-01 | 12/600 |
|
|
[考研] 286求調(diào)劑 +5 | Sa67890. 2026-04-01 | 7/350 |
|
|
[考研] 350求調(diào)劑 +7 | 阿佳~ 2026-03-31 | 7/350 |
|
|
[基金申請] esi高被引論文是不是能對中標有所加分和幫助呢 +3 | redcom 2026-04-01 | 3/150 |
|
|
[考研] 324求調(diào)劑 +4 | 想上學求調(diào) 2026-04-01 | 5/250 |
|
|
[考研] 358求調(diào)劑 +3 | 王向陽花 2026-03-31 | 3/150 |
|
|
[考研] 335求調(diào)劑 +3 | 321* 2026-03-31 | 4/200 |
|
|
[考研] 080500-315分復(fù)試調(diào)劑 +9 | 上岸3821 2026-03-31 | 9/450 |
|
|
[考研] 340求調(diào)劑 +4 | 希望如此i 2026-03-31 | 4/200 |
|
|
[考研]
|
小羊36 2026-03-30 | 4/200 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學材料與化工方向336分 +13 | 辰沐5211314 2026-03-26 | 13/650 |
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +4 | 崔wj 2026-03-31 | 4/200 |
|
|
[考研] 11408總分309,一志愿東南大學求調(diào)劑,不挑專業(yè) +5 | 天賦帶到THU 2026-03-29 | 6/300 |
|
|
[考研] 284求調(diào)劑 +14 | junqihahaha 2026-03-26 | 15/750 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|