求助;如何得到时频分布图

eaglezxw

我做时频分布,但是想得到WVD的时频分布时碰到问题,请各位高手指点!

[ 本帖最后由 eight 于 2007-4-11 21:30 编辑 ]

yangzj

怎么都没说是什么问题

eaglezxw

就是如何得到WVD时频分布的图像

songzy41

如果想调用MATLAB的工具箱函数来做，在时频分析工具箱中有tfrwv函数。

eaglezxw

我知道可以使用tfrwv,像
[tfr,t,f]=tfrwv(sig)
用contour(t,f,tfr)可得到WV分布的等高线图,单我想得到时频分布平面图是不是用imshow(tfr),感觉不对啊
谢谢

zhlong

单边谱图要用解析信号形式进行分析，即tfrwv(hilbert(sig));

songzy41

原帖由 eaglezxw 于 2007-4-12 11:32 发表
我知道可以使用tfrwv,像
=tfrwv(sig)
用contour(t,f,tfr)可得到WV分布的等高线图,单我想得到时频分布平面图是不是用imshow(tfr),感觉不对啊
谢谢

在[tfr,t,f]=tfrwv(sig)后,不妨用
imagesc(tfr); axis('xy');
来作图,看能否满足楼主的要求.

石头王石头

songzy41 发表于 2007-4-14 09:02

                               
登录/注册后可看大图

在[tfr,t,f]=tfrwv(sig)后,不妨用
imagesc(tfr); axis('xy');
来作图,看能否满足楼主的要求.

N=1024;
n=1:N;
fs=4000;
t=n/fs;
fx=1000;
x=cos(2*pi*fx*t);
y=hilbert(x);
[tfr,t1,f]=tfrwv(y');
f=fs*(n-1)/N/2;%%在WVD变换中其分辨率不是fs/N（不同于FFT的变换），而变为fs/（2*N）
mesh(t,f,abs(tfr));view(-75,30);
figure;
imagesc(tfr);
axis('xy');
按你说的做出来的频率是500，而不是1000，怎么回事，还有横坐标是采样点数，怎么变成时间？
不好意思，不知道怎么传图片

songzy41

石头王石头 发表于 2013-5-14 11:29

                               
登录/注册后可看大图

N=1024;
n=1:N;
fs=4000;

把程序改为
N=1024;
n=1:N;
fs=4000;
t=n/fs;
fx=1000;
x=cos(2*pi*fx*t);
y=hilbert(x);
[tfr,t1,f]=tfrwv(y');
%f=fs*(n-1)/N/2;%%在WVD变换中其分辨率不是fs/N（不同于FFT的变换），而变为fs/（2*N）
mesh(t,f*fs,abs(tfr));view(-75,30);
figure;
imagesc(t1/fs,f*fs,tfr);
axis('xy');
得图中频率为1000Hz，横轴也以时间为单位。

石头王石头

songzy41 发表于 2013-5-14 19:38

                               
登录/注册后可看大图

把程序改为
N=1024;
n=1:N;

非常感谢，再帮忙看看这个
clear;
clc;
T=1,fs=1024;
%-------------------频率突变信号------------------------
t1=0:1/fs:T/2-1/fs;
t2=T/2:1/fs:T-1/fs;
t=[t1,t2];
x1=[sin(2*pi*50*t1),zeros(1,length(t2))];
x2=[zeros(1,length(t1)),sin(2*pi*100*t1)];
x=x1+x2;
x=x';
tfrwv(hilbert(x),1:length(t));
是不是频率归一化了？
怎么变成实际的频率

石头王石头

songzy41 发表于 2013-5-14 19:38

                               
登录/注册后可看大图

把程序改为
N=1024;
n=1:N;

还有一个问题，怎么把上图中的背景变成透明的。

谢谢了

chybeyond

石头王石头 发表于 2013-5-14 20:23

                               
登录/注册后可看大图

还有一个问题，怎么把上图中的背景变成透明的。

谢谢了

set(gca,'color','white')