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')
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')
El umbral de la frecuencia resalta aún más las regiones de frecuencia común.