% accéleromètre + filtre RC
close all
f = 6000 % Hz
zeta = 0.01  % amortissement

% filtre
R = 100e3
C = 4.55e-9

w = 2*pi*f
TF_acc = tf (1, [1/w^2  2*zeta/w  1])
figure; bode (TF_acc)

TF_RC = tf (1, [R*C 1])

TF = TF_acc * TF_RC
figure; bode (TF)

f_3db = 10^2.54

% réponse indicielle
figure; step (TF_acc); 
title ('Accéleromètre seul')

figure; step (TF); 
title ('Accéleromètre + filtre')

[u,t] = gensig ('sin', 1/f_3db, 0.01, 10e-6)
figure;   lsim( TF, u,t)