关于滤波

subtropical

易小凯 发表于 2015-9-17 14:35
是的，要检验A是否显著。可以参考统计学的相关分析。

大侠，请问有木有高通滤波的程序，脚本。。。或者函数什么的，想做不知道怎么做。。。。

易小凯

subtropical 发表于 2016-3-6 20:50
大侠，请问有木有高通滤波的程序，脚本。。。或者函数什么的，想做不知道怎么做。。。。

一般都是低通滤波，高频波一般都是干扰信号。想要得到高频波，用低通滤波滤出低频波，然后去掉低频部分的就是高频部分了。

subtropical

易小凯 发表于 2016-3-7 07:43
一般都是低通滤波，高频波一般都是干扰信号。想要得到高频波，用低通滤波滤出低频波，然后去掉低频部分的 ...

嗯，谢谢你，我昨晚在官网上搜了然后滤低频波，再把低频减掉，但是我不知道哪儿出了问题，图只出来了一半，可以麻烦大侠帮我看下么？也没提示错误，所以我也不知道怎么回事。。

脚本的重要部分如下，再次感谢您了，您有时间的话帮忙看一下吧

易小凯

33行的wgt_runave用错了，你这行表示的是在纬向上进行滤波，前面一段自然就没有了。
正确的应该使用函数wgt_runave_n_Wrap, s3=wgt_runave_n_Wrap(s2,wgts,0,0)

易小凯

subtropical 发表于 2016-3-7 09:29
嗯，谢谢你，我昨晚在官网上搜了然后滤低频波，再把低频减掉，但是我不知道哪儿出了问题，图只出来了一半 ...

33行的wgt_runave用错了，你这行表示的是在纬向上进行滤波，前面一段自然就没有了。
正确的应该使用函数wgt_runave_n_Wrap, s3=wgt_runave_n_Wrap(s2,wgts,0,0)

subtropical

易小凯 发表于 2016-3-8 08:00
33行的wgt_runave用错了，你这行表示的是在纬向上进行滤波，前面一段自然就没有了。
正确的应该使用函数 ...

非常感谢您的回复，但是我的ncl不知道是不是版本低了，输入s3=wgt_runave_n_Wrap(s2,wgts,0,0)这个没用
那应该怎么解决呢？谢谢您

subtropical

易小凯 发表于 2016-3-8 08:00
33行的wgt_runave用错了，你这行表示的是在纬向上进行滤波，前面一段自然就没有了。
正确的应该使用函数 ...

我明白您的意思了，我用wgt_runave_leftdim表示对时间的滤波，函数是起作用的

但是结果却是这样的。。。

图就更奇怪了，竟然有一张是对的，其他的全是空白，麻烦大神帮忙看下，这是为什么啊？

易小凯

subtropical 发表于 2016-3-8 17:28
非常感谢您的回复，但是我的ncl不知道是不是版本低了，输入s3=wgt_runave_n_Wrap(s2,wgts,0,0)这个没用
...

那你就自行变换一下维度，s3=wgt_runave(s2(lat|:,lon|:,time|:)
然后在变回来就行

baggio

{:lxm_20:}{:lxm_20:}{:lxm_20:}{:lxm_20:}{:lxm_20:}

subtropical

易小凯 发表于 2016-3-9 08:05
那你就自行变换一下维度，s3=wgt_runave(s2(lat|:,lon|:,time|:)
然后在变回来就行

非常感谢您的回复，我其实一直都有一个疑问就是怎么把原来改了顺序的变量参数再变回来，求大神指教。谢谢您