Main Content

Espectrograma cruzado de señales complejas

Genere dos señales, cada una de ellas muestreada a 3 kHz durante 1 segundo. La primera señal es un chirp cuadrático cuya frecuencia aumenta desde 300 Hz hasta 1300 Hz durante la medición. El chirp está integrado en un ruido blanco gaussiano. La segunda señal, también integrada en ruido blanco, es un chirp con contenido de frecuencia que varía sinusoidalmente.

fs = 3000;
t = 0:1/fs:1-1/fs;

x1 = chirp(t,300,t(end),1300,'quadratic')+randn(size(t))/100;

x2 = exp(2j*pi*100*cos(2*pi*2*t))+randn(size(t))/100;

Calcule y represente el espectrograma cruzado de las dos señales. Divida las señales en segmentos de 256 muestras con 255 muestras de solapamiento entre segmentos contiguos. Utilice una ventana de Kaiser con un factor de forma β = 30 para acotar los segmentos. Utilice el número de puntos DFT predeterminado. Centre el espectrograma cruzado en la frecuencia cero.

nwin = 256;

xspectrogram(x1,x2,kaiser(nwin,30),nwin-1,[],fs,'centered','yaxis')

Figure contains an axes object. The axes object with xlabel Time (ms), ylabel Frequency (kHz) contains an object of type image.

Calcule el espectro de potencia en lugar de la densidad espectral de potencia. Ponga a cero los valores inferiores a -40 dB. Centre la gráfica en la frecuencia de Nyquist.

xspectrogram(x1,x2,kaiser(nwin,30),nwin-1,[],fs, ...
    'power','MinThreshold',-40,'yaxis')
title('Cross-Spectrogram of Quadratic Chirp and Complex Chirp')

Figure contains an axes object. The axes object with title Cross-Spectrogram of Quadratic Chirp and Complex Chirp, xlabel Time (ms), ylabel Frequency (kHz) contains an object of type image.

El umbral de la frecuencia resalta aún más las regiones de frecuencia común.

Consulte también

|

Temas relacionados