% Paramètres du système
m = 1; % Masse (kg)
wn = 2*pi*1; % Fréquence naturelle (rad/s), exemple pour 1 Hz
zeta = 0.2; % Coefficient d'amortissement (sous-amorti: 0 < zeta < 1)
% Calcul de la fréquence amortie
wd = wn*sqrt(1 - zeta^2);
% Définition du temps
t = linspace(0, 5, 1000); % De 0 à 5 secondes avec 1000 points
% Calcul de la réponse impulsionnelle
h = (exp(-zeta*wn*t) ./ (m*wd)) .* sin(wd*t);
% Tracé de la réponse
figure;
plot(t, h, 'b', 'LineWidth', 1.5);
grid on;
xlabel('Temps (s)');
ylabel('Amplitude h(t)');
title('Réponse impulsionnelle d''un système du second ordre sous-amorti');
legend(['ζ = ' num2str(zeta) ', ω_n = ' num2str(wn) ' rad/s']);
% Ajout des lignes d'enveloppe exponentielle si souhaité
hold on;
env_pos = exp(-zeta*wn*t) ./ (m*wd);
env_neg = -env_pos;
plot(t, env_pos, 'r--', t, env_neg, 'r--');
hold off;