Fortran experiment_6 实验报告

6.1 ex0601.f90

有一个4维数组a(2,2,2,2), 编写程序，根据a的元素在内存中的存储顺序，依次输出其下标。

代码：

program ex0601
implicit none

integer :: i, j, k, l
integer :: a(2,2,2,2)
a = reshape((/(i,i=1,16)/),(/2,2,2,2/))

do l=1,2
    do k=1,2
        do j=1,2
            do i=1,2
                write(*,10) a(i,j,k,l),":(",i,",",j,",",k,",",l,")" 
            end do
        end do
    end do
end do
10 format (i2,a,i1,a,i1,a,i1,a,i1,a)

End

结果：

分析：

这个编程题较为简单，主要在于理解数组在内存中的存储方式，与C语言等一些语言里按行存储不同，Fortran中数组是按列存储的，因此对于一个多维数组a(2,2,2,2)，其下标依次应为(1,1,1,1)，(2,1,1,1)等等。为了方便观察结果，将数组中的元素按下标的顺序编了号（Fortran中数组下标是从1开始的），另外，还要注意输出的格式。

6.2 ex0602.f90

设有10个测试数据，存放在一维数组a中，编写程序计算其方差。

代码：

program ex0602
implicit none

integer :: i, j, k
real ave, sum1, sum2, var
sum1 = 0
sum2 = 0

write(*,*) "comfirm the data"

read(*,*) (a(i),i=1,10)

do j=1, 10
    sum1 = sum1 + a(j)
end do

ave = sum1 / 10

do k=1, 10
    sum2 = sum2 + (a(k)-ave)**2
end do

var = sum2 / 10

write(*,*) var

end

结果：

分析：

这道题要求方差，方差等于每个元素减去平均数的平方再除以个数，因此按部就班的做，首先建立一个数组并读入数据，然后利用一个循环求出平均数，最后再求出方差，注意由于会出现小数所以不要全部用整型变量。