b********2
发帖数: 5 | 1 我现在对如下的时滞微分方程组写了个小code:
y1'=-a*y1^2-pai*S*(c-y1-y2）+k*y2*y2(t-\tau);
y2'=ip*y2*(c-y1-y2)-k*y2*y2(t-\tau)
(这里\tau是时滞。 y1,y2是变量， a， pai， k， c, ip是参数。
\tau=2, t in [0, 5].)
这个程序运行起来没有问题，但是我不知道我的code是不是正确的描述了这个方程组。
1。 特别是我的ylag1(2）的表达
2。 还有在figure(2) 和figure（4）中，我想分别画 时间t与 y1, y1与y2的图形，
不知
道我的code的表达对不对？ 我用了 sol.y(1,:), sol.y(2,:).这样写是正确的吗？
谢了
。
function v=manuf(t, Z, y)
v=zeros(2, 1);
a=2; pai=3; c=20; k=3; ip=2;
ylag1=Z(:,1);
v(1)=-a*y(1)*y(1)-pai*y(1)*(c-y(1)-y(2))+k*y(2)*ylag1(2);
v(2)=ip*y(2)*(c-y(1)-y(2))-k*y(2)*ylag1(2);
function v=manuh(t)
v(1)=10;
v(2)=14;
opts = ddeset('RelTol',1e-5,'AbsTol',1e-8);
sol=dde23(@manuf, 2, @manuh, [0, 5], opts);
figure(1)
plot(sol.x, sol.y);
figure(2)
plot(sol.x, sol.y(1,:));
figure(3)
plot(sol.x, sol. y(2,:));
figure(4)
plot(sol.y(1,:), sol.y(2,:)) |
|
