短期气候预测实习报告二(2)

enddo

enddo

enddo

do it=1,61

do j=1,21

do i=1,41

f(21*(i-1)+j,it)=h(i,j,it)

enddo

enddo

enddo

2. 计算协方差矩阵

3. 用Jacobi方法计算协方差矩阵的特征值与特征向量

4. 将特征值从大到小排列

5. 计算特征向量的时间系数

6. 计算每个特征向量的方差贡献

7. 结果输出

输出的eof.txt：

EIGENVALUE AND ANALYSIS ERROR

H LAMDA SLAMDA PH SPH

1 11844.33105 11844.33105 0.22552 0.22552

2 11124.51270 22968.84375 0.21181 0.43733

3 4 5 6 9327.01270 32295.85547 5934.51562 38230.37109 3124.63110 41355.00391 2463.64844 43818.65234 0.17759 0.61492 0.11299 0.72791 0.05949 0.78740 0.04691 0.83431

7 2322.31592 46140.96875 0.04422 0.87853

8 1282.84216 47423.81250 0.02443 0.90295

实习结果：

第一空间场和时间序列：