常微分方程的数值解
如
一般地,凡表示未知函数,,叫常微分方程;未知函数是多元函数的,叫做偏微分方程.
Matlab实现
[t,x]=ode45(***@f,ts,x0,options,p1,p2,......)
function dx=f(t,x)
dx=[f1;f2;…;fn]
options=odeset(‘reltol’,rt,’abstol’,at)
默认时rt=10^-3,at=10^-6
食饵—捕食者模型
表示时刻食饵的密度, 表示捕食者的密度;
表示食饵独立生存时的增长率;
表示捕食者独立生存时的死亡率;
表示捕食者的存在对食饵增长的影响系数,反映捕食者对食饵的捕获能力;
表示食饵的存在对捕食者增长的促进系数,反映食饵对捕食者的喂养能力
高阶常微分方程的解法
高阶常微分方程
可以将以上高阶微分方程化为如下一阶常微分方程组
令
常微分方程数值解法 来自淘豆网www.taodocs.com转载请标明出处.