| 查看: 1073 | 回复: 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 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有149人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有186人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有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