Matlab使用前向欧拉方法

Question

DT / DX = X（1.15-0.02x-0.02y） DT / DY = Y（1.34-0.02y-0.04X）

T=15;
dt=1;
N=floor(T/dt);
t(1)=1;
x(1)=21;
y(1)=20;

for i=1:N   
x(i+1)=x(i)+dt*x(i)*(1.15-0.02*x(i)-0.02*y(i));
y(i+1)=y(i)+dt*y(i)*(1.34-0.02*y(i)-0.04*x(i));
t(i+1)=i*dt;
end
title('state variable')
subplot(2,1,1)
plot(t,x,'r*-')
subplot(2,1,2)
plot(t,y,'bo-')

使用前进euler方法 这是对的吗？

Answer 1

我是这么认为的，但是你的等式中有一个错误

dx/dt=...
dy/dt=...