声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2151|回复: 3

[滤波] 如何画STFT能量,时间,频率三维图

[复制链接]
发表于 2012-9-16 11:06 | 显示全部楼层 |阅读模式

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

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

x
n=length(x);
fs=12000;
t=0:1/fs:(n-1)/fs;
figure(1);
plot(t,x');
title('原始信号');
xlabel('时间 t');
ylabel('幅值 A');
grid on;
nfft=nextpow2(n);
figure(2);
[tfr,t,f]=tfrstft(x,1:nfft,nfft);
contour(t/fs,f(1:length(f)/2)*fs,abs(tfr(1:length(f)/2,:)));
title('短时傅里叶变换布时频图');
xlabel('时间 t');
ylabel('频率 f');
grid on;
figure(3);
mesh(t/fs,f(1:length(f)/2)*fs,abs(tfr(1:length(f)/2,:)));
请问我最后的三维图是能量,时间,频率的三维图吗?
如果不是应该怎样画?
求高手解答,万分感谢!
回复
分享到:

使用道具 举报

发表于 2012-9-30 11:58 | 显示全部楼层
也是,画这种三维图很多种方法,是幅值吧,好像不是能量
发表于 2012-12-24 21:02 | 显示全部楼层
楼上说得对,因为你用的是tfrstft线性时频分析,所以计算出来的tfr的abs是幅值,不是能量。
发表于 2012-12-24 21:39 | 显示全部楼层
个人观点,大家轻拍。
其实,STFT做出来的幅值在一定程度上可以认为是能量在这一频率上的分布情况(不管是加速度信号,还是速度信号,还是位移信号,平方都有能量的意义,另一个方面从数学上说绝对值和平方一样都可以用来衡量距离,也就是说在意义上绝对值和平方是一致的。这是在时间域上的讨论,而STFT也是将时域信号转换到频域上,对应的上述的能量意义在频域上也是成立的)。
补充的是:STFT一般做出来的二维图,但是颜色深浅也是一维,你可以试着旋转你做出来的结果,这样就能看到三维图了。
另外:请把你的结果图贴上来。
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-29 05:13 , Processed in 0.073697 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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