实验报告二 经验正交分解

气象统计分析与预报方法 课程实验报告

实验名称 实验二 经验正交函数分解

系 别 大气科学 姓 名 学 号 班 级

应气101

实验地点 机房 实验日期 11月13日

评 分

指导老师

肖国杰

同组其他成员

一、实验内容（含实验原理介绍）：实验所提供的资料为NCEP/NCAR 59年（1948年-2006年）逐年1～12月的

850hPa 高度场资料，资料范围为（90

N

-90S ，0E -360E ），网格距为2.5*2.5，纬向格点数为144，经向格点

数为73。资料为NC 格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放，注意读取方式以及记录长度。

对（0N -90N ，60E -120W ）850hPa 高度场进行经验正交展开（EOF.FOR ），输出分析主要参数指标；绘制环流型图和相应的时间系数序列图，并加以分析。 本实验运用EOF 方法：

EOF （经验正交函数分解）是针对气象要素场进行的，其基本原理是把包含p 个空间点

（变量）的场随时间变化进行分解。设抽取样本容量为n 的资料.则场中任一空间点i 和任一时间点j 的距平观测值

ij x 可看成由p 个空间函数ik v 和时间函数kj y (k=1,2,…,p)的线性组合,表示成

11221

p

ij ik

kj i j i j ip pj k x v

y v y v y v y ==

=+++∑

EOF 功能是从一个气象场多次观测资料中识别出主要空间型及其时间演变规律。 EOF 展开就是将气象变量场分解为空间函数（V ）和时间函数（T ）两部分的乘积之和： X=VT 。

应用步骤：

资料预处理(距平或标准化处理)

计算协方差矩阵、用Jacobi 方法或迭代法计算协方差矩阵的特征值与特征向量、将特征值从大到小排列、计算特征向量的时间系数、计算每个特征向量的方差贡献、结果输出

二、实验目的：经验正交函数分解（EOF）是统计天气分析中气象要素场最基础的研究模型，是必须理解和掌握的方法之一，是后续课程中许多气象要素场的计算结果的理解的基础理论，也是毕业设计和论文中的基本分析方法。该方法用个数较少的几个空间分布模态来描述环流形势，而且基本涵盖环流场的信息，既能作为天气分析模型，其方法的延拓又能作为天气预报模型，在实际工作中也有极强的实用意义。通过该实验，深刻理解气象要素场的统计模型的意义，掌握气象要素场分析的基本方法，为实际预报业务和科研工作打下一定的基础。

三、涉及实验的相关情况介绍（包含使用软件或实验设备等情况）：

本实验要求运用Fortran6.5编译软件以及grads 计算机一台。

四、实验结果（含程序、数据记录及分析和实验总结等，可附页）：

要求编写主程序，其中包括资料读入，范围截取，子程序调用。注意：EOF的资料输入，时间场一维，空间场一维。

*********************（附程序，对关键部分标志出）**********************

EOF程序

C**********************************************************************

C *

C PROGRAM NOTES *

C *

C THIS PROGRAM USES EOF TO ANALYSIS TIME SERIES *

C OF METEOROLOGICAL FIEL

D *

C *

C**********************************************************************

C *

C ******** Parameter Table ********* *

C *

C Mt===>LENTH OF TIME SERIES *

C N ===>NUMBER OF GRID-POINTS ( or STATIONS ) *

C KS=-1, SELF; KS=0, DEPATURE; KS=1, STANDERDLIZE

D DEPATUR

E *

C KV = NUMBER OF EIGENVALUES WILL BE OUTPUT *

C KVT = NUMBER OF EIGENVECTORS AN

D TIM

E SERIES WILL BE OUTPUT *

C MNH = Minimum(Mt,N) *

C EGVT===>EIGENVECTORS, ECOF===>TIME COEFFICIENTS FOR EGVT *

C ER(KV,1)====>LAMDA; LAMDA===>EIGENVALUE *

C ER(KV,2)====>ACCUMULATE LAMDA *

C ER(KV,3)====>THE SUM OF COMPONENTS VECTORS PROJECTE

D ONTO *

C EIGENVACTOR. *

C ER(KV,4)====>ACCUMULATE ER(KV,3) *

C *

C**********************************************************************

PARAMETER(N=73*37, MT=58, MNH=58)

PARAMETER(KS=1, KV=10, KVT=10)

REAL F(N,MT),AVF(N),DF(N),ER(MNH,4)

REAL A(MNH,MNH),S(MNH,MNH),V(MNH)

c**************************************************************************

c INFN－输入数据文件名；OUTERA－输出特征值及方差贡献、累积方差贡献的文件名（文本）;

c OUTTC1－输出时间系数文件（文本）；OUTTC2－输出时间系数文件（二进制）；

c OUTTEVT－输出特征向量文件（二进制）；

c**************************************************************************

CHARACTER*50 INFN,OUTERA,OUTTC1,OUTTC2,OUTEVT

DATA INFN/'hgt8501948-2005july.grd'/

DATA OUTERA/'hgt_XT03ER3.DAT'/

DATA OUTTC1/'hgt_XT03TC13.DAT'/

DATA OUTTC2/'hgt_XT03TC23.DAT'/

DATA OUTEVT/'hgt_XT03VT3.DAT'/

C---------------- Read ORIGINAL DATA ----------------------------