% per confrontare i diagrammi di Bode (moduli) dei tre sistemi 
% che compaiono nello schema Simulink LegamiTransitorioW

s=tf('s')

sys1=1/(s^2+s+1)

sys2=1/(s^2+0.2*s+1)

sys3=9/(s^2+3*s+9)

w=logspace(-1,1,1000);

[mag1,phase1]=bode(sys1,w);
[mag2,phase2]=bode(sys2,w);
[mag3,phase3]=bode(sys3,w);

semilogx(w,20*log10(mag1(1,:)),'b*',w,20*log10(mag2(1,:)),'r*',w,...
    20*log10(mag3(1,:)),'g*'),grid
axis([0.1,10,-10,15])
xlabel('pulsazione (rad/sec)')
ylabel('modulo (db)')
title('blu=1; rosso=2; verde=3')
