Question

用matlab画出微分方程组dx/dt=-5*x+y+1; dy/dt=x-0.6*y+1;的x-y相图？

微分方程组dx/dt=-5*x+y+1;
dy/dt=x-0.6*y+1;
x0=1，y0=1

Answer 1

题主的问题，可以这样解决：

第一步：自定义微分方程组函数

func=@(t,x)[-5*x(1)+x(2)+1;x(1)-0.6*x(2)+1];

第二步：确定x0的初值及时间范围

x0=[1;1];

tspan=[0,10];

第三步：使用ode45函数，求解t、x、y

[t,x]=ode45(@(t,x)func(t,x),tspan,x0);

第四步：使用plot函数，绘制x-y相位图

plot(x(:,1),x(:,2),'.-')

xlabel('x'),ylabel('y')

第五步：运行程序，得到如下结果