| 5 | 1/1 | 返回列表 |
| 查看: 2137 | 回復(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),謝謝! 有的前輩可能說把v及u的導(dǎo)數(shù)具體算出來,然后利用ode45來求解,可是這樣來算是特別麻煩,并且結(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ā)表 | |
|---|---|---|---|---|
|
[考研]
|
nnnnnnn5 2026-03-25 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
[考研] 一志愿華東理工大學(xué)081700,初試分?jǐn)?shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +4 | @taotao 2026-03-26 | 5/250 |
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
[考研] 求調(diào)劑 一志愿 本科 北科大 化學(xué) 343 +6 | 13831862839 2026-03-24 | 7/350 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +4 | 柳溪邊 2026-03-26 | 4/200 |
|
|
[考研] 打過很多競(jìng)賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 考研調(diào)劑 +6 | 來好運(yùn)來來來 2026-03-21 | 7/350 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|