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.

Mida el poder de una señal

El poder de una señal es la suma de los cuadrados absolutos de sus muestras de dominio del tiempo dividido por la longitud de la señal, o, equivalentemente, el cuadrado de su nivel RMS. La función le permite estimar la potencia de la señal en un solo paso.bandpower

Considere una unidad Chirp incrustada en ruido Gaussiano blanco y muestreada a 1 kHz durante 1,2 segundos. La frecuencia del Chirp aumenta en un segundo de un valor inicial de 100 Hz a 300 Hz. El ruido tiene varianza

<math display="block">
<mrow>
<mn>0</mn>
<mo>.</mo>
<mn>0</mn>
<msup>
<mrow>
<mn>1</mn>
</mrow>
<mrow>
<mn>2</mn>
</mrow>
</msup>
</mrow>
</math>
. Restablezca el generador de números aleatorios para obtener resultados reproducibles.

N = 1200; Fs = 1000; t = (0:N-1)/Fs;  sigma = 0.01; rng('default')  s = chirp(t,100,1,300)+sigma*randn(size(t));

Compruebe que la estimación de potencia dada por es equivalente a la definición.bandpower

pRMS = rms(s)^2
pRMS = 0.5003 
powbp = bandpower(s,Fs,[0 Fs/2])
powbp = 0.5005 

Utilice la función para estimar el ancho de la banda de frecuencias que contiene el 99% de la potencia de la señal, los límites inferior y superior de la banda, y la potencia en la banda.obw La función también traza la estimación del espectro y anota el ancho de banda ocupado.

obw(s,Fs);

 [wd,lo,hi,power] = obw(s,Fs); powtot = power/0.99
powtot = 0.5003 

A un amplificador de potencia no lineal se le da una sinusoide de 60 Hz como entrada y emite una señal ruidosa con distorsión de tercer orden. La frecuencia de muestreo es de 3,6 kHz. Reste el componente de frecuencia cero (DC) para concentrarse en el contenido espectral.

load(fullfile(matlabroot,'examples','signal','AmpOutput.mat')) Fs = 3600; y = y-mean(y);

Debido a que el amplificador presenta distorsión de tercer orden, se espera que la señal de salida tenga

  • Un componente con la misma frecuencia que la entrada, 60 Hz;fundamental

  • Dos--componentes de frecuencia en dos y tres veces la frecuencia de la entrada, 120 y 180 Hz.harmonics

Se utiliza para determinar la potencia almacenada en el fundamental y los armónicos.bandpower Exprese cada valor como un porcentaje de la potencia total y en decibelios. Mostrar los valores como una tabla.

pwrTot = bandpower(y,Fs,[0 Fs/2]);  Harmonic = {'Fundamental';'First';'Second'};  Freqs = [60 120 180]';  Power = zeros([3 1]); for k = 1:3     Power(k) = bandpower(y,Fs,Freqs(k)+[-10 10]); end  Percent = Power/pwrTot*100;  inDB = pow2db(Power);  T = table(Freqs,Power,Percent,inDB,'RowNames',Harmonic)
T=3×4 table
                   Freqs     Power      Percent      inDB  
                   _____    ________    _______    ________

    Fundamental      60       1.0079    12.563     0.034136
    First           120      0.14483    1.8053      -8.3914
    Second          180     0.090023    1.1221      -10.456

Consulte también

| | |

Temas relacionados