声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2452|回复: 11

[综合] 求助;如何得到时频分布图

[复制链接]
发表于 2007-4-11 20:48 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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

[ 本帖最后由 eight 于 2007-4-11 21:30 编辑 ]
回复
分享到:

使用道具 举报

发表于 2007-4-11 21:45 | 显示全部楼层
怎么都没说是什么问题
 楼主| 发表于 2007-4-11 22:16 | 显示全部楼层
就是如何得到WVD时频分布的图像
发表于 2007-4-12 06:41 | 显示全部楼层
如果想调用MATLAB的工具箱函数来做,在时频分析工具箱中有tfrwv函数。

评分

1

查看全部评分

 楼主| 发表于 2007-4-12 11:32 | 显示全部楼层
我知道可以使用tfrwv,像
[tfr,t,f]=tfrwv(sig)
用contour(t,f,tfr)可得到WV分布的等高线图,单我想得到时频分布平面图是不是用imshow(tfr),感觉不对啊
谢谢
发表于 2007-4-13 11:56 | 显示全部楼层
单边谱图要用解析信号形式进行分析,即tfrwv(hilbert(sig));
发表于 2007-4-14 09:02 | 显示全部楼层
原帖由 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');
来作图,看能否满足楼主的要求.

评分

1

查看全部评分

发表于 2013-5-14 11:29 | 显示全部楼层

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,怎么回事,还有横坐标是采样点数,怎么变成时间?
不好意思,不知道怎么传图片



发表于 2013-5-14 19:38 | 显示全部楼层
石头王石头 发表于 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,横轴也以时间为单位。


发表于 2013-5-14 20:04 | 显示全部楼层
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));
是不是频率归一化了?
怎么变成实际的频率
发表于 2013-5-14 20:23 | 显示全部楼层
songzy41 发表于 2013-5-14 19:38
把程序改为
N=1024;
n=1:N;

未命名.jpg
还有一个问题,怎么把上图中的背景变成透明的。

谢谢了
发表于 2013-5-14 21:03 | 显示全部楼层
石头王石头 发表于 2013-5-14 20:23
还有一个问题,怎么把上图中的背景变成透明的。

谢谢了

set(gca,'color','white')
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 11:31 , Processed in 0.066335 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表