How do you plot the direction (vector) field of a second-order homogeneous ode using Matlab?
I've already used MATLAB to check the solution to the ode and I've tried to use tutorials online to plot the direction (vector) field, but haven't had any luck. Here's what I have done in MATLAB:
eqn1 = 'D2x + 5*Dx + 4*x = 0';
x = dsolve(eqn1, 't')
The above gives me the correct solution to the second-order ode, but isn't helpful for plotting the direction (vector) field. I'm new to MATLAB, so any help would be greatly appreciated.
Best Answer
The second order ODE $X'' +5X'+4X =0$ rewrites as a first-order ODE system. Indeed, setting $(x,y)=(X,X')$, one has \begin{aligned} \frac{\text{d}x}{\text{d}t} &= y \, ,\\ \frac{\text{d}y}{\text{d}t} &= -5y-4x \, . \end{aligned} The phase-space plot can then be obtained in a similar manner as described in this post, e.g. using Matlab's
quiver
function.