一个固定读入4列数据并进行简单计算的代码
此代码固定读入四列数据，读入规则如下：输入文件名为TEST.txt，首行为空行或列代号。数据内容从第二行开始。
数据行为动态数组，可改变。列数通过修改也可以改为任意。
需要注意的是，涉及动态数组在子程序中的调用，则子程序必须用module封装，或在主程序中采用interface进行声明。本例采用了前者所述方法。
感谢fortran群2338021，群友：臭石头雪球，马克钛师傅，小爱的讲解帮助。
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
不是特别理解这个帖子需要讨论的问题是什么？ 这不是讨论问题的，只是写出来，就发上来了，以后有人需要可以直接用
页:
[1]