主要内容gydF4y2Ba

互Wigner-Ville分布gydF4y2Ba

交叉维格纳-维尔分布和交叉平滑伪维格纳-维尔分布gydF4y2Ba

描述gydF4y2Ba

例子gydF4y2Ba

dgydF4y2Ba=xwvd(gydF4y2BaxgydF4y2Ba,gydF4y2BaygydF4y2Ba)gydF4y2Ba返回的交叉维格纳-维尔分布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,gydF4y2Ba弗温gydF4y2Ba)gydF4y2Ba指定时间窗口,gydF4y2Ba成双的gydF4y2Ba,频率窗口,gydF4y2Ba弗温gydF4y2Ba,用于平滑。要使用默认窗口进行时间平滑或频率平滑,请指定相应的参数为空,gydF4y2Ba[]gydF4y2Ba.gydF4y2Ba

dgydF4y2Ba=xwvd(gydF4y2Ba___gydF4y2Ba,‘smoothedPseudo’,‘NumFrequencyPoints’,gydF4y2Ba核因子gydF4y2Ba)gydF4y2Ba使用以下公式计算交叉平滑的伪Wigner-Ville分布:gydF4y2Ba核因子gydF4y2Ba频率点。您可以指定gydF4y2Ba成双的gydF4y2Ba和gydF4y2Ba弗温gydF4y2Ba在这个语法中,或者你可以省略它们。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

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

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));gydF4y2Ba

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

项(x + y, fs)gydF4y2Ba

Figure包含axes对象。标题为Wigner Ville Distribution的axes对象包含image类型的对象。gydF4y2Ba

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

xwvd(x,y,fs)gydF4y2Ba

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

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

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

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

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

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

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

创建由高斯原子组成的参考信号,该原子以1 kHz的频率采样1秒。高斯原子是由高斯原子调制的正弦波。指定50 Hz的正弦波频率。高斯原子以64毫秒为中心,其方差为gydF4y2Ba 0gydF4y2Ba .gydF4y2Ba 01gydF4y2Ba 2gydF4y2Ba .gydF4y2Ba

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);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“平滑假想”gydF4y2Ba); wvd(xa+xr,fs,gydF4y2Ba“平滑假想”gydF4y2Ba)gydF4y2Ba

图中包含一个轴对象。标题为平滑伪维格纳-维尔分布的轴对象包含一个类型为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的轴对象包含一个类型为image的对象。gydF4y2Ba

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

D = w + c;d = abs(真正的(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的对象。gydF4y2Ba

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

测试= 2*tc - mu;= 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指定以1 Hz采样4秒的随机变量。gydF4y2Ba

数据类型:gydF4y2Ba仅有一个的gydF4y2Ba|gydF4y2Ba双重的gydF4y2Ba
复数的支持:金宝appgydF4y2Ba对gydF4y2Ba

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

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

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

  • 的默认长度gydF4y2Ba成双的gydF4y2Ba最小的奇数是否大于或等于gydF4y2Ba轮gydF4y2Ba(gydF4y2Ba长度gydF4y2Ba(gydF4y2BaxgydF4y2Ba) / 10)gydF4y2Ba.gydF4y2Ba

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

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

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

频率点的数量,指定为整数。此参数控制频率过采样的程度。频率点的数量必须至少为gydF4y2Ba(gydF4y2Ba长度gydF4y2Ba(gydF4y2Ba弗温gydF4y2Ba)+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

交叉项能量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分布采用独立窗口在时间和频率上进行平滑: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. 恩格尔伍德悬崖,新泽西州:普伦蒂斯大厅,1995年。gydF4y2Ba

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

Malnar, Damir, Victor Sucic, Boualem Boashash。“一种基于交叉项几何的方法,利用交叉维格纳-维尔分布进行组件瞬时频率估计。”在gydF4y2Ba第十一届国际信息科学、信号处理及其应用会议(ISSPA)gydF4y2Ba, 1217 - 1222页。蒙特利尔:IEEEgydF4y2Ba®gydF4y2Ba,2012年。gydF4y2Ba

扩展能力gydF4y2Ba

另请参阅gydF4y2Ba

功能gydF4y2Ba

主题gydF4y2Ba

介绍了R2018bgydF4y2Ba