Phase Plot - Velocity vs Position

6 views (last 30 days)
Martin Romero
Martin Romero on 21 Apr 2019
Answered: Star Strider on 21 Apr 2019
How can I create a graph of velocity vs position(x' vs x) (phase plot).
Here is the current code I am working with..
function LAB08ex1
m = 1; %mass [kg]
k = 4; %spring constant [N/m]
omega0 = sqrt(k/m);
y0 = 0.1; v0 = 0; %initial conditions
[t,Y] = ode45(@f,[0,10],[y0,v0],[],omega0); % solve for 0<t<10
y = Y(:,1); v = Y(:,2); % time series for y and v
E = (1/2*m*v.^2)+(1/2*k*y.^2); %calculating the energy
figure(1); plot(t,y, 'b+-',t,v,'ro-', t,E, 'gro-'); % retrieve y, v from Y
%-----------------------------------------x
function dYdt = f(t,Y,omega0)
y = Y(1); v = Y(2);
dYdt = [ v ; -omega0^2*y ];

Answers (1)

Star Strider
Star Strider on 21 Apr 2019
Unless I am not understanding your code, that would appear to be:
figure(2)
plot(y, v)
grid
xlabel('Position')
yleabel('Velocity')

Categories

Find more on MATLAB 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!