EOF分析及其應用最新版本

1、.1 組員：組員： 中國氣象科學研究院中國氣象科學研究院 .2一、引一、引 言言 經(jīng)驗正交函數(shù)（經(jīng)驗正交函數(shù)（EOFEOF）方法：最早由統(tǒng)計學家）方法：最早由統(tǒng)計學家pearsonpearson（19021902）提出，由）提出，由LorenzLorenz（19561956）引入氣象）引入氣象問題分析中。該方法以場的時間序列為分析對象，問題分析中。該方法以場的時間序列為分析對象，由于對計算條件要求甚高，直到由于對計算條件要求甚高，直到2020世紀世紀6060年代后期年代后期才在實際工作中得到廣泛應用。才在實際工作中得到廣泛應用。近近3030年來，出現(xiàn)了適合于各種分析目的的年來，出現(xiàn)了適合于各種

2、分析目的的EOFEOF分析方分析方法，如法，如擴展擴展EOFEOF（EEOFEEOF）方法）方法，旋轉(zhuǎn)旋轉(zhuǎn)EOFEOF（REOFREOF）方）方法，風場法，風場EOFEOF（EOFWEOFW）方法，復變量）方法，復變量EOFEOF（CEOFCEOF）方）方法。法。.3l1，沒有固定的函數(shù)；沒有固定的函數(shù)；l2 2，能在有限區(qū)域?qū)Σ灰?guī)則分布的站點進行分解；，能在有限區(qū)域?qū)Σ灰?guī)則分布的站點進行分解；l3 3，展開收斂快，很容易將變量場的信息集中在幾個模，展開收斂快，很容易將變量場的信息集中在幾個模態(tài)上；態(tài)上；l4 4，分離出的空間結(jié)構(gòu)具有一定的物理意義；，分離出的空間結(jié)構(gòu)具有一定的物理意義；l E

3、OFEOF方法不但用于觀測資料的分析，還用方法不但用于觀測資料的分析，還用于于GCMGCM資料的分析和數(shù)值模式的設計?，F(xiàn)在，資料的分析和數(shù)值模式的設計?，F(xiàn)在，EOFEOF方法已作為一種基本的分析手段頻繁地出方法已作為一種基本的分析手段頻繁地出現(xiàn)在大氣科學研究的文獻中。現(xiàn)在大氣科學研究的文獻中。.4二、二、EOF分析方法原理分析方法原理 將某氣候變量場的觀測資料以矩陣形將某氣候變量場的觀測資料以矩陣形式給出式給出m是空間點（觀測站或網(wǎng)格點），是空間點（觀測站或網(wǎng)格點），n是時間序列長度（觀測次數(shù)）。是時間序列長度（觀測次數(shù)）。 mnmmnnxxxxxxxxx212222111211X.5氣象場的

4、自然正交展開，是將氣象場的自然正交展開，是將X X分解為時間函分解為時間函數(shù)數(shù)Z Z和空間函數(shù)和空間函數(shù)V V兩部分兩部分, ,即即或者或者含義：含義：場中第場中第i i個格點上的第個格點上的第t t次觀測值，可以次觀測值，可以看作是看作是p p個空間函數(shù)個空間函數(shù) 和時間函數(shù)和時間函數(shù) 的線性組的線性組合合 。VZX mi,2, 1nt,2,1pkptiptitiktikitzvzvzvzvx12211ikvkizpk,2, 1.6其中，其中， 是第是第j個典型場，只是空間的函數(shù)。個典型場，只是空間的函數(shù)。 mmmmmmvvvvvvvvv212222111211V mnmmnnzzzzzzz

5、zz212222111211ZTmjjjjvvvv),(21.7第第t個空間場可表示為個空間場可表示為或者或者mtmtttzzzvvvx2211mtmmmmtmtmmtttzvvvzvvvzvvvxxx2122221211211121.8 上式表明，第上式表明，第t t個場可以表示為個場可以表示為m個空間個空間典型場，按照不同的權(quán)重線性疊加而成。典型場，按照不同的權(quán)重線性疊加而成。V的每一列表示一個空間典型場，由于這個場的每一列表示一個空間典型場，由于這個場由實際資料確定，故又叫經(jīng)驗正交函數(shù)。由實際資料確定，故又叫經(jīng)驗正交函數(shù)。 上述分解要求滿足下列兩個條件：上述分解要求滿足下列兩個條件：.9

6、jijizzintjtit01TjiZZmji,21jijivvpkkjkij101vvTi性性質(zhì)質(zhì).10三、分解方法三、分解方法A為實對稱矩陣，根據(jù)實對稱矩陣分解原理，為實對稱矩陣，根據(jù)實對稱矩陣分解原理，一定有一定有 或者或者TTVVZZXXTTXXA AVVTTVVA.11 V的列就是的列就是A的特征向量的特征向量， 是是A的的特征值特征值組成的對角矩陣。組成的對角矩陣。Z就是時間系數(shù)矩陣就是時間系數(shù)矩陣， 第第i個格點上的個格點上的第第t 個時間系數(shù)。個時間系數(shù)。XVZTmkktkiitxvz1itz.12四、誤差估計和計算四、誤差估計和計算 是擬合場是擬合場.可以證明誤差可以證明誤差

7、nppmZ.VXXX21111()pmnmititiiitiiQxx.13第第i個特征向量對個特征向量對X場的貢獻率場的貢獻率前前p個特征向量對個特征向量對X場的貢獻率場的貢獻率miiii1miiipiii11五、重要參數(shù)五、重要參數(shù).14六、計算步驟六、計算步驟1）根據(jù)分析目的，對原始資料矩陣）根據(jù)分析目的，對原始資料矩陣X作距作距平或者標準化處理；平或者標準化處理；2）由）由X求協(xié)方差矩陣求協(xié)方差矩陣 ；3）求實對稱矩陣）求實對稱矩陣A的全部特征值的全部特征值 、特征、特征向量向量 ，h=1H（通常使用（通常使用Jacobi法）；法）；hVhTAXX.15 4）將特征值作非升序排列（通常使

8、用）將特征值作非升序排列（通常使用沉浮法），并對特征向量序數(shù)作相應變動；沉浮法），并對特征向量序數(shù)作相應變動； 5）根據(jù)）根據(jù) , ,h=1H和和X總方差，求出全總方差，求出全部部 、 , , h=1H； 6 6）由）由X及主要及主要 求其時間系數(shù)求其時間系數(shù) 、h h=1=1H，主要的數(shù)量由分析目的及分析對象，主要的數(shù)量由分析目的及分析對象定；定； 7 7）輸出主要計算結(jié)果。）輸出主要計算結(jié)果。 hhVhZhhP.16七七 經(jīng)驗正交函數(shù)的物理意義經(jīng)驗正交函數(shù)的物理意義 特征向量以及時間系數(shù)的分析。特征向量以及時間系數(shù)的分析。 vvvv第一特征向量（第一空間典型場）是與第一特征向量（第一空間典

9、型場）是與n n張張X X圖平均最相似的，或者說具有與所要展開的圖平均最相似的，或者說具有與所要展開的資料矩陣的資料矩陣的n n個樣本最相似的特征。比如：若原個樣本最相似的特征。比如：若原始資料矩陣是始資料矩陣是7 7月份月份5050年實測將水場（非距平年實測將水場（非距平場），則第一特征向量就可以解釋為這場），則第一特征向量就可以解釋為這5050年的年的平均場，其相應的時間系數(shù)基本對應我國大尺平均場，其相應的時間系數(shù)基本對應我國大尺度旱澇年。但當降水場由距平組成，第一特征度旱澇年。但當降水場由距平組成，第一特征向量就解釋為與向量就解釋為與5050年夏季距平場最相似的特征年夏季距平場最相似的特

10、征場，它指出了我國夏季經(jīng)常出現(xiàn)的大尺度澇區(qū)場，它指出了我國夏季經(jīng)常出現(xiàn)的大尺度澇區(qū)和旱區(qū)。和旱區(qū)。.17八、時空轉(zhuǎn)換問題八、時空轉(zhuǎn)換問題 當當 時，先求出時，先求出 的特征值，的特征值，然后求然后求 的特征向量，這種方法叫時的特征向量，這種方法叫時空轉(zhuǎn)換。空轉(zhuǎn)換。 令令 的特征值為的特征值為 ，其特征向量，其特征向量為為 ， 的特征值也為的特征值也為 ，其，其特征向量為特征向量為 nm XXTTXXXXTiiuTXXiiv.18iiiXuvmvvvV,21XVZT.19l分析表明，南亞夏季風的爆發(fā)主要體現(xiàn)在降水的突然增加和季風雨帶的快速推進上，雨帶的時空分布有突變的特點。l第1 模態(tài)降水量的突

11、然增加。l第2 模態(tài)從南向北的快速推進過程。l第3模態(tài)東西分布型態(tài)，及在季風爆發(fā)后印度半島降水快速增加的過程。l第4模態(tài)印度次大陸東海岸降水的準雙周振蕩型態(tài)。.20.21.22.23.24l將東亞500 hPa 風場,溫度場,降水量場進行經(jīng)驗正交分解,得到它們的主要模態(tài)的時空變化特征。結(jié)果表明東亞風場EOF 的主要模態(tài)與我國溫度,降水量的EOF 的主要模態(tài)對應,其第一EOF 模態(tài)與盛夏溫度,降水量的關(guān)系密切。.25.26圖中,淡灰(灰,黑色) 是0. 05 (0. 01 ,0. 001) 信度的相對區(qū).27圖a 圖b .28l魏鳳英，現(xiàn)代氣候統(tǒng)計診斷與預測技術(shù)，氣象出版社，北京，2007;lE

12、OF在大氣科學研究中的新進展；丁裕國；氣象科技199303期 ；l近年來中國統(tǒng)計氣象學的新進展；周家斌 黃嘉佑 ；氣象學報 1997年03期 ；l我國盛夏500 hPa 風場的EOF 分析及其與大尺度氣候異常的關(guān)系。顧澤，封國林，顧駿強，施能；l南亞夏季風爆發(fā)前后降水量時空變化特征。朱 敏，張 銘；.29Thank you!.30lKarl Pearson（18571936） 統(tǒng)統(tǒng)計學計學之父之父。 lK. Pearson 1879年畢業(yè)于劍橋大學數(shù)學系；1884年進入倫敦大倫敦大學學院學學院 。 lK. Pearson 最重要的學術(shù)成就，是為現(xiàn)代統(tǒng)計學打下基礎。K. Pearson 在1893-1912年間寫出18篇在演化論上的數(shù)學貢獻的文章，而這門算術(shù)，也就是今日的統(tǒng)計。許多熟悉的統(tǒng)計名詞如標準差，成分分析，卡方檢定都是他提出的。.31l混沌理論之父、MIT教授lOne of the great modern science stories is the so-