主要内容gydF4y2Ba

xwvdgydF4y2Ba

交叉Wigner-Ville分布和交叉平滑伪Wigner-Ville分布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

生成2个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

计算信号和的维格纳-维尔分布。gydF4y2Ba

项(x + y, fs)gydF4y2Ba

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

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

fs xwvd (x, y)gydF4y2Ba

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

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

fs = 3000;t = (0:1 / fs: 1 - 1 / fs) ';x1 =唧唧喳喳(800 t, 1400 t(结束),);x2 = 2 *唧唧喳喳(1000 t, 200 t(结束),);gydF4y2Ba

将信号存储为时间表。计算和绘制两个渠道的交叉维格纳-维尔分布。gydF4y2Ba

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

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

利用已知的参考信号和交叉维格纳-维尔分布计算信号的瞬时频率。gydF4y2Ba

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

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

创建“未知”信号进行分析,包括一个啁啾。信号在0.4秒时突然开始,半秒后突然结束。在这个过程中,啁啾的频率从400hz线性下降到100hz。gydF4y2Ba

f0 = 400;f1 = 100;xa = 0(大小(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 =项(xa + xr, fs,gydF4y2Ba“smoothedPseudo”gydF4y2Ba);项(xa + xr, fs,gydF4y2Ba“smoothedPseudo”gydF4y2Ba)gydF4y2Ba

图中包含一个axes对象。标题为平滑伪Wigner-Ville分布的axis对象包含一个类型为image的对象。gydF4y2Ba

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

画出交叉维格纳-维尔分布的实部。gydF4y2Ba

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

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

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

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

图中包含一个axes对象。坐标轴对象包含一个image类型的对象。gydF4y2Ba

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

ff = tfridge (c, fc);电视=和(c) > 0;ff = ff(电视);tc = tc(电视);持有gydF4y2Ba在gydF4y2Ba情节(tc、ff,gydF4y2Ba“r——”gydF4y2Ba,gydF4y2Ba“线宽”gydF4y2Ba, 2)gydF4y2Ba从gydF4y2Ba

图中包含一个axes对象。坐标轴对象包含两个类型为image、line的对象。gydF4y2Ba

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

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

图中包含一个axes对象。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
复数的支持:gydF4y2Ba是的gydF4y2Ba

抽样速率,指定为正数值标量。gydF4y2Ba

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

用于平滑的时间和频率窗口,指定为奇数长度的向量。默认情况下,gydF4y2BaxwvdgydF4y2Ba使用具有形状因子的Kaiser窗口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

交叉维格纳-维尔分布,作为矩阵返回。的列的时间增加了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

交叉项能量gydF4y2Ba

连续信号gydF4y2BaxgydF4y2Ba(gydF4y2BatgydF4y2Ba)gydF4y2Ba而且gydF4y2BaygydF4y2Ba(gydF4y2BatgydF4y2Ba)gydF4y2Ba,gydF4y2Ba交叉项能量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分布使用独立窗口在时间和频率上进行平滑: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

[1]科恩,莱昂。gydF4y2Ba时频分析:理论与应用gydF4y2Ba.恩格尔伍德悬崖,新泽西州:Prentice-Hall, 1995。gydF4y2Ba

[2] Mallat,史蒂芬。gydF4y2Ba信号处理的小波之旅gydF4y2Ba.第二版。圣地亚哥,加州:学术出版社,1999年。gydF4y2Ba

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

扩展功能gydF4y2Ba

版本历史gydF4y2Ba

介绍了R2018bgydF4y2Ba

另请参阅gydF4y2Ba

功能gydF4y2Ba

主题gydF4y2Ba

Baidu
map