如何在MATLAB中绘制非线性微分方程系统？ [重复]

Question

这个问题在这里已有答案：

实际上，我又问了同样的问题

DX = Y

的Dy = -k * YX ^ 3个+ 9.8 * COS（t）的

inits =（ 'X（0）= 0，Y（0）= 0'）

这些是我想绘制的微分方程。

首先，我试图解决微分方程，然后绘制图形。

Dsolve（'Dx = y'，'Dy = -k * yx ^ 3 + 9.8 * cos（t）'，inits）

但是，像这样，这个系统没有明确的解决方案。

现在我卡住了:(

如何在不解决方程式的情况下绘制该系统？

克里斯泰勒通过绘制y vs t和x vs t图来解决我的问题。

但是我想知道绘制y对x图的方法（这是对的吗？我的意思是垂直轴y和水平轴x的图形）

请帮我

Answer 1

使用Chris Taylor的答案如下：

k = 1;
f = @(t,x) [x(2); -k * x(2) - x(1)^3 + 9.8 * cos(t)];
tspan = [0, 10];
xinit = [0, 0];
[t, y] = ode45(f, tspan, xinit);

给出值x和y作为t的函数，从0到10，分别存储在y的第一列和第二列中。 因此，您可以通过简单地调用来绘制y（x）

plot(y(:,2),y(:,1));