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

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

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