主要内容GydF4y2Ba

wvdGydF4y2Ba

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

描述GydF4y2Ba

例子GydF4y2Ba

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

例子GydF4y2Ba

D.GydF4y2Ba=wvd(GydF4y2BaXGydF4y2Ba那GydF4y2BaFS.GydF4y2Ba)GydF4y2Ba返回Wigner Ville分布,当GydF4y2BaXGydF4y2Ba以一定的速率进行采样GydF4y2BaFS.GydF4y2Ba。GydF4y2Ba

例子GydF4y2Ba

D.GydF4y2Ba=wvd(GydF4y2BaXGydF4y2Ba那GydF4y2BatsGydF4y2Ba)GydF4y2Ba返回Wigner Ville分布,当GydF4y2BaXGydF4y2Ba以时间间隔进行采样GydF4y2BatsGydF4y2Ba在样本之间。GydF4y2Ba

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

例子GydF4y2Ba

D.GydF4y2Ba=wvd(GydF4y2Ba___GydF4y2Ba,'smoothedpseudo',GydF4y2Ba成双的GydF4y2Ba那GydF4y2Ba弗温GydF4y2Ba)GydF4y2Ba指定时间窗口,GydF4y2Ba成双的GydF4y2Ba和频率窗口,GydF4y2Ba弗温GydF4y2Ba,用于平滑。若要使用时间或频率平滑的默认窗口,请将相应的参数指定为空,GydF4y2Ba[]GydF4y2Ba。GydF4y2Ba

例子GydF4y2Ba

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

例子GydF4y2Ba

D.GydF4y2Ba=wvd(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

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

例子GydF4y2Ba

全部收缩GydF4y2Ba

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

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

绘制Wigner-Ville分布。GydF4y2Ba

图像SC(t、f、d)轴GydF4y2BaxyGydF4y2Ba色条GydF4y2Ba

图中包含一个轴。轴包含图像类型的对象。GydF4y2Ba

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

wvd(x+y)GydF4y2Ba

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

生成一个信号,该信号由在1 kHz下采样1.5秒的200 Hz正弦波组成。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分布的轴包含类型为image的对象。GydF4y2Ba

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

x=x+vco(cos(2*pi*t),[250450],fs);xt=时刻表(秒(t),x);wvd(xt)GydF4y2Ba

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

将振幅小于0的分布元素设置为零。GydF4y2Ba

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

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

生成一个以1kHz采样1秒的信号。信号的一个成分是啁啾,在测量过程中,啁啾的频率从100Hz二次增加到400Hz。信号的另一个成分是啁啾,在同一时间段内频率从350 Hz线性降低到50 Hz。GydF4y2Ba

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

fs=1000;t=0:1/fs:1;x=chirp(t,100,1400,GydF4y2Ba“二次”GydF4y2Ba)+啁啾(t,350,1,50);GydF4y2Ba

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

wvd(x,fs)GydF4y2Ba

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

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

wvd(x,fs,GydF4y2Ba“平滑假想”GydF4y2Ba那GydF4y2Ba“NumFrequencyPoints”GydF4y2Ba,501,GydF4y2Ba'numtimepoints'GydF4y2Ba,502)GydF4y2Ba

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

增加时间点的数量,以便二次啁啾变得可见。GydF4y2Ba

wvd(x,fs,GydF4y2Ba“平滑假想”GydF4y2Ba那GydF4y2Ba“NumFrequencyPoints”GydF4y2Ba,501,GydF4y2Ba'numtimepoints'GydF4y2Ba,522)GydF4y2Ba

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

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

wvd(x,fs,GydF4y2Ba“平滑假想”GydF4y2Ba那GydF4y2Ba“NumFrequencyPoints”GydF4y2Ba,1000,GydF4y2Ba'numtimepoints'GydF4y2Ba,1502)GydF4y2Ba

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

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

fs=3000;t=0:1/fs:1-1/fs;dt=s(t(2)-t(1));x1=chirp(t,300,t(end),1300,GydF4y2Ba“二次”GydF4y2Ba); x2=exp(2j*pi*100*cos(2*pi*2*t));x=x1+x2+randn(尺寸(t))/10;GydF4y2Ba

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

wvd(x,dt,GydF4y2Ba“平滑假想”GydF4y2Ba,哈明(601),雷克特温(305),GydF4y2Ba...GydF4y2Ba“NumFrequencyPoints”GydF4y2Ba,600,GydF4y2Ba“MinThreshold”GydF4y2Ba,-50)GydF4y2Ba

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

生成由四个高斯原子组成的信号。每个原子都由一个高斯调制的正弦波组成。正弦频率为100 Hz和400 Hz。高斯分布以150毫秒和350毫秒为中心,方差为GydF4y2Ba 0.GydF4y2Ba 。GydF4y2Ba 01GydF4y2Ba 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=gaussFun([1 1],t[mu1 mu1 mu2 mu2],[f1 f2]);GydF4y2Ba

计算并显示信号的Wigner-Ville分布。可以具有负值的干扰术语在每对自动术语之间存在负值。GydF4y2Ba

WVD(S,FS)GydF4y2Ba

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

计算并显示信号的平滑伪Wigner-Ville分布。时间和频率平滑可衰减干扰项。GydF4y2Ba

WVD(S,FS,GydF4y2Ba“平滑假想”GydF4y2Ba)GydF4y2Ba

图中包含一个轴。标题为平滑伪Wigner-Ville分布的轴包含一个image类型的对象。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
复数支持:万博1manbetxGydF4y2Ba对GydF4y2Ba

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

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

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

  • 的默认长度GydF4y2Ba成双的GydF4y2Ba最小奇数整数是否大于或等于GydF4y2Ba圆形的GydF4y2Ba(GydF4y2Ba长度GydF4y2Ba(GydF4y2BaXGydF4y2Ba)/ 10)GydF4y2Ba。GydF4y2Ba

  • 的默认长度GydF4y2Ba弗温GydF4y2Ba最小奇数整数是否大于或等于GydF4y2Ba核因子GydF4y2Ba/4GydF4y2Ba哪里GydF4y2Ba核因子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必须出现在引号内。您可以以任何顺序指定多个名称和值对参数GydF4y2Ba名称1,值1,…,名称,值GydF4y2Ba。GydF4y2Ba

例子:GydF4y2Ba“NumFrequencyPoints”,201,“NumTimePoints”,300GydF4y2Ba计算201个频率点和300个时间点的Wigner-Ville分布。GydF4y2Ba

频率点的数量,指定为逗号分隔对,由GydF4y2Ba“NumFrequencyPoints”GydF4y2Ba和整数。此参数控制频率过采样的程度。频率点的数量必须至少为GydF4y2Ba(GydF4y2Ba长度GydF4y2Ba(GydF4y2Ba弗温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并且频率随着行增加。矩阵的大小GydF4y2BaN.GydF4y2BaFGydF4y2Ba×GydF4y2BaN.GydF4y2BaT.GydF4y2Ba哪里GydF4y2BaN.GydF4y2BaFGydF4y2Ba是的长度GydF4y2BaFGydF4y2Ba和GydF4y2BaN.GydF4y2BaT.GydF4y2Ba是的长度GydF4y2BaT.GydF4y2Ba。GydF4y2Ba

频率,作为向量返回。GydF4y2Ba

  • 如果输入有时间信息,则GydF4y2BaFGydF4y2Ba包含以Hz表示的频率。GydF4y2Ba

  • 如果输入没有时间信息,则GydF4y2BaFGydF4y2Ba包含以rad/sample表示的标准化频率。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

对于具有GydF4y2BaN.GydF4y2Ba样本,则分布变为GydF4y2Ba

WVD.GydF4y2Ba XGydF4y2Ba (GydF4y2Ba NGydF4y2Ba 那GydF4y2Ba K.GydF4y2Ba )GydF4y2Ba =GydF4y2Ba ∑GydF4y2Ba mGydF4y2Ba =GydF4y2Ba −GydF4y2Ba N.GydF4y2Ba N.GydF4y2Ba XGydF4y2Ba (GydF4y2Ba NGydF4y2Ba +GydF4y2Ba mGydF4y2Ba /GydF4y2Ba 2GydF4y2Ba )GydF4y2Ba XGydF4y2Ba *GydF4y2Ba (GydF4y2Ba NGydF4y2Ba −GydF4y2Ba mGydF4y2Ba /GydF4y2Ba 2GydF4y2Ba )GydF4y2Ba E.GydF4y2Ba −GydF4y2Ba jGydF4y2Ba 2GydF4y2Ba πGydF4y2Ba K.GydF4y2Ba mGydF4y2Ba /GydF4y2Ba N.GydF4y2Ba 。GydF4y2Ba

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

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

SPWVDGydF4y2Ba 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恩格尔伍德悬崖,新泽西州:普伦蒂斯大厅,1995年。GydF4y2Ba

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

[3] 计算二次时频分布的快速且内存有效的算法GydF4y2Ba应用与计算谐波分析GydF4y2Ba。卷。35,第2,2013号,第350-358页。GydF4y2Ba

扩展能力GydF4y2Ba

也可以看看GydF4y2Ba

职能GydF4y2Ba

话题GydF4y2Ba

在R2018B中介绍GydF4y2Ba