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
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