| 查看: 1084 | 回复: 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 |
回复此楼» 猜你喜欢
散金币祈福
已经有187人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有208人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有118人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：如何用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