爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
楼主: chenshadow2014

[分享资料] 带通滤波画图

[复制链接]
发表于 2015-7-4 13:54:18 | 显示全部楼层
chenshadow2014 发表于 2015-7-3 22:14
解决了,我是用ncl画的,直接用的滤波函数,气候态函数和方差公式。 你可以试试看。

麻烦发一下NCL相应的函数吧!想看看是啥?
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2015-7-5 14:52:44 | 显示全部楼层
学而不厌 发表于 2015-7-4 13:54
麻烦发一下NCL相应的函数吧!想看看是啥?

请看楼下.....................
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2015-7-5 14:54:48 | 显示全部楼层
你自己看看吧,可能和你的不一样。。。。

bandpass-VAR.ncl

2.15 KB, 下载次数: 25, 下载积分: 金钱 -5

售价: 5 贡献  [记录]

密码修改失败请联系微信:mofangbao
发表于 2015-10-24 20:39:31 | 显示全部楼层
楼主你好~我现在也遇到了感觉和你差不多的问题,要求是对某一年的资料进行时间上的谐波滤波,滤出2-7天,8-20天,30-60天不同频段上的温,压,风场,比较年平均和四季的各不同时间尺度波的分布特征及传播特征。
目前的问题是:
1.这个不同的滤波时间没有说清楚,所以不确定到底是要几天的,这个是有固定的说法还是根据大气的平均态再做定论?
2.请问这个谐波滤波就是做滑动平均的意思么?是的话可以直接用Matlab做的;
3.如果对第一个条件理解为滑动平均的话,滤出的意思是,对所有数据做一个5天的滑动,把高频的信号平滑掉,然后需要再减去什么东西么还是可以直接做图了?谢谢你~
密码修改失败请联系微信:mofangbao
发表于 2015-10-24 20:46:56 | 显示全部楼层
言深深 发表于 2015-6-15 17:36
不知道是不是我理解的问题:
1、多年平均是年资料,3-10天,这种属于天尺度的,怎么滤波?
2、是不是逐日 ...

深深老师你好~我的问题和这位楼楼很相近,也想问问你的想法。先谢谢你了~
要求是对某一年的资料进行时间上的谐波滤波,滤出2-7天,8-20天,30-60天不同频段上的温,压,风场,比较年平均和四季的各不同时间尺度波的分布特征及传播特征。
目前的问题是:
1.这个不同的滤波时间没有说清楚,所以不确定到底是要几天的,这个是有固定的说法还是根据大气的平均态再做定论?
2.请问这个谐波滤波就是做滑动平均的意思么?是的话可以直接用Matlab做的;
3.如果对第一个条件理解为滑动平均的话,滤出的意思是,对所有数据做一个5天的滑动,把高频的信号平滑掉,然后需要再减去什么东西么还是可以直接做图了?谢谢你~
密码修改失败请联系微信:mofangbao
 成长值: 0
发表于 2015-10-25 13:14:59 | 显示全部楼层
小萌萌 发表于 2015-10-24 20:46
深深老师你好~我的问题和这位楼楼很相近,也想问问你的想法。先谢谢你了~
要求是对某一 ...

你说的这个都是天气尺度上面的,30天当然也可以认为是短期气候。
滤波当然是可以的,我没搞懂,你滤波之后打算做什么?考察不同分量的特征?仅仅看特征,那么你绘图就是了。
第二点,平滑相当于滤波,意思是这样的,对于几点平滑的话,相当于滤掉高频的东西,至于你说到底滤掉了几个波,这个是没办法定量去描述的。平滑之后,直接绘图,不需要再减去什么。
最后,至于你用matlab,还是excel,都是无所谓的,能出来结果就行。
密码修改失败请联系微信:mofangbao
发表于 2015-10-25 14:34:10 | 显示全部楼层
本帖最后由 小萌萌 于 2015-10-25 14:37 编辑
言深深 发表于 2015-10-25 13:14
你说的这个都是天气尺度上面的,30天当然也可以认为是短期气候。
滤波当然是可以的,我没搞懂,你滤波之 ...

昂,好的,谢谢您的回复~我是想对物理量场做了不同时间的滤波之后比较年平均和四季的各不同时间尺度波的分布特征及传播特征,所以应该是需要画图就可以了~昂 那我比如要做滤波的话,就可以用Matlab里的函数做了,但是它里面的函数很明确,都没有一个范围,就说我能做3点,5点或者9点的滑动平均,不像是我想得到的一个范围内(2-7天)的天数,想问问您这个怎么弄~还有一个小问题,请问怎么去掉微博和家园账号的关联,就比如不让微博动态在签名栏里显示?
密码修改失败请联系微信:mofangbao
 成长值: 0
发表于 2015-10-26 09:51:38 | 显示全部楼层
小萌萌 发表于 2015-10-25 14:34
昂,好的,谢谢您的回复~我是想对物理量场做了不同时间的滤波之后比较年平均和四季的各不同时间尺度波的 ...

当然不是3-5点平滑,这个只是去掉高频而已,不能等同于滤波。你这个还得用带通滤波才行。

去掉微博,你不要绑定微博就可以了。
密码修改失败请联系微信:mofangbao
发表于 2015-10-27 12:52:46 | 显示全部楼层
言深深 发表于 2015-10-26 09:51
当然不是3-5点平滑,这个只是去掉高频而已,不能等同于滤波。你这个还得用带通滤波才行。

去掉微博, ...

昂原来是要做带通滤波,之前还一直都不明白怎么回事。。现在知道了,原来这种要求时间范围的滤波已经不适用了。谢谢深深老师~微博我解除绑定了,但现在还有不知道怎么回事。这个不是很着急。刚刚在搜grads的时候看见您的一条回复,我一起在这里问了吧,是这样的,我在画图时出现了“Contouring: -2.5e+38 to 3e+38 interval 5e+37”的错误,但是在ctl里面我已经加了“options byteswapped”,而且,这个情况是有时好有时差,想问问这是什么原因呢?
密码修改失败请联系微信:mofangbao
发表于 2016-8-31 21:12:38 | 显示全部楼层
言深深 发表于 2015-10-25 13:14
你说的这个都是天气尺度上面的,30天当然也可以认为是短期气候。
滤波当然是可以的,我没搞懂,你滤波之 ...

请教用ncl滤波滤出8-20天,是要对365天(每一年分开)滤波后再合起来,还是可以直接对1825天(5年)这整时间序列滤???????
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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