zhyuer 发表于 2006-5-17 09:04

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

最好给个程序例子，谢谢

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

roltem 发表于 2008-4-13 08:35

用时频工具箱很容易就画出来了
给一个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分布
=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 发表于 2009-3-23 21:12

？？？

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

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

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

rxs1980 发表于 2009-11-16 09:44

你真强大！

wuyifeng 发表于 2010-3-4 23:55

请教

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

fadeblack 发表于 2010-3-6 23:23

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

suxh 发表于 2010-8-30 20:44

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

龙真诚 发表于 2011-5-7 09:37

值得收藏

sjh2100 发表于 2011-10-19 09:19

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

miwueshine 发表于 2012-2-8 10:32

{:{39}:}

elena 发表于 2012-4-11 16:29

不错{:{39}:}

大海无边 发表于 2012-12-17 10:46

太好了，又学习了！
页: [1]
查看完整版本: 求Wigner-Ville分布算法的时频图