请教滑动平均值与滑动均方差的求法

duan2008

请教高手们滑动平均值和滑动方差的求法，上面这一公式中21天的滑动平均和滑动方差平均怎么求? 数据是1981年至2010年的30年的hgt的nc 文件数据， 首先我把这30个nc 文件合并成一个NC文件。 然后应用以下程序求得：'set t 11 10947''set lev 500''define daimean=ave(hgt.1, t-10, t+10)''set sdfwrite D:\Writing\disaster_management\daily_hgt_21days.nc''sdfwrite daimean'滑动平均文件并保存为一个nc文件，然后再应用以下程序：'define centmean=(sum(daimean.1, t=1, t=10937)/21)''set sdfwrite D:\Writing\disaster_management\centered_hgt_21days.nc''sdfwrite centmean'求的21天的滑动平均的nc文件。第一个问题：求得的这个平均文件，我想问一下高手是不是正确的，但是我觉得数值太大，因为如果用前面的X去减这个值，那么永远都是负数，所以感觉应该是错了，但不知道怎么理解这个21-day running mean values.第二个问题：而对于21-day running standard deviations， 感觉又很复杂的去求解，不知道各位高手能不能给点建议，在这个合并的30的nc文件中怎么求的这个值？
不胜感激！

duan2008

呃，高手们帮帮忙啊，怎么计算这两个值，谢谢了。

lqouc

不是很明白你要做什么，所以我说的不一定是你想要的
1.这个滑动平均不是可以用函数做么，例如100个时次
set t 11 90
d tloop(ave(hgt,t-10,t+10))
详细用法参考LASG的手册，或者你可以用while循环做。
2.滑动方差你的文献写的很清楚了啊，就是你做完了这个滑动平均，在某一时刻的观测值与滑动平均之差再平方不就是方差么？这个麻烦么？

duan2008

lqouc 发表于 2013-1-19 19:15

                               
登录/注册后可看大图

不是很明白你要做什么，所以我说的不一定是你想要的
1.这个滑动平均不是可以用函数做么，例如100个时次
s ...

tloop求出来的是这30年中每一天的一个滑动平均值吗？当我用tloop(ave(hgt, t-10,t+10),得到的这30年每一天的一个值，是不是还要平均？如果一平均那怎么不是就是平均值了麽？标准差是一个值吧。

lqouc

duan2008 发表于 2013-1-21 12:38

                               
登录/注册后可看大图

tloop求出来的是这30年中每一天的一个滑动平均值吗？当我用tloop(ave(hgt, t-10,t+10),得到的这30年每一天 ...

我不知道你是要接下来做什么，但是就从数学角度看当然不要再做平均了，否则你还做什么滑动平均啦。标准差是多少个值都无所谓，看你取多少的样本集，你把所有的数据作为一个样本集那就是一个值，你要是把每一年的单独加起来那就是每年的标准差，怎么取的话看你分析的具体问题啦。

duan2008

lqouc 发表于 2013-1-21 13:17

                               
登录/注册后可看大图

我不知道你是要接下来做什么，但是就从数学角度看当然不要再做平均了，否则你还做什么滑动平均啦。标准差 ...

谢谢你，知道怎么做了，只是数据量大，估计要处理很长时间。谢谢高手了。

blui

有没有滑动平均的fortran程序呢？

wpp851118

有没有滑动平均的matlab程序呢？

yangwawa

超前滞后相关+滑动平均
http://bbs.06climate.com/forum.p ... 2&fromuid=21881
(出处: 气象家园)
这个算吗？

shiyueling

沙发