Main Content

Estimaciones del espectro Welch

Cree una señal compuesta por tres sinusoides con ruido y un chirp, muestreados a 200 kHz durante 0,1 segundos. Las frecuencias de las sinusoides son 1 kHz, 10 kHz y 20 kHz. Las sinusoides presentan diferentes amplitudes y niveles de ruido. El chirp sin ruido presenta una frecuencia que comienza a 20 kHz y aumenta linealmente a 30 kHz durante el muestreo.

Fs = 200e3; 
Fc = [1 10 20]'*1e3; 
Ns = 0.1*Fs;

t = (0:Ns-1)/Fs;
x = [1 1/10 10]*sin(2*pi*Fc*t)+[1/200 1/2000 1/20]*randn(3,Ns);
x = x+chirp(t,20e3,t(end),30e3);

Calcule la estimación de PSD de Welch y los espectros de retención máxima y mínima de la señal. Represente los resultados.

[pxx,f] = pwelch(x,[],[],[],Fs);
pmax = pwelch(x,[],[],[],Fs,'maxhold');
pmin = pwelch(x,[],[],[],Fs,'minhold');

plot(f/1000,pow2db(pxx))
hold on
plot(f/1000,pow2db([pmax pmin]),':')
hold off
xlabel('Frequency (kHz)')
ylabel('PSD (dB/Hz)')
legend('pwelch','maxhold','minhold')
grid

Figure contains an axes object. The axes object contains 3 objects of type line. These objects represent pwelch, maxhold, minhold.

Repita el procedimiento, esta vez calculando las estimaciones del espectro de potencia centrado.

[pxx,f] = pwelch(x,[],[],[],Fs,'centered','power');
pmax = pwelch(x,[],[],[],Fs,'maxhold','centered','power');
pmin = pwelch(x,[],[],[],Fs,'minhold','centered','power');

plot(f/1000,pow2db(pxx))
hold on
plot(f/1000,pow2db([pmax pmin]),':')
hold off
xlabel('Frequency (kHz)')
ylabel('Power (dB)')
legend('pwelch','maxhold','minhold')
title('Centered Power Spectrum Estimates')
grid

Figure contains an axes object. The axes object with title Centered Power Spectrum Estimates contains 3 objects of type line. These objects represent pwelch, maxhold, minhold.

Consulte también

| |