%FY1001/TFY4145 Mekanisk fysikk hoesten 2012.
%Oving 12, oppgave 4: Plotting av resonanstopp for
%svakt dempet dreven oscillator.
%
clear all;
%F0 = kraftamplitude (N)
F0 = 10;
%m = masse (kg)
m = 0.1;
%k = fjaerkonstant (N/m)
k = 1000;
%b = dempingskonstant (Ns/m)
b = 0.2;
w0 = sqrt(k/m);
gamma = b/(2*m);
%w = vinkelfrekvensintervallet (1/s)
w = 95:0.01:105;
%Regner ut og plotter A^2(w)
A2 = (F0/m)^2./((w.*w-w0*w0).^2+(2.*gamma.*w).^2);
plot(w,A2);
%grid on;
xlabel('Vinkelfrekvens (1/s)');
ylabel('Kvadrert amplitude');
