| 6 | 1/1 | 返回列表 |
| 查看: 3260 | 回復(fù): 5 | |||
mathsxdg金蟲 (正式寫手)
|
[求助]
如何用MATLAB求解脈沖微分方程并畫圖
|
|
有以下脈沖微分方程,請高手賜教如何用Matlab求解并畫圖,多謝了。! l: 其中,最后一個(gè)式子表示脈沖線l, b0 為一常數(shù)。 [ Last edited by mathsxdg on 2013-10-15 at 21:23 ] |

|
你這里 l 是一點(diǎn)吧,可以這樣求: 編一個(gè)循環(huán) 每個(gè)循環(huán)長度為 l ,在其中求解ode45的方程,然后將解得最后一點(diǎn)作為下一個(gè)循環(huán)的初值。 即可。 程序: function draw() l=10; x10=1;x20=2;b0=0.2; for j=1:10 [t,x]=ode45(@test,[(j-1)*l,j*l],[x10,x20]); [mm,nn]=size(x); plot(t,x);hold on; x10=x(nn,1)+(sqrt(3)/2-1)*x(nn,1)-1/2*b0*x(nn,2); y10=x(nn,2)+1/2*b0*x(nn,1)+(sqrt(3)/2-1)*x(nn,2); x=[]; t=[]; end plot(t,x); function dydt = test(t,x) %%以下為參數(shù),可修改。 b0=1; %%以上為參數(shù),可修改。 x1=x(1);x2=x(2); dydt =[2*x1-x2 x1+2*x2]; |



金蟲 (正式寫手)

金蟲 (小有名氣)
開心就好!
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化工學(xué)碩306求調(diào)劑 +10 | 42838695 2026-03-12 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8 | Ncdx123456 2026-03-13 | 9/450 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +7 | rare12345 2026-03-18 | 7/350 |
|
|
[考博] 環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|