爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4706|回复: 3

【短时傅里叶变换】【加窗傅里叶变换】【gabor变换】的MATLAB实现问题

[复制链接]

新浪微博达人勋

发表于 2018-12-15 19:37:00 | 显示全部楼层 |阅读模式
MATLAB
问题截图: -
问题概况: lz 最近在做加窗傅里叶变换(也称短时傅里叶变换)的问题。在MATLAB和网页上找了半天的资料,找到spectrogram这个相关的函数,无奈网上相关资料太少,论坛上也没见到有成熟的代码。比如对于一组连续时间的水位,对它做短时傅里叶变换,画出图来是:横坐标为时间,纵坐标是频率的二维图,就像小波分析的图一样。


lz还有一个问题,【gabor变换】是【短时傅里叶变换】的一种,怎么使用MATLAB做gabor变换
我看过提问的智慧: 看过
自己思考时长(天): 1
系统平台:

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 13top 于 2018-12-15 20:00 编辑

【短时傅里叶变换】【加窗傅里叶变换】【gabor变换】的MATLAB实现问题
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-12-15 19:39:07 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-12-15 19:49:55 | 显示全部楼层
%%------------------------------------------------------------------------
%%功能:使用MATLAB对特定信号做短时傅里叶分析
%%------------------------------------------------------------------------
fs=1024;
t=0:1/fs:2;
x=cos(2*pi*50*t).*(t>=0&t<1)+0.2*cos(2*pi*150*t).*(t>=1&t<2);
f=-500:0.1:0;
a=size(f)
subplot(2,1,1);
spectrogram(x,kaiser(60,5),20,f,fs,'yaxis');
title('Kaiser窗,时域宽度为60');
%%或者用下面的语句绘图
%[s,f,t,p]=spectrogram(x,kaiser(60,5),20,f,fs,'yaxis');
%surf(t,f,10*log10(abs(p)),'EdgeColor','none');
%axis xy;
%axis tight;
%colormap(jet);
%view(0,90);
%xlabel('Time');
%ylabel('Frequency (Hz)');
%%以上语句与用无输出参数的spectrogram相同
subplot(2,1,2);
spectrogram(x,kaiser(260,5),20,f,fs,'yaxis');
title('Kaiser窗,时域宽度为260');

http://www.ilovematlab.cn/thread-338529-1-1.html 刚才在ilovematlab找的文件,不知道对不对
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-12-15 20:07:33 | 显示全部楼层
希望有童靴解答
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

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