爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 42904|回复: 51

[源代码] 分享一个九点二次平滑的程序(相当于低通滤波)

  [复制链接]

新浪微博达人勋

发表于 2013-4-12 10:24:15 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
九点二次平滑(不仅起到低通滤波器的作用,并且克服了滑动平均削弱过多波幅的缺点),程序很简单,希望能对大家有帮助

ps附件上传不上,我只能复制粘贴了
implicit none
integer m,n,yr
parameter(m=61,n=2)
integer i,j
real a(m,n),c(m)

open(1,file='E:\2\160\jpwinall.txt')
open(2,file='E:\2\160\jpwinall9h.txt')
do i=1,m
  read(1,*) (a(i,j),j=1,n)
end do


do i=5,m-4
yr=1950
c(i)=(-21*a(i-4,2)+14*a(i-3,2)+39*a(i-2,2)+54*a(i-1,2)+59*a(i,2)+54*a(i+1,2)+39*a(i+2,2)+14*a(i+3,2)-21*a(i+4,2))/231
yr=yr+i
write(2,"(i5,f8.1)") yr,c(i)
enddo

end

评分

参与人数 2金钱 +15 贡献 +2 收起 理由
lgg00828 + 5
mofangbao + 10 + 2

查看全部评分

本帖被以下淘专辑推荐:

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-4-12 10:42:33 | 显示全部楼层
程序不错,简洁支持一下!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-4-12 11:23:34 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-7-20 19:30:04 | 显示全部楼层
这个的原理有没有出处?你采用的是气象统计上的非平均平滑?还是哪里的原理?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-7-20 18:27:11 | 显示全部楼层
程序很实用,赞一个~~~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-5-8 10:30:50 | 显示全部楼层
楼主真个好人啊,很有用的东西
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-5-10 09:22:20 | 显示全部楼层
简洁但实用,谢谢.
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-5-10 15:14:13 | 显示全部楼层
多谢分享!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-5-24 15:07:54 | 显示全部楼层
LZ给力!真心赞一个~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-6-21 09:41:58 | 显示全部楼层
真好,大谢!
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表