24小时热门版块排行榜    

查看: 418  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 孤鸿飘逸 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 gr4vd950ka 2026-07-18 4/200 2026-07-19 21:10 by mxyes0ckns
[硕博家园] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 bjm2e4yior 2026-07-18 4/200 2026-07-19 20:50 by mxyes0ckns
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 ufgqpyvpi6 2026-07-18 4/200 2026-07-19 20:43 by mxyes0ckns
[基金申请] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 ufgqpyvpi6 2026-07-18 5/250 2026-07-19 20:39 by mxyes0ckns
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 ufgqpyvpi6 2026-07-18 6/300 2026-07-19 20:34 by mxyes0ckns
[找工作] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 ufgqpyvpi6 2026-07-18 4/200 2026-07-19 20:23 by mxyes0ckns
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 ufgqpyvpi6 2026-07-18 5/250 2026-07-19 20:19 by mxyes0ckns
[基金申请] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 ufgqpyvpi6 2026-07-18 5/250 2026-07-19 20:14 by mxyes0ckns
[基金申请] 面上的水很深 +17 布布和一二 2026-07-18 22/1100 2026-07-19 20:03 by botar
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 ufgqpyvpi6 2026-07-18 6/300 2026-07-19 16:03 by 0u3sh5hhfi
[基金申请] 国社科通讯评审 +5 潜水冠军123 2026-07-18 6/300 2026-07-19 12:28 by lfy8008
[基金申请] 生命口面上会评结束了吗大家有知道结果吗?现在还没有结果正常不? +6 碧水00 2026-07-16 9/450 2026-07-17 15:21 by jaasxulu
[有机交流] 伯胺上甲基求助 100+4 002tmumu 2026-07-14 7/350 2026-07-16 19:47 by 002tmumu
[论文投稿] 跨出版社商投稿 20+5 倪好520 2026-07-13 5/250 2026-07-16 17:24 by 无视军团
[基金申请] 生生命学部会评 +4 wuke100666 2026-07-15 4/200 2026-07-16 00:12 by 啦啦小妮子
[基金申请] H口,面上,时间戳7月7日。 +3 plumage5 2026-07-15 3/150 2026-07-15 21:06 by 6543yes
[基金申请] E口会评 +11 布布和一二 2026-07-13 13/650 2026-07-15 16:34 by Kittylucky
[论文投稿] 农业机械学报投稿周期 +3 G0DLIN 2026-07-14 7/350 2026-07-15 09:26 by xs74101122
[基金申请] 国自然面上D口祈祷 +4 monkeynana 2026-07-14 4/200 2026-07-15 09:14 by 白水煮
[基金申请] 没收到消息,再次陪跑 +3 天空星辰fly 2026-07-14 3/150 2026-07-14 19:14 by 6543yes
信息提示
请填处理意见