Main Content

Medir frecuencia media, potencia, ancho de banda

Genere 1024 muestras de un chirp muestreado a 1024 kHz. El chirp presenta una frecuencia inicial de 50 kHz y alcanza los 100 kHz al finalizar el muestreo. Añada ruido blanco gaussiano hasta que la relación señal-ruido sea de 40 dB.

nSamp = 1024;
Fs = 1024e3;
SNR = 40;

t = (0:nSamp-1)'/Fs;

x = chirp(t,50e3,nSamp/Fs,100e3);
x = x+randn(size(x))*std(x)/db2mag(SNR);

Calcule el 99% del ancho de banda ocupado de la señal y anótelo en una gráfica de densidad del espectro de potencia (PSD).

obw(x,Fs);

Figure contains an axes object. The axes object with title 99% Occupied Bandwidth: 55.377 kHz contains 4 objects of type line, patch.

Calcule la potencia en la banda y verifique que es el 99% del total.

[bw,flo,fhi,powr] = obw(x,Fs);

pcent = powr/bandpower(x)*100
pcent = 99.0000

Genere otro chirp. Especifique una frecuencia inicial de 200 kHz, una frecuencia final de 300 kHz y una amplitud que sea el doble que la de la primera señal. Añada ruido blanco gaussiano.

x2 = 2*chirp(t,200e3,nSamp/Fs,300e3);
x2 = x2+randn(size(x2))*std(x2)/db2mag(SNR);

Añada los dos chirps para formar una nueva señal. Represente la PSD de la señal y anote su frecuencia mediana.

medfreq([x+x2],Fs);

Figure contains an axes object. The axes object with title Median Frequency Estimate: 237.563 kHz contains 2 objects of type line.

Represente la PSD y anote la frecuencia media.

meanfreq([x+x2],Fs);

Figure contains an axes object. The axes object with title Mean Frequency Estimate: 214.962 kHz contains 2 objects of type line.

Ahora, piense que cada chirp representa un canal independiente. Calcule la frecuencia media de cada canal. Anote las frecuencias medias en una gráfica de las PSD.

meanfreq([x x2],Fs)

Figure contains an axes object. The axes object with title Mean Frequency Estimate contains 4 objects of type line.

ans = 1×2
105 ×

    0.7503    2.4999

Calcule el ancho de banda de media potencia de cada canal. Anote los anchos de banda de 3 dB en una gráfica de las PSD.

powerbw([x x2],Fs)

Figure contains an axes object. The axes object with title 3 dB Bandwidth contains 8 objects of type line, patch.

ans = 1×2
104 ×

    4.4386    9.2208

Consulte también

| | | |