求助，Butterworth滤波末端序列变得很小

兰先生

本人之前在对数据进行bw滤波处理的时候，发现自己数据序列的最后几个值都很小，看程序也似乎没有什么问题，想不通问题究竟出在脚本还是什么其他。希望有经验的朋友帮忙解答。

具体过程是使用ERSST的195001-201602的月平均资料先去线性趋势，然后求出距平值，最后做的6-96个月的滤波，各个点所得结果在最后几个都会变小。以东太平洋（0°，120W）为例，201501-201602的SSTA的值见图，最后5个月接近0，然而事实上那会儿是El nino峰值期应该是大值才对。本人换用其他数据和变量也出现相似的结果。

附上脚本。

begin
f1 = addfile("/cygdrive/d/test/data/sst.mnmean.v4.nc","r")
f2 = addfile("/cygdrive/d/test/data/precip.mon.mean.nc","r")
;sst:1854年1月-2016年2月   prec:197901-201602  气候态取1979-2008
t1 = f1->time  ;sst's timeline
x1 =f1->sst
sst = dtrend_msg_n(t1,x1,False,False,0)
copy_VarMeta(x1,sst)
yyyymm1 = cd_calendar(t1,-1)
;计算气候态
clstart1 = ind(yyyymm1.eq.197901)
clend1 = ind(yyyymm1.eq.200812)
clsst = clmMonTLL(sst(clstart1:clend1,:,:))
;计算异常
ostart1 = ind(yyyymm1.eq.195001)
oend1 = ind(yyyymm1.eq.201512)
asst = calcMonAnomTLL(sst(ostart1:oend1,:,:),clsst)
;printVarSummary(asst)
at = new((/794,89,180/),"float")
dstart1 = ind(yyyymm1.eq.201601)
dend1 = ind(yyyymm1.eq.201602)
dt = sst(dstart1:dend1,:,:)
at(0:791,:,:)=asst
do i = 0,1
at(792+i,:,:) = dt(i,:,:)-clsst(i,:,:)
end do
at&time = yyyymm1(ostart1:dend1)
;printVarSummary(at)
  ;bandpass-----------------
  ca    = 96.0        ; band start (longer period)
  cb    = 6.0        ; band end
  fca   = 1.0/ca      ; 'left'  frequency
  fcb   = 1.0/cb      ; 'right' frequency
  opt   = False
  bat = bw_bandpass_filter (at,fca,fcb,opt,0)   ; (ua,fca,fcb,)
  copy_VarMeta(at,bat)
  printVarSummary(bat)

泌泌

就是这么小的，末端数据不能用的。因为滤波都会损失数据，而butterworth滤波虽然得到的序列是完整的，但是两端数据也是它估计出来的是不能用的。具体多少个数不能用取决于你滤波的频段

兰先生

泌泌 发表于 2016-5-19 09:26
就是这么小的，末端数据不能用的。因为滤波都会损失数据，而butterworth滤波虽然得到的序列是完整的，但是 ...

受教了，对BW滤波了解不多，之前只知道有损但两端会保留，没想到还是不能用。多谢~

泌泌

兰先生 发表于 2016-5-19 11:33
受教了，对BW滤波了解不多，之前只知道有损但两端会保留，没想到还是不能用。多谢~

建议把序列选长一点

qq469015280

泌泌 发表于 2016-5-19 09:26
就是这么小的，末端数据不能用的。因为滤波都会损失数据，而butterworth滤波虽然得到的序列是完整的，但是 ...

请移步~~~~~

求助 滤波  急~
http://bbs.06climate.com/forum.p ... 2&fromuid=22714
(出处: 气象家园)

泌泌

qq469015280 发表于 2017-1-3 14:43
请移步~~~~~

求助 滤波  急~

链接不完整

qq469015280

泌泌 发表于 2017-1-3 16:51
链接不完整

http://bbs.06climate.com/forum.p ... 2&fromuid=22714

顾十一

泌泌 发表于 2016-5-19 09:26
就是这么小的，末端数据不能用的。因为滤波都会损失数据，而butterworth滤波虽然得到的序列是完整的，但是 ...

请问不能用的个数应该怎么计算呢，求教

小懒C

不知道这个和阶数有没有关系，李建平老师用的2阶的Butterworth，ncl默认好像是6阶，而且我也觉得最后几个数据不能要