| 查看: 1043 | 回复: 1 | ||
mathda新虫 (小有名气)
|
[求助]
matlab 求脉冲微分方程
|
|
程序代码运行出错，求教 function predprey a=0.25;b=-0.01;c=-1.0;d=0.01; impcount=0; Y=-1.2*c/d; p=0.2; tstar=0; tfinal=100; while true if (impcount ==0 ) options = ddeset('Events',@events,'AbsTol',1e-9,... 'RelTol',1e-9); [t,y] = ode15s(@ddes,[tstar,tfinal],[80,30], options); plot(y(:,1),y(:,2)); else % Specify the new solution at impulse times... options = ddeset(options, 'InitialY',... [(1-p).*y(:,1); y(:,2)]); [t,y] = ode15s(@ddes, [tstar,tfinal],... [(1-p).*y(:,1);y(:,2)],options); plot(y(:,1),y(:,2)); end tstar = max(t); impcount = impcount + 1; if (tstar >= tfinal) break; end end %===Nested functions============================== % Evaluate the DDES... function dydt = ddes(t,y) dydt = [a*y(1)+b*y(1)*y(2); ... c*y(2)+d*y(1)*y(2)]; end % Evaluate the event function residuals... function [value,isterminal,direction] = events(t,y) value = y(1)-Y; isterminal = 1; direction = 1; end end |
回复此楼» 猜你喜欢
大豆异黄酮分离
已经有0人回复
-
湖南大学材料学院急招2026年博士生，临时增加一名博士联培指标
已经有10人回复
-
物理学I论文润色/翻译怎么收费?
已经有146人回复
-
天津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生1-2名
已经有1人回复
-
中国科学院物理研究所谌志国研究员团队招收2027年博士研究生
已经有5人回复
-
2026年中德博士后交流项目 - 新型量子和磁性材料：材料制备表征和中子散射研究
已经有12人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有1人回复
-
如何从铁电相到顺电相。
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：如何用MATLAB的偏微分方程工具箱画出不同区域内的温度分布图？
已经有5人回复
- matlab求解偏微分方程组
已经有3人回复
- 求助matlab求解复杂微分方程
已经有7人回复
- 高手进来用matlab求解常微分方程，十分感谢！
已经有17人回复
- matlab微分方程组参数拟合，以周为单位求解，汇总后以年为单位进行数值比较
已经有4人回复
- matlab 线性常微分方程求解
已经有5人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 高金求助matlab解微分方程组
已经有12人回复
- matlab解微分方程组
已经有15人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
2楼2012-11-01 16:30:44