Сгенерируйте два сигнала, каждый произведенный на уровне 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')
Пороговая обработка дальнейшие подсветки области общей частоты.