生成两个信号,每个信号在3 kHz上采样1秒。第一信号是二次啁啾,其在测量期间从300 Hz到1300Hz增加。啁啾嵌入白高斯噪音中。第二信号,也嵌入白噪声,是具有正弦变化频率内容的啁啾。
FS = 3000;T = 0:1 / FS:1-1 / FS;X1 = Chirp(T,300,T(端),1300,'二次')+ Randn(尺寸(t))/ 100;X2 = EXP(2J * PI * 100 * COS(2 * PI * 2 * T))+ RANDN(尺寸(t))/ 100;
计算并绘制两个信号的跨谱图。将信号分成256个样本段,在相邻的区段之间具有255个重叠样本。使用带状因子的kaiser窗口β= 30窗口段。使用默认的DFT点数。以零频率为中心跨谱图。
nwin = 256;XspectRogre(X1,X2,Kaiser(NWIN,30),NWIN-1,[],FS,'中心'那'yaxis')
计算功率谱而不是功率谱密度。设置为零,值小于-40 dB。以奈奎斯特频率为中心。
XspectRogre(X1,X2,Kaiser(NWIN,30),NWIN-1,[],FS,......'力量'那'minthreshold',-40,'yaxis') 标题(“二次啁啾和复杂啁啾的跨谱图”)
阈值化进一步突出了公共频率的区域。