求Wigner-Ville分布算法的时频图

zhyuer

最好给个程序例子，谢谢

[ 本帖最后由 eight 于 2008-4-14 20:55 编辑 ]

roltem

用时频工具箱很容易就画出来了
给一个LFM的例子：

%产生线性调频信号
N=64; sig=fmlin(N,0,0.3);
%加噪声
sig=sigmerge(sig,noisecg(N),0);
%绘制时域波形
figure(1); plot(real(sig),'LineWidth',2); xlabel('时间 t'); ylabel('幅值 A');
nfft=1024; dsp=abs(fft(sig,nfft)).^2;
%绘制频谱图
figure(2); plot((0:nfft/2-1)/nfft,dsp(1:nfft/2)); xlabel('频率 f'); ylabel('功率谱 PSD');
%计算信号的Wigner-Ville分布
[tfr,t,f]=tfrwv(sig);
%绘图, 绘制等高线
figure(3); contour(t,f,abs(tfr).^2); xlabel('时间 t'); ylabel('频率 f'); %%绘制三维图
figure(4); mesh(t,f,abs(tfr).^2); xlabel('时间 t'); ylabel('频率 f'); zlabel('幅值 A');

同时也想请教各位大侠：时频工具箱的函数都是归一化处理的结果，如果要设置载频、采样频率、脉宽等参数之后又该如何进行分析呢，最好附上程序，谢谢！

[ 本帖最后由 ChaChing 于 2010-3-10 11:29 编辑 ]

xielna

刚学matlab，上面程序输入后，没有出现图像，
显示??? Undefined function or method 'fmlin' for input arguments of type 'double'.，
能指点一下是什么问题吗？

晓得了晓得是怎么回事了，O(∩_∩)O~~~~~~~~~~:@)

[ 本帖最后由 ChaChing 于 2010-3-5 09:21 编辑 ]

rxs1980

你真强大！

wuyifeng

各位大侠，为什么我用这个贴子里的程序执行后得到以下的图形。可是为什么频率坐标只有最到到0.5呢？我的信号明明是比较高的频率但是用这个程序执行后却只有这么小的频率？

fadeblack

因为调用tfrwv程序中已经对最终的频率f进行了归一化，这个可以到tfrwv源程序中改动

suxh

也可以将归一化的频率乘以抽样频率得出

龙真诚

值得收藏

sjh2100

2楼的程序原来是个”例子“啊。没看出来，厉害厉害。

miwueshine

elena

不错

大海无边

太好了，又学习了！