- 积分
- 102
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-8-17
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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
|
评分
-
查看全部评分
|