爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 14689|回复: 22

[源程序] 本人原创用FFT看nc文件的主频主周期的matlab程序

[复制链接]

新浪微博达人勋

发表于 2015-6-28 12:58:57 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 purplejack 于 2015-6-29 12:49 编辑

近期对信号处理有点感兴趣,上学时没学过但研究中常用到,于是自学了一些信号处理的东西,其中基础的当然就是Fourier变换了,不过信号处理通常针对的是一个时间序列,而像我们气象数据乃是时空数据,也就是每个格点上都有一个时间序列,于是根据气象数据的特点自己编了一个小程序,主要目的是“读取nc文件后对每个格点时间序列作Fourier变换得到频谱序列后找出每个格点的主频次主频和主周期次主周期”,特此分享出来,供需要做相关东西的朋友使用,也请大家检验和指正。


gridFFT.m

1.8 KB, 下载次数: 41, 下载积分: 金钱 -5

售价: 1 贡献  [记录]

purplejack原创

评分

参与人数 1金钱 +10 贡献 +5 体力 +50 收起 理由
二爷名声在外 + 10 + 5 + 50

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-6-28 13:10:57 | 显示全部楼层
先说明一点的是:可能由于这个nc文件(为CMIP5中下载的中国两个模式FGOALS和BCC_CSM的个别输出文件)选择得不好,Fourier变换后的频谱序列都是前三个幅值明显太高,考虑到倍频和可能频谱泄漏,我把前三个频率去掉了认为是虚高的频率,也不知道这样对不对?还请高手指正。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-6-28 14:25:25 | 显示全部楼层
请附上数据或图,大家才能判断啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-6-28 15:12:33 | 显示全部楼层
结果图是什么呢?      
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-6-28 21:42:21 | 显示全部楼层
附件中是对主周期分别用matlab和grads画图的结果


mainPeriod_hfls2007.fig

574.26 KB, 下载次数: 36, 下载积分: 金钱 -5

mainPeriod_hfls2007.gmf

934.31 KB, 下载次数: 18, 下载积分: 金钱 -5

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-8-1 12:24:10 | 显示全部楼层
你好,请问下你有对矢量进行过傅里叶变换吗
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-11-28 10:06:01 | 显示全部楼层
谢谢楼主分享
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-4-29 12:15:57 | 显示全部楼层
感谢分享、
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-4-29 13:19:23 | 显示全部楼层
{:5_231:}{:5_231:}{:5_231:}{:5_231:}{:5_231:}{:5_231:}{:5_231:}{:5_231:}{:5_231:}
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2017-4-29 19:21:02 | 显示全部楼层
谢谢楼主分享
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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