Question

如何用matlab解微分方程组 Dx=-ay(-bx) Dy=-cx(-dy) x(0)=x0 y(0)=y0 并在x、y相平面内作图

Answer 1

1
将下面的代码保存为M文件：

function ccc

x0 = 0.1;
y0 = 0.2;

z0 = [x0;y0]; % 初值
[T,Z]=ode45(@cc,[0,2],z0);
plot(Z(:,1),Z(:,2))
title('运动轨迹')
xlabel('X')
ylabel('Y')

function dz = cc(t,z)

a = 1;
b = 2;
c = 1.3;
d = 1.4;

x = z(1);
y = z(2);

dz = [-a*y*(-b*x);-c*x*(-d*y)];

2
将下面的代码复制到命令窗口，按回车键，运行即可
ccc