滑动平均小程序

无敌三脚猫

program huadong
implicit none
integer,parameter::m=2000,n=47           !n为滑动步数
real ser(m),have(m)
real a,b,c
integer i,j,k,len1,len2                 !len1为原序列长,len2为滑动后的序列长
open(1,file='',status='old')
do i=1,m
read(1,*)ser(i)
if(eof(1)) goto 100
end do
100 continue
len1=i
len2=len1-n
close(1)
do j=1,len2
a=0
do i=j,j+n
  a=a+ser(i)
end do
have(j)=a/(n+1)
print*,j,have(j)
end do
open(2,file='',status='replace')
write(2,200)(have(j),j=1,len2)
close(2)
200 format(f5.1)
end

天道酬勤

是不是经过滑动之后，会缩减数据长度啊~

mofangbao

自己动手写写这种程序挺重要的

Aires

天道酬勤 发表于 2013-6-28 14:15

                               
登录/注册后可看大图

是不是经过滑动之后，会缩减数据长度啊~

滑动主要用于体现趋势，低通滤波，两端数据会根据滑动平均的长度有所减少，但不影响分析

善人/jw

天道酬勤 发表于 2013-6-28 14:15

                               
登录/注册后可看大图

是不是经过滑动之后，会缩减数据长度啊~

是有公式可以把两端缺少的数据补齐的

残月

善人/jw 发表于 2013-6-28 14:58

                               
登录/注册后可看大图

是有公式可以把两端缺少的数据补齐的

是什么公式啊？

随缘

直接贴上来，楼主好人啊

天道酬勤

Aires 发表于 2013-6-28 14:54

                               
登录/注册后可看大图

滑动主要用于体现趋势，低通滤波，两端数据会根据滑动平均的长度有所减少，但不影响分析

是不影响分析，但如果想要低通滤波的结果来进行运算，就有影响了。

天道酬勤

善人/jw 发表于 2013-6-28 14:58

                               
登录/注册后可看大图

是有公式可以把两端缺少的数据补齐的

同求！贴上来更好。。。

善人/jw

残月 发表于 2013-6-28 15:34

                               
登录/注册后可看大图

是什么公式啊？

我也记不清了，只是本科上课的时候，老师讲过，可以看看魏凤英老师的书，应该会有