%Ćwiczenie nr 2
%Korzystając ze wzorów na amplitudę drgań tłumionych narysuj wykres
%zależności amplitudy od częstości siły wymuszającej dla 3 różnych
%współczynników tłumienia b dla podanych danych.

%F0=1
%b0=0.1
%b1=1
%b2=10
%wwym: 0:0.05:100
%LpKM=3
%MASA (m=290[g]=0,290[kg]
%Wsp SPREŻYSTOŚCI: (k=98[N/m])
k=98;
f0=1;
m=0.290;
b0=0.1;
b1=1;
b2=10;
%w0=k/m=98/0.290;
w0=sqrt(k/m);
wwym=0:0.05:100;
w1=(w0^2-b0.^2)^(1/2);
w2=(w0^2-b1.^2)^(1/2);
w3=(w0^2-b2.^2)^(1/2);

A1=(f0/m).*(1./((wwym.^2-w1.^2).^2+4.*(b0.^2).*(w1.^2)).^(1/2));
A2=(f0/m).*(1./((wwym.^2-w2.^2).^2+4.*(b1.^2).*(w2.^2)).^(1/2));
A3=(f0/m).*(1./((wwym.^2-w3.^2).^2+4.*(b2.^2).*(w3.^2)).^(1/2));

hold;

plot(wwym,A1);
plot(wwym,A2);
plot(wwym,A3);
title('A zależne od częstości siły wymuszającej');
xlabel('omega[rad/s]');
ylabel('Amplituda[m]');
%====
legend('b0','b1','b2')