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