From 4313ccbb4770c66b97d12bb7343a32291476cef7 Mon Sep 17 00:00:00 2001
From: PERBEN Anatole <anatole.perben@imt-atlantique.net>
Date: Wed, 28 May 2025 23:28:52 +0000
Subject: [PATCH] Upload New File

---
 src/wav/question_3_4.m | 23 +++++++++++++++++++++++
 1 file changed, 23 insertions(+)
 create mode 100644 src/wav/question_3_4.m

diff --git a/src/wav/question_3_4.m b/src/wav/question_3_4.m
new file mode 100644
index 0000000..0ce1067
--- /dev/null
+++ b/src/wav/question_3_4.m
@@ -0,0 +1,23 @@
+
+% É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
-- 
GitLab