close all
clear all
G=1;
w=logspace(-3,3);
w0=0.1*2*pi;
K=0.0039;
m=0.1;
D=1;

x=(G*w.^2)./(sqrt((w.^2-w0^2).^2+4*w.^2*w0.^2*1));
X=real(x);
f=w./(2*pi);
figure(1)
loglog(f,X)
%%%
G2=1;
w2=logspace(-3,3);
w02=0.19748;
K=0.0039;
m2=0.01;
D=1;

x2=(G2*w2.^2)./(sqrt((w2.^2-w02^2).^2+4*w2.^2*w02.^2*1));
X2=real(x2);
f2=w2./(2*pi);
figure(2)
loglog(f2,X2)


figure(3)
loglog(f,X,'r')
hold on
loglog(f2,X2,'b')
hold on

ylim([0.01 2])
xlim([0.001 10])
legend('#1','#2')
title('Amplitude response')
xlabel('Linear frequency - f')

