一个固定读入4列数据并进行简单计算的代码

Felix804665 · 发表于 2014-3-4 19:15:20

      此代码固定读入四列数据，读入规则如下：
      输入文件名为TEST.txt，首行为空行或列代号。数据内容从第二行开始。
      数据行为动态数组，可改变。列数通过修改也可以改为任意。
      需要注意的是，涉及动态数组在子程序中的调用，则子程序必须用module封装，或在主程序中采用interface进行声明。本例采用了前者所述方法。
      感谢fortran群2338021，群友：臭石头雪球，马克钛师傅，小爱的讲解帮助。

[Fortran] 纯文本查看 复制代码

module TEST
contains

subroutine mean(A,N,O,P,Q,R)
implicit none

    real :: O,P,Q,R
    integer :: N,i
    real,intent(in) :: A(:,:)

    !记录ABCD的平均数，分别对应OPQR    
    O=0.0
    P=0.0
    Q=0.0
    R=0.0

    do i=1,N
        O=O+A(i,1)
        P=P+A(i,2)
        Q=Q+A(i,3)
        R=R+A(i,4)
    enddo

    O=O/real(N)
    P=P/real(N)
    Q=Q/real(N)
    R=R/real(N)

return
end subroutine mean

subroutine FANGCHA(A,N,O,P,Q,R,K,L,M,H)
implicit none

    real,intent(in) :: A(:,:)
    integer :: i,N
    real :: O,P,Q,R,K,L,M,H
    real :: B(N),C(N),D(N),E(N)

    K=0.0
    L=0.0
    M=0.0
    H=0.0

    do i=1,N
        B(i)=(A(i,1)-O)**2
        C(i)=(A(i,2)-P)**2
        D(i)=(A(i,3)-Q)**2
        E(i)=(A(i,4)-R)**2
    enddo        

    do i=1,N
        K=K+B(i)
        L=L+C(i)
        M=M+D(i)
        H=H+E(i)
    enddo

        K=(K)/real(N)
        L=(L)/real(N)
        M=(M)/real(N)
        H=(H)/real(N)

return
end subroutine FANGCHA

subroutine BIAOZHUNCHA(K,L,M,H,bzc1,bzc2,bzc3,bzc4)
implicit none

    real :: K,L,M,H,bzc1,bzc2,bzc3,bzc4

    bzc1=sqrt(K)
    bzc2=sqrt(L)
    bzc3=sqrt(M)
    bzc4=sqrt(H)

return
end subroutine BIAOZHUNCHA

subroutine BZCS(bzc1,bzc2,bzc3,bzc4,bzc5,bzc6,bzc7,bzc8,bzc9,bzc10,bzc11,bzc12)
implicit none

    real :: bzc1,bzc2,bzc3,bzc4,bzc5,bzc6,bzc7,bzc8,bzc9,bzc10,bzc11,bzc12

    bzc5=bzc1/3.0
    bzc6=bzc2/3.0
    bzc7=bzc3/3.0
    bzc8=bzc4/3.0

    bzc9=bzc1/4.0
    bzc10=bzc2/4.0
    bzc11=bzc3/4.0
    bzc12=bzc4/4.0

return
end subroutine BZCS

end module TEST

program main
use TEST
implicit none

    real,allocatable :: A(:,:)
    integer :: i,j,N
    real :: O,P,Q,R,K,L,M,H,bzc1,bzc2,bzc3,bzc4
    real :: bzc5,bzc6,bzc7,bzc8,bzc9,bzc10,bzc11,bzc12

    write(*,*) "输入ABCD每列元素个数(注意元素个数必须相等):"
    read(*,*) N

    j=4

    open(1,file='TEST.txt',status='old')
    open(2,file='TEST01.txt')

    allocate (A(N,j))

    read(1,*)
    write(2,*) "A,B,C,D分别为:"
    write(2,*)

    do i=1,N
        read(1,*) (A(i,j),j=1,4)
        write(2,*) (A(i,j),j=1,4)
    enddo

    write(2,*)

!计算平均值    
    call mean(A,N,O,P,Q,R)

!计算方差
    call FANGCHA(A,N,O,P,Q,R,K,L,M,H)

!计算标准差
    call BIAOZHUNCHA(K,L,M,H,bzc1,bzc2,bzc3,bzc4)

!标准差再计算
    call BZCS(bzc1,bzc2,bzc3,bzc4,bzc5,bzc6,bzc7,bzc8,bzc9,bzc10,bzc11,bzc12)

    write(2,*) "A,B,C,D的平均数分别为:"
        write(2,*)
    write(2,*) O,P,Q,R
        write(2,*)
    write(2,*) "A,B,C,D的方差分别为:"
        write(2,*)
    write(2,*) K,L,M,H
        write(2,*)
    write(2,*) "A,B,C,D的标准差为:"
        write(2,*) 
    write(2,*) bzc1,bzc2,bzc3,bzc4
        write(2,*)
    write(2,*) "A,B,C,D标准差的三分之一:"
        write(2,*)
    write(2,*) bzc5,bzc6,bzc7,bzc8
        write(2,*)
    write(2,*) "A,B,C,D标准差的四分之一:"
        write(2,*)
    write(2,*) bzc9,bzc10,bzc11,bzc12

    DEALLOCATE(A)

pause
end