diff --git a/SAR3Q4.m b/SAR3Q4.m
new file mode 100644
index 0000000000000000000000000000000000000000..56501dc3340c96ed3e04acd9c1713fd40d2dfafd
--- /dev/null
+++ b/SAR3Q4.m
@@ -0,0 +1,19 @@
+load('signal_excitation.mat');  % xe1, xe2, Fe
+x = xe1;                         % On choisit xe1
+x = x ;             % Normalisation
+
+load('empreinte_ac');            % fe
+
+ys = simule_piece(x, fe);        % Convolution avec la vraie empreinte
+
+% Corrélation
+[h_est, lags] = xcorr(ys, x, 'none');
+t = lags / fe; 
+
+% Affichage amplitude réelle
+figure;
+plot(t, h_est);
+xlabel('Temps (s)');
+ylabel('Amplitude (V²)');
+title('Réponse impulsionnelle estimée de la pièce (non normalisée)');
+grid on;