eaglezxw 发表于 2007-4-11 20:48

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

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

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

yangzj 发表于 2007-4-11 21:45

怎么都没说是什么问题

eaglezxw 发表于 2007-4-11 22:16

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

songzy41 发表于 2007-4-12 06:41

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

eaglezxw 发表于 2007-4-12 11:32

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

zhlong 发表于 2007-4-13 11:56

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

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

原帖由 eaglezxw 于 2007-4-12 11:32 发表
我知道可以使用tfrwv,像
=tfrwv(sig)
用contour(t,f,tfr)可得到WV分布的等高线图,单我想得到时频分布平面图是不是用imshow(tfr),感觉不对啊
谢谢
在=tfrwv(sig)后,不妨用
imagesc(tfr); axis('xy');
来作图,看能否满足楼主的要求.

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

songzy41 发表于 2007-4-14 09:02 static/image/common/back.gif
在=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);
=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 19:38

石头王石头 发表于 2013-5-14 11:29 static/image/common/back.gif
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);
=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，横轴也以时间为单位。


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

songzy41 发表于 2013-5-14 19:38 static/image/common/back.gif
把程序改为
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=;
x1=;
x2=;
x=x1+x2;
x=x';
tfrwv(hilbert(x),1:length(t));
是不是频率归一化了？
怎么变成实际的频率

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

songzy41 发表于 2013-5-14 19:38 static/image/common/back.gif
把程序改为
N=1024;
n=1:N;



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

谢谢了

chybeyond 发表于 2013-5-14 21:03

石头王石头 发表于 2013-5-14 20:23 static/image/common/back.gif
还有一个问题，怎么把上图中的背景变成透明的。

谢谢了

set(gca,'color','white')
页: [1]
查看完整版本: 求助;如何得到时频分布图