xspectrogram
利用短时傅立叶变换交叉频谱
句法
描述
[___] = xSpectRogre（___那返回短期交叉功率谱密度估计Spectrumtype.）Spectrumtype.被指定为'psd'并返回短期交叉功率谱估计Spectrumtype.被指定为“权力”．
xspectrogram (___）不带输出参数绘制在目前的数字窗口中的交叉频谱。
xspectrogram (___那指定绘制频率的轴。指定freqloc）freqloc因为无论是“x轴”或“桠溪”．
例子
线性啁啾的交叉谱图
产生两个线性啁啾采样在1 MHz为10毫秒。
第一啁啾具有150 kHz的初始频率，通过测量结束时增加到350 kHz。
第二个啁啾的初始频率为200 kHz，在测量结束时增加到300 kHz。
加入高斯白噪声，使信噪比为40 dB。
nSamp = 10000;Fs = 1000年e3;信噪比= 40;t = (0: nSamp-1) / Fs;x1 =唧唧声(t, 150年e3, t(结束),350年e3);x1 = x1 + randn(大小(x1)) *性病(x1) / db2mag(信噪比);x2 =唧唧声(t, 200年e3, t(结束),300年e3);x2 = x2 + randn(大小(x2)) *性病(x2) / db2mag(信噪比);
计算并绘制两个啁啾的交叉谱图。将信号分成200个样本段，每个样本段有一个汉明窗口。指定80个相邻片段之间的重叠样本和1024个样本的DFT长度。
xspectrogram（X1，X2，汉明（200），80,1024，FS，“桠溪”）
修改第二啁啾，使频率在测量期间从50 kHz到350 kHz上升。使用带有形状因子的500样品kaiser窗口 到窗口的段。指定的重叠和256计算的DFT长度的450个样本并画出交叉谱图。
x2 =唧唧声(t, 50 e3, t(结束),350年e3);x2 = x2 + randn(大小(x2)) *性病(x2) / db2mag(信噪比);xspectrogram (x1, x2,凯瑟(500 5),450256年,Fs,“桠溪”）
在这两种情况下，该功能突出显示两个信号共同的频率内容。
语音信号的跨谱图
加载包含在44,100 Hz上采样的两个语音信号的文件。
第一信号是拍摄“变换功能”的雌性语音的记录。
第二个信号是相同的女声说的记录“改革正义”。
绘制两个信号。垂直偏移第二信号，因此两者都是可见的。
负载('Voice.Mat'）% To hear, type soundsc(transform,fs)，pause(2)，soundsc(reform,fs)t =（0：长度（改革）-1）/ fs;情节（T，变换，T，改革+ 0.3）传奇（'“变换功能”'那“改革司法”的）
计算两个信号的交叉谱图。将信号分成1000个样本段，并用汉明窗对其加窗。指定800个相邻段重叠的样本。只包括4千赫以下的频率。
nwin = 1000;nvlp = 800;弗林特= 0:4000;f (s, t) = xspectrogram(变换,改革,汉明(nwin) nvlp,弗林特,fs);网格(t, f, 20 * log10 (s))视图(2)轴紧的
跨谱图突出了信号具有更多共同频率内容的时间间隔。音节“form”是特别明显的。
两个二次啁啾之间的相移
生成两个二次啁啾，各在1 kHz取样2秒钟。两个线性调频脉冲具有100Hz的初始频率增大到通过测量200赫兹中途。第二线性调频脉冲具有23的相位差°相比于第一。
fs = 1 e3;t = 0:1 / fs: 2;日元=唧唧声(t, 100, 1200,“二次”，0）;Y2 =啁啾（T，100,1,200，“二次”，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,“桠溪”）
提取并显示跨谱图的最大能量时频脊。
[TFR，〜，Lridge] = TFRIDGE（C，F）;抓住在图（XT，TFR，数k那'行宽'，2）持有从
相移是沿山脊的时变交叉谱的虚部与实部之比。计算相移并以度数表示。显示其平均值。
pshft =夹角（c（lridge））* 180 / PI;平均（pshft）
ans = -23.0000.
复杂信号的交叉谱图
产生两个信号，每个信号在3千赫取样用于1秒。该第一信号是二次啁啾的频率在测量期间从300赫兹到1300赫兹的增加。啁啾是嵌入在高斯白噪声。所述第二信号，也嵌入在白噪声，是用正弦变化的频率成分的线性调频脉冲。
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;xspectrogram (x1, x2, kaiser (nwin 30), nwin-1, [], fs,“中心”那“桠溪”）
计算功率谱而不是功率谱密度。设置为零，值小于-40 dB。以奈奎斯特频率为中心。
xspectrogram (x1, x2, kaiser (nwin 30), nwin-1, [], fs,．..“权力”那'MinThreshold', -40,“桠溪”） 标题（二次啁啾和复啁啾的交叉谱图）
阈值进一步突出了公共频率区域。
两个序列的交叉谱图
计算和绘制两个序列的跨谱图。
指定每个序列为4096个样本长。
n = 4096;
为了创建第一个序列，生成一个嵌入高斯白噪声的凸二次啁啾，并对其进行带通滤波。
啁啾的初始归一化频率为0.1π，在测量结束时增加到0.8π。
16阶滤波器的归一化频率在0.2 ~ 0.4π rad/sample之间，阻带衰减为60 dB。
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π rad/sample范围内抑制归一化频率，通带纹波为1 dB。
一=唧唧声(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，．..'passbandropple'，1）;Y =过滤器（Dy，Ry）;
绘制这两个序列。垂直偏移第二个序列，使两者都可见。
绘图（[x y + 2]）
计算和绘制的横频谱X和y．使用512样本汉明窗口。指定500个相邻段和2048个DFT点之间的重叠样本。
xspectrogram (x, y,汉明(512),500年,2048年,“桠溪”）
将跨谱图值设置为零，小于-50 dB。
xspectrogram (x, y,汉明(512),500年,2048年,'MinThreshold', -50,“桠溪”）
谱图显示了滤波器增强或抑制的频率区域。
输入参数
输出参数
参考
Sanjit K. Mitra数字信号处理：基于计算机A方法．第二次。纽约：麦格劳 - 希尔，2001。
[2]奥本海姆，艾伦五，和Ronald W.谢弗，约翰·R·巴克。离散时间信号处理．第二次。上部鞍河，NJ：Prentice Hall，1999。
扩展功能
你也可以从以下列表中选择一个网站:
