数学建模之长江水质监测问题

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　长江水质监测

　　摘要

　　本文解决的是长江水质的评价与监测问题，通过分析过去十年不同监测站收集到的长江水质数据，运用不同的理论建立不同的模型，对长江过去十年的水质情况作出评价，然后再预测未来十年长江水质的变化情况。

　　针对问题一：考虑到问题一中需要对长江水质情况作出定量的评价，并分析各地区水质的污染状况，为此，建立模糊综合评价模型确定了其隶属度函数，建立评判因子的权重矩阵，求得最终结果为：水质最差的地方是江西南昌滁槎（ 15号），其次水质差的地方为四川乐山岷江大桥（8号）、湖南长沙新港（12号）以及四川泸州沱江二桥（10号），此四处水质污染严重；水质最好的地方是湖北丹江口胡家岭 （11号）。

　　针对问题二：根据长江的降解系数，可得到污染物随时间的变化量。由于污染源的污染物排放量等于本地区污染物的流量与上游流下的污染物流量之差。因此，建立污染物流量随时间变化的微分方程模型。最后求得：高锰酸钾指数和氨氮的污染源主要集中在宜昌至岳阳之间。

　　针对问题三：根据已知的过去10年的主要统计数据，建立了灰色预测模型。 在相对误差较小的情况下对未来10年的水质情况作出了预测，分析得出结论：未来10年可饮用水所占的比例越来越低，排污量有明显的上升趋势。

　　针对问题四：在问题四中建立多元线性回归方程，利用最小二乘法求解系数，在满足问题四要求的前提下，求出未来10年的允许最大相对排污量，继而求得未来10年每年的相应排污量，后者与前者的差值与未来10年的长江水总流量的乘积，求得最终结果如下表：

　　针对问题五：分析总结前几个问题的结果，找出水质污染的根本原因。结合考察团的调查结果，给出合理的建议和意见。

　　最后，对模型中运用的方法进行了优、缺点评价，在模型的推广中提出了可以建立类似模型解决生活中的一类问题。

　　关键词：模糊评价 微分方程 灰色预测 线性回归

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　1问题重述

　　1.1问题背景

　　长江是我国第一、世界第三大河流，其水质的污染程度日趋严重，已引起了相关政府部门和专家们的高度重视。2004年10月，由全国政协与中国发展研究院联合组成“保护长江万里行”考察团，从长江上游宜宾到下游上海，对沿线21个重点城市做了实地考察，揭示了一幅长江污染的真实画面，其污染程度让人触目惊心。为此，专家们提出“若不及时拯救，长江生态10年内将濒临崩溃”。

　　长江水是许多人赖以生存和发展的资源，保护水资源就是保护我们自己，就长江近年来的水质情况，采取合理的保护和治理措施刻不容缓。

　　1.2解决问题

　　为了制定出合理的治理长江水质污染的方案，根据长江地区近两年多主要水质指标的检测数据，现讨论以下几个问题：

　　(1) 对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的 污染状况。

　　(2) 研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源 的主要分布地区。

　　(3) 假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江 未来水质污染的发展趋势做出预测分析。

　　(4) 根据预测分析，如果未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的 比例控制在20%以内，且没有劣Ⅴ类水,那么每年需要处理多少污水？

　　(5) 发表对解决长江水质污染问题切实可行的建议和意见。

　　2问题分析

　　问题需要对长江水质作出评价和预测，根据附件中已知的数据，建立相应的评价预测模型，分析得出长江水质过去十年以及未来十年的污染状况。

　　2.1问题一的分析

　　问题提供的数据主要反映了三方面的内容即：水质划分等级标准、17个测站点在28个月份中主要污染物的浓度和水质等级。而观测点水质等级的确定是由主要污染物的浓度决定的，因此为了对长江及各观测点作出定量的综合评价我们需要将主要污染物的浓度归一到一个比较量中。这个比较量要能够快速而准确的反映出各测站点的水质污染状况，因此我们评价整个长江流域的水质状况采用模糊综合评价模型。

　　虽然PH值是影响水质等级的因素之一，但是给出的所有PH值都在6~9这个范围内，即PH值对于水质等级几乎没有影响，所以在以后的问题讨论中不再考虑PH值的影响。

　　2.2问题二的分析

　　因为污染物在时间和空间上是动态变化的，为简便起见，不考虑支流的因素，只从纵向的角度，结合本地与上游污水进行水质的污染分析。根据常识可知，污

　　1

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　染物的污染源就是新增污染物较多的地方。因此，考虑到降解系数，求得污染物 随时间的变化量，通过相邻两主要观测点之间每千米的相对排污量判定污染源，即本地新增污染物的量等于本地区现有的污染物减去上游通过降解后到达该处的污染物的总量，再比上两地之间的距离，得到相对排污量。建立污染物的量随时间变化的微分方程模型，最后通过新增污染物的相对量大小关系得出主要的污染源地区。

　　2.3问题三的分析测

　　第三问根据过去10年长江的总体水污染状况的监测数据，可以看出长江总体水流量变化不大，但年排污的总量在增加，这使得污染河段比例增加，污染的严重程度呈现快速增长的趋势，即每年污染情况主要与当年的排污量和总水流量等因素有关。为此，首先利用回归分析方法确定出可饮用水的比例与总排污量和总水流量的关系式，然后根据过去10 年排污量，利用灰色预测方法对未来的年排污量做出预测，最后根据总排污量的增长趋势来推断出可饮用水比例的变化趋势，从而可以预测出未来10 年长江水质的变化情况。

　　2.4问题四的分析

　　在问题四中建立多元线性回归方程，利用最小二乘法求解系数，在满足未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水的前提下，求出未来10年的允许最大相对排污量，继而求得未来10年每年的相应排污量，后者与前者的差值与未来10年的长江水总流量的乘积，即可求得未来10年要处理的排污量。

　　2.5问题五的分析

　　分析总结前几个问题的结果，找出水质污染的原因。结合考察团的调查结果，给出合理的建议和意见。

　　3模型假设与符号说明

　　3.1模型假设

　　假设一： 降解系数在一定时间段固定不变；

　　假设二：长江干流的自然净化能了近似均匀；

　　假设三：干流的相邻两观测站点之间的排污口主要集中在下游处；

　　假设四：干流污染物的富集主要受上游影响，支流的影响忽略不计。

　　3.2符号说明

　　2

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　4模型的建立与求解

　　4.1问题一的解答

　　4.1.1模型一的建立与求解

　　（1）确立评判指标：根据《地表水环境质量标准》（GB3838—2002）中4个主要项目标准限值设置表，结合28个月的各观测点的测量结果，选取溶氧量（DO）、高锰酸盐指数(CODMn)和氨氮（NH3-N）为评价因子，设置评价因子集：

　　u??x1,x2,x3?，求出17个观测点28个月的个评价因子的均值，得到17个点的评价因子集。

　　（2）建立评价集

　　3

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　4个指标将水质分为6个等级。

　　（3）建立隶属函数，进行单因素评价：

　　由于水质污染程度和水质分级标准都是模糊的，所以用隶属度来描述分级界限较为合理，先根据各指标的6级标准，做出6个级别的隶属函数，其中DO的评价指标以数值大为优，其余两个指标以数值小为优，函数如下：

　　DO（A）以升半梯形分布建立隶属函数

　　?1?x?6?

　　Au1?x???

　　?7.5?6??0??0?

　　Au3?x???x?5

　　?x?3??5?3??0?

　　Au5?x???x?2

　　?x??2

　　?x?7.5??6?x?7.5??x?6??x?6orx?3??5?x?6??3?x?5??x?3orx?0??2?x?3??0?x?2?

　　?x?7.5orx?5??0

　　?x?6?

　　?6?x?7.5?Au2?x???

　　?7.5?6??x?5?5?x?6?

　　?0?x?3?

　　Au4?x???

　　?5?3??x?2?0?x?

　　Au6?x???

　　?2??1

　　?x?5orx?2??3?x?5??2?x?3??x?2??0?x?2??x?0?

　　COOMn(B)和NH3-N（C）以升半梯形分布建立隶属函数：

　　4

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　?1?x?2?

　　Bu1?x???

　　?4?2??0??0?

　　?x?4

　　Bu3?x???

　　?6?4?x?6??10?6?0

　　?x?10?

　　Bu5?x???

　　?15?10??1

　　?x?2??2?x?4??x?4??x?10orx?4??4?x?6??6?x?10??x?0??10?x?15??x?15?

　　?

　　?0?

　　?x?2

　　Bu2?x???

　　?4?2?x?4??6?4??0?

　　?x?3

　　Bu4?x???

　　?10?6?x?10??15?10Bu6?x??1

　　?x?2orx?6??2?x?4??4?x?6??x?15orx?6??6?x?10??10?x?15?

　　?x???

　　?x?0.15??1

　　?x?0.15?

　　?0.15?x?0.5?Cu1?x???

　　0.5?0.15?

　　?x?0.5???0

　　?

　　?0?

　　?x?0.5

　　Cu3?x???

　　?1?0.5?x?1??1.5?1?0

　　?x?1.5?

　　Cu5?x???

　　?2?1.5??1

　　?

　　?0?x?0.15orx?1??

　　?x?0.15

　　?0.15?x?0.5?Cu2?x???

　　0.5?0.15?

　　?x?0.5

　　?0.5?x?1???1?0.5?

　　?0?

　　?x?1

　　Cu4?x???

　　?1.5?1?x?1.5??2?1.5Cu6?x??1

　　?x?1.5orx?0.5??0.5?x?1??1?x?1.5??x?1.5??1.5?x?2??x?2?

　　?x?2orx?1??1?x?1.5??1.5?x?2?

　　?x???

　　有了各指标的隶属函数，就可以进行单因素评价，将各观测点的3个指标带入相应的隶属函数，计算出指标的隶属度，得到各观测点模糊关系矩阵Rj（j=1，2，3...，17），具体数据见附录一。

　　5

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　（4）建立评价因素的权重集：

　　由于DO,COOMn和NH3-N等污染指标对水质的影响不同，因此对各指标应赋予不同的权重，根据污染物对水质的污染大、权重的的原则来决定权重的大小。

　　定义：对于DO越大越优型：ai?

　　xi

　　，对于COOMn和NH3-N越小越优型：mi

　　ai?

　　ni

　　，其中，ai为第i种指标的权重，xi是第i种指标的实时均值，mi,ni分xi

　　ai

　　别为多级浓度标准的最大值和最小值。再将权重归一化处理，Qi?

　　?a

　　i?1

　　3

　　，评价

　　i

　　因素权重的集合为A=[0.145,0.226,0.629]。

　　（5）模糊综合评价：

　　评价因素的权重集A乘以单因子矩阵Rj得到模糊综合评价结果，

　　Bj?ARj。将6个等级看作一种相对位置使其连续化，设各等级分别用1，2，3，4，5，6表示，并称为各等级的秩，然后用B中对应的量将各级的秩加权求积，得到被评等级的相对位置，即：

　　B*?

　　?b

　　k?16k?1

　　6

　　k

　　?k

　　k

　　?b

　　式中bk为隶属于第k级的隶属度。具体计算结果见下表：

　　4.1.2问题一的结果分析

　　由表1可知，各观测站中水质等级最大的为江西南昌滁槎 （15号），可知此处水质情况最差，其次水质差的地方为四川乐山岷江大桥（8号）、湖南长沙新港（12号）以及四川泸州沱江二桥（10号），此四处水质污染严重；水质最

　　6

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　好的地方是湖北丹江口胡家岭 （11号）。

　　4.2问题二的解答

　　因为污染物在时间和空间上是动态变化的，为简便起见，不考虑支流的因素，只从纵向的角度，结合本地与上游污水进行水质的污染分析。建立污染物流量随时间变化的微分方程模型，最后通过新增污染物的流量大小关系得出主要的污染源地区。

　　4.2.1模型二的建立与求解

　　（1）某一污染物扩散所满足的微分方程是一个抛物线方程，结合实际问题的假设，常假设其水流近似的处于稳定状态，断面均匀，普通对流扩散方程为

　　?C?C?2C?u?D2?kC ?t?x?x

　　上式是解决污染物的一般方程，结合问题二的情况，假设观测站干流的污染

　　?C?0，又因为D?0,所以上式进一步简化为：物浓度在一个月内保持不变，则 ?t

　　?Cu??kC ?x

　　其中，u为断面水流平均流速，C为某组分子在x处的浓度，k为降解常数取0.2。由此解出长江干流污染物浓度Cx与据观测站x米处的函数关系：

　　Cx?Cie?kxui

　　（2）对于长江上任一段干流AB，A为第i个观测点，B为第i?1个观测点，两观测点的距离为dii?1，假设该段干流之间有n个排污口（包括支流入口和直排口），为计算的方便，假设排污口聚集在一起，其总的排污口的流量为q，平均 流量为u，污染物的浓度为c，qc即为AB段的排污量,则

　　?kxi?1

　　ui?1

　　Ci?qc?Qi?1Ci?1e

　　Qi

　　?kdAB

　　U CB?qc?QACAe

　　QB

　　?kdAB

　　UW?qc?CBQB?CAQAe

　　U?UA?UB 2

　　对于任意一个江段AB，污染物的浓度CA,CB，水流量QA,QB，流速UA,UB，距

　　7

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　离dAB和降解系数k都已知，则可以算出AB段污染物的排放量Wj。

　　（3）根据所给数据可以算出，对于每个月每个江段的排污区间Wj

　　113

　　（j=1,2,...,13个月），取13个月的平均值W??Wj,表示平均排污量，则

　　13j?1该江段每千米的排污量w?

　　W

　　为相对排污量，每个江段都有一个相对排污量，dAB

　　以此为指标可以判断长江干流的主要污染源。

　　4.2.2模型二的求解

　　按步骤运用MATLAB编求解（程序见附录3），结果如下表：

　　4.2.3问题二的结果分析

　　问题二用长江干流7个观测站点将长江分为6个江段，通过污染物流量随时

　　间的变化关系建立了微分方程模型，先计算出每月每段的高锰酸盐和氨氮的量，再对六个污染源近一年所排放的污染物的数量求期望，分析得出污染物高锰酸盐指数和氨氮的污染源主要集中在宜昌至岳阳之间。

　　4.3问题三的解答

　　4.3.1模型三的建立与求解

　　当系统内部信息和特性是部分已知的，另一部分是未知时，人们往往难以建立客观的物理原形，内部因素难以辨识，以及相互之间的关系较为隐蔽，难已准确了解这类系统的行为特征，因此，对于这类问题进行定量描述，即建立模型难度较大，便选择建立灰色系统，对问题求解。

　　1. 灰色动态模型

　　8

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　附件4中给出了1995-2004年长江各等级在丰水期，枯水期，水文年的百分比，为了预测此趋势之下未来十年长江水质的变化，通过对水文年全流域水质的变化研究来对长江未来水质污染的发展趋势作出预测分析。建立灰色预测模型，首先舍去异常数据：1995年IV类、劣V类数据，然后从附件4中提取出水文年全流域随时间变化的原始数据。

　　x?0??x?0??1?,x?0??2?....x?0??n?n?10

　　对于原始数据进行累加处理，即

　　x?k???x?0??i??k?1,2....10? 1

　　i?1k??

　　加生成数列为

　　???????x?1??x?1??1?,x?1??2?....x?1??n?n?10?

　　采用一阶单变量微分方程进行拟合，得到白化方程的GM（1,1）模型

　　dx?1?

　　?ax?1??t??u ? ?dt??

　　其中式中的a,u为待定系数。

　　综上所述，灰色动态模型为：

　　?1??k???0.5x?1??k??0.5x?1??k?1??z ??0??1???x?k??az?k??u

　　???????

　　2. 模型求解

　　取水文年全流域的饮用水预测为例，可通过对其变化研究对长江未来水质污 染的发展趋势作出预测分析。

　　(1) 求解步骤：

　　第一步：从附表中提取出数据

　　x?0???93.1,85.3,80.7,88.4,80.2,74,73.7,76.7,77.5,68?

　　第二步：利用公式对x?0?做累加生成得

　　x?1???93.1178.4259.1347.5427.7501.7575.4652.1729.6797.6?

　　?其中 第三步：构造矩阵B和数据向量，x?0?与x?1?满足关系Y?n??ba

　　Yn?x?0??2?...x?0??n???Tn?10

　　T??85.380.788.480.27473.776.777.568?

　　9

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　?1??2?1B???

　　?2?1????2

　　???

　　??135.751?

　　??218.751?????303.31????1??1??x?1??x?1?,1???387.61???

　　?464.71??1??1??x?2??x?3?,1??? ???538.551?????1??1?

　　x?n?1??x?n?,1???613.751?

　　????690.851???690.851??????763.61??

　　?

　　?

　　?

　　第四步：求解常系数?

　　?a?

　　?a????BTB?u?

　　??

　　?1

　　?0.0231?

　　?????????????? BYn???

　　?88.8436?

　　T

　　第五步：得出表达式：

　　?1u??a?ku

　　Y?1????3753?e?0.0231k?3846 ?x?1????ek?1?

　　a?a?

　　?? 令

　　11

　　??? Y?0?Y?Yk??k?1??k?

　　所以Y?0k?表示还原后的值。

　　(2) 结果求解

　　按步骤运用MATLAB编求解（程序见附录4），运用此程序根据过去十年的 数据可以预测出未来10年的排污量和可饮用水如下表：

　　表4：未来10年排污量

　　表5：未来10年可饮用水所占比例

　　4.3.2问题三的结果分析

　　从表5可以看出，可饮用水所占的比例越来越低，所以，可以看出水质越来越差，从2012年起，可饮用水的百分比甚至低于60%,如此严重的水污染，应该引起足够重视，采用合理的措施来减少对长江水的污染。

　　10

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　4.4问题四的解答

　　4.4.1模型的建立与求解

　　（1）长江干流的Ⅳ类和Ⅴ类水以及劣V类水的比例与每年长江总水流量x1?t?和总排污量x2?t?有关，因此分别以Ⅳ类和Ⅴ类水所占的比例y1和劣V类水所占比例y2为响应变量,以历年来的

　　? ,9513,9171.26,9513,9924,8892.8 ,10210,9980,9405x1?t?=?9205

　　x2?t?=?174,179,183,207,234,220.5,256,270,285?

　　为解释变量，做出多元线性回归，得到Ⅳ类和Ⅴ类水的总和及劣V类水与长江总水流量及总排污量的关系，做出一般的多元线性回归模型：

　　y1?a1x1?t??b1x2?t??c1 y2?a2x1?t??b2x2?t??c2

　　运用matlab用最小二乘法求解得到回归系数

　　?a1,b1,c1????0.0032,0.079,15.9276?

　　?a2,b2,c2????0.0077,0.1485,56.2805?

　　就求得了Ⅳ类和Ⅴ类水的总和及劣V类水与长江总水流量及总排污量的关系。

　　问题四要解决问题是：未来10年内每年都要求长江干流的Ⅳ类和Ⅴ类水的比例控制在20%以内，且没有劣Ⅴ类水。则约束条件为：

　　?y1??0.0077x1?t??0.1485x2?t??56.2805?20??y2??0.0032x1?t??0.079x2?t??15.9276?0

　　综上所述，所建立的模型为：

　　?y1??0.0032x1?t??0.079x2?t??15.9276??y2??0.0077x1?t??0.1845x2?t??c

　　?y??0.0077x1?t??0.1485x2?t??56.2805?20s.t.?1

　　?y2??0.0032x1?t??0.079x2?t??15.9276?0

　　对上述回归方程进行回归检验时，发现回归效果并不是太好，考虑到长江干流的Ⅳ类和Ⅴ类水的比例与江总水流量及总排污量不是直接的线性关系，因此对模型进行改进。

　　（2）优化模型

　　长江干流的Ⅳ类和Ⅴ类水以及劣V类水的比例与每年x2?t?总排污量和长江x1?t? 11

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　x?t?总水流量的比值Z?t??2相对排污量有关，因此分别以Ⅳ类和Ⅴ类水所占的比x1t例y1和劣V类水所占比例y2为解释变量，Z?t?为响应变量，做多元线性回归。对表四的数据进行观测，得到Z?t?与时间t成线性关系，因此列出如下线性回归方程：

　　?Z1?t???t????Z2?t??ay1?t??by2?t??c

　　运用matlab用最小二乘法进行求解，得到回归系数分别为

　　?? ? ,????0.0012,0.0165

　　?a,b,c???0.0002,0.0005,0.0185?

　　计算时剔除了2003年这个突变点，且检验到回归性较好。得到的回归方程如下：

　　t?0.0165?Z1?t??0.0012 ???????Zt?0.0002yt?0.0005yt?0.016512?2

　　将y1?t??20,y2?t??0带入Z2?t?的表达式中，得到最大允许相对排污量

　　将t?11,12,13,14,...,20带入Z1?t?的表达式中得到未来10年每年预Z2?t??0.0225，

　　测的相对排污量，则污水处理量w?t??max?Z1?t??Z2?t?,0??L，因为每年的总水流量（除特殊年份如1998年）变化不大，故取前10年的总水流量均值为每年的110

　　总水流量L??x?t?，得到每年的污水处理量如下表：

　　10t?1

　　4.4.2问题四的结果分析

　　从表6可以看出，未来10年所需处理的污水越来越多，所以，可以看出水质越来越差，如此严重的水污染，应该引起足够重视，采用合理的措施来减少对长江水的污染。

　　4.5解决长江水质问题的建议

　　通过问题一的计算可以看出，长江水质的枯水期、丰水期与平水期都由水流量的不同导致水质情况也不同。从问题三、四的求解结果可以看出，在未来的十

　　12

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　年里，长江水质不断恶化，废水排放量逐年递增，需要处理的污水也越来越多。

　　根据模型假设，在不发生特大洪涝灾害时长江的水流量是一定的，总结得出解决长江水质问题的关键，在于减少废水的排放量的同时增强废水的处理能力。鉴于此，我们认为可以从以下几个方面着手解决长江水质问题：

　　（1）加大处罚力度。

　　（2） 对污染严重地区进行重点治理。

　　（3） 加强宣传，增强民众保护长江的意识。

　　（4） 严格进行检测，处理违规企业，严惩腐败。

　　（5） 保护长江立法，引起政府和人民的高度重视。

　　5模型的评价、改进与推广

　　5.1模型的评价

　　5.1.1模型的优点

　　（1）问题一中建立了模糊综合评价模型，通过精确的数字手段处理模糊的评价对象，对蕴藏信息呈现模糊性的资料作出了比较科学、合理、贴近实际的量化评价；

　　（2）问题二的模型，简化了许多因素，忽略了支流的影响，但这种理想化还是合理的，较符合实际情况;

　　（3）问题三采用灰色预测模型，具有少数据性、良好的时效性、较强的系统性和关联性等特征，可以合理的对数据做出预测。

　　5.1.2模型的缺点

　　在求解模型时，为了使问题得到方便的解决，往往采用简化的手段进行求解，因此求出的结果与真实值会存在一定的偏差。

　　5.2模型的改进

　　在问题三中运用了灰色动态模型对未来10年长江的水质作出了预测，但并未对预测的结果进行检验，不知预测结果的准确度有多高，在此，最好能将模型三的基础上对灰色模型的预测结果进行检验，以保证预测结果有较高的可信度。

　　5.3模型的推广

　　模型一中，除了用模糊算法评价外，还可以用BP神经网络来求长江水质的在综合评价而；而在求解问题三和问题四时采用的灰色预测，还可以用回归分析预测、时间序列预测和神经网络预测等代替。

　　6参考文献

　　[1] 宋晓秋编著，模糊数学原理与方法，北京，中国矿业大学出版社，1999

　　[2] 张俊福等编著，应用模糊数学，北京，地质出版社，2001

　　[3] 岳超源，决策理论与方法，北京：科学出版社，2003

　　[4] 邓聚龙，灰理论基础,武汉：华中科技大学出版社，2002

　　[5] 韩中庚，数学建模方法及其应用，北京：高等教育出版社，2005

　　13

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　附录一

　　各观测点模糊关系矩阵：

　　00000?00000??1?1?2号?0.04820.04820000? 0.21610.216100001号????????0.0940.0940000???0.51940.51940000??

　　00000?00000??1?1?4号?0.89290.89290000? 0.43750.437500003号???????0.90000??0.32660.32660000???0.9?

　　00000??1?0.97030.97030000??6号?0.28750.28750000? 0.21430.214300005号????????0.02970.02970000???0.2280.2280000??

　　?0.99410.99410000??00.55860.5586000??8号?00.62150.6215000? 0.04650.046500007号???????00000??1??00.84860.8486000??

　　00000?0000??1?0.57670.5767?10号?0.66970.6697? 0.36790.3679000000009号???????0.62360.6236000??0.80110.80110000???0?

　　0000??100000??0.07330.0733?12号?0.24290.2429? 100000000011号???????0.83280.8328000??100000???0?

　　00000??1?0.94760.94760000??14号?0.66250.66250000? 00.09650.096500013号???????0000??0.67340.6734??0.13570.13570000??

　　0.69820.6982000?00000??0?1?16号?0.87150.87150000? 0.16200.1620000015号???????00010??0??0.38970.38970000??

　　00000??1?0.51070.5107000017号?????0.39170.39170000??

　　附录二

　　问题一的程序（以17号观测站为例）：

　　14

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　clc;clear

　　a=[1 0 0 0 0 0;0.5107 0.5107 0 0 0 0;0.3917 0.3917 0 0 0 0];

　　b=[0.145 0.226 0.629];

　　c=b*a;

　　d=0;

　　e=0;

　　for i=1:6

　　d=d+c(1,i)*i;

　　end

　　for i=1:6

　　e=e+c(1,i);

　　end

　　B=d/e

　　问题二的程序：

　　clc;clear

　　[num11]=xlsread('11.xls');

　　[num29]=xlsread('29.xls');

　　for i=1:6

　　Ca=num11(i,6);

　　Cb=num11(i+1,6);

　　Qa=num29(2,i+2);

　　Qb=num29(2,i+3);

　　d1=num29(1,i+1);

　　d2=num29(1,i+2);

　　d=(d2-d1)*1000;

　　Ua=num29(3,i+2);

　　Ub=num29(3,i+3);

　　U=(86400*(Ua+Ub))/2;

　　m=exp(-0.2*d/U);

　　W=Cb*Qb-Ca*Qa*m;

　　W

　　end

　　问题三的程序：

　　clc;clear

　　y=input('请输入数据 ');

　　n=length(y);

　　yy=ones(n,1);

　　yy(1)=y(1);

　　15 附录三附录四

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　for i=2:n

　　yy(i)=yy(i-1)+y(i);

　　end

　　B=ones(n-1,2);

　　for i=1:(n-1)

　　B(i,1)=-(yy(i)+yy(i+1))/2;

　　B(i,2)=1;

　　end

　　BT=B';

　　for j=1:n-1

　　YN(j)=y(j+1);

　　end

　　YN=YN';

　　A=inv(BT*B)*BT*YN;

　　a=A(1);

　　u=A(2);

　　t=u/a;

　　t_test=input('请输入需要预测个数：');

　　i=1:t_test+n;

　　yys(i+1)=(y(1)-t).*exp(-a.*i)+t;

　　yys(1)=y(1);

　　for j=n+t_test:-1:2

　　ys(j)=yys(j)-yys(j-1);

　　end

　　x=1:n;

　　xs=2:n+t_test;

　　yn=ys(2:n+t_test);

　　plot(x,y,'^r',xs,yn,'*-b');

　　det=0;

　　for i=2:n

　　det=det+abs((yn(i)-y(i))/xs(i));

　　end

　　det=det/(n-1);

　　disp(['百分相对误差为：', num2str(det),'%']);

　　disp(['预测值为： ', num2str(ys(n+1:n+t_test))]);

　　附录五

　　问题四的程序：

　　clc;clear

　　a=[9205 9513 9171.26 13127 9513 9924 8892.8 10210 c=mean2(a);

　　z=0.0002*20+0.0185;

　　b=[11 12 13 14 15 16 17 18 19 20];

　　for i=1:10

　　16 9980 9405];

　　本文解决的是长江水质的评价与监测问题,通过分析过去十年不同监测站收集到的长江水质数据,运用不同的理论建立不同的模型,对长江过去十年的水质情况作出评价,然后再预测未来十年长江水质的变化情况。

　　Z=0.0012*b(1,i)+0.0165; M=(Z-z)*c;

　　M

　　end

　　17