[b,a] = zp2tf(0.1,[-1 -30 -500],1);
figure(1);
bode(b,a);
grid;
G = tf(b,a);
omeg = 1;
S = tf(omeg,[1, 0, omeg^2]) % laplace d'un sinus sin(omeg t)
figure(2);
clf
[out2,t2]=impulse(S*G);
plot(t2(1:800),out2(1:800));
hold on;
gain = 4.73*10^(-5); %valeur maximum dans la figure 2
20*log10(gain)
[out1,t1]=impulse(gain*S); % j'avais pas mis les bons temps t1 et t2 qui sont retournés par impulse dans le zoom live
plot(t1(1:200),out1(1:200));
% on remarque le dephasage qui correspond à la partie angulaire du
% diagramme de Bode. Le temps entre deux pics multiplié par omeg donne
% l'angle. Grosso modo 45 degres donc 45*3.14/180 = 0.78;
% temps associé au déphasage omeg*t, et comme omeg = 1, on a environ 0,78
% secondes, c'est ce que l'on remarque entre les maximums des deux courbes
% de la figure 