实际上,我再次提出同样的问题
DX = Y
的Dy = -k * Y-X ^ 3 + 9.8 * COS(t)的
inits =( 'X(0)= 0,Y(0)= 0')
这些是我想绘制的微分方程。
首先,我尝试解决微分方程,然后绘制图形。
Dsolve('Dx = y','Dy = -k * y-x ^ 3 + 9.8 * cos(t)',inits)
然而,像这样,这个系统没有明确的解决方案。现在我被卡住了:(
如何在不解决方程式的情况下绘制该系统?
Chris Taylor通过绘制y vs t和x vs t图来解决我的问题。
但是我想知道绘制y对x图的方法(这是对的吗?我的意思是垂直轴y和水平轴x的图)
请帮帮我
答案 0 :(得分:3)
使用Chris Taylor的answer代码如下:
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));