How does one plot phase portraits for systems of differential equations?

Question

Aaron Graf 2020-4-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/519014-how-does-one-plot-phase-portraits-for-systems-of-differential-equations

回答： Brian Hong 2022-1-26

I'm sort of new to this whole process and a lof of my homework for diff eq's asks for us to use technology to plot some of the solutions.

We just started last week on system of diff eq's and I wanted to learn how to do this with MATLAB. Be easy on me as I have very limited knowledge of the language.

The example I have on hand of something I've been asked to plot is the following:

Where x is a function of t.

Any help would be greayly appreciated. Thank in advance!

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

请先登录，再进行评论。

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

Answer 1

Ameer Hamza 2020-4-19

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/519014-how-does-one-plot-phase-portraits-for-systems-of-differential-equations#answer_426893

编辑：Ameer Hamza 2020-4-19

在 MATLAB Online 中打开

See this: http://matlab.cheme.cmu.edu/2011/08/09/phase-portraits-of-a-system-of-odes/ to understand how to make phase portraits in MATLAB. The following shows an example without using for-loop for your equations.

[X1,X2] = meshgrid(-5:0.5:5);
xs = arrayfun(@(x,y) {odeFun([],[x,y])}, X1, X2);
x1s = cellfun(@(x) x(1), xs);
x2s = cellfun(@(x) x(2), xs);
quiver(x1s, x2s)
xlabel('x_1')
ylabel('x_2')
axis tight equal;
function dxdt = odeFun(t,x)
    dxdt(1) = 5*x(1)-9*x(2);
    dxdt(2) = 1*x(1)-5*x(2);
end