马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
需要模拟一个脉冲微分方程组,无脉冲方程简要为: dx/dt=r(K-x)-axy, dy/dt=qaxy-(m+d)y..
现在需要进行模拟,其中的r,K,a,q,d均为常数,而m的情况比较复杂,当x/y小于1时,m=0.1*(y/x-1);当x/y大于等于1时,m=0. 求助如何在编写DHTI时候体现参数m的变化,我写了一个,但有些问题,如下:
function DHTI=Model20120326000(t,x)
global r m a K q d;
if x(1)/x(2)<1;
m=0.1*((x(2)/x(1))-1);
else
m=0;
end
DHTI=[r*(K-x(1))-a*x(1)*x(2);
q*a*x(1)*x(2)-x(2)*(m+d); |