xspectrogram
使用短时傅里叶变换的交叉谱图
句法
描述
[___] = xspectrogram（___，返回短期互功率谱密度估计，如果spectrumtype）spectrumtype被指定为'PSD'返回短期互相功率谱估计，如果spectrumtype被指定为'功率'。
xspectrogram（___）在没有输出参数的情况下，在当前图窗口中绘制交叉光谱图。
xspectrogram（___，指定要在其上绘制频率轴。指定freqloc）freqloc因为无论是“x轴”要么“Y轴”。
例子
线性啁啾的交叉谱图
产生在1MHz取样用于10毫秒的两个线性调频脉冲。
第一线性调频脉冲具有150千赫的初始频率由测量结束增加到350千赫。
第二线性调频脉冲具有200千赫的初始频率，通过测量的端增加至300kHz。
加性高斯白噪声，使得信噪比为40分贝。
nSamp = 10000;FS = 1000e3;SNR = 40;T =（0：nSamp-1）'/ FS;X1 =啁啾（T，150e3，T（端部），350e3）;X1 = X1 + randn（尺寸（X1））* STD（X1）/ db2mag（SNR）;X2 =啁啾（T，200e3，T（端部），300e3）;X2 = X2 + randn（大小（×2））* STD（×2）/ db2mag（SNR）;
计算和情节两个啁啾的交叉频谱。划分信号转换成200-样品段和窗口中的每个段具有汉明窗口。指定邻接段和1024个样本的DFT长度之间的80个样本的重叠。
xspectrogram（X1，X2，汉明（200），80,1024，FS，“Y轴”）
修改第二啁啾使频率在测量期间从50千赫至350千赫上升。使用具有形状因子500样品Kaiser窗 到窗口的段。指定的重叠和256计算的DFT长度的450个样本并画出交叉谱图。
X2 =啁啾（T，50e3，T（端部），350e3）;X2 = X2 + randn（大小（×2））* STD（×2）/ db2mag（SNR）;xspectrogram（X1，X2，凯瑟（500,5），450256，FS，“Y轴”）
在这两种情况下，该函数突出了频率内容，这两个信号具有共同点。
语音信号的互谱图
加载包含在44100赫兹采样的两个语音信号的文件。
第一个信号是女声的录音说：“变换函数”。
第二个信号是同一个女人说“改革司法”的录音。
画出两个信号。偏移所述第二信号，以便上下两个是可见的。
加载（'voice.mat'）％听到，类型soundsc（变换，FS），暂停（2），soundsc（改革，FS）T =（0：长度（改革）-1）/ FS;情节（T，变换，吨，改革+ 0.3）图例（“‘变换函数’”，“‘改革司法’”）
计算这两个信号的横谱图。划分信号转换成1000个采样分段和窗口它们与汉明窗口。指定邻接节段之间的重叠800个样品。仅包含频率最高为4千赫。
nwin = 1000;NVLP = 800;FINT = 0：4000;[S，F，T] = xspectrogram（变换，改革，汉明（nwin），NVLP，FINT，FS）;目（T，F，20 *日志10（S））的视图（2）轴紧
交叉频谱图突出了时间间隔，其中信号具有更多的共同频率内容。音节“形式”是格外引人注目。
相移两二次啁啾之间
产生两个二次啁啾，每个采样在1千赫为2秒。两种啁啾的初始频率都是100hz，在测量过程中增加到200hz。与第一个相比，第二个啁啾的相位差为23°。
FS = 1E3;t = 0时：1 / FS：2;Y1 =啁啾（T，100,1,200，“二次”,0);y2 =唧唧声(t, 100, 1200,“二次”、23);
计算啁啾的复杂交叉谱图，提取它们之间的相移。将信号分成128个样本段。指定相邻段之间重叠的120个样本。使用带有形状因子的Kaiser窗口β= 18，并指定128个样本的DFT长度。使用的绘图功能xspectrogram来显示交叉光谱图。
[〜，F，XT，C] = xspectrogram（Y1，Y2，凯瑟（128,18），120128，FS）;xspectrogram（Y1，Y2，凯瑟（128,18），120128，FS，“Y轴”）
提取和显示交叉频谱的最大能量的时间 - 频率脊。
[TFR，〜，lridge] = tfridge（C，F）;保持上情节(xt,总和生育率,“k”，“线宽”2）保持离
该相移是虚部的沿脊随时间变化的交叉谱的实部的比率。计算的相移，并表示它的度。显示其平均值。
pshft =角(c (lridge)) * 180 /π;意思是(pshft)
ANS = -23.0000
复杂信号的互谱图
生成两个信号，每个信号以3 kHz采样1秒。第一个信号是二次啁啾，其频率在测量过程中从300hz增加到1300hz。啁啾信号被嵌入到高斯白噪声中。第二个信号也包含在白噪声中，它是一个啁啾，其频率内容呈正弦变化。
FS = 3000;t = 0时：1 / FS：1-1 / FS;X1 =啁啾（吨，300，T（端部），1300，“二次”）+ randn（大小（T））/ 100;X2 = EXP（2J * PI * 100 * COS（2 * PI * 2 * T））+ randn（大小（T））/ 100;
计算和情节两个信号的交叉频谱。划分信号与邻接节段之间的重叠255个样品256样品段。使用具有形状因子Kaiser窗β= 30〜窗口的段。使用DFT点的默认号码。中心的十字，频谱在零频率。
nwin = 256;xspectrogram（X1，X2，凯瑟（nwin，30），nwin-1，[]，FS，“中心”的，“Y轴”）
计算功率谱，而不是功率谱密度。设置为零值小于-40 dB的小。中心奈奎斯特频率的情节。
xspectrogram（X1，X2，凯瑟（nwin，30），nwin-1，[]，FS，...'功率'，“MinThreshold”，-40，“Y轴”）标题（“二次啁啾的交叉谱图和复杂的啁啾声”）
阈值处理进一步凸显公共频率的区域。
两个序列的交叉谱图
计算和绘图两个序列的交叉频谱。
指定每个序列为4096个样本长。
N = 4096;
为了创建第一序列，生成一个凸二次啁啾嵌入在白高斯噪声和带通滤波器它。
该线性调频脉冲具有0.1π的初始标准化频率由测量结束增加到0.8π。
第16阶滤波器传递0.2π和0.4π弧度/样品之间归一化频率，并具有60dB的阻带衰减。
RX =啁啾（0：N-1,0.1 / 2，N，0.8 / 2，“二次”，[]“凸”)+ randn (N - 1) / 100;dx = designfilt ('bandpassiir'，'FilterOrder'16，...“StopbandFrequency1”，0.2，“StopbandFrequency2”，0.4％，...'StopbandAttenuation'，60）;X =过滤器（DX，RX）;
以创建所述第二序列，生成嵌入在白高斯噪声和带阻滤波器其线性啁啾。
该线性调频脉冲具有0.9π的初始标准化频率其降低由测量结束到0.1π。
第16阶滤波器停止0.6π和0.8π弧度/样品之间归一化频率，并且具有1dB的通带波纹。
RY =啁啾（0：N-1,0.9 / 2，N，0.1 / 2）'+ randn（N，1）/ 100;DY = designfilt（“bandstopiir”，'FilterOrder'16，...'PassbandFrequency1'，0.6，'PassbandFrequency2'，0.8％，...'PassbandRipple'，1）;Y =滤波器（DY，RY）;
画出两个序列。偏移第二序列垂直因此两者是可见的。
图（[X Y + 2]）
计算和绘制的交叉光谱图X和ÿ。使用512样本汉明窗。指定邻接区段和2048 DFT点之间的重叠的500个样本。
xspectrogram（X，Y，汉明（512），500,2048，“Y轴”）
设为零交叉频谱值小于-50dB小。
xspectrogram（X，Y，汉明（512），500,2048，“MinThreshold”，-50，“Y轴”）
谱图表明，由滤波器增强或抑制的频率区域。
输入参数
输出参数
参考文献
[1]米特拉，Sanjit K.数字信号处理：基于计算机A方法。第二版。纽约：麦格劳 - 希尔，2001。
[2]奥本海姆，艾伦五，罗纳德·W·谢弗，和John R.巴克。离散时间信号处理。第二版。上马鞍河，NJ：Prentice Hall出版社，1999年。
扩展功能
介绍了在R2017a
您还可以选择从下面的列表中的网站：
