Esta página aún no se ha traducido para esta versión. Puede ver la versión más reciente de esta página en inglés.

Welch Spectrum Estimates

Cree una señal que consista en tres sinusoides ruidosos y un chirrido, muestreado a 200 kHz durante 0,1 segundos. Las frecuencias de los sinusoides son 1 kHz, 10 kHz y 20 kHz. Los sinusoides tienen diferentes amplitudes y niveles de ruido. El chirrido silencioso tiene una frecuencia que comienza en 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);

Calcular la estimación de Welch PSD y los espectros de retención máxima y de retención mínima de la señal. Trazar 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 estimaciones de espectro de potencia centradas.

[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

Consulte también

| |