作业帮 > 综合 > 作业

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

来源:学生作业帮 编辑:搜搜做题作业网作业帮 分类:综合作业 时间:2024/06/22 13:53:05
如何用matlab解微分方程组 Dx=-ay(-bx) Dy=-cx(-dy) x(0)=x0 y(0)=y0 并在x、y相平面内作图
如何用matlab解微分方程组 Dx=-ay(-bx) Dy=-cx(-dy) x(0)=x0 y(0)=y0 并在x、
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