24小时热门版块排行榜    

查看: 444  |  回复: 2
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 孤鸿飘逸 的 58 个金币

孤鸿飘逸

金虫 (正式写手)

[交流] 【求助】跪求MATLAB关于脉冲微分方程程序

跪求MATLAB关于脉冲微分方程程序
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

★ ★
孤鸿飘逸(金币+2): 2010-03-27 15:06
adu886886(金币+2):感谢应助 2010-03-27 20:59
Differential equations:
dx1/dt=-x1/t1+x2+x3*(x2-x1)
dx2/dt=r(t)-x2-x3*(x2-x1)
dx3/dt=A*x2+x3*(x2-x1)-x3/t0

initial conditions:
x1(0)=x2(0)=x3(0)=0;
r(t) is as follow:
mr=mod(t,tp+ti);
r=0;
if mr<2;
    r=4;
end
here is :
tp=2;ti=50;t0=10;t1=0.5;A=1.5e-8;
range of time: from 0 to 200.
the program is as follow:

程序为:

function dx=cyj(t,x);
% 脉冲作用下的解微分方程
% Example:
%     [t,y]=ode45('cyj',[0,200],[0,0,0]);
%     plot(t,y(:,1),t,y(:,2),t,y(:,3));
%     legend('x1','x2','x3',0);
%
tp=2;
ti=50;
mr=mod(t,tp+ti);
r=0;
if mr<2;
    r=4;
end
t0=10;
t1=0.5;
A=1.5e-8;
Ji=x(3)*(x(2)-x(1));
dx=[-x(1)/t1+x(2)+Ji;r-x(2)-Ji;A*x(2)+Ji-x(3)/t0];
2楼2010-03-27 14:56:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孤鸿飘逸

金虫 (正式写手)

引用回帖:
Originally posted by change0618 at 2010-03-27 14:56:54:
Differential equations:
dx1/dt=-x1/t1+x2+x3*(x2-x1)
dx2/dt=r(t)-x2-x3*(x2-x1)
dx3/dt=A*x2+x3*(x2-x1)-x3/t0

initial conditions:
x1(0)=x2(0)=x3(0)=0;
r(t) is as follow:
mr=mod(t,tp+ti ...

能详细点吗?谢谢
3楼2010-03-27 15:07:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 孤鸿飘逸 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 18:07 by hZiFeudZyoGR
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 18:07 by hZiFeudZyoGR
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 17:25 by hZiFeudZyoGR
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ASdOkHsho7FD 2026-08-28 14/700 2026-08-30 17:01 by hZiFeudZyoGR
[基金申请] 有没有仍没收到信息的 +6 德尚中行 2026-08-27 7/350 2026-08-30 12:53 by Bxndjjd
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:03 by l0VvVHGBGRLv
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +8 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:40 by l0VvVHGBGRLv
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 06:57 by ZPa0EcMwuECS
[基金申请] 中青基了要发朋友圈吗? +4 349506619 2026-08-28 4/200 2026-08-29 22:41 by alongwaytogo
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
信息提示
请填处理意见