一元线性回归

horizon

Program main   
Parameter my=33,mo=12,nx=144,ny=72
Real x(mo,my),y(nx,ny,mo,my),avex(4,my),avey(nx,ny,4,my)
        integer i,j,it,im
        !!!!!!!!!!it=year,im=month!!!!!!!!!!
        real a(nx,ny,mo),b(nx,ny,mo),r(nx,ny,mo),qx(mo),xx(mo),yy(nx,ny,mo),dx(nx,ny,mo),dy(nx,ny,mo),dxy(nx,ny,mo),t(nx,ny,mo)
        !!!!!!!!!xx=ave_x,yy=ave_y,dx,dy,dxy=zhongjianliang!!!!!!!!!!!!!
    open(11,file='c:\2013hg\precip.grd',form='binary',access='direct',recl=nx*ny)
    Open(22,file='c:\2013hg\aoi.grd',form='binary')
Print*,"开始读入数据："
     Do j=1,my
  do i=1,mo
            Read(22),x(i,j)
     End do
  end do
  print*,'aoi is ok'
     print*,x(1,1)
  print*,'precip'
irec=1
       do  it=1,my
    do  im=1,mo
       read(11,rec=irec)((y(i,j,im,it),i=1,nx),j=1,ny)
       irec=irec+1
      enddo
   enddo
print*,'read chang-var is ok!'
     Close(11)
     Close(22)
     Print*,"数据读入完毕！"
   do j=1,my  
  avex(1,j)=(x(3,j)+x(4,j)+x(5,j))/3
  avex(2,j)=(x(6,j)+x(7,j)+x(8,j))/3
  avex(3,j)=(x(9,j)+x(10,j)+x(11,j))/3
     avex(4,j)=(x(12,j)+x(1,j)+x(2,j))/3
   enddo

do j=1,ny
do i=1,nx
do it=1,my
     avey(i,j,1,it)=(y(i,j,3,it)+y(i,j,4,it)+y(i,j,5,it))/3
     avey(i,j,2,it)=(y(i,j,6,it)+y(i,j,7,it)+y(i,j,8,it))/3
     avey(i,j,3,it)=(y(i,j,9,it)+y(i,j,10,it)+y(i,j,11,it))/3
     avey(i,j,4,it)=(y(i,j,12,it)+y(i,j,1,it)+y(i,j,2,it))/3
enddo
enddo
enddo
!????????????????????????????????????????????????????????????????????

s7udY

大哥，这里是matlab的地盘，fortran在隔壁....

kongfeng0824

发错版块了

zengxiaoqing

s7udY 发表于 2014-4-16 12:21
大哥，这里是matlab的地盘，fortran在隔壁....

你这算忽悠人吗

zengxiaoqing

kongfeng0824 发表于 2014-4-16 13:11
发错版块了

你还真的被忽悠了

lqouc

zengxiaoqing 发表于 2014-4-16 13:37
你还真的被忽悠了

仔细看一下吧，这个帖子是被清风移动过来的。这种小事情没必要这么纠结。

天马不行空

正需要！！！

ypliu2015

我也在找线性回归的，(*^__^*) 嘻嘻……

qxylnuist

毕业论文需要，感谢分享

海洋里的鱼

感谢分享，~~~~楼主好人