水資源短缺風(fēng)險綜合評價模型-數(shù)學(xué)建模競賽參賽論文

PAGEPAGE30封一答卷編號（參賽學(xué)校填寫）：THSY06答卷編號（競賽組委會填寫）：論文題目：B水資源短缺風(fēng)險綜合評價模型組別：本科生參賽隊員信息(必填)：姓名專業(yè)班級及學(xué)號聯(lián)系電話參賽隊員1數(shù)學(xué)與應(yīng)用數(shù)學(xué)200906010242參賽隊員2數(shù)學(xué)與應(yīng)用數(shù)學(xué)200906010220參賽隊員3數(shù)學(xué)與應(yīng)用數(shù)學(xué)200906010145 參賽學(xué)校：報名序號：（可以不填）封二答卷編號（參賽學(xué)校填寫）：THSY06答卷編號（競賽組委會填寫）：評閱情況（學(xué)校評閱專家填寫）：學(xué)校評閱1.學(xué)校評閱2.學(xué)校評閱3.評閱情況（省賽評閱專家填寫）：省賽評閱1.省賽評閱2.省賽評閱水資源短缺風(fēng)險綜合評價模型摘要本文通過對水資源系統(tǒng)影響因素的分析，確立了水資源短缺的主要風(fēng)險因子：風(fēng)險率、脆弱性、重現(xiàn)期、可恢復(fù)性、風(fēng)險度，并以這些因子為衡量指標(biāo)建立了數(shù)學(xué)模型對北京市水資源短缺風(fēng)險進(jìn)行了綜合評價。在問題1中，建立了衡量指標(biāo)模型、概率模型，運用方差、標(biāo)準(zhǔn)差對風(fēng)險因子進(jìn)行了檢驗。在問題2中，綜合問題1的風(fēng)險因子建立了一個描述北京市水資源短缺風(fēng)險程度的模糊綜合評價模型，在模型建立過程中運用了分類、劃分方法將水資源短缺風(fēng)險因子分為5個等級，并用層次分析法確定了風(fēng)險因子的權(quán)系數(shù)，構(gòu)造隸屬函數(shù)確立了水資源短缺風(fēng)險綜合評價模型，最終劃分了水資源短缺風(fēng)險的5個等級。在問題3中，運用圖形分析法，根據(jù)北京市水資源情況圖像及問題2中建立的水資源短缺風(fēng)險綜合評價模型對北京市未來兩年的水資源短缺風(fēng)險進(jìn)行了預(yù)測，預(yù)測結(jié)果表明：盡管加大再生水利用量、南水北調(diào)工程在一定程度上緩解了北京市水資源短缺的緊張局面，但未來兩年北京市水資源短缺風(fēng)險仍將處于高風(fēng)險水平。問題4則對北京市水行政主管部門提出了解決水資源短缺的可行性方案。在整個模型建立過程中還用到了EXCEL、MATLAB等工具。 關(guān)鍵詞北京、水資源短缺風(fēng)險、模糊數(shù)學(xué)、概率、風(fēng)險因子問題重述改革開放以來，我國工業(yè)和城鎮(zhèn)生活用水持續(xù)增長，由于氣候條件、水利工程設(shè)施、工業(yè)污染、人口規(guī)模等因素，近年來我國北方地區(qū)水資源短缺問題日益嚴(yán)重，而北京人均水資源占有量不足300立方米，為全國人均的1/8，世界人均的1/30，已成為世界上水資源嚴(yán)重缺乏的大都市之一?，F(xiàn)根據(jù)《北京2009統(tǒng)計年鑒》及市政統(tǒng)計資料提供的有關(guān)1979年至2000年北京市水資源短缺的狀況表（見附錄）討論如下問題：1、北京市水資源短缺風(fēng)險的主要風(fēng)險因子。 （影響水資源的因素，如：氣候條件、水利工程設(shè)施、工業(yè)污染、農(nóng)業(yè)用水、管理制度，人口規(guī)模等。）2、對北京市水資源短缺風(fēng)險進(jìn)行綜合評價并根據(jù)綜合評價標(biāo)準(zhǔn)作出風(fēng)險等級劃分，并考慮如何對主要風(fēng)險因子進(jìn)行調(diào)控才能使得風(fēng)險降低？3、對北京市未來兩年水資源的短缺風(fēng)險進(jìn)行預(yù)測，并提出應(yīng)對措施。4、以北京市水行政主管部門為報告對象，寫一份建議報告。問題分析近年來，受氣候變化和經(jīng)濟(jì)社會不斷發(fā)展的影響，水資源短缺問題日趨嚴(yán)重，對水資源短缺風(fēng)險的研究已引起了廣泛的重視。要對北京市水資源短缺風(fēng)險進(jìn)行綜合評價，首先要明確北京市水資源短缺風(fēng)險的影響因子及其影響程度的大??；其次，要建立模型考察各風(fēng)險因子對水資源短缺風(fēng)險的具體影響。在影響北京市水資源短缺風(fēng)險的各個主要因子中又有其它因素對其影響，在這里我們忽略了影響較小的因素。判斷北京市水資源短缺風(fēng)險問題要綜合考慮，既要考慮各因素對風(fēng)險因子的影響，又要考慮各風(fēng)險因子對水資源短缺風(fēng)險的影響，然后再對其進(jìn)行綜合評價。根據(jù)此題特征，我們很容易就想到了模糊概率的相關(guān)知識，通過計算概率來實現(xiàn)風(fēng)險的綜合評價。為此，先構(gòu)造各風(fēng)險因子的評價指標(biāo)，然后，再根據(jù)各風(fēng)險因子來構(gòu)造水資源短缺風(fēng)險的評價函數(shù)。最后，根據(jù)模型函數(shù)進(jìn)行風(fēng)險等級劃分并提出相應(yīng)的解決方案。符號說明符號說明S水資源失事狀態(tài)Z水資源正常狀態(tài)M總用水量G水資源總量Xt水資源系統(tǒng)狀態(tài)變量It水資源系統(tǒng)狀態(tài)量ui（i=1、2、3、4、5）分別表示衡量風(fēng)險率、脆弱性、重現(xiàn)期、可恢復(fù)性、風(fēng)險度的指標(biāo)T水資源系統(tǒng)工作總歷時T1水資源系統(tǒng)失事總時間QZt第t年水資源短缺量SZt第t年水資源總量X失事狀態(tài)下缺水量年均值Y水資源總量年均值Zt水資源可恢復(fù)狀態(tài)量THF水資源可恢復(fù)狀態(tài)總量TSX水資源狀態(tài)總量D(x)衡量風(fēng)險度的方差E(x)衡量風(fēng)險度的標(biāo)準(zhǔn)差d(q,n)第n與n+1次失事時間間隔Vi(i=1,2,3,4,5)水資源短缺風(fēng)險等級U、V定義的兩個有限域i（i=1,2,3,4,5）各風(fēng)險的權(quán)系數(shù)bi（i=1,2,3,4,5）水資源短缺風(fēng)險等級判別指標(biāo)A、B分別為U、V的模糊子集RU變換關(guān)系矩陣Uvi()Ui對Vi的隸屬函數(shù)模型假設(shè)北京市人均年用水量一定；北京市單位產(chǎn)業(yè)用水量一定；不考慮水資源自身凈化能力；假設(shè)水資源短缺風(fēng)險的各風(fēng)險因子間影響程度可以忽略；在計算過程中忽略微小數(shù)值對結(jié)果的影響；6、忽略信息的有限性和不完全性對模型的影響。模型建立、求解與分析水資源系統(tǒng)是一個復(fù)雜的大系統(tǒng)，廣泛存在著隨機(jī)性和模糊性，由于隨機(jī)性是因果率的破缺、模糊性是排中率的破缺，所以在水資源短缺風(fēng)險評價模型的設(shè)計中應(yīng)同時考慮這兩種因素的影響。為了緩解北京市水資源短缺問題，建立數(shù)學(xué)模型對北京市水資源短缺風(fēng)險進(jìn)行綜合評價。根據(jù)北京市水資源短缺風(fēng)險的主要風(fēng)險因子：風(fēng)險率、脆弱性、重現(xiàn)期、可恢復(fù)性、風(fēng)險度建立如下模型：問題1水資源是人類賴以生存的能源。影響水資源的因素有很多，例如：氣候條件、自然災(zāi)害、水利工程、工農(nóng)業(yè)污染、污水利用率、管理制度、人口數(shù)量及規(guī)模、入境水量等，而這些因素都影響著水資源系統(tǒng)處于失事狀態(tài)時所經(jīng)歷的時間、損失的嚴(yán)重程度、失事的次數(shù)、恢復(fù)到正常狀態(tài)經(jīng)歷的時間等。綜合以上原因可以將水資源的風(fēng)險因子概括為：1、風(fēng)險率考慮到水資源系統(tǒng)的不確定性，如果供水系統(tǒng)處于失事狀態(tài)，且水資源系統(tǒng)的工作有長期的記錄，風(fēng)險率可以定義為水資源系統(tǒng)不能正常工作的時間與整個工作歷時之比，即水資源系統(tǒng)風(fēng)險率衡量指標(biāo)為：其中 2、脆弱性脆弱性是指水資源系統(tǒng)處于失事狀態(tài)的平均損失嚴(yán)重程度。為定量描述水資源系統(tǒng)的脆弱性，我們假設(shè)水資源系統(tǒng)第t年水資源短缺量為QZt水資源總量為SZt，T1為水資源系統(tǒng)失事總時間,則水資源系統(tǒng)脆弱性衡量指標(biāo)為：上式中,失事狀態(tài)下缺水量年均值X，水資源總量年均值Y為：3、重現(xiàn)期重現(xiàn)期是指再次進(jìn)入失事狀態(tài)S時所經(jīng)歷的時間。則水資源短缺重現(xiàn)期衡量指標(biāo)為：其中，重現(xiàn)期4、可恢復(fù)性可恢復(fù)性即水資源系統(tǒng)從失事狀態(tài)到正常狀態(tài)的可能性。系統(tǒng)的可恢復(fù)性越高，表明該系統(tǒng)從失事狀態(tài)恢復(fù)到正常狀態(tài)的能力越強(qiáng)，反之，則越弱。根據(jù)已知條件它可以用如下指標(biāo)來衡量：其中， 5、風(fēng)險度標(biāo)準(zhǔn)差反映了變量的穩(wěn)定與波動、集中與離散的程度，標(biāo)準(zhǔn)差越小，波動越小，反之，則波動越大。利用標(biāo)準(zhǔn)差的特征可以定義以下指標(biāo)來衡量風(fēng)險度的大?。浩渲?，方差在類似分析中，一般假定P1=P2=……=PT=1/(T-1)標(biāo)準(zhǔn)差，問題2針對問題1中的各種風(fēng)險因子，建立一個基于模糊綜合評價的水資源短缺風(fēng)險模型。如果北京市年水資源總量小于年需水量，即水資源處于短缺狀態(tài)。則基于水資源系統(tǒng)的模糊不確定性，可以構(gòu)造一個合適的模糊綜合評價模型：設(shè)U={U1,U2,U3,U4,U5}和V={V1,V2,V3,V4,V5}為兩個有限域，A={1，2，3，4，5},0≤i≤1,B={b1,b2,b3,b4,b5}分別為U、V的模糊子集,則模糊綜合評價可表示為下列模糊變換：B=ARU 由上述假設(shè)可知，關(guān)系矩陣RU與水資源短缺的主要風(fēng)險因子有關(guān)。關(guān)系矩陣RU中的元素即為因素Ui對應(yīng)等級Vi的隸屬度，其值可根據(jù)各級評價因子的實際數(shù)值對照各因子的分級指標(biāo)推求。下面我們將評語分為5個等級：表一：各風(fēng)險因子指標(biāo)分級水資源短缺風(fēng)險等級風(fēng)險率指標(biāo)U1脆弱性指標(biāo)U2可恢復(fù)性指標(biāo)U3重現(xiàn)期指標(biāo)U4風(fēng)險度指標(biāo)U5V1低≤0.20≤0.20≥0.80≥9.00≤0.20V2較低0.20-0.400.20-0.400.60-0.806.00-9.000.20-0.60V3中0.40-0.600.40-0.600.40-0.603.00-6.000.60-1.00V4較高0.60-0.800.60-0.800.20-0.401.00-3.001.00-2.00V5高≥0.80≥0.80≤0.20≤1.00≥2.00由于水資源風(fēng)險率、脆弱性、風(fēng)險度是越小越優(yōu)性指標(biāo)，而可恢復(fù)性、重現(xiàn)期是越大越優(yōu)性指標(biāo)，所以，對于U1、U2、U3、U4、U5各級評語構(gòu)造如下隸屬函數(shù)： 從而，得到RU的具體向量表示：利用層次分析法確定水資源短缺風(fēng)險各因子的權(quán)重系數(shù)，即A={1，2，3，4，5}={1/2，1/4,1/8,1/12,1/24}于是可得到如下綜合評判向量：在綜合評判中，我們可以將上述的向量轉(zhuǎn)化為一般實數(shù)的加法，即最后，選取max{bj}對應(yīng)的評語作為水資源短缺風(fēng)險的綜合評價結(jié)果。根據(jù)上述評價標(biāo)準(zhǔn)，我們可以將水資源短缺風(fēng)險分為如下級別： 表二：水資源短缺風(fēng)險級別評價 水資源短缺風(fēng)險等級風(fēng)險級別水資源短缺風(fēng)險的特征V1(1<b)低風(fēng)險可以忽略的風(fēng)險V2(0.9<b<1)較低風(fēng)險可以接受的風(fēng)險V3(0.8<b<0.9)中風(fēng)險邊緣風(fēng)險V4(0.7<b<0.8)較高風(fēng)險不可接受風(fēng)險V5(b<0.7)高風(fēng)險災(zāi)變風(fēng)險經(jīng)過Matlab程序運行得到的水資源短缺風(fēng)險等級指標(biāo)顯示，近年來北京市水資源短缺風(fēng)險一直處于高風(fēng)險水平。由問題1建立的模型及本問題中建立的模糊綜合評價模型可知，當(dāng)風(fēng)險率指標(biāo)、脆弱性指標(biāo)、風(fēng)險度指標(biāo)越低時，水資源短缺風(fēng)險程度越低；當(dāng)重現(xiàn)期指標(biāo)、可恢復(fù)性指標(biāo)越高時，水資源短缺風(fēng)險程度越低；因此，要想降低水資源短缺風(fēng)險程度就要調(diào)控主要分險因子。我們可以通過高新技術(shù)、水利工程、人口控制等方法降低水資源系統(tǒng)的風(fēng)險率、脆弱性、風(fēng)險度或提高水資源系統(tǒng)的可恢復(fù)性與重現(xiàn)期，從而使水資源短缺風(fēng)險程度降低。問題3

根據(jù)北京市歷年主要用水分布圖，我們知道北京市近年來農(nóng)業(yè)用水在逐年減少，工業(yè)用水也呈現(xiàn)出負(fù)增長趨勢，但隨著社會的發(fā)展第三產(chǎn)業(yè)等用水卻逐年顯著增加。且根據(jù)北京市歷年用水量與水資源總量比較分布圖及北京市政統(tǒng)計資料提供的數(shù)據(jù)顯示，1999年—2008年由于干旱北京市一直處于水資源短缺狀態(tài)，而2008年基本處于平衡是南水北調(diào)工程開始投入使用從而提供了大量水資源。在研究北京市水資源短缺現(xiàn)狀的基礎(chǔ)上，采用概率模型對北京市水資源短缺風(fēng)險進(jìn)行預(yù)測：在對水資源短缺風(fēng)險因子進(jìn)行分析后，將他們的變化作為參數(shù)代入水資源短缺風(fēng)險綜合評價模型中，通過計算機(jī)數(shù)據(jù)分析，得到未來北京市水資源短缺的風(fēng)險數(shù)據(jù)，對應(yīng)數(shù)據(jù)就可以得到其水資源短缺風(fēng)險程度。結(jié)果顯示，在保持現(xiàn)有政策和技術(shù)水平下，加大再生水回用，南水北調(diào)工程雖然在一定程度上可以緩解北京市水資源短缺狀況。但研究表明，水資源短缺風(fēng)險因子對北京市水資源供需平衡有很大影響，如果沒有南水北調(diào)工程，未來兩年內(nèi)整個首都圈的水資源短缺風(fēng)險將會處于高風(fēng)險水平，水資源供需狀況極度危險，對水資源采取有效的風(fēng)險管理措施已刻不容緩。應(yīng)對措施：1、控制外來人口入京數(shù)量，加大計劃生育的實施力度，綜合調(diào)控北京市人口數(shù)量，抑制水資源需求過度膨脹。2、實施法律制度式節(jié)水：硬性規(guī)定用水單位的用水量，對違規(guī)單位進(jìn)行適當(dāng)?shù)膽土P措施以鼓勵市民節(jié)約用水。3、解決全國城市均衡發(fā)展問題，調(diào)整各城市間的空間結(jié)構(gòu)和產(chǎn)業(yè)布局，達(dá)到資源共享。4、提高科學(xué)技術(shù)，改進(jìn)完善污水處理設(shè)備，從而提高污水凈化率。5、合理進(jìn)行人工降雨，增加地表水和地下水的來源。6、采用高新技術(shù)，進(jìn)行更深層的地下水位挖掘。7、對大型水庫進(jìn)行維修處理，增加水庫的蓄水量。8、增加海水利用率。問題4關(guān)于水資源短缺風(fēng)險的建議報告北京市水行政主管部門：近年來，我國，特別是北方地區(qū)水資源短缺問題日趨嚴(yán)重。根據(jù)北京市人口普查相關(guān)數(shù)據(jù)顯示，1979年——2008年北京市人口年均增長率逐年增加，且從2001年北京市人口總量就已經(jīng)大大超過了北京市水資源的承載力，使得北京成為世界上嚴(yán)重缺水的大城市之一?，F(xiàn)對以上問題提出如下報告：一、北京市地處海河流域，有豐富的地下水資源，因此，北京一直以地下水為飲用水源。但是，隨著經(jīng)濟(jì)快速增長，人口增加，長期超量開采，導(dǎo)致地下水位下降。而地表水主要靠降雨補(bǔ)給，北京的湖泊都很小，水量有限。從而，北京水資源總量長期處于短缺狀態(tài)。北京市水行政主管部門應(yīng)采取措施進(jìn)行水資源的深度開發(fā)和優(yōu)化配置，開發(fā)新水源，實行水資源聯(lián)合調(diào)度。另外，可以通過調(diào)節(jié)和重新分配地表徑流、淡化海水和高礦化地下水以及把一部分地表徑流轉(zhuǎn)化為地下徑流的方法來解決水資源短缺問題。二、北京沒有天然湖泊，其大量用水靠水庫提供，但北京的密云水庫、官廳水庫、懷柔水庫、海子水庫四座大型水庫的水源入境量都有不同程度的衰減，近年來官廳水庫尤為嚴(yán)重。由于各年入境水量不穩(wěn)定，所以，應(yīng)該大力改進(jìn)修建水庫，擴(kuò)大蓄水量以保證水資源的合理利用。三、根據(jù)所給數(shù)據(jù)顯示，北京市近年來的工、農(nóng)業(yè)用水雖然都呈現(xiàn)下降趨勢，但第三產(chǎn)業(yè)用水卻呈現(xiàn)明顯增長趨勢。為此，必須從各方面狠抓節(jié)約用水，以節(jié)約用水支持社會經(jīng)濟(jì)發(fā)展。首先，要調(diào)整北京的工業(yè)結(jié)構(gòu)布局，使高新技術(shù)產(chǎn)業(yè)成為首都經(jīng)濟(jì)發(fā)展的主導(dǎo)。其次，要科學(xué)調(diào)整農(nóng)業(yè)種植結(jié)構(gòu)，建節(jié)水灌溉工程，并加強(qiáng)城鎮(zhèn)生活用水的管理。最后，實行用水總量控制和定額管理相結(jié)合的制度，盡快制定行業(yè)綜合用水定額、居民生活用水定額，在實踐中不斷完善和改進(jìn)用水方案，并研究科學(xué)可行的評估方法對各個用水單位進(jìn)行用水評估，且對評估結(jié)果進(jìn)行一定的獎懲措施。四、污水按性質(zhì)分為自然污染水和人為污染水。北京當(dāng)前對水體危害較大的是人為污染。北運河是北京天然河道五大水系中一條常年有水的河流，但由于城市污水污染、工業(yè)廢水污染、農(nóng)業(yè)回流水污染、固體廢物污染等使得北運河將在七年內(nèi)才能完成綜合治理。對于污水處理問題，可以加大污水處理廠的建設(shè)，也可以投入建設(shè)中水設(shè)施，以洗浴、盥洗等日常雜排水為水源，經(jīng)過處理達(dá)到中水水質(zhì)標(biāo)準(zhǔn)后，可以用于沖廁、洗車、綠化等。因此，建立科學(xué)合理的水資源利用效率指標(biāo)體系，研究科學(xué)可行的評估方法是很有必要的，對于節(jié)水型社會建設(shè)具有重要的社會意義。針對北京市的特殊背景和情況，上述建議為地方和行業(yè)水資源利用效率、調(diào)整產(chǎn)業(yè)結(jié)構(gòu)布局提供了信息支撐。以上報告，如無不妥，請批轉(zhuǎn)各地執(zhí)行。通化師范學(xué)院二0一一年五月五日模型評價與改進(jìn)本文提出的是使用模糊綜合評價的方法解決水資源短缺問題，與其它計算方法相比，最主要的特點是：1、模型的優(yōu)點：本模型基于多個指標(biāo)對北京市水資源短缺風(fēng)險進(jìn)行了綜合評價，具有一定的可操作性和實用性。本模型成功地分析了北京市水資源短缺風(fēng)險問題，其中采用了計算機(jī)數(shù)據(jù)對實際理論情況的分析，分析數(shù)據(jù)的結(jié)果表明此模型的建立有效合理。本模型與實際聯(lián)系較為密切，綜合實際情況進(jìn)行了合理分析，為北京市未來水資源短缺風(fēng)險的預(yù)測提供了有效的理論依據(jù)，也為解決水資源短缺問題提供了切實可行的方案。2、模型的缺點：由于現(xiàn)實生活中水資源系統(tǒng)的不確定性，所以，還有很多難以預(yù)料的不定因素的影響，這使得運算結(jié)果具有一定的誤差；同時，我們的理論體系還不夠完善，需要進(jìn)一步的改進(jìn)。3、模型的改進(jìn)：在計算概率時還可以應(yīng)用Logistic回歸方法進(jìn)一步提高模型的精度，由于時間關(guān)系，一些改進(jìn)的思路還來不及實現(xiàn)。模型推廣水資源是人類賴以生存和發(fā)展的重要資源之一，是不可或缺、不可代替的特殊資源。通過定性分析影響水資源短缺風(fēng)險因子，可以了解影響水資源短缺的因素，從而，從這些因素入手制定合理的方案降低水資源短缺現(xiàn)狀。模型推廣之后，可以對北京市未來水資源短缺風(fēng)險進(jìn)行預(yù)測并提出相應(yīng)的應(yīng)對措施，具有一定的實用性和可操作性。也可以借此方法來研究我國其他缺水地區(qū)水資源短缺風(fēng)險。 參考文獻(xiàn)[1]北京市統(tǒng)計局.北京2009年統(tǒng)計年鑒[M].中國統(tǒng)計,2006.[2]劉衛(wèi)國.MATLAB程序設(shè)計與應(yīng)用.北京：高等教育出版社，2002.[3]馮平.供水系統(tǒng)干旱期的水資源風(fēng)險管理[J].自然資源學(xué)報，1998,13（2）：139—144.[4]葉其孝.大學(xué)生數(shù)學(xué)建模競賽輔導(dǎo)教材.長沙：湖南教育出版社，1993.[5]姜啟源.數(shù)學(xué)模型（第二版）.高等教育出版社，1993. 附錄附表1979年至2008年北京市水資源短缺的狀況年份總用水量(億立方米)農(nóng)業(yè)用水(億立方米)工業(yè)用水(億立方米)第三產(chǎn)業(yè)及生活等其它用水(億立方米)水資源總量（億方）197942.9224.1814.374.3738.23198050.5431.8313.774.9426198148.1131.612.214.324198247.2228.8113.894.5236.6198347.5631.611.244.7234.7198440.0521.8414.3764.01739.31198531.7110.1217.24.3938198636.5519.469.917.1827.03198730.959.6814.017.2638.66198842.4321.9914.046.439.18198944.6424.4213.776.4521.55199041.1221.7412.347.0435.86199142.0322.711.97.4342.29199246.4319.9415.5110.9822.44199345.2220.3515.289.5919.67199445.8720.9314.5710.3745.42199544.8819.3313.7811.7730.34199640.0118.9511.769.345.87199740.3218.15199840.4317.3910.8412.237.7199941.7118.4510.5612.714.22200040.416.4910.5213.3916.86200138.919.2200234.615.57.511.616.1200335.813.88.413.618.4200434.613.57.713.421.4200534.513.26.814.523.2200634.324.5200734.812.45.816.623.8200835.112.05.217.934.2程序1（計算水資源風(fēng)險率衡量指標(biāo)）clcclearM=[42.9250.5448.1147.2247.5640.0531.7136.5530.9542.4344.6441.1242.0346.4345.2245.8744.8840.0140.3240.4341.7140.438.934.635.834.634.534.334.835.1];%總用水量G=[38.23262436.634.739.313827.0338.6639.1821.5535.8642.2922.4419.6745.4230.3445.8722.2537.714.2216.8619.216.118.421.423.224.523.834.2];%水資源總量C=M>G;D=M<GM(C)=1;M(D)=0It=M%水資源系統(tǒng)狀態(tài)量Tsx=0;fork=MTsx=Tsx+kendclcTsx%水資源狀態(tài)總量T=30%水資源系統(tǒng)工作的總歷時U1=Tsx/T%水資源風(fēng)險率衡量指標(biāo)運行結(jié)果：Tsx=26T=30U1=0.8667%總用水量與水資源總量分布圖%M=[42.9250.5448.1147.2247.5640.0531.7136.5530.9542.4344.6441.1242.0346.4345.2245.8744.8840.0140.3240.4341.7140.438.934.635.834.634.534.334.835.1];plot(M,'^')%總用水量分布圖%plot(G,'*')%水資源總量分布圖%程序2（計算水資源脆弱性衡量指標(biāo)）clcclearM=[42.9250.5448.1147.2247.5640.0531.7136.5530.9542.4344.6441.1242.0346.4345.2245.8744.8840.0140.3240.4341.7140.438.934.635.834.634.534.334.835.1];%總用水量G=[38.23262436.634.739.313827.0338.6639.1821.5535.8642.2922.4419.6745.4230.3445.8722.2537.714.2216.8619.216.118.421.423.224.523.834.2];%水資源總量C=M>G;D=M<GM(C)=1;M(D)=0T1=0;fork=MT1=T1+kendclcT1M=[42.9250.5448.1147.2247.5640.0531.7136.5530.9542.4344.6441.1242.0346.4345.2245.8744.8840.0140.3240.4341.7140.438.934.635.834.634.534.334.835.1];%總用水量G=[38.23262436.634.739.313827.0338.6639.1821.5535.8642.2922.4419.6745.4230.3445.8722.2537.714.2216.8619.216.118.421.423.224.523.834.2];%水資源總量QZt1=M-GB=QZt1<0;QZt1(B)=0;clcQZt1M=[42.9250.5448.1147.2247.5640.0531.7136.5530.9542.4344.6441.1242.0346.4345.2245.8744.8840.0140.3240.4341.7140.438.934.635.834.634.534.334.835.1];%總用水量G=[38.23262436.634.739.313827.0338.6639.1821.5535.8642.2922.4419.6745.4230.3445.8722.2537.714.2216.8619.216.118.421.423.224.523.834.2];%水資源總量QZt2=G-MB=QZt2<0;QZt2(B)=0;clcQZt2QZt=QZt1+QZt2QZ=sum(QZt1+QZt2)clcQZt%第t年水資源缺水量QZ%水資源短缺總量T1%水資源系統(tǒng)失事總時間X=QZ/T1%失事狀態(tài)下缺水量年均值SZt=G;%第t年水資源總量SZ=sum(SZt);%水資源總量T=30%水資源系統(tǒng)工作的總歷時Y=SZ/T%水資源總量年均值U2=X/Y%水資源脆弱性衡量指標(biāo)運行結(jié)果：QZt=Columns1through64.690024.540024.110010.620012.86000.7400Columns7through126.29009.52007.71003.250023.09005.2600Columns13through180.260023.990025.55000.450014.54005.8600Columns19through2418.07002.730027.490023.540019.700018.5000Columns25through3017.400013.200011.30009.800011.00000.9000QZ=376.9600T1=26X=14.4985T=30Y=29.2327U2=0.4960%第t年水資源缺水量與第t年水資源總量%plot(QZt,'g-')%第t年水資源缺水plot(SZt)%第t年水資源總量程序3（計算水資源重現(xiàn)期衡量指標(biāo)）clcclearM=[42.9250.5448.1147.2247.5640.0531.7136.5530.9542.4344.6441.1242.0346.4345.2245.8744.8840.0140.3240.4341.7140.438.934.635.834.634.534.334.835.1];%總用水量G=[38.23262436.634.739.313827.0338.6639.1821.5535.8642.2922.4419.6745.4230.3445.8722.2537.714.2216.8619.216.118.421.423.224.523.834.2];%水資源總量C=M>G;D=M<GM(C)=1;M(D)=0T1=0;fork=MT1=T1+kendclcT1It=M%水資源系統(tǒng)狀態(tài)量fori=1:29dqnm(i)=It(i+1)-It(i)endclcdqnmdqnn=dqnm<0;dqnm(dqnn)=1;dqn=sum(dqnm)/2clcdqn%第n與n+1次失事的總時間間隔U3=dqn/(T1-1)%水資源重現(xiàn)期衡量指標(biāo)運行結(jié)果：dqn=4U3=0.1600%失事的時間間隔%plot(dqnm,'x')%失事的時間間隔程序4（計算水資源可恢復(fù)性衡量指標(biāo)）clcclearM=[42.9250.5448.1147.2247.5640.0531.7136.5530.9542.4344.6441.1242.0346.4345.2245.8744.8840.0140.3240.4341.7140.438.934.635.834.634.534.334.835.1];%總用水量G=[38.23262436.634.739.313827.0338.6639.1821.5535.8642.2922.4419.6745.4230.3445.8722.2537.714.2216.8619.216.118.421.423.224.523.834.2];%水資源總量C=M>G;D=M<GM(C)=1;M(D)=0It=M%水資源系統(tǒng)狀態(tài)量fort=2:30ifIt(t-1)==1&&It(t)==0;Zt(t-1)=1;elseZt(t-1)=0;endendclcZt%水資源可恢復(fù)狀態(tài)量Thf=sum(Zt)%水資源可恢復(fù)狀態(tài)總量M=[42.9250.5448.1147.2247.5640.0531.7136.5530.9542.4344.6441.1242.0346.4345.2245.8744.8840.0140.3240.4341.7140.438.934.635.834.634.534.334.835.1];%總用水量G=[38.23262436.634.739.313827.0338.6639.1821.5535.8642.2922.4419.6745.4230.3445.8722.2537.714.2216.8619.216.118.421.423.224.523.834.2];%水資源總量C=M>G;D=M<GM(C)=1;M(D)=0Tsx=0;fork=MTsx=Tsx+kendclcTsx%水資源狀態(tài)總量ifTsx==0;U4=0;elseU4=Thf/Tsx;endU4%水資源可恢復(fù)性衡量指標(biāo)運行結(jié)果：Tsx=26U4=0.1538%水資源可恢復(fù)狀態(tài)量%plot(Zt,'+')%水資源可恢復(fù)狀態(tài)量程序5（計算水資源風(fēng)險度衡量指標(biāo)）clcclearM=[42.9250.5448.1147.2247.5640.0531.7136.5530.9542.4344.6441.1242.0346.4345.2245.8744.8840.0140.3240.4341.7140.438.934.635.834.634.534.334.835.1];%總用水量G=[38.23262436.634.739.313827.0338.6639.1821.5535.8642.2922.4419.6745.4230.3445.8722.2537.714.2216.8619.216.118.421.423.224.523.834.2];%水資源總量QZt1=M-GB=QZt1<0;QZt1(B)=0;clcQZt2=G-MB=QZt2<0;QZt2(B)=0;clcQZt=QZt1+QZt2%第t年水資源缺水量SZt=G;%第t年水資源總量SZ=sum(SZt);%水資源總量T=30%水資源系統(tǒng)工作的總歷時Y=SZ/T%水資源總量年均值QZc=QZt-Y%第t年缺水量與水資源總量年均值的差值QZcc=[QZc(1)^2QZc(2)^2QZc(3)^2QZc(4)^2QZc(5)^2QZc(6)^2QZc(7)^2QZc(8)^2QZc(9)^2QZc(10)^2QZc(11)^2QZc(12)^2QZc(13)^2QZc(14)^2QZc(15)^2QZc(16)^2QZc(17)^2QZc(18)^2QZc(19)^2QZc(20)^2QZc(21)^2QZc(22)^2QZc(23)^2QZc(24)^2QZc(25)^2QZc(26)^2QZc(27)^2QZc(28)^2QZc(29)^2QZc(30)^2];QZz=sum(QZcc);Dx=QZz/(T-1)%衡量風(fēng)險度的方差Ex=sqrt(Dx)%衡量風(fēng)險度的標(biāo)準(zhǔn)差U5=Ex/Y%水資源風(fēng)險度衡量指標(biāo)運行結(jié)果：QZt=Columns1through64.690024.540024.110010.620012.86000.7400Columns7through126.29009.52007.71003.250023.09005.2600Columns13through180.260023.990025.55000.450014.54005.8600Columns19through2418.07002.730027.490023.540019.700018.5000Columns25through3017.400013.200011.30009.800011.00000.9000T=30Y=29.2327QZc=Columns1through6-24.5427-4.6927-5.1227-18.6127-16.3727-28.4927Columns7through12-22.9427-19.7127-21.5227-25.9827-6.1427-23.9727Columns13through18-28.9727-5.2427-3.6827-28.7827-14.6927-23.3727Columns19through24-11.1627-26.5027-1.7427-5.6927-9.5327-10.7327Columns25through30-11.8327-16.0327-17.9327-19.4327-18.2327-28.3327Dx=362.5391Ex=19.0405U5=0.6513程序6（計算水資源短缺風(fēng)險級別）U1=0.8667;%水資源風(fēng)險率衡量指標(biāo)U2=0.4695;%水資源脆弱性衡量指標(biāo)U3=0.1600;%水資源重現(xiàn)期衡量指標(biāo)U4=0.1538;%水資源可恢復(fù)性衡量指標(biāo)U5=0.6513;%水資源風(fēng)險度衡量指標(biāo)U=[U1U2U3U4U5]%風(fēng)險因子衡量指標(biāo)的集合miuv1phi11=0.2miuv2phi11=0.3miuv3phi11=0.5miuv4phi11=0.7miuv5phi11=0.8%風(fēng)險率、脆弱性、重現(xiàn)期、可恢復(fù)性、風(fēng)險度低風(fēng)險的指標(biāo)miuv1phi12=0.2miuv2phi12=0.3miuv3phi12=0.5miuv4phi12=0.7miuv5phi12=0.8%風(fēng)險率、脆弱性、重現(xiàn)期、可恢復(fù)性、風(fēng)險度較低風(fēng)險的指標(biāo)miuv1phi13=0.8miuv2phi13=0.7miuv3phi13=0.5miuv4phi13=0.3miuv5phi13=0.2%風(fēng)險率、脆弱性、重現(xiàn)期、可恢復(fù)性、風(fēng)險度中風(fēng)險的指標(biāo)miuv1phi14=9miuv2phi14=7.5miuv3phi14=4.5miuv4phi14=2miuv5phi14=1%風(fēng)險率、脆弱性、重現(xiàn)期、可恢復(fù)性、風(fēng)險度較高風(fēng)險的指標(biāo)miuv1phi15=0.2miuv2phi15=0.4miuv3phi15=0.8miuv4phi15=1.5miuv5phi15=2%風(fēng)險率、脆弱性、重現(xiàn)期、可恢復(fù)性、風(fēng)險度高風(fēng)險的指標(biāo)W1=[miuv1phi11miuv2phi11miuv3phi11miuv4phi11miuv5phi11]%風(fēng)險率U1衡量指標(biāo)W2=[miuv1phi12miuv2phi12miuv3phi12miuv4phi12miuv5phi12]%脆弱性U2衡量指標(biāo)W3=[miuv1phi13miuv2phi13miuv3phi13miuv4phi13miuv5phi13]%可恢復(fù)性U3衡量指標(biāo)W4=[miuv1phi14miuv2phi14miuv3phi14miuv4phi14miuv5phi14]%重現(xiàn)期U4衡量指標(biāo)W5=[miuv1phi15miuv2phi15miuv3phi15miuv4phi15miuv5phi15]%風(fēng)險度U5衡量指標(biāo)clcRu0=[W1;W2;W3;W4;W5]%變換關(guān)系矩陣ifU1<=miuv1phi11miuv1phi11=1;elseifU1<miuv2phi11&&U1>=miuv1phi11miuv2phi11=(miuv2phi11-U1)/(miuv2phi11-miuv1phi11);elseifU1<miuv3phi11&&U1>=miuv2phi11miuv3phi11=0;elseifU1<miuv4phi11&&U1>=miuv3phi11miuv4phi11=0;elseifU1>=miuv1phi11miuv5phi11=0;endW1=[miuv1phi11miuv2phi11miuv3phi11miuv4phi11miuv5phi11]W2=[miuv1phi12miuv2phi12miuv3phi12miuv4phi12miuv5phi12]W3=[miuv1phi13miuv2phi13miuv3phi13miuv4phi13miuv5phi13]W4=[miuv1phi14miuv2phi14miuv3phi14miuv4phi14miuv5phi14]W5=[miuv1phi15miuv2phi15miuv3phi15miuv4phi15miuv5phi15]clcRu1=[W1;W2;W3;W4;W5]%只對風(fēng)險率進(jìn)行了衡量的變換關(guān)系矩陣ifU2<=miuv1phi12miuv1phi12=U2/miuv1phi12;elseifU2<miuv2phi12&&U2>=miuv1phi12miuv1phi12=1;elseifU2<miuv3phi12&&U2>=miuv2phi12miuv3phi12=(miuv3phi12+0.1-U2)/(miuv3phi12-miuv2phi12);elseifU2<miuv4phi12&&U2>=miuv3phi12miuv4phi12=0;elseifmiuv1phi12<=U2miuv5phi12=0;endW1=[miuv1phi11miuv2phi11miuv3phi11miuv4phi11miuv5phi11]W2=[miuv1phi12miuv2phi12miuv3phi12miuv4phi12miuv5phi12]W3=[miuv1phi13miuv2phi13miuv3phi13miuv4phi13miuv5phi13]W4=[miuv1phi14miuv2phi14miuv3phi14miuv4phi14miuv5phi14]W5=[miuv1phi15miuv2phi15miuv3phi15miuv4phi15miuv5phi15]clcRu2=[W1;W2;W3;W4;W5]%對風(fēng)險率、脆弱性進(jìn)行了衡量的變換關(guān)系矩陣clcifU3>=miuv1phi13miuv1phi13=0;elseifU3<miuv2phi13&&U3>=miuv1phi13miuv2phi13=miuv2phi13/U3;elseifU3<miuv3phi13&&U3>=miuv2phi13miuv3phi13=1;elseifU3<miuv4phi13&&U3>=miuv3phi13miuv4phi13=(miuv3phi13-U3)/(miuv3phi13-miuv4phi13);elseifU3<=miuv4phi13miuv5phi13=0;endW1=[miuv1phi11miuv2phi11miuv3phi11miuv4phi11miuv5phi11]W2=[miuv1phi12miuv2phi12miuv3phi12miuv4phi12miuv5phi12]W3=[miuv1phi13miuv2phi13miuv3phi13miuv4phi13miuv5phi13]W4=[miuv1phi14miuv2phi14miuv3phi14miuv4phi14miuv5phi14]W5=[miuv1phi15miuv2phi15miuv3phi15miuv4phi15miuv5phi15]clcRu3=[W1;W2;W3;W4;W5]%對風(fēng)險率、脆弱性、重現(xiàn)期進(jìn)行了衡量的變換關(guān)系矩陣clcifU4>=miuv2phi14miuv1phi14=0;elseifU4<miuv2phi14&&U4>=miuv3phi14miuv3phi14=miuv3phi14/U4;elseifU4<miuv3phi14&&U4>=miuv4phi14miuv3phi14=1;elseifU4<miuv4phi14miuv5phi14=(U4)/miuv5phi14;endW1=[miuv1phi11miuv2phi11miuv3phi11miuv4phi11miuv5phi11]W2=[miuv1phi12miuv2phi12miuv3phi12miuv4phi12miuv5phi12]W3=[miuv1phi13miuv2phi13miuv3phi13miuv4phi13miuv5phi13]W4=[miuv1phi14miuv2phi14miuv3phi14miuv4phi14miuv5phi14]W5=[miuv1phi15miuv2phi15miuv3phi15miuv4phi15miuv5phi15]clcRu4=[W1;W2;W3;W4;W5]%對風(fēng)險率、脆弱性、重現(xiàn)期、可恢復(fù)性進(jìn)行了衡量的變換關(guān)系矩陣clcifU5<=miuv3phi15miuv3phi15=0;elseifU5<miuv3phi15&&U5>=miuv4phi15miuv4phi15=miuv4phi15/U5;elseifU5>miuv4phi15miuv5phi15=1;endW1=[miuv1phi11miuv2phi11miuv3phi11miuv4phi11miuv5phi11]W2=[miuv1phi12miuv2phi12miuv3phi12miuv4phi12miuv5phi12]W3=[miuv1phi13miuv2phi13miuv3phi13miuv4phi13miuv5phi13]W4=[miuv1phi14miuv2phi14miuv3phi14miuv4phi14miuv5phi14]W5=[miuv1phi15miuv2phi15miuv3phi15miuv4phi15miuv5phi15]clcRu5=[W1;W2;W3;W4;W5]%對風(fēng)險率、脆弱性、重現(xiàn)期、可恢復(fù)性、風(fēng)險度進(jìn)行了衡量的變換關(guān)系矩陣clcRu=Ru5A=[1/21/41/81/121/24]%水資源短缺風(fēng)險各因素的權(quán)重系數(shù)B=A*Ru%水資源短缺風(fēng)險評價等級運行結(jié)果：Ru0=0.20000.30000.50000.70000.80000.20000.30000.50000.70000.80000.80000.70000.50000.30000.20009.00007.50004.50002.00001.00000.20000.40000.80001.50002.0000Ru1=0.20000.30000.50000.700000.20000.30000.50000.70000.80000.80000.70000.50000.30000.20009.00007.50004.50002.00001.00000.20000.40000.80001.50002.0000Ru2=0.20000.30000.50000.700000.20000.30000.65250.70000.80000.80000.70000.50000.30000.20009.00007.50004.50002.00001.00000.20000.40000.80001.50002.0000Ru3=0.20000.30000.50000.700000.20000.30000.65250.70000.80000.80000.70000.50000.300009.00007.50004.50002.00001.00000.20000.40000.80001.50002.0000Ru4=0.20000.30000.50000.700000.20000.30000.65250.70000.80000.80000.70000.50000.300009.00007.50004.50002.00000.15380.20000.40000.80001.50002.0000Ru5=0.20000.30000.50000.700000.20000.30000.65250.70000.80000.80000.70000.50000.300009.00007.50004.50002.00000.15380.20000.400001.50002.0000A=0.50000.25000.12500.08330.0417B=1.00830.95420.85060.79170.2962%水資源短缺風(fēng)險各因素的權(quán)重系數(shù)、變換關(guān)系矩陣與水資源短缺風(fēng)險級別%plot(A,'s')%水資源短缺風(fēng)險各因素的權(quán)重系數(shù)plot(B,'r-d')%水資源短缺風(fēng)險級別plot(Ru)%變換關(guān)系矩陣plot(U,'c-*')%風(fēng)險因子衡量指標(biāo)的集合程序7（計算水資源短缺風(fēng)險衡量值）U1=0.8667;%水資源風(fēng)險率衡量指標(biāo)U2=0.4695;%水資源脆弱性衡量指標(biāo)U3=0.1600;%水資源重現(xiàn)期衡量指標(biāo)U4=0.1538;%水資源可恢復(fù)性衡量指標(biāo)U5=0.6513;%水資源風(fēng)險度衡量指標(biāo)clcB=[1.00830.95420.85060.79170.2962]%水資源短缺風(fēng)險評價等級A=[1/21/41/81/121/24]%水資源短缺風(fēng)險各因素的權(quán)重系數(shù)U=[U1U2U3U4U5]%風(fēng)險因子衡量指標(biāo)的集合UZ=A*U'%水資源短缺風(fēng)險衡量值運行結(jié)果：A=0.50000.25000.12500.08330.0417U=0.86670.46950.16000.15380.6513UZ=0.6107%水資源短缺風(fēng)險衡量值在水資源短缺風(fēng)險評價等級中的顯示%plot(B,'y-<')holdonplot(UZ,'c-*')%水資源短缺風(fēng)險衡量值在水資源短缺風(fēng)險評價等級中的顯示基于C8051F單片機(jī)直流電動機(jī)反饋控制系統(tǒng)的設(shè)計與研究基于單片機(jī)的嵌入式Web服務(wù)器的研究MOTOROLA單片機(jī)MC68HC（8）05PV8/A內(nèi)嵌EEPROM的工藝和制程方法及對良率的影響研究基于模糊控制的電阻釬焊單片機(jī)溫度控制系統(tǒng)的研制基于MCS-51系列單片機(jī)的通用控制模塊的研究基于單片機(jī)實現(xiàn)的供暖系統(tǒng)最佳啟停自校正（STR）調(diào)節(jié)器單片機(jī)控制的二級倒立擺系統(tǒng)的研究基于增強(qiáng)型51系列單片機(jī)的TCP/IP協(xié)議棧的實現(xiàn)基于單片機(jī)的蓄電池自動監(jiān)測系統(tǒng)基于32位嵌入式單片機(jī)系統(tǒng)的圖像采集與處理技術(shù)的研究基于單片機(jī)的作物營養(yǎng)診斷專家系統(tǒng)的研究基于單片機(jī)的交流伺服電機(jī)運動控制系統(tǒng)研究與開發(fā)基于單片機(jī)的泵管內(nèi)壁硬度測試儀的研制基于單片機(jī)的自動找平控制系統(tǒng)研究基于C8051F040單片機(jī)的嵌入式系統(tǒng)開發(fā)基于單片機(jī)的液壓動力系統(tǒng)狀態(tài)監(jiān)測儀開發(fā)模糊Smith智能控制方法的研究及其單片機(jī)實現(xiàn)一種基于單片機(jī)的軸快流CO〈,2〉激光器的手持控制面板的研制基于雙單片機(jī)沖床數(shù)控系統(tǒng)的研究基于CYGNAL單片機(jī)的在線間歇式濁度儀的研制基于單片機(jī)的噴油泵試驗臺控制器的研制基于單片機(jī)的軟起動器的研究和設(shè)計基于單片機(jī)控制的高速快走絲電火花線切割機(jī)床短循環(huán)走絲方式研究基于單片機(jī)的機(jī)電產(chǎn)品控制系統(tǒng)開發(fā)基于PIC單片機(jī)的智能手機(jī)充電器基于單片機(jī)的實時內(nèi)核設(shè)計及其應(yīng)用研究基于單片機(jī)的遠(yuǎn)程抄表系統(tǒng)的設(shè)計與研究基于單片機(jī)的煙氣二氧化硫濃度檢測儀的研制基于微型光譜儀的單片機(jī)系統(tǒng)單片機(jī)系統(tǒng)軟件構(gòu)件開發(fā)的技術(shù)研究基于單片機(jī)的液體點滴速度自動檢測儀的研制基于單片機(jī)系統(tǒng)的多功能溫度測量儀的研制基于PIC單片機(jī)的電能采集終端的設(shè)計和應(yīng)用基于單片機(jī)的光纖光柵解調(diào)儀的研制氣壓式線性摩擦焊機(jī)單片機(jī)控制系統(tǒng)的研制基于單片機(jī)的數(shù)字磁通門傳感器基于單片機(jī)的旋轉(zhuǎn)變壓器-數(shù)字轉(zhuǎn)換器的研究基于單片機(jī)的光纖Bragg光柵解調(diào)系統(tǒng)的研究單片機(jī)控制的便攜式多功能乳腺治療儀的研制基于C8051F020單片機(jī)的多生理信號檢測儀基于單片機(jī)的電機(jī)運動控制系統(tǒng)設(shè)計Pico專用單片機(jī)核的可測性設(shè)計研究基于MCS-51單片機(jī)的熱量計基于雙單片機(jī)的智能遙測微型氣象站MCS-51單片機(jī)構(gòu)建機(jī)器人的實踐研究基于單片機(jī)的輪軌力檢測基于單片機(jī)的GPS定位儀的研究與實現(xiàn)基于單片機(jī)的電液伺服控制系統(tǒng)用于單片機(jī)系統(tǒng)的MMC卡文件系統(tǒng)研制基于單片機(jī)的時控和計數(shù)系統(tǒng)性能優(yōu)化的研究基于單片機(jī)和CPLD的粗光柵位移測量系統(tǒng)研究單片機(jī)控制的后備式方波UPS提升高職學(xué)生單片機(jī)應(yīng)用能力的探究基于單片機(jī)控制的自動低頻減載裝置研究基于單片機(jī)控制的水下焊接電源的研究基于單片機(jī)的多通道數(shù)據(jù)采集系統(tǒng)基于uPSD3234單片機(jī)的氚表面污染測量儀的研制基于單片機(jī)的紅外測油儀的研究96系列單片機(jī)仿真器研究與設(shè)計基于單片機(jī)的單晶金剛石刀具刃磨設(shè)備的數(shù)控改造基于單片機(jī)的溫度智能控制系統(tǒng)的設(shè)計與實現(xiàn)基于MSP430單片機(jī)的電梯門機(jī)控制器的研制基于單片機(jī)的氣體測漏儀的研究基于三菱M16C/6N系列單片機(jī)的CAN/USB協(xié)議轉(zhuǎn)換器基于單片機(jī)和DSP的變壓器油色譜在線監(jiān)測技術(shù)研究基于單片機(jī)的膛壁溫度報警系統(tǒng)設(shè)計基于AVR單片機(jī)的低壓無功補(bǔ)償控制器的設(shè)計基于單片機(jī)船舶電力推進(jìn)電機(jī)監(jiān)測系統(tǒng)基于單片機(jī)網(wǎng)絡(luò)的振動信號的采集系統(tǒng)基于單片機(jī)的大容量數(shù)據(jù)存儲技術(shù)的應(yīng)用研究基于單片機(jī)的疊圖機(jī)研究與教學(xué)方法實踐基于單片機(jī)嵌入式Web服務(wù)器技術(shù)的研究及實現(xiàn)基于AT89S52單片機(jī)的通用數(shù)據(jù)采集系統(tǒng)基于單片機(jī)的多道脈沖幅度分析儀研究機(jī)器人旋轉(zhuǎn)電弧傳