| 查看: 424 | 回复: 2 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 孤鸿飘逸 的 58 个金币 | |||
孤鸿飘逸金虫 (正式写手)
|
[交流]
【求助】跪求MATLAB关于脉冲微分方程程序
|
||
| 跪求MATLAB关于脉冲微分方程程序 |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有13人回复
-
关于代码变化问题，想知道的进来
已经有20人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
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
孤鸿飘逸(金币+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
孤鸿飘逸
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 5861.4
- 散金: 499
- 红花: 17
- 帖子: 617
- 在线: 432.5小时
- 虫号: 864807
- 注册: 2009-10-07
- 性别: GG
- 专业: 机器人学及机器人技术
3楼2010-03-27 15:07:46
