% capteur de force
w = 1e4  % rad/s
zeta = 0.1  % amortissement

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

f_lim = 10^2.54

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