Question

这道matlab用ODE45怎么写？

Answer 1

该常微分方程用ode45函数求解。可以这样处理：

y0=[100,100];

[t,x]=ode45(@odefun,[0 100],y0); %odefun（）是微分方程组的自定义函数

figure(1)

plot(t,x),grid on

legend('x1(t)','x2(t)','Location','northwest')

xlabel('t'),ylabel({'x1(t)','x2(t)'})

figure(2)

plot(x(:,1),x(:,2)),grid on

xlabel('x1(t)'),ylabel('x2(t)')

title('x1(t) 与x2(t) 相平面图')

运行结果