diff --git a/src/comparaison_rep_freq.m b/src/comparaison_rep_freq.m new file mode 100644 index 0000000000000000000000000000000000000000..6b4270df22fb3a247fc87e734d70f8589c2fdbaf --- /dev/null +++ b/src/comparaison_rep_freq.m @@ -0,0 +1,33 @@ +clear; clc; + +% Paramètres +tau = 100; +g = 0.6; +N = 2048; + +nu = linspace(0, 0.1, N/10); + +% Réponse théorique +H_th = abs(1 ./ (1 + g * exp(-1j * 2 * pi * nu * tau))); +H_th_db = 20 * log10(H_th); + +% Réponse impulsionnelle numérique (on a recopié analyse_delay.m) +a = [1, zeros(1, tau - 1), g]; +b = 1; +x = zeros(1, N); +x(1) = 1; +h = filter(b, a, x); + +% FFT de h et conversion en dB +H_fft = abs(fft(h, N)); +H_fft_db = 20 * log10(H_fft(1:N/10) + eps); + +% Tracé +figure; +plot(nu, H_th_db, 'r-', 'LineWidth', 1.5); hold on; +plot(nu, H_fft_db, 'b--'); +xlabel('Fréquence réduite \nu'); +ylabel('|H(\nu)| en dB'); +title('Réponse en fréquence (module) sur plage réduite [0, 0.1]'); +legend('Théorique', 'Numérique'); +grid on;