| 查看: 995 | 回复: 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人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有136人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：如何用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