| 查看: 364 | 回复: 2 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 孤鸿飘逸 的 58 个金币 | |||
孤鸿飘逸金虫 (正式写手)
|
[交流]
【求助】跪求MATLAB关于脉冲微分方程程序
|
||
| 跪求MATLAB关于脉冲微分方程程序 |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有6人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
论文终于录用啦！满足毕业条件了
已经有28人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
-
实验室接单子
已经有3人回复
-
磺酰氟产物，毕不了业了！
已经有8人回复
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