画小波变换的时频图-面包板社区

http://www.chinavib.com/forum/thread-46138-1-3.html

jzc1983:

一调频信号：x=sin(2*pi*(10*t+1).*t)+sin(2*pi*(30*t+1).*t);
用tfrscalo怎么只得到一个频率的图线，tfrscalo(xd',1:k1,N,0.001,0.50);这个N是根据信号得的么？N不同得到的时频就不同，小弟今天才接触小波，谢谢了

采样频率为1000HZ
-----------------------------------------------

zhlong:

给你一个例子，按图示信息就可以做出来了。
t=1:512;
t=t*1/2000;
f1=600;
f2=100;
f3=200;
a=1000;
x=sin(2*pi*(f1-a*t).*t)+sin(2*pi*(f2+a*t).*t)+sin(2*pi*f3*t);
tfrscalo(hilbert(x'));
回车之后会出现第一副图，即信号的能量谱图。

如图2按matlab提示输入分析频率范围的最小值和最大值，还有频率点数（Number of frequency samples）（这个数取512就可以了，它默认值为信号长度，如果信号很长的话，那么频率点数也很长，这样时频图的矩阵就非常大，很难显示，会提示“out of memory"）。

回车后出现图3，图3菜单各项的功能自己可以试试。其中一个便是改变采样频率，tfrscalo(hilbert(x'));这样的调用默认采样频率为1。

点击 changing the sampling frequency 输入采样频率：

点击close后，输入命令colormap(flipud(gray)); 得到白底黑图的图片。

--------------------------------------------------------

按照你说的方法，我怎么只求出一个频率信息啊，zhlong你再帮我看看，辛苦了，而且我怎么看这个图，这么怪啊一半白一半蓝的

--------------------------------------------------

zhlong:

回复 #6 jzc1983 的帖子

这是分析频率范围设置问题。
蓝色底的那一部分就是你分析的范围。
你可以这样设置，
lowwer frequency bound: 0.001
upper frequency bound:0.5
如上面我贴的第二个图。

--------------------------------------------------

问题已经解决了，但是我想把横坐标改成1000ms到2000ms之间，知道在哪个函数里么？我找了几个函数都没有找到相应的语句，谢谢

--------------------------------------------------------

eight:

引用:

原帖由 jzc1983 于 2007-6-13 11:12 发表
问题已经解决了，但是我想把横坐标改成1000ms到2000ms之间，知道在哪个函数里么？我找了几个函数都没有找到相应的语句，谢谢

help xlim

--------------------------------------------------------

好像这样不行啊，xlim([1000 2000])之后图形显示了从1000~2000之间的但是好像只是向右平移了1000个单位，得到的图是空白啊

----------------------------------------------------

zhlong:

回复 #10 jzc1983 的帖子

你的意思是只改变x轴的坐标刻度显示值吗？

-------------------------------------------------

eight:

引用:

原帖由 jzc1983 于 2007-6-13 11:19 发表
好像这样不行啊，xlim([1000 2000])之后图形显示了从1000~2000之间的但是好像只是向右平移了1000个单位，得到的图是空白啊

你的目的不是看 1000-2000 之间的图形吗？这个命令就是了，要不我理解错了，要不你表述得不清楚。

--------------------------------------------------

是我表述不太清楚，我想改Xlabel的值，我用图形框EDIT下面的工具条给改了，谢谢了8同志和zhlong同志，你们是好同志，我顺便问一下，小波多尺度分解后得到的分解系数和原信号的长度不同，用什么命令给重构的相同？小波工具箱中得到的系数是和信号长度相同的
---------------------------------------------------

zhlong:\

help wrcoef

------------------------------------------------

这个好像是根据小波系数重构信号的吧，我想知道小波分解得到的小波系数长度不同，但小波工具箱中的小波系数是相同的且等于信号长度，是怎么弄得？或许我理解的有误区

------------------------------------------------

hitsniq:

不好意思，我用了方法画小波的时频图，是没有能成功。出了以下问题：
??? Error using ==> tfrview
unsupported matlab version. please send an email.

Error in ==> tfrqview at 183
tfrview(tfr,sig,t,method,param,p1,p2,p3);

Error in ==> tfrscalo at 183
tfrqview(tfr,hilbert(real(X)),t,'tfrscalo',wave,N,f);
没能出现蓝色的时频图，我用的Matlab的版本是7.3 R2006b.

--------------------------------------------------------

zhlong:

回复 #17 hitsniq 的帖子

没碰到过你的问题，版面里面还有pengzk版友发的有关小波时频图的程序，可以参考。
我的matlab版本还是7.1

--------------------------------------------------

eight:

引用:

原帖由 hitsniq 于 2007-10-28 21:49 发表
不好意思，我用了方法画小波的时频图，是没有能成功。出了以下问题：
??? Error using ==> tfrview
unsupported matlab version. please send an email.

Error in ==> tfrqview at 183
tfrview(tfr,si ...

估计是作者编写这个函数的时候作了版本限制，自己 edit tfrview，然后搜索 unsupported matlab version. please send an email，最后修改一下吧

原创画小波变换的时频图

文章评论（1条评论）

用户377235 2013-8-27 11:28

用户391779 2008-8-1 15:39

原创 画小波变换的时频图

文章评论（1条评论）

用户377235 2013-8-27 11:28

用户391779 2008-8-1 15:39

原创画小波变换的时频图