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