平滑伪Wigner-Ville分布的时频图

beefyy

一频率为20Hz的周期信号，计算其平滑伪Wigner-Ville分布，然后画其频谱图！
结果频谱图中出现20和30两条频率线！   采用频率为100Hz！
若换为10Hz，则出现10和40两条，规律为相加都为50!
这是为什么画谱图就不会出现这样的问题！

[ 本帖最后由 zhangnan3509 于 2007-10-27 12:03 编辑 ]

破凰

应该是你编写的程序出了问题。
下面我写的代码：
clear;
clc;
f=20;
fs=100;
t=1/fs:1/fs:1;
y=sin(2*pi*f*t);
z=hilbert(y');
tfrpwv(z);
得到的时频图如下

beefyy

谢谢你的帮助，是我自己编写错了！

beefyy

想再向你请教下，Threshold=5%在哪里体现出来？

jingrenzhi

你好!你在程序里的z=hilbert(y');是啥意思?为啥还要做这个变换?谢谢

beefyy

这个是指求y的解析信号啊！

jingrenzhi

是不是关于横轴不对称的信号,不能进行z=hilbert(y')这个变换?

破凰

原帖由 beefyy 于 2007-10-28 10:23 发表
想再向你请教下，Threshold=5%在哪里体现出来？

好像就是在时频图上体现。我的理解是：阈值越大，将能量较小的那些点上的系数置0，中间那根线就越接近理想的时频分布。
下面是Threshold=30%时得到的时频图

beefyy

再请教下，这个阈值在程序里怎么体现，我没见到你的程序里有关于阈值的语句啊。
谢谢！

破凰

原帖由 beefyy 于 2007-10-29 11:12 发表
再请教下，这个阈值在程序里怎么体现，我没见到你的程序里有关于阈值的语句啊。
谢谢！

运行程序后出现界面



点黑框中的选项，再在命令窗口输入阈值即可。

zhlong

呵呵，楼主是想看程序里面实现该功能的语句。

破凰

时频工具箱的函数好像没有提供接口，只有自己编了。

zhlong

里面的函数都可以打开编辑的呀，我以前也找过这个功能的语句，不过里面程序太多，没仔细去找。

破凰

我的意思是用tfrpwv不行，:lol
应该是通过这个函数tfrview改的阈值

beefyy

谢谢破凰，我又学了一招，以前我都是直接用[tfr,t,f]=tfrpwv(z);这语句，然后就画图了。