24小时热门版块排行榜    

查看: 3269  |  回复: 5

mathsxdg

金虫 (正式写手)

[求助] 如何用MATLAB求解脉冲微分方程并画图

有以下脉冲微分方程,请高手赐教如何用Matlab求解并画图,多谢了!!!
如何用MATLAB求解脉冲微分方程并画图


l: 如何用MATLAB求解脉冲微分方程并画图-1

其中,最后一个式子表示脉冲线l, b0 为一常数。

[ Last edited by mathsxdg on 2013-10-15 at 21:23 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

生命科学
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

niubenhit

金虫 (正式写手)

【答案】应助回帖

你这里 l 是一点吧,可以这样求:
编一个循环
每个循环长度为 l ,在其中求解ode45的方程,然后将解得最后一点作为下一个循环的初值。
即可。

程序:
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];
有眼大如天,山高月更阔
2楼2013-10-26 08:21:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

niubenhit

金虫 (正式写手)

引用回帖:
2楼: Originally posted by niubenhit at 2013-10-26 08:21:08
你这里 l 是一点吧,可以这样求:
编一个循环
每个循环长度为 l ,在其中求解ode45的方程,然后将解得最后一点作为下一个循环的初值。
即可。

程序:
function draw()
l=10;
x10=1;x20=2;b0=0.2;
for j=1 ...

没看到最后一行,不好意思,楼主,上面程序有误,删了吧。
有眼大如天,山高月更阔
3楼2013-10-26 08:22:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

niubenhit

金虫 (正式写手)

引用回帖:
3楼: Originally posted by niubenhit at 2013-10-26 08:22:26
没看到最后一行,不好意思,楼主,上面程序有误,删了吧。...

有个想法,可以用定步长R-K差分格式去写,这样脉冲线的位置程序中比较好确定。试一下吧。
有眼大如天,山高月更阔
4楼2013-10-26 08:24:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathsxdg

金虫 (正式写手)

引用回帖:
4楼: Originally posted by niubenhit at 2013-10-26 08:24:06
有个想法,可以用定步长R-K差分格式去写,这样脉冲线的位置程序中比较好确定。试一下吧。...

多谢您的回复,让您费心了。这是一个状态依赖的脉冲微分方程,脉冲线为一条过原点的半直线。如果您知道如何编写程序,麻烦您再帮我一下,再次感谢!
生命科学
5楼2013-10-26 08:53:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Master121

金虫 (小有名气)

开心就好!

老师您好!这个程序写出来了吗?能否借鉴学习下?
6楼2017-08-25 13:22:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mathsxdg 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 306求调剂 +4 chuanzhu川烛 2026-03-18 4/200 2026-03-21 08:25 by laoshidan
[考研] 南昌大学材料专硕311分求调剂 +6 77chaselx 2026-03-20 6/300 2026-03-21 07:24 by JourneyLucky
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 083200学硕321分一志愿暨南大学求调剂 +3 innocenceF 2026-03-17 3/150 2026-03-21 02:35 by JourneyLucky
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 271材料工程求调剂 +8 .6lL 2026-03-18 8/400 2026-03-21 00:58 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
信息提示
请填处理意见