%   FY1001/TFY4145 Mekanisk fysikk november 2011
%   Oving 13
%   Faseplott for udempet og dempet harmonisk oscillator
%
%   Date: 22.11.2011
%   Author: Jon Andreas Stovneng

% Tabell med tidspunkter fra 0 til 100, i alt 1001 verdier
tid = linspace(0,100,1001);

figure(1);

subplot(1,3,1);
% Udempet oscillator, A = 1 og 0.5, omega = 1, m = 1
% slik at x(t)=sin(t) resulterer i p(t)=mv(t)=cos(t):
x11 = sin(tid);
p11 = cos(tid);
x12 = 0.5.*x11;
p12 = 0.5.*p11;
plot(x11,p11,x12,p12);
xlabel('x');
ylabel('p');
axis([-1 1 -1 1]);
title('Udempet');

subplot(1,3,2);
% Dempet oscillator, underkritisk. Her er omega=1 og beta=0.15:
beta = 0.15;
x21 = exp(-beta.*tid).*cos(tid);
p21 = exp(-beta.*tid).*(-beta.*cos(tid)-sin(tid));
x22 = 0.5*exp(-beta.*tid).*cos(tid);
p22 = 0.5*exp(-beta.*tid).*(-beta.*cos(tid)-sin(tid));
plot(x21,p21,x22,p22);
xlabel('x');
ylabel('p');
axis([-1 1 -1 1]);
title('Underkritisk dempet');

subplot(1,3,3);
% Dempet oscillator, overkritisk. Her er beta=1 og omega=0.2:
omega = 0.2;
x3 = exp(-tid).*cosh(omega.*tid);
p3 = exp(-tid).*(omega.*sinh(omega.*tid)-cosh(omega.*tid));
plot(x3,p3);
xlabel('x');
ylabel('p');
axis([-1 1 -1 1]);
title('Overkritisk dempet');
