% Exercice page 18
% Simulation d'un accéléromètre à poutre
close all

w = 1e4
zeta = 0.01

tf_acc = tf (1, [1/w^2  2*zeta/w  1])
figure
bode (tf_acc)
grid

w_vec = [0:10:20e3];  % on défini une plage de fréquences
[mag,phase] = bode (tf_acc, w_vec);
i = 1;
while (mag(:,:,i) < 1.05)  % voir le format de mag >> whos mag
    i = i+1;
end
freq_05 = w_vec(i-1)/2/pi  % on trouve la fréquence recherchée
    
figure; step (tf_acc)

[u,t] = gensig ('sin', 1/350, 0.02, 1/100000);
figure; lsim(tf_acc, u, t)
grid on