How to plot phase portrait given a system of equations
29 views (last 30 days)
Show older comments
Missael Hernandez
on 12 Feb 2021
Commented: prajith samraj
on 18 Aug 2021
Suppose I have the system represented as
How would I plot the phase portriats for A=0.5,1, and 3? I know how to do it if I was given the equations, but I am confused how I would do it given matrix notation. Thank you.
0 Comments
Accepted Answer
Alan Stevens
on 12 Feb 2021
Do you mean something like this
A = 0.5;
M = [A 0; 1 A];
dXdt = @(t,X) M*X;
tspan = [0 1];
X0 = [1; 1];
[t, X] = ode45(dXdt, tspan, X0);
x1 = X(:,1); x2 = X(:,2);
V = dXdt(t,X')';
v1 = V(:,1); v2 = V(:,2);
figure
plot(t,x1,t,x2), grid
xlabel('t'),ylabel('x')
legend('x1','x2')
figure
subplot(2,1,1)
plot(x1,v1),grid
xlabel('x1'),ylabel('v1')
subplot(2,1,2)
plot(x2,v2),grid
xlabel('x2'),ylabel('v2')
3 Comments
prajith samraj
on 18 Aug 2021
x'' + a x' +bx+cx^3 +c x5= f coswt
how would you plot phase portrait . i have some parameter from above equations and i have choose the control parameter f. how to write phase portrait code in matlab.
More Answers (0)
See Also
Categories
Find more on Ordinary Differential Equations in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!