- 积分
- 48
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-9-4
- 最后登录
- 1970-1-1
![[如果我不开拖拉机你还会爱我吗] 粉丝数:120 微博数:157 新浪微博达人勋](source/plugin/sina_login/img/light.png)
|
10金钱
对500hPa逐日位势高度场进行10-25天带通滤波,输出看了一下发现数值变得好小,只有±50左右,这是为什么呀?我理解的脚本里wgt_runave_n只是做滑动平均,得到的应该也是5500左右的数值,附脚本
;-----------------------------------设置滤波器------------------------------
nWgt = 21 ; loose 10d each end
ihp = 2 ; band pass
sigma = 1.0 ; Lanczos sigma
fca = 1./25. ; start freq
fcb = 1./10. ; last freq
wgt = filwgts_lanczos (nWgt, ihp, fca, fcb, sigma ) ;得到滤波权重
printVarSummary(wgt)
BP=new((/NY,dsum,nlat,nlon/),"float") ;dsum是120天,对120天滤波
do YY=1979,2021
pc1=unew({YY},:,:,:)
bp = wgt_runave_n (pc1, wgt, 0,0)
BP(YY-1979,:,:,:)=bp
end do
copy_VarMeta(unew,BP)
printVarSummary(BP)
print(BP({1979},:,{40},{120}))
|
最佳答案
查看完整内容
因为滤波会先减去平均值,也就是虽然你输入是原始场,但是实际上是距平进行滤波,50多正常的
|