| 5 | 1/1 | 返回列表 |
| 查看: 3261 | 回復: 5 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
mathsxdg金蟲 (正式寫手)
|
[求助]
如何用MATLAB求解脈沖微分方程并畫圖
|
|
|
有以下脈沖微分方程,請高手賜教如何用Matlab求解并畫圖,多謝了!! l: 其中,最后一個式子表示脈沖線l, b0 為一常數。 [ Last edited by mathsxdg on 2013-10-15 at 21:23 ] |


|
你這里 l 是一點吧,可以這樣求: 編一個循環(huán) 每個循環(huán)長度為 l ,在其中求解ode45的方程,然后將解得最后一點作為下一個循環(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) %%以下為參數,可修改。 b0=1; %%以上為參數,可修改。 x1=x(1);x2=x(2); dydt =[2*x1-x2 x1+2*x2]; |


金蟲 (正式寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 295求調劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 08工科 320總分 求調劑 +5 | 梨花珞晚風 2026-03-17 | 5/250 |
|
|
[考研] 070300化學319求調劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
[考研] 304求調劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
[考研] 296求調劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
[考研] 301求調劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 材料工程專碩調劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 318求調劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 321求調劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 0856求調劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
[考研] 288求調劑 +4 | 奇點0314 2026-03-14 | 4/200 |
|
|
[考研] 255求調劑 +3 | 李嘉慧, 2026-03-12 | 4/200 |
|
|
[考研] 330求調劑 +3 | ?醬給調劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 求材料調劑 085600英一數二總分302 前三科235 精通機器學習 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
[考研] 311求調劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[考研] 工科材料085601 279求調劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
[考研] 290求調劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
[考研] 274求調劑 +3 | S.H1 2026-03-12 | 3/150 |
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|