资源描述：
《实验报告二 经验正交分解》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、气象统计分析与预报方法课程实验报告实验名称实验二经验正交函数分解系别大气科学姓名学号班级应气101实验地点机房实验日期11月13日评分指导老师肖国杰同组其他成员一、实验内容（含实验原理介绍）：实验所提供的资料为NCEP/NCAR59年（1948年-2006年）逐年1～12月的850hPa高度场资料，资料范围为（-，-），网格距为2.5*2.5，纬向格点数为144，经向格点数为73。资料为NC格式，资料从南到北、自西向东排列，每月为一个记录，按年逐月排放，注意读取方式以及记录长度。对（-，-）850hPa高度场进行经验正交展开（EOF.FOR）
2、，输出分析主要参数指标；绘制环流型图和相应的时间系数序列图，并加以分析。本实验运用EOF方法：EOF（经验正交函数分解）是针对气象要素场进行的，其基本原理是把包含p个空间点（变量）的场随时间变化进行分解。设抽取样本容量为n的资料.则场中任一空间点i和任一时间点j的距平观测值可看成由p个空间函数和时间函数(k=1,2,…,p)的线性组合,表示成EOF功能是从一个气象场多次观测资料中识别出主要空间型及其时间演变规律。EOF展开就是将气象变量场分解为空间函数（V）和时间函数（T）两部分的乘积之和：X=VT。应用步骤：资料预处理(距平或标准化处理)计
3、算协方差矩阵、用Jacobi方法或迭代法计算协方差矩阵的特征值与特征向量、将特征值从大到小排列、计算特征向量的时间系数、计算每个特征向量的方差贡献、结果输出二、实验目的：经验正交函数分解（EOF）是统计天气分析中气象要素场最基础的研究模型，是必须理解和掌握的方法之一，是后续课程中许多气象要素场的计算结果的理解的基础理论，也是毕业设计和论文中的基本分析方法。该方法用个数较少的几个空间分布模态来描述环流形势，而且基本涵盖环流场的信息，既能作为天气分析模型，其方法的延拓又能作为天气预报模型，在实际工作中也有极强的实用意义。通过该实验，深刻理解气象要
4、素场的统计模型的意义，掌握气象要素场分析的基本方法，为实际预报业务和科研工作打下一定的基础。三、涉及实验的相关情况介绍（包含使用软件或实验设备等情况）：本实验要求运用Fortran6.5编译软件以及grads计算机一台。四、实验结果（含程序、数据记录及分析和实验总结等，可附页）：要求编写主程序，其中包括资料读入，范围截取，子程序调用。注意：EOF的资料输入，时间场一维，空间场一维。*********************（附程序，对关键部分标志出）**********************EOF程序C*******************
5、***************************************************C*CPROGRAMNOTES*C*CTHISPROGRAMUSESEOFTOANALYSISTIMESERIES*COFMETEOROLOGICALFIELD*C*C**********************************************************************C*C********ParameterTable**********C*CMt===>LENTHOFTIMESERIES*CN===>N
6、UMBEROFGRID-POINTS(orSTATIONS)*CKS=-1,SELF;KS=0,DEPATURE;KS=1,STANDERDLIZEDDEPATURE*CKV=NUMBEROFEIGENVALUESWILLBEOUTPUT*CKVT=NUMBEROFEIGENVECTORSANDTIMESERIESWILLBEOUTPUT*CMNH=Minimum(Mt,N)*CEGVT===>EIGENVECTORS,ECOF===>TIMECOEFFICIENTSFOREGVT*CER(KV,1)====>LAMDA;LAMDA===>E
7、IGENVALUE*CER(KV,2)====>ACCUMULATELAMDA*CER(KV,3)====>THESUMOFCOMPONENTSVECTORSPROJECTEDONTO*CEIGENVACTOR.*CER(KV,4)====>ACCUMULATEER(KV,3)*C*C**********************************************************************PARAMETER(N=73*37,MT=58,MNH=58)PARAMETER(KS=1,KV=10,KVT=10)RE
8、ALF(N,MT),AVF(N),DF(N),ER(MNH,4)REALA(MNH,MNH),S(MNH,MNH),V(MNH)c*****************
