关于去掉序列月振荡——兼答坛友(若兰)求助

言深深

@若兰在帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=11812中提出来这个问题：
如何去点月平均序列中的{气候态}【错误表达：年内振荡，也就是季节振荡（感谢@游子[/url]指点）】，这个程序本不是太难，真好我曾经也写过这样的小程序，顺便编辑一下贴上来。首先用下面一张图来表示上面这句话的意思：

其中红色时间序列是原序列，存在12这样的周期性在里面，我们要去掉这样的趋势，于是想到的办法是“距平值要求是1月的减去1月的多年平均，2月的减2月的以此类推(见原帖)”，这样我们用程序
parameter(m=100)
real x(m)
open(1,file="test.txt")
do i=1,m
read(1,*) x(i)
enddo
close(1)
do mon=0,11
  s=0
  nc=0
  do i=1,m
  if(mod(i,12).eq.mon)then
  s=s+x(i)
  nc=nc+1
  endif
  enddo
   do i=1,m
   if(mod(i,12).eq.mon)then
   x(i)=x(i)-s/nc
   endif
   enddo
enddo
open(2,file="result by 言深深.txt")
do i=1,m
write(2,*)x(i)
enddo
close(2)
end
DeTrend.f90 (381 Bytes, 下载次数: 33)

2012-12-21 20:37 上传

点击文件名下载附件 处理一下就得到蓝色圈圈表示的序列。
附上我的测试数据 test.txt (324 Bytes, 下载次数: 29)

2012-12-21 20:37 上传

点击文件名下载附件

尽头的尽头

深深辛苦了~~

言深深

尽头的尽头 发表于 2012-12-21 21:11

                               
登录/注册后可看大图

深深辛苦了~~

汗，一个小程序，正好有就顺手发上来了。不值这么多钱啦···不过还是很开心哈···哈哈哈哈

river

程序猿的帖子必须顶！

游子

这个不叫季节震荡，叫做气候态。季节震荡不是这样定义的

若兰

万分感谢深深指导！！！

言深深

若兰 发表于 2012-12-22 15:27

                               
登录/注册后可看大图

万分感谢深深指导！！！

不客气，欢迎论坛多多讨论交流哈

157041245

这个用来去处季节趋势的吧？

言深深

157041245 发表于 2012-12-26 20:47

                               
登录/注册后可看大图

这个用来去处季节趋势的吧？

差不多就是这个意思呢······我不知道这个要怎么表达才好，呵呵

ayzqs

看了你好多帖子，真是太赞了！不简单，但易懂