2005高教社杯全國大學(xué)生數(shù)學(xué)建模競賽

1、裝 訂 線論 文 封 面論文格式、用紙及論文提交要求1論文用紙、格式要求所有參賽隊(duì)必須首先從北京市教育委員會(huì)高教處網(wǎng)站上下載論文封面，本封面必須放在你的論文的上面，與你的論文一同裝訂。本封面為參賽隊(duì)參賽的合法標(biāo)志，無此封面的論文將被視為無效論文。北京市教育委員會(huì)高教處網(wǎng)址為：論文封面及正式論文一律使用A4復(fù)印紙，論文用紙由參賽校自行解決。紙的具體使用方法及要求如下：每一頁紙的上方留出30mm的空白即對(duì)應(yīng)裝訂線的上部，空白內(nèi)不允許書寫任何文字。在距離上方45mm處開始書寫正文。每一頁的下方中間部位要注明頁號(hào)。如需繪制圖形，應(yīng)盡量繪制在論文紙上。如果使用計(jì)算機(jī)繪制圖形，也應(yīng)將圖形繪制在論文的相應(yīng)部

2、位，或?qū)D形剪下，粘貼在論文紙的相應(yīng)部位。如使用計(jì)算機(jī)打印自己編寫的程序，可將打印的程序紙附在論文后面。2論文試卷的密封論文完成之后，按頁的順序整理好，并將學(xué)校名稱，參賽隊(duì)員姓名，指導(dǎo)教師姓名寫在本論文封面紙的第3頁上方的指定位置。然后，沿裝訂線用訂書器將論文裝訂好。裝訂好后，將論文封面的第一張紙即本頁紙沿裝訂線向上折起，折向論文的反面，并與論文的最后一頁的反面用膠水粘牢。此時(shí)第2頁紙成為首頁。3論文提交要求必須以兩種格式提交論文：論文的完全打印本，及內(nèi)容與打印本完全一致的相應(yīng)的電子版。兩種版本缺一不可，否那么將被視為無效參賽論文。打印本需按要求進(jìn)行裝訂密封；電子版用軟盤提供，文件名為三名參賽

3、隊(duì)員的姓名+AB、C或D。A、B、C、D為所選題目的題號(hào)，電子版論文可用Word，或LaTex可從網(wǎng)站下載模版編寫。注意：本統(tǒng)一論文封面紙共三頁，不能缺頁。試題下載 全國大學(xué)生數(shù)學(xué)建模競賽網(wǎng)站：高等教育出版社 以及網(wǎng)站：北京大學(xué)生數(shù)學(xué)建模與計(jì)算機(jī)應(yīng)用競賽競 賽 規(guī) 那么1數(shù)學(xué)建模競賽是一項(xiàng)開放式競賽，學(xué)生在競賽過程中可以在不違反競賽章程的前提下，查閱任何書籍或資料包括使用互聯(lián)網(wǎng)查閱資料；使用任何實(shí)驗(yàn)手段；做必要的社會(huì)調(diào)查；使用計(jì)算機(jī)及任何必要的軟件等。同一參賽隊(duì)的三名同學(xué)之間可以進(jìn)行隨意的討論。數(shù)學(xué)建模競賽以最大程度地發(fā)揮學(xué)生自身的創(chuàng)造力為競賽第一宗旨。每個(gè)參賽隊(duì)至多由三名隊(duì)員組成。2. 引用

4、的所有資料包括：雜志、書籍、網(wǎng)絡(luò)等一切公開發(fā)行，具有版權(quán)的資料。必須按公開發(fā)表論文的要求，注明：作者，書籍名稱與頁數(shù)或雜志名稱、期號(hào)、頁數(shù)，出版地，出版單位，出版日期。如引用的是譯著，同時(shí)還應(yīng)注明原著者，否那么將被視為剽竊。3如在建模的過程中使用了自己編寫的程序，需將程序全部打印，并將程序拷貝到一張軟盤中，在軟盤上注明學(xué)校、參賽隊(duì)員的姓名，并在論文最后的備注中注明“有程序軟盤字樣。然后由學(xué)校統(tǒng)一密封在一個(gè)小信封內(nèi)，并在信封上注明學(xué)校名稱，放入論文試卷密封袋。4每個(gè)參賽隊(duì)必須獨(dú)立完成論文的全部工作。參賽學(xué)生不得在任何地方、以任何方式、任何手段包括利用網(wǎng)絡(luò)與本參賽隊(duì)之外的任何人討論與競賽題目有關(guān)的

5、問題。一經(jīng)發(fā)現(xiàn)立即取消參賽資格，并通知學(xué)校處理。5密封后的論文的任何部位都不能出現(xiàn)可辨識(shí)參賽隊(duì)員和參賽隊(duì)所在學(xué)校的標(biāo)記。否那么將被取消參賽資格。6各參賽隊(duì)?wèi)?yīng)模范遵守?全國大學(xué)生數(shù)學(xué)建模競賽暨北京大學(xué)生數(shù)學(xué)建模與計(jì)算機(jī)應(yīng)用競賽章程?。任何參賽隊(duì)具有參賽隊(duì)間相互監(jiān)督與檢查競賽違紀(jì)現(xiàn)象的權(quán)利和義務(wù)。對(duì)發(fā)現(xiàn)其他隊(duì)的違章問題有向本校教務(wù)處、北京賽區(qū)組委會(huì)、北京市教育委員會(huì)高教處直接反映的權(quán)利和義務(wù)。7所有參賽隊(duì)必須使用統(tǒng)一的封面獲得方法見前頁，否那么被視為無效論文。8競賽設(shè)立爭議期制度。自初評(píng)結(jié)果公布之日起的15天內(nèi)，為初評(píng)結(jié)果爭議期。在爭議期內(nèi)競賽組委會(huì)受理一切與本次競賽中的違紀(jì)行為，及評(píng)卷異議有關(guān)的

6、投訴。投訴必須為書面形式。投訴信可發(fā)至北京市教育委員會(huì)高教處或北京大學(xué)生數(shù)學(xué)建模與計(jì)算機(jī)應(yīng)用競賽組委會(huì)。投訴信必須在爭議期內(nèi)發(fā)出以郵戳日期為限，過時(shí)不再受理。來信必須寫出真實(shí)學(xué)校名稱與個(gè)人姓名，匿名信恕不受理。初評(píng)結(jié)果不作為最終競賽成績。9競賽初評(píng)結(jié)果公布后，正式結(jié)果公布前，競賽專家評(píng)審組將根據(jù)實(shí)際需要對(duì)參賽隊(duì)進(jìn)行論文辯論。被約請參加論文辯論的參賽隊(duì)三名隊(duì)員，必須在指定的時(shí)間到達(dá)指定的地點(diǎn)參加論文辯論，否那么將被取消競賽資格。論文辯論的成績將被作為評(píng)定最終成績的重要參考依據(jù)。請?jiān)谕瓿勺珜懻撐牡娜抗ぷ髦?，在所選題號(hào)后面的括號(hào)內(nèi)畫“，不得涂改。選作題目代號(hào)：甲組：A ；B ；乙組：C；D 省市

7、代號(hào)：01 學(xué)校名稱： 中央民族大學(xué) 隊(duì)員姓名： 吳鋼 、 余思維 、 唐碧艷 指導(dǎo)教師姓名：指導(dǎo)組 裝 訂 線2005高教社杯全國大學(xué)生數(shù)學(xué)建模競賽編 號(hào) 專 用 頁賽區(qū)評(píng)閱編號(hào)由賽區(qū)組委會(huì)評(píng)閱前進(jìn)行編號(hào)：全國統(tǒng)一編號(hào)由賽區(qū)組委會(huì)送交全國前編號(hào)：全國評(píng)閱編號(hào)由全國組委會(huì)評(píng)閱前進(jìn)行編號(hào)： 雨量預(yù)報(bào)方法的評(píng)價(jià)模型摘要本文利用夾角余弦計(jì)算客觀性權(quán)重，用線性加權(quán)公式計(jì)算總評(píng)價(jià)值，建立了雨量預(yù)測方法優(yōu)劣的評(píng)價(jià)模型，對(duì)兩種雨量預(yù)測方案進(jìn)行綜合評(píng)價(jià)。在問題一中，首先我們確定站點(diǎn)有效影響范圍，找到該范圍內(nèi)的網(wǎng)格點(diǎn)，然后，處理該范圍內(nèi)網(wǎng)格點(diǎn)的預(yù)測數(shù)據(jù)，將其與該站點(diǎn)的實(shí)測數(shù)據(jù)進(jìn)行比擬。在此根底上，求出該站點(diǎn)有效

8、范圍內(nèi)網(wǎng)格點(diǎn)的預(yù)測數(shù)據(jù)與站點(diǎn)相應(yīng)的實(shí)測數(shù)據(jù)的綜合誤差，并采用模擬投票方法，評(píng)價(jià)出兩個(gè)雨量預(yù)測方法的優(yōu)劣。結(jié)論是方法一優(yōu)于方法二。對(duì)于問題二，我們將預(yù)測數(shù)據(jù)按照所給的等級(jí)標(biāo)準(zhǔn)轉(zhuǎn)化為7個(gè)等級(jí)，與站點(diǎn)的實(shí)際等級(jí)進(jìn)行了比擬。在此根底上，設(shè)立參加公眾損失度表征公眾對(duì)雨量預(yù)報(bào)過失的感受，進(jìn)而將公眾損失度與第一問的綜合誤差結(jié)合，并賦予權(quán)重，得出兩種方法的最終評(píng)價(jià)指標(biāo)，從而判斷出方法一優(yōu)于方法二。關(guān)鍵詞 : 夾角余弦 線性加權(quán) 損失度 Matlab問題重述雨量預(yù)報(bào)對(duì)農(nóng)業(yè)生產(chǎn)和城市工作和生活有重要作用，但準(zhǔn)確，及時(shí)地對(duì)雨量做出預(yù)報(bào)是一個(gè)十分困難的問題，廣受世界各國關(guān)注。我國某地氣象臺(tái)和氣象研究所正在研究6小時(shí)雨

9、量預(yù)報(bào)方法，即每天晚上20點(diǎn)預(yù)報(bào)從21點(diǎn)開始的4個(gè)時(shí)段 21點(diǎn)到次日3點(diǎn)，次日3點(diǎn)到9點(diǎn)，9點(diǎn)到15點(diǎn)，15點(diǎn)至21點(diǎn)在某些位置的雨量，這些位置位于東經(jīng)120度，北緯32度附近的53*47的等距網(wǎng)格點(diǎn)上同時(shí)設(shè)立91個(gè)不均勻的觀測站點(diǎn)實(shí)測這些時(shí)段的實(shí)際雨量?，F(xiàn)氣象部門提供了用兩種不同方法的預(yù)報(bào)數(shù)據(jù)和相應(yīng)的實(shí)測數(shù)據(jù)。請對(duì)雨量預(yù)測的兩種不同方法建立數(shù)學(xué)模型，具體要求如下：(1)請建立數(shù)學(xué)模型來評(píng)價(jià)兩種6小時(shí)雨量預(yù)報(bào)方法的準(zhǔn)確性；60毫米為大暴雨，大于60.1毫米為特大暴雨。假設(shè)按此分級(jí)向公眾預(yù)報(bào)，如何在評(píng)價(jià)方法中考慮公眾的感受？問題分析根據(jù)氣象部門提供的網(wǎng)格點(diǎn)和站點(diǎn)位置的經(jīng)緯度，可知該地區(qū)大致位于我

10、國浙江杭州到江蘇鹽城，安徽合肥到上海沿海的區(qū)域內(nèi)具體地圖如下，時(shí)間為6月中旬至7月底，正值該地區(qū)的梅雨時(shí)期。在此背景下，我們根據(jù)站點(diǎn)星號(hào)與網(wǎng)格點(diǎn)藍(lán)點(diǎn)的經(jīng)緯度用MATLAB軟件各點(diǎn)二維平面圖如下，不難看出該地區(qū)的地形因素是造成站點(diǎn)設(shè)置不均勻的主要原因。算出對(duì)于各站點(diǎn)的兩種不同方法的預(yù)測數(shù)據(jù)的綜合評(píng)價(jià)指標(biāo)，使得網(wǎng)格點(diǎn)的數(shù)據(jù)和各站點(diǎn)的數(shù)據(jù)能夠比擬，對(duì)于兩者的精確度，可以進(jìn)行模擬投票，比擬出兩種方法的優(yōu)劣；對(duì)于問題二，我們對(duì)降雨量的七個(gè)等級(jí)進(jìn)行數(shù)值轉(zhuǎn)換，即晴天為0，小雨為1，中雨為2，以此類推，相應(yīng)的也將各預(yù)測數(shù)據(jù)和預(yù)報(bào)數(shù)據(jù)進(jìn)行數(shù)值轉(zhuǎn)化，并引入損失度，再將二者進(jìn)行比擬，得出優(yōu)劣度。根本假設(shè)1. 各個(gè)站

11、點(diǎn)所得的實(shí)際數(shù)據(jù)都是精確的；2. 各點(diǎn)間的地理距離理想化為在二維平面距離；3. 緯度相差一單位的距離與經(jīng)度相差一單位的距離在二維平面上相等符號(hào)說明wc ： 代表某站點(diǎn)的有效范圍內(nèi)各網(wǎng)格點(diǎn)的綜合測量值與實(shí)際值之間的誤差 ： 代表各網(wǎng)格點(diǎn)的預(yù)測降雨量 ： 各站點(diǎn)的實(shí)測雨量 ： 網(wǎng)格點(diǎn)的橫坐標(biāo) ： 站點(diǎn)的橫坐標(biāo) ： 網(wǎng)格點(diǎn)的縱坐標(biāo) ：站點(diǎn)的縱坐標(biāo)pc ： 站點(diǎn)預(yù)測與實(shí)測差距給公眾帶來的綜合損失度f ： 個(gè)預(yù)測方法最后的評(píng)價(jià)結(jié)果q ： 各網(wǎng)格點(diǎn)對(duì)站點(diǎn)的權(quán)重模型建立與求解問題一：一 有效站點(diǎn)確實(shí)定本文定義有效范圍是實(shí)測站點(diǎn)能夠有效地測定的范圍。確定站點(diǎn)有效影響范圍，找到該范圍內(nèi)的網(wǎng)格點(diǎn)。有效范圍是以該站

12、點(diǎn)為中心以r為半徑的圓，其中r的計(jì)算公式為：其中dmin表示與該站點(diǎn)最近的網(wǎng)格點(diǎn)到該站點(diǎn)的距離；理由：用dmin加一個(gè)數(shù)是為了保證該站點(diǎn)有效范圍內(nèi)有網(wǎng)格點(diǎn)；網(wǎng)格點(diǎn)間的最小距離和最大距離的平方和的開方小數(shù)點(diǎn)后取一位有效數(shù)字，加上這樣一個(gè)數(shù)可保證有效范圍內(nèi)網(wǎng)格點(diǎn)數(shù)量比擬均勻。二 站點(diǎn)雨量預(yù)測值確實(shí)定處理該范圍內(nèi)網(wǎng)格點(diǎn)的預(yù)測數(shù)據(jù)，將其與該站點(diǎn)的實(shí)測數(shù)據(jù)進(jìn)行比擬，可以得到一個(gè)站點(diǎn)雨量預(yù)測的誤差?？紤]到該站點(diǎn)有效影響范圍內(nèi)網(wǎng)格點(diǎn)的個(gè)數(shù)不同且到該站點(diǎn)的距離也不同，我們通過各網(wǎng)格點(diǎn)到站點(diǎn)的距離計(jì)算出其對(duì)站點(diǎn)的權(quán)重，使距離越近的點(diǎn)權(quán)重越大，其計(jì)算公式為： 這里 ;在此根底上，求出該站點(diǎn)有效范圍內(nèi)網(wǎng)格點(diǎn)的預(yù)測數(shù)

13、據(jù)與站點(diǎn)相應(yīng)的實(shí)測數(shù)據(jù)的綜合誤差： 其中n表示該站點(diǎn)有效范圍內(nèi)網(wǎng)格點(diǎn)的個(gè)數(shù).三 雨量預(yù)測方法綜合評(píng)價(jià)模型的建立和求解 模型建立步驟如下：1 通過一天的4時(shí)段91站點(diǎn)的有效范圍內(nèi)網(wǎng)格點(diǎn)的預(yù)測數(shù)據(jù)與站點(diǎn)相應(yīng)的實(shí)測數(shù)據(jù)的綜合誤差構(gòu)成矩陣，為一天中第i個(gè)站點(diǎn)第j個(gè)時(shí)段的w值。2 確定該天各站點(diǎn)的兩個(gè)相對(duì)偏差矩陣： ；其中 3. 計(jì)算U，V的對(duì)應(yīng)列向量的夾角余弦得出初始權(quán)重即未歸一化的權(quán)重。歸一化后得到客觀性權(quán)重,它表示該天4個(gè)時(shí)段的客觀性權(quán)重。4建立本錢型矩陣， 5算出綜合評(píng)價(jià)值 6 結(jié)果 我們對(duì)6月的數(shù)據(jù)6月18日到6月28日進(jìn)行計(jì)算，結(jié)果可用Matlab求得 f1=; f2=; f1<f2,

14、所以認(rèn)為第一種預(yù)報(bào)6小時(shí)雨量的方法更為準(zhǔn)確；7模型檢驗(yàn)本文采用虛擬投票法對(duì)模型結(jié)果進(jìn)行檢驗(yàn)，虛擬各站點(diǎn)每天每時(shí)段都有工作人員根據(jù)站點(diǎn)實(shí)測數(shù)據(jù)與有效影響范圍內(nèi)網(wǎng)格點(diǎn)的預(yù)測數(shù)據(jù)對(duì)兩方法進(jìn)行投票，投票遵循的原那么是投數(shù)據(jù)最接近的方法，如兩者一樣那么對(duì)兩方法都投票。用Matlab進(jìn)行模擬6月6月18日到6月28日的情況，結(jié)果是： 方法一 2079票 方法二 1925票可見方法一確實(shí)比方法二更優(yōu)，問題一種的綜合評(píng)價(jià)指標(biāo)是合理的符合實(shí)際的。問題二一 原始數(shù)據(jù)轉(zhuǎn)換按照氣象部門的降雨量等級(jí)劃分，將原始數(shù)據(jù)都轉(zhuǎn)換為相應(yīng)的等級(jí)；轉(zhuǎn)換遵循的原那么為：二 等級(jí)比擬在問題一模型的根底上，處理有效范圍內(nèi)網(wǎng)格點(diǎn)的預(yù)測等級(jí)，

15、將其與該站點(diǎn)實(shí)等級(jí)進(jìn)行比擬，結(jié)果為：三 損失度與損失原那么確定引入氣象預(yù)報(bào)錯(cuò)會(huì)給公眾帶來的損失度作為評(píng)價(jià)標(biāo)準(zhǔn)來考慮公眾的感受，為此將p值轉(zhuǎn)換為損失度，轉(zhuǎn)換原那么如下：四 結(jié)合損失度計(jì)算該站點(diǎn)預(yù)測與實(shí)測差距給公眾帶來的綜合損失度：五 評(píng)價(jià)模型確實(shí)定直接引入問題一模型中各天時(shí)段的權(quán)重，通過某一天的4時(shí)段91站點(diǎn)的有效范圍內(nèi)網(wǎng)格點(diǎn)的預(yù)測等級(jí)與站點(diǎn)相應(yīng)的實(shí)測等級(jí)的綜合損失度構(gòu)成的矩陣可直接建立本錢型矩陣并算出綜合評(píng)價(jià)值， 其中， 六 結(jié)果我們對(duì)6月的數(shù)據(jù)6月18日到6月28日進(jìn)行計(jì)算，結(jié)果可用Matlab求得f1=; f2=;f1<f2,所以認(rèn)為第一種預(yù)報(bào)6小時(shí)雨量的方法更為準(zhǔn)確：七 模型檢驗(yàn)虛

16、擬各站點(diǎn)每天每時(shí)段都有公眾根據(jù)站點(diǎn)實(shí)測等級(jí)與有效影響范圍內(nèi)網(wǎng)格點(diǎn)的預(yù)測等級(jí)對(duì)兩方法進(jìn)行投票，投票遵循的原那么是投損失度最小的方法，如兩者一樣那么棄權(quán)。用Matlab進(jìn)行模擬6月6月18日到6月28日的情況，結(jié)果是： 方法一 613票 方法二 526票可見方法一確實(shí)比方法二更優(yōu)，在考慮公眾感受的情況下，用我們的模型依舊可以保證其穩(wěn)定有效性。模型評(píng)價(jià) 在模型的建立過程中，我們參考了大量有關(guān)地理知識(shí)的書籍和網(wǎng)站，因此我們的信息更廣泛更準(zhǔn)確；在對(duì)大量數(shù)據(jù)的處理方面，我們編寫的程序簡單，同時(shí)借助Matlab軟件，因此程序運(yùn)行速度快，可操作性強(qiáng)，能處理龐大的數(shù)據(jù)，模型易推廣。在模型假設(shè)中，我們假設(shè)中的2，

17、3點(diǎn)，即：各站點(diǎn)與被測點(diǎn)之間的地理距離理想化為在二維平面坐標(biāo)圖上的距離；緯度相差一單位的實(shí)際距離與經(jīng)度相差一單位的實(shí)際距離相等，雖然與實(shí)際情況有所差異，可能使結(jié)果不夠精確，但并不影響評(píng)價(jià)結(jié)果。參考文獻(xiàn)1 胡守信 李伯年，基于MATLAB的數(shù)學(xué)實(shí)驗(yàn),北京：科學(xué)出版社，20042 王建國 ，中國地圖冊，成都：成都地圖出版社編著，20034 互聯(lián)網(wǎng)三維地圖 ，5 中國數(shù)學(xué)建模網(wǎng) . 附件%畫出觀測站點(diǎn)與預(yù)測網(wǎng)格點(diǎn)位置圖x=A020618(:,3);y=A020618(:,2);plot(x,y,'rp') grid on xlabel('東經(jīng)度') ylabel(&#

18、39;北緯度')title('觀測站點(diǎn)與預(yù)測網(wǎng)格點(diǎn)位置圖') hold onplot(lon,lat,'b.') %將觀測站點(diǎn)與預(yù)測網(wǎng)格點(diǎn)的位置變?yōu)閺?fù)數(shù)形式 ZWZ=A020618(:,3)+A020618(:,2)*i; DWZ=lon+lat*i;%問題一模型求解%尋找網(wǎng)格點(diǎn)間的最小和最大距離Lr=zeros(52,47);Lc=zeros(53,46);for i=1:52 Lr(i,:)=abs(DWZ(i,:)-DWZ(i+1,:);endfor j=1:46 Lc(:,j)=abs(DWZ(:,j)-DWZ(:,j+1);endL1=min(

19、min(Lc);L2=max(max(Lc);%比擬方法(一)各天各時(shí)段各站點(diǎn)實(shí)測值與其有效范圍內(nèi)網(wǎng)格點(diǎn)預(yù)測值for day=618:628 wc=zeros(91,4); for i=1:4 disp('i=',i); for n=1:91 S=DWZ-ZWZ(n); R=abs(S); Rmin=min(min(R)+0.2; ok=0; sd=0; wcd=zeros(20,10); for r=1:53 for c=1:47 if R(r,c)<=Rmin ok=ok+1; M=eval('f',num2str(day),num2str(i),&#

20、39;_dis',num2str(1); N=eval('A020',num2str(day); dst=abs(M(r,c)-N(n,i+3); wcd(ok,1)=dst; wcd(ok,2)=1/(R(r,c)+0.001);%為防止出現(xiàn)站點(diǎn)與網(wǎng)格點(diǎn)重合而出現(xiàn)分母為零，對(duì)Rr,c進(jìn)行加0.001處理下同 wcd(ok,3)=r; wcd(ok,4)=c; wcd(ok,5)=sum(wcd(:,2); end end end for k=1:ok wcd(k,6)=1/(R(wcd(k,3),wcd(k,4)+0.001)/wcd(ok,5); wcd(k,7)=

21、wcd(k,1)*wcd(k,6); wcd(k,8)=sum(wcd(:,7); wc(n,i)=wcd(ok,8); end endenddayc=day-617;gwc1(:,:,dayc)=wc;end%問題一模型求解%比擬方法(二)各天各時(shí)段各站點(diǎn)實(shí)測值與其有效范圍內(nèi)網(wǎng)格點(diǎn)預(yù)測值for day=618:628 wc=zeros(91,4); for i=1:4 disp('i=',i); for n=1:91 S=DWZ-ZWZ(n); R=abs(S); Rmin=min(min(R)+0.2; ok=0; sd=0; wcd=zeros(20,10); for r

22、=1:53 for c=1:47 if R(r,c)<=Rmin ok=ok+1; M=eval('f',num2str(day),num2str(i),'_dis',num2str(2); N=eval('A020',num2str(day); dst=abs(M(r,c)-N(n,i+3); wcd(ok,1)=dst; wcd(ok,2)=1/(R(r,c)+0.001); wcd(ok,3)=r; wcd(ok,4)=c; wcd(ok,5)=sum(wcd(:,2); end end end for k=1:ok wcd(k,6)

23、=1/(R(wcd(k,3),wcd(k,4)+0.001)/wcd(ok,5); wcd(k,7)=wcd(k,1)*wcd(k,6); wcd(k,8)=sum(wcd(:,7); wc(n,i)=wcd(ok,8); end endenddayc=day-617;gwc2(:,:,dayc)=wc;end%計(jì)算兩方法在各天的綜合評(píng)價(jià)值s1=zeros(11,2);for way=1:2 for i=1:11 E=eval('gwc',num2str(way); m1=ones(91,1)*max(E(:,:,i); m2=ones(91,1)*min(E(:,:,i);

24、n=ones(91,1)*range(E(:,:,i); u=abs(E(:,:,i)-m1)./n; v=abs(E(:,:,i)-m2)./n; d=dot(u,v); nu=norm(u(:,1),norm(u(:,2),norm(u(:,3),norm(u(:,4); nv=norm(v(:,1),norm(v(:,2),norm(v(:,3),norm(v(:,4); w1=d./(nu.*nv); w=w1./sum(w1); B=E(:,:,i)./m1; F=B*w' s1(i,way)=sum(F); endend%問題一模擬工作人員投票檢驗(yàn)?zāi)Ｐ秃侠硇詔p1=0;tp

25、2=0;for day=1:11 for i=1:91 for j=1:4 if gwc1(i,j,day)>gwc2(i,j,day) tp2=tp2+1; elseif gwc1(i,j,day)<gwc2(i,j,day) tp1=tp1+1; else tp1=tp1+1; tp2=tp2+1; end end endend%問題二模型求解%對(duì)實(shí)測數(shù)據(jù)進(jìn)行等級(jí)替換for day=618:628 N=eval('A020',num2str(day); for i=1:4 for n=1:91 if N(n,i+3)=0; elseif N(n,i+3)=1;

26、elseif N(n,i+3)=2; elseif N(n,i+3)=3; elseif N(n,i+3)=4; elseif N(n,i+3)=5; else N(n,i+3)=6; end end end NNN(:,:,day-617)= N;end%對(duì)網(wǎng)格點(diǎn)預(yù)測數(shù)據(jù)進(jìn)行等級(jí)替換for day=618:628 for i=1:4 disp('i=',i); for way=1:2 M=eval('f',num2str(day),num2str(i),'_dis',num2str(way); for r=1:53 for c=1:47 if

27、M(r,c)=0; elseif M(r,c M(r,c)=1; elseif M(r,c)=2; elseif M(r,c)=3; elseif M(r,c)=4; elseif M(r,c)=5; else M(r,c)=6; end end end szh=(day-618)*8+(i-1)*2+way; MMM(:,:,szh)=M; end endend %比擬方法(一)各天各時(shí)段各站點(diǎn)實(shí)測值等級(jí)與其有效范圍內(nèi)網(wǎng)格點(diǎn)預(yù)測值等級(jí)差異給公眾的損失度for day=618:628 wc2=zeros(91,4); for i=1:4 disp('i=',i); for n=

28、1:91 S=DWZ-ZWZ(n); R=abs(S); Rmin=min(min(R)+0.2; ok=0; wcd2=zeros(20,10); for r=1:53 for c=1:47 if R(r,c)<=Rmin ok=ok+1; szh=(day-618)*8+(i-1)*2+1; djc=MMM(r,c,szh)-NNN(n,i+3,day-617); if djc<0 qz=-djc*2; else qz=djc; end wcd2(ok,1)=qz; wcd2(ok,2)=1/(R(r,c)+0.001); wcd2(ok,3)=r; wcd2(ok,4)=c; wcd2(ok,5)=sum(wcd2(:,2); end end end for k=1:ok wcd2(k,6)=1/(R(wcd2(k,3),wcd2(k,4)+0.001)/wcd2(ok,5); wcd2(k,7)=wcd2(k,1)*wcd2(k,6); wcd2(k,8)=sum(wcd2(:,7); wc2(n,i)=wcd2(ok,8); end endenddayc=day-617;gwc3(:,:,dayc)=wc2;end%比擬方法(二)各天各時(shí)段各站點(diǎn)實(shí)測值等級(jí)與其有效范圍內(nèi)網(wǎng)格點(diǎn)預(yù)測值等級(jí)差異給公眾的損失度for day=6