大气统计预报期末考已结束lecture.pptxVIP

第七章气象变量场时空结构分离;7.1 主分量分析 7.2 经验正交分解 7.3 扩展经验正交分解 7.4 旋转经验正交分解 7.5 主振荡型分析;某一区域的气候变量场通常由许多个观测站点或网格点构成，这给直接研究其时空变化特征带来困难。如果能用个数较少的几个空间分布模态来描述原变量场，且又能基本涵盖原变量场的信息，这将十分有助于我们分析气候要素场的时间和空间特征。这就要寻找某种数学表达式将变量场的主要空间分布结构有效的分离出来。气候统计诊断中，应用最广的办法就是把原变量场分解为正交函数的线性组合，构成为数较少的不相关典型模态，代替原始的气候变量场;7.1 主分量分析;主成分的性质;7.2 经验正交分解 EOF;EOF 分解技术介绍;利用经验正交函数展开，就是把X 分解成正交的空间函数F与正交的时间函数T的乘积。即： 根据正交性，F和T应满足下面条???： 当k=l时 当k≠l时 ;那我们现在的任务就是求出F与T 的形式，下面看如何求它们。先观察一下X场的协方差矩阵S。有 , 根据线性代数的知识，S是一个m 乘m 的实对称矩阵。 将 带入，并利用 ,及 得: 右乘F得， 即： ;或: 这实际上是关于 的m方程组，若只考虑第k 个方程组，则有： ;或: ……….. 这是一个线性齐次方程组，根据线性代数的知识，它有非零解的充要条件是：系数行列式为零。此时我们把当作未知数，则它们有非零解的条件是： 这其实是矩阵S的特征方程， 是其特征值，而 是对应的特征向量。 ;由此可得到m 个特征值 ，及其对应的m个特征向量，将这些特征向量放到一起就是我们所要求得的F。它所代表的意义就是：将X 分成了m 个互相独立、互不相关的空间模态. 这m 个模态中，哪个是最主要的模态呢？哪个是次要模态呢？ 可将所求得的m 个特征值 按由大到小的顺序排列，即： 则由 所对应的特征向量是最主要的模态， 对应的则是第二位重要， 对应的则是最不重要的模态。 ;那末， 所反映的模态占X 的分量有多少呢？由下式r决定： （由于S是实对称的，所以， 均为正值。） 得到了F，对应的T 就很容易得到了： ;每个特征向量所对应的时间向量（称为时间系数）就反映了X 的空间区域中由此特征向量所表示的空间型的时间变化特征。系数绝对值越大表明对应时刻（某月、某年）这种空间分布型（模态）越明显（典型）。 由于EOF的收敛速度很快，一般特征值较大的前几个模态就能反映出气候场X的主要特征了。;(1) 它没有固定的函数，不像有些分解需要有一种特殊的函数作基函数，如球谐函数。 (2) 它能在有限的区域内对不规则分布的站点进行分解。 (3) 它的展开收敛速度快，很容易将变量场的信息集中在几个模态上。 (4) 分离出的空间模态具有一定的物理意义; 计算步骤 EOF的一般计算步骤是： (1) 对原始资料矩阵X作距平或标准化处理。然后计算其协方差矩阵 用求实对称矩阵的特征值和特征向量的方法（最常用的是Jacobi 方法）求出S 的特征值 和特征向量F (3) 利用 的关系求出T。 (4) 计算每个特征向量所占的方差贡献 ,前p 个特征向量所占的累积方差贡献 . (5) 凭借气候学知识对前几项有意义的模态（mode）作分析。;摘自Dai等（2004）;The first mode of EOF of SLP is called AO. The corresponding time series is AO index;资料：SST； 格点资料，空间格点数5670； 时间长度：1953-2003年月资料，612个月； 采用标准化数据进行PCA。;;;;;我国盛夏500 hPa 风场的EOF 分析及其与大尺度气候异常的关系;1951 - 2004 年500 hPa 7 - 8 月平均风场EOF 分解的第一特征向量图(a)和标准化的时间权重系数(b) (图中斜直线是回归线);　盛夏500 hPa 第一时间权重系数与我国同期平均温度(a) 和月总降水量(b) 的相关系数图;　1951 - 2004 年我国7 - 8 月月平均温度(a) 和月总降水量(b) 的第一特征向量图;7.3 扩展经验正交分解(EEOF) ;EEOF和EOF的方法类似，关键是计算协方差矩阵的资料矩阵是由几个连续的时间上的观测值构成的，相当于