Calcul du Périodogramme
Le programme MATLAB suivant permet de calculer le périodogramme d'un signal x généré par la fonction sig_bruit ainsi développée.
Simulation : Script MATLAB N_1
clear all; close all;
N = 1024; % nombre d'échantillons
% Générer le signal à analyser en utilisant sig_bruit soit 250 Hz sin plus bruit blanc; N échantillons [x,t] =sig_bruit(250,N);
fs = 1000; % la fréquence d'échantillonnage est de 1 kHz Y = fft(x); %Calculer la TFR
PS = abs(Y).ˆ2; % calculer PS comme magnitude au carré
freq = (1:N)/fs; % vecteur fréquence à utiliser avec le tracé plot(freq,20*log10(PS)); tracer PS en échelle logarithmique title('spectre en puissance (symétrique de fs/2)'); xlabel('Fréquence (Hz)');
ylabel('spectre en puissance (dB)');
