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.

bandpower

Descripción

ejemplo

p = bandpower(x) devuelve la potencia media en la señal de entrada, .x Si es una matriz, calcula la potencia media en cada columna de forma independiente.xbandpower

ejemplo

p = bandpower(x,fs,freqrange) devuelve la potencia media en el rango de frecuencia, , especificada como un vector de dos elementos.freqrange Debe introducir la frecuencia de muestreo, , para devolver la potencia en un intervalo de frecuencia especificado. utiliza un periodograma modificado para determinar la potencia media en .fsbandpowerfreqrange

ejemplo

p = bandpower(pxx,f,'psd') devuelve la potencia media calculada mediante la integración de la estimación de densidad espectral de potencia (PSD), .pxx La integral se aproxima mediante el método rectángulo. La entrada, , es un vector de frecuencias correspondiente a las estimaciones PSD en .fpxx La opción indica que la entrada es una estimación PSD y no datos de series temporales.'psd'

ejemplo

p = bandpower(pxx,f,freqrange,'psd') devuelve la potencia media contenida en el intervalo de frecuencia, .freqrange Si las frecuencias en no coinciden con los valores en , se utilizan los valores más cercanos.freqrangef La potencia media se calcula integrando la estimación de densidad espectral de potencia (PSD), .pxx La integral se aproxima mediante el método rectángulo. La opción indica que la entrada es una estimación PSD y no datos de series temporales.'psd'

Ejemplos

contraer todo

Cree una señal que consista en una onda sinusoidal de 100 Hz en ruido gaussiano blanco (0,1).N La frecuencia de muestreo es de 1 kHz. Determinar la potencia media y compararla con la

<math display="block">
<mrow>
<msub>
<mrow>
<mo></mo>
</mrow>
<mrow>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</math>
Norma.

t = 0:0.001:1-0.001; x = cos(2*pi*100*t)+randn(size(t));  p = bandpower(x)
p = 1.5264 
l2norm = norm(x,2)^2/numel(x)
l2norm = 1.5264 

Determine el porcentaje de la potencia total en un intervalo de frecuencia especificado.

Cree una señal que consista en una onda sinusoidal de 100 Hz en ruido gaussiano blanco (0,1).N La frecuencia de muestreo es de 1 kHz. Determine el porcentaje de potencia total en el intervalo de frecuencia entre 50 Hz y 150 Hz. Restablezca el generador de números aleatorios para obtener resultados reproducibles.

rng('default')  t = 0:0.001:1-0.001; x = cos(2*pi*100*t)+randn(size(t));  pband = bandpower(x,1000,[50 150]); ptot = bandpower(x,1000,[0 500]); per_power = 100*(pband/ptot)
per_power = 51.9591 

Determine la potencia media calculando primero una estimación de PSD utilizando el periodograma. Introduzca la estimación PSD en .bandpower

Cree una señal que consista en una onda sinusoidal de 100 Hz en ruido gaussiano blanco (0,1).N La frecuencia de muestreo es de 1 kHz. Obtenga el periodograma y utilice la bandera para calcular la potencia media utilizando la estimación PSD.'psd' Compare el resultado con la potencia media calculada en el dominio de tiempo.

t = 0:0.001:1-0.001; Fs = 1000; x = cos(2*pi*100*t)+randn(size(t));  [Pxx,F] = periodogram(x,rectwin(length(x)),length(x),Fs); p = bandpower(Pxx,F,'psd')
p = 1.5264 
avpow = norm(x,2)^2/numel(x)
avpow = 1.5264 

Determine el porcentaje de la potencia total en un intervalo de frecuencia especificado utilizando el periodograma como entrada.

Cree una señal que consista en una onda sinusoidal de 100 Hz en ruido gaussiano blanco (0,1).N La frecuencia de muestreo es de 1 kHz. Obtenga el periodograma y el vector de frecuencia correspondiente. Utilizando la estimación PSD, determine el porcentaje de la potencia total en el intervalo de frecuencia entre 50 Hz y 150 Hz.

Fs = 1000; t = 0:1/Fs:1-0.001; x = cos(2*pi*100*t)+randn(size(t));  [Pxx,F] = periodogram(x,rectwin(length(x)),length(x),Fs); pBand = bandpower(Pxx,F,[50 150],'psd'); pTot = bandpower(Pxx,F,'psd'); per_power = 100*(pBand/pTot)
per_power = 49.1798 

Cree una señal multicanal que consista en tres sinusoides en ruido gaussiano blanco (0,1) aditivo.N Las frecuencias de los sinusoidos son 100 Hz, 200 Hz y 300 Hz. La frecuencia de muestreo es de 1 kHz, y la señal tiene una duración de 1 s.

Fs = 1000;  t = 0:1/Fs:1-1/Fs;  f = [100;200;300];  x = cos(2*pi*f*t)'+randn(length(t),3);

Determinar la potencia media de la señal y compararla con la

<math display="block">
<mrow>
<msub>
<mrow>
<mo></mo>
</mrow>
<mrow>
<mn>2</mn>
</mrow>
</msub>
</mrow>
</math>
Norma.

p = bandpower(x)
p = 1×3

    1.5264    1.5382    1.4717

 l2norm = dot(x,x)/length(x)
l2norm = 1×3

    1.5264    1.5382    1.4717

Argumentos de entrada

contraer todo

Introduzca datos de series temporales, especificados como vector de fila o columna o como matriz. Si es una matriz, sus columnas se tratan como canales independientes.x

Ejemplo: es una señal de vector de columna de un solo canal.cos(pi/4*(0:159))'+randn(160,1)

Ejemplo: es un sinusoides ruidoso de dos canales.cos(pi./[4;2]*(0:159))'+randn(160,2)

Tipos de datos: double | single
Soporte de números complejos:

Frecuencia de muestreo para los datos de la serie temporal de entrada, especificada como un escalar positivo.

Tipos de datos: double | single

Rango de frecuencia para el cálculo de potencia de banda, especificado como un vector de fila o columna de valor real de dos elementos. Si la señal de entrada, , contiene muestras, debe estar dentro de los siguientes intervalos:xNfreqrange

  • [0, /2] si es de valor real e inclusofsxN

  • [0, ( – 1) /(2 )] si tiene valor real y es imparNfsNxN

  • [–( – 2) /(2 ), /2] si tiene un valor complejo e inclusoNfsNfsxN

  • [–( – 1) /(2 ), ( – 1) /(2 )] si tiene un valor complejo y es imparNfsNNfsNxN

Tipos de datos: double | single

Estimaciones PSD de una o dos caras, especificadas como un vector o matriz de columna de valor real con elementos no negativos.

La densidad espectral de potencia debe expresarse en unidades lineales, no en decibelios. Se utiliza para convertir valores de decibelios en valores de potencia.db2pow

Ejemplo: especifica la estimación de la DSP de periodograma de un sinusoides de dos canales ruidoso muestreado a 2 oS y las frecuencias a las que se calcula.[pxx,f] = periodogram(cos(pi./[4;2]*(0:159))'+randn(160,2))

Tipos de datos: double | single

Vector de frecuencia, especificado como vector de columna. El vector de frecuencia, , contiene las frecuencias correspondientes a las estimaciones PSD en .fpxx

Tipos de datos: double | single

Argumentos de salida

contraer todo

Poder promedio de la banda, devuelto como un escalar no negativo.

Tipos de datos: double | single

Consulte también

|

Introducido en R2013a