2018東三省數(shù)學(xué)建模聯(lián)賽E優(yōu)秀論文

《2018東三省數(shù)學(xué)建模聯(lián)賽E優(yōu)秀論文》由會(huì)員分享，可在線閱讀，更多相關(guān)《2018東三省數(shù)學(xué)建模聯(lián)賽E優(yōu)秀論文（36頁(yè)珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。

封一 答卷編號(hào) 參賽學(xué)校填寫 答卷編號(hào) 競(jìng)賽組委會(huì)填寫 論文題目 E 組 別 本科生 參賽隊(duì)員信息 必填 姓 名 學(xué) 號(hào) 聯(lián)系電話 參賽學(xué)校 封二 答卷編號(hào) 參賽學(xué)校填寫 答卷編號(hào) 競(jìng)賽組委會(huì)填寫 評(píng)閱情況 學(xué)校評(píng)閱專家填寫 學(xué)校評(píng)閱 1 學(xué)校評(píng)閱 2 學(xué)校評(píng)閱 3 評(píng)閱情況 聯(lián)賽評(píng)閱專家填寫 聯(lián)賽評(píng)閱 1 聯(lián)賽評(píng)閱 2 聯(lián)賽評(píng)閱 3 E 我國(guó)水資源現(xiàn)狀的分析發(fā)展及預(yù)測(cè) 摘要 本文通過(guò)近五年我國(guó)水資源相關(guān)信息與數(shù)據(jù) 進(jìn)行科學(xué)性的分析研究 判斷出我 國(guó)水資源現(xiàn)狀及其存在的問(wèn)題 如何科學(xué)有效的管理規(guī)劃未來(lái)水資源的管理 以及和 水資源指標(biāo)之間的內(nèi)在聯(lián)系 為今后我國(guó)水資源的可持續(xù)發(fā)展起到積極作用 針對(duì)問(wèn)題一 我們收集了 2012 2016 年的水資源公報(bào)數(shù)據(jù)信息 通過(guò)列表做折線 圖 我們發(fā)現(xiàn)我國(guó)用水量呈上升趨勢(shì) 并且預(yù)測(cè)在 2030 年會(huì)超過(guò)水資源開(kāi)發(fā)利用中全 國(guó)用水總量紅線 7000 億立方米 為此我們分析了各個(gè)指標(biāo) 并且結(jié)合我國(guó)國(guó)情和特點(diǎn) 發(fā)現(xiàn)了在水資源問(wèn)題中影響因素較大的因子 我國(guó)是農(nóng)業(yè)大國(guó) 全年總用水量約有近 50 用來(lái)農(nóng)業(yè)發(fā)展 并且這個(gè)數(shù)字仍在不斷上升 我們發(fā)現(xiàn) 盡管我國(guó)每年的農(nóng)業(yè)用水 量極大 但是農(nóng)田灌溉水有效利用系數(shù)卻很低 這導(dǎo)致了極大的浪費(fèi) 與我國(guó)短缺的 水資源問(wèn)題相悖 因此 我們著手分析農(nóng)田灌溉水有效利用系數(shù)這一數(shù)據(jù)進(jìn)行分析 假設(shè)農(nóng)田灌溉用水有效利用率系數(shù)函數(shù)為 y aebt 經(jīng) Excel 擬合得到農(nóng)田灌溉用水有 效利用率系數(shù)函數(shù) y 0 510e0 012x 從而預(yù)測(cè) 2030 年農(nóng)田灌溉用水有效利用率系數(shù)偏 差為 2 6 針對(duì)問(wèn)題二 首先我們明確的農(nóng)田灌溉水有效利用系數(shù) 萬(wàn)元工業(yè)增加用水量以 及用水效率的定義函數(shù)表達(dá)式 聯(lián)系分析水資源開(kāi)發(fā)利用數(shù)據(jù) 得到農(nóng)田灌溉水有效 利用系數(shù) 萬(wàn)元工業(yè)增加用水量以及用水效率相關(guān)模型 4 5 我們收集廣東省 2008 2016 年的數(shù)據(jù) 使用最小二乘法進(jìn)行擬合得到工業(yè)和農(nóng)業(yè)用水占比系數(shù)與其他指 標(biāo)的函數(shù)關(guān)系表達(dá)式 通過(guò)使用 MATLAB 待定系數(shù)法進(jìn)行在保證目標(biāo)的前提下的擬合 直觀得到其余各個(gè)指標(biāo)的變化趨勢(shì) 解決多變量不能一一求解的問(wèn)題 得到要確保每 年要以約 1 57 的萬(wàn)元工業(yè)增加值用水量才能達(dá)到目標(biāo)這一結(jié)論 針對(duì)問(wèn)題三 我們利用灰色關(guān)聯(lián)度矩陣 利用 MATLAB 軟件的處結(jié)果圖形與關(guān)聯(lián)矩 陣 通過(guò)圖像直觀分析該關(guān)聯(lián)矩陣 R 初得出萬(wàn)元工業(yè)增加值用水量 萬(wàn)元國(guó)內(nèi)生產(chǎn)總 值用水量即農(nóng)田灌溉水有效利用系數(shù)的內(nèi)在規(guī)律 再利用 MATLAB 的 K means 函數(shù) 以 萬(wàn)元工業(yè)增加值用水量 萬(wàn)元國(guó)內(nèi)生產(chǎn)總值用水量即農(nóng)田灌溉水有效利用系數(shù)為影響 因子 將全國(guó) 31 省按各省相對(duì)于 K means 函數(shù)中聚點(diǎn)的相對(duì)位置分層 選擇較為具有 代表性的省 導(dǎo)入數(shù)據(jù)到 SPSS 中進(jìn)行差分 進(jìn)一步分析得出結(jié)果 問(wèn)題四則是對(duì)以上結(jié)果的總結(jié)與概括 對(duì)每個(gè)問(wèn)題進(jìn)行了評(píng)價(jià)分析 結(jié)合我國(guó)國(guó) 情和當(dāng)今發(fā)展現(xiàn)狀 給國(guó)家水利部提交一份建議書 為今后我國(guó)水資源發(fā)展做出貢獻(xiàn) 關(guān)鍵字 MATLAB 用水效率 待定系數(shù)法 最小二乘法 灰色關(guān)聯(lián)度矩陣 1 一 問(wèn)題重述 1 1 問(wèn)題背景 水是生命之源 新中國(guó)成立以來(lái)特別是改革開(kāi)放以來(lái) 我國(guó)的水資源開(kāi)發(fā)利用 配置保護(hù)和管理等工作取得顯著成績(jī) 為社會(huì)發(fā)展做出了巨大貢獻(xiàn) 但是我國(guó)的基本 國(guó)情和水情產(chǎn)生了當(dāng)前人多水少 水資源時(shí)空分布不均勻的問(wèn)題 并且還有水資源短 缺 污染 水資源惡化等一系列問(wèn)題十分突出 這已成為當(dāng)前我國(guó)實(shí)現(xiàn)可持續(xù)發(fā)展戰(zhàn) 略的主要問(wèn)題 隨著社會(huì)發(fā)展 水資源需求將在較長(zhǎng)一段時(shí)間內(nèi)持續(xù)增長(zhǎng) 供需矛盾 將更加尖銳 形勢(shì)也將十分嚴(yán)峻 為解決水資源問(wèn)題 2011 年中央 1 號(hào)文件和中央水利工作會(huì)議明確要求實(shí)行最嚴(yán) 格水資源管理制度 2012 年國(guó)務(wù)院發(fā)布 關(guān)于實(shí)行最嚴(yán)格水資源管理制度的意見(jiàn) 每年發(fā)布 中國(guó)水資源公報(bào) 確立水資源開(kāi)發(fā)利用 用水效率控制和水功能區(qū)限制 納污 三條紅線 1 確立水資源開(kāi)發(fā)利用控制紅線 到 2030 年全國(guó)用水總量控制在 7000 億立方米 以內(nèi) 2 確立用水效率控制紅線 到 2030 年用水效率達(dá)到或接近世界先進(jìn)水平 萬(wàn)元工 業(yè)增加值降低到 40 億立方米以下 農(nóng)田灌溉水有效利用系數(shù)提高到 0 6 以上 3 確立水功能區(qū)限制納污紅線 到 2030 年主要污染物入河湖總量控制在水功能區(qū) 納污能力范圍之內(nèi) 水功能區(qū)水質(zhì)達(dá)標(biāo)率提高到 95 以上 1 2 問(wèn)題概述 現(xiàn)要求收集近五年相關(guān)信息數(shù)據(jù)完成問(wèn)題 1 根據(jù) 2012 2016 我國(guó)水資源公報(bào) 分析我國(guó)水資源開(kāi)發(fā)利用現(xiàn)狀和問(wèn)題并與國(guó) 家提出的 2020 年和 2030 年目標(biāo)偏差并進(jìn)行對(duì)比 2 構(gòu)建相關(guān)數(shù)學(xué)模型 如何保證在不超過(guò)水資源開(kāi)發(fā)利用量和用水效率控制紅線 的前提下實(shí)現(xiàn)萬(wàn)元工業(yè)增加用水量呈下降趨勢(shì)以及農(nóng)田灌溉水利用系數(shù)提高到 0 6 以 上的目標(biāo) 3 構(gòu)建相關(guān)數(shù)學(xué)模型 分析不同省的萬(wàn)元國(guó)內(nèi)生產(chǎn)總值用水量 耕地實(shí)際灌溉畝 均用水量和萬(wàn)元工業(yè)增加值的內(nèi)在規(guī)律并預(yù)測(cè) 2020 至 2030 年各省水管理控制發(fā)展趨 勢(shì) 4 根據(jù)研究成果寫一封建議信 2 二 問(wèn)題分析 問(wèn)題一需要分析 2012 2016 這五年的水資源公報(bào)數(shù)據(jù)信息 分析我國(guó)水資源各類 類別數(shù)據(jù) 進(jìn)行列表做折線圖 分析實(shí)際情況和存在的主要問(wèn)題 并與國(guó)家提出的 2020 和 2030 年的目標(biāo)進(jìn)行對(duì)比 同時(shí)進(jìn)行線性擬合農(nóng)田灌溉水有效利用系數(shù)這一重要 數(shù)據(jù)的線性函數(shù) 分析到 2030 年的數(shù)據(jù)計(jì)算偏差 對(duì)此我們收集相關(guān)數(shù)據(jù)并制作表格 和折線圖 以明確清晰地表示近幾年數(shù)據(jù)的變化趨勢(shì) 以進(jìn)行對(duì)比分析 問(wèn)題二則需要先明確各個(gè)變量的函數(shù)關(guān)系表達(dá)式 并綜合到一起聯(lián)立得到個(gè)指標(biāo) 與工農(nóng)業(yè)用水占比之間的函數(shù)關(guān)系式 首先需收集一個(gè)省份的數(shù)據(jù) 先得到此省近年 來(lái)工業(yè)與農(nóng)業(yè)用水占比系數(shù)的 建立散點(diǎn)圖 并進(jìn)行線性擬合預(yù)測(cè) 由此確立 2030 年 的工農(nóng)業(yè)用水占比系數(shù) 由此可以確立出在工農(nóng)業(yè)用水占比系數(shù)確定的情況下其他指 標(biāo)的函數(shù)關(guān)系式 為明確直觀分析出其他標(biāo)量的變化趨勢(shì) 我們利用 MATLAB 來(lái)進(jìn)行此 函數(shù)的三維圖像擬合 進(jìn)而分析得出結(jié)論 問(wèn)題三需要我們根據(jù)不同省份的特點(diǎn)進(jìn)行分析研究萬(wàn)元國(guó)內(nèi)生產(chǎn)總值用水量 耕 地實(shí)際灌溉畝均用水量和萬(wàn)元工業(yè)增加值用水量之間存在的內(nèi)在規(guī)律 我們需收集各 省份的近幾年的數(shù)據(jù) 利用灰色預(yù)測(cè)中關(guān)聯(lián)度模型講各省數(shù)據(jù)聚合到一起 利用特定 函數(shù) 以題中所需的萬(wàn)元國(guó)內(nèi)生產(chǎn)總值用水量 耕地實(shí)際灌溉畝均用水浪和萬(wàn)元工業(yè) 增加值用水量三個(gè)量為影響因子 選出具有明顯代表性特點(diǎn)的省份進(jìn)行擬合 進(jìn)一步 分析從而得出結(jié)果 預(yù)測(cè)分析 2020 2030 年各省的水資源管理及控制的發(fā)展趨勢(shì) 3 模型假設(shè) 1 假設(shè)從 2016 起之后 不出現(xiàn)自然天災(zāi) 瘟疫等重大影響自然資源事件 2 假設(shè)我國(guó)在未來(lái)的年份內(nèi)不發(fā)生重大經(jīng)濟(jì)轉(zhuǎn)型 能源戰(zhàn)略維持不變 3 假設(shè)所給水資源公報(bào)各個(gè)數(shù)據(jù)準(zhǔn)確無(wú)誤 4 假設(shè)全國(guó)總用水量年均增長(zhǎng)率為 1 01 5 假設(shè)我國(guó)未來(lái)幾年無(wú)任何重大水資源調(diào)遣戰(zhàn)略變化 6 假設(shè)我國(guó)各用水指標(biāo)仍呈有利趨勢(shì)發(fā)展 3 4 模型的建立與求解 4 1 問(wèn)題一 4 1 1 符號(hào)說(shuō)明 4 1 2 問(wèn)題概述 根據(jù) 2012 2016 我國(guó)水資源公報(bào) 分析我國(guó)水資源開(kāi)發(fā)利用現(xiàn)狀和問(wèn)題并與國(guó)家 提出的 2020 年和 2030 年目標(biāo)偏差并進(jìn)行對(duì)比 4 1 3 模型的建立與求解 以下數(shù)據(jù)是根據(jù) 2012 2016 年中華人民共和國(guó)水利部水資源公報(bào)所收集的部分?jǐn)?shù)據(jù) 根據(jù)數(shù)據(jù)做出部分對(duì)比圖像如下 1 用水量與供水量占當(dāng)年水資源總量百分比折線圖 t 時(shí)間 年份 a 待定系數(shù) b 待定系數(shù) A 偏差 y 2030 年的函數(shù)值 4 2012 2013 2014 2015 2016 2020 2030 0 1000 2000 3000 4000 5000 6000 7000 8000 19 00 19 50 20 00 20 50 21 00 21 50 22 00 22 50 23 00 用 水 量 億 供 水 量 占 當(dāng) 年 水 資 源 總 量 由此圖我們可以看出 2012 到 2016 年全國(guó)供水量主要呈下降趨勢(shì) 但是通過(guò)計(jì) 算往年增長(zhǎng)率得知 現(xiàn)我國(guó)總用水量年均增長(zhǎng)率為 1 01 由此可預(yù)測(cè)出 2020 年 與 2030 年我國(guó)總用水量分別為 6634 億 m 和 7335 4m 再次做出折線圖 6131 2 6183 4 6095 6103 2 6040 2 6634 7335 4 2012 2013 2014 2015 2016 2020 2030 0 1000 2000 3000 4000 5000 6000 7000 8000 用 水 量 億m 用 水 量 億m 我們發(fā)現(xiàn)隨著人口增長(zhǎng)和經(jīng)濟(jì)的持續(xù)高速發(fā)展 到 2020 年和 2030 年 全國(guó)用水 量較往年有很大增長(zhǎng) 并且在 2030 年總用水量超過(guò)了水資源開(kāi)發(fā)利用紅線即 2030 年 控制在 7000 億立方米以內(nèi) 這可以說(shuō)明 我國(guó)仍面臨著水資源供需嚴(yán)重緊迫的問(wèn)題 我國(guó)現(xiàn)水資源總量為 2 8 萬(wàn)億立方米 居世界第六位 但人均占有量為 2240 立方 5 米 約為世界人均水平的 1 4 在世界銀行連續(xù)統(tǒng)計(jì)的 153 個(gè)國(guó)家中居第 88 位 總量 并不豐富 人均占有量低 水資源分配不均是我國(guó)水資源現(xiàn)狀的一部分 并且隨著社 會(huì)發(fā)展人口在不斷增長(zhǎng)我國(guó)人均占有水資源量仍面臨著挑戰(zhàn) 2 農(nóng)田灌溉水有效利用系數(shù)與萬(wàn)元工業(yè)增加值用水量折線圖 0 516 0 523 0 53 0 536 0 542 2012 2013 2014 2015 2016 0 51 0 515 0 52 0 525 0 53 0 535 0 54 0 545 69 67 59 5 58 3 52 8 0 10 20 30 40 50 60 70 80 農(nóng) 田 灌 溉 水 有 效 利 用 系 數(shù) 萬(wàn) 元 工 業(yè) 增 加 值 用 水 量 工 業(yè) 用 水 量 立 方 米 工 業(yè) 增 加 值 萬(wàn) 元 我們由此圖看到 近幾年的數(shù)據(jù)表示我國(guó)萬(wàn)元工業(yè)增加值用水量呈總體下降趨勢(shì) 我國(guó)近兩年工業(yè)增加值隨著國(guó)家經(jīng)濟(jì)實(shí)力的不斷進(jìn)步定會(huì)不斷增長(zhǎng) 但這不意味著萬(wàn) 元工業(yè)增加值用水量也在下降 由分子與分母的關(guān)系即可得出 據(jù)了解 目前我國(guó)工 業(yè)萬(wàn)元產(chǎn)值用水量約為 80 億立方米 是發(fā)達(dá)國(guó)家的 10 20 倍 我國(guó)水的重復(fù)利用率為 40 左右 而發(fā)達(dá)國(guó)家為 75 85 左右 工業(yè)用水量巨大且重復(fù)利用率低 是我們水資 源開(kāi)發(fā)利用存在的一大問(wèn)題 由圖還知 我國(guó)近五年農(nóng)田灌溉水有效利用系數(shù)逐年上漲 但仍在 0 5 0 55 之間 部分地區(qū)可能只有 0 3 0 4 左右 可見(jiàn)我國(guó)農(nóng)田灌溉水有效利用率之低 發(fā)達(dá)國(guó)家早 在 40 50 年代就開(kāi)始采用節(jié)水灌溉 現(xiàn)在很多國(guó)家實(shí)現(xiàn)了輸水渠道防滲化 管道化 大田噴灌 滴灌化 灌溉科學(xué)化 自動(dòng)化 灌溉水利用系數(shù)達(dá)到 0 7 0 8 經(jīng)對(duì)比我國(guó) 水資源的浪費(fèi)很嚴(yán)重 這是目前我國(guó)水資源現(xiàn)狀之一 3 廢污水排放量折線圖 6 785 775 771 770 765 2012 2013 2014 2015 2016 760 765 770 775 780 785 790 廢 污 水 排 放 量 t 這是我國(guó)近五年廢污水排放量的總趨勢(shì) 呈下降狀態(tài) 年平均污水占總比約 7 據(jù)統(tǒng)計(jì) 我國(guó)每年的污水約 80 未經(jīng)任何處理直接排入江河湖庫(kù) 90 以上的城市地表 水體 97 以的城市地下含水層收到污染 由此 水資源水環(huán)境惡化 也是我國(guó)水資源 現(xiàn)存在的問(wèn)題 從這幾方面因素來(lái)看 我國(guó)農(nóng)業(yè)用水量占主要地位 而農(nóng)業(yè)用水效率較低 大部 分都被浪費(fèi)掉 所以我們主要對(duì)農(nóng)田灌溉水有效利用系數(shù)進(jìn)行分析 假設(shè)農(nóng)田灌溉用 水有效利用率系數(shù)函數(shù)為 y aebt 使用 Excel 進(jìn)行參數(shù)擬合得到農(nóng)田灌溉水有效利用系 數(shù) 0 1 2 3 4 5 6 0 51 0 515 0 52 0 525 0 53 0 535 0 54 0 545 農(nóng) 田 灌 溉 有 效 利 用 系 數(shù) 擬 合 經(jīng) Excel 線性擬合所得出的函數(shù)為 y 0 510e0 012t 由此可知 當(dāng)?shù)?2030 年時(shí) t 取 18 y 0 5844 計(jì)算與國(guó)家 2030 年的目標(biāo)偏差 A A 7 得目標(biāo)偏差 A 2 6 綜上所述 我國(guó)現(xiàn)在水資源開(kāi)發(fā)利用現(xiàn)狀及存在的主要問(wèn)題如下 1 水資源供需矛盾日益尖銳 我國(guó)總水資源總量占世界的 6 而人口卻占全球的 23 左右 人均水資源量很低 只有世界平均值的四分之一 甚至曾被聯(lián)合國(guó)可持續(xù)發(fā)展委員會(huì)列為世界人均水資源 最貧乏的國(guó)家之一 并且我國(guó)人口基數(shù)大 隨著社會(huì)發(fā)展人口數(shù)仍將保持上漲趨勢(shì) 面臨著重大挑戰(zhàn) 2 水資源浪費(fèi)嚴(yán)重 我國(guó)雖擁有著較多的水資源 卻也是水資源浪費(fèi)最嚴(yán)重的國(guó)家之一 生產(chǎn)同樣多的 糧食 我國(guó)要比其他發(fā)達(dá)國(guó)家多用一倍的水 我國(guó)萬(wàn)元產(chǎn)值的耗水量為 225m 發(fā)達(dá)國(guó) 家卻僅有 110m 此外 城市生活用水量雖遠(yuǎn)低于農(nóng)業(yè)用水和工業(yè)用水量 但生活中人 們對(duì)水資源的毫不吝嗇和肆無(wú)忌憚的浪費(fèi)卻與前兩者相差無(wú)幾 這反映了人們節(jié)約水 資源的意識(shí)的淡薄 與當(dāng)前我國(guó)水資源現(xiàn)狀形成極大地反差 3 水污染嚴(yán)重 近幾年我國(guó)的污水排放量有所下降 但仍存在很大問(wèn)題 下降并不能說(shuō)明情況好轉(zhuǎn) 由于用廉價(jià)淡水稀釋污水從而達(dá)到排污標(biāo)準(zhǔn)成為眾多排污企業(yè)的慣用伎倆 大量未經(jīng) 處理不達(dá)標(biāo)的污水廢水北直街排入江河湖海庫(kù)等公用水體 甚至一些地區(qū)陷入 越污 染 越缺水 越污染 的惡性循環(huán) 此外 濫用化肥 農(nóng)藥 水土流失對(duì)水資源帶來(lái) 了嚴(yán)重的污染 由于水資源污染日益加劇 部分公共水體的承載能力被打破 出現(xiàn)嚴(yán) 重的水質(zhì)退化 導(dǎo)致水資源利用量進(jìn)一步減少 4 2 問(wèn)題二 4 2 1 符號(hào)說(shuō)明 8 h 用水效率 p 城市一年總產(chǎn)值 w 城市全年總用水 s 萬(wàn)元工業(yè)增加值用水量 u 工業(yè)增加值 f 農(nóng)田灌溉用水利用系數(shù) a1 工業(yè)用水占比 a2 農(nóng)業(yè)用水占比 w2 城市農(nóng)田用水總量 w3 農(nóng)田實(shí)際吸水總量 R2 擬合優(yōu)度 4 2 2 問(wèn)題概述 保證在不超過(guò)水資源開(kāi)發(fā)利用量和用水效率控制紅線的前提下實(shí)現(xiàn)萬(wàn)元工業(yè)增加 用水量呈下降趨勢(shì)以及農(nóng)田灌溉水利用系數(shù)提高到 0 6 以上的目標(biāo) 4 2 3 模型的建立與求解 模型建立 為了分析在農(nóng)田灌溉水有效利用系數(shù)確定和工農(nóng)業(yè)用水占比系數(shù)可預(yù)測(cè)的情況下 各指標(biāo)的變化情況需要構(gòu)建以下變量指標(biāo)的函數(shù)關(guān)系式 因此需要明確各個(gè)指標(biāo)的定 義 1 用水效率 所用單位體積水產(chǎn)生的生產(chǎn)力結(jié)果或收益結(jié)果 在一個(gè)省所產(chǎn)生的效 9 果 直觀表示為每立方米的水所產(chǎn)生的經(jīng)濟(jì)效益 即該省全民生產(chǎn)總值與該省一年內(nèi) 總用水量的比值 記做 h 1 2 萬(wàn)元工業(yè)增加值用水量 工業(yè)用水量 工業(yè)增加值 記做 S 2 3 農(nóng)田灌溉水有效利用系數(shù) 在一次灌水期間被農(nóng)作物利用的凈水量與水源渠首處 總引進(jìn)的水量的比值 在一個(gè)省中即可表示為一年內(nèi)該省農(nóng)田實(shí)際吸收水的總量與一 年內(nèi)該省農(nóng)業(yè)用水總量的比值 記做 f 3 2 3 由 1 2 3 式可得 P 4 1 P 5 3 2 10 其中 a1 a2 分別為表格中所算工業(yè)用水占比和農(nóng)業(yè)用水占比 由 4 5 式 得到用水效率 萬(wàn)元工業(yè)增加值用水量和農(nóng)田灌溉用水有效利用系 數(shù)之間的關(guān)系 即可建立其相應(yīng)關(guān)系模型 模型求解 為確保水資源開(kāi)發(fā)利用量和用水效率不超過(guò)控制紅線的前提下保證萬(wàn)元工業(yè)增值用 水量呈下降趨勢(shì)以及農(nóng)田灌溉水有效利用系數(shù)提高到 0 6 以上 我們以廣東省為例進(jìn) 行分析 收集了歷年廣東省水資源信息數(shù)據(jù) 2008 年至 2016 年的數(shù)據(jù)如下 年 份 w總 用 水 量 億 m 農(nóng) 業(yè) 用 水 億 m 工 業(yè) 用 水 億 m 工 業(yè) 用 水 占 比 農(nóng) 業(yè) 用 水 占 比2008 462 5 224 8 141 1 0 3051 0 4861 2009 461 5 227 7 137 2 0 2973 0 49342010 463 4 228 7 136 2 0 2939 0 4935 2011 469 227 5 138 8 0 2959 0 48512012 469 227 5 138 8 0 2959 0 4851 2013 451 227 6 121 6 0 2696 0 50472014 443 2 223 7 119 6 0 2699 0 5047 2015 442 5 224 3 117 0 2644 0 50692016 433 1 227 112 5 0 2598 0 5241 我們用工業(yè)和農(nóng)業(yè)的用水占比系數(shù)繪制散點(diǎn)圖 并進(jìn)一步分析在農(nóng)田灌溉水有效利 用系數(shù)紅線確定限制的條件下其他指標(biāo)的變化趨勢(shì) 來(lái)分析并說(shuō)明如何調(diào)節(jié)相關(guān)參數(shù) 來(lái)確保其他紅線能夠不被超過(guò) 下面是工業(yè)與農(nóng)業(yè)近幾年的散點(diǎn)圖像 11 0 3051 0 2973 0 2939 0 2959 0 2959 0 2696 0 2699 0 2644 0 2598 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 0 2550 0 2600 0 2650 0 2700 0 2750 0 2800 0 2850 0 2900 0 2950 0 3000 0 3050 0 3100 工 業(yè) 用 水 占 比 工 業(yè) 用 水 占 比 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 0 4800 0 4850 0 4900 0 4950 0 5000 0 5050 0 5100 0 5150 0 5200 0 5250 0 5300 農(nóng) 業(yè) 用 水 占 比 農(nóng) 業(yè) 用 水 占 比 隨后利用 Excel 軟件來(lái)擬合函數(shù) 工業(yè)用水占比函數(shù)式為 12 y 0 0003 x 2 0 0027 x 0 3072 R2 0 8923 農(nóng)業(yè)用水占比函數(shù)式為 y 0 0008 x 2 0 0045 x 0 494 R2 0 8289 擬合的函數(shù)圖像如下 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 0 2500 0 2600 0 2700 0 2800 0 2900 0 3000 0 3100 工 業(yè) 用 水 占 比 Poly 工 業(yè) 用 水 占 比 工 業(yè) 用 水 占 比 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 0 4800 0 4900 0 5000 0 5100 0 5200 0 5300 農(nóng) 業(yè) 用 水 占 比 Poly 農(nóng) 業(yè) 用 水 占 比 農(nóng) 業(yè) 用 水 占 比 13 R2 均大于 0 8 說(shuō)明擬合程度比較高 我們預(yù)測(cè) 2030 年的數(shù)據(jù) 將 t 2030 代入到上述函數(shù)式 可預(yù)測(cè)到 2030 年的時(shí)候 廣東省工業(yè)占比系數(shù) a1 和農(nóng)業(yè)占比系數(shù) a2 分別為 a1 0 1026 a2 0 7822 此時(shí) 我們讓 2030 年農(nóng)業(yè)灌溉水有效有效利用系數(shù) f 0 6 代入到 5 式中 進(jìn)而 確定工業(yè)總產(chǎn)值 P 172250 00 由 1 5 式我們可以得到如下工農(nóng)業(yè)占比系數(shù)在 農(nóng)業(yè)灌溉水有效利用系數(shù) f 0 6 條件限制下的函數(shù) P hw sufa2 W3a1 我們知道了了 f 每年的變化趨勢(shì) 為了更直觀表達(dá)農(nóng)業(yè)灌溉水有效利用率 f 確定的情況 下其他變量的變化趨勢(shì) 我們利用 MATLAB 做出相應(yīng)的三維網(wǎng)格圖 經(jīng) MATLAB 繪制次函數(shù)三維圖像如下 14 其中 Z 軸表示全年生產(chǎn)總值 Y 軸表示 一年該省萬(wàn)元工業(yè)增加值 X 軸即我們要分析 的萬(wàn)元工業(yè)增加值用水量的變化趨勢(shì) 我們可以看出 在我們所觀測(cè)的年份范圍內(nèi) 要想使 Z 量呈增大趨勢(shì) 同時(shí) Z 是與總用水量有直接關(guān)系 要確保每年要以約 1 57 的萬(wàn)元工業(yè)增加值用水量才能達(dá)到目標(biāo) 同時(shí)通過(guò)占比分析保證了全國(guó)用水總量不超 過(guò) 7000 億立方米的紅線 4 3 問(wèn)題三 4 3 1 符號(hào)說(shuō)明 Y 萬(wàn)元工業(yè)增加值用水量 15 X1 萬(wàn)元國(guó)內(nèi)生產(chǎn)總值用水量 X2 耕地實(shí)際灌溉畝均用水量 b 待定系數(shù) R 灰色關(guān)聯(lián)度 4 3 2 問(wèn)題概述 分析不同省的萬(wàn)元國(guó)內(nèi)生產(chǎn)總值用水量 耕地實(shí)際灌溉畝均用水量和萬(wàn)元工業(yè)增 加值的內(nèi)在規(guī)律并預(yù)測(cè) 2020 至 2030 年各省水管理控制發(fā)展趨勢(shì) 4 3 3 模型建立與求解 對(duì) 2015 年和 2016 年萬(wàn)元國(guó)內(nèi)生產(chǎn)總值用水量 m3 耕地實(shí)際灌溉畝均用水量 m3 和萬(wàn)元工業(yè)增加值用水量 m3 數(shù)據(jù)建立非線性回歸模型 Y b 1 b 2 X1 X1 b 3 X1 b 4 X2 X2 b 5 X2 b 6 X1 X2 數(shù)據(jù)收集 16 17 分別將 2015 年和 2016 年各省數(shù)據(jù)以 x y z 的形式代入 利用 MATLAB 建立三維散 點(diǎn)圖 并擬合出三維曲面圖 求解系數(shù) b 2015 年 b 5 28 0861 0 0008 0 4900 0 0002 0 1759 0 00015 18 2016 年 b 5 23 5716 0 0007 0 4340 0 0001 0 1738 0 0000 由圖形和回歸方程可得萬(wàn)元國(guó)內(nèi)生產(chǎn)總值用水量 耕地實(shí)際灌溉畝均用水量和萬(wàn)元 工業(yè)增加值用水量的相關(guān)變化趨勢(shì) 隨著萬(wàn)元國(guó)內(nèi)生產(chǎn)總值用水量的增加與耕地實(shí)際 官高畝均用水量的增加 萬(wàn)元工業(yè)增加值用水量普遍也會(huì)增加 現(xiàn)建立灰色關(guān)聯(lián)矩陣模型 對(duì)于一個(gè)參考數(shù)列 X0 有若干個(gè)比較數(shù)列 X1 X2 Xn 各比較數(shù)列與參考數(shù)列在各個(gè)時(shí)刻 即曲線中的各點(diǎn) 的關(guān)聯(lián)系數(shù) Xi 可由 下列公式算出 為分辨系數(shù) 一般在 0 1 之間 通常取 0 5 第二級(jí)最小差 記為 min 兩級(jí)最大差 記為 max 為各比較數(shù)列 Xi 曲線上的每一個(gè)點(diǎn)與參考數(shù)列 X0 曲線上的每一個(gè)點(diǎn)的絕對(duì)差值 記為 oi k 關(guān)聯(lián)系數(shù) Xi 簡(jiǎn)化如下列公式 19 將各個(gè)時(shí)刻 即曲線中的各點(diǎn) 的關(guān)聯(lián)系數(shù)集中為一個(gè)值 即求其平均值 作為比 較數(shù)列與參考數(shù)列間關(guān)聯(lián)程度的數(shù)量表示 關(guān)聯(lián)度 公式如下 利用此公式結(jié)合 MATLAB 編程 計(jì)算出 2015 年的數(shù)據(jù)的灰色關(guān)聯(lián)度 此處用 R 表示 Rx1 0 6252 Rx2 0 7114 2016 年 Rx1 0 6180 Rx2 0 7222 下標(biāo) x1 x2 代表因子 我們發(fā)現(xiàn) R x2均大于 Rx1 表明耕地實(shí)際灌溉畝均用水量 對(duì)萬(wàn)元工業(yè)增加值用水量關(guān)聯(lián)程度大 X 1萬(wàn)元國(guó)內(nèi)生產(chǎn)總值用水量灰色關(guān)聯(lián)度與 X2表 示耕地實(shí)際灌溉畝均用水量灰色關(guān)聯(lián)度相差不大 但實(shí)際上來(lái)說(shuō) 我國(guó)終是農(nóng)業(yè)大國(guó) 隨著今后的不斷發(fā)展 農(nóng)業(yè)的相關(guān)系數(shù)占比只會(huì)越來(lái)越大 并且國(guó)內(nèi)生產(chǎn)總值也與農(nóng) 業(yè)生產(chǎn)有著密不可分的關(guān)系 因此 耕地實(shí)際灌溉畝均用水量是極為重要的因子 由 此預(yù)測(cè) 2020 年至 2030 年各省的水資源管理及控制的發(fā)展趨勢(shì) 將以耕地實(shí)際灌溉畝 均用水量提高為主 這樣才能最大化的科學(xué)有效的提升國(guó)內(nèi)生產(chǎn)總值等一系列指標(biāo) 4 4 問(wèn)題四 致國(guó)家水利部門一封信 尊敬的各位領(lǐng)導(dǎo) 您們好 鑒于地球水資源極為有限 但需求量卻每日劇增這一情況 根據(jù)我們對(duì)水資源的研 究 隨著社會(huì)發(fā)展 工業(yè)水平不斷進(jìn)步 工業(yè)全國(guó)工業(yè)用水呈下降趨勢(shì) 每立方米工 業(yè)用水產(chǎn)生的經(jīng)濟(jì)效益不斷增加 我國(guó)還處于經(jīng)濟(jì)發(fā)展階段 農(nóng)業(yè)用水占大部分 且 20 農(nóng)田有效用水利用系數(shù)較低 和世界上發(fā)達(dá)國(guó)家相比較有較大差距 我國(guó)水資源豐富 但是人均水資源占有量較低 用水效率有待提高 若不加以調(diào)整水資源利用的方式 制定相關(guān)計(jì)劃 科學(xué)合理的利用水資源 中國(guó)總用水將按照每年 1 01 的增長(zhǎng)率持續(xù)增 長(zhǎng) 到 2030 年 我國(guó)全國(guó)用水總量將達(dá)到 7335 4 億立方米 超過(guò)國(guó)家計(jì)劃用水量控 制紅線 7000 億立方米 有效處理好產(chǎn)值的增加與用水量的減小的矛盾 企業(yè)用水支出 廢水處理支出與企 業(yè)利潤(rùn)之間的矛盾 提高人民群眾對(duì)有效用水的積極性 合理運(yùn)用水資源這一行動(dòng)刻不容緩 在科技方面 我們應(yīng)當(dāng)加以完善 通過(guò)完善工 業(yè)設(shè)備 提高我國(guó)先進(jìn)的工業(yè)節(jié)約用水技術(shù) 高效利用 合理地開(kāi)發(fā)和利用現(xiàn)有的水 資源 并且 對(duì)于工業(yè)廢水的排放處理 應(yīng)當(dāng)嚴(yán)格達(dá)到標(biāo)準(zhǔn)后再進(jìn)行排放 使用再生 水 另外 在完善科技的同時(shí) 還需要通過(guò)教育來(lái)加強(qiáng)人們的節(jié)約用水的意識(shí) 使其內(nèi)化 于心 外化于行 因?yàn)檫@是全民意識(shí)的良好體現(xiàn) 生活上尤其要做到一水多用 因?yàn)?這是節(jié)約用水最行之有效的辦法之一 我們應(yīng)當(dāng)有效地達(dá)到與自然和諧的目的 真正實(shí)現(xiàn)自然之水的循環(huán)利用 全面貫徹 落實(shí)可持續(xù)發(fā)展戰(zhàn)略 目前 生態(tài)文明建設(shè)已經(jīng)被國(guó)家給予了高度重視 我們不能一 味地向自然索取 應(yīng)踐行 十二理念 中的天人和諧 對(duì)于水資源的合理利用 是我 們每一個(gè)人義不容辭的責(zé)任和使命 讓我們齊心協(xié)力 共同實(shí)現(xiàn)這一目標(biāo) 絕不能讓 地球最后一滴水成為我們的眼淚 謝謝各位領(lǐng)導(dǎo) 5 模型改進(jìn) 評(píng)價(jià)與推廣 5 1 模型改進(jìn) 利用 K means 函數(shù)以三個(gè)目標(biāo)用水量 即萬(wàn)元國(guó)內(nèi)生產(chǎn)總值用水量 耕地實(shí)際灌溉 畝均用水量 萬(wàn)元工業(yè)增加值用水量為影響因子 將全國(guó) 31 個(gè)省份按各省相對(duì)于 K means 中聚點(diǎn)相對(duì)位置分層 并挑選出各層次最有代表性的省份 分析數(shù)據(jù) 得到差分 序列圖 進(jìn)而分析得出結(jié)果 5 2 模型評(píng)價(jià) 本模型多用預(yù)測(cè)模型來(lái)解決問(wèn)題 可通過(guò)舊數(shù)據(jù)來(lái)預(yù)測(cè)新數(shù)據(jù) 來(lái)擬合并了解未來(lái) 的發(fā)展趨勢(shì) 可以較好的為相應(yīng)問(wèn)題提供信息并制定相關(guān)戰(zhàn)略 模型在使用的時(shí)候應(yīng) 盡量收集多的數(shù)據(jù) 這樣才能更好地?cái)M合函數(shù)曲線和涵蓋函數(shù)的真實(shí)情況 但對(duì)于不 好收集數(shù)據(jù)的問(wèn)題 比如說(shuō)政治 人文或是真實(shí)技術(shù)的方面 此模型并不能提供較好 的解決方法 21 總體來(lái)說(shuō) 預(yù)測(cè)總歸是預(yù)測(cè) 不能 100 的預(yù)言未來(lái) 我們只能盡可能的去擬合目 標(biāo) 所以 存在一定的誤差也是此模型缺點(diǎn)的一方面 5 3 模型推廣 此模型在多方面領(lǐng)域都有涉足 具有較為穩(wěn)定的性質(zhì) 一般來(lái)說(shuō) 若具有較多數(shù)據(jù) 并且這些數(shù)據(jù)都具有相應(yīng)的關(guān)系 此模型都可以較為準(zhǔn)確地解決問(wèn)題 此模型還可以用在人口增長(zhǎng)預(yù)測(cè) 區(qū)域經(jīng)濟(jì)指標(biāo)預(yù)測(cè) 房地產(chǎn)價(jià)格指數(shù)和基于大數(shù) 據(jù)背景下的問(wèn)題研究等等 通過(guò)建立數(shù)據(jù)集合 建立相關(guān)指標(biāo)的關(guān)系 可進(jìn)一步分析 該問(wèn)題所需要了解并實(shí)行的發(fā)展規(guī)劃和趨勢(shì) 參考文獻(xiàn) 1 楊啟帆 何勇 談之奕 數(shù)學(xué)建模競(jìng)賽 杭州 浙江大學(xué)出版社 2005 5 2 張曉宇 竇世卿 我國(guó)水資源管理現(xiàn)狀及對(duì)策 第 15 卷 3 期 2006 3 中華人民共和國(guó)水利部 水資源公報(bào) 2008 2016 2008 2016 4 吳芳 張新鋒 崔學(xué)鋒 中國(guó)水資源利用特征及未來(lái)發(fā)展趨勢(shì)分析 2017 34 1 30 39 5 姜兵 中國(guó)水資源利用效率趨于差異及影響因素研究 第 8 卷 8 期 2017 4 6 云舟工作室 MATLAB 數(shù)學(xué)建?；A(chǔ)教程 M 人民郵電出版社 2001 7 尚月強(qiáng) 楊一都 Matlab 及其在數(shù)學(xué)建模中的應(yīng)用 J 貴州師范大學(xué)學(xué)報(bào) 自然 科學(xué)版 2005 23 1 77 81 8 郎連和 區(qū)域水資源可持續(xù)利用評(píng)價(jià)指標(biāo)與評(píng)價(jià)方法研究 自然資報(bào) 2014 8 15 9 姜啟源等 數(shù)學(xué)模型 第三版 北京 高等教育出版社 2003 10 周凱 宋軍全等 數(shù)學(xué)建模競(jìng)賽入門與提高 杭州 浙江大學(xué)出版社 2011 10 22 附件一 灰色模型預(yù)測(cè) Matlab 程序 clc clear syms a b c a b A 69 67 59 5 58 3 52 8 B cumsum A 原始數(shù)據(jù)累加 n length A for i 1 n 1 C i B i B i 1 2 生成累加矩陣 end 計(jì)算待定參數(shù)的值 D A D 1 D D E C ones 1 n 1 c inv E E E D c c a c 1 b c 2 預(yù)測(cè)后續(xù)數(shù)據(jù) F F 1 A 1 for i 2 n 14 只推測(cè)后 14 個(gè)數(shù)據(jù) 可以從此修改 F i A 1 b a exp a i 1 b a end G G 1 A 1 for i 2 n 14 只推測(cè)后 14 個(gè)數(shù)據(jù) 可以從此修改 23 G i F i F i 1 得到預(yù)測(cè)出來(lái)的數(shù)據(jù) end t1 2012 2016 t2 2012 2030 多 14 組數(shù)據(jù) G h plot t1 A o t2 G 原始數(shù)據(jù)與預(yù)測(cè)數(shù)據(jù)的比較 set h LineWidth 1 5 G 1 至 12 列 69 0000 66 1349 61 4069 57 0170 52 9409 49 1562 45 6420 42 3791 39 3494 36 5364 33 9244 31 4992 13 至 19 列 29 2473 27 1564 25 2150 23 4124 21 7387 20 1846 18 7416 附件二 灰色關(guān)聯(lián)度 Matlab 程序 x 90 17 16 63 57 103 49 94 236 42 82 43 131 77 147 34 60 102 114 61 178 124 50 88 93 109 300 50 175 111 242 619 394 227 218 213 186 327 389 351 442 431 427 355 282 617 547 177 165 430 517 753 873 995 339 406 384 393 612 282 497 565 753 617 58 3 10 5 7 7 17 8 31 3 23 6 18 4 35 9 58 8 90 9 85 4 30 0 96 8 66 0 88 2 11 4 32 6 80 9 81 3 37 3 87 5 66 7 58 5 45 8 76 9 58 5 206 6 18 6 65 1 32 0 44 4 43 9 24 s size x len s 2 num s 1 for i 1 num x i x i x i 1 end dx num len 1 0 for i 1 num for j 1 len 1 dx i j x i j 1 x i j end end r 1 1 len 1 1 for i 2 num for k 1 len 1 r i k 1 1 abs dx 1 k dx i k end end r1 sum r 2 num 2 len 1 s 25 3 32 r1 0 6252 0 7114 x 81 16 15 57 58 102 61 89 229 38 76 39 121 66 134 32 57 87 106 55 159 111 44 82 85 101 271 47 165 103 206 588 380 225 236 217 188 326 358 325 404 439 410 347 280 553 551 173 166 320 513 748 840 990 332 394 389 377 587 287 487 565 688 608 52 8 9 9 7 6 16 6 32 0 22 4 29 2 33 6 56 0 90 2 83 7 26 9 94 0 59 5 81 6 11 5 29 9 74 6 79 6 34 2 73 7 65 5 50 9 48 3 69 7 52 7 164 3 18 3 64 5 28 4 42 1 44 4 s size x len s 2 num s 1 for i 1 num x i x i x i 1 end dx num len 1 0 for i 1 num for j 1 len 1 dx i j x i j 1 x i j end 26 end r 1 1 len 1 1 for i 2 num for k 1 len 1 r i k 1 1 abs dx 1 k dx i k end end r1 sum r 2 num 2 len 1 s 3 32 r1 0 6180 0 7222 附件三 非線性回歸方程 Matlab 程序 clear clc data 90 17 16 63 57 103 49 94 236 42 82 43 131 77 147 34 60 102 114 61 178 124 50 88 93 109 300 50 175 111 242 619 394 227 218 213 186 327 389 351 442 431 427 355 282 617 547 177 165 430 517 753 873 995 339 406 384 393 612 282 497 565 753 617 27 58 3 10 5 7 7 17 8 31 3 23 6 18 4 35 9 58 8 90 9 85 4 30 0 96 8 66 0 88 2 11 4 32 6 80 9 81 3 37 3 87 5 66 7 58 5 45 8 76 9 58 5 206 6 18 6 65 1 32 0 44 4 43 9 x1 data 1 x2 data 2 y data 3 X ones size x1 x1 x1 x1 x2 x2 x2 x1 x2 b bint r rint stats regress y X scatter3 x1 x2 y filled hold on x1fit min x1 0 5 max x1 x2fit min x2 0 5 max x2 X1FIT X2FIT meshgrid x1fit x2fit YFIT b 1 b 2 X1FIT X1FIT b 3 X1FIT b 4 X2FIT X2FIT b 5 X2FIT b 6 X1FIT X2FIT mesh X1FIT X2FIT YFIT xlabel x1 ylabel x2 zlabel Y view 140 30 b 28 0861 0 0008 28 0 4900 0 0002 0 1759 0 0001 bint 91 5312 35 3590 0 0014 0 0001 0 2534 1 2335 0 0004 0 0001 0 1025 0 4542 0 0013 0 0014 r 0 4908 1 9169 3 2500 13 3480 5 8312 34 0032 22 4542 20 4043 40 0697 29 50 8412 28 7424 5 9395 34 1762 6 8611 4 9844 3 0558 8 5618 16 6230 9 3367 11 4702 4 6839 12 4885 20 6742 11 9880 18 6204 6 0701 95 3076 14 3564 24 2742 39 8736 58 7484 7 0007 rint 30 63 4205 64 4022 61 4065 57 5728 62 3486 55 8485 74 7076 48 0116 54 3658 66 0281 95 5291 27 5228 84 2480 39 3397 83 6959 42 8873 91 9602 11 8207 6 4934 108 1758 33 6700 91 1548 68 4176 56 5386 21 3410 89 6933 53 6216 67 3438 57 7698 67 7386 62 5173 56 4058 49 6934 66 8170 46 7019 79 9479 53 6774 72 3509 64 1926 41 2521 48 8265 58 1944 23 1044 48 0814 41 9803 83 3287 75 6242 51 6482 44 8507 82 0915 69 6353 57 4951 55 8451 134 7700 31 77 4349 48 7221 85 3185 36 7701 100 4235 20 6764 101 1035 16 3934 13 9492 0 0523 stats 0 4472 4 2064 0 0062 972 1002