24小时热门版块排行榜    

查看: 3462  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见