主要内容GyD.F4y2Ba

互Wigner-Ville分布GyD.F4y2Ba

交叉维格纳-维尔分布和交叉平滑伪维格纳-维尔分布GyD.F4y2Ba

折叠所有页面GyD.F4y2Ba

句法GyD.F4y2Ba

d = xwvd（x，y）GyD.F4y2Ba
d = xwvd（x，y，fs）GyD.F4y2Ba
d = xwvd（x，y，ts）GyD.F4y2Ba
d = xwvd（GyD.F4y2Ba___GyD.F4y2Ba，'smoothedpseudo'）GyD.F4y2Ba
d = xwvd（GyD.F4y2Ba___GyD.F4y2Ba“smoothedPseudo”,双胞胎,fwin)GyD.F4y2Ba
d = xwvd（GyD.F4y2Ba___GyD.F4y2Ba,‘smoothedPseudo’,‘NumFrequencyPoints’,nf)GyD.F4y2Ba
d = xwvd（GyD.F4y2Ba___GyD.F4y2Ba、“MinThreshold”打)GyD.F4y2Ba
[D，F，T] = XWVD（GyD.F4y2Ba___GyD.F4y2Ba）GyD.F4y2Ba
XWVD（GyD.F4y2Ba___GyD.F4y2Ba）GyD.F4y2Ba

描述GyD.F4y2Ba

例子GyD.F4y2Ba

D.GyD.F4y2Ba=xwvd(GyD.F4y2BaXGyD.F4y2Ba那GyD.F4y2BayGyD.F4y2Ba）GyD.F4y2Ba返回十字架Wigner-Ville分销GyD.F4y2BaXGyD.F4y2Ba和GyD.F4y2BayGyD.F4y2Ba．GyD.F4y2Ba

例子GyD.F4y2Ba

D.GyD.F4y2Ba=xwvd(GyD.F4y2BaXGyD.F4y2Ba那GyD.F4y2BayGyD.F4y2Ba那GyD.F4y2BaFS.GyD.F4y2Ba）GyD.F4y2Ba在以下情况下返回交叉Wigner-Ville分布：GyD.F4y2BaXGyD.F4y2Ba和GyD.F4y2BayGyD.F4y2Ba抽样率是多少GyD.F4y2BaFS.GyD.F4y2Ba．GyD.F4y2Ba

D.GyD.F4y2Ba=xwvd(GyD.F4y2BaXGyD.F4y2Ba那GyD.F4y2BayGyD.F4y2Ba那GyD.F4y2BatsGyD.F4y2Ba）GyD.F4y2Ba在以下情况下返回交叉Wigner-Ville分布：GyD.F4y2BaXGyD.F4y2Ba和GyD.F4y2BayGyD.F4y2Ba用时间间隔进行采样GyD.F4y2BatsGyD.F4y2Ba之间的样本。GyD.F4y2Ba

D.GyD.F4y2Ba=xwvd(GyD.F4y2Ba___GyD.F4y2Ba，'smoothedpseudo'）GyD.F4y2Ba返回交叉平滑的伪Wigner-Ville分销GyD.F4y2BaXGyD.F4y2Ba和GyD.F4y2BayGyD.F4y2Ba. 该函数使用输入信号的长度来选择用于时间和频率平滑的窗口长度。此语法可以包括来自以前语法的输入参数的任意组合。GyD.F4y2Ba

D.GyD.F4y2Ba=xwvd(GyD.F4y2Ba___GyD.F4y2Ba，'smoothedpseudo'，GyD.F4y2Ba成双的GyD.F4y2Ba那GyD.F4y2Ba弗温GyD.F4y2Ba）GyD.F4y2Ba指定时间窗口，GyD.F4y2Ba成双的GyD.F4y2Ba和频率窗口，GyD.F4y2Ba弗温GyD.F4y2Ba，用于平滑。要使用默认窗口进行时间平滑或频率平滑，请指定相应的参数为空，GyD.F4y2Ba［］GyD.F4y2Ba．GyD.F4y2Ba

D.GyD.F4y2Ba=xwvd(GyD.F4y2Ba___GyD.F4y2Ba,‘smoothedPseudo’,‘NumFrequencyPoints’,GyD.F4y2Ba核因子GyD.F4y2Ba）GyD.F4y2Ba使用以下公式计算交叉平滑的伪Wigner-Ville分布：GyD.F4y2Ba核因子GyD.F4y2Ba频率点。您可以指定GyD.F4y2Ba成双的GyD.F4y2Ba和GyD.F4y2Ba弗温GyD.F4y2Ba在此语法中，或者您可以省略它们。GyD.F4y2Ba

D.GyD.F4y2Ba=xwvd(GyD.F4y2Ba___GyD.F4y2Ba“MinThreshold”，GyD.F4y2Ba打GyD.F4y2Ba）GyD.F4y2Ba设置为零的那些元素GyD.F4y2BaD.GyD.F4y2Ba其振幅小于GyD.F4y2Ba打GyD.F4y2Ba．该语法既适用于交叉Wigner-Ville分布，也适用于交叉平滑伪Wigner-Ville分布。GyD.F4y2Ba

[GyD.F4y2BaD.GyD.F4y2Ba那GyD.F4y2BaFGyD.F4y2Ba那GyD.F4y2BaT.GyD.F4y2Ba] = XWVD（GyD.F4y2Ba___GyD.F4y2Ba）GyD.F4y2Ba还返回频率向量，GyD.F4y2BaFGyD.F4y2Ba，和一个时间向量，GyD.F4y2BaT.GyD.F4y2Ba，在GyD.F4y2BaD.GyD.F4y2Ba计算。GyD.F4y2Ba

例子GyD.F4y2Ba

XWVD（GyD.F4y2Ba___GyD.F4y2Ba）GyD.F4y2Ba没有输出争论绘制十字架的真正的一部分，或者在当前数字中交叉平滑的伪Wigner-Ville分布。GyD.F4y2Ba

例子GyD.F4y2Ba

全部收缩GyD.F4y2Ba

打开实时脚本GyD.F4y2Ba

产生两个信号采样在1 kHz为1秒，并嵌入白噪声。一个信号是频率为150赫兹的正弦波。另一种信号是频率在200hz和400hz之间正弦变化的啁啾。噪声的方差为GyD.F4y2Ba 0.GyD.F4y2Ba ．GyD.F4y2Ba 1GyD.F4y2Ba 2GyD.F4y2Ba ．GyD.F4y2Ba

fs = 1000;t = (0:1 / fs: 1) ';X = cos(2*pi*t*150) + 0.1*randn(size(t));y = vco (cos(3 *π* t) (200 400), fs) + 0.1 * randn(大小(t));GyD.F4y2Ba

计算信号总和的Wigner-Ville分布。GyD.F4y2Ba

项(x + y, fs)GyD.F4y2Ba

图包含轴。标题Wigner-Ville分布的轴包含类型图像的对象。GyD.F4y2Ba

计算并绘制信号的交叉Wigner-Ville分布。交叉分布对应于Wigner-Ville分布的交叉项。GyD.F4y2Ba

xwvd（x，y，fs）GyD.F4y2Ba

图中包含一个轴。标题为Cross Wigner Ville分布的轴包含一个image类型的对象。GyD.F4y2Ba

打开实时脚本GyD.F4y2Ba

生成由两个啁啾组成的双通道信号。信号以3 kHz的频率采样1秒。第一个啁啾的初始频率为400 Hz，在采样结束时达到800 Hz。第二次啁啾从500 Hz开始，最后达到1000 Hz。第二个啁啾的振幅是第一个啁啾的两倍。GyD.F4y2Ba

FS = 3000;t =（0：1 / fs：1-1 / fs）';x1 = chirp（t，1400，t（端），800）;x2 = 2 * chirp（t，200，t（end），1000）;GyD.F4y2Ba

将信号存储为时间表。计算并绘制两个通道的交叉Wigner-Ville分布。GyD.F4y2Ba

xt =时间表（秒（t），x1，x2）;XWVD（XT（：，1），XT（:,2））GyD.F4y2Ba

图中包含一个轴。标题为Cross Wigner Ville分布的轴包含一个image类型的对象。GyD.F4y2Ba

打开实时脚本GyD.F4y2Ba

通过使用已知的参考信号和交叉Wigner-Ville分布来计算信号的瞬时频率。GyD.F4y2Ba

创建一个参考信号组成的高斯原子采样在1 kHz 1秒。高斯原子是由高斯调制的正弦波。指定正弦波频率为50hz。高斯分布以64毫秒为中心，方差为GyD.F4y2Ba 0.GyD.F4y2Ba ．GyD.F4y2Ba 01GyD.F4y2Ba 2GyD.F4y2Ba ．GyD.F4y2Ba

fs=1e3；t=（0:1/fs:1-1/fs）；（mu=0.064；sigma=0.01；fsin=50；xr=exp（-（t-mu）。^2/（2*sigma^2））*sin（2*pi*fsin*t）；GyD.F4y2Ba

创建“未知”信号以分析，由啁啾组成。该信号突然在0.4秒时开始，突然半秒突然结束。在这种流逝中，啁啾的频率从400Hz到100Hz线性降低。GyD.F4y2Ba

F0 = 400;F1 = 100;xa =零（尺寸（t））;Xa（t> 0.4＆t <= 0.9）= chirp（（0：1 / fs：0.5-1 / fs）'，f0,0.5，f1）;GyD.F4y2Ba

创建一个由未知信号和参考信号之和组成的双分量信号。结果的光滑伪Wigner-Ville分布提供了一个“理想”时频表示。GyD.F4y2Ba

计算和显示平滑的伪Wigner-Ville分布。GyD.F4y2Ba

w =项(xa + xr, fs,GyD.F4y2Ba“平滑假想”GyD.F4y2Ba); wvd（xa+xr，fs，GyD.F4y2Ba“平滑假想”GyD.F4y2Ba）GyD.F4y2Ba

图中包含一个轴。标题为平滑伪Wigner-Ville分布的轴包含类型为image的对象。GyD.F4y2Ba

计算未知信号和参考信号的交叉Wigner-Ville分布。取分布的绝对值，将振幅小于10的元素设置为零。交叉Wigner-Ville分布等于两分量信号的交叉项。GyD.F4y2Ba

绘制交叉Wigner-Ville分布的真正部分。GyD.F4y2Ba

[c，fc，tc]=xwvd（xa，xr，fs）；c=abs（c）；c（c<10）=0；xwvd（xa，xr，fs）GyD.F4y2Ba

图中包含一个轴。标题为Cross Wigner Ville分布的轴包含一个image类型的对象。GyD.F4y2Ba

通过向交叉Wigner-Ville分布添加理想的时频表示来增强Wigner-Ville横向。Wigner-Ville分布的横向于参考信号和未知信号之间的中途发生。GyD.F4y2Ba

D = w + c;d = abs(真正的(d));显示亮度图像(tc、fc、d)轴GyD.F4y2BaxyGyD.F4y2BacolorbarGyD.F4y2Ba

图中包含一个轴。这些轴包含一个image类型的对象。GyD.F4y2Ba

识别并绘制与交叉项相对应的高能脊。为了隔离山脊，找出交叉分布能量非零的时间值。GyD.F4y2Ba

ff=tfridge（c，fc）；tv=sum（c）>0；ff=ff（tv）；tc=tc（tv）；保持GyD.F4y2Ba在…上GyD.F4y2Ba绘图（tc、ff、，GyD.F4y2Ba“r——”GyD.F4y2Ba那GyD.F4y2Ba“线宽”GyD.F4y2Ba，2）持有GyD.F4y2Ba离开GyD.F4y2Ba

图中包含一个轴。轴包含2个类型为image、line的对象。GyD.F4y2Ba

使用脊和参考功能重建未知信号的瞬时频率。绘制瞬时频率作为时间的函数。GyD.F4y2Ba

测试= 2*tc - mu;= 2*ff - fsin;情节(测试、节日)GyD.F4y2Ba

图包含轴。轴包含类型线的对象。GyD.F4y2Ba

输入参数GyD.F4y2Ba

全部收缩GyD.F4y2Ba

输入信号，指定为vectors或matlabGyD.F4y2Ba®GyD.F4y2Ba每个时间表只包含一个矢量变量。GyD.F4y2BaXGyD.F4y2Ba和GyD.F4y2BayGyD.F4y2Ba必须两者都是vector或两者都是时间表，并且必须具有相同的长度。GyD.F4y2Ba

如果输入信号的长度为奇数，则函数会附加一个0使其长度为偶数。GyD.F4y2Ba

例子:GyD.F4y2BaCOS（PI / 8 *（0：159））'+ RANDN（160,1）/ 10GyD.F4y2Ba指定嵌入白噪声的正弦曲线。GyD.F4y2Ba

例子:GyD.F4y2Ba时间表（秒（0：5）'，兰特（6,1））GyD.F4y2Ba指定以1 Hz采样4秒的随机变量。GyD.F4y2Ba

数据类型：GyD.F4y2Ba仅有一个的GyD.F4y2Ba|GyD.F4y2Ba双GyD.F4y2Ba
复数支持：万博1manbetxGyD.F4y2Ba对GyD.F4y2Ba

采样率，指定为正数字标量。GyD.F4y2Ba

采样时间，指定为aGyD.F4y2Ba期间GyD.F4y2Ba标量。GyD.F4y2Ba

用于平滑的时间和频率窗口，指定为奇数长度的向量。默认情况下，GyD.F4y2Ba互Wigner-Ville分布GyD.F4y2Ba使用带有形状因子的kaiser窗口GyD.F4y2BaβGyD.F4y2Ba= 20.GyD.F4y2Ba．GyD.F4y2Ba

  • 的默认长度GyD.F4y2Ba成双的GyD.F4y2Ba最小的奇数是否大于或等于GyD.F4y2Ba轮GyD.F4y2Ba（GyD.F4y2Ba长度GyD.F4y2Ba（GyD.F4y2BaXGyD.F4y2Ba）/ 10）GyD.F4y2Ba．GyD.F4y2Ba

  • 的默认长度GyD.F4y2Ba弗温GyD.F4y2Ba最小的奇数是否大于或等于GyD.F4y2Ba核因子GyD.F4y2Ba/4GyD.F4y2Ba．GyD.F4y2Ba

每个窗口的长度必须小于或等于GyD.F4y2Ba2*GyD.F4y2BaCEIL.GyD.F4y2Ba（GyD.F4y2Ba长度GyD.F4y2Ba（GyD.F4y2BaXGyD.F4y2Ba)/2)GyD.F4y2Ba．GyD.F4y2Ba

例子:GyD.F4y2Ba凯撒GyD.F4y2Ba(65,0.5)GyD.F4y2Ba指定具有0.5的形状因子的65样本kaiser窗口。GyD.F4y2Ba

频率点的数量，指定为整数。此参数控制频率过采样的程度。频率点的数量必须至少为GyD.F4y2Ba（GyD.F4y2Ba长度GyD.F4y2Ba（GyD.F4y2Ba弗温GyD.F4y2Ba)+1)/2GyD.F4y2Ba并且不能大于默认值。GyD.F4y2Ba

最小非零值，指定为实标量。函数将这些元素设置为零GyD.F4y2BaD.GyD.F4y2Ba谁的幅度小于GyD.F4y2Ba打GyD.F4y2Ba．GyD.F4y2Ba

输出参数GyD.F4y2Ba

全部收缩GyD.F4y2Ba

交叉Wigner-Ville分布，作为矩阵返回。时间会随着时间的推移而增加GyD.F4y2BaD.GyD.F4y2Ba并且频率随着行增加。矩阵的大小GyD.F4y2BaNGyD.F4y2BaFGyD.F4y2Ba×GyD.F4y2BaNGyD.F4y2BaT.GyD.F4y2Ba哪里GyD.F4y2BaNGyD.F4y2BaFGyD.F4y2Ba是的长度GyD.F4y2BaFGyD.F4y2Ba和GyD.F4y2BaNGyD.F4y2BaT.GyD.F4y2Ba是的长度GyD.F4y2BaT.GyD.F4y2Ba．GyD.F4y2Ba

频率，作为向量返回。GyD.F4y2Ba

  • 如果输入有时间信息，则GyD.F4y2BaFGyD.F4y2Ba包含以Hz表示的频率。GyD.F4y2Ba

  • 如果输入没有时间信息，则GyD.F4y2BaFGyD.F4y2Ba包含以rad/sample表示的标准化频率。GyD.F4y2Ba

时间瞬间，作为向量返回。GyD.F4y2Ba

  • 如果输入有时间信息，则GyD.F4y2BaT.GyD.F4y2Ba包含以秒为单位的时间值。GyD.F4y2Ba

  • 如果输入没有时间信息，则GyD.F4y2BaT.GyD.F4y2Ba包含样本编号。GyD.F4y2Ba

时间点的数量固定为GyD.F4y2Ba4*GyD.F4y2BaCEIL.GyD.F4y2Ba（GyD.F4y2Ba长度GyD.F4y2Ba（GyD.F4y2BaXGyD.F4y2Ba)/2)GyD.F4y2Ba．GyD.F4y2Ba

更多关于GyD.F4y2Ba

全部收缩GyD.F4y2Ba

交叉项能量GyD.F4y2Ba

连续信号GyD.F4y2BaXGyD.F4y2Ba（GyD.F4y2BaT.GyD.F4y2Ba）GyD.F4y2Ba和GyD.F4y2BayGyD.F4y2Ba（GyD.F4y2BaT.GyD.F4y2Ba）GyD.F4y2Ba， 这GyD.F4y2Ba交叉Wigner-Ville分布GyD.F4y2Ba定义为GyD.F4y2Ba

XWVDGyD.F4y2Ba XGyD.F4y2Ba 那GyD.F4y2Ba yGyD.F4y2Ba （GyD.F4y2Ba T.GyD.F4y2Ba 那GyD.F4y2Ba FGyD.F4y2Ba ）GyD.F4y2Ba =GyD.F4y2Ba ∫GyD.F4y2Ba −GyD.F4y2Ba ∞GyD.F4y2Ba ∞GyD.F4y2Ba XGyD.F4y2Ba （GyD.F4y2Ba T.GyD.F4y2Ba +GyD.F4y2Ba τGyD.F4y2Ba 2GyD.F4y2Ba ）GyD.F4y2Ba yGyD.F4y2Ba *GyD.F4y2Ba （GyD.F4y2Ba T.GyD.F4y2Ba −GyD.F4y2Ba τGyD.F4y2Ba 2GyD.F4y2Ba ）GyD.F4y2Ba E.GyD.F4y2Ba −GyD.F4y2Ba jGyD.F4y2Ba 2GyD.F4y2Ba πGyD.F4y2Ba FGyD.F4y2Ba τGyD.F4y2Ba D.GyD.F4y2Ba τGyD.F4y2Ba ．GyD.F4y2Ba

对于具有GyD.F4y2BaNGyD.F4y2Ba样本，则分布变为GyD.F4y2Ba

XWVDGyD.F4y2Ba XGyD.F4y2Ba 那GyD.F4y2Ba yGyD.F4y2Ba （GyD.F4y2Ba NGyD.F4y2Ba 那GyD.F4y2Ba K.GyD.F4y2Ba ）GyD.F4y2Ba =GyD.F4y2Ba ∑GyD.F4y2Ba mGyD.F4y2Ba =GyD.F4y2Ba −GyD.F4y2Ba NGyD.F4y2Ba NGyD.F4y2Ba XGyD.F4y2Ba （GyD.F4y2Ba NGyD.F4y2Ba +GyD.F4y2Ba mGyD.F4y2Ba /GyD.F4y2Ba 2GyD.F4y2Ba ）GyD.F4y2Ba yGyD.F4y2Ba *GyD.F4y2Ba （GyD.F4y2Ba NGyD.F4y2Ba −GyD.F4y2Ba mGyD.F4y2Ba /GyD.F4y2Ba 2GyD.F4y2Ba ）GyD.F4y2Ba E.GyD.F4y2Ba −GyD.F4y2Ba jGyD.F4y2Ba 2GyD.F4y2Ba πGyD.F4y2Ba K.GyD.F4y2Ba mGyD.F4y2Ba /GyD.F4y2Ba NGyD.F4y2Ba ．GyD.F4y2Ba

对于的奇数值GyD.F4y2BamGyD.F4y2Ba，定义需要评估半整数样本值的信号。因此，它需要插值，这使得必须零填充离散的傅里叶变换以避免混叠。GyD.F4y2Ba

交叉Wigner-Ville分布包含通常使其解释复杂化的干扰术语。为了锐化分布，可以使用低通窗过滤定义。交叉平滑的伪Wigner-Ville分配使用独立的窗口在时间和频率下平滑：GyD.F4y2Ba

XSPWVDGyD.F4y2Ba XGyD.F4y2Ba 那GyD.F4y2Ba yGyD.F4y2Ba GGyD.F4y2Ba 那GyD.F4y2Ba HGyD.F4y2Ba （GyD.F4y2Ba T.GyD.F4y2Ba 那GyD.F4y2Ba FGyD.F4y2Ba ）GyD.F4y2Ba =GyD.F4y2Ba ∫GyD.F4y2Ba −GyD.F4y2Ba ∞GyD.F4y2Ba ∞GyD.F4y2Ba GGyD.F4y2Ba （GyD.F4y2Ba T.GyD.F4y2Ba ）GyD.F4y2Ba HGyD.F4y2Ba （GyD.F4y2Ba FGyD.F4y2Ba ）GyD.F4y2Ba XGyD.F4y2Ba （GyD.F4y2Ba T.GyD.F4y2Ba +GyD.F4y2Ba τGyD.F4y2Ba 2GyD.F4y2Ba ）GyD.F4y2Ba yGyD.F4y2Ba *GyD.F4y2Ba （GyD.F4y2Ba T.GyD.F4y2Ba −GyD.F4y2Ba τGyD.F4y2Ba 2GyD.F4y2Ba ）GyD.F4y2Ba E.GyD.F4y2Ba −GyD.F4y2Ba jGyD.F4y2Ba 2GyD.F4y2Ba πGyD.F4y2Ba FGyD.F4y2Ba τGyD.F4y2Ba D.GyD.F4y2Ba τGyD.F4y2Ba ．GyD.F4y2Ba

参考文献GyD.F4y2Ba

[1]科恩，莱昂。GyD.F4y2Ba时频分析：理论与应用GyD.F4y2Ba. 恩格尔伍德悬崖，新泽西州：普伦蒂斯大厅，1995年。GyD.F4y2Ba

[2] 玛拉特，圣菲。GyD.F4y2Ba信号处理的小波之旅GyD.F4y2Ba．第二版。加州圣地亚哥:学术出版社，1999。GyD.F4y2Ba

Malnar, Damir, Victor Sucic, Boualem Boashash。“一种基于交叉项几何的方法，利用交叉维格纳-维尔分布进行组件瞬时频率估计。”在GyD.F4y2Ba第11届信息科学国际会议，信号处理及其申请（ISSPA）GyD.F4y2Ba, 1217 - 1222页。蒙特利尔:IEEEGyD.F4y2Ba®GyD.F4y2Ba, 2012年。GyD.F4y2Ba

扩展能力GyD.F4y2Ba

也可以看看GyD.F4y2Ba

职能GyD.F4y2Ba

话题GyD.F4y2Ba

在R2018B中介绍GyD.F4y2Ba

信号处理工具箱文档GyD.F4y2Ba

万博1manbetx

基于MATLAB的信号处理深度学习GyD.F4y2Ba

基于MATLAB的信号处理深度学习GyD.F4y2Ba

下载白皮书GyD.F4y2Ba