How To Plot Directional Field of 2nd Order Differential Equation IVP

Question

Jordan Stanley 2022-4-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1693950-how-to-plot-directional-field-of-2nd-order-differential-equation-ivp

评论： Jordan Stanley 2022-4-11

采纳的回答： Sam Chak

在 MATLAB Online 中打开

Hello,

I have the second order differential equation initial value problem, y'' + 2y' + y = 0, y(-1) = 0, y'(0) = 0.

In MATLAB, I need to plot the directional field of the solution to the equation without the initial conditions.

I have used the meshgrid() command so far and know that I have to use the quiver() command but I don't know how to enter what I need as parameters to plot the solution.

Here is what I have so far...

% Finds solution to the DE
syms y(x) 
Dy = diff(y); 
D2y = diff(y,2); 
ode = D2y + 2*Dy + y == 0; 
ySol = dsolve(ode) 
 
% Sets up directional field 
[x,y]=meshgrid(-3:0.3:3,-2:0.3:2); 
quiver() %Not sure what to include here.

Any help is appreciated!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Sam Chak 2022-4-11

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1693950-how-to-plot-directional-field-of-2nd-order-differential-equation-ivp#answer_940310

编辑：Sam Chak 2022-4-11

在 MATLAB Online 中打开

@Jordan Stanley

Basically, it should look something like this:

[X, Y] = meshgrid(-3:6/14:3, -3:6/14:3);
U = Y;              % x1' = y'
V = - 2*Y - X;      % x2' = y'' = - 2*y' - y
quiver(X, Y, U, V)
% quiver(X, Y, U, V, 1.5)   % can adjust arrow size
xlabel('x')
ylabel('y')

3 个评论
显示 1更早的评论隐藏 1更早的评论

Sam Chak 2022-4-11

编辑：Sam Chak 2022-4-11

Hi @Jordan Stanley

Well, the ODE can be rewritten in the form of a state-space representation.

Begin by defining

Taking the time derivative yields

.

Rewritting the dynamics in system states

.

Back to quiver function, this quiver(X, Y, U, V) command plots arrows with directional components U and V at the Cartesian coordinates specified by the grid of X and Y values.

The directional components U and V mean the motion of the the point

that extends horizontally according to U vector, and extends vertically according to V vector. Naturally, these imply the first-order equations. So, if we assign

then

.

Hope this explanation is helpful for you to plot the direction field of the desired ODEs.

Jordan Stanley 2022-4-11