水資源短缺風(fēng)險綜合評價

2023高教社杯全國大學(xué)生數(shù)學(xué)建模競賽承諾書我們仔細(xì)閱讀了中國大學(xué)生數(shù)學(xué)建模競賽的競賽規(guī)則.我們完全明白，在競賽開始后參賽隊員不能以任何方式（包括、電子郵件、網(wǎng)上咨詢等）與隊外的任何人（包括指導(dǎo)教師）研究、討論與賽題有關(guān)的問題。我們知道，抄襲別人的成果是違反競賽規(guī)則的,如果引用別人的成果或其他公開的資料（包括網(wǎng)上查到的資料），必須按照規(guī)定的參考文獻(xiàn)的表述方式在正文引用處和參考文獻(xiàn)中明確列出。我們鄭重承諾，嚴(yán)格遵守競賽規(guī)則，以保證競賽的公正、公平性。如有違反競賽規(guī)則的行為，我們將受到嚴(yán)肅處理。我們授權(quán)全國大學(xué)生數(shù)學(xué)建模競賽組委會，可將我們的論文以任何形式進(jìn)行公開展示（包括進(jìn)行網(wǎng)上公示，在書籍、期刊和其他媒體進(jìn)行正式或非正式發(fā)表等）。我們參賽選擇的題號是（從A/B/C/D中選擇一項填寫）：我們的參賽報名號為（如果賽區(qū)設(shè)置報名號的話）：所屬學(xué)校（請?zhí)顚懲暾娜簠①愱爢T(打印并簽名)：1.2.3.指導(dǎo)教師或指導(dǎo)教師組負(fù)責(zé)人(打印并簽名)：稠鹽半僻籌蘆禍痛同日期：券他年惰緣月慶蓬日許賽區(qū)評閱編炭號（由賽區(qū)訓(xùn)組委會評閱殿前進(jìn)行編號擴(kuò)）：語2023杯高教社杯全紹國大學(xué)生數(shù)教學(xué)建模競賽賴編號專線用頁洗賽區(qū)評閱編艙號（由賽區(qū)仙組委會評閱孝前進(jìn)行編號翼）：松賽區(qū)評閱記子錄（可供賽丹區(qū)評閱時使壺用）：到評攜閱海人授評側(cè)分芽備辰注跳全國統(tǒng)一編共號（由賽區(qū)剪組委會送交穿全國前編號稠）：慌全國評閱編腸號（由全國騎組委會評閱此前進(jìn)行編號綁）：警水資源短缺賞風(fēng)險綜合評唯價摘要留水資源短缺鈴問題是影響糕我國發(fā)展的自重大問題，晉本文針對水三資源短缺風(fēng)楊險問題找出嘆了主要風(fēng)險摘因子，建立雞了水資源短檢缺風(fēng)險評價過模型，對水商資源短缺風(fēng)魄險進(jìn)行等級泥劃分，并提芹出淚相應(yīng)的有效妨措施規(guī)避風(fēng)裕險陜。右對于問題一折，我們建立幸主成分和灰眠色關(guān)聯(lián)度分葛析模型，分窯析附表和相革關(guān)資料，先件確立了閃北京市水資禿源短缺風(fēng)險書的風(fēng)險因素賞主要包括自授然因素，即誕降雨量和常超住人口，和幫社會因素，炮即農(nóng)業(yè)用水逝，工業(yè)用水狼，第三產(chǎn)業(yè)炕及生活其他臟用水，污水牽處理率，城漠市綠化覆蓋茅率。然后利麗用主成分分仇析撇得到個各個天因子的貢獻(xiàn)每率碎，再利用灰社色關(guān)聯(lián)度分盲析，得到各籠個因子與缺活水量的槍關(guān)聯(lián)度的大勺小申，基本與主晚成分分析一末致，最后得屑到主要風(fēng)險客因子。冊對于問題二替，仁我們用綜合板評價蜻的編模型,遭對模型一所呈確定的主要閱風(fēng)險因子冰做相應(yīng)卸高優(yōu)和低優(yōu)什指標(biāo)處理邁，工并塵對北京市水號資源短缺進(jìn)普行風(fēng)險等級奔劃分。最后蝴對抵主要風(fēng)險因示子潤進(jìn)行城調(diào)控，刻來餡降低風(fēng)險等井級變。沉對于問題三銀，我們建立綱休鎮(zhèn)模型，要對各北京市未來材兩年水資源某的短缺風(fēng)險彩進(jìn)行預(yù)測，樂我們通過對晨主要風(fēng)險因朝子付進(jìn)行色預(yù)測，溪并對預(yù)測模斧型進(jìn)行麗后驗差檢驗飾，掩然后再千用脊模型身，給未來的乓兩年劃分風(fēng)捉險等級。通對于問題四忠，我們通過京分析上面的冬數(shù)據(jù)和查找冠相關(guān)資料，竹給北京市水贈行政主管部掩門寫一份建隸議報告?？瓣P(guān)鍵詞：主慚成分分析閘灰色關(guān)聯(lián)度舒分析啄模型械后模型怎后驗差檢驗壘劃一、隨問題重述鮮新中國成立拒以來，黨和抖政府領(lǐng)導(dǎo)全性國人民進(jìn)行傲了大規(guī)模的逃水利建設(shè)，間取得了巨大萄成就，但水勵的問題仍然啦困擾著我國拌經(jīng)濟(jì)社會的走發(fā)展?！?器998年中拴國水資源公位報》中的數(shù)演據(jù)表明：中鏈國面臨的水礙資源形勢依頭然嚴(yán)峻。水顯的問題已經(jīng)圾成為制約我舒國絕大部分津地區(qū)經(jīng)濟(jì)社樓會可持續(xù)發(fā)失展的重要因奉素。默水資源，是捐指可供人類孩直接利用，備能夠不斷更運(yùn)新的天然水摟體。主要包泥括陸地上的剪地表水和地比下水。艙風(fēng)險，是指詳某一特定危茄險情況發(fā)生毒的可能性和集后果的組合茅。剖水資源短缺娛風(fēng)險，泛指貸在特定的時皺空環(huán)境條件冒下，由于來石水和用水兩獄方面存在不械確定性，使段區(qū)域水資源適系統(tǒng)發(fā)生供偶水短缺的可假能性以及由蒼此俱可能輸產(chǎn)生種相關(guān)的快的損失。傲北京是世界喘上水資源嚴(yán)鉗重缺乏的大孤都市之一，蹈目前，北京敵水資源年均您用水總量達(dá)魂36億立方亭米，而年均厲水資源總量島僅有21億鬧立方米，缺致口達(dá)15億醬立方米。北擇京的年均水赴資源量平攤奶到每個人身色上，甚至不害足蝦100立方寨米舊。擔(dān)為全國人均副的1/8，且世界人均的慕1/30，聽屬重度缺水紡地區(qū)，進(jìn)為了維持到擦現(xiàn)在，北京哭用盡了各種穿方法，除了擁尋找多種水攻源悟“免開源塞”含，還通過調(diào)胞整用水結(jié)構(gòu)罩、再生水利蠢用等辦法督“嶄節(jié)流?！鳖D。用三句話他概括，就是營農(nóng)業(yè)用水負(fù)湊增長，工業(yè)仔用水零增長腹，生活用水拉適度增長。的附表中所列朋的數(shù)據(jù)給出堡了1979沸年至200退0年北京市隱水資源短缺煮的狀況。北灰京市水資源總短缺已經(jīng)成陳為影響和制閣約首都社會傭和經(jīng)濟(jì)發(fā)展哀的主要因素仙。政府采取海了一系列措僚施,如南癥水北調(diào)工程破建設(shè),建暢立污水處理皆廠,產(chǎn)業(yè)結(jié)剃構(gòu)調(diào)整等。妥但是，氣候顆變化和經(jīng)濟(jì)倡社會不斷發(fā)替展，水資源翠短缺風(fēng)險始磚終存在。如偉何對水資源范風(fēng)險的主要抗因子進(jìn)行識姿別，對風(fēng)險膀造成的危害皺等級進(jìn)行劃蒼分，對不同煤風(fēng)險因子采門取相應(yīng)的有敬效措施規(guī)避咱風(fēng)險或減少蔥其造成的危堆害，這對社畫會經(jīng)濟(jì)的穩(wěn)升定、可持續(xù)西發(fā)展戰(zhàn)略的名實施具有重宜要的意義。胖《北京想2023崗統(tǒng)計年鑒》遇及市政統(tǒng)計烏資料墓提供了北京吐市水資源的鄰有關(guān)信息。浮利用這些資狐料和我們可如獲得的其他裳資料，討論矮了以下問題謠：變問題一：匙影響水資源念的因素很多孫,例如：氣菊候條件、水攻利工程設(shè)施河、工業(yè)污染侍、農(nóng)業(yè)用水嶄、管理制度鍋，人口規(guī)模咐等。我們評魂價判定北京謠市水資源短忙缺風(fēng)險的主拌要風(fēng)險因子睬是什么。領(lǐng)問題二：建扶立一個數(shù)學(xué)涌模型對北京撒市水資源短貪缺風(fēng)險進(jìn)行角綜合評價，足作出風(fēng)險狗等級劃分并春陳述理由。株并通過相關(guān)培數(shù)據(jù)對主要騰風(fēng)險因子,派如何進(jìn)行調(diào)喉控，使得風(fēng)尚險降低。膛問題三：建蠶立預(yù)測模型詠，對北京市終未來兩年水源資源的短缺愛風(fēng)險進(jìn)行預(yù)扇測，并提出鳳應(yīng)對措施。避問題四：以餃北京市水行獄政主管部門在為報告對象載，寫一份建太議報告。愁二柳、票基本假設(shè)贈（耽1慮）噴用水量是農(nóng)腰業(yè)用水，工沒業(yè)用水，第釋三產(chǎn)業(yè)及生薯活其他用水袋的總和。蹲模型的分析救與建立均以耍北京市為基撓礎(chǔ)。逼（2）切由于水資源腸的不穩(wěn)定性峽，在計算中捕排除某一差不異較大的數(shù)血，對模型無捧影響。模型價的分析與建火立均以北京蒸市為基礎(chǔ)。篩（3）誦引起水資源沿短缺的因素弊除我們選擇厲的7個因子院外，其他的者因子影響都盆很小，可以滋忽略不計。械（4）莖在一年時間湊內(nèi)，各風(fēng)險怨因子對北京農(nóng)市水資源短律缺的影響是儉一定的。懂（5）珠模型的分析憲與建立均以造北京市為基鑄礎(chǔ)。略三鳳、腫符號說明：分辨率；傍：滑矩陣相關(guān)系叔數(shù)且；愈:評價對展象；墻:織傘評價指標(biāo)排刻成；鄰:燭各組數(shù)據(jù)壇的秩休；脹:懂相應(yīng)權(quán)重貢；蒸:暢畏稱為發(fā)展灰斑數(shù)氏；琴:直道稱撤為內(nèi)生控制煩灰數(shù)鴉四象、盟問題分析偵對于問題一兇，市我們建立的天是用主成分貨和灰色關(guān)聯(lián)寒度分析法對朵水資源短缺遇的風(fēng)險因子矮進(jìn)行綜合評秋價。最通過收集資內(nèi)料和分析數(shù)謎據(jù)，得到影躍響水資源短斥缺風(fēng)險因素膝，主要有社赴會因素和自姜然因素，社揚(yáng)會因素又可述以細(xì)分為螺農(nóng)業(yè)用水，舅工業(yè)用水，虛第三產(chǎn)業(yè)及抓生活其他用些水，污水處換理率，城市健綠化覆蓋率需。自然因素亦可以細(xì)分為捉降雨量和常享住人口。染主成分分析反中得到每個執(zhí)因子在各個損主成分中的哨貢獻(xiàn)率，根嬸據(jù)貢獻(xiàn)率的茅大小提取主包要因子。然鍬后再結(jié)合灰路色關(guān)聯(lián)度分辱析來求得各另個因子與缺假水量的關(guān)聯(lián)冒度，對主成秧分分析中得藥到的主要風(fēng)薦險因子進(jìn)行堡檢驗。情對于問題二金，著基于問題一滔所確立的主粘要風(fēng)險因子川，和對知低優(yōu)指標(biāo)相感應(yīng)處理，利脈用暑模型觸對哲北京市水資遺源短缺風(fēng)險恢進(jìn)行綜合評哨價。根據(jù)奸對應(yīng)的概率鄉(xiāng)單位值進(jìn)行秀相應(yīng)低等級劃分置。通過控制痛影響塑北京市水資沉源短缺遷的主要思風(fēng)險嚴(yán)因子的權(quán)重沒，來改變風(fēng)略險等級。皺對于問題三那，我們用掀顯恢模型握，分別預(yù)測跨了未來兩年罩主要風(fēng)險因異子，即污水能處理率，農(nóng)沙業(yè)用水，第從三產(chǎn)業(yè)及生煎活其他用水風(fēng)，降雨量，療常住人口，貿(mào)在利用問題潤二的對模型抬，得到未來如兩年的風(fēng)險棋等級。殘對于問題四耽，通過鉤對水資源短應(yīng)缺風(fēng)險因子事的分析，和胞相關(guān)的預(yù)測渾，我們可以忍采取一些措臣施，降低水喝資源短缺的什風(fēng)險等級。理五從、吸模型的建立疊與求解抵5控.1敗問題一模型畢的建立與求膏解絲通過查找相智關(guān)資料，得防到影響水資盲源短缺風(fēng)險剪因素，主要孕有社會因素江即為萬農(nóng)業(yè)用水，限工業(yè)用水，戲第三產(chǎn)業(yè)及楊生活其他用怎水，污水處藏理率，城市紐綠化覆蓋率式。和自然因賄素即降雨量煉和常住人口灘。如圖1所鑰示。水資源短缺風(fēng)險因素社會因素水資源短缺風(fēng)險因素社會因素自然因素農(nóng)業(yè)用水污水處理率常住人口城市綠化覆蓋率降雨量第三產(chǎn)業(yè)及其他用水圖1風(fēng)險因子分布圖工業(yè)用水插漂主成分分析下的研究原理啊在數(shù)據(jù)分析慰工作中，?，F(xiàn)常需要把很乞復(fù)雜的數(shù)據(jù)拔集簡化，即賭將飽評個指標(biāo)所構(gòu)豈成的徹叔維簡化為一惑位系統(tǒng)。主雙成分分析就復(fù)是多個指標(biāo)做化為少數(shù)幾膚個綜合指標(biāo)掉的一種統(tǒng)計噸分析方法。茂在多指標(biāo)（痕變量）研究序中，往往由陳于變量個數(shù)瓦太多，且彼停此之間存在要著一定的相則關(guān)性，因而汽使得所觀測父的數(shù)據(jù)在一編定程度上有聞信息的重疊塊。當(dāng)變量較興多時，在高閣維空間中研狂究樣本的分死布規(guī)律就更吐復(fù)雜。主成齡分分析采取躁一種降維的訊方法，找出廢幾個綜合因窩子盡可能地孫反映原來變渴量的信息量公，而且彼此超之間互不相懲關(guān)，從而達(dá)輕到簡化的目油的。責(zé)主成分分析扣的計算步驟欲如下：設(shè)觀失測樣本矩陣咽為（新為樣本數(shù)，倍娘為變量數(shù)）共：商為使該樣本扶集在降維中競所引起的平妹方誤差最小賤，必須進(jìn)行環(huán)兩方面的工楚作：一是進(jìn)扒行坐標(biāo)變換讓，即用雅可暮比方法求解剃正交變換矩尸陣；二是選乏取忠個主成分.戲?qū)⒃紨?shù)據(jù)虎進(jìn)行標(biāo)準(zhǔn)化患處理，即對孔樣本集中元丟素：作變換襲幻（i=1,淺2，茶…濕，n；k=咬1,2，涉…嬌，供蝴）慌作變換，即，其中：椒為參與評價漆的指標(biāo)個數(shù)模。主成分分稼析的明顯特穿征是每個主普分量依賴于械測量初始變藥量所用的尺章度，當(dāng)尺度汪改變時，會爬得到不同的俯特征值?？藛璺@個困難陸的方法是對限初始變量進(jìn)液行以上標(biāo)準(zhǔn)偵化處理，使蹦其方差為1診.洞計算樣本矩訊陣的相關(guān)系精數(shù)矩陣。暢對應(yīng)相關(guān)系倉數(shù)矩陣，計續(xù)算特征值以糞及各個主成橡分的貢獻(xiàn)率剖。再把貢獻(xiàn)嫌率大的因子咱找出來，從癥而起到篩選右的作用。貫宇主成分分奶析法求解過獅程斥在附錄宗1籮中給出了1扣979到2否010年北旋京水資源短疫缺風(fēng)險評價逝的原始數(shù)據(jù)蛋，欣我們利用深釣將這些數(shù)據(jù)壯進(jìn)行標(biāo)準(zhǔn)差驚標(biāo)準(zhǔn)化處理距后計算其相活關(guān)系數(shù)矩陣縱（見表砍1妥）盜酸斑瘋表1相皂關(guān)系數(shù)矩陣豈污水處理率合（%）陽農(nóng)業(yè)用水（染億立方米）腰工業(yè)用水（嗽億立方米）青第三產(chǎn)業(yè)及穿生活其他用憑水（億立方盡米）榨降雨量（毫嚇米）晚城市綠化覆屯蓋率（%）澡常住人口（犬萬人）散污水處理率殊（%）沾1.000蟲-0.64幻8沖-0.92弄4保0.877鐵-0.38剝0民0.867高0.946高農(nóng)業(yè)用水（潤億立方米）旱-0.64升8醫(yī)1.000烘0.506物-0.71片0日-0.01艇9鄙-0.71敢2色-0.72額2燒工業(yè)用水（逮億立方米辯-0.92遮4醒0.506址1.000碰-0.81指3壩0.397趕-0.81軍4妄-0.87介9迷第三產(chǎn)業(yè)及見生活其他用嘴水（億立方稼米環(huán)0.877絨-0.71敬0嫁-0.81偶3飲1.000堵-0.29譽(yù)2采0.964襪0.958任降雨量（毫害米）豐-0.38結(jié)0斥-0.01活9慰0.397首-0.29帶2珠1.000槳-0.31黃8講-0.30塑0該城市綠化覆硬蓋率（%）嫂0.867論-0.71懇2且-0.81殃4變0.964膽-0.31棍8趙1.000乎0.946枯常住人口（泥萬人）登0.946尤-0.72版2辱-0.87尿9墾0.958啦-0.30棚0潤0.946傘1.000惡再由相關(guān)系慧數(shù)矩陣計算爬特征值和各辰個主成分的傾貢獻(xiàn)率和累暫積貢獻(xiàn)率如靈表2，表3宴所示。昂芳員爹茶表2特征萍值和貢獻(xiàn)率訓(xùn)和累積貢獻(xiàn)涼率水成分橡初始特征值獅提取平方和技載入拼合計嘗方差的%滋累積%理合計錦方差的%乖累積%饅1鈔5.243杏74.90姑0榴74.90鑰0樣5.243貨74.90役0挨74.90勾0挎2排1.050哨14.99蝕4洲89.89舌4昏1.050轉(zhuǎn)14.99反4電89.89淺4斑3漿0絹.381連5.440撈95.33我4紐0撐.381懇5.440技95.33維4化4呢0美.207撒2.958稅98.29噴2箱0陶.207瘋2.958喚98.29苦2原5岔0斃.065縫0救.933竭99.22究4漿0擠.065躺0活.933銜99.2緩24悼6困0監(jiān).034飛0宴.482玉99.70擔(dān)6金7嚴(yán)0箏.021贏0資.294戒100.0堅00撲嶼匯腥晚表3成分暢矩陣粗1挨2委3肢4簽5傷污水處理率慘（%）叼.958升-.067劈.134怠-.167寬.157欣農(nóng)業(yè)用水（浮億立方米）聲-.742芝-.486桑.417舊.194通.042異工業(yè)用水（疼億立方米）收-.905耗.164丸-.308凳.186她.155制第三產(chǎn)業(yè)及散生活其他用護(hù)水（億立方希米）脆.962略.075祥-.008史.228那.021身降雨量（毫雷米）因-.374帥.878犯.298給.023諸.004畝城市綠化覆綁蓋率（%）匯.960昂.054朗-.039教.228傳-.085孩常住人口（掘萬人）絞.985奉.058蚊.060廣.044突.085油由表廊2、3湯可以看出，桑第一、科第二主成分惠的累計貢獻(xiàn)父率已達(dá)到了潑89.89瘦4%，壯我們在第一駐和第二主成垂分中選擇因樸子，各個因怠子系數(shù)的大稻小可以反映沿因子對主成鄙分的貢獻(xiàn)率雅，所以我們?nèi)邕x擇了農(nóng)業(yè)濁用水，第三將產(chǎn)業(yè)及生活屈其他用水，馬降雨量，令城市綠化覆線蓋率，籮常住人口，聞作為水資源仔短缺的主要駛風(fēng)險因子。轎季任灰色關(guān)聯(lián)屢度的原理拒和求解原灰色關(guān)聯(lián)度折分析是分析駱系統(tǒng)中各因騙素關(guān)聯(lián)的程聰度的方法，紛在計算關(guān)聯(lián)甩度之前，需戚先計算關(guān)聯(lián)碎系數(shù)。關(guān)聯(lián)系數(shù)：游則關(guān)聯(lián)系數(shù)別定義為任：懲依其中：宿織竟分豐為第在個點融剝和閣挨跑的絕對誤差郊；燦為兩級最小嫌差；陸為兩級最大子差；悄洗稠街稱為分辨率符，0<秩桂<1,若衫集越小，關(guān)聯(lián)磚系數(shù)間差異低越大，區(qū)分宮能力越強(qiáng)。況一般取該升=0.5控。澇對單位不一疏，初值不同歪的序列，在郊計算相關(guān)系拜數(shù)前應(yīng)首先者進(jìn)行初始化槳，即將該序鈴列所有的數(shù)松據(jù)分別除以知第一個數(shù)據(jù)遲。生利用上述方寧法，我們?nèi)鳖i水量來衡量板水資源短缺夸程度，所以況把缺水量作外為參考序列寶，各個因子找作為比較序顯列，得到各沉個因子與缺她水量的關(guān)聯(lián)凍度QUOTEr翠導(dǎo)（如表蜂4典所示）。仰表好4定各因子關(guān)固聯(lián)度值門因子典污水處理率擋農(nóng)業(yè)用水教工業(yè)供水純第三產(chǎn)業(yè)及滴生活其他用哪水弱降雨量棗城市綠化率緞常住人口鹽0.550蹦1昆0.652并0炸0.589榆2壺0.693露7洪0.650業(yè)2汽0.699疼1格0.704縫2無根據(jù)關(guān)聯(lián)度退的大小，我販們同樣可以調(diào)選擇農(nóng)業(yè)用油水，第三產(chǎn)撤業(yè)及生活其恨他用水，降荒雨量，城市握綠化率，常歡住人口。梅傅結(jié)果說明遷用主成分分卸析評價影響拳水資源短缺從的主要風(fēng)險毛因子和用灰單色關(guān)聯(lián)度分赤析方法對影歇響水資源短狂缺的主要風(fēng)媽險因子，其丘結(jié)果基本一觀致，說明對吳于北京水資拖源短缺的主泥要風(fēng)險因子潔為將農(nóng)業(yè)用水，松第三產(chǎn)業(yè)及離生活其他用殼水，降雨量演，城市綠化匯率，常住人犧口。炒5.2問題添一模型的建漢立與求解議候QUOTERSR賞抬原理炎秩和比(測)指的是表獅中行(或列泊)秩次的平亭均值，是一靠個非參數(shù)計饞量的綜合指涌數(shù)，具有惡0橡～種1區(qū)間連續(xù)捆變量的特征刪，其所有比竟較組秩和比默之和為嚇子暴；如果編秩枝不按照經(jīng)典方的秩變換方長法，各組秩能和比的合計劍可能不為裳編漲張。墻其基本思想銜是在一個乞行（醫(yī)評價對象）螞列（滑個評價指標(biāo)增或等級）矩扣陣中，通過勿秩轉(zhuǎn)換，獲揪得無量綱的戴統(tǒng)計量QUOTERSR同吸，以QUOTERSR晴夸值對評深價對象的優(yōu)草劣進(jìn)行排序聽，進(jìn)而根據(jù)灑比較組數(shù)的錦多少，進(jìn)行慣分檔處理（談比濫組夫數(shù)較多）或按進(jìn)行錘QUOTERSR循村平方根反正姜弦變換值可桂信區(qū)間處理辛（比較組數(shù)飾較少）請。卷霜悶?zāi)Ｐ偷那蠼鈰穑?）確定慶評價對象和義綜合評價因撫子和編秩貪將布個評價對象般和遵個評價指標(biāo)許排成策行粗列的原始數(shù)樂據(jù)矩陣。高術(shù)優(yōu)指標(biāo)從小撕到大編秩,擇低優(yōu)指標(biāo)從掩大到小編秩卻,相同者取薦平均值。謎本題我們確狂定5個主要顛影響因素：泥農(nóng)業(yè)用水、赴第三產(chǎn)業(yè)及劃其他用水、棗常住人口、只污水處理率未、降雨量。降從附錄表1燦中提取它們霞對應(yīng)的32低年原始數(shù)據(jù)通。膏（2）計算功和口通過血軟件計算各腹組數(shù)據(jù)的秩乏（QUOTER沃攻）和相應(yīng)權(quán)糊重(土)，公式如漸下：宗斜鈴馬欄QUOTERSR=1ｍRm?n愁垮拉團(tuán)另給（獨膛挖錯弱臥QUOTEWRSR=1ｍR?Wm?n跌顛遭連另必通奴索（2舞床應(yīng)用上一步舒的5組原始季數(shù)據(jù)，通過頌公式（1）冊、（2）得蒜到如下結(jié)果悟：=[槳玩大0.345弟2刷蜘0.508吩0溜0.436史9太0.415后6損0.465鼓5性0.390著1役想圍淋雜0.200蠟9糠行0.436苦9璃靈0.303首0精膀桃忘0.483買1呆0.574逝2償0.497袍2酷0.545規(guī)8帽晝0.58翁13嶺0.58其33訪0.53神34樹0.56備18勵0.49狂98境絹0.531缺4劇益0.482殊9尚0.619辯9蜻0.595哨9花0.606冒1蛾0.558徹3怖抗0.591圈0少足0.573澇1稿0.599齊7覺0.621擺1學(xué)0.590儲6亭0.567碼5同崗0.608貿(mào)7叛士0.591篩4盡籍盾蔑]背（3）QUOTEWRSR單習(xí)的分布QUOTEWRSR各遙的分布是指穩(wěn)用概率單位搶讓QUOTEprobit煙暴表達(dá)的QUOTEWRSR礎(chǔ)丙定值的向下練累計頻率湊。編制QUOTEWRSR譯其頻數(shù)分布表板,列出各些組頻數(shù)曾、累計頻數(shù)碑QUOTEf臺掌；確定各組QUOTEWRSR要渠的秩次Ｒ三及平均秩次蔬專計算向下累緒計頻率常賢QUOTERn×100%閑喜，釀修正最后一迎個累計頻率博，最后一個塞累積頻率按味摸拌QUOTE1-14n×100%徐暗計算；喚將百分率培區(qū)換算為概率皇單位我午QUOTEprobit剝售,妙它邀涉故勇鍋枯撓喪不貢屋QUOTEprobit瓜唯為百分率柜壁對應(yīng)的標(biāo)準(zhǔn)斥正態(tài)離差u招加5。接表5QUOTEWRSR渾繳的分布渾累計真頻數(shù)打0.209丘1筒1打1打0.031享3酒3.137術(shù)3擺0.300福1懸2飼2射0.秧0625毀3.465國9旗0.342我1圖3斥3蜓0.093皇8無3.682腔0統(tǒng)0.391鹿1養(yǎng)4罵4脖0.125譽(yù)0攤3.849巷7吃0.416濱1鄙5拋5系0.156勤3籮3.99故00俘0.439澆1皇6蝴6賄0.187脾5南4.112僻9爭0.439寇1肯7海7拴0.218拋8近4.223申6皂0.465悲1鬼8悔8配0.250宣0芝4.325梁5南0.489嫂1圍9嶄9灰0.281慚3害4.420歷9蜓0.481持1套10懂10礙0.312層5踩4.511防2煮0.492世1妨11剛11遇0.343依8匙4.597脊7循0.498鐵1暴12巧12最0.375攪0撇4.681膚4棵0.500版1俊13贈13崗0.406鳴3氣4.762懸8斜0.534菌1摸14辜14遼0.437魂5廢4.842膝7逐0.534扁1鼠15是15柴0.468間8涌4.921驚6壇0.548宮1腰16貞16篇0.500充0萬5簽.0000方0.553掩1州17算17扭0.531帖3瘡5.078致4嶺0.568希1姑18掏18塘0.562騾5累5.157攝3競0.565盈1哄19嚴(yán)19摔0.593疊8池5.237愛2臭0.571強(qiáng)1琴20儲20耕0.625巖0續(xù)5.318蒜6鏈0.572隙1脖21眾21升0.656籮3刮5.402貝3聰0.583物1訓(xùn)22已22巧0.687企5具5.488銀8著0.583決1侮23艦23漲0.718缸8邁5.579犯1為0.596籍1陳24亭24葬0.750柿0蛛5.674傳5瀉0.590袍1置25匆25掉0.781過5恩5.776森4音0.594匯1應(yīng)26恨26乘0.812銷5架5.887飛1厭0.599瀉1權(quán)27充27杏0.843未8嗎6.01免00根0.59悼7孕1魄28康28璃0.875跡0技6.150朵3雞0.601賄1追29污29把0.906債3廢6.318口0糧0.607半1肺30的30莖0.937怠5辨6.534擋1梁0.619貧1勿31烤31庸0.968贈8享6.862冰7猛0.621妻1突32盞32碌0.992鼓2訪7.417寸6疾（4）計算濁回歸方程屠以累計頻率告對應(yīng)的概率四單位值益為自變量鎮(zhèn),哨值為變量徒,計算回歸刮方程掀。敞根據(jù)表5，川利用鈔軟件計算回腿歸方程結(jié)果盞為牛方程擬合優(yōu)賊度為0.9支147，獎?wù)f明通過檢嘴驗。隙（5）分檔謝排序然根據(jù)表5常辦用分檔數(shù)與犯百分位數(shù)突及概率單位斃Y,以匙值對評價對撥象進(jìn)行分檔齊排序株表6年常讓用分檔數(shù)與蜘百分位數(shù)績及概率單位必?fù)鯏?shù)撥百分位數(shù)保擋數(shù)與百分位數(shù)彎3彎15.86和6以下燈7椅10.02執(zhí)7過15.86企6簽33.36論84.13援4屑67.00帳3撤4百6.681蹈以下閣89.97蓄3兩6.681駕98.35埋2訓(xùn)50標(biāo)8捷1.222肝以下竊93.31持9太1.222務(wù)5含3.593論以下勞6.681渡3.593面22.66型3霧27.42欣5概50論72.57姜5分77.33素7蝴96.40脈7懷93.31舌9韻6錫2.275辟以下直98.67雹8橋2.275白9撿0.990穩(wěn)以下聚15.86油6圍0.99追50推4.746貞84.13來4菠15.86匹6箭97.72俊5誠37.07求7缺1.168咳以下身62.93壤1.168爪84.13魔4碼根據(jù)表5、叼6，結(jié)合本杰題需要，我腳們確定采用銅5檔的分法懶吃，結(jié)果如表蒼7所示：渴表7結(jié)果牌分析小等級黨等級對應(yīng)年斧份孕對應(yīng)年份總貸數(shù)飾的范圍客一級脊1985團(tuán)1勞<3.20墓二級璃1979、成1982、掀1984、心1986、禁1987砌5呢3.20~息4.40窩三級衛(wèi)1981、弦1983、醋1988、藥1990、箏1998農(nóng)5寬4.40~睡5.60隸四級訪1996、酸1980、公1989、核1990、松1991、圓1992、狂1993蚊1994、邁1995、司1997、臣2002、房2004、智2023瓣13渴5.60~蛛6.80奧五級乒1999、漿2000、離2001、嚇2003、限2005、例2006、盆2007、腔2023白8惡>6.80隨（6）秩和朋檢驗耀秩和檢驗用姓于檢驗分檔太是否合理，退同等檔級之房間會有不和小諧現(xiàn)象，用愿秩和檢驗可片以排除此類須問題的干擾輩。職燒秩和檢驗值疑，當(dāng)?shù)騋UOTEp>0.05蔑在，說明分檔運(yùn)和諧，否則倡不和諧。即我們利用賭程序翼（見程序4瓶）辜計算出各等粱級秩和檢驗棉值旅如下表俗表8各等醋級對應(yīng)秩和坐檢驗值販等級椅檢驗值仇是否和諧蜻一級疤1矮和諧抄二級徒0.982盟5施和諧外三級燥0.919仙9誰和諧男四級香0.780分9嗚和諧輛五級駕0.977叔9赤和諧膝通過表8我直們可以知道鏡各等級是和趣諧的，因此賽我們的風(fēng)險奔等級劃分合車?yán)?。?.3.幅問題三的模殺型建立與求早解糧柏模型原理緒在對興模型對實際梯問題進(jìn)行預(yù)案測時發(fā)現(xiàn)，悄影響預(yù)測精鼓度的因素有賺很多時，鐘其中發(fā)展灰駕數(shù)第增長率較小顛時，淋預(yù)測精度越慌高,肥當(dāng)原始數(shù)據(jù)表增長率變化濟(jì)較大時，即裳模型預(yù)測方售法中地發(fā)展紫灰數(shù)較大時蜜，預(yù)測精度剖低珍。嘴模型包預(yù)測，具體任步驟如下：株根據(jù)原始數(shù)央據(jù)，建立時匹間序列：浙通過將原始訪數(shù)據(jù)進(jìn)行一伯次性累加生鍵成新序列：其中：才建立廁模型相應(yīng)的巾微分方程為招：恥構(gòu)造累加數(shù)及據(jù)矩陣欠和常數(shù)向量譽(yù)用最小二乘何法求得灰色醒參數(shù)陽其中：逃稱為發(fā)展灰站數(shù)，QUOTEμ并宜稱為內(nèi)生控壟制灰數(shù)。當(dāng)設(shè)玩為待估參數(shù)投向量，解的虹：距求解微分方騎程，即可得取預(yù)測模型：嘩能模型的求解粱通過建立的涂模型，因式為降雨量存丈在著隨機(jī)性應(yīng)，所以我們規(guī)對目農(nóng)業(yè)用水，褲第三產(chǎn)業(yè)及歇生活其他用哪水，城市綠叼化率，常住患人口，拳進(jìn)行他預(yù)測畜擬合芝，利用駝編程（見附橡錄程序3）憤，畫出了擬萌合散點圖，焦并得到了擬窗合方程。它圓圈產(chǎn)點為原始值揪，菱形點為干預(yù)測值。勿圖2農(nóng)業(yè)趟用水的擬合宰圖脈從圖中可以催看到原始值設(shè)和預(yù)測值的形都是以遞減親的方向變動誰，原始值在即預(yù)測值的上泰下均勻分布拍。楚對蝕農(nóng)業(yè)用水懼的面擬合方程：披寧省圖3第三態(tài)產(chǎn)業(yè)及生活警其他用水傍的退擬合圖磁從圖中可以牌看到原始值演和預(yù)測值的痛都是以遞增盈的方向變動常，原始值在咱預(yù)測值的上廣下均勻分布末。泥對甲第三產(chǎn)業(yè)及搜生活其他用轎水決的GM插(1,1)強(qiáng)浸擬合方程燭圖4城市袍綠化率篇的房擬合圖坊從圖中可以陶看到原始值景和預(yù)測值的輩都是以遞增罷的方向變動釘，原始值在答預(yù)測值的上銹下均勻分布論。鞠對昨城市綠化率攏的GM奶(1,1)巖擬合方程：民圖5常住致人口由的驅(qū)擬合圖谷從圖中可以直看到原始值凝和預(yù)測值的倒都是以遞增摸的方向變動遵，原始值在貨預(yù)測值的上含下均勻分布菠。鄭對夾常住人口足的莖擬合方程：繼巨模型的檢驗約對睡模型進(jìn)行后艇驗差檢驗橡（1）某計算原始序被列標(biāo)準(zhǔn)差：副（2）光計算絕對誤貞差序列的標(biāo)水準(zhǔn)差：蕩（3）癢計算后驗差啟比值：厘（4）缺計算小誤差懸概率：糕等級高>0.95物<0.35敗好沃>0.80瓣<0.5共合格膛>0.70虛<0.65駕勉強(qiáng)合格槍<=0.7影0支>=0.6童5憑不合格廢表鄰9妖筒檢驗值對照嚷表閣用上述檢驗狂方法對主要迎風(fēng)險因子進(jìn)異行檢驗，結(jié)暴果如表塊10刪所示黎表史10降風(fēng)險因子砌檢驗形風(fēng)險因子翅農(nóng)業(yè)用水脹0.652禿8壯0.843吉8紡勉強(qiáng)合格亭第三產(chǎn)業(yè)及癢生活其他用偏水齊0.234瞇3處1同好貪城市綠化率粒0.178藥4網(wǎng)1吸好喘常住人口姨0.188泄1致0.968喬8拋好柔主要風(fēng)險因肝子基本上通村過了亞后驗差檢驗燙，所以可以伐用吃擬合方程來定預(yù)測未來兩汁年的主要風(fēng)加險的值，結(jié)浴果如表1疼1叮所示。勢表1員1齡悔2023年爹和2023送年風(fēng)險因子籌預(yù)測值梁風(fēng)險因子悄農(nóng)業(yè)用水懲第三產(chǎn)業(yè)及調(diào)其它用水而常住人口消污水處理率寒2023坊11.3弟21.09劇41宵1848墊103.5躺2023驅(qū)10.99調(diào)22.08片10蔽1893槳115.1涉再莖根據(jù)的MA秘TLAB程敏序絹（如附錄程怨序2）理計算出20負(fù)11年和2省012年的吩,烤查附件表2說百分?jǐn)?shù)與概這率單位對應(yīng)續(xù)值可以得到甩如下結(jié)論：被表1晉2忠2023段和2023勻風(fēng)險等級預(yù)筒測絮年份框風(fēng)險等級傅2023傍0.570區(qū)7餐5.178拔9抄第三等級燈2023塊0.606捧9浩5.271潮5紹第三等級秤六、模型的團(tuán)評價盆質(zhì)合比法因勺為數(shù)據(jù)太少嶼，主觀因素增太大，可能漲會造成誤差課。而且質(zhì)合狡檢驗仰缺點挽是嶄相同秩次較村多時，統(tǒng)計聲量要流經(jīng)過多次彎校正。薦關(guān)聯(lián)度分析當(dāng)?shù)娜毕菰谟谖浕舅枷肭凡腿睖?zhǔn)確和規(guī)撕范內(nèi)涵欠全搏面兩個方面耍，而本文也討無法避開這昂兩個方面，抽這一點有待酬改進(jìn)。嬸莊逐方法本身是概一種確定的閱外推，在處惜理數(shù)據(jù)、擬紛合曲線和得底到的擬合曲瞞線的過程中片都不考慮隨腸機(jī)誤差。參考文獻(xiàn)攜[1]蘇恨金明，王永央利，MAT烘LAB7.中0使用指南承[M]，北伐京：電子工填業(yè)出版社，菠2004.皆[2袍]菠江世宏，席MATLA劍B語言與數(shù)車學(xué)實驗[M縮]，北京：柄科學(xué)出版社范，2007背.輛[3]趙去靜,但琦,牛數(shù)學(xué)建模與歸數(shù)學(xué)實驗[危M],高等即教育出版社叛,2023循.1.惠[4]冉于啟康,張振援宇,張立柱曬,常用數(shù)學(xué)斬軟件教程,創(chuàng)人民郵電出榴版社,蘭2023.脈10電.摩[5]若朋毛吉，申應(yīng)用主成分堵和灰色關(guān)聯(lián)京度分析法評燃價共和盆地守天然草地產(chǎn)域草量的影響冶，躺:比//www與.doci在n貍/p-45我01029濫58.ht逮ml#do蛛cumen贈tinfo購，竿2023年指8月22日墾.玻[6]植亦吳清平宜,促張丹屋,藥秩和比法和薄幾種常用評坊價方法在醫(yī)早療質(zhì)量評價編中的比較敘[J],估中國醫(yī)院態(tài)統(tǒng)計訓(xùn),價2003眉.寇03車.工25撐.沙[7]吊馮長春望,纖公侯玉亭慧,幟城鎮(zhèn)土地評墾價中主成分蟻分析法的應(yīng)價用眼[J],拐質(zhì)中國國土資倚源經(jīng)濟(jì)伐,缸2007槳.軌07界.北15捕.增[8]屑榨何國華秩,險區(qū)域物流需舞求預(yù)測及灰質(zhì)色預(yù)測模型炸的應(yīng)用峰,粱身北京交通大意學(xué)學(xué)報組[J]癢,蛇2023以.奧01妄.儉05盟.豈[9]漿田志紅垮，擾我國水資源草保護(hù)存在的峰問題與對策白研究[J]屆，添科技情報開壁發(fā)與經(jīng)濟(jì)吐，州2023年罩02期脫.附錄押表1初1979到恐2023的梢相關(guān)數(shù)據(jù)害年份版污水處理率全（%）文農(nóng)業(yè)用水

適（億立方米亮）概工業(yè)用水

樓（億立方米眠）盤第三產(chǎn)業(yè)及獻(xiàn)生活他

降用水（億立踩方米）葬降雨量（毫勺米）謙1979禮10.2有24.18牲14.37磁4.37驅(qū)718.4門1980炊9.4偶31.83寒13.77鬼4.94酸3夢80.7皆1981件10.8嚷31.6常12.21狼4.3花393.2棕1982晶10.9竭28.81閥13.89沖4.52禁544.4岡1983繳10.2袋31.6墓11.24暈4.72踩489.9較1984僻10毒21.84埋14.37畫6儀4.017貴488.8癢1985耳10拳10.12琴17.2海4.39絹721糠1986易8.9捧19.46廟9.91晨7.18鎮(zhèn)665.3船1987蘆7.7便9.68塊14.01糕7.26碧683.9橋1988踢7.4胸21.99愁14.04勢6.4初673.3客1989幼6.6鋸24.42嫌13.77貓6.45嘩442.2芝1990速7.3障21.74卻12.34滿7.04主697.3劃1991互6.6辛22.7強(qiáng)11.9江7.43待747.9味1992哥1.2義19.94濾15.51遣10.98臉541.5處1993粘3.1柄20.35練15.28螺9.59頁506.7隙1994緞9.6喘20.93叔14.57墳10.37染813.2晴1995田19.4愛19.33階13.78妥11.77對572.5文1996順21.2若18.95夜11.76央9.3聰700.9翁1997望22揭18.12場11.1炮11.1盜430.9辛1998旅22.5連17.39遷10.84銳12.2燥731.7木1999桶25鴨18.45煮10.56欣12.7欠266.9靈2000殖39.4稀16.49淘10.52烏13.39匯371.1跪2001乒42押17.4煌9.2隔12.3徑338.9底2002常45浪15.5膽7.5幅11.6串370.4販2003假50.1押13.8打8.4甘13.6挽444.9初2004疾53.9卵13.5座7.7旬13.4捎483.5彎2005序62.4刑13.2序6.8句14.5評410.7列2006錦73.8探12.8襪6.2員15.3使318宜2007受76.2遇12.4幕5.8隱16.6穴483.9林2023亮78.9廉12齒5.2沾17.9浸626.3貓2023杏80.3售12更5.2巷18.3治480.6普2023來81繳11.4滲5.1爸18.8素522.5響表污2孤百分?jǐn)?shù)與概辜率單位對照都表伯%圣0.0錢0.1詠0.2季0.3愿0.4絡(luò)0.5爽0.6眨0.7撐0.8燒0.9掠0蔥-罪1.909春8忌2.121兆8桑2.252流2暖2.347偽9費2.424厘2都2.487槐9駛2.542滾7游2.591島1服2.634個4誰1蘇2.673于7部2.709黨6醒2.742促9嘴2.773吊8躍2.802裂7椅2.829羊9岡2.855株6尸2.879標(biāo)9宵2.903洞4抬2.925熔1可2菊2.946孝3泰2.966景5惑2.985淋9鑰3.004戰(zhàn)6槐3.022泛6壘3.040厘0教3.056糕9溉3.073墾2誼3.089醫(yī)0勝3.104抬3常3后3.119側(cè)2趁3.133療7惱3.147療8房3.161擠6多3.175斷9隊3.188厲1棋3.200桌9喝3.213恒4禍3.225帥6繼3.237池6賤4謙3.249鬼3把3.260您8非3.272恩1慶3.283障1臂3.294社0御3.304鈴6謙3.315物1跡3.325類3砌3.335喊4夫3.345足4術(shù)5掘3.355祖1促3.364兔8灘3.374貪2貝3.383凡6慕3.392優(yōu)8羨3.401貞8睡3.410折7綁3.419拍5難3.428突2惱3.426宅8來6關(guān)3.445墊2殲3.453砍6漆3.461芳8閑3.469夕9宏3.478這0浸3.485慢9踢3.493學(xué)7取3.501偉5孝3.509領(lǐng)1糧3.516挖7師7東3.524列2翼3.531殃6薄3.538豈9墨3.546約2廈3.553我4響3.560尺6包3.567曾5籮3.574待5遷3.581伸3蠢3.588膝2寄8息3.594蛾9旬3.601戲6模3.608喉3蘆3.614向8相3.621傅3旁3.627獸8純3.634冊2嚷3.640必5業(yè)3.646男8鋤3.653膝1拘9待3.659五2負(fù)3.665錦4搖3.671檔5炒3.677涉5刊3.683殘5塑3.689嘩4體3.695周3佛3.701荒2輪3.707焦0仁3.712戶7豬掛苦責(zé)愉廊悲夜膜丸舊三10駕3.718農(nóng)4動3.724摸1霸3.729押8辦3.735孝4裝3.740貿(mào)9變3.746非4淡3.751伯9風(fēng)3.754躍7毫3.762浪5夠3.768參1纏11殊3.773塔5負(fù)3.778索8林3.784頑0末3.789緞3啞3.794秀5藏3.799胡6瞇3.804嘩8泄3.809稻9撕3.815范0過3.820傍0耀12灘3.825前0歉3.830禾0他3.835斗0斷3.839支9異3.844蓋8軟3.849框7孤3.854儀5近3.859灣3寇3.864撒1和3.868讀9爆13腎3.873禿6泊3.878糧3仇3.883敬0畜3.887梁7廈3.892抵3黨3.896造9殲3.901誓5促3.906蘿1狠3.910辰7釋3.915呈2擁14畝3.919亭7準(zhǔn)3.924睛2庫3.926魚8至3.933揚(yáng)1諷3.937腫5頸3.941先9龜3.946干3諷3.950晝6飲3.955總0床3.959丹3冊定沿薪惕劉污當(dāng)搭略步居15弱3.963誤6軋3.967著8搜3.972專1綿3.976服3俊3.980射6轎3.984熄8愧3.989棟0遵3.993摧1妨3.997至3兩4.001蛇4許16報4.005客5零4.009饑6簽4.013濫7他4.017技8側(cè)4.021糕8醋4.025繡9找4.029餡9灑4.033撥9企4.037星9鼠4.041屠9泛17殊4.045浪8策4.049螞8誓4.053餃7搶4.057希6樹4.061寸5下4.065霸4才4.069桿3疾4.073挖1鞏4.077巡0冰4.080簽8誼18機(jī)4.084涼6剖4.088體4般4.092欲2幫4.096桐0列4.099襪8軟4.103飄5傻4.107譽(yù)3倍4.111腦0肺4.114訂7縮4.118醒4盟19聚4.122衡1增4.125驗8運(yùn)4.129俗5醬4.133錯1從4.136含7患4.140齡4之4.144值0抗4.147耐6雁4.151肌2堆4.154汽8棋構(gòu)抵乏浮坡孕傲女回族誓20晶4.158村4走4.161請9營4.165莫5寶4.169報0穿4.172常6儀4.176底1疼4.179叮6愚4.183汽1育4.186行6柴4.190辭1控21仔4.193伶6夏4.197離0證4.200腦5用4.203旬9照4.207未4驅(qū)4.210絲8界4.214囑2肯4.217循6資4.221吩0統(tǒng)4.224界4鑒22還4.227辯8杜4.231賽2己4.234熟5盯4.237什9值4.241燭2至4.244余6薯4.247卻9專4.251晨2增4.254息6文4.257刻9拖23祥4.261病2森4.264床4價4.267執(zhí)7施4.271耍0婚4.274模3瓦4.277泄5這4.280筋8爽4.284搏0扇4.287忠2在4.290床5填24估4.293學(xué)7會4.296魂9逃4.300麻1坡4.303矛3萄4.306拿5獨4.309帥7棗4.312晝9容4.316鹽0裙4.319甚2哲4.322鼓4偶癢鵲蓄訓(xùn)泳萌緣抱皂間扣25皆4.325愿5注4.328志7筆4.331限8流4.334周9盡4.338輪0撲4.341懷2帝4.344撥3滴4.347近4水4.350冶5番4.353墻6屈26存4.356怠7育4.359掩7隊4.362翼8箏4.365擁9土4.368返9泡4.372帖0夏4.375踏0式4.378田1豆4.381份1榮4.384濕2奸27差4.387劫2悲4.390拾8校4.393澆2慣4.396嚇2沒4.399栽2威4.402旺2巖4.405桶2宮4.408茅2春4.411個2撲4.414怒2慨28農(nóng)4.417號2悟4.420籮1饅4.423侍1街4.426閉0蛛4.429扛0乖4.431北9蘋4.434把9辰4.437氧8柱4.440鐵8慶4.443顯7示29犁4.446副6撞4.449川5芳4.452但4恐4.455情4蛛4.458網(wǎng)3腫4.461畝2當(dāng)4.464咐1亮4.467苗0任4.469存8篇4.472浮7敲挪鋼狀他盼牙鴿陳鄭贊號30期4.475躲6虹4.478念5患4.481遭3雀4.484遼2吳4.487駐1三4.489毅9澡4.498跌2湯4.495酬6您4.498畝5君4.501冒3儲31壓4.504走1芽4.505員0攜4.509餡8收4.512廚9角4.515云5臘4.518慘3修4.521導(dǎo)1婚4.523橋9曬4.526塌7月4.529刷5飄32著4.532炸3鋼4.535罵1尼4.537蝦9裙4.540錢7株4.543拼5驗4.546籮2陷4.549伴0夸4.551旁8池4.554吧6垃4.557攜3很33慰4.560注1沖4.562貢8咳4.565恩6姐4.568序4扣4.571屈1聽4.573勻9掏4.576碑6超4.579擋3崗4.582率1衣4.584渠5旺34職4.587壺5嚇4.590猾3脾4.593豪0徒4.595賽7潑4.598哭4斜4.601撈1扛4.603按9近4.606償6暗4.609紛3拌4.612拿0學(xué)無厘照姓底繁她垂狠辣民35豈4.614早7辛4.617焰4屯4.620狐1寬4.622懇8黑4.625畫5累4.628循1送4.630哈8珠4.633話5追4.636船2質(zhì)4.638姐9糾36滴4.641鼓5卻4.644姨2符4.646鞋9灣4.649染5秀4.652蓋2界4.654們9恭4.657相5沒4.660訴2盟4.662型8纖4.665應(yīng)5顏37造4.668霜1密4.670戲8粱4.673備4屆4.676衡1概4.678界7盈4.681只4快4.684警0狀4.686壞6政4.689計3夸4.691武9魯38誘4.694鴨5杠4.697貓1撐4.699五2狀4.702老4脆4.705弟0旨4.707吸6碌4.710回2腎4.712披9盈4.715嘗5怠4.718篩1澤39狗4.720你7弊4.723口3沈4.725厚9凡4.728兄5織4.731斧1去4.733艘7膚4.736師3削4.738刊9制4.741施5搬4.744瞎1曲論豆旗受褲隸辭鍬株稠俘40治4.746挨7遷4.749杰2縱4.751斜8悟4.754駱4返4.757牛0孫4.759慎6漢4.762并2響4.764乞7鑼4.767皂3安4.769鞋9鵲41役4.772雀5偉4.775節(jié)0濃4.777可6碑4.780候2叨4.782她7說4.785景3矩4.787砍9乏4.790豬4婆4.793疊0黑4.795技5詞42罪4.798集1長4.800沒7庫4.803慣2漏4.805滑8仍4.808俱3營4.810槳9跨4.813剖4谷4.816同0消4.818蠢5譜4.821摟1輩43次4.823溉6烘4.826穴2稱4.828偶7奏4.831庫3而4.833昌8肢4.836旋3稠4.838繪9棒4.841肢4坦4.844矛0棕4.846估5管44僅4.849晉0塌4.851伙6透4.854番1敬4.856桂6令4.859纏2禁4.861型7糕4.864劍2闖4.866襲8互4.869狐3怒4.871讀8律位纏煎雞申翼績貼錦茄彎45概4.874偉3達(dá)4.876亂9晴4.879譜4蠢4.881慈9浪4.884梳4長4.887按0差4.889尿5灘4.892斗0黨4.894接5沃4.897欲0烈46居4.899勞5頁4.902揀1神4.904余6息4.907粒1窮4.909番6續(xù)4.912廳2懸4.914蜻7彩4.917站2堅4.919絞7頁4.922柴2襪47嫂4.924幕7零4.927筋2椒4.929致8失4.932廳3催4.935龍8今4.937傳3籮4.939顆8熄4.942赴3娘4.944士8偉4.947困3佩48幟4.949暖8床4.952詳4焦4.954質(zhì)9莊4.957治4易4.959鞏9肉4.962機(jī)4和4.964僚9追4.967都4剝4.969徹9犁4.972綿4染49淺4.974梁9移4.977輸4敬4.979于9雁4.982婦5班4.985給0膨4.987輪5藏4.990快0京4.992峽5孝4.995唇0秩4.997膝5撇槽炮犬廉攀賀碰紹章丈罰50寒5.000赴0蝶5.002耳5統(tǒng)5.005費0升5.007的5服5.010拒0紫5.012駱5負(fù)5.015落0貪5.017脊5必5.020仗1繭5.022娛6喊51艦5.025饑1照5.027遇6御5.030遍1護(hù)5.032口6果5.035范1繪5.037魯6洲5.040槐1幕5.042巷6謝5.045迷1遍5.047太6算52帝5.050鬧2族5.052服7您5.055男2雜5.057暢7攔5.060葛2院5.062勸7事5.065掘2抬5.067抄7略5.070田2季5.072清8傅53止5.075罪3桿5.077商8脊5.080榜3美5.082池8佛5.085絕3缺5.087兼8襲5.090鑒4泊5.092蘭9轎5.095成4燭5.097饑9風(fēng)54尤5.100姻4千5.103綢0盡5.105爭5潔5.108更0匠5.110拾5是5.113灶0姿5.115記6澇5.118鉤1茂5.120逃6爺5.123由1兼壩允怕凳傅涼奶需縮皮剖55塔5.125售7丈5.128拉2殼5.130嘴7赴5.133斷2警5.135橡8過5.138幼3死5.140羅8齊5.143繩4屠5.145稱9身5.148腦4也56單5.151思0潛5.153枕5末5.156軟0娃5.158固6特5.161老1井5.163儉7驚5.166剖2導(dǎo)5.168礦7絹5.171咳3配5.173悶8資57海5.176搶4饅5.178冊9鎖5.181勻5磚5.184斗0偽5.186保6編5.189奉1旱5.191艦7霜5.194嶺2艦5.196內(nèi)8爸5.199伏3像58乳5.201善9議5.204倒5肢5.207售0避5.209話6且5.212臣1娛5.214轟7伶5.217依3呈5.219惜8譯5.222牛4陳5.225估0儀59掙5.227徐5肌5.230污1勉5.232女7違5.235磨3文5.237揮8壽5.240揚(yáng)4脅5.243寄0蟻5.245肝6若5.248番2高5.250詢8蟻牧兄熟爭吵辯透渠魔招濫60遇5.253帥3漆5.255約9熔5.258駛5跳5.261常1攪5.262垮7廳5.266排3皆5.268抖9睡5.271息5綢5.274蟻1糊5.276抄7惰61租5.279調(diào)3詠5.281秀9壤5.284投5炮5.287堤1耍5.289凱8彩5.292甚4朽5.295把0關(guān)5.297影6兇5.300英2愛5.302籌9客62瞇5.305導(dǎo)5丑5.308掏1啄5.310淹7必5.313艦4糕5.316奪0忠5.318椒6瀉5.321雙3增5.323關(guān)9擾5.326熱6瘦5.329賽2壟63純5.331東9尾5.334極5喪5.337嚼2據(jù)5.339課8符5.342同5或5.345苦1哄5.347啟8鹽5.350嫁5畏5.353忌1早5.355沉8肅64炭5.358芳5制5.361扇1悠5.363仰8豆5.366搜5渴5.369謎2菠5.371騾9扛5.374陽5茶5.377府2淡5.379蛋9踏5.382垂6該辮皂榮叉趨歉凍質(zhì)卵盜酒65鍛5.385裂3腔5.388隙0瘦5.390匪7倍5.393渴4茶5.396煎1錦5.398忽9料5.401肺6賢5.404冊3式5.407獨0猜5.409語7裳66干5.412湯5投5.415來2轟5.417壞9棋5.420店7鵲5.423科4外5.426罵1譽(yù)5.428牢9猶5.431狂0稠5.434埋4誓5.437脾2淡67潛5.439以9斧5.442慮7例5.445功4檔5.448誼2體5.451堪0巾5.453猜8蚊5.456貓5鉆5.459蜘3昆5.462啊1項5.464根9愉68烈5.467卷7棍5.470璃5清5.473乎3甲5.476猾1楊5.468滴9牌5.481停7芒5.484多5站5.487切4當(dāng)5.490惰2鴿5.493嘴0陰69眠5.485夠8們5.498寇7改5.501絡(luò)5諸5.504斬4總5.507遵2哄5.510麗1悄5.512嫁9探5.515憑8餐5.518慮7屢5.521胳5追蠻扎千字嬌星腥熟搭技租70離5.524犁4攀5.527陶3菠5.530乓2撓5.533澆0宗5.535胃9懸5.538燕8廁5.541哪7戰(zhàn)5.544坐5題5.547悟6傘5.550步5奔71檢5.553嗓4逝5.556溜3溪5.559渠2裝5.562小2戚5.565相1規(guī)5.568羽1隆5.571徐0擦5.574希0炮5.576所9赴5.579塞9守72四5.582乖8擔(dān)5.585搞8粥5.588閣8小5.591溝8及5.594杏8仆5.597叛8灣5.600妥8煮5.603縮8怖5.606話8叛5.609輝8秘73第5.612叉8服5.615提8格5.618蓬9咳5.621三9承5.625狀0象5.628柜0獅5.631竹1父5.634鉛1抗5.637禿2悠5.640斗3御74探5.643聲3飲5.646難4蠟5.649國5找5.652擁6臺5.655稼7模5.658遼8到5.662哲0生5.565和1請5.668膀2泛5.671伯3施影樣字驢加垃柜領(lǐng)仰農(nóng)族75晴5.674鍋5友5.677最6加5.680他8容5.684待0腸5.687桑1驅(qū)5.690責(zé)3筐5.693塊5傷5.696荷7蝶5.699還9盈5.703幸1晴76譯5.706練3主5.709妥5篇5.712骨8監(jiān)5.716助0鈔5.719蹈2賺5.722跑5議5.725側(cè)7姨5.729辜0顫5.732臥3延5.735牌6緊77吧5.738特8床5.742離1俘5.745熟4偽5.748濁8超5.752匆1衰5.755滋4貨5.758凳8類5.762段1賊5.765貌5伶5.768回8末78棵5.772黃2乳5.775哀6匙5.779緩0賄5.782蛙4篇5.785煤8劑5.789勉2御5.792導(dǎo)6翁5.796兇1疤5.799番5唯5.803磚0紹79悠5.806遙4沈5.809箭9毫5.813練4區(qū)5.816灶9脅5.820澤4漠5.823袖9愿5.827洗4圾5.831鴉0臺5.834新5減5.833喇1勇奸鏈招勵梁嬌甚煌討叫據(jù)80胳5.841晝6擊5.845蛙2甩5.848姥8燕5.852視4恥5.856彎0損5.859賤6保5.863裙3滔5.866窄9律5.870恭5炎5.874僚2襪81路5.877衡9境5.881蘿6丹5.885高3莊5.889掉0過5.892渡7陸5.896妖5斬5.900膨2寺5.904慎0振5.907卷8救5.911遺6猴82磨5.915折4杜5.919銀2蓋5.923咳0浮5.926叫9飼5.930輸7域5.934劉6埋5.938夫5議5.942宿4酷5.946網(wǎng)3鹿5.950反2掌83挪5.954理2肚5.958衣1帆5.962賴1落5.966隆1端5.970冰1啦5.974朽1緒5.978富2衡5.982煌2肉5.986域3慮5.990紛4施84盜5.994膛5健5.998芒5聯(lián)5.002踢7丸6.006歲9浙6.011聽0箱6.015謊2津6.019屯4測6.023剃7倒6.027稿9招6.032翅2暗萍柳風(fēng)器充稀攜再削塔癢85免6.036夾4圈6.040梨7烘6.045六0侄6.049縫4超6.053鐵7悶6.058管1漸6.062侄5志6.066短9挑6.071垂4測6.075您8偏86兼6.080肚3狀6.084猶8益6.089北3德6.092斯9錫6.098補(bǔ)5愁6.103葡1臘6.107井7臘6.112頸3雖6.117焦0患6.121紐7貢87寫6.126償4多6.131堪1克6.135預(yù)9逗6.140使7培6.145鋼5杏6.150惕3碌6.155么2赴6.160徐1贊6.165經(jīng)0百6.170干0泡88鼠6.175紐0放6.180吸0裳6.185虜0刻6.190氏1德6.195退2禽6.200淚4眉6.205須5逃6.210宣7刻6.216蘭0單6.221房2貨89迅6.226銹5澇6.231犬9暴6.237川2金6.242奪6盾6.243受1熱6.253紋6法6.259鋤1值6.264均6抱6.270期2拘6.275紗9宇漫偏碼瀉食規(guī)史端苗熱株90瓦6.281餡6菊6.267敘3婆6.293寺0熟6.298碼8飛6.304摘7咸6.310似6噸6.316甩5蟲6.322票5躲6.328者5軋6.334坡6訓(xùn)91湯6.340歸8哄6.346零9呈6.353授2疾6.359泉5擠6.365文8籮6.372烏2屢6.378臉7濫6.385攪2貸6.391迷7糕6.398區(qū)4賣92傘6.505泥1作6.411裁8仿6.418吵7創(chuàng)6.425短5兇6.432騎5皇6.439丈5涉6.446畏6準(zhǔn)6.453強(qiáng)8幫6.461均1窗6.458敘4坡93捏6.575蓄8活6.483皮3鄭6.490垂9圍6.498賺5狗6.506暑3令6.514即1誰6.522款0句6.530衛(wèi)1春6.532拒8萍6.548歸4鴉94育6.554完8孤6.563分2諷6.571吧8伴6.580魄5允6.589雞3邪6.598施2屋6.607弓2悄6.616黨4錯6.625彩8稠6.635繩2舞顏浴懸沾崖戴闊舉小三妹95幕6.644吸9雜6.654載6議6.664扒6進(jìn)6.674臺7激6.684壤9茂6.695鈔4暴6.705散0幼6.716悼9型6.727衫9嬌6.739骨2表96憲6.750銳7漁6.762拘4刃6.774事4進(jìn)6.786覺6盤6.799鳥1允6.811局9臂6.825甩0辯6.838遷4礎(chǔ)6.852降2斑6.866值3姥97熊6.880像8抵6.895雖7茶6.911貌0碑6.926稿8蘭6.943泥1指6.960瘋0畫6.977尼4呀6.995攏4彩7.014降1撿7.033仙5露98何7.053煎7肚7.074汗9德7.096掩9死7.120顧1疏7.144鞭4雅7.170霧1磁7.197繩3謝7.226曬2澤7.257錄1恭7.290蠶4刺99俯7.326職3笑7.365棋6遞7.408戴9廊7.457惰3酬7.512露1宇7.575前8覆7.652姥1啟7.747掀8郵7.878六2知8.090鋼2程序1沃aw=[2紹4.18備4.3型7肌897.1鞋10賀.2喉718.戶4滋31.83見4.環(huán)94警904.柔39陰.438正0.7厘31.6烏4.拾3919饞.2站10.8弟39抹3.2賣28.81啟4.私52斬935掉10.9征54慰4.4限31.6賣4.機(jī)72廢950室10.2聞48腹9.9鍛21.84也4.拒017涂965別104掌88.8體10.12住4.輔39套981珍107版21與19.46手7.呼18椒1028竿8每.966哈5.3劃9.68匠7.付26鉤1047鹽7核.768喪3.9州21.99頓6.衛(wèi)4106通1多7.46騾73.3限24.42秧6.穗45止1075閥6瀉.644械2.2籍21.74另7.燃04棕1086洗7仗.369耳7.3圾22.7咐7.充43著1094廟6叫.674程7.9跟19.94舊10艷.98綢1102殼1富.254寸1.5艷20.35業(yè)9.它59換1112緒3階.150糕6.7蜂20.93造10淡.37噴1125歪9痕.681石3.2藍(lán)19.33擇11列.77寸1251盆.11棚9.4文572聽.5振18.95蘿9.勸3125束9.4灣21.2歉70敬0.9腹18.12嚇11既.1濟(jì)1240屆2肝243尾0.9瘦17.39些12族.2駝1245川.62跨2.5運(yùn)731擱.7型18.45廈12窯.7曠1257惜.22經(jīng)526凳6.9粗16.49腔13話.39老1363墾.63寨9.4給371禿.1使17.4論12駐.3待1385愚.14己233輸8.9貍15.5逝11應(yīng).6昂1423癢.24勸537壇0.4骨13.8鬧13鐮.6屬1456食.45唇0.1常444顫.9于13.5較13禁.4門1492想.75棒3.9熔483傅.5李13.2繞14鐮.5萬1538點6增2.4盼410無.7夜12.8狡15探.3志1581鹽7剪3.8遭318都12.4筆16蝴.6辰1633蠶7艦6.2楚483膽.9起121銜7.9歉169奸5受78.9或62道6.3蘋121簡8.3崖175允5筍80.3獅48亦0.6半11.4逗18暢.8劣1961士.98航152它2.5莊];錄%導(dǎo)入數(shù)據(jù)三a=aw(條1:end室,:);附%提取指標(biāo)跑數(shù)據(jù)糾a(:,[煤4,5])淚=-a(:幣,[4,5易]);補(bǔ)%低優(yōu)度處脅理以ra=ti耕edran頃k(a);獨%對每個指徒標(biāo)值分別編掌秩，即對a健的每一列分權(quán)別編秩沖[n,m]勉=size蜓(a);拿%計算矩陣鉗的維數(shù)仰b=[0.注2573,削0.253廳7,0.2晌533,0叢.1215泰,0.09魔01];敢%權(quán)重趕w=b/s某um(b)享;廉%歸一化處橋理體WRSR=壩ra*w'殖/n;瘡[WRSR形i]=s下ort(W贊RSR,'扯ascen穗d')陸%升序排列刃WRSRn=32;呢p=[1:惱n]/n;旱%計算累計強(qiáng)頻率所p(end政)=1-1伴/(4*n貓);魚%修正最后縣一個累計頻匆率，最后一比個累積頻率待按1-1/危(4n)估載計距p=[3.降1337牛3.461顏8嚴(yán)3.677源5惹3.849書7屬3.989揉0鐵4.114輪7帝4.224順4落4.325篇5暢4.420顫1禍4.512哲9敲4.598棕4套4.681倉4魚4.762潑2虎4.844羅0盲4.922幸2瓣5.000峽0蜻5.077恐8劍5.158給6熟5.237個8快5.318純6砌5.401椒6哀5.490替2耗5.579坑9狹5.674小5砌5.775而6侄5.889夜0培6.011職0重6.150較3側(cè)6.316害5脈6.506航3唯6.866礙3佩7.408巷9];歷X=[on插es(n,棟1),p]稻;逝%構(gòu)造一元牙線性回歸分尚析的數(shù)據(jù)矩咐陣填[ab,a乳bint,隙r,rin佳t,sta灑ts]=r全egres頂s(WRS薯R,X)灶%一元線性月回歸分析粥WRSRf撲it=ab披(1)+a腥b(2)*努probi出t才%計算WR想SR的估計抓值程序2：脅aw=[1亡1.3石21屈.0941跑184畫810衰3.5菠1桑0.99既22組.081洞1893碗11傳5.1朽1披0.397角024移.2炕1987垃142.制4吵1使0.11哥25秒.33亡2036揚(yáng)158.董4]才；%導(dǎo)入數(shù)據(jù)犧a=aw(痰1:end郵,:)畢；%提取指分標(biāo)數(shù)據(jù)葛a(:,[紗4,4])煤=-a(:蓄,[4,4員])掘；%低優(yōu)度匠處理席ra=ti眾edran覺k(a)陸%對每個指飄標(biāo)值分別編環(huán)秩，即對a俊的每一列分憑別編秩率[n,m]斗=size色(a)粱%計算矩陣惑的維數(shù)子b=[0.孤2573,獲0.253沸7,0.2陷533,0駱.1215狼]做%權(quán)重關(guān)w=b/s飾um(b)循%歸一化處縫理柔WRSR=撲ra*w'潤/n砌[WRSR儉i]=s宿ort(W既RSR,'僻ascen牌d')竊%升序排列東WRSR繼程序3：唱灰色預(yù)測的旨程序襯x=[飄4.6924.5424.1110.6212.860.923-6.299.52-7.713.2523.095.26-0.2623.9925.550.4514.54-5.8618.072.7327.4923.5419.718.517.413.211.312.2110.913.712.2]'艘sizex當(dāng)d2=歪size(章x,2);舌%求數(shù)組長農(nóng)度k=0;追fory觸1=x撿k耳=k+1;抹i甚fk>1超顆x1排(k)=x殃1(k-1續(xù))+x(k鋸);敵顯兼%累加生成在德z1宅(k-1)們=-0.5典*(x1(忠k)+x1唯(k-1)蹤);千脆輔%z1維數(shù)俘減1，用于處計算B鞭姐yn預(yù)1(k-1擠)=x(k冷);各e皇lse叢俘x1除(k)=x憶(k);甲e儉ndend牌%x1,z汽1,k,y硬n1蝦sizez書1=siz沫e(z1,宮2);驚%size葡(yn1)銷;賊z2=棍z1';犬z季3=o鼻nes(1膝,size立z1)';照YN=職yn1';梁滅%轉(zhuǎn)置%YN豐B=[z2罵z3];團(tuán)au0=i聯(lián)nv(B'秋*B)*B上'*YN;賀au=層au0';扶%B,au淺0,au久afor藍(lán)=au(徹1);責(zé)ufor繡=au(鄭2);歲ua=疼au(2)甩./au(賀1);梢%afor泉,ufor戲,ua濟(jì)%輸出預(yù)測遙的a舒u和u腿