Upload New File

src/wav/question_3_4.m

+
+% Étape 0 : charger les signaux et fréquence d'échantillonnage
+load("signal_excitation.mat");  % contient xe1, xe2, fe
+
+% Étape 1 : Choisir le signal d'excitation
+x = xe1;              % tu peux tester aussi avec xe2 si tu veux comparer
+N = length(x);
+% Étape 2 : Simuler la sortie de la pièce
+y = simule_piece(x,fe);  % réponse mesurée simulée
+
+% Étape 3 : Estimer la réponse impulsionnelle via la corrélation croisée
+R_yx = xcorr(y, x, 'none');  % estimation de h
+L = length(R_yx);
+lags = -floor((L-1)/2):ceil((L-1)/2);
+t = lags / fe;                 % axe temporel en secondes
+
+% Étape 4 : Afficher la réponse impulsionnelle estimée
+figure;
+plot(t, R_yx);
+xlabel("Temps (s)");
+ylabel("Amplitude");
+title("Réponse impulsionnelle estimée de la pièce");
+grid on;
\ No newline at end of file