2005年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学) 蔡虎,卢亚军。许超指导教师汪晓银 长江水质评价及预测模型的建立与分析 摘要: 水资源污染日趋严重,本文仅就长江流域若干观测点的最近10年的观测数据进行 建模分析,探讨水资源污染的预测与控制的问题。 对于问题一、我们首先采用模糊综合评判法进行评价,我们以溶解氧、髙锰酸盐指 数、氨氮组成的评价参数集为U={u1,u2,u3},水质等级集为{vl,v2,v3,v4,v5,v6},用隶 属度来刻画长江沿岸17个测点水质在三个参数方面属于不同水质级别的情况,由观测 值和隶属度确定模糊评判矩阵,由最大隶属原则,求出28个时点平均值得到17个测点 的2年多的平均水质情况,其中江西南昌滁槎污染最为严重(测点15),其次是四川乐 山岷江大桥(测点8)。 我们还采用了另外一种一般的综合指数评价方法对水质情况作综合评价。选取Ⅲ等 水的4项指标值PH=7,DO=5,COD=6,NH=1为评价标准。先计算分指数K,然后再把 分指数加权(n)求和得到综合指标爬∑nK。因为17个测点分别有28个时点的监测 数据,我们把这28个月按年分为枯水期、丰水期和平水期。分类标准参照题中附件4 的说明(3)。然后用平均的方法求出按年按期分类的综合指标矩阵。用 Excel分别画出 2003、2004、2005年各测点的水质综合评价柱状图(见附件中图、图、图)。其中 Ⅲ等水的量化指标为0.1,值越低水质越好。最后得出:江西南昌滁槎污染最为严重(测 点15),其次是四川乐山岷江大桥(测点8)。可以看出两种方法评价的结果是一致的。 对于问题二、考虑到江水的流动性和自降解性,测点浓度不能与污染源等同起来, 所以我们把每个监测值看作是由上游对它的影响和本地污染两部分组成的。所以我们考 虑水流速度、测点间距计算出降解次数,进而计算出一个测点对下一个测点的影响度, 通过减去来自上游的影响进而求出每个测点的本地污染度。这样就可以通过比较本地污 染度找出主要污染源了。具体计算出主要污染源有:江西南昌滁槎、四川泸州沱江二桥 四川乐山岷江大桥、湖南长沙新港、湖南岳阳城陵矶、重庆朱沱、湖南岳阳楼等地区 对于问题三、问题要求预测未来十年的水质污染趋势,我们用预测未来10年中六 种水的比例来评价水质污染趋势。具体用灰色系统方法进行预测。又由于每年的总评价 河长不一样,所以不能直接预测,但是比例是不会变的,我们把10年的六种水的比例 用同一个总评价河长(我们假定为10000)转化。再进行预测,得出未来10年的六种水 的河长,然后再转化为比例值。可以看出未来10年长江劣V类水越来多,污染越来越 严重。 对于问题四、我们把6种水重新分为3类:清洁水(I、Ⅱ和Ⅲ)、污染水(Ⅳ和V) 劣质水(V)。用灰色系统方法预测未来10年原6种水的比例,折算出新3类水的比例, 建立废水率与3类水比例的逐步回归模型,确定未来10年每一年废水率与各种水比例 之间的回归方程。然后根据回归方程得到在控制状态下(Ⅳ和Ⅴ类不超过20%、劣Ⅴ类 没有)的废水率,然后与不加控制的预测值进行比较,得到需要处理的废水量(亿吨)。 2005年2006年2007年2008年2009年2010年2011年2012年2013年2014年 79497.51189141.9166.11926231.3261.32948330.1 对于问题五、根据我们所做模型的分析和预测结果,以及目前我国治理长江水质污 染的现状给出:合理制定工业发展规划、污水处理设施的利用、控制排污总量、加大执 法力度等四点建议和意见 关键字:水质评价模糊综合评价综合指数评价灰色系统逐步回归
2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学) 蔡虎，卢亚军。许超 指导教师 汪晓银 1 长江水质评价及预测模型的建立与分析 摘要： 水资源污染日趋严重，本文仅就长江流域若干观测点的最近 10 年的观测数据进行 建模分析，探讨水资源污染的预测与控制的问题。 对于问题一、我们首先采用模糊综合评判法进行评价，我们以溶解氧、高锰酸盐指 数、氨氮组成的评价参数集为 U={u1,u2,u3},水质等级集为{v1,v2,v3,v4,v5,v6}，用隶 属度来刻画长江沿岸 17 个测点水质在三个参数方面属于不同水质级别的情况，由观测 值和隶属度确定模糊评判矩阵，由最大隶属原则，求出 28 个时点平均值得到 17 个测点 的 2 年多的平均水质情况，其中江西南昌滁槎污染最为严重（测点 15），其次是四川乐 山岷江大桥（测点 8）。 我们还采用了另外一种一般的综合指数评价方法对水质情况作综合评价。选取Ⅲ等 水的 4 项指标值 PH0=7，DO0=5，COD0=6，NH0=1 为评价标准。先计算分指数 K，然后再把 分指数加权（η）求和得到综合指标 W=∑ηKj。因为 17 个测点分别有 28 个时点的监测 数据，我们把这 28 个月按年分为枯水期、丰水期和平水期。分类标准参照题中附件 4 的说明（3）。然后用平均的方法求出按年按期分类的综合指标矩阵。用 Excel 分别画出 2003、2004、2005 年各测点的水质综合评价柱状图（见附件中图[1]、图[2]、图[3]）。其中 Ⅲ等水的量化指标为 0.1，值越低水质越好。最后得出：江西南昌滁槎污染最为严重（测 点 15），其次是四川乐山岷江大桥（测点 8）。可以看出两种方法评价的结果是一致的。 对于问题二、考虑到江水的流动性和自降解性，测点浓度不能与污染源等同起来， 所以我们把每个监测值看作是由上游对它的影响和本地污染两部分组成的。所以我们考 虑水流速度、测点间距计算出降解次数，进而计算出一个测点对下一个测点的影响度， 通过减去来自上游的影响进而求出每个测点的本地污染度。这样就可以通过比较本地污 染度找出主要污染源了。具体计算出主要污染源有：江西南昌滁槎、四川泸州沱江二桥 、 四川乐山岷江大桥、湖南长沙新港、湖南岳阳城陵矶、重庆朱沱、湖南岳阳楼等地区。 对于问题三、问题要求预测未来十年的水质污染趋势，我们用预测未来 10 年中六 种水的比例来评价水质污染趋势。具体用灰色系统方法进行预测。又由于每年的总评价 河长不一样，所以不能直接预测，但是比例是不会变的，我们把 10 年的六种水的比例 用同一个总评价河长（我们假定为 10000）转化。再进行预测，得出未来 10 年的六种水 的河长，然后再转化为比例值。可以看出未来 10 年长江劣Ⅴ类水越来多，污染越来越 严重。 对于问题四、我们把 6 种水重新分为 3 类：清洁水（I、Ⅱ和Ⅲ）、污染水（Ⅳ和Ⅴ）、 劣质水（Ⅴ）。用灰色系统方法预测未来 10 年原 6 种水的比例，折算出新 3 类水的比例， 建立废水率与 3 类水比例的逐步回归模型，确定未来 10 年每一年废水率与各种水比例 之间的回归方程。然后根据回归方程得到在控制状态下（Ⅳ和Ⅴ类不超过 20%、劣Ⅴ类 没有）的废水率，然后与不加控制的预测值进行比较，得到需要处理的废水量（亿吨）。 2005 年 2006 年 2007 年 2008 年 2009 年 2010 年 2011 年 2012 年 2013 年 2014 年 79.4 97.5 118.9 141.9 166.1 192.6 231.3 261.3 294.8 330.1 对于问题五、根据我们所做模型的分析和预测结果，以及目前我国治理长江水质污 染的现状给出：合理制定工业发展规划、污水处理设施的利用、控制排污总量、加大执 法力度等四点建议和意见。 关键字：水质评价 模糊综合评价 综合指数评价 灰色系统 逐步回归
2005年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学) 蔡虎,卢亚军。许超指导教师汪晓银 正文 问题的提出 黄河脏了,长江可以救黄河,长江完了,拿什么救长江?”。我们的生命线一 长江正在倍受煎熬,保护长江、保护水资源就是保护我们自己。而如今长江污染日趋严 重。那么长江现在到底现况如何?污染源在哪里?如果再不采取有效措施,长江的命运 将会如何?最主要的污染物一一污水的排放对长江产生什么样的影响?怎么样规划、采 取怎么样的措施才能使长江在保持“生命力”的前提下达到环境与经济和谐发展?等等 这些都是我们亟待解决的问题 二、问题的分析 2.1对水质情况的定量综合评价 水质的评价是对水中各种物质的浓度进行测定,分析是否超过某个标准的过程。题 目中只给出了4种物质作为评价指标。通过查阅资料我们知道一般分析水质情况的方法 有指数评价法、模糊评价法、灰色评价法、物元分析法、人工神经网络评价法等诸多方 法叫。我们在这里先用模糊评价法作出评价,然后用比较简单常用的指数评价法进行验 22寻找主要污染源 污染源是指污染的源头,即污染物的排放点。但是题目只给出了各个测点的污染物 的浓度值,我们很容易理解一个测点的监测值高不一定说明该点就是主要污染源。因为 江水是流动的,一个测点的数据应该包括两个部分,即来自上游的污染物的影响和本地 的污染物的浓度。各测点本地污染物的浓度才是我们考察的目标。那么上游污染物对下 测点的影响度该怎么计算呢?这就要考虑两测点间距离、水流速度及污染物降解与上游 污染物影响度之间的关系 23水质污染趋势的预测 预测需要有历史数据支持,根据题目所给数据的特点:过去10年每一年中6类水 的比例,以及水流量总量及排污量等不确定因素。我们可以灰色系统方法进行预测。且 用一年中6类水的比例也可充分表明水质状况 24水质污染的控制 在2.3的预测中,我们知道6类水的比例反映水质污染的程度,且题目也是要求我 们通过对污水的处理达到对6类水的比例进行控制的目的,即6类水的比例是我们的控 制目标,污水量是我们的控制点。所以我们要先找到污水量与6类水比例之间的关系 我们采用多元回归对污水排放比例和6类水的比例进行逐步回归分析,这样我们可以得 到一个回归方程,它就是我们要找排污比例与6类水比例之间的关系式,通过它我们就 可以在给定控制标准下的污水排放上限,然后跟我们预测的不加控制的排污量比较就得 到了需要处理的污水量。 、模型的假设 (1)假设溶解氧(DO)浓度越髙水质越好,不考虑过含氧情况。 (2)假设各监测指标之间无相互作用。 (3)假设我们研究的长江是一条平直的河流。 (4)假设所给数据真实可靠。 (5)假设水质状况只与题目给我们的4个项目有关,不考虑其他项目 2
2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学) 蔡虎，卢亚军。许超 指导教师 汪晓银 2 正文： 一、问题的提出 “黄河脏了，长江可以救黄河，长江完了，拿什么救长江？”。我们的生命线—— 长江正在倍受煎熬，保护长江、保护水资源就是保护我们自己。而如今长江污染日趋严 重。那么长江现在到底现况如何？污染源在哪里？如果再不采取有效措施，长江的命运 将会如何？最主要的污染物——污水的排放对长江产生什么样的影响？怎么样规划、采 取怎么样的措施才能使长江在保持“生命力”的前提下达到环境与经济和谐发展？等等 这些都是我们亟待解决的问题。 二、问题的分析 2.1 对水质情况的定量综合评价 水质的评价是对水中各种物质的浓度进行测定，分析是否超过某个标准的过程。题 目中只给出了 4 种物质作为评价指标。通过查阅资料我们知道一般分析水质情况的方法 有指数评价法、模糊评价法、灰色评价法、物元分析法、人工神经网络评价法等诸多方 法[1]。我们在这里先用模糊评价法作出评价，然后用比较简单常用的指数评价法进行验 证。 2.2 寻找主要污染源 污染源是指污染的源头，即污染物的排放点。但是题目只给出了各个测点的污染物 的浓度值，我们很容易理解一个测点的监测值高不一定说明该点就是主要污染源。因为 江水是流动的，一个测点的数据应该包括两个部分，即来自上游的污染物的影响和本地 的污染物的浓度。各测点本地污染物的浓度才是我们考察的目标。那么上游污染物对下 测点的影响度该怎么计算呢？这就要考虑两测点间距离、水流速度及污染物降解与上游 污染物影响度之间的关系。 2.3 水质污染趋势的预测 预测需要有历史数据支持，根据题目所给数据的特点：过去 10 年每一年中 6 类水 的比例，以及水流量总量及排污量等不确定因素。我们可以灰色系统方法进行预测。且 用一年中 6 类水的比例也可充分表明水质状况。 2.4 水质污染的控制 在 2.3 的预测中，我们知道 6 类水的比例反映水质污染的程度，且题目也是要求我 们通过对污水的处理达到对 6 类水的比例进行控制的目的，即 6 类水的比例是我们的控 制目标，污水量是我们的控制点。所以我们要先找到污水量与 6 类水比例之间的关系， 我们采用多元回归对污水排放比例和 6 类水的比例进行逐步回归分析，这样我们可以得 到一个回归方程，它就是我们要找排污比例与 6 类水比例之间的关系式，通过它我们就 可以在给定控制标准下的污水排放上限，然后跟我们预测的不加控制的排污量比较就得 到了需要处理的污水量。 三、模型的假设 （1）假设溶解氧(DO)浓度越高水质越好，不考虑过含氧情况。 （2）假设各监测指标之间无相互作用。 （3）假设我们研究的长江是一条平直的河流。 （4）假设所给数据真实可靠。 （5）假设水质状况只与题目给我们的 4 个项目有关，不考虑其他项目
2005年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学) 蔡虎,卢亚军。许超指导教师汪晓银 四、号的定义与说明 41模糊评判模型符号定义 符号 符号说明 第i个参数在第j级别上的标准值 Y 第i个参数在某一级别上的监测值的第k个 i参数监测值介于Ay到A,之间的个数 L 第i个参数污染物监测值的个数 Z 参数污染物监测值介于AP到A.之间的1,个 监测值的平均值 z对第j级水质的隶属度 计算得到的i参数的权重 参数而言,介于A.到AJ之间的监测值发生在 J水质下的概率 归一化处理后i参数的权重 参数发在j水质下的模糊概率 水域水体出现j级水质的模糊综合概率 42综合指数评价模型符号定义 符号 符号说明 W 水质综合指标 K 分指数(监测数据与评价标准之比) 评价指标的权重 全流域17个测点COD的监测值 全流域17个测点NH3的监测值 P 全流域17个测点PH的监测值 D 全流域17个测点DO的监测值 支流10个测点COD的监测值 支流10个测点NH3的监测值 干流7个测点COD的监测值 N 干流7个测点NH3的监测值 五、模型的建立与求解 5.1长江水质的综合评价 5.1.1模糊综合评判模型 根据水域情况的质量标准我们把水污染监测浓度看成是一个离散的随机变量,用概 率统计方法进行统计可以得到水域属于某个标准的概率,因为可以拟定不同的水域标 准,评价参数集为U={u1,u2,u3},水质分级集为{v1,V2,v3,V45,v6},其中u1,u2,u3分别表 示为溶解氧,高锰酸盐指数,氨氮(NH3-N),因为PH值对水域影响不大,所以对其不
2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学) 蔡虎，卢亚军。许超 指导教师 汪晓银 3 四、 号的定义与说明 4.1 模糊评判模型符号定义 符号 符号说明 Aij 第 i 个参数在第 j 级别上的标准值 Xi,k 第 i 个参数在某一级别上的监测值的第 k 个 li,j i 参数监测值介于 Ai,j-1到 Ai,j之间的个数 Li 第 i 个参数污染物监测值的个数 Zi,j i 参数污染物监测值介于 Ai,j-1到 Ai,j之间的 li,j个 监测值的平均值 ri,j,j zi,j对第 j 级水质的隶属度 wi 计算得到的 i 参数的权重 pi,j i 参数而言，介于 Ai,j-1到 Ai,j之间的监测值发生在 j 水质下的概率 ai 归一化处理后 i 参数的权重 qi,j i 参数发在 j 水质下的模糊概率 qj 水域水体出现 j 级水质的模糊综合概率 4.2 综合指数评价模型符号定义 符号 符号说明 W 水质综合指标 K 分指数（监测数据与评价标准之比） η 评价指标的权重 Cj 全流域 17 个测点 COD 的监测值 Nj 全流域 17 个测点 NH3 的监测值 Pj 全流域 17 个测点 PH 的监测值 Dj 全流域 17 个测点 DO 的监测值 Czj 支流 10 个测点 COD 的监测值 Nzj 支流 10 个测点 NH3 的监测值 Cgj 干流 7 个测点 COD 的监测值 Ngj 干流 7 个测点 NH3 的监测值 五、模型的建立与求解 5.1 长江水质的综合评价 5.1.1 模糊综合评判模型[2] 根据水域情况的质量标准我们把水污染监测浓度看成是一个离散的随机变量，用概 率统计方法进行统计可以得到水域属于某个标准的概率，因为可以拟定不同的水域标 准，评价参数集为 U={u1,u2,u3},水质分级集为{v1,v2,v3,v4,v5,v6},其中 u1,u2,u3 分别表 示为溶解氧，高锰酸盐指数，氨氮（NH3-N），因为 PH 值对水域影响不大，所以对其不
2005年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学) 蔡虎,卢亚军。许超指导教师汪晓银 予考虑,-V6分别表示为I类到劣Ⅴ类,设i参数污染物监测值共有L个,其中介 于Ar到A之间的监测值有l个, 高锰酸盐指数,氨氮(NH3-N),的监测值为(i=2,3) 监测值Xk≤A1 l=〈lA1≤Hk≤A=2,3456) 溶解氧的监测值为(i=1) l,1监测值Xk≥A1 A≤Mk≤AA(=2……n) 那么,对于i参数而言,介于A1到A之间的监测值Xk发生在A1下或A下的 概率为p=l/L1 我们用隶属度来刻画水质分级界限: 已知水质等级标准为A1,i参数污染物监测值k介于A到A之间的l个监测值 的平均值为=(∑Xkl(k=1,2.bk为第i个参数污染物位于A到A之间的监 测值的第k个。则对;参数而言,对第j级水质和j-1级水质的隶属度m和r,可由 下列各式求得(其中产=2,3),即: r12:≤Acru=A}4j/A1-A Aij-l<<Aij(Ii,n+1,n=0 ≤A ri,1,1=0z,1=0 b≈00 而仁=1时对于溶解氧的隶属度的求法与上面方法相反 对于评价参数的权重的确定: 对于溶解氧权重按如下确定w=(×0-x1)x0-s)而高锰酸盐指数,氨氮的权重分别 为w=xs1,其中ⅹ-—-第ⅰ种污染物的实测浓度算术平均值,x0-—溶解氧在某条件下 的饱和浓度(标准浓度),s-—-第ⅰ种污染物各级标准的算术平均值。对其进行归 化处理得到a=w/∑;(r1,2,3) 3个参数构成权重矩阵即为A=(a1,a2,a3) 三个指标的权重(见表1):
2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学) 蔡虎，卢亚军。许超 指导教师 汪晓银 4 予考虑，v1-- v6 分别表示为Ⅰ类到劣Ⅴ类，设 i 参数污染物监测值共有 Li个，其中介 于 Ai,j-1到 Ai,j之间的监测值有 li,j 个， 高锰酸盐指数，氨氮（NH3-N），的监测值为(i=2,3) li,1 监测值 Xi,k ≤Ai,1 li,j= li,j Ai,j-1≤Xi,k≤Ai,,j(j=2,3,4,5,6) li,n+1 Xi,k ≥Ai,n(k=1,2,… … , Li) 溶解氧的监测值为(i=1)： li,1 监测值 Xi,k≥Ai,1 li,j= li,j-1 Ai,j-1≤Xi,k≤Ai,j(j=2,… … ,n) li,n Xi,k≥Ai,n 那么，对于 i 参数而言，介于 Ai,j-1 到 Ai,j 之间的监测值 Xi,k 发生在 Ai,j-1 下或 Ai,j 下 的 概率为 pi,j= li,j / Li 我们用隶属度来刻画水质分级界限： 已知水质等级标准为 Ai,j, i 参数污染物监测值 Xi,k 介于 Ai,j-1 到 Ai,j 之间的 li,j 个监测值 的平均值为 zi,j=(∑Xi,k)/ li,j(k=1 ,2,…li,j), Xi,k 为第 i个参数污染物位于 Ai,j-1 到 Ai,j 之间的监 测值的第 k 个。则对 i 参数而言，zi,j 对第 j 级水质和 j-1 级水质的隶属度 ri,j,j 和 ri,j,j-1,可由 下列各式求得(其中 i=2,3)，即： ri,1,1=1 zi,1≤Ai,j ri,j,j-1=Ai,j-zi,j／Ai,j- Ai,j-1 Ai,j-1≤zi,j≤Ai,j ri,n+1,n=0 ri,j,j=zi,j 一 Ai,j/／Ai,j 一 Ai,j-1 Ai,j-1≤zi,j≤Ai,j (zi,n+1=0) ri,j,j-1=0 zi,j=0 ri,n+1,n=1 ri,1,1=0 zi,1=0 ri,j,j=0 (zi,n+1> Ai,n 而 i=1 时对于溶解氧的隶属度的求法与上面方法相反 对于评价参数的权重的确定： 对于溶解氧权重按如下确定 w1=(x0-x1)/(x0-s1)， 而高锰酸盐指数，氨氮的权重分别 为 wi=xi/si ，其中 xi---第 i 种污染物的实测浓度算术平均值，x0---溶解氧在某条件下 的饱和浓度（标准浓度） ，si---第 i 种污染物各级标准的算术平均值。对其进行归一 化处理得到 ai= wi/∑wi（i=1,2,3） 3 个参数构成权重矩阵即为 A=（a1,a2,a3） 三个指标的权重（见表 1）：
2005年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学) 蔡虎,卢亚军。许超指导教师汪晓银 指标溶解氧 高锰酸盐指数氨氮 0.2523 0.6854 0.053 i参数发生各级水质的模糊概率为qi=F·pu+r1+1·P+(=1,23j=1,23456) 则模糊概率关系矩阵为q3×6于是每个单指标模糊概率评价矩阵为q×6=A*q36(见附表 1)其中q=Σa*q即为水域水体出现j级水质的模糊综合概率,并按最大隶属原则求 得q判断属于第几个级别。 在我们已知条件中我们对28个月中17个长江流域主要城市的水质监测项目求出平均值按以 上说明做出模糊综合评判,由 matlab程序求出17个城市28个月的平均水质情况,得到长江沿岸 17个观测点近两年多的水质平均情况(见表2): 表2 测点1234567891011121314151617 水质等级ⅡIIⅡⅡⅡⅡⅡIvⅢⅢIⅢⅡⅡ 劣 IlⅡ 5.1.2综合指数模型 为了验证我们评价的客观性,我们现在用综合指数评价模型对模糊综合评价模型进 行验证。 ①首先求分指数K 4个指标的K分别为:P/PH;D/DO;c/CODb;N/NH 其中选取第三类等级水的4项标准P4=7,DO=5,COD=6,NH=1为评价标准。由 于D0是逆指标,P是非单调指标(不是值越大表示水质越好或越坏)。所以综合指标计 算公式如下: P-/ P (不加权的) COD NH ②求权数n 权重的确定是参照附件3中的数据,通过统计28组数据中“主要污染指标中”各 指标出现次数,依据次数越多,权重越大的思想计算各权重(见表3) 表3 PH DO COD NHo 1.1% 25.9% 16.5% 56.5% W= n P n n (加权) Nh 因为有28个月,这样我们就得到了一个17×28的综合指标矩阵,数据过于庞大
2005 年全国大学生数学建模竞赛二等奖获奖论文(华中农业大学) 蔡虎，卢亚军。许超 指导教师 汪晓银 5 表 1 指标 溶解氧 高锰酸盐指数 氨氮 a 0.2523 0.6854 0.053 i 参数发生各级水质的模糊概率为 qi,j=ri,j,j * pi,j+ ri,j+1,j * pi,j+1(i=1,2,3,j=1,2,3,4,5,6), 则模糊概率关系矩阵为 q3×6,于是每个单指标模糊概率评价矩阵为 q1×6=A* q3×6（见附表 1）,其中 qj=∑ai *qi,j，即为水域水体出现 j 级水质的模糊综合概率，并按最大隶属原则求 得 qj 判断属于第几个级别。 在我们已知条件中我们对 28 个月中 17 个长江流域主要城市的水质监测项目求出平均值按以 上说明做出模糊综合评判，由 matlab 程序求出 17 个城市 28 个月的平均水质情况，得到长江沿岸 17 个观测点近两年多的水质平均情况（见表 2）： 表 2 测点 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 水质等级 Ⅱ Ⅱ Ⅱ Ⅱ Ⅱ Ⅱ Ⅱ IV Ⅱ Ⅱ Ⅰ Ⅲ Ⅱ Ⅱ 劣 Ⅴ Ⅱ Ⅱ 5.1.2 综合指数模型[3] 为了验证我们评价的客观性，我们现在用综合指数评价模型对模糊综合评价模型进 行验证。 ①首先求分指数 K 4 个指标的 K 分别为：Pj / PH0 ； Dj / DO0 ； Cj / COD0 ；Nj / NH0 其中选取第三类等级水的 4 项标准 PH0=7，DO0=5，COD0=6，NH0=1 为评价标准。由 于 DO 是逆指标，PH 是非单调指标（不是值越大表示水质越好或越坏）。所以综合指标计 算公式如下： W= 1 0 - PH Pj + D j DO0 + COD0 C j + NH0 N j （不加权的） ②求权数η 权重的确定是参照附件 3 中的数据，通过统计 28 组数据中“主要污染指标中”各 指标出现次数，依据次数越多，权重越大的思想计算各权重（见表 3） 表 3 PH DO0 COD0 NH0 η 1.1% 25.9% 16.5% 56.5% W=η1 1 0 - PH Pj +η2 D j DO0 +η3 COD0 C j +η4 NH0 N j （加权） 因为有 28 个月，这样我们就得到了一个 17×28 的综合指标矩阵，数据过于庞大