pwelch
Estimación de densidad espectral de potencia de Welch
Sintaxis
Descripción
devuelve la estimación de densidad espectral de potencia (PSD), pxx
= pwelch(x
)pxx
, de la señal de entrada, x
, obtenida utilizando el estimador de promediado de segmentos superpuestos de Welch. Cuando x
es un vector, se trata como un solo canal. Cuando x
es una matriz, la PSD se calcula de forma independiente en cada columna y se almacena en la columna correspondiente de pxx
. Si x
tiene valor real, pxx
es una estimación de PSD unilateral. Si x
tiene valor complejo, pxx
es una estimación de PSD bilateral. De forma predeterminada, x
se divide en los segmentos más largos posibles para acercarse lo máximo posible, sin sobrepasar, a ocho segmentos con un 50% de solapamiento. Cada segmento se dispone en ventanas con una ventana Hamming. Los periodogramas modificados se promedian para obtener la estimación de PSD. Si no puede dividir la longitud de x
exactamente en un número entero de segmentos con 50% de solapamiento, x
se trunca en consecuencia.
usa el vector de entrada o entero, pxx
= pwelch(x
,window
)window
, para dividir la señal en segmentos. Si window
es un vector, pwelch
divide la señal en segmentos de la misma longitud que la longitud de window
. Los periodogramas modificados se calculan mediante los segmentos de señal multiplicados por el vector, window
. Si window
es un entero, la señal se divide en segmentos de longitud window
. Los periodogramas modificados se calculan mediante una ventana de Hamming de longitud window
.
utiliza pxx
= pwelch(x
,window
,noverlap
) noverlap
muestras de solapamiento de segmento a segmento. noverlap
debe ser un entero positivo menor que window
si window
es un entero. noverlap
debe ser un entero positivo menor que la longitud de window
si window
es un vector. Si no especifica noverlap
, o especifica noverlap
como vacío, el número predeterminado de muestras solapadas es el 50% de la longitud de la ventana.
[
devuelve un vector de frecuencia, pxx
,f
] = pwelch(___,fs
)f
, en ciclos por unidad de tiempo. La tasa de muestreo, fs
, es el número de muestras por unidad de tiempo. Si la unidad de tiempo es el segundo, f
estará en ciclos/segundo (Hz). Para señales de valor real, f
abarca el intervalo [0,fs
/2] cuando nfft
es par y [0,fs
/2) cuando nfft
es impar. Para señales de valor complejo, f
abarca el intervalo [0,fs
). fs
debe ser la quinta entrada a pwelch
. Para introducir una tasa de muestreo y seguir utilizando los valores predeterminados de los argumentos opcionales anteriores, especifique estos argumentos como vacíos, []
.
[
devuelve las estimaciones de PSD de Welch bilaterales en las frecuencias especificadas en el vector, pxx
,f
] = pwelch(x
,window
,noverlap
,f
,fs
)f
. El vector f
debe contener al menos dos elementos porque, de lo contrario, la función lo interpreta como nfft
. Las frecuencias de f
están en ciclos por unidad de tiempo. La tasa de muestreo, fs
, es el número de muestras por unidad de tiempo. Si la unidad de tiempo es el segundo, f
estará en ciclos/segundo (Hz).
[___,
devuelve los intervalos de confianza de pxxc
] = pwelch(___,'ConfidenceLevel',probability
)probability
× 100% para la estimación de PSD en pxxc
.
[___] = pwelch(___,
devuelve la estimación de PSD si spectrumtype
)spectrumtype
está especificado como 'psd'
y devuelve el espectro de potencia si spectrumtype
está especificado como 'power'
.
pwelch(___)
sin argumentos de salida representa la estimación de PSD de Welch en la ventana de figura actual.
Ejemplos
Argumentos de entrada
Argumentos de salida
Más acerca de
Referencias
[1] Hayes, Monson H. Statistical Digital Signal Processing and Modeling. New York: John Wiley & Sons, 1996.
[2] Stoica, Petre, and Randolph Moses. Spectral Analysis of Signals. Upper Saddle River, NJ: Prentice Hall, 2005.