主要内容gydF4y2Ba

xwvdgydF4y2Ba

交叉Wigner-Ville分布和交叉平滑伪Wigner-Ville分布gydF4y2Ba

折叠所有页面gydF4y2Ba

语法gydF4y2Ba

D = xwvd(x,y)gydF4y2Ba
D = xwvd(x,y,fs)gydF4y2Ba
D = xwvd(x,y,ts)gydF4y2Ba
D = xwvd(gydF4y2Ba___gydF4y2Ba, ' smoothedPseudo ')gydF4y2Ba
D = xwvd(gydF4y2Ba___gydF4y2Ba“smoothedPseudo”,双胞胎,fwin)gydF4y2Ba
D = xwvd(gydF4y2Ba___gydF4y2Ba,‘smoothedPseudo’,‘NumFrequencyPoints’,nf)gydF4y2Ba
D = xwvd(gydF4y2Ba___gydF4y2Ba、“MinThreshold”打)gydF4y2Ba
[d,f,t] = xwvd(gydF4y2Ba___gydF4y2Ba）gydF4y2Ba
xwvd (gydF4y2Ba___gydF4y2Ba）gydF4y2Ba

描述gydF4y2Ba

例子gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2BaxgydF4y2Ba，gydF4y2BaygydF4y2Ba）gydF4y2Ba的交叉Wigner-Ville分布gydF4y2BaxgydF4y2Ba而且gydF4y2BaygydF4y2Ba．gydF4y2Ba

例子gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2BaxgydF4y2Ba，gydF4y2BaygydF4y2Ba，gydF4y2BafsgydF4y2Ba）gydF4y2Ba返回交叉Wigner-Ville分布gydF4y2BaxgydF4y2Ba而且gydF4y2BaygydF4y2Ba以一定的速率采样gydF4y2BafsgydF4y2Ba．gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2BaxgydF4y2Ba，gydF4y2BaygydF4y2Ba，gydF4y2BatsgydF4y2Ba）gydF4y2Ba返回交叉Wigner-Ville分布gydF4y2BaxgydF4y2Ba而且gydF4y2BaygydF4y2Ba是否有时间间隔进行采样gydF4y2BatsgydF4y2Ba之间的样本。gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2Ba___gydF4y2Ba, ' smoothedPseudo ')gydF4y2Ba的交叉平滑伪Wigner-Ville分布gydF4y2BaxgydF4y2Ba而且gydF4y2BaygydF4y2Ba．该函数使用输入信号的长度来选择用于时间和频率平滑的窗口长度。此语法可以包括来自以前语法的输入参数的任何组合。gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2Ba___gydF4y2Ba“smoothedPseudo”,gydF4y2Ba双胞胎gydF4y2Ba，gydF4y2BafwingydF4y2Ba）gydF4y2Ba指定时间窗口，gydF4y2Ba双胞胎gydF4y2Ba，频率窗口，gydF4y2BafwingydF4y2Ba，用于平滑。要使用默认窗口进行时间或频率平滑，指定相应的参数为空，gydF4y2Ba［］gydF4y2Ba．gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2Ba___gydF4y2Ba,‘smoothedPseudo’,‘NumFrequencyPoints’,gydF4y2BanfgydF4y2Ba）gydF4y2Ba计算交叉平滑的伪Wigner-Ville分布gydF4y2BanfgydF4y2Ba频率点。你可以指定gydF4y2Ba双胞胎gydF4y2Ba而且gydF4y2BafwingydF4y2Ba在这个语法中，或者你可以省略它们。gydF4y2Ba

dgydF4y2Ba= xwvd (gydF4y2Ba___gydF4y2Ba“MinThreshold”,gydF4y2Ba打gydF4y2Ba）gydF4y2Ba的元素设置为零gydF4y2BadgydF4y2Ba谁的振幅小于gydF4y2Ba打gydF4y2Ba．这种语法既适用于交叉Wigner-Ville分布，也适用于交叉平滑的伪Wigner-Ville分布。gydF4y2Ba

［gydF4y2BadgydF4y2Ba，gydF4y2BafgydF4y2Ba，gydF4y2BatgydF4y2Ba= xwvd(gydF4y2Ba___gydF4y2Ba）gydF4y2Ba也返回一个频率向量，gydF4y2BafgydF4y2Ba为时间向量，gydF4y2BatgydF4y2Ba，在gydF4y2BadgydF4y2Ba计算。gydF4y2Ba

例子gydF4y2Ba

xwvd (gydF4y2Ba___gydF4y2Ba）gydF4y2Ba在没有输出参数的情况下，绘制当前图中交叉Wigner-Ville分布或交叉平滑伪Wigner-Ville分布的实部。gydF4y2Ba

例子gydF4y2Ba

全部折叠gydF4y2Ba

打开实时脚本gydF4y2Ba

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

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

计算信号和的Wigner-Ville分布。gydF4y2Ba

项(x + y, fs)gydF4y2Ba

图中包含一个轴对象。标题为Wigner-Ville Distribution的axes对象包含一个image类型的对象。gydF4y2Ba

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

fs xwvd (x, y)gydF4y2Ba

图中包含一个轴对象。标题为Cross Wigner-Ville Distribution的axes对象包含一个image类型的对象。gydF4y2Ba

打开实时脚本gydF4y2Ba

产生一个由两个啁啾组成的双通道信号。信号以3千赫频率采样一秒。第一个啁啾的初始频率为400hz，在采样结束时达到800hz。第二个啁啾从500hz开始，最后达到1000hz。第二个啁啾的振幅是第一个啁啾的两倍。gydF4y2Ba

Fs = 3000;T = (0:1/fs:1-1/fs)';X1 = chirp(t,1400,t(end)，800);X2 = 2*chirp(t,200,t(end)，1000);gydF4y2Ba

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

Xt =时间表(秒(t)，x1,x2);xwvd (xt (: 1), xt (:, 2))gydF4y2Ba

图中包含一个轴对象。标题为Cross Wigner-Ville Distribution的axes对象包含一个image类型的对象。gydF4y2Ba

打开实时脚本gydF4y2Ba

利用已知的参考信号和交叉Wigner-Ville分布计算信号的瞬时频率。gydF4y2Ba

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

Fs = 1e3;T = (0:1/fs:1-1/fs)';Mu = 0.064;σ = 0.01;Fsin = 50;(xr = exp () - tμ。^ 2 /(2 *σ^ 2))。*罪(2 *π*国际* t);gydF4y2Ba

创建要分析的“未知”信号，由啁啾组成。信号在0.4秒突然开始，半秒后突然结束。在这个过程中，啁啾的频率从400hz线性下降到100hz。gydF4y2Ba

F0 = 400;F1 = 100;Xa = 0 (size(t));xa (t > 0.4 & t < = 0.9) =唧唧声((0:1 / fs: 0.5 - 1 / fs)”,f0, 0.5, f1);gydF4y2Ba

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

计算并显示平滑的伪Wigner-Ville分布。gydF4y2Ba

W = wvd(xa+xr,fs，gydF4y2Ba“smoothedPseudo”gydF4y2Ba）;项(xa + xr, fs,gydF4y2Ba“smoothedPseudo”gydF4y2Ba）gydF4y2Ba

图中包含一个轴对象。标题为Smoothed Pseudo Wigner-Ville Distribution的axes对象包含一个image类型的对象。gydF4y2Ba

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

绘制交叉Wigner-Ville分布的实部。gydF4y2Ba

[c,fc,tc] = xwvd(xa,xr,fs);C = abs(C);C (C <10) = 0;xwvd (xa xr, fs)gydF4y2Ba

图中包含一个轴对象。标题为Cross Wigner-Ville Distribution的axes对象包含一个image类型的对象。gydF4y2Ba

通过在交叉Wigner-Ville分布中添加理想时频表示来增强Wigner-Ville交叉项。Wigner-Ville分布的交叉项出现在参考信号和未知信号之间的中间。gydF4y2Ba

D = w + c;D = abs(real(D));显示亮度图像(tc、fc、d)轴gydF4y2BaxygydF4y2BacolorbargydF4y2Ba

图中包含一个轴对象。axis对象包含一个image类型的对象。gydF4y2Ba

识别并绘制出交叉项对应的高能脊线。为了隔离脊，找到交叉分布具有非零能量的时间值。gydF4y2Ba

Ff = tfridge(c,fc);TV = sum(c)>0;Ff = Ff (tv);Tc = Tc (tv);持有gydF4y2Ba在gydF4y2Ba情节(tc、ff,gydF4y2Ba“r——”gydF4y2Ba，gydF4y2Ba“线宽”gydF4y2Ba, 2)gydF4y2Ba从gydF4y2Ba

图中包含一个轴对象。axis对象包含image、line类型的2个对象。gydF4y2Ba

利用脊线和参考函数重构未知信号的瞬时频率。画出瞬时频率作为时间的函数。gydF4y2Ba

tEst = 2*tc - mu;fEst = 2*ff - fsin;情节(测试、节日)gydF4y2Ba

图中包含一个轴对象。axis对象包含一个line类型的对象。gydF4y2Ba

输入参数gydF4y2Ba

全部折叠gydF4y2Ba

输入信号，指定为矢量或MATLABgydF4y2Ba®gydF4y2Ba每个时间表都包含一个矢量变量。gydF4y2BaxgydF4y2Ba而且gydF4y2BaygydF4y2Ba必须都是向量或者都是时间表，而且长度必须相同。gydF4y2Ba

如果输入信号长度为奇数，则该函数附加一个0使长度为偶数。gydF4y2Ba

例子:gydF4y2Ba因为(π/ 8 * (0:159))' + randn (160 1) / 10gydF4y2Ba指定嵌入白噪声中的正弦波。gydF4y2Ba

例子:gydF4y2Ba时间表(秒(0:5)”,兰德(6,1))gydF4y2Ba以1hz频率采样4秒的随机变量。gydF4y2Ba

数据类型:gydF4y2Ba单gydF4y2Ba|gydF4y2Ba双gydF4y2Ba
复数支持:金宝appgydF4y2Ba是的gydF4y2Ba

采样率，指定为正数值标量。gydF4y2Ba

采样时间，指定为agydF4y2Ba持续时间gydF4y2Ba标量。gydF4y2Ba

用于平滑的时间和频率窗口，指定为奇数长度的向量。默认情况下,gydF4y2BaxwvdgydF4y2Ba使用具有形状因子的凯撒窗gydF4y2BaβgydF4y2Ba= 20gydF4y2Ba．gydF4y2Ba

  • 的默认长度。gydF4y2Ba双胞胎gydF4y2Ba最小的奇数是否大于等于gydF4y2Ba轮gydF4y2Ba（gydF4y2Ba长度gydF4y2Ba（gydF4y2BaxgydF4y2Ba) / 10)gydF4y2Ba．gydF4y2Ba

  • 的默认长度。gydF4y2BafwingydF4y2Ba最小的奇数是否大于等于gydF4y2BanfgydF4y2Ba/ 4gydF4y2Ba．gydF4y2Ba

每个窗口的长度必须小于或等于gydF4y2Ba2 *gydF4y2Ba装天花板gydF4y2Ba（gydF4y2Ba长度gydF4y2Ba（gydF4y2BaxgydF4y2Ba) / 2)gydF4y2Ba．gydF4y2Ba

例子:gydF4y2Ba凯撒gydF4y2Ba》(65,0.5)gydF4y2Ba指定形状因子为0.5的65个样本Kaiser窗口。gydF4y2Ba

频率点的数目，指定为整数。这个参数控制频率的过采样程度。频率点个数至少为gydF4y2Ba（gydF4y2Ba长度gydF4y2Ba（gydF4y2BafwingydF4y2Ba) + 1) / 2gydF4y2Ba且不能大于默认值。gydF4y2Ba

最小非零值，指定为实标量。函数将的元素设为零gydF4y2BadgydF4y2Ba谁的振幅小于gydF4y2Ba打gydF4y2Ba．gydF4y2Ba

输出参数gydF4y2Ba

全部折叠gydF4y2Ba

交叉Wigner-Ville分布，作为矩阵返回。的列中时间增加gydF4y2BadgydF4y2Ba，频率沿行递减。矩阵有大小gydF4y2BaNgydF4y2BafgydF4y2Ba×gydF4y2BaNgydF4y2BatgydF4y2Ba,在那里gydF4y2BaNgydF4y2BafgydF4y2Ba的长度gydF4y2BafgydF4y2Ba而且gydF4y2BaNgydF4y2BatgydF4y2Ba的长度gydF4y2BatgydF4y2Ba．gydF4y2Ba

频率，作为向量返回。gydF4y2Ba

  • 如果输入有时间信息，那么gydF4y2BafgydF4y2Ba包含以Hz表示的频率。gydF4y2Ba

  • 如果输入没有时间信息，则gydF4y2BafgydF4y2Ba包含以rad/sample表示的归一化频率。gydF4y2Ba

时间瞬间，以矢量形式返回。gydF4y2Ba

  • 如果输入有时间信息，那么gydF4y2BatgydF4y2Ba包含以秒表示的时间值。gydF4y2Ba

  • 如果输入没有时间信息，则gydF4y2BatgydF4y2Ba包含样本编号。gydF4y2Ba

时间点的数量固定为gydF4y2Ba4 *gydF4y2Ba装天花板gydF4y2Ba（gydF4y2Ba长度gydF4y2Ba（gydF4y2BaxgydF4y2Ba) / 2)gydF4y2Ba．gydF4y2Ba

更多关于gydF4y2Ba

全部折叠gydF4y2Ba

交叉Wigner-Ville分布gydF4y2Ba

对于连续信号gydF4y2BaxgydF4y2Ba（gydF4y2BatgydF4y2Ba）gydF4y2Ba而且gydF4y2BaygydF4y2Ba（gydF4y2BatgydF4y2Ba）gydF4y2Ba,gydF4y2Ba交叉Wigner-Ville分布gydF4y2Ba定义为gydF4y2Ba

XWVDgydF4y2Ba xgydF4y2Ba ，gydF4y2Ba ygydF4y2Ba （gydF4y2Ba tgydF4y2Ba ，gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba ∫gydF4y2Ba −gydF4y2Ba ∞gydF4y2Ba ∞gydF4y2Ba xgydF4y2Ba （gydF4y2Ba tgydF4y2Ba +gydF4y2Ba τgydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba ygydF4y2Ba ＊gydF4y2Ba （gydF4y2Ba tgydF4y2Ba −gydF4y2Ba τgydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba egydF4y2Ba −gydF4y2Ba jgydF4y2Ba 2gydF4y2Ba πgydF4y2Ba fgydF4y2Ba τgydF4y2Ba dgydF4y2Ba τgydF4y2Ba ．gydF4y2Ba

对于一个离散信号gydF4y2BaNgydF4y2Ba样本，分布变成gydF4y2Ba

XWVDgydF4y2Ba xgydF4y2Ba ，gydF4y2Ba ygydF4y2Ba （gydF4y2Ba ngydF4y2Ba ，gydF4y2Ba kgydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba ∑gydF4y2Ba 米gydF4y2Ba ＝gydF4y2Ba −gydF4y2Ba NgydF4y2Ba NgydF4y2Ba xgydF4y2Ba （gydF4y2Ba ngydF4y2Ba +gydF4y2Ba 米gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba ygydF4y2Ba ＊gydF4y2Ba （gydF4y2Ba ngydF4y2Ba −gydF4y2Ba 米gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba egydF4y2Ba −gydF4y2Ba jgydF4y2Ba 2gydF4y2Ba πgydF4y2Ba kgydF4y2Ba 米gydF4y2Ba /gydF4y2Ba NgydF4y2Ba ．gydF4y2Ba

的奇数值gydF4y2Ba米gydF4y2Ba，定义要求在半整数采样值处计算信号。因此，它需要插值，这使得有必要对离散傅里叶变换进行零垫以避免混叠。gydF4y2Ba

交叉Wigner-Ville分布包含干扰项，这往往使其解释复杂化。为了锐化分布，可以使用低通窗口过滤定义。交叉平滑伪Wigner-Ville分布使用独立窗口在时间和频率上进行平滑:gydF4y2Ba

XSPWVDgydF4y2Ba xgydF4y2Ba ，gydF4y2Ba ygydF4y2Ba ggydF4y2Ba ，gydF4y2Ba HgydF4y2Ba （gydF4y2Ba tgydF4y2Ba ，gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba ∫gydF4y2Ba −gydF4y2Ba ∞gydF4y2Ba ∞gydF4y2Ba ggydF4y2Ba （gydF4y2Ba tgydF4y2Ba ）gydF4y2Ba HgydF4y2Ba （gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba xgydF4y2Ba （gydF4y2Ba tgydF4y2Ba +gydF4y2Ba τgydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba ygydF4y2Ba ＊gydF4y2Ba （gydF4y2Ba tgydF4y2Ba −gydF4y2Ba τgydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba egydF4y2Ba −gydF4y2Ba jgydF4y2Ba 2gydF4y2Ba πgydF4y2Ba fgydF4y2Ba τgydF4y2Ba dgydF4y2Ba τgydF4y2Ba ．gydF4y2Ba

参考文献gydF4y2Ba

科恩，里昂。gydF4y2Ba时频分析:理论与应用gydF4y2Ba．恩格尔伍德悬崖，新泽西州:Prentice-Hall, 1995。gydF4y2Ba

[2]马拉特，Stéphane。gydF4y2Ba信号处理小波之旅gydF4y2Ba．第二版。圣地亚哥，加州:学术出版社，1999年。gydF4y2Ba

马纳尔，达米尔，维克多·苏契奇，博阿什。“一种基于交叉项几何的方法，用于使用交叉Wigner-Ville分布的组件瞬时频率估计。”在gydF4y2Ba第十一届信息科学、信号处理及其应用国际会议(ISSPA)gydF4y2Ba，第1217-1222页。蒙特利尔:IEEEgydF4y2Ba®gydF4y2Ba, 2012年。gydF4y2Ba

扩展功能gydF4y2Ba

另请参阅gydF4y2Ba

功能gydF4y2Ba

主题gydF4y2Ba

在R2018b中引入gydF4y2Ba

小波工具箱文档gydF4y2Ba

金宝app

试试MATLAB、Sim金宝appulink和其他产品下载188bet金宝搏gydF4y2Ba

现在就去审判吧gydF4y2Ba