24小时热门版块排行榜

返回列表

当前主题已经存档。

【有奖交流】积极回复本帖子，参与交流，就有机会分得作者孤鸿飘逸的 58 个金币

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

孤鸿飘逸

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 5861.4
散金: 499
红花: 17
帖子: 617
在线: 432.5小时
虫号: 864807
注册: 2009-10-07
性别: GG
专业: 机器人学及机器人技术

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

跪求MATLAB关于脉冲微分方程程序

回复此楼

» 猜你喜欢

参与限项已经有3人回复
假如你的研究生提出不合理要求已经有7人回复
实验室接单子已经有4人回复
全日制（定向）博士已经有4人回复
对氯苯硼酸纯化已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有12人回复
不自信的我已经有12人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复
北核录用已经有3人回复

1楼 2010-03-27 13:09:31

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

孤鸿飘逸

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 5861.4
散金: 499
红花: 17
帖子: 617
在线: 432.5小时
虫号: 864807
注册: 2009-10-07
性别: GG
专业: 机器人学及机器人技术

引用回帖:

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 ...

能详细点吗？谢谢

赞一下(1人)

回复此楼

3楼2010-03-27 15:07:46

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

查看全部 3 个回答

change0618

铁杆木虫 (著名写手)

方丈大师

应助: 44 (小学生)
金币: 17724.5
红花: 17
帖子: 2413
在线: 546.7小时
虫号: 496517
注册: 2008-01-19
专业: 化学反应工程

★ ★
孤鸿飘逸(金币+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];

赞一下(1人)

回复此楼

2楼2010-03-27 14:56:54

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

查看全部 3 个回答