滤波出现问题，求助

曦尘

我要翻一下C. Duchon（1979）的《Lanczos Filtering in One and Two Dimensions》才能给你答案了，nWgt的长度应该是有个参考范围的，你选的值越大，虽然滤波结果更准确，但是首尾损失的数据就更多，上面官网的例子说明里有反映，用的时候要权衡一下。具体参数，我晚点回复你~

subtropical

曦尘 发表于 2016-3-11 15:12
Lanczos滤波程序有个地方是值得注意的，比如说你的滤波窗口nWgt选择121时步时，你的时间序列前后各有0.5*（ ...

我想做的就是在130年里把elnino年份的春季海温异常年份提取出来看看，但是刚开始做的时候发现海温异常有年代际和线性趋势，然后我就去除线性趋势，再用滤波滤去年代际信号，您觉得我应该用哪种滤波好一点，我已开始用的是十一点滑动平均，但是有人说这个方法不好，所以我就用滤波了。。。麻烦您指教，非常感谢

subtropical

曦尘 发表于 2016-3-11 15:32
我要翻一下C. Duchon（1979）的《Lanczos Filtering in One and Two Dimensions》才能给你答案了，nWgt的长 ...

非常非常感谢你，大侠，我也好好学学官网例子，把这些弄明白{:eb511:}{:eb511:}

曦尘

nWgt的选取应该是参考2*n+1，官网里是低通滤波两年（24个月），所以最后选了2*24+1=49，所以你要是滤掉10年以上的成分，建议你用nWgt=2*10+1=21。另外一个，就是建议你快速傅里叶分解FTT和Lanczos滤波都用一下，互相对比，快速傅里叶FTT首尾两端都是有数据的，不会影响数据的长度，可以参考官网的例子4。其实十一点滑动平均也能做，但是响应函数估计没有上述两种方法好，滤波后信息丢失会厉害一点。响应函数就是你给的例子里右边的图，越接近那个折线图的形状就越完美。

subtropical

曦尘 发表于 2016-3-11 16:06
nWgt的选取应该是参考2*n+1，官网里是低通滤波两年（24个月），所以最后选了2*24+1=49，所以你要是滤掉10年 ...

好的，非常感谢，您说的对。我刚试了nwgt选取不同，出来的结果不太一样，好像是nwgt大一些的信号更强一点如图：左边是nwgt=3.右边是nwgt=21，但是因为我提取的年份比较靠前和靠后，导致那些只参与计算的年份画不出来，那我试试您说的FTT，另外一个问题如果我用带通滤波的话滤出我需要的年际变化，那nWgt怎么设置呢？还是21么？

曦尘

嗯，如果是过滤10年，还是21吧。但是这里可能还有个小问题，我们说的年际和年代际的划分不是以10年为标准，而是以8年来作为阈值划分的，这个问题还需要斟酌一下。还有，你可以按照原来的思路，先低通滤波后，再用原序列去减，估计和高通滤波的结果应该是一样的。

右手边的小星星

我最近也在做滤波，想把高频和低频信号分离开来，分别试了一下高通滤波，以及原序列减掉低通滤波，得到的结果并不一样，差别还比较明显，楼主有对比过吗？

subtropical

曦尘 发表于 2016-3-11 17:15
嗯，如果是过滤10年，还是21吧。但是这里可能还有个小问题，我们说的年际和年代际的划分不是以10年为标准， ...

嗯是的，我做了发现两个是一样的

subtropical

右手边的小星星 发表于 2016-3-11 17:47
我最近也在做滤波，想把高频和低频信号分离开来，分别试了一下高通滤波，以及原序列减掉低通滤波，得到的结 ...

您好，我做了发现我的结果是一样的。。。

subtropical

曦尘 发表于 2016-3-11 16:06
nWgt的选取应该是参考2*n+1，官网里是低通滤波两年（24个月），所以最后选了2*24+1=49，所以你要是滤掉10年 ...

您好，我又遇到问题了，在算FFT时候，老是出错，想请教您一下 真的谢谢您了