声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1272|回复: 0

[编程技巧] 关于参数变化的常微分方程求解

[复制链接]
发表于 2012-3-27 04:20 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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);
回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-24 19:30 , Processed in 0.052292 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表