氣象資料的分析與預測建模摘要：本文建立了用于氣象資料的分析與預測的數(shù)學模型。經(jīng)比照該城市與北京的海拔、氣候等極為相似，因此，我們以北京的標準氣象指數(shù)為參照建立模型。首先針對問題一：對該城市兩年來的總體氣象進行整體評價，并對該城市氣候走勢進行中長期預測。我們仿照科學家對環(huán)境空氣質(zhì)量綜合指數(shù)評價的數(shù)學模型，以第一年每個月的平均氣壓、平均氣溫、平均相對濕度、平均風速、最高氣壓、最高氣溫、最高相對濕度、最高風速和北京的標準氣候指數(shù)為參數(shù)，通過matlab建立與之相關的方程來確定該城市當月的氣象質(zhì)量指數(shù)，按照指數(shù)數(shù)值的大小分為優(yōu)、良、差三大類，從而評價每個月的氣候質(zhì)量。運用第二年的數(shù)據(jù)進行檢驗模型的正確性：隨機選取幾個月的氣象因素數(shù)據(jù)，并各自與對應的北京標準氣象數(shù)據(jù)做差，數(shù)值越小那么氣象質(zhì)量越好，將分析結(jié)果與通過權(quán)重綜合指數(shù)法計算得出的結(jié)論做比擬。跟據(jù)建立的氣象質(zhì)量評價數(shù)學模型和第一、第二兩年數(shù)據(jù)比照趨勢圖，對該城市氣候進行整體評價和中長期的分析預測。然后針對問題二：對影響極端天氣發(fā)生的主要指標，比方：降水、溫度等建立監(jiān)控預報體系的數(shù)學模型，并用兩年內(nèi)的累積氣象資料進行驗證。我們運用多元線性回歸分析的數(shù)學方法，建立了監(jiān)控預報最高溫度的數(shù)學模型。該模型中我們先假設了最高溫度的主要影響因素是平均氣壓、平均氣溫、平均濕度、日照時數(shù)、地面平均溫度、降水量等，通過matlab編寫程序驗證取舍得出平均氣壓、平均氣溫、平均濕度、日照時數(shù)、地面平均溫度是影響降水和溫度的主要影響因素；然后，檢驗多元線性回歸方程的擬合優(yōu)度、相關性；最后，帶入兩年內(nèi)的累積氣象資料進行驗證。最后我們評價了模型的優(yōu)缺點，并對模型的缺乏之處進行了改良。關鍵詞：權(quán)重綜合氣象質(zhì)量指數(shù)；多元線性回歸；正態(tài)分布。1.問題重述近年來，我國極端天氣呈現(xiàn)出發(fā)生頻率加大、致災性加重等新特點，極端天氣趨于常態(tài)化。雖然局部地方加大防災減災建設并取得一些成效，但相比現(xiàn)實需求，對極端天氣監(jiān)測預警手段仍然缺乏，防御應對體系建設仍存在明顯短板。附件中是某城市兩年內(nèi)連續(xù)的日氣象資料，包括氣壓、溫濕度、降水量、風力風向等多項氣象資料指標。請你完成以下任務：(1)對該城市兩年來的總體氣象進行整體評價，并對該城市氣候走勢進行中長期預測；請詳細給出評價的指標體系以及評級和預測的數(shù)學模型；(2)對影響極端天氣發(fā)生的主要指標，比方：降水、溫度等建立監(jiān)控預報體系的數(shù)學模型，并用兩年內(nèi)的累積氣象資料進行驗證。注意：這里的主要指標并不限于降水和溫度等指標，你們也可根據(jù)實際需求自行選擇。注：該城市的海拔約為30-50米。2.問題的背景與分析雖然我國幅員遼闊，地形復雜，但各地的氣象在空間分布上仍有一定規(guī)律。我國分布著世界上最大的溫帶季風區(qū)，秦嶺淮河以北是溫帶季風氣候，這里夏季高溫多雨，冬季寒冷枯燥。冬冷夏熱，雨熱同期；秦嶺淮河以南是亞熱帶季風氣候，這里夏季高溫多雨，冬季溫和少雨，熱量充足，氣溫年較差較小，降水豐富，但季節(jié)變化較大；西部的新疆、寧夏、內(nèi)蒙古、青海、甘肅等多是溫帶大陸性氣候，這里冬冷夏熱，年溫差大，降水集中，四季清楚，年雨量較少，大陸性強；面積廣闊的青藏高原等地是高原山地氣候，這里海拔高，氣溫低，但輻射強，日照豐富，降水少，冬半年風力強勁，氣溫的年較差小，日較差大。近年來，我國極端天氣呈現(xiàn)出發(fā)生頻率加大、致災性加重等新特點，極端天氣趨于常態(tài)化。雖然局部地方加大防災減災建設并取得一些成效，但相比現(xiàn)實需求，對極端天氣監(jiān)測預警手段仍然缺乏，防御應對體系建設仍存在明顯短板。向社會提供準確及時的天氣監(jiān)控預測是我們的宗旨；滿足人民對氣象信息的多種需求是我們的目標。因此，準確的對極端天氣監(jiān)測預警，有著十分重要的意義。我們建立的這個數(shù)學模型就是本著對極端天氣監(jiān)測預警的角度出發(fā)，以城市多年的歷史數(shù)據(jù)為參照，建立氣象評價體系數(shù)學模型和對極端溫度、風速的監(jiān)控預報體系數(shù)學模型。模型的假設與符號說明模型的假設：1．假設氣象部門提供的實測數(shù)據(jù)是準確的，能較真實地反映該城市的氣象情況。2．假設北京氣象與該城市氣象相似。3．假定網(wǎng)上所給的北京標準氣候指標可靠。4.假定最高氣溫與與平均氣壓、平均氣溫、平均濕度、日照時數(shù)、地面平均溫度成線性函數(shù)。5.假定最高氣溫時隨機變量，服從均值為零的正態(tài)分布。符號說明：Ii——第i項氣象因素指數(shù)；N——參數(shù)項數(shù)；Cimax——第i項氣象因素(月)均最大值；Ci——第i項氣象因素(月)平均值；Si——第i頂氣象因素標準值；Qi——第i項氣象因素指數(shù)的權(quán)數(shù)；Pi——第i項氣象因素指數(shù)的修正；I——綜合氣象質(zhì)量指數(shù)；K——權(quán)重綜合氣象質(zhì)量指數(shù)。Y——最高氣溫；x1——平均氣壓；x2——平均氣溫；x3——平均濕度；x4——日照時數(shù)；x5——地面平均氣溫；bi——Xi的回歸系數(shù)；3.氣象評價體系模型的建立與修正為了能夠更客觀地評價和預測某地的氣象，我們仿照科學家對環(huán)境空氣質(zhì)量綜合指數(shù)評價的數(shù)學模型，以每個月的平均氣壓C1、平均氣溫C2、平均相對濕度C3、平均風速C4、最高氣壓C1極、最高氣溫C2極、最高相對濕度C3極、最高風速C4極和北京的標準氣候指數(shù)Si為參數(shù)，通過matlab建立與之相關的函數(shù)方程來確定該城市當月的綜合氣象質(zhì)量指數(shù)I，然后在評級列表中查出氣象質(zhì)量等級。3.1綜合氣象質(zhì)量指數(shù)法計算公式:〔1〕〔2〕其中第i頂氣象因素評價標準Si從表1—1北京氣象標準指數(shù)表中獲取。表3—1北京市標準氣象值平均氣壓平均氣溫極端最高氣溫極端最低氣溫降水量平均相對濕度平均風速(%)1月1024.2-3.712.9-18.32.7442.62月1022-0.719.8-164.9442.83月1017.45.826.4-158.3463.14月101014.233-3.221.2463.25月1005.719.936.82.634.2532.96月1001.224.439.29.878.1612.57月999.726.241.916.6185.2752.18月1003.724.936.111.4159.7771.89月1010.52034.44.345.568210月1016.713.129.3-3.521.8612.111月1021.34.622-10.67.4572.412月1023.8-1.519.5-15.62.8492.6運用公式〔1〕〔2〕得出該市綜合氣象質(zhì)量指數(shù)I，然后從表1—2中對應查出氣象質(zhì)量等級。表3—2氣象質(zhì)量分級標準綜合指數(shù)0.90~1.101.10—1.301.10—1.30評價等級優(yōu)良差可見,綜合氣象質(zhì)量指數(shù)法計算簡單,但綜合氣象質(zhì)量指數(shù)采用算術(shù)平均值計算,所參加評價的氣象環(huán)境因素指數(shù)在計算中權(quán)重值相同,主要污染物對環(huán)境空氣質(zhì)量的影響評價結(jié)果偏低。3.2綜合指數(shù)法的修正為能更加客觀地評價氣象質(zhì)量,在綜合指數(shù)法的根底上,給各單項指數(shù)根據(jù)一定條件賦予一個權(quán)重值,各單項指數(shù)與權(quán)重值的乘積之和為評價綜合氣象質(zhì)量指數(shù),并將此法稱為權(quán)重綜合氣象質(zhì)量指數(shù)法,計算公式如下:設(i=1,2,3…n)〔3〕假設S為的個數(shù),且,t為的個數(shù),,分指數(shù)權(quán)重賦值公式:0.05,〔4〕那么且?！?〕修正綜合指數(shù)法說明:對分指數(shù)太小或太大分指數(shù)賦值,即或時,對分指數(shù)取值偏差不大,時,用加權(quán)算術(shù)平均數(shù)計算分指數(shù)權(quán)數(shù),用公式(4)求出分指數(shù)權(quán)數(shù),再由公式(5)求出綜合指數(shù)。3.3權(quán)重綜合指數(shù)法的應用表3—3該市第一年氣象質(zhì)量指數(shù)對照表月份氣壓指數(shù)氣溫指數(shù)濕度指數(shù)風力指數(shù)綜合指數(shù)法質(zhì)量等級權(quán)重綜合指數(shù)質(zhì)量等級一0.99171.06511.41950.84271.0798優(yōu)1.1214良二0.99331.03391.36521.27781.1676良1.1888良三0.99031.02881.35780.70671.0209優(yōu)1.0731優(yōu)四1.00021.02311.53800.67081.0580優(yōu)1.1490良五0.99751.01871.21310.70670.9840優(yōu)1.0173優(yōu)六1.00181.01211.28500.57800.9692優(yōu)1.0351優(yōu)七0.99881.00771.11200.60330.9305優(yōu)0.9710優(yōu)八0.99921.00791.09180.85450.9883優(yōu)0.9957優(yōu)九0.99781.01481.14400.60990.9416優(yōu)0.9839優(yōu)十1.00171.01311.11920.62380.9395優(yōu)0.9771優(yōu)十一1.00141.02341.25660.63530.9792優(yōu)1.0226優(yōu)十二0.99901.02341.66100.58011.0659優(yōu)1.1332良從表3—3可以看出,使用兩種方法計算的綜合污染指數(shù)略有不同,權(quán)重綜合氣象質(zhì)量指數(shù)法計算結(jié)果略高于綜合氣象質(zhì)量指數(shù)法,該市第一年的三月、五月、六月、七月、八月、九月、十月、十一月的權(quán)重氣象質(zhì)量指數(shù)均為優(yōu)。3.4權(quán)重綜合指數(shù)法的檢驗我們常用單個氣候因素的真實值與標準值作差或商，來判斷氣候質(zhì)量。所以可以用氣候的差或商值來檢驗模型的正確性。將該市第二年每個月的平均氣壓、平均氣溫、平均相對濕度、平均風速、與北京市標準氣象值做比照，畫出折線圖如圖：表3—5該市第二年月平均氣壓與標準氣壓比照折線圖表3—6該市第二年月平均氣溫與標準氣溫比照折線圖表3—7該市第二年月平均濕度與標準濕度比照折線圖風速表隨機選取第二年中數(shù)月〔如二月份、七月份、十一月〕的氣象因素數(shù)值進行做差檢驗，差值較小者的為優(yōu)，其次為良，最后為差。做差得：二月份平均氣壓差值：1018.2-1022=-3.8平均氣溫差值：-1.4+0.7=-0.7平均濕度差值：0.58-0.44=0.14平均風力差值：1.4-2.8=-1.4七月份平均氣壓差值：993.6-999.7=-6.1平均氣溫差值：28.7-26.2=2.5平均濕度差值：0.62-0.75=-0.13平均風力差值：1.5-2.1=-0.6十一月份平均氣壓差值：1011.5-1021.3=-9.8平均氣溫差值：9.9-4.6=5.3平均濕度差值：0.45-0.57=-0.12平均風力差值：1.3-2.4=-1.1分別運用綜合氣象質(zhì)量指數(shù)法和權(quán)重綜合氣象質(zhì)量指數(shù)法計算得二月、七月、十一月質(zhì)量等級如表3—4表3—4該市第二年氣象質(zhì)量指數(shù)對照表月份綜合指數(shù)質(zhì)量等級權(quán)重綜合指質(zhì)量等級二月1.0791優(yōu)1.1901良七月0.9973優(yōu)1.0161優(yōu)十一月1.39335差1.803差經(jīng)做差分析知七月的各個氣候因素差值相對較小，氣象質(zhì)量等級應最高，十一月的各個氣候因素差值相對較大，氣象質(zhì)量等級應最低，此結(jié)果正與運用權(quán)重綜合氣象質(zhì)量指數(shù)法計算結(jié)果相同。3.5結(jié)論：權(quán)重綜合污染指數(shù)法對極端值進行了處理,根據(jù)分指數(shù)的大小不同分別賦予其不同的權(quán)重值；與綜合氣象質(zhì)量指法相比，對主要氣象因素在權(quán)重上有所側(cè)重,運用此法對該市第二年的數(shù)據(jù)進行檢驗，評價結(jié)果也與實際情況能較好吻合。4.對該市氣候的整體評價及中長期預測4.1對該市氣候的整體評價：該市的氣候整體表現(xiàn)為夏季高溫多雨，冬季寒冷枯燥，冬冷夏熱，雨熱同期，屬于典型的溫帶季風氣候。4.2對該市氣候的中長期預測：表4—1預測對照數(shù)據(jù)日期綜合氣象評估指數(shù)檢測數(shù)據(jù)權(quán)重氣象評估檢測數(shù)據(jù)1月份1.07981.12142月份1.16761.07911.18881.19013月份1.02091.07314月份1.0581.1495月份0.9841.01736月份0.96921.03517月份0.93050.99730.9711.01618月份0.98830.99599月份0.94160.983910月份0.93950.977111月份0.97921.39341.02261.80312月份1.06591.1332表4—2預測對照圖根據(jù)建立的氣象質(zhì)量評價數(shù)學模型和第一、第二兩年數(shù)據(jù)比照趨勢圖，分析預測得綜合氣象質(zhì)量會逐漸變差,平均氣壓有緩慢降低趨勢，平均氣溫有緩慢升高趨勢，最高氣溫緩慢降低趨勢，最低氣溫緩慢升高趨勢，平均濕度緩慢降低趨勢，降水量緩慢降低，平均風力緩慢增大等。5．極端天氣監(jiān)測預報體系模型的建立5.1建立模型：為了能夠更客觀地對極端天氣監(jiān)測預警，我們對影響極端天氣發(fā)生的主要指標之一——最高氣溫Y，建立監(jiān)測預報的數(shù)學模型。由假設可知，最高氣溫Y是隨機變量，它服從均值為零的正態(tài)分布，平均氣壓x1、平均氣溫x2、平均濕度x3、日照時數(shù)x4、地面平均溫度x5等是影響最高氣溫Y的因素，所以可以建立多元線性回歸模型：建模所需的知識點及其在Matlab中實現(xiàn)的方法：求解多元線性回歸方程用Matlab中的regress函數(shù)，函數(shù)命令形式為：[b,bint,r,rint,stats]=regress(Y,x,alpha)其中因變量數(shù)據(jù)向量Y和自變量數(shù)據(jù)向量均按以下方式輸入：Y=[y1y2…yn]x=[x1x2…xn]Alpha為顯著性水平〔缺省時設定為0.05〕Regress函數(shù)的功能：b為回歸系數(shù)的估計值b1，b2，b3，b4，b5，bint為回歸系數(shù)估計值b的置信區(qū)間；r、rint為殘差向量及其置信區(qū)間；stats是用于檢驗回歸模型的統(tǒng)計量，第一個是R2,稱為決定系數(shù)，R2越接近1說明引入方程的自變量與因變量的相關程度越高，xi與Y的回歸效果越好。R是相關系數(shù)，第二個是F統(tǒng)計量，第三個是與統(tǒng)計量相對應的概率p，當p<α時，說明回歸模型假設成立〔即Y與x的關系〕。模型求解在Matlab中的實現(xiàn)：在Matlab中求解該模型的程序代碼見附錄。輸入原始數(shù)據(jù)后的結(jié)果如下：b=-44.07830.04941.1526-0.04790.0052-0.0545bint=-103.522815.3661-0.00890.10761.00691.2982-0.0741-0.0218-0.00010.0106-0.17380.0649stats=1.0e+003*0.00102.95260從上述結(jié)果可知：回歸方程如果出現(xiàn)某一項因素的系數(shù)過小，那么可以刪除這個因子，重新計算回歸方程。多元線性回歸方程中的相關性檢驗：由運行結(jié)果可知：F=2952.6，查F分布的上α分位數(shù)Fα〔m，n〕表可得：F0.05(5,23)=2.53F>F0.05(5,23)，可以認為Y與x1、x2、x3、x4、x5之間顯著地有線性關系。多元線性回歸方程中的擬合優(yōu)度檢驗：由運行結(jié)果可知：R2=1，多元線性回歸方程擬合度很高。多元線性回歸方程中的方差檢驗：δ2==0.11864其中Q為5.2驗證模型：多元線性回歸方程大致反映了變量Y與xi之間的變化規(guī)律。但是，由于Y與x之間的關系不是確定的，所以對于一組x的值，只能得到Y(jié)i的估計值。它的精確性如何？對此，我們應當對Yi進行區(qū)間估計，即對于給定的置信度1-α，求出Yi的置信區(qū)間。對于x的任一值，我們得到置信度為95%的預測區(qū)間是〔-1.96s,+1.96s〕其中當n充分大時，xi與Yi近似的服從正態(tài)分布N〔i，s2〕方差為S2=Se=標準差S=根據(jù)兩年內(nèi)的累積氣象資料可得到相應的置信區(qū)間：表5—1置信區(qū)間的分布第一年氣象資料月份平均氣壓/hpa平均氣溫/oC平均濕度/%日照時數(shù)/hour地面平均溫度/oC最高氣溫/oC/oC置信區(qū)間11009.19.647220.910.315.415.17(14.50，15.85)210115.541291.22.811.611.60(10.93，12.28)31009.19.647220.910.315.415.17(14.50，15.85)41002.717.155240.919.622.922.71(22.04，23.39)5998.721.653277.925.327.527.68(27.01，28.36)699824.964213.427.930.630.45(29.77，31.12)7994.426.474185.427.831.531.38(30.70，32.05)8998.824.776152.126.328.929.45(28.77，30.12)9100521.768204.423.12727.13(26.45，27.80)101012.714.959185.415.221.320.43(19.76，21.11)111014.4857197.56.913.413.17(12.50，13.85)121018.40.26997.9-143.72(3.04，4.39)第二年氣象資料月份平均氣壓/hpa平均氣溫/oC平均濕度/%日照時數(shù)/hour地面平均溫度/oC最高氣溫/oC/oC置信區(qū)間11017.1-1.954167.5-7.42.42.66〔1.99，3.34〕21018.2-1.45898.5-0.82.12.38〔1.71，3.06〕31012.383426610.113.914.35〔13.68，15.03〕41002.317.836250.120.824.324.39〔23.72，25.07〕5998.820.753247.825.126.226.50〔25.83，27.18〕6991.327.850203.531.834.433.86〔33.19，34.54〕7993.628.762144.932.134.134.12〔33.44，34.79〕8997.52670170.427.830.931.18〔30.51，31.86〕91006.1226416823.827.127.49〔26.82，28.17〕101012.215.851189.915.621.621.83〔21.15，22.50〕111011.59.945195.4815.515.73〔15.05，16.40〕由以上兩表可知，95%以上的Yi值都在置信區(qū)間內(nèi)，即可驗證極端天氣監(jiān)測預報模型的正確合理性。5.3模型的評價：極端天氣監(jiān)測預報模型最大優(yōu)點在于對原始數(shù)據(jù)擬合時,采用多種方法進行,使之愈來愈完善,具有很高的擬合精度和適度性在此根底上,對模型作進一步討論便可得到一系列可靠而實用的信息。并且所得結(jié)論與客觀事實很好地吻合,從而進一步說明模型是合理的。不過，此模型的影響因素有限，而且僅用每年的月平均值計算，不能全面預測到每一天。6.參考文獻[1]王麗霞概率論與數(shù)理統(tǒng)計——理論、歷史及應用大連理工大學出版社，2010.10[2]費浦生羿旭明數(shù)學建模及其根底知識詳解武漢大學出版社，2006.57.附錄模型一Matlab命令：函數(shù)式functionIi=f4(Ci,Di,Si)

Ii=(sqrt(Ci.*Di))./Sifunction=f6(Ai,Bi,Ci,Ei)Mi=Ai+Bi+Ci+Ei

氣壓均值Di=[1022.4,1019.3,1006,1017.7,1007.7,1008,1002.7,1007.1,1011.6,1024.1,1031.1,1027.2]Si=[1024.2,1022,1017.4,1010,1005.7,1001.2,999.7,1003.7,1010.5,1016.7,1021.3,1023.8]Ci=[1009.1,1011,1009.1,1002.7,998.7,998,994.4,998.8,1005,1012.7,1014.4,1018.4 ];Ii=f4(Ci,Di,Si)Ii=0.99170.99330.99031.00020.99751.00180.99880.99920.99781.00171.00140.9990氣溫均值Ci=[9.6,5.5,9.6,17.1,21.6,24.9,26.4,24.7,21.7,14.9,8,0.2];Si=[-3.7,-0.7,5.8,14.2,19.9,24.4,26.2,24.9,20,13.1,4.6,-1.5];Di=[18.1,11.6,18.1,24.6,29.2,31.1,30.6,29.8,27,18.8,14.2,9.6];Ii=f4(Ci,Di,Si)Ii=Ii=Columns1through11-3.5627-11.41072.27271.44441.26201.14051.08481.08961.21031.27762.3170Column12-0.9238華氏溫度Ci=282.6000278.5000282.6000290.1000294.6000297.9000299.4000297.7000294.7000287.9000281.0000273.2000>>Si=269.3000272.3000278.8000287.2000292.9000297.4000299.2000297.9000293.0000286.1000277.6000271.5000Di=291.1000284.6000291.1000297.6000302.2000304.1000303.6000302.8000300.0000291.8000287.2000282.6000Ci=[282.6000278.5000282.6000290.1000294.6000297.9000299.4000297.7000294.7000287.9000281.0000273.2000];Di=[291.1000284.6000291.1000297.6000302.2000304.1000303.6000302.8000300.0000291.8000287.2000282.6000];Si=[269.3000272.3000278.8000287.2000292.9000297.4000299.2000297.9000293.0000286.1000277.6000271.5000];Ii=f4(Ci,Di,Si)Ii=1.06511.03391.02881.02311.01871.01211.00771.00791.01481.01311.02341.0234溫濕Ci=[0.47,0.41,0.47,0.55,0.53,0.64,0.74,0.76,0.68,0.59,0.57,0.69];Di=[0.83,0.88,0.83,0.91,0.78,0.96,0.94,0.93,0.89,0.79,0.90,0.96];Si=[0.44,0.44,0.46,0.46,0.53,0.61,0.75,0.77,0.68,0.61,0.57,0.49];Ii=f4(Ci,Di,Si)Ii=Ii=Columns1through111.41951.36521.35781.53801.21311.28501.11201.09181.14401.11921.2566Column121.6610風力Ci=[1.5,1.6,1.5,1.44,1.4,1.16,1.07,0.91,0.93,0.78,0.93,0.91];Si=[2.6,2.8,3.1,3.2,2.9,2.5,2.1,1.8,2.0,2.1,2.4,2.6];Di=[3.2,8,3.2,3.2,3,1.8,1.5,2.6,1.6,2.2,2.5,2.5];Ii=f4(Ci,Di,Si)Ii=Ii=Columns1through110.84271.27780.70670.67080.70670.57800.60330.85450.60990.62380.6353Column120.5801求和Ai=[0.99170.99330.99031.00020.99751.00180.99880.99920.99781.00171.00140.9990];Bi=[1.06511.03391.02881.02311.01871.01211.00771.00791.01481.01311.02341.0234];Ci=[1.41951.36521.35781.53801.21311.28501.11201.09181.14401.11921.25661.6610];Ei=[0.84271.27780.70670.67080.70670.57800.60330.85450.60990.62380.63530.5801];Mi=f6(Ai,Bi,Ci,Ei)Mi=4.31904.67024.08364.23213.93603.87693.72183.95343.76653.75783.91674.2635矩陣Iij=0.99170.99330.99031.00020.99751.00180.99880.99920.99781.00171.00140.99901.06511.03391.02881.02311.01871.01211.00771.00791.01481.01311.02341.02341.41951.36521.35781.53801.21311.28501.11201.09181.14401.11921.25661.66100.84271.27780.70670.67080.70670.57800.60330.85450.60990.62380.63530.5801Iij的反轉(zhuǎn)矩陣0.99171.06511.41950.84270.99331.03391.36521.27780.99031.02881.35780.70671.00021.02311.53800.67080.99751.01871.21310.70671.00181.01211.28500.57800.99881.00771.11200.60330.99921.00791.09180.85450.99781.01481.14400.60991.00171.01311.11920.62381.00141.02341.25660.63530.99901.02341.66100.5801Gi=[0.99171.06511.41950.8427];Q1=Gi./4.3190Q1=0.22960.24660.32870.1951>>Gi=[0.99331.03391.36521.2778];Q2=Gi./4.6702Q2=0.21270.22140.29230.2736>>Gi=[0.99031.02881.35780.7067];Q3=Gi./4.0836Q3=0.24250.25190.33250.1731>>Gi=[1.00021.02311.53800.6708];Q4=Gi./4.2321Q4=0.23630.24170.36340.1585>>Gi=[0.99751.01871.21310.7067];Q5=Gi./3.936Q5=0.25340.25880.30820.1795>>Gi=[1.00181.01211.28500.5780];Q6=Gi../3.8769Q6=0.25840.26110.33150.1491>>Gi=[0.99881.00771.11200.6033];Q7=Gi./3.7218Q7=0.26840.27080.29880.1621>>Gi=[0.99921.00791.09180.8545];Q8=Gi./3.9534Q8=0.25270.25490.27620.2161>>Gi=[0.99781.01481.14400.6099];Q9=Gi./3,7665Q9=0.33260.33830.38130.2033>>Gi=[1.00171.01311.11920.6238];Q10=Gi./3.7578Q10=0.26660.26960.29780.1660>>Q11=[1.00141.02341.25660.6353];Q11=Gi./3.9167Q11=0.25580.25870.28580.1593>>Q12=[0.99901.02341.66100.5801];Q12=Gi./4.2635Q12=0.23490.23760.26250.1463Qij=0.22960.24660.32870.19510.21270.22140.29230.27360.24250.25190.33250.17310.23630.24170.36340.15850.25340.25880.30820.17950.25840.26110.33150.14910.26840.27080.29880.16210.25270.25490.27620.21610.33260.33830.38130.20330.26660.26960.29780.1660.25580.25870.28580.15930.23490.23760.26250.1463Qij的反轉(zhuǎn)矩陣Aij=0.22960.21270.24250.23630.25340.25840.26840.25270.33260.26660.25580.23490.24660.22140.25190.24170.25880.26110.27080.25490.33830.26960.25870.23760.32870.29230.33250.36340.30820.33150.29880.27620.38130.29780.28580.26250.19510.27360.17310.15850.17950.14910.16210.21610.20330.16600.15930.1463Qij矩陣列項求和Ri=1.00001.00001.00000.99990.99991.00011.00010.99991.25551.00000.95960.8813Qi=[0.22960.24660.32870.1951];Ii=[0.99171.06511.41950.8427];P1=(Qi./1).*IiQi=[0.21270.22140.29230.2736];Ii=[0.99331.03391.36521.2778];P2=(Qi./1).*IiQi=[0.24250.25190.33250.1731];Ii=[0.99031.02881.35780.7067];P3=(Qi./1).*IiQi=[0.23630.24170.36340.1585];Ii=[1.00021.02311.53800.6708];P4=(Qi./0.9999).*IiQi=[0.25340.25880.30820.1795];Ii=[0.99751.01871.21310.7067];P5=(Qi./0.9999).*IiQi=[0.25840.26110.33150.1491];Ii=[1.00181.01211.28500.5780];P6=(Qi./1.0001).*IiQi=[0.26840.27080.29880.1621];Ii=[0.99881.00771.11200.6033];P7=(Qi./1.0001).*IiQi=[0.25270.25490.27620.2161];Ii=[0.99921.00791.09180.8545];P8=(Qi./0.9999).*IiQi=[0.33260.33830.38130.2033];Ii=[0.99781.01481.14400.6099];P9=(Qi./1.2555).*IiQi=[0.26660.26960.29780.166];Ii=[1.00171.01311.11920.6238];P10=(Qi./1).*IiQi=[0.25580.25870.28580.1593];Ii=[1.00141.02341.25660.6353];P11=(Qi./0.9596).*IiQi=[0.23490.23760.26250.1463];Ii=[0.99901.02341.66100.5801];P12=(Qi./0.8813).*IiPij=0.22770.26270.46660.16440.21130.22890.39900.34960.24010.25920.45150.12230.23640.24730.55900.10630.25280.26370.37390.12690.25880.26420.42590.08620.26810.27290.33220.09780.25250.25690.30160.18470.26430.27340.34740.09880.26710.27310.33330.10360.26690.27590.37430.10550.26630.27590.49470.0963Pij反轉(zhuǎn)矩陣0.22770.21130.24010.23640.25280.25880.26810.25250.26430.26710.26690.26630.26270.22890.25920.24730.26370.26420.27290.25690.27340.27310.27590.27590.46660.39900.45150.55900.37390.42590.33220.30160.34740.33330.37430.49470.16440.34960.12230.10630.12690.08620.09780.18470.09880.10360.10550.0963權(quán)重氣象評估Ki=1.12141.18881.07311.14901.01731.03510.97100.99570.98390.97711.02261.1332綜合氣象評估指數(shù)Ii〔Pij反轉(zhuǎn)矩陣列項求和〕=1.07981.16761.02091.05800.98400.96920.93050.98830.94160.93950.97921.06591.10—1.30《070、》1.30優(yōu)良良差氣象資料分級檢驗1第二年二月份Ci=[1018.2,272.6,0.58,1.44];Di=[1027.8,277.5,0.89,2.5];Si=[1022,273.7,0.44,2.8];Ii=f4(Ci,Di,Si)Ii=1.00101.00491.63290.6776I=1.0791-----綜合氣象評估優(yōu)Ii=[1.00101.00491.63290.6776];Qi=Ii./4.3164Qi=0.23190.23280.37830.1570Qi求和為1Pi=[0.23190.23280.37830.1570];Ii=[1.00101.00491.63290.6776];Ni=Pi.*IiNi=0.23210.23390.61770.1064K=1.1901權(quán)重氣象評估良第二年七月份Ci=[993.6,28.7,0.62,1.45];Di=[1001,34.2,0.88,2];Si=[999.7,26.2,0.75,2.1];Ii=f4(Ci,Di,Si)Ii=0.99761.19580.98490.8109I=0.9973-------綜合氣象評估優(yōu)Ii=[0.99761.19580.98490.8109];Qi=Ii./3.9892Qi=0.25010.29980.24690.2033Q=1.0001