TP 2 Estimation spectrale non-paramétrique

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.

SimulationScript 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)');

Périodogramme
PrécédentPrécédentSuivantSuivant
AccueilAccueilImprimerImprimerRéalisé avec Scenari (nouvelle fenêtre)