动态加权综合评价、灰色关联度分析、BP神经网络模型.doc

大气污染预报问题摘要本文针对大气污染问题，采用动态加权综合评价方法建立了合理的空气质量评价模型，同时，采用灰色关联度分析方法和BP神经网络模型较好地研究了空气质量和气象参数之间的关系。问题一中，考虑到污染物浓度这一评价指标的“质的差异”和“量的差异”，采用动态加权综合评价方法建立评价模型。首先对评价指标数据进行归一化处理，然后选取偏大型正态分布函数作为动态加权函数建立评价模型，从而对评价指标每天的观测值进行排序，最后用决策分析中的Borda数方法对四个城市的空气质量综合排序。得到的最终排序结果为：空气质量最差的是城市，其次是城市，排在第三位是城市，而城市的空气质量最好。问题二中，对于空气质量与气象参数关系的问题，采用灰色关联度分析方法和BP神经网络模型进行探讨。首先，通过灰色关联度分析确定了大气污染物浓度与气象参数强弱主次关系，然后针对其复杂非线性关系建立BP神经网络预测模型，预测2009年7月26日至30日的污染物浓度。最后用实际值对预测值进行了误差分析，结果表明预测值与实际值的误差较小，即BP神经网络模型的预测值具有较高的精度。 本文最大的特色是采用了动态加权综合评价方法建立评价模型，增大了评价结果的客观性，比定常加权模型更科学合理。其次，鉴于空气质量与气象参数复杂的非线性关系，建立了BP神经网络模型，较好地讨论了大气污染物浓度与气象参数的关系，经过检验分析知此模型是解决非线性问题的有力工具。关键词：动态加权综合评价、灰色关联度分析、BP神经网络模型、一、问题提出大气是指包围在地球外围的空气层，是地球自然环境的重要组成部分之一。人类生活在大气里，洁净大气是人类赖于生存的必要条件。一个人在五个星期内不吃饭或5天内不喝水，尚能维持生命，但超过5分钟不呼吸空气，便会死亡。随着地球上人口的急剧增加，人类经济增长的急速增大，地球上的大气污染日趋严重，其影响也日趋深刻，如由于一些有害气体的大量排放，不仅造成局部地区大气的污染，而且影响到全球性的气候变化。因此，加强大气质量的监测和预报是非常必要。目前对大气质量的监测主要是监测大气中、悬浮颗粒物（主要为）等的浓度。附件一给出了城市A、B、C、D从2009年6月1日至2009年7月25日测量的污染物含量及城市A的气象参数的数据；附件二给出了城市A从2009年7月26日至2009年7月30日测量的污染物含量及气象参数的数据。请解决下面两个问题：（1）建立由污染物浓度评价空气质量的数学模型，然后利用附件一中的数据对四个城市的空气质量进行排序。（2）分析城市A的空气质量(指、PM10的浓度)与气象参数之间的关系，并利用附件二中的数据进行检验。二、基本假设1、假设评价空气质量的各指标间相互作用关系忽略不计；2、假设空气质量仅与附件中的四个气象参数有关；3、假设题中数据为每天的统计平均值，能客观反映当天空气污染物浓度的实际情况；4、假设在预测模型中，在未来一年没有发生重大自然突变；三、符号说明符号意义备注表示第个评价指标的第个观测值的原始值表示第个评价指标的第个观测值归一后的值表示动态加权函数表示第个评价对象的第个综合评价指标值表示第个排序方案中排在第个被评价对象后的个数表示第个气象参数对评价指标的第个点的关联系数表示第个气象参数与评价指标的关联度四、问题分析大气污染问题愈加严重，加强大气质量的监测和预报十分必要。问题中要求解决两个问题：一是建立由污染物浓度评价空气质量的数学模型，据此对四个城市的空气质量进行排序；二是需要分析空气质量与气象参数之间关系，并利用数据检验分析。针对问题一，查阅国标（GB 3095-1996）规定，环境空气质量标准分为三级。每一等级对每一项指标都有相应的标准值，且同一等级的空气在污染物的含量上也有差别。这种既有“质的差异”又有“量的差异”的问题可采用动态加权评价方法建立评价模型，并利用决策分析中的Borda函数方法确定最终的排序方案。针对问题二，通过观察附件一、二中的数据可知，需要利用附件一中数据建立预测模型，并用问题二中数据检验分析。故可先利用灰色关联度分析方法讨论四种气象参数分别对三种污染物浓度的影响程度，为了进一步讨论空气质量与气象参数的关系，建立BP神经网络预测模型，并利用附件二中数据检验分析。 具体各问问题分析见本文（五）。五、模型的建立与求解5.1 问题一模型建立与求解5.1.1 问题一的分析问题一要求建立由污染物浓度评价空气质量的数学模型，并对四个城市空气质量进行排序。查阅国标（GB 3095-1996）规定知，环境空气质量标准分为三级，每一个级别对每一项指标都有相应的标准值（相关数据见表1）。也就是说对于每一个评价指标（即、的浓度），既有同级别的差异，同级别又有不同量值的差异。对于这种既有“质差”，又有“量差”的问题，采用定常加权法显然是不合理的，故合理有效地方法是动态加权综合评价方法。表1：国标GB 3095-1996中3个主要项目标准限值 （单位：）污染物名称日平均浓度限值范围一级标准二级标准三级标准(0,0.05(0.05,0.15(0.15,0.25(0,0.08(0,0.08(0.08,0.12(0,0.05(0.05,0.15(0.15,0.25 5.1.2 问题一模型的建立评价指标的规范化处理因为评价指标可能有极大型的、极小型的、中间型、或区间型四种情况，且可能各有不同的量纲，故需要对不同类型的指标变换成统一的、无量纲的标准化指标。1评价指标类型的一致化处理通过判断可知、的浓度这三类指标均为极小型指标，即总是期望指标的取值越小越好。故在此不需要将进行指标类型一致化处理。2评价指标无量纲化处理一般来说，数据的无量纲化处理有标准差方法、极值差方法和功效系数方法等。在此，选取极值差方法对三类指标进行标准化处理。处理方法为： （1.1）其中，则是无量纲的指标观测值。经过此处理我们可以得到评价指标经标准化处理后的三级区间。动态加权函数的确定根据空气质量问题的实际问题，各项指标对综合评价的影响比较符合随着类别的增加呈现先缓慢增加，中间快速增长，最后平缓增加趋于最大。于是不妨选取呈正态分布曲线的偏大型正态分布函数作为动态加权函数，即： （1.2）其中在这里取指标的第一级浓度标准区间的中值，即，由确定。由标准化处理后的实际数据经计算可得， ，则代入上式可以得到、三项指标的动态加权函数。图1：加权函数的图像建立综合评价模型为了给每次的检测值进行排序，在基于上述模型的同时，取综合评价模型为各评价指标的动态加权和，即： （1.3）由此综合评价指标函数可以求出每个评价对象的个综合评价指标值，且据此大小排序，可得到个评价对象的个排序方案。利用决策分析中的Borda函数方法来确定综合排序方案。若在第个排序方案中排在第个被评价对象后的个数为，令，则被评价对象的Borda数为 （1.4）根据此式的计算结果大小进行排序，便可得到个被评价对象的总排序结果。综合以上的数据处理、权的求法及排序函数的数学模型，建立以下由污染物浓度评价空气质量的数学模型：5.1.3 问题一模型的求解与分析1.算法（1）运用极值差法，将数据先归一化化成可比较的区间上的数值（2）根据偏大型正态分布函数，确定三类指标的动态加权函数（3）从1开始到3，从1到4，从1到55，对（2）得到的新数值矩阵进行加权求和，得到不同的的分数（4）对分数进行从大到小的排序，得到4个评价对象的55个排序方案。（5）利用Borda函数计算4个被评价对象的Borda数，并根据Borda数从大到小进行排序，得到最终排序结果2.求解及分析（1）求解运用软件编程（见附录1）对各次检测值进行加权求和，得到各综合评价指标值：部分结果见下表;表2 :四个城市的综合评价指标值城市得分日期2009/6/10.034030.000640.098080.348542009/6/20.002760.103050.237950.239212009/6/30.000390.000200.059980.061322009/6/40.052170.123630.240070.666142009/7/240.005030.014140.006640.011622009/7/250.000390.017340.006520.00379根据表2的结果进行排序得到55个排序方案后，利用Borda函数编程计算可得的Borda数及总排序结果如表3 所示表 3 ：按各城市空气污染物浓度总排序结果城市排序Borda数68928981总排序4123由上表知四个城市的空气污染物浓度排序为：，可见空气污染最严重的是城市，其次是城市，排在第三位的是城市，而空气质量最好的是城市。5.2 问题二模型建立与求解5.2.1 问题二的分析问题二要求解决以下两个问题：一是分析城市A的空气质量与气象参数之间的关系，二是利用附件二中的数据进行检验分析。根据常识，大气污染物与气象条件有着紧密联系，特别是风速对其影响非常明显，而大气压、温度、湿度对大气污染物浓度的影响并不明确。故可以通过灰色关联分析方法分析大气污染物、的浓度与气象参数之间的强弱主次关系。通过观察附件一和附件二中数据的关系，可知需要建立空气质量与气象参数相关的预测模型。而空气质量与气象参数之间存在着复杂的非线性关系，对于非线性关系研究，可以通过神经网络模型，二次多项式逼近等方法解决。故在此刻建立BP神经网络预测模型，并利用附件二中数据对预测值进行检验分析。5.2.2 问题二模型的建立灰色关联分析方法模型的建立灰色关联度分析法是灰色系统理论中一种定量描述因素间发展势态的相似或相异程度的量化比较方法。它的基本思想是根据序列曲线几何形状的相似程度来判断其联系是否紧密。一般地，曲线越接近，相应序列之间的关联度就越大，反之就越小。灰色关联度分析法的步骤与模型如下：1分别选取大气污染物、PM10的日平均浓度数列为参考数列，气象参数数列为比较数列，其中参考数列记为，比较数列记为。2由于气象参数数列中变量的量纲不同，为消除量纲对分析结果的影响，需要进行无量纲化处理。常用的有标准化、初始化、极差法、最大值化等，在此也选用极差法进行处理，如下： 令 （2.1）其中，则是无量纲的指标观测值。3计算关联度系数，公式为 （2.2）其中：为和第个点的绝对误差；为两级最小差；为分辨率，一般取；越大，分辨率越小，越小，分辨率越大。4计算关联度，其公式为 （2.3）其中，即为对的关联度BP神经网络进行预测模型的建立通过上述灰色关联度分析，明确了大气污染物、的浓度与气象参数之间的强弱主次关系。但是大气污染物浓度与气象参数存在复杂的非线性关系，故在此建立BP神经网络模型进一步研究。本文采用的是多层感知器的误差反向传播算法(Error Back Proragation),即BP算法。其基本思想是，学习过程由信号的正向传播与误差的反向传播两个过程组成。正向传播时，输入样本从输入层传入，经各隐层逐层处理后，传向输出层。若输出层的实际输出与期望的输出不符，则转入误差的反向传输阶段。误差反传是将输出误差以某种形式通过隐层向输入曾逐层反传，并将误差分摊给各层的所有单元，从而获得各层单元的误差信号，该误差信号即作为修正各单元权值的依据。信号正向传播与误差反向传播的各层权值的调整是反复进行的，直至网络输出的误差减少到可接受的程度，或是进行到预先设定的学习次数，结构图如下。 符号说明： 输入变量向量组（四个变量） 隐层输出变量（为隐变量个数）向量组 输出层输出向量组（大气污染物浓度指标） 期望输出向量组（为常向量组） 输入层到隐层之间的权值矩阵 隐层到输出层之间的权值矩阵 比例常数，表示学习速率用表示四个气象参数变量（大气压、温度、湿度、风速），隐层输出变量设为7个，用表示，输出层输出向量（大气污染物浓度）用表示，而期望输出为常向量根据BP神经网络原理，对于输出层，有 （2.4） 对于隐层，有 （2.5） 其中 （2.6）5.2.3 问题二模型的求解5.2.3.1灰色关联度分析模型求解首先，将代入公式，用编程（见附录3 ）可计算得到大气污染物、浓度与各气象要素的关联系数，部分结果见表4 表 4 ：大气污染物与各个气象参素的关联系数0.753930.774740.648150.583680.924970.960510.528390.659860.768560.746990.449730.909360.928320.695180.730090.754710.949250.744310.640380.822280.899680.79220.629470.875150.453590.641190.877920.49840.491350.76280.679730.551680.433810.615360.900560.477060.712480.955440.62740.681240.969750.679350.479290.514150.681160.994590.631950.689850.532030.454120.57640.530780.745290.586930.847960.742550.554280.46640.605480.552860.382710.612890.464980.420540.438930.849960.568430.496430.389040.651110.480070.43055然后将表4 中的关联系数代入（2.3）式计算得大气污染物、浓度与各气象要素的关联度，如下表表 5：大气污染物浓度与各气象参数的关联度0.621830.705600.639750.665090.649060.701040.662890.631070.626540.698980.670280.69972通过对表5中的数据分析可知：影响浓度的气象参数排序为：温度风速湿度大气压影响浓度的气象参数排序为：温度湿度大气压风速影响浓度的气象参数排序为：风速温度湿度大气压综上，可知四种气象参数对大气污染物、浓度的影响程度不同，其中对和影响程度最大的是温度，而风速则是影响浓度的主要因子。5.2.3.2 BP神经网络模型求解算法Step1：初始化。对权值矩阵、赋随机数，将样本模式计数器和训练次数计数置为1，误差置0，学习率设为内的小数，网络训练后达到的精度设为一个正的小数。Step2：输入训练样本对，计算各层输出。用当前样本对向量数组赋值，并计算各分量。Step3：计算网络输出误差。共有6对训练样本，网络对于第个样本具有误差，总输出误差（均方根误差）为 Step4：计算各层误差信号Step5：调整各层权值。权值调整量为Step:6：检查是否对所有样本完成一次轮训。若，计数器加1，返回Step2，否则进行下一步。Step7：检查网络总误差是否达到精度要求。若，训练结束，否则置0，置1，返回Step2结果使用Matlab编程计算（程序见附录4），训练后的神经网络模拟的指标与实际数据几乎完全吻合，最终得到的、浓度预测值如下：表6：大气污染物浓度预测值2009/7/260.03460.0279220.0656572009/7/270.0323850.0288860.0639062009/7/280.0276930.0257160.0573132009/7/290.0304560.0285720.0617492009/7/300.0325950.0298260.0646695.2.4 问题二结果的分析及验证由于通过神经网络模型预测出了2009年7月26至2009年7月30日的大气污染物浓度，故可通过附录二中这五天的实际值进行检验分析。实际值与预测值对比如下：表7：浓度实际值与预测值的统计比较结果实际值预测值绝对误差2009/7/260.0310.03460.00359992009/7/270.0210.0323850.0113852009/7/280.0250.0276930.00269292009/7/290.0240.0304560.00645572009/7/300.0260.0325950.0065946表8：浓度实际值与预测值的统计比较结果实际值预测值绝对误差2009/7/260.0370.027922-0.0090782009/7/270.0220.0288860.00688612009/7/280.0340.025716-0.00828422009/7/290.0350.028572-0.0064282009/7/300.0330.029826-0.0031743表9：浓度实际值与预测值的统计比较结果实际值预测值绝对误差2009/7/260.0470.0656570.0186572009/7/270.030.0639060.0339062009/7/280.0340.0573130.0233132009/7/290.0350.0617490.0267492009/7/300.0810.064669-0.016331由上述所得的绝对误差，用编程绘制出误差曲线如下图：图2：误差曲线图从图2中可见5天内、浓度的预测值与实际值的绝对误差波动不大，通过计算得浓度的平均绝对误差为0.0061，浓度的平均绝对误差为-0.004，浓度的平均绝对误差为0.017。故BP神经网络模型预测值的误差较小，且波动相对稳定，预测效果理想。而通过图中残差平方和、检验值以及组内离均差平方和判断可知此模型的拟合优度较高，预测精度较高。六、模型的评价与推广7.1 模型的评价问题一中，考虑到大气污染物浓度这一评价指标的“质差”和“量差”的关系，采用了动态加权综合评价的方法建立模型，最后利用Borda函数得到了四个城市空气质量的总排序。这一模型不仅充分考虑到了每一因素每一属性的所有“差异”的影响和作用，使得评价结果科学合理；而且增加了综合评价的客观性，与定常加权法相比大大淡化了评价人主观因素的影响。问题二中，灰色关联度分析模型清晰明了地反映了大气污染物、的浓度与气象参数之间的强弱主次关系。而BP神经网络预测模型则更进一步对污染物浓度与气象参数间复杂的非线性关系进行了研究，用实际值对预测值检验分析后知该模型精度较高，预测结果理想。但是，BP神经网络模型也存在一定的缺点：隐层节点数需要经过多次尝试；学习速度慢；易于陷入局部极小点。7.2 模型的推广问题一中的动态加权综合评价模型不仅可以科学地评价空气质量问题，还可以推广到水质评价问题以及经济、军事等领域的许多综合评价问题，具有广泛的应用价值。问题二中的BP神经网络模型是目前最为广泛应用的神经网络模型之一，在气象预报中得到了广泛应用。该模型是一种解决非线性问题的有力工具，可以推广的其他领域应用。但由于BP神经网络模型自身的局限性，在实际应用中需要进行一定的改进和修正，本文中的BP神经网络模型若进行改进效果会更好。七、参考文献1 中华人民共和国国家标准/环境空气质量（GB3095-1996），http:/www.stats.gov.cn/tjsj/qtsj/hjtjzl/hjtjsj1999.htm，1996年2 姜启源等, 数学模型（第三版），高等教育出版社，2003年8月3 王庚、王敏生，现代数学建模方法 ，科学出版社，2008年2月4 苏静芝等，神经网络在空气污染预报中的应用研究，四川环境，2008年4月5 张德丰，MATLAB神经网络应用设计，机械工业出版社，2009年1月八、附录8.1 附录清单附录1：求解问题一的程序附录2：问题一的完整数据结果附录3：求解问题二中灰色关联度分析的程序附录4：求解问题二中BP神经网络模型的程序8.2 附录正文附录1：求解问题一的程序首先是求解动态加权函数中的值x=0.6,0.67,0.6;ca=0.1,0.335,0.1;ca1=x-ca./sqrtm(-log(0.1);从而求出结果ca1 = 0.53410000000000 0.44920000000000 0.53410000000000然后编程求解评价指标值和最后总得分A=0.037 0.031 0.098 0.030 0.032 0.058 0.029 0.028 0.028 0.025 0.034 0.108 0.034 0.027 0.067 0.034 0.032 0.060 0.029 0.024 0.059 0.028 0.028 0.094 0.027 0.025 0.039 0.018 0.028 0.052 0.014 0.020 0.031 0.011 0.023 0.061 0.040 0.035 0.067 0.039 0.042 0.065 0.047 0.033 0.033 0.036 0.045 0.090 0.035 0.030 0.076 0.036 0.041 0.068 0.011 0.036 0.040 0.025 0.039 0.106 0.027 0.021 0.056 0.025 0.030 0.041 0.000 0.000 0.000 0.010 0.025 0.080 0.026 0.023 0.047 0.023 0.029 0.048 0.011 0.019 0.024 0.013 0.027 0.081 0.026 0.023 0.061 0.017 0.028 0.050 0.007 0.021 0.037 0.014 0.027 0.070 0.027 0.019 0.051 0.010 0.030 0.034 0.005 0.008 0.015 0.010 0.020 0.050 0.025 0.018 0.038 0.009 0.023 0.024 0.005 0.014 0.011 0.009 0.015 0.037 0.036 0.024 0.061 0.013 0.028 0.031 0.011 0.025 0.019 0.015 0.026 0.061 0.028 0.017 0.055 0.016 0.027 0.034 0.006 0.013 0.017 0.011 0.021 0.062 0.032 0.025 0.050 0.023 0.033 0.042 0.005 0.023 0.022 0.014 0.025 0.060 0.035 0.019 0.049 0.051 0.040 0.046 0.014 0.028 0.029 0.016 0.025 0.061 0.037 0.027 0.074 0.081 0.054 0.052 0.017 0.037 0.034 0.019 0.041 0.082 0.044 0.031 0.128 0.100 0.062 0.094 0.013 0.036 0.077 0.019 0.047 0.136 0.029 0.029 0.102 0.019 0.041 0.071 0.003 0.048 0.067 0.005 0.040 0.111 0.033 0.026 0.085 0.019 0.031 0.059 0.006 0.031 0.060 0.005 0.030 0.094 0.031 0.023 0.095 0.009 0.023 0.054 0.004 0.032 0.053 0.003 0.022 0.078 0.045 0.023 0.129 0.017 0.029 0.084 0.017 0.023 0.057 0.022 0.034 0.158 0.032 0.022 0.105 0.011 0.026 0.066 0.009 0.015 0.050 0.014 0.026 0.123 0.033 0.021 0.068 0.020 0.022 0.047 0.021 0.025 0.029 0.018 0.022 0.085 0.033 0.024 0.075 0.018 0.029 0.043 0.016 0.020 0.028 0.016 0.027 0.071 0.056 0.035 0.100 0.036 0.038 0.048 0.042 0.033 0.033 0.033 0.040 0.125 0.059 0.041 0.101 0.045 0.051 0.089 0.077 0.046 0.054 0.032 0.047 0.118 0.065 0.049 0.148 0.081 0.052 0.127 0.054 0.055 0.084 0.037 0.000 0.157 0.051 0.030 0.068 0.029 0.039 0.052 0.013 0.032 0.054 0.010 0.032 0.079 0.044 0.021 0.101 0.010 0.026 0.081 0.030 0.018 0.062 0.014 0.014 0.128 0.066 0.024 0.062 0.011 0.032 0.049 0.017 0.018 0.022 0.017 0.027 0.079 0.035 0.018 0.053 0.008 0.023 0.049 0.008 0.009 0.021 0.009 0.004 0.045 0.052 0.023 0.038 0.046 0.025 0.042 0.019 0.013 0.020 0.024 0.014 0.060 0.045 0.025 0.030 0.056 0.028 0.038 0.025 0.017 0.020 0.023 0.016 0.050 0.057 0.034 0.063 0.061 0.041 0.050 0.063 0.031 0.030 0.019 0.027 0.082 0.076 0.043 0.093 0.094 0.055 0.081 0.030 0.047 0.049 0.027 0.058 0.035 0.073 0.036 0.089 0.000 0.067 0.085 0.020 0.021 0.058 0.020 0.025 0.000 0.056 0.028 0.077 0.000 0.057 0.073 0.009 0.020 0.040 0.013 0.019 0.000 0.070 0.031 0.051 0.035 0.037 0.045 0.031 0.029 0.025 0.024 0.031 0.092 0.049 0.029 0.061 0.019 0.028 0.044 0.013 0.024 0.033 0.008 0.028 0.056 0.034 0.026 0.041 0.019 0.031 0.041 0.004 0.032 0.018 0.007 0.031 0.051 0.037 0.040 0.096 0.038 0.050 0.073 0.013 0.052 0.032 0.020 0.057 0.130 0.026 0.027 0.058 0.028 0.030 0.044 0.000 0.000 0.000 0.012 0.035 0.047 0.023 0.029 0.065 0.066 0.027 0.050 0.011 0.042 0.035 0.009 0.034 0.045 0.034 0.031 0.079 0.091 0.038 0.066 0.023 0.035 0.031 0.028 0.044 0.100 0.022 0.028 0.050 0.071 0.033 0.040 0.003 0.046 0.029 0.004 0.033 0.038 0.032 0.043 0.063 0.049 0.046 0.056 0.011 0.032 0.015 0.017 0.048 0.070 0.033 0.035 0.078 0.011 0.031 0.057 0.006 0.034 0.034 0.005 0.041 0.064 0.024 0.031 0.057 0.004 0.032 0.055 0.004 0.019 0.037 0.003 0.036 0.032 0.029 0.030 0.040 0.006 0.032 0.040 0.004 0.015 0.024 0.004 0.029 0.021 0.024 0.052 0.109 0.010 0.050 0.078 0.004 0.031 0.033 0.006 0.060 0.099 0.022 0.026 0.052 0.027 0.028 0.046 0.003 0.012 0.031 0.013 0.027 0.035 0.030 0.043 0.068 0.063 0.045 0.069 0.013 0.049 0.031 0.023 0.056 0.000 0.041 0.037 0.075 0.039 0.036 0.044 0.018 0.037 0.020 0.021 0.047 0.064 0.068 0.036 0.060 0.061 0.031 0.047 0.013 0.038 0.032 0.014 0.038 0.054 0.024 0.028 0.041 0.059 0.022 0.033 0.011 0.018 0.022 0.013 0.022 0.042 0.022 0.025 0.037 0.057 0.021 0.027 0.004 0.018 0.013 0.012 0.018 0.029 ;ca=0.1,0.335,0.1;%中值ca1=0.5341 0.4492 0.5341;%的值B=0.2 0.6 1;0.67 0.67 1;0.2 0.6 1;%标准分级a=size(A,1);%A的行数b=size(A,2);%A的列数MX=max(A);%A的每列的最大值MN=min(A);%A的每列的最小值f1=1;%B的行变量f2=1;%B的列变量for j=1:b%将A里面的值标准化 for i=1:a A(i,j)=(A(i,j)-MN(j)/(MX(j)-MN(j); endendX=zeros(a,4);%建立放置每个城市每天空气质量的值矩阵for i=1:a%计算每个城市每天的空气质量 h=1; k=1; flag0=1;%算完每个城市每天的空气质量后，跳到下一个城市（h=h+1） for j=1:b for flag=1:3 if A(i,j)ca(k)&A(i,j)A(i,j) C(k)= C(k)+1; else end end h=h+1; endend附录2：问题一的完整数据结果四个城市的综合评价指标得分完整结果（取小数点后五位）为：ABCD2009/6/10.03403 0.00064 0.09808 0.34854 2009/6/20.00276 0.10305 0.23795 0.23921 2009/6/30.00039 0.00020 0.05998 0.06132 2009/6/40.05217 0.12363 0.24007 0.66614 2009/6/50.05878 0.11685 0.16207 0.38856 2009/6/60.00061 0.00010 0.00000 0.04958 2009/6/70.00031 0.00003 0.00882 0.05714 2009/6/80.00004 0.00024 0.05339 0.06882 2009/6/90.00002 0.00001 0.00027 0.00704 2009/6/100.00000 0.00000 0.00156 0.00100 2009/6/110.00003 0.00061 0.03163 0.12279 2009/6/120.00064 0.00054 0.00057 0.05521 2009/6/130.00044 0.00144 0.02172 0.06898 2009/6/140.00011 0.04811 0.10128 0.11553 2009/6/150.05116 0.43618 0.18570 0.20942 2009/6/160.34321 0.49469 0.50670 0.58250 2009/6/170.10857 0.06252 0.57531 0.30074 2009/6/180.03632 0.05183 0.24770 0.05998 2009/6/190.02065 0.00000 0.19665 0.04202 2009/6/200.33182 0.01540 0.16726 0.09518 2009/6/210.11603 0.04535 0.01590 0.25513 2009/6/220.00061 0.00003 0.08252 0.07344 2009/6/230.04746 0.00004 0.01114 0.11512 2009/6/240.17774 0.06186 0.16145 0.71678 2009/6/250.39631 0.27731 0.40810 0.74042 2009/6/260.62218 0.40253 0.16286 0.47511 2009/6/270.03529 0.01596 0.20446 0.08952 2009/6/280.13965 0.02077 0.24596 0.26102 2009/6/290.26530 0.00130 0.00636 0.10280 2009/6/300.00001 0.00000 0.00081 0.00002 2009/7/10.01978 0.03610 0.00057 0.15256 2009/7/20.04181 0.01954 0.00456 0.01610 2009/7/30.16545 0.12568 0.39322 0.09036 2009/7/40.22075 0.62081 0.36871 0.60661 2009/7/50.50593 0.01380 0.16867 0.05067 2009/7/60.16225 0.26703 0.04510 0.00526 2009/7/70.36332 0.05960 0.10230 0.17295 2009/7/80.03751 0.00020 0.07232 0.08562 2009/7/90.00125 0.00028 0.08389 0.05371 2009/7/100.14931 0.19347 0.46092 0.72179 2009/7/110.00276 0.00002 0.00000 0.08471 2009/7/120.00820 0.13261 0.25892 0.07600 2009/7/130.06068 0.42639 0.16508 0.49390 2009/7/140.00509 0.17135 0.33876 0.06782 2009/7/150.19731 0.09581 0.08418 0.33270 2009/7/160.09633 0.00085 0.14921 0.19348 2009/7/170.01771 0.00064 0.04774 0.09399 2009/7/180.01282 0.00080 0.00199 0.04096 2009/7/190.13893 0.13753 0.12046 0.10999 2009/7/200.00125 0.00019 0.04904 0.03002 2009/7/210.19754 0.20714 0.39681 0.41284 2009/7/220.17641 0.05260 0.14126 0.30287 2009/7/230.36866 0.10004 0.20234 0.16576 2009/7/240.00503 0.01414 0.00664 0.01162 2009/7/250.00039 0.01734 0.00652 0.00379 附录3：求解问题二中灰色关联度分析的程序A=0.0370.0340.0270.040.0350.0270.0260.0260.0270.0250.0360.0280.0320.0350.0370.0440.0290.0330.0310.0450.0320.0330.0330.0560.0590.0650.0510.0440.0660.0350.0520.0450.0570.0760.0730.0560.070.0490.0340.0370.0260.0230.0340.0220.0320.0330.0240.0290.0240.0220.030.0410.0680.0240.0220.0310.0270.0250.0350.030.0210.0230.0230.0190.0180.0240.0170.0250.0190.0270.0310.0290.0260.0230.0230.0220.0210.0240.0350.0410.0490.030.0210.0240.0180.0230.0250.0340.0430.0360.0280.0310.0290.0260.040.0270.0290.0310.0280.0430.0350.0310.030.0520.0260.0430.0370.0360.0280.0250.0980.0670.0390.0670.0760.0560.0470.0610.0510.0380.0610.0550.050.0490.0740.1280.1020.0850.0950.1290.1050.0680.0750.10.1010.1480.0680.1010.0620.0530.0380.030.0630.0930.0890.0770.0510.0610.0410.0960.0580.0650.0790.050.0630.0780.0570.040.1090.0520.0680.0750.060.0410.037663.73660.13668.62665.63664.15669.17672.15671.5663.92666.15665.3664.25666.82669.17668.67673.52671.75675.26663.89662.38659.02662.81661.79663.45658.96666.58668.88657.65654.79659.69663.71664.6664.33671.08668.92667.85666.42667.63665.21663.5668.06671.65663.