气象统计实习5多元回归方程~仅供参考~

最爱珍珠奶茶

      program main
implicit none
parameter m=13,n=4,k=n-1
real da(m,n),nda(m,n),nnda(k,n),ave(n),ss(n)
real x(k,k),y(k),b(k),ba(k),xa(k,k),e
integer i,j,p
!!!m为时间序列，n为预测因子和预报量的总变量数
!!!da(m,n)为原始数据，最后一列为预报量，nda(m,n)为标准化后的数据,nnda(k,n)为处理后变量矩阵
!!!ave(n)为n个变量的平均值,ss(n)为n个变量的方差
!!!x(k,k)为因子矩阵,y(k)为预报量,b(k)是回归系数,ba和xa是为保存x和y,赋值给ba,xa,再进行后续运算
open(1,file='g:\sx5\1.txt')
      read(1,*) ((da(j,i),i=1,n),j=1,m)
      !!!标准化
do i=1,n
ave(i)=sum(da(1:m,i))/m
enddo
      
do i=1,n
ss(i)=sum((da(1:m,i)-ave(i))**2)/m
enddo

do j=1,n
do i=1,m
nda(i,j)=(da(i,j)-ave(j))/sqrt(ss(j))
enddo;enddo
      !!!求因子矩阵和常数矩阵的协方差
do j=1,k
do p=1,n
      nnda(j,p)=sum(nda(1:m,j)*nda(1:m,p))
enddo;enddo
!!!Gauss Jordan消元法求回归系数b
      do i=1,k
  y(i)=nnda(i,n)
       ba(i)=y(i)
do j=1,k
  x(i,j)=nnda(i,j)
  xa(i,j)=x(i,j)
enddo;enddo
      !!!将变量处理成上三角矩阵
do i=1,k-1
do j=i+1,k
e=xa(j,i)/xa(i,i)
xa(j,1:k)=xa(j,1:k)-xa(i,1:k)*e
ba(j)=ba(j)-ba(i)*e
      enddo;enddo
!!!继续将上面上三角矩阵处理成下三角矩阵
do i=k,2,-1
do j=i-1,1,-1
e=xa(j,i)/xa(i,i)
xa(j,1:i)=xa(j,1:i)-xa(i,1:i)*e
      ba(j)=ba(j)-ba(i)*e
enddo;enddo
      
do i=1,k
b(i)=ba(i)/xa(i,i)
enddo
print*,'b1=',b(1),'b2=',b(2),'b3=',b(3)
end program main

漂洋过海

感谢分享~

maomaomomo

谢谢，谢谢！

若卿宸

楼主有心啦~~~不过这排版看着好眼花啊

笨笨兔

谢谢楼主！正好用上可以好好参考下~

mlittlej

这个程序好像有问题

wallace

感谢分享~

watermelon11111

谢谢楼主！

ThuRonaldo

感谢分享，有收获