| 5 | 1/1 | 返回列表 |
| 查看: 3273 | 回復(fù): 5 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
mathsxdg金蟲 (正式寫手)
|
[求助]
如何用MATLAB求解脈沖微分方程并畫圖
|
|
|
有以下脈沖微分方程,請(qǐng)高手賜教如何用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]; |



| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 07化學(xué)280分求調(diào)劑 +3 | 722865 2026-03-23 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +3 | 番茄頭--- 2026-03-22 | 4/200 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +3 | WW.' 2026-03-23 | 5/250 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-17 | 10/500 |
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
[考研]
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
[考研] 319求調(diào)劑 +4 | 小力氣珂珂 2026-03-20 | 4/200 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|