关于长江水质的综合评价与预测问题.doc

2014高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了全国大学生数学建模竞赛章程和全国大学生数学建模竞赛参赛规则（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的参赛报名号为（如果赛区设置报名号的话）： 所属学校（请填写完整的全名）： 参赛队员 (打印并签名) ：1. 郝琪琪 2. 高盈超 3. 指导教师或指导教师组负责人 (打印并签名)： 日期： 年 月 日赛区评阅编号（由赛区组委会评阅前进行编号）：2013高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：- 2 -关于长江水质的综合评价与预测问题摘要作为世界第三大河流的长江，面临着前所未有的水资源污染问题，由于污染严重，长江岸边形成了许多污染带，在干流21个城市中，重庆、岳阳、武汉、南京、镇江、上海六大城市累计污染带长度占长江干流污染带总长的73%。本文针对长江水质的水污染问题，进行长江水质的综合评价与其污染程度进行预测。针对问题一，对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况。通过对附表: 地表水环境质量标准（GB38382002）中4个主要项目标准限值的分析，利用层次分析法，确定溶解氧、高锰酸盐指数、氨氮、PH值四个因素对长江水质的影响程度，也即是各个所对应得权重。同时，再结合模糊综合评价模型，利用最大隶属度的原则，对长江水质的四个影响指标合理准确的评估，从而对长江近两年的水质情况作出定量的综合评价。并利用模糊模型的结果对各地区水质污染状况进行分析。针对问题三，假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析。利用附件3中的各种数据，借助MATLAB软件建立灰色预测GM(1,1)模型，对其残差检验、关联度检验、后验差检验进行验算，较为科学的预测出长江水未来十年的水质状况。利用灰色模型处理数据不仅对数据没有很强的限制，而且精度高，计算简便，能科学的预测出未来十年水质污染的发展趋势。针对问题四，通过建立非线性规划模型，借助附表四的数据和问题三已经解决的未来十年的预测问题，可以明确长江处理污水排放量的函数关系。针对问题五，提出治理水质的切实可行的办法这一要求，我们从分析的数据与问题出发，分别就问题的根源、影响因素与较高权重等方面，较为科学、全面、系统的提出了问题的解决方案。综上所述：近两年在长江区域主要污染物起主导作用的是溶解氧，次之是高锰酸盐指数，高锰酸盐浓度主要对四川乐门沱江二桥影响最为显著，尤其是在2013年12月到2014年六月影响最为显著。对四川泸州沱江二桥影响次之，在这两个地区出现较高的高锰酸盐浓度，对其水质影响严重。氨氮浓度主要对江西南昌滁槎影响最为显著，在个别月里严重超标，影响尤为恶劣。关键词：层次分析法 模糊综合评价 GM(1,1)模型 非线性规划模型一、问题重述作为世界第三大河流的长江，面临着前所未有的水资源污染问题，由于污染严重，长江岸边形成了许多污染带，在干流21个城市中，重庆、岳阳、武汉、南京、镇江、上海六大城市累计污染带长度占长江干流污染带总长的73%。本文针对长江水质的水污染问题，进行长江水质的综合评价与其污染程度进行预测。要解决的问题如下：（1）对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况。建立基于层次分析法的迷糊综合评价模型，来解决实际问题。（2）研究、分析长江干流近一年多主要污染物高锰酸盐指数和氨氮的污染源主在哪要些地区。建立统计模型，来分析长江干流的主要污染物的污染状况。（3）假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析，建立GM(1，1)灰色预测模型，科学的对长江污染发展趋势作出预测。（4）根据预测分析，如果未来10年内每年都要求长江干流的类和类水的比例控制在20%以内，且没有劣类水(约束条件),那么每年需要处理的污水。建立非线性回归方程，借助软件行合理规划。二、问题的分析针对问题一，对长江近两年多的水质情况做出定量的综合评价，并分析各地区水质的污染状况。通过对附表: 地表水环境质量标准（GB38382002）中4个主要项目标准限值的分析，利用层次分析法，确定溶解氧、高锰酸盐指数、氨氮、PH值四个因素对长江水质的影响程度，也即是各个所对应得权重。同时，再结合模糊综合评价模型，利用最大隶属度的原则，对长江水质的四个影响指标合理准确的评估，从而对长江近两年的水质情况作出定量的综合评价。并利用模糊模型的结果对各地区水质污染状况进行分析。针对问题三，假如不采取更有效的治理措施，依照过去10年的主要统计数据，对长江未来水质污染的发展趋势做出预测分析。利用附件3中的各种数据，借助MATLAB软件建立灰色预测GM(1,1)模型，对其残差检验、关联度检验、后验差检验进行验算，较为科学的预测出长江水未来十年的水质状况。利用灰色模型处理数据不仅对数据没有很强的限制，而且精度高，计算简便，能科学的预测出未来十年水质污染的发展趋势。针对问题四，通过建立非线性规划模型，借助附表四的数据和问题三已经解决的未来十年的预测问题，可以明确长江处理污水排放量的函数关系。针对问题五，提出治理水质的切实可行的办法这一要求，我们从分析的数据与问题出发，分别就问题的根源、影响因素与较高权重等方面，较为科学、全面、系统的提出了问题的解决方案。三、对问题的合理假设（1）假设各个测量数据不存在人为的影响因素。（2）一个观测站的水质污染主要来自本地区的排污和上游的污水，河水在通过观测点后立即被混入下段区域的污水。（3）不考虑洪水干旱等特殊气候对水质的影响。 （4）两观测点之间江水的流速是所测的两点流速的平均值。四、符号说明 符号表示含义Q长江水质量影响指标W判断矩阵权重归一化向量t判断矩阵最大特征根表示溶解氧,高锰酸盐指数,氨氮,PH值v1,v2,v6各个水质类别R=(rig)hxm评价矩阵R1，R2本月，上月评价S模型评价率五、问题一的解决方案（1）通过附表提供的数据利用层次分析法把问题层次结构分两层：第一层为目标层(O);第二层为准则层（C），相关条件，共有四个因素，依次是A溶解氧(DO)、B高锰酸盐指数(CODMn)、 C氨氮（NH3-N）、D PH值（无量纲）。构造的判断矩阵如表一所示。Q长江水质影响影响A溶解氧B高锰酸盐指数（CODMn）D PH值C氨氮（NH3-N）表一 比较判断矩阵WABCDA1234B1/2123C1/31/212D1/41/31/21借助MATLAB软件计算出该判断矩阵的权重归一化向量w和最大特征根t, w= 0.4673 0.2772 0.1601 0.0954 t= 4.0310 若想使结果具有科学性需对归一化向量进行一致性检测，如果满足CR0.1，则结果具有说服力，否则需从新构造判断矩阵，再进行一致性检验，知道满足条件为止。对w向量进行一致性检验，一致性指标RI的取值如表二所示123456789101100058090112124132141145149151也即RI取值为0.90，则一致性指标CI= ,CI=0.01033，一致性比率指标CR=,CR=0.011480.6的。第七步 后验差检验（1）计算： （2） 计算 序列的均方差 （3） 计算残差的均值： （4） 计算残差的均方差: （5） 计算C: （6） 计算小残差概率： 所有的都小于，故小残差概率Pp i=i+1; x(:,i)=A*y(:,i-1); m(i)=max(x(:,i); y(:,i)=x(:,i)/m(i); k=abs(m(i)-m(i-1);enda=sum(y(:,i);w=y(:,i)/a;t=m(i);disp(w);disp(t);请输入判断矩阵A(n阶)A=1 2 3 4;1/2 1 2 3;1/3 1/2 1 2;1/4 1/3 1/2 1 w= 0.4673 0.2772 0.1601 0.0954 t= 4.0310(b)模糊因子算法：function ab=synt2(a,b);m=size(a,1);n=size(b,2);for i=1:m for j=1:n ab(i,j)=max(a(i,:)b(:,j) endend（c）数据拟合程序：x=1995:1:2004; z=174 179 183 189 207 234 220.5 256 270 285; a=polyfit(x,z,2);Warning: Polynomial is badly conditioned. Remove repeated data points or try centering and scaling as described in HELP POLYFIT. In polyfit at 81 plot(x,z,+,x,s,r) ylabel(废水排放总量/亿吨); xlabel(年份); title(1995-2004长江流域废水排放总量) grid on（d）GM(1,1)灰色模型源程序：x0=174 179 183 189 207 234 220.5 256 270 285; n=length(x0); lamda=x0(1:n-1)./x0(2:n); range=minmax(lamda); %计算及比的范围 x1=cumsum(x0); %累加运算 B=-0.5*(x1(1:n-1)+x1(2:n),ones(n-1,1); Y=x0(2:n); u=BY; x=