主要内容GydF4y2Ba

WVD.GydF4y2Ba

Wigner-Ville分销和平滑的伪Wigner-Ville分销GydF4y2Ba

描述GydF4y2Ba

例子GydF4y2Ba

D.GydF4y2Ba=项(GydF4y2BaXGydF4y2Ba)GydF4y2Ba返回Wigner-Ville分销GydF4y2BaXGydF4y2Ba.GydF4y2Ba

例子GydF4y2Ba

D.GydF4y2Ba=项(GydF4y2BaXGydF4y2Ba那GydF4y2BaFS.GydF4y2Ba)GydF4y2Ba返回Wigner-Ville分销GydF4y2BaXGydF4y2Ba以速度进行抽样GydF4y2BaFS.GydF4y2Ba.GydF4y2Ba

例子GydF4y2Ba

D.GydF4y2Ba=项(GydF4y2BaXGydF4y2Ba那GydF4y2BaTS.GydF4y2Ba)GydF4y2Ba返回Wigner-Ville分销GydF4y2BaXGydF4y2Ba用时间间隔进行采样GydF4y2BaTS.GydF4y2Ba之间的样本。GydF4y2Ba

D.GydF4y2Ba=项(GydF4y2Ba___GydF4y2Ba,'smoothedpseudo')GydF4y2Ba的平滑伪Wigner-Ville分布GydF4y2BaXGydF4y2Ba.该功能使用输入信号的长度来选择用于时间和频率平滑的窗口的长度。此语法可以包括来自先前语法的输入参数的任何组合。GydF4y2Ba

例子GydF4y2Ba

D.GydF4y2Ba=项(GydF4y2Ba___GydF4y2Ba,'smoothedpseudo',GydF4y2Ba双胞胎GydF4y2Ba那GydF4y2BaFWIN.GydF4y2Ba)GydF4y2Ba指定时间窗口,GydF4y2Ba双胞胎GydF4y2Ba和频率窗口,GydF4y2BaFWIN.GydF4y2Ba,用于平滑。要使用默认窗口进行时间平滑或频率平滑,请指定相应的参数为空,GydF4y2Ba[]GydF4y2Ba.GydF4y2Ba

例子GydF4y2Ba

D.GydF4y2Ba=项(GydF4y2Ba___GydF4y2Ba,'smoothedpseudo',GydF4y2Ba名称,价值GydF4y2Ba)GydF4y2Ba使用名称值对参数指定Smooteed Pseudo Wigner-Ville分发的其他选项。您可以指定GydF4y2Ba双胞胎GydF4y2Ba和GydF4y2BaFWIN.GydF4y2Ba在此语法中,或者您可以省略它们。GydF4y2Ba

例子GydF4y2Ba

D.GydF4y2Ba=项(GydF4y2Ba___GydF4y2Ba,'minthreshold',GydF4y2Ba打GydF4y2Ba)GydF4y2Ba设置为零的那些元素GydF4y2BaD.GydF4y2Ba谁的幅度小于GydF4y2Ba打GydF4y2Ba.这种语法适用于Wigner-Ville分布和平滑的伪Wigner-Ville分布。GydF4y2Ba

例子GydF4y2Ba

[GydF4y2BaD.GydF4y2Ba那GydF4y2BaFGydF4y2Ba那GydF4y2BaT.GydF4y2Ba] = WVD(GydF4y2Ba___GydF4y2Ba)GydF4y2Ba还返回频率向量,GydF4y2BaFGydF4y2Ba,和一个时间向量,GydF4y2BaT.GydF4y2Ba,在那里GydF4y2BaD.GydF4y2Ba计算。GydF4y2Ba

项GydF4y2Ba___GydF4y2Ba)GydF4y2Ba在没有输出参数的情况下,绘制当前图中Wigner-Ville或平滑伪Wigner-Ville分布。GydF4y2Ba

例子GydF4y2Ba

全部收缩GydF4y2Ba

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

x = 0 (1001 1);x = 10 (500);y =罪(π* (0:1000)/ 2)';[d, f t] =项(x + y);GydF4y2Ba

绘制Wigner-Ville分布。GydF4y2Ba

ImagesC(T,F,D)轴GydF4y2BaxyGydF4y2BacolorbarGydF4y2Ba

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

通过调用GydF4y2BaWVD.GydF4y2Ba没有输出参数。GydF4y2Ba

wvd(x + y)GydF4y2Ba

图中包含一个坐标轴。标题Wigner-Ville分布的轴包含类型图像的对象。GydF4y2Ba

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

fs = 1000;t =(0:1 / fs:1.5)';x = cos(2 * pi * t * 200);GydF4y2Ba

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

WVD(X,FS)GydF4y2Ba

图中包含一个坐标轴。标题Wigner-Ville分布的轴包含类型图像的对象。GydF4y2Ba

添加到信号啁啾,其频率变化正弦,在250Hz和450 Hz之间。将信号转换为MATLAB®时间表。计算Wigner-Ville分发。GydF4y2Ba

X = X + vco(cos(2*pi*t),[250 450],fs);xt =时间表(x秒(t));项(xt)GydF4y2Ba

图中包含一个坐标轴。标题Wigner-Ville分布的轴包含类型图像的对象。GydF4y2Ba

设置为零的分配元素,幅度小于0。GydF4y2Ba

项(xt,GydF4y2Ba“MinThreshold”GydF4y2Ba, 0)GydF4y2Ba

图中包含一个坐标轴。标题Wigner-Ville分布的轴包含类型图像的对象。GydF4y2Ba

产生一个信号采样在1 kHz 1秒。信号的一个组成部分是一个啁啾,在测量过程中频率从100hz二次增加到400hz。信号的另一个组成部分是啁啾,它的频率在相同的延时下从350hz线性下降到50hz。GydF4y2Ba

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

fs = 1000;t = 0:1 / fs: 1;x =唧唧声(t, 100, 1400,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

在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 (2 j *π* 100 * cos(2 *π* 2 * t));X = x1 + x2 + randn(size(t))/10;GydF4y2Ba

计算并绘制信号的平滑伪Wigner Ville。窗口使用601样本汉明窗口及时使用305样本矩形窗口及时窗口分布。使用600个频率点显示。设置为零的分布的零组件,幅度小于GydF4y2Ba -GydF4y2Ba 50.GydF4y2Ba .GydF4y2Ba

WVD(X,DT,GydF4y2Ba'smoothedpseudo'GydF4y2Barectwin汉明(601),(305),GydF4y2Ba......GydF4y2Ba“NumFrequencyPoints”GydF4y2Ba600,GydF4y2Ba“MinThreshold”GydF4y2Ba,-50)GydF4y2Ba

图中包含一个坐标轴。具有标题平滑的伪Wigner-Ville分布的轴包含类型图像的对象。GydF4y2Ba

生成由四个高斯原子组成的信号。每个原子由高斯调制的正弦曲线组成。正弦波具有100Hz和400 Hz的频率。高斯以150毫秒为中心,350毫秒为中心,并具有差异GydF4y2Ba 0.GydF4y2Ba .GydF4y2Ba 01.GydF4y2Ba 2GydF4y2Ba .所有原子都有单位幅度。信号在1 kHz下对半秒进行采样。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

WVD(S,FS)GydF4y2Ba

图中包含一个坐标轴。标题Wigner-Ville分布的轴包含类型图像的对象。GydF4y2Ba

计算并显示信号的平滑伪Wigner-Ville分布。在时间和频率上的平滑可以减弱干涉项。GydF4y2Ba

WVD(S,FS,GydF4y2Ba'smoothedpseudo'GydF4y2Ba)GydF4y2Ba

图中包含一个坐标轴。具有标题平滑的伪Wigner-Ville分布的轴包含类型图像的对象。GydF4y2Ba

输入参数GydF4y2Ba

全部收缩GydF4y2Ba

输入信号,指定为向量或MATLABGydF4y2Ba®GydF4y2Ba包含单个矢量变量的时间表。GydF4y2Ba

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

例子:GydF4y2BaCOS(PI / 8 *(0:159))'+ RANDN(160,1)/ 10GydF4y2Ba指定嵌入白噪声的正弦曲线。GydF4y2Ba

例子:GydF4y2Ba时间表(秒(0:5)',兰特(6,1))GydF4y2Ba指定以1 Hz采样的随机变量5秒。GydF4y2Ba

数据类型:GydF4y2Ba单身的GydF4y2Ba|GydF4y2Ba双GydF4y2Ba
复数支持:金宝appGydF4y2Ba是的GydF4y2Ba

采样率,指定为正数字标量。GydF4y2Ba

采样时间,指定为aGydF4y2Ba期间GydF4y2Ba标量。GydF4y2Ba

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

  • 默认长度GydF4y2Ba双胞胎GydF4y2Ba最小的奇数是否大于或等于GydF4y2Ba轮GydF4y2Ba(GydF4y2Ba长度GydF4y2Ba(GydF4y2BaXGydF4y2Ba)/ 10)GydF4y2Ba.GydF4y2Ba

  • 默认长度GydF4y2BaFWIN.GydF4y2Ba最小的奇数是否大于或等于GydF4y2BaNF.GydF4y2Ba/ 4.GydF4y2Ba,在那里GydF4y2BaNF.GydF4y2Ba使用GydF4y2Banumfrequencypoints.GydF4y2Ba.GydF4y2Ba

每个窗口的长度必须小于或等于GydF4y2Ba2 *GydF4y2BaCEIL.GydF4y2Ba(GydF4y2Ba长度GydF4y2Ba(GydF4y2BaXGydF4y2Ba)/ 2)GydF4y2Ba.GydF4y2Ba

例子:GydF4y2Ba凯泽GydF4y2Ba(65,0.5)GydF4y2Ba指定具有0.5的形状因子的65样本kaiser窗口。GydF4y2Ba

最小非零值,指定为真正的标量。该功能设置为零的元素GydF4y2BaD.GydF4y2Ba谁的幅度小于GydF4y2Ba打GydF4y2Ba.GydF4y2Ba

名称值对参数GydF4y2Ba

指定可选的逗号分离对GydF4y2Ba名称,价值GydF4y2Ba参数。GydF4y2Ba姓名GydF4y2Ba是参数名称和GydF4y2Ba价值GydF4y2Ba是相应的价值。GydF4y2Ba姓名GydF4y2Ba必须出现在引号内。您可以以任何顺序指定多个名称和值对参数GydF4y2BaName1, Value1,…,的家GydF4y2Ba.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.时间点的数量必须至少是GydF4y2Ba2 *GydF4y2Ba长度GydF4y2Ba(GydF4y2Ba双胞胎GydF4y2Ba)GydF4y2Ba并且不能大于默认值。GydF4y2Ba

提示GydF4y2Ba

如果输入信号较大,可以减少时间点的数量,以降低内存要求,加快计算速度。GydF4y2Ba

输出参数GydF4y2Ba

全部收缩GydF4y2Ba

Wigner-Ville分布,作为矩阵返回。时间增加了跨越的列GydF4y2BaD.GydF4y2Ba并且频率随着行增加。矩阵的大小GydF4y2BaNGydF4y2BaFGydF4y2Ba×GydF4y2BaNGydF4y2BaT.GydF4y2Ba,在那里GydF4y2BaNGydF4y2BaFGydF4y2Ba是的长度GydF4y2BaFGydF4y2Ba和GydF4y2BaNGydF4y2BaT.GydF4y2Ba是的长度GydF4y2BaT.GydF4y2Ba.GydF4y2Ba

频率,作为向量返回。GydF4y2Ba

  • 如果输入有时间信息,那么GydF4y2BaFGydF4y2Ba包含Hz中表达的频率。GydF4y2Ba

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

时间瞬间,返回为矢量。GydF4y2Ba

  • 如果输入有时间信息,那么GydF4y2BaT.GydF4y2Ba包含以秒为单位的时间值。GydF4y2Ba

  • 如果输入没有时间信息,那么GydF4y2BaT.GydF4y2Ba包含样品号。GydF4y2Ba

更多关于GydF4y2Ba

全部收缩GydF4y2Ba

能量分布GydF4y2Ba

这GydF4y2BaWigner-Ville分销GydF4y2Ba提供信号的高分辨率时频表示。该分布在信号可视化、检测和估计等方面都有应用。GydF4y2Ba

对于连续信号GydF4y2BaXGydF4y2Ba(GydF4y2BaT.GydF4y2Ba)GydF4y2Ba,Wigner-Ville分布定义为GydF4y2Ba

WVD.GydF4y2Ba XGydF4y2Ba (GydF4y2Ba T.GydF4y2Ba 那GydF4y2Ba FGydF4y2Ba )GydF4y2Ba =GydF4y2Ba ∫GydF4y2Ba -GydF4y2Ba ∞GydF4y2Ba ∞GydF4y2Ba XGydF4y2Ba (GydF4y2Ba T.GydF4y2Ba +GydF4y2Ba τGydF4y2Ba 2GydF4y2Ba )GydF4y2Ba XGydF4y2Ba *GydF4y2Ba (GydF4y2Ba T.GydF4y2Ba -GydF4y2Ba τGydF4y2Ba 2GydF4y2Ba )GydF4y2Ba E.GydF4y2Ba -GydF4y2Ba jGydF4y2Ba 2GydF4y2Ba πGydF4y2Ba FGydF4y2Ba τGydF4y2Ba D.GydF4y2Ba τGydF4y2Ba .GydF4y2Ba

用于离散信号GydF4y2BaNGydF4y2Ba样品,分布变为GydF4y2Ba

WVD.GydF4y2Ba XGydF4y2Ba (GydF4y2Ba N.GydF4y2Ba 那GydF4y2Ba K.GydF4y2Ba )GydF4y2Ba =GydF4y2Ba σ.GydF4y2Ba mGydF4y2Ba =GydF4y2Ba -GydF4y2Ba NGydF4y2Ba NGydF4y2Ba XGydF4y2Ba (GydF4y2Ba N.GydF4y2Ba +GydF4y2Ba mGydF4y2Ba /GydF4y2Ba 2GydF4y2Ba )GydF4y2Ba XGydF4y2Ba *GydF4y2Ba (GydF4y2Ba N.GydF4y2Ba -GydF4y2Ba mGydF4y2Ba /GydF4y2Ba 2GydF4y2Ba )GydF4y2Ba E.GydF4y2Ba -GydF4y2Ba jGydF4y2Ba 2GydF4y2Ba πGydF4y2Ba K.GydF4y2Ba mGydF4y2Ba /GydF4y2Ba NGydF4y2Ba .GydF4y2Ba

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

Wigner-Ville分布包含通常使其解释复杂的干扰术语。为了锐化分布,可以使用低通窗过滤定义。Smoothed Pseudo Wigner-Ville分配使用独立的窗口在时间和频率下平滑:GydF4y2Ba

spwvd.GydF4y2Ba XGydF4y2Ba GGydF4y2Ba 那GydF4y2Ba HGydF4y2Ba (GydF4y2Ba T.GydF4y2Ba 那GydF4y2Ba FGydF4y2Ba )GydF4y2Ba =GydF4y2Ba ∫GydF4y2Ba -GydF4y2Ba ∞GydF4y2Ba ∞GydF4y2Ba GGydF4y2Ba (GydF4y2Ba T.GydF4y2Ba )GydF4y2Ba HGydF4y2Ba (GydF4y2Ba FGydF4y2Ba )GydF4y2Ba XGydF4y2Ba (GydF4y2Ba T.GydF4y2Ba +GydF4y2Ba τGydF4y2Ba 2GydF4y2Ba )GydF4y2Ba XGydF4y2Ba *GydF4y2Ba (GydF4y2Ba T.GydF4y2Ba -GydF4y2Ba τGydF4y2Ba 2GydF4y2Ba )GydF4y2Ba E.GydF4y2Ba -GydF4y2Ba jGydF4y2Ba 2GydF4y2Ba πGydF4y2Ba FGydF4y2Ba τGydF4y2Ba D.GydF4y2Ba τ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

在R2018B中介绍GydF4y2Ba