| 查看: 1010 | 回复: 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 |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有188人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有13人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助：如何用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
5