Перекрестная спектрограмма комплексных сигналов

Сгенерируйте два сигнала, каждый из которых дискретизирован на частоте 3 кГц в течение 1 секунды. Первый сигнал является квадратичным щебетом, частота которого увеличивается с 300 Гц до 1300 Гц во время измерения. Щебет встроен в белый Гауссов шум. Второй сигнал, также встроенный в белый шум, является щебетом с синусоидально изменяющимся содержимым частоты.

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;

Вычислите и постройте график перекрестной спектрограммы двух сигналов. Разделите сигналы на 256 сегментов с 255 выборками перекрытия между смежными сегментами. Используйте окно Кайзера с масштабным фактором β = 30, чтобы открыть сегменты. Используйте количество точек ДПФ по умолчанию. Центрируйте поперечную спектрограмму на нулевой частоте.

nwin = 256;

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

Вычислите спектр степени вместо степени спектральной плотности. Установите таким образом, чтобы нуль значения меньше -40 дБ. Центрируйте график на частоте Найквиста.

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

Пороговое значение дополнительно выделяет области общей частоты.

См. также

spectrogram | xspectrogram