回归程序

mhx_rong

别人写好的，我改进了一下，可以求回归，超前滞后回归。

cz---整型变量，输入参数。超前滞后参数

n——整型变量，输入参数。观测数列长度

x——长度为n的一维实型数组，输入参数。存放自变量x的n个值。

y——长度为n的一维实型数组，输入参数。存放与自变量x相对应的n个y观测值。

a——实型变量，输出参数。回归系数，即线性回归方程的一次项系数。

b——实型变量，输出参数。回归系数，即线性回归方程的常数项。

q——实型变量，输出参数。偏差平方和。

s——实型变量，输出参数。平均标准偏差。

p——实型变量，输出参数。回归平方和。

umax——实型变量，输出参数。最大偏差。

umin——实型变量，输出参数。最小偏差。

u——实型变量，输出参数。偏差的平均值。

!*****************************

! 回归子程序

   subroutine yxhg(cz,n,x,y,a,b,q,s,p,umax,umin,u)

    implicit none

    integer n,i,j,k

    integer nn

    real xv,yv,a,b,s,q,p,umax,umin,u,dxy,dxx,dx

    integer cz

real(4),dimension(n)::x,y

real,dimension(n-abs(cz))::xx,yy

if(cz>=0) then

do i=1,n-cz

    xx(i)=x(i)

    yy(i)=y(i+cz)

end do

else

do i=1,n+cz

    xx(i)=x(i+abs(cz))

    yy(i)=y(i)

end do

end if

nn=n-abs(cz)

xv=0         

yv=0         

do i=1,nn

  xv=xv+xx(i)

  yv=yv+yy(i)

end do

xv=xv/n

yv=yv/n

dxx=0.0

dxy=0.0

do i=1,nn

  q=xx(i)-xv

  dxx=dxx+q*q

  dxy=dxy+q*(yy(i)-yv)

end do

a=dxy/dxx

b=yv-a*xv

q=0

u=0

p=0

umax=-1.0e10

umin=1.0e20

do i=1,nn

  s=a*xx(i)+b

  q=q+(yy(i)-s)**2

  p=p+(s-yv)**2

  dx=abs(yy(i)-s)

  if(dx.gt.umax)then

  umax=dx

  endif

  if(dx.lt.umin)then

  umin=dx

  endif

  u=u+dx/n

end do

s=sqrt(q/n)

end

!*****************************************

mofangbao

有说明看起来会轻松不少 谢谢楼主共享

hubb0119

谢谢楼主分享！

Choo_Choo

多谢楼主，辛苦了

漂洋过海

谢谢了~~

husaian

楼主给力，谢谢

麻花

学习学习~~~~~~

letaotao22

学习学习~~~！

气象小超人

学习学习

linLIN880407

谢谢