主要内容gydF4y2Ba

WVD.gydF4y2Ba

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

页面上倒塌gydF4y2Ba

语法gydF4y2Ba

d = wvd（x）gydF4y2Ba
d = WVD（X，FS）gydF4y2Ba
d = wvd（x，ts）gydF4y2Ba
d = wvd（gydF4y2Ba＿＿＿gydF4y2Ba, ' smoothedPseudo ')gydF4y2Ba
d = wvd（gydF4y2Ba＿＿＿gydF4y2Ba，'smoothedpseudo'，双胞胎，fwin）gydF4y2Ba
d = wvd（gydF4y2Ba＿＿＿gydF4y2Ba、“smoothedPseudo”名称,值)gydF4y2Ba
d = wvd（gydF4y2Ba＿＿＿gydF4y2Ba，'minthreshold'，thresh）gydF4y2Ba
[D，F，T] = WVD（gydF4y2Ba＿＿＿gydF4y2Ba）gydF4y2Ba
WVD（gydF4y2Ba＿＿＿gydF4y2Ba）gydF4y2Ba

描述gydF4y2Ba

例子gydF4y2Ba

dgydF4y2Ba= WVD（gydF4y2BaxgydF4y2Ba）gydF4y2Ba返回的Wigner-Ville分布gydF4y2BaxgydF4y2Ba．gydF4y2Ba

例子gydF4y2Ba

dgydF4y2Ba= WVD（gydF4y2BaxgydF4y2Ba，gydF4y2BafsgydF4y2Ba）gydF4y2Ba返回Wigner-Ville分销gydF4y2BaxgydF4y2Ba以速度进行抽样gydF4y2BafsgydF4y2Ba．gydF4y2Ba

例子gydF4y2Ba

dgydF4y2Ba= WVD（gydF4y2BaxgydF4y2Ba，gydF4y2BaTS.gydF4y2Ba）gydF4y2Ba返回Wigner-Ville分销gydF4y2BaxgydF4y2Ba用时间间隔进行采样gydF4y2BaTS.gydF4y2Ba样品之间。gydF4y2Ba

dgydF4y2Ba= WVD（gydF4y2Ba＿＿＿gydF4y2Ba, ' smoothedPseudo ')gydF4y2Ba返回Smoothed Pseudo Wigner-Ville分销gydF4y2BaxgydF4y2Ba．该函数使用输入信号的长度来选择用于时间和频率平滑的窗口的长度。该语法可以包含以前语法中的输入参数的任何组合。gydF4y2Ba

例子gydF4y2Ba

dgydF4y2Ba= WVD（gydF4y2Ba＿＿＿gydF4y2Ba“smoothedPseudo”,gydF4y2Ba双胞胎gydF4y2Ba，gydF4y2BaFWIN.gydF4y2Ba）gydF4y2Ba指定时间窗口，gydF4y2Ba双胞胎gydF4y2Ba，频率窗口，gydF4y2BaFWIN.gydF4y2Ba，用于平滑。要使用默认窗口进行时间或频率平滑，请将相应的参数指定为空，gydF4y2Ba[]gydF4y2Ba．gydF4y2Ba

例子gydF4y2Ba

dgydF4y2Ba= WVD（gydF4y2Ba＿＿＿gydF4y2Ba“smoothedPseudo”,gydF4y2Ba名称,值gydF4y2Ba）gydF4y2Ba使用名称-值对参数指定平滑伪Wigner-Ville发行版的其他选项。您可以指定gydF4y2Ba双胞胎gydF4y2Ba和gydF4y2BaFWIN.gydF4y2Ba在这个语法中，或者你可以省略它们。gydF4y2Ba

例子gydF4y2Ba

dgydF4y2Ba= WVD（gydF4y2Ba＿＿＿gydF4y2Ba，'minthreshold'，gydF4y2Ba阈值gydF4y2Ba）gydF4y2Ba将这些元素置零gydF4y2BadgydF4y2Ba谁的幅度小于gydF4y2Ba阈值gydF4y2Ba．这种语法适用于Wigner-Ville分布和平滑的伪Wigner-Ville分布。gydF4y2Ba

例子gydF4y2Ba

［gydF4y2BadgydF4y2Ba，gydF4y2BafgydF4y2Ba，gydF4y2BatgydF4y2Ba] =项gydF4y2Ba＿＿＿gydF4y2Ba）gydF4y2Ba还返回频率向量，gydF4y2BafgydF4y2Ba和时代的矢量，gydF4y2BatgydF4y2Ba，在那里gydF4y2BadgydF4y2Ba被计算。gydF4y2Ba

WVD（gydF4y2Ba＿＿＿gydF4y2Ba）gydF4y2Ba没有输出参数在当前数字中绘制Wigner-Ville或平滑的伪Wigner-Ville分布。gydF4y2Ba

例子gydF4y2Ba

全部折叠gydF4y2Ba

打开直播脚本gydF4y2Ba

产生一个1000样本脉冲和一个具有归一化频率的1000样本音调gydF4y2Ba πgydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ．计算两个信号之和的Wigner-Ville分布。gydF4y2Ba

x =零（1001,1）;x（500）= 10;Y = SIN（PI *（0：1000）/ 2）';[d，f，t] = wvd（x + y）;gydF4y2Ba

绘制Wigner-Ville分布。gydF4y2Ba

ImagesC（T，F，D）轴gydF4y2BaXY.gydF4y2Ba彩色栏gydF4y2Ba

图包含轴。轴包含类型图像的对象。gydF4y2Ba

通过呼叫重现结果gydF4y2BaWVD.gydF4y2Ba没有输出参数。gydF4y2Ba

wvd（x + y）gydF4y2Ba

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

打开直播脚本gydF4y2Ba

生成由在1kHz上采样的200 Hz正弦曲面组成的信号1.5秒。gydF4y2Ba

fs = 1000;t = (0:1 / fs: 1.5);x = cos(2 *π* t * 200);gydF4y2Ba

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

WVD（X，FS）gydF4y2Ba

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

给信号加上一个频率在250hz和450hz之间正弦变化的啁啾。转换信号到MATLAB®时间表。计算Wigner-Ville分布。gydF4y2Ba

x = x + vco（cos（2 * pi * t），[250 450]，fs）;xt =时间表（秒（t），x）;WVD（XT）gydF4y2Ba

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

设置为零的分配元素，幅度小于0。gydF4y2Ba

WVD（XT，gydF4y2Ba'minthreshold'gydF4y2Ba，0）gydF4y2Ba

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

打开直播脚本gydF4y2Ba

生成在1 kHz上采样的信号1秒。信号的一个组件是啁啾，在测量期间从100Hz到400Hz等于频率增加。信号的其他组件是啁啾，其在同一流逝中从350Hz到50Hz的频率下降。gydF4y2Ba

将信号存储在时间表中。gydF4y2Ba

fs = 1000;t = 0：1 / FS：1;x = chirp（t，100,1,400，gydF4y2Ba'二次'gydF4y2Ba）+ Chirp（T，350,1,50）;gydF4y2Ba

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

WVD（X，FS）gydF4y2Ba

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

计算信号的平滑伪Wigner-Ville分布。指定501个频率点和502个时间点。gydF4y2Ba

WVD（X，FS，gydF4y2Ba'smoothedpseudo'gydF4y2Ba，gydF4y2Ba'numfrequencypoints'gydF4y2Ba，501，gydF4y2Ba“NumTimePoints”gydF4y2Ba，502）gydF4y2Ba

图包含轴。具有标题平滑的伪Wigner-Ville分布的轴包含类型图像的对象。gydF4y2Ba

增加时间点的数量，使二次啁啾变得可见。gydF4y2Ba

WVD（X，FS，gydF4y2Ba'smoothedpseudo'gydF4y2Ba，gydF4y2Ba'numfrequencypoints'gydF4y2Ba，501，gydF4y2Ba“NumTimePoints”gydF4y2Ba, 522)gydF4y2Ba

图包含轴。具有标题平滑的伪Wigner-Ville分布的轴包含类型图像的对象。gydF4y2Ba

增加频率点和时间点以获得更清晰的图像。gydF4y2Ba

WVD（X，FS，gydF4y2Ba'smoothedpseudo'gydF4y2Ba，gydF4y2Ba'numfrequencypoints'gydF4y2Ba，1000，gydF4y2Ba“NumTimePoints”gydF4y2Ba, 1502)gydF4y2Ba

图包含轴。具有标题平滑的伪Wigner-Ville分布的轴包含类型图像的对象。gydF4y2Ba

打开直播脚本gydF4y2Ba

在3 kHz上采样1秒的双组分信号产生。第一组件是二次啁啾，其频率在测量期间从300 Hz增加到1300 Hz。第二组分是具有正弦变化频率内容的啁啾。该信号嵌入在白色高斯噪声中。将连续样本之间的时间表示为agydF4y2Ba持续时间gydF4y2Ba标量。gydF4y2Ba

FS = 3000;T = 0：1 / FS：1-1 / FS;dt =秒（t（2）-t（1））;X1 = Chirp（T，300，T（端），1300，gydF4y2Ba'二次'gydF4y2Ba）;x2 = exp（2j * pi * 100 * cos（2 * pi * 2 * t））;x = x1 + x2 + randn（尺寸（t））/ 10;gydF4y2Ba

计算并绘制信号的平滑伪维格纳维尔。使用601样本的汉明窗和使用305样本的矩形窗的频率分布。显示使用600个频率点。将振幅小于的分布分量设为零gydF4y2Ba -gydF4y2Ba 50.gydF4y2Ba ．gydF4y2Ba

WVD（X，DT，gydF4y2Ba'smoothedpseudo'gydF4y2Ba，汉明（601），矫正器（305），gydF4y2Ba......gydF4y2Ba'numfrequencypoints'gydF4y2Ba600，gydF4y2Ba'minthreshold'gydF4y2Ba, -50)gydF4y2Ba

图包含轴。具有标题平滑的伪Wigner-Ville分布的轴包含类型图像的对象。gydF4y2Ba

打开直播脚本gydF4y2Ba

生成由四个高斯原子组成的信号。每个原子由高斯调制的正弦曲线组成。正弦波具有100Hz和400 Hz的频率。高斯以150毫秒为中心，350毫秒为中心，并具有差异gydF4y2Ba 0gydF4y2Ba ．gydF4y2Ba 01.gydF4y2Ba 2gydF4y2Ba ．所有原子都有单位振幅。信号以1khz采样半秒。gydF4y2Ba

fs = 1000;t =（0：1 / fs：0.5）';F1 = 100;F2 = 400;mu1 = 0.15;mu2 = 0.35;Gaussfun = @（a，x，mu，f）exp（ - （x-mu）。^ 2 /（2 * 0.01 ^ 2））。* sin（2 * pi * f。* x）* a';S =高度（[1 1 1 1]，T，[MU1 MU1 MU2]，[F1 F2 F1 F2]）;gydF4y2Ba

计算并显示信号的Wigner-Ville分布。干扰项可以是负值，出现在每对自动项之间。gydF4y2Ba

项(s, fs)gydF4y2Ba

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

计算和显示信号的平滑伪Wigner-Ville分布。在时间和频率下平滑衰减干扰术语。gydF4y2Ba

项,fs,gydF4y2Ba'smoothedpseudo'gydF4y2Ba）gydF4y2Ba

图包含轴。具有标题平滑的伪Wigner-Ville分布的轴包含类型图像的对象。gydF4y2Ba

输入参数gydF4y2Ba

全部折叠gydF4y2Ba

输入信号，指定为向量或MATLABgydF4y2Ba®gydF4y2Ba包含单一向量变量的时间表。gydF4y2Ba

如果输入信号具有奇数长度，则该函数附加零以使长度均匀。gydF4y2Ba

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

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

数据类型：gydF4y2Ba单身的gydF4y2Ba|gydF4y2Ba双倍的gydF4y2Ba
复数的支持:万博1manbetxgydF4y2Ba是的gydF4y2Ba

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

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

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

  • 默认长度gydF4y2Ba双胞胎gydF4y2Ba是最小的奇数整数大于或等于gydF4y2Ba圆形的gydF4y2Ba（gydF4y2Ba长度gydF4y2Ba（gydF4y2BaxgydF4y2Ba) / 10)gydF4y2Ba．gydF4y2Ba

  • 默认长度gydF4y2BaFWIN.gydF4y2Ba是最小的奇数整数大于或等于gydF4y2BaNF.gydF4y2Ba/ 4.gydF4y2Ba， 在哪里gydF4y2BaNF.gydF4y2Ba使用gydF4y2BaNumFrequencyPointsgydF4y2Ba．gydF4y2Ba

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

例子：gydF4y2Ba凯泽gydF4y2Ba（65,0.5）gydF4y2Ba指定一个65个样本的Kaiser窗口，形状因子为0.5。gydF4y2Ba

最小非零值，指定为真正的标量。该功能设置为零的元素gydF4y2BadgydF4y2Ba谁的振幅小于gydF4y2Ba阈值gydF4y2Ba．gydF4y2Ba

名称-值对的观点gydF4y2Ba

指定可选的逗号分隔的对gydF4y2Ba名称,值gydF4y2Ba论点。gydF4y2Ba姓名gydF4y2Ba参数名和gydF4y2Ba价值gydF4y2Ba是相应的价值。gydF4y2Ba姓名gydF4y2Ba必须出现在引号内。可以以任意顺序指定多个名称和值对参数gydF4y2Baname1，value1，...，namen，valuengydF4y2Ba．gydF4y2Ba

例子：gydF4y2Ba'numfrequencypoints'，201，'numtimepoints'，300gydF4y2Ba在201频点和300个时间点计算Wigner-Ville分布。gydF4y2Ba

频率点数，指定为逗号分隔对gydF4y2Ba'numfrequencypoints'gydF4y2Ba和一个整数。此参数控制频率过采样程度。频率点的数量必须至少是gydF4y2Ba（gydF4y2Ba长度gydF4y2Ba（gydF4y2BaFWIN.gydF4y2Ba）+1）/ 2gydF4y2Ba且不能大于默认值。gydF4y2Ba

指定为逗号分隔对的时间点数gydF4y2Ba“NumTimePoints”gydF4y2Ba和一个偶数整数。这个论点控制了时间过采样的程度gydF4y2Ba［3］gydF4y2Ba（信号处理工具箱）gydF4y2Ba．时间点的数量必须至少是gydF4y2Ba2 *gydF4y2Ba长度gydF4y2Ba（gydF4y2Ba双胞胎gydF4y2Ba）gydF4y2Ba且不能大于默认值。gydF4y2Ba

提示gydF4y2Ba

如果输入信号大，则减少时间点的数量，以降低存储器要求并加快计算。gydF4y2Ba

输出参数gydF4y2Ba

全部折叠gydF4y2Ba

Wigner-Ville分布，以矩阵形式返回。时间在列中增加gydF4y2BadgydF4y2Ba，频率沿行向下增加。这个矩阵是有大小的gydF4y2BaNgydF4y2BafgydF4y2Ba×gydF4y2BaNgydF4y2BatgydF4y2Ba， 在哪里gydF4y2BaNgydF4y2BafgydF4y2Ba为长度gydF4y2BafgydF4y2Ba和gydF4y2BaNgydF4y2BatgydF4y2Ba为长度gydF4y2BatgydF4y2Ba．gydF4y2Ba

频率，作为向量返回。gydF4y2Ba

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

  • 如果输入没有时间信息，那么gydF4y2BafgydF4y2Ba包含Rad /样品中表达的归一化频率。gydF4y2Ba

时间瞬间，作为向量返回。gydF4y2Ba

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

  • 如果输入没有时间信息，那么gydF4y2BatgydF4y2Ba包含样品号。gydF4y2Ba

更多关于gydF4y2Ba

全部折叠gydF4y2Ba

Wigner-Ville分销gydF4y2Ba

的gydF4y2Ba能量分布gydF4y2Ba提供信号的高分辨率时间频率表示。分布在信号可视化，检测和估计中具有应用。gydF4y2Ba

对于连续信号gydF4y2BaxgydF4y2Ba（gydF4y2BatgydF4y2Ba）gydF4y2Ba，Wigner-Ville分布定义为gydF4y2Ba

项gydF4y2Ba xgydF4y2Ba （gydF4y2Ba tgydF4y2Ba ，gydF4y2Ba fgydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba ∫gydF4y2Ba -gydF4y2Ba ∞gydF4y2Ba ∞gydF4y2Ba xgydF4y2Ba （gydF4y2Ba tgydF4y2Ba +gydF4y2Ba τ.gydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba xgydF4y2Ba ＊gydF4y2Ba （gydF4y2Ba tgydF4y2Ba -gydF4y2Ba τ.gydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba egydF4y2Ba -gydF4y2Ba jgydF4y2Ba 2gydF4y2Ba πgydF4y2Ba fgydF4y2Ba τ.gydF4y2Ba dgydF4y2Ba τ.gydF4y2Ba ．gydF4y2Ba

用于离散信号gydF4y2BaNgydF4y2Ba样品，分布变为gydF4y2Ba

项gydF4y2Ba xgydF4y2Ba （gydF4y2Ba ngydF4y2Ba ，gydF4y2Ba kgydF4y2Ba ）gydF4y2Ba ＝gydF4y2Ba σ.gydF4y2Ba 米gydF4y2Ba ＝gydF4y2Ba -gydF4y2Ba NgydF4y2Ba NgydF4y2Ba xgydF4y2Ba （gydF4y2Ba ngydF4y2Ba +gydF4y2Ba 米gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba xgydF4y2Ba ＊gydF4y2Ba （gydF4y2Ba ngydF4y2Ba -gydF4y2Ba 米gydF4y2Ba /gydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba egydF4y2Ba -gydF4y2Ba jgydF4y2Ba 2gydF4y2Ba πgydF4y2Ba kgydF4y2Ba 米gydF4y2Ba /gydF4y2Ba NgydF4y2Ba ．gydF4y2Ba

对于奇数值gydF4y2Ba米gydF4y2Ba，该定义要求在半整数采样值处对信号进行评估。因此需要插值，这就需要对离散傅里叶变换进行零填充以避免混叠。gydF4y2Ba

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

spwvd.gydF4y2Ba xgydF4y2Ba 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 xgydF4y2Ba ＊gydF4y2Ba （gydF4y2Ba tgydF4y2Ba -gydF4y2Ba τ.gydF4y2Ba 2gydF4y2Ba ）gydF4y2Ba egydF4y2Ba -gydF4y2Ba jgydF4y2Ba 2gydF4y2Ba πgydF4y2Ba fgydF4y2Ba τ.gydF4y2Ba dgydF4y2Ba τ.gydF4y2Ba ．gydF4y2Ba

参考gydF4y2Ba

[1]科恩，莱昂。gydF4y2Ba时频分析：理论与应用gydF4y2Ba．Englewood Cliffs，NJ：1995年Prentice-Hall。gydF4y2Ba

[2] Mallat，Stéphane。gydF4y2Ba信号处理的小波巡回赛gydF4y2Ba．第二版。圣地亚哥，加利福尼亚州：1999年学术出版社。gydF4y2Ba

[3] O'Toole，John M.和Boualem Boashash。“用于计算二次时频分布的快速和内存高效算法。”gydF4y2Ba应用和计算谐波分析gydF4y2Ba．第35卷，第2期，2013年，350-358页。gydF4y2Ba

扩展能力gydF4y2Ba

另请参阅gydF4y2Ba

功能gydF4y2Ba

主题gydF4y2Ba

介绍了R2018bgydF4y2Ba

小波工具箱文档gydF4y2Ba

万博1manbetx

尝试matlab，sim万博1manbetxulink等产品s manbetx 845gydF4y2Ba

得到审判现在gydF4y2Ba