数学模型论文

模型论文 空气污染监测问题 摘 要 本文是一个如何有效对空气污染进行监测，并以此监测结果作为事实依据来解决空气污染的问题。以基本的代数、集合、图论的知识等数学工具建立了层次分析模型、模糊聚类分析模型。本文利用全市均匀分布的多个监测点定时、定点监测大气中二氧化硫、氮氧化物和漂尘的含量，对城市空气污染状况进行科学的量化归类，使人们对环境污染程度有了更为直观的了解。 针对问题(1)，通过对本问题的分析，提出两种解决此类问题的可行性数学模型——层次分析模型和模糊聚类分析模型。通过对两种模型在解决此类问题中的可行性对比，认为层次分析模型在解决本文问题中更具有优势。由于影响空气质量的三大因素污染状况的不确定性及监测数据的无规律性，层次分析模型可以对数据量化并能对其简便分类。 针对问题(2)，本文应用问题(1)提出的层次分析模型，首先通过分析国家环境质量标准规定并依据国家空气污染指数得出空气污染状况划分指标。再对监测点监测到空气污染三大因素的数据进行了量化处理。根据上述提出的划分指标及处理数据，有效的对十八个监测点进行了归类划分。监测点9、13、16属于 、2、3、6、7、8、10、14、18属于一般污染地区，监严重污染地区，监测点1 测点4、5、11、12、15、17属于基本没有污染地区 针对问题(3)，利用问题(2)得出的空气污染状况划分指标及数据处理方法，首先对三个采集点污染元素数据进行处理，依据划分指标，对其进行了归属分类。采集点19、21属于一般污染地区，采集点20属于基本没有污染地区。 针对问题(4)，在解决上述三个问题的过程中，利用层次分析模型及划分指标在对空气污染状况分类中，发现题中所给空气污染因素不够全面及监测点数目不足。因此有必要设计问卷调查表格，适当增加必要的空气质量评价因素，提高评价结果的准确性及划分的细致性，从而达到反映和评价空气质量的最佳效果。 关键词:层次分析模型 模糊聚类分析模型 划分指标 一、 问题重述 1(1 问题背景 改革开放以来随着工业生产和经济建设的高速发展，汽车数量的迅速增加，煤炭和石油等燃料消耗量的剧增，导致了空气污染的加剧。我们都知道当前世界的能源结构中，煤炭占主导地位，这就使得空气污染物表现为二氧化硫和漂尘即可吸入颗粒物污染比例较大。空气中含有多种氮氧化物，其主要来源于化石燃料的燃烧以及生产、使用硝酸的过程，它们排放的一氧化氮和进一步形成的二氧化氮直接严重影响到空气的质量。 二氧化硫是造成空气污染的主要物质之一，我国已将二氧化硫列为一种主要的法规控制空气污染物，并将大气中二氧化硫的浓度水平作为评价空气质量的一项重要指标。二氧化硫对人体健康的危害表现为它易溶于人体的体液和其他黏性液中，其长期的影响会导致多种疾病，甚至在氧化剂、光的作用下，会生成使人致病、甚至增加病人死亡率的硫酸盐气溶胶;二氧化硫对植物的危害表现为在低浓度二氧化硫的情况下，植物的生长机能受到影响，造成产量下降，品质变坏，在高浓度二氧化硫的情况下，植物产生急性危害或直接使植物叶片枯萎脱落;二 氧化硫对金属的腐蚀主要表现在对钢结构的腐蚀，据统计，金属腐蚀造成的直接损失远大于水灾、风灾、火灾、地震造成损失的总和，且金属腐蚀直接威胁到工业设施、生活设施和交通设施的安全;二氧化硫对生态环境的影响表现为其形成的酸雨和酸雾对湖泊、地下水、建筑物、森林、古文物以及人的衣物构成腐蚀。 通常把粒径在10微米以下的颗粒物称为PM10，又称为可吸入颗粒物或漂尘。由于漂尘粒径小，能被人直接吸入呼吸道造成危害，又由于它能在大气中长期漂浮，易将污染物带到很远的地方，导致污染范围扩大，同时在大气中还可为化学反应提供反应床，因此，漂尘是从事环境科学工作者所注目的研究对象之一。国 -1996)》中将漂尘家环保总局1996年颁布修订的《环境空气质量标准(GB3095改称为可吸入颗粒物，作为正式大气环境质量标准。漂尘通常来自未铺沥青、水泥路面上行驶的机动车、材料的破碎碾磨处理过程以及被风扬起的尘土。漂尘被人吸入后，会累积在呼吸系统中，引发许多疾病。粗颗粒物的暴露可侵害呼吸系统，诱发哮喘病;细颗粒物可能引发心脏病、肺病、呼吸道疾病，降低肺功能等。 大气中的氮氧化物主要来源于化石燃料的燃烧和植物体的焚烧以及农田土壤和动物排泄物中含氮化合物的转化，随着这些来源的增加，氮氧化物的浓度将会迅速升高。氮氧化物的危害是多方面的，包括对人体健康、生态系统和建筑设施都有直接和潜在的危害。它可刺激肺部，产生使人较难抵抗感冒之类的呼吸系统疾病，同时其在空气中的水结合最终转化成的硝酸是酸雨的成因之一，它与其他污染物在一定条件下能产生化学烟雾污染。 由上述研究结果得出，这三种空气污染物对农作物、森林、水产及陆地动物都有严重危害，是空气污染物中主要因素。它们对全球大气环境的影响目前已明显表现在三个方面:臭氧层的消耗、酸雨、全球变暖，这些问题如不及时控制将对整个地球造成灾难性的危害。所以，深入研究空气污染监测问题对于人类社会的发展有着至关重要的作用。 空气污染指数(AIR POLLUTION INDEX，简称API)是一种反映和评价空气质量的方法，就是将常规监测的集中空气污染物的浓度简化成为单一的概念性数值形式，并分级表征空气质量与空气污染的程度，其结果简明直观，使用方便，适用于表示城市的短期空气质量状况和变化趋势。所以在本题中，我们将空气污染指数作为衡量监测污染的标准。 1.2 目标任务 问题一:提出对该问题分类的一般数学方法，即通过分析表格1所给数据结合题目中对空气污染评价指标，寻找适合解决此问题的数学方法并建立合理的数学模型。 问题二:假设将城市污染状况分为严重污染、一般污染和基本没有污染三类地区，试建立数学模型，对18个监测点的污染程度进行排序以及分类。 问题三:根据表二数据，按照问题二中所解决的数学方法对19,20,21三个地区污染程度进行归属划分，给出解决这类问题的一般方法。 问题四:自行设计出一个量化评价公民对环境的关心程度的问卷表格，通过该表格，要达到能够体现设计目的的要求。 一、 模型假设 1) 在这段时间内，该地区不发生沙尘暴，不会影响污染物的测量; 2) 市内化工厂能够严格按照国家规定指标排放化学物质; 3) 温度、湿度等其他因素不会对测试仪造成影响; 4) 忽略一氧化碳、臭氧等空气污染物; 5) 各监测点测量标准统一; 6) 问卷调查人员能够积极配合，依据事实填写表格; 二、 符号说明 符号 意义 AP 空气污染指数(评估空气质量状况的一组数字) I W 表示二氧化硫的权值 1 W 表示氮氧化物的权值 2 W 表示漂尘的权值 3 C 表示二氧化硫对空气的污染作用 1 C 表示氮氧化物对空气的污染作用 2 C 表示漂尘对空气污染的作用 3 A 表示成对比较矩阵 a表示Ci/Cj(其中i=1,2,3;j=1,2,3)ij Q 表示污染程度 D 表示18个监测点的数据矩阵(D=D，D,D,„,D) 12318 D表示18个监测点(i=1,2,3，„,18),D=(d，d，d) i ii1i2i3 d代表二氧化硫的含量 i1 d代表氮氧化物的含量 i2 d代表漂尘的含量 i3 I表示污染指数 n ’Q 采集点19，20,21对应的污染程度矩阵 ’D 采集点19,20,21的数据矩阵 ’q采集点19的污染程度值 1 ’q采集点20的污染程度值 2 ‘q采集点21的污染程度值 3 四、模型建立与求解 4.1 问题一 对于问题一，给出对空气污染监测问题进行分类的一般数学方法，经研究对比发现，层次分析法和模糊聚类分析法可以解决这个问题。 4.1.1模型1:层次分析模型 4.1.1.1模型的概述 层次分析法(The analytic hierarchy process)简称AHP，在20世纪70年代中期由美国运筹学家托马斯?塞蒂(T.L.Saaty)正式提出。它是一种定性和定量相结合的、系统化、层次化的分析方法。由于它在处理复杂的决策问题上的实用性和有效性，很快在世界范围得到重视。它的应用已遍及经济和管理、能源政策和分配、行为科学、军事指挥、运输、农业、教育、人才、医疗和环境等领域。 层次分析法的基本思路与人对一个复杂的决策问题的思维、判断过程大体上是一样的。不妨用假期旅游为例:假如有3个旅游胜地A、B、C供你选择，你会根据诸如景色、费用和居住、饮食、旅途条件等一些准则去反复比较这3个候选地点(首先，你会确定这些准则在你的心目中各占多大比重，如果你经济宽绰、醉心旅游，自然分别看重景色条件，而平素俭朴或手头拮据的人则会优先考虑费用，中老年旅游者还会对居住、饮食等条件寄以较大关注。其次，你会就每一个准则将3个地点进行对比，譬如A景色最好，B次之;B费用最低，C次之;C居住等条件较好等等。最后，你要将这两个层次的比较判断进行综合，在A、B、C中确定哪个作为最佳地点。 4.1.1.2层次分析法的基本步骤 1、建立层次结构模型。在深入分析实际问题的基础上，将有关的各个因素按照不同属性自上而下地分解成若干层次，同一层的诸因素从属于上一层的因素或对上层因素有影响，同时又支配下一层的因素或受到下层因素的作用。最上层为目标层，通常只有1个因素，最下层通常为或对象层，中间可以有一个或几个层次，通常为准则或指标层。当准则过多时(譬如多于9个)应进一步分解出子准则层。 2、构造成对比较阵。从层次结构模型的第2层开始，对于从属于(或影响)上一层每个因素的同一层诸因素，用成对比较法和1—9比较尺度构追成对比较阵，直到最下层。 比较第 i 个元素与第 j 个元素相对上一层某个因素的重要性时，使用数量化的相对权重a来描述。设共有 n 个元素参与比较，则称为成ij 对比较矩阵。 成对比较矩阵中a的取值可参考 Satty 的提议，按下述标度进行赋值。aijij在 1-9 及其倒数中间取值。 a = 1元素 i 与元素 j 对上一层次因素的重要性相同; ij a = 3元素 i 比元素 j 略重要; ij a = 5元素 i 比元素 j 重要; ij a = 7 元素 i 比元素 j 重要得多; ij a = 9元素 i 比元素 j 的极其重要; ij a = 2n，n=1,2,3,4元素 i 与 j 的重要性介于a = 2n ? 1与a = 2n + 1ijijij 之间; ，n=1,2,...,9 当且仅当a = n。 ij 成对比较矩阵的特点: 。 3、计算权向量并做一致性检验。对于每一个成对比较阵计算最大特征根及对应特征向量，利用一致性指标、随机一致性指标和一致性比率做一致性检验。 若检验通过，特征向量(归一化后)即为权向量:若不通过，需重新构追成对比较阵。 从理论上分析得到:如果A是完全一致的成对比较矩阵，应该有 aa = a。 ijjkik 但实际上在构造成对比较矩阵时要求满足上述众多等式是不可能的。因此退而要求成对比较矩阵有一定的一致性，即可以允许成对比较矩阵存在一定程度的不一致性。 由分析可知，对完全一致的成对比较矩阵，其绝对值最大的特征值等于该矩阵的维数。对成对比较矩阵 的一致性要求，转化为要求: 的绝对值最大的特征值和该矩阵的维数相差不大。 检验成对比较矩阵 A 一致性的步骤如下: , 计算衡量一个成对比矩阵 A (n>1 阶方阵)不一致程度的指标CI: 其中λ是矩阵 A 的最大特征值。 max 从有关资料查出检验成对比较矩阵 A 一致性的标准RI:RI称为平均随机一致性指标，它只与矩阵阶数 有关。 按下面公式计算成对比较阵 A 的随机一致性比率 CR: 。 判断方法如下: 当CR<0.1时，判定成对比较阵 A 具有满意的一致性，或 其不一致程度是可以接受的;否则就调整成对比较矩阵 A，直到达到满意的 一致性为止。 4、计算组合权向量并做组合一致性检验。计算最下层对目标的组合权向量，并根据公式做组合一致性检验，若检验通过，则可按照组合权向量表示的结果进行决策，否则需要重新考虑模型或重新构造那些一致性比率较大的成对比较阵。 4.1.2模型2:模糊聚类分析模型 4.1.2.1模型的概述 在科学技术、经济管理中常常要按一定的标准(相似程度或亲疏关系)进行分类。例如，根据生物的某些性状可对生物分类，根据土壤的性质可对土壤分类等。对所研究的事物按一定标准进行分类的数学方法称为聚类分析，它是多元统计“物以类聚”的一种分类方法。由于科学技术、经济管理中的分类界限往往不分明，因此采用模糊聚类方法通常比较符合实际。 4.1.2.2模糊聚类分析的一般步骤 [9] 1、第一步:数据标准化 (1) 数据矩阵 Uxxx,{,,,}m 设论域为被分类对象，每个对象又有个指标表示其性状，12n 即 xxxx,{,,,} ， (1,2,,)in,iiiim12 于是，得到原始数据矩阵为 xxx,,11121m,,xxx21222m,,。 ,, ,,xxx,,nnnm12 xnm其中表示第个分类对象的第个指标的原始数据。 nm (2) 数据标准化 在实际问题中，不同的数据一般有不同的量纲，为了使不同的量纲也能进行比较，通常需要对数据做适当的变换。但是，即使这样，得到的数据也不一定在 [0,1]区间上。因此，这里说的数据标准化，就是要根据模糊矩阵的要求，将数据 [0,1]压缩到区间上。通常有以下几种变换: ? 平移标准差变换 xx,ikk,x,iks(1,2,,;1,2,,)inkm,,k nn112sxx,,()xx,,kikk,kiknn,1ii,1其中 ， 。 经过变换后，每个变量的均值为0，标准差为1，且消除了量纲的影响。但 ,x[0,1]是，再用得到的还不一定在区间上。 ik ? 平移极差变换 ,,xx,min{}ikik1,,in,,x, ， (1,2,,)km,ik,,max{}min{}xx,ikik1,,in1,,in ,,01,,x显然有，而且也消除了量纲的影响。 ik ? 对数变换 ,xx,lg (1,2,,;1,2,,)inkm,,ikik 取对数以缩小变量间的数量级。 2、第二步:标定(建立模糊相似矩阵) Uxxx,{,,,}xxxx,{,,,}12niiiim12 设论域，，依照传统聚类方法确定相似 xrRxx,(,)rRxx,(,)xijijijijji系数，建立模糊相似矩阵，与的相似程度。确定的 方法主要借用传统聚类的相似系数法、距离法以及其他方法。具体用什么方法， 可根据问题的性质，选取下列公式之一计算。 (1) 相似系数法 ? 夹角余弦法 m xx,ikjkk,1 。 ,rijmm22xx,,ikjkkk,,11 ? 最大最小法 m ()xx,,ikjkk,1r, 。 ijm ()xx,,ikjkk,1 ? 算术平均最小法 m 2()xx,,ikjkk,1r, 。 ijm ()xx，,ikjkk,1 ? 几何平均最小法 m 2()xx,,ikjkk,1r, 。 ijm xx,ikjkk,1 x,0以上3种方法中要求，否则也要做适当变换。 ij ? 数量积法 1,,ij,,,m， r,1,ijxxij,,,,ikjk,M,1k, m Mxx,max()其中 。 ,ikjkij,k,1? 相关系数法 m xxxx,,,ikijkj,1k ， r,ijmm22()()xxxx,,,,ikijkj,,11kk mm11xx,xx,其中 ，。 ,,jjkiikmmkk,1,1? 指数相似系数法 2m()xx,13ikjk ， r,,exp[],ij2ms4,k1k n12sxx,,()其中 ， ,kikikni1, n1xx,,kikkm,(1,2,,)ni,1而 。 (2) 距离法 ? 直接距离法 rcdxx,,1(,)ijij ， 01,,rdxx(,)ijijc其中为适当选取的参数，使得，表示他们之间的距离。经常 用的距离有: 海明距离 m dxxxx(,),,,ijikjk,1k 欧几里得距离 m2 dxxxx(,)(),,,ijikjk,1k 切比雪夫距离 m 。 dxxxx(,),,,ijikjk,1k ? 倒数距离法 1,,ij,, ,M 。 r,,ij,,ij,,dxx(,)ij, 01,,r其中为适当选取的参数，使得。 Mij? 指数距离法 rdxx,,exp[(,)]ijij 。 3、第三步:聚类(求动态聚类图) (1)基于模糊等价矩阵聚类方法 传递闭包法 ? *RR 根据标定所得的模糊矩阵还要将其改造称模糊等价矩阵。用二次方法 *tR(),RR求的传递闭包，即=。再让由大变小，就可形成动态聚类图。 ? 布尔矩阵法 布尔矩阵法的理论依据是下面的定理: Uxxx,{,,,} 定理4.1.2.2 设是上的一个相似的布尔矩阵，则具有传RR12n递性(当是等价布尔矩阵时)矩阵在任一排列下的矩阵都没有形如RR,11111001,,,,,,,,的特殊子矩阵。 ,,,,,,,,,,,10011111,,,,,,,, 布尔矩阵法的具体步骤如下: R ? 求模糊相似矩阵的截矩阵. ,,, RR ? 若按定理4.1.2.2判定为等价的，则由可得在水平上的分类，U,,,RR若判定为不等价，则在某一排列下有上述形式的特殊子矩阵，此时只要将,, 其中特殊子矩阵的0一律改成1直到不再产生上述形式的子矩阵即可。如此得到 **RR的为等价矩阵。因此，由可得水平上的分类 ,,, (2) 直接聚类法 所谓直接聚类法，是指在建立模糊相似矩阵之后，不去求传递闭包，也tR()不用布尔矩阵法，而是直接从模糊相似矩阵出发求得聚类图。其步骤如下: ,,1x[]x ? 取(最大值)，对每个作相似类，且 1iiR {|1}xr,[]x =， jijiR r,1xx即将满足的与放在一类，构成相似类。相似类与等价类的不同之处是，ijji 不同的相似类可能有公共元素，即可出现 []{,}xxx,[][]xx,,,[]{,}xxx, ,,。 iRjkijiRik ,,1此时只要将有公共元素的相似类合并，即可得水平上的等价分类。 1 (,)xx,, ? 取为次大值，从中直接找出相似度为的元素对(即Rij22 r,,x,,1x)，将对应于的等价分类中所在的类与所在的类合并，将所有的ij2j1i ,这些情况合并后，即得到对应于的等价分类。 2 (,)xx,, ? 取为第三大值，从中直接找出相似度为的元素对(即Rij33 r,,x,x)，将对应于的等价分类中所在的类与所在的类合并，将所有的这ij3j2i ,些情况合并后，即得到对应于的等价分类。 3 ? 以此类推，直到合并到成为一类为止。 U ,4.1.2.3最佳阈值的确定 在模糊聚类分析中对于各个不同的,,[0,1]，可得到不同的分类，许多实际问题需要选择某个阈值，确定样本的一个具体分类，这就提出了如何确定阈值, 的问题。一般有以下两个方法: , ? 按实际需要，在动态聚类图中，调整的值以得到适当的分类，而不需, 要事先准确地估计好样本应分成几类。当然，也可由具有丰富经验的专家结合专业知识确定阈值，从而得出在水平上的等价分类 ,, ? 用F统计量确定最佳值, Uxxx,{,,,}nxm 设论域为样本空间(样本总数为)，而每个样本有个12ni xxxx,{,,,}特征:，。于是得到原始数据矩阵，如下表所示，(1,2,,)in,iiiim12 n1xxkm,,(1,2,,)其中，称为总体样本的中心向量。 x,kikn,1i 样 指 标 本 1 2 k m xxxxx111211km1 xxxxx212222km2 xxxxxiiikim12i xxxxxnnnknm12n x()xxxx12km rn设对应于值的分类数为，第类的样本数为，第类的样本记为: jj,j ()()()()jjjj()j()()()jjjxxxx,(,,,)xxxx,,,，第类的聚类中心为向量，其中为第jk112m12nj 个特征的平均值，即 k nj1()()jj (1,2,,)km,，， xx,,kikn,1ij 作统计量 F r()jnxxr,,(1),jj,1 ， F,njr()()jjxxnr,,(),,i,,11ji m()()2jj其中 xxxx,,,(),kk,1k ()()jj()j()j()jxx,为x与间的距离，为第j类中第个样本x与其中心x间的距xii nr,离。称为统计量，它是遵从自由度为，的分布。它的分子表征类Fr,1F 与类之间的距离，分母表征类内样本间的距离。因此，值越大，说明类与类之F 间的距离越大;类与类间的差异越大，分类就越好。 4.1.3问题的综合分析与进一步研究 4.1.3.1层次分析法的优点 运用层次分析法有很多优点，其中最重要的一点就是简单明了。层次分析法不仅适用于存在不确定性和主观信息的情况，还允许以合乎逻辑的方式运用经验、洞察力和直觉。也许层次分析法最大的优点是提出了层次本身，它使得买方能够认真地考虑和衡量指标的相对重要性。 4.1.3.2模糊聚类分析法的优点 聚类分析是将事物根据一定的特征，并按某种特定要求或规律分类的方法。由于聚类分析的对象必定是尚未分类的群体，而且现实的分类问题往往带有模糊性，对带有模糊特征的事物进行聚类分析，分类过程中不是仅仅考虑事物之间有无关系，而是考虑事物之间关系的深浅程度，显然用模糊数学的方法处理更为自然，因此称为模糊聚类分析。 4.1.3.3层次分析法与模糊聚类分析法的比较 通过层次分析法与模糊聚类分析法的优点分析与比较，针对空气污染监测问题，层次分析法对空气污染的三大因素污染程度的不确定行和主观成分允许以合乎逻辑的方式进行量化运用，使得问题的运算结果并不局限于本文，更能发挥其优势。 4.2 问题二 对于问题二，采用层次分析法，依据表1中的数据给出18个监测点的污染程度的排序以及分类结果。 4.2.1模型:层次分析模型 1、 构造空气污染问题的层次结构 将空气污染问题分解为2个层次，最上层为目标层，即污染程度，最下层为 准则层，有二氧化硫w、氮氧化物w、漂尘w三个准则。 123 结构图如下图所示: 污染程度 目标层 准则层 氮氧化物w漂 尘w二氧化硫w1 2 3 2、构造成对比较矩阵 参考附录1中国家环境质量标准规定，假设当二氧化硫日平均浓度为0.15毫克/立方米、氮氧化物日平均浓度为0.10毫克/立方米、漂尘日平均浓度为0.3毫克/立方米时，三种污染因素对空气的污染作用(或者说对环境造成的危害)是一样的。则在含量相同的情况下，二氧化硫对空气的污染作用(C)、氮氧化1物对空气的污染作用(,)、漂尘对空气的污染作用(,)三者之间的比值为23 ,:,:,,2:3:1(可得 12 3 1 2/3 2 A= 3/2 1 3 (1) 1/2 1/3 1 其中，a表示C/C。 ijij 3、确定权向量 一般地，如果一个正互反阵A满足 a*a=a i,j,k=1,2,3,…,n (2) ijjkik 则A称为一致性矩阵，简称一致阵. (1)式给出的A显然是一致阵。 如果得到的成对比较阵是一致阵，像(1)式的A，自然应取对应于特征根的，归一化的特征根向量(即分量之和为1)表示诸因素O的权重，这个向量称为权向量。如果成对比较阵A不是一致阵，但在不一致的容许范围内(下面将说明如何确定这个范围)Saaty 等人建议用对应于A最大特征根(记作λ)的特征向量(归一化后)作为权向量w,即w满足 A *w=λ* w (3) 由于本题中的成对比较阵A是一致阵，经演算可得A的特征根A: 1 11/3 A= 11/2 ( 4 ) 1 11/6 因为11/3+11/2+11/6=11,经归一化的特征向量w作为权向量: 11/3 1/3 w= 11/2 * 1/11 = 1/2 (5) 11/6 1/6 即有w=1/3 , w=1/2 , w=1/6 123 4、模型的一致性检验 前面已经给出n阶一致性的特征根是n，n阶正互反阵A的最大特征根λ?n，而当λ=n时，A是一致阵。 根据这个定理和λ连续地依赖于a的事实可知，λ比n大得越多，A的不ij 一致程度越严重，用特征向量作为权向量引起的判断误差越大。因而可以用λ- n数值的大小来衡量A的不一致程度。T.L.Saaty将 CI=(λ- n)/(n –1) (6) 定义为一致性指标。CI=0时A为一致阵;CI越大A的不一致程度越严重。 为了确定A的不一致程度的容许范围，需要找出衡量A的一致性指标CI的标准。Satty又引入随机一致性指标RI(见表5)。 对于n?3的成对比较阵A，将它的一致性指标CI与同阶(指n相同)的随机一致性指标RI之比称为一致性比率CR，当 CR=CI/RI,0.1 (7) 时认为A的不一致程度在容许范围之内，可用其特征向量作为权向量。(上式中0.1的选取是带有一定主观信度的。) 由于本题所得成对比较阵A是一致性矩阵，一致性程度已达最高标准，并不需要进行一致性检验，这里仅作介绍，不再引用具体事例演算。 另外，对于组合一致性检验同样如此，具体演算方法见4.1.1.2中第4部分。 5、污染程度确定 设污染程度用Q表示，18个监测点的数据矩阵为D=(D,D,D,…,D)，对12318 于每个监测点有D=(d,d,d)。其中，d,d,d分别代表D二氧化硫、氮氧化物ii1i2i3i1i2i3i 和漂尘的含量。则有 0.045 0.066 0.094 0.003 0.048 0.210 0.086 0.196 0.187 D= 0.043 0.039 0.061 0.003 0.015 0.066 0.072 0.072 0.082 0.265 0.264 0.194 0.102 0.106 0.263 0.274 0.211 0.301 0.053 0.020 0.035 0.205 0.088 0.058 0.225 0.002 0.091 0.060 0.008 0.015 0.068 0.058 0.025 0.091 0.001 0.070 (8) 0.209 0.112 0.170 0.284 0.215 0.147 0.342 0.101 0.199 0.045 0.066 0.094 0.003 0.048 0.210 0.086 TQ= w)D= 1/3 1/2 1/6 0.043 0.039 0.061 0.003 0.015 0.066 0.072 0.265 0.264 0.194 0.102 0.106 0.263 0.274 0.196 0.187 0.053 0.020 0.035 0.205 0.088 0.058 0.225 0.002 0.091 0.072 0.082 0.060 0.008 0.015 0.068 0.058 0.025 0.091 0.001 0.070 0.211 0.301 0.209 0.112 0.170 0.284 0.215 0.147 0.342 0.101 0.199 = ( 0.081 0.087 0.094 0.020 0.041 0.147 0.110 0.137 0.154 0.083 0.029 0.048 0.150 0.094 0.056 0.178 0.018 0.098 ) (9) 6、分类结果确定 污染情况分为严重污染、一般污染和基本没有污染三类。分类依据就是污染 程度Q的矩阵: 严重污染 Q?I 2 分类结果= 一般污染 I,Q?I (10) 21 基本没有污染 I,Q?0 1 采用附录2中的国家空气污染指数，对应本题中的18个监测点，设 基本没有污染为API : 0~50 一般污染为API : 50~100 为API : 100以上 严重污染 由表3 污染指数 污染物浓度 I TSP SO NO n2X 100 0.300 0.150 0.100 50 0.120 0.050 0.050 可得 0.050 0.150 0.050 0.150 TI I =w* 0.050 0.100 = 1/3 1/2 1/6 0.050 0.100 12 0.120 0.300 0.120 0.300 = 0.062 0.150 (11) 即分类条件设置如下: 严重污染 Q?0.150 分类结果 = 一般污染 0.150,Q?0.062 (12) 基本没有污染 0.062,Q?0 用以上条件对18个监测点进行排序分类，结果如下 表5: 4.3 问题三 对于问题三，采集到19,20,21三个地区的污染数据(参见表2)，确定分类，并给出这类问题的一般方法。 4.3.1模型:层次分析法 ’ 1、设定给定三个采集点的污染元素数据矩阵为D,则有 0.101 0.045 0.015 ’ D= 0.052 0.005 0.011 (13) 0.181 0.122 0.323 ’用(7)式的方法，求出对应的污染程度矩阵，设为Q，则有 TQ’=w)D’= 0.090 0.038 0.064 (14) 2、采样点污染状况分类结果 由q’=0.090 , 0.150>q’>0.062 , 故采集点19属于一般污染; 11 由q’=0.038 , q’<0.062 , 故采集点20属于基本没有污染; 22 由q’=0.064, 0.150>q’>0.062 , 故采集点21属于一般污染。 33 4.3.2解决这类问题的一般方法步骤: 1、将决策问题分解为目标层和准则层2~3个层次。 2、通过相互比较确定各准则对于目标的权重，及各个方案对于每一准则的权重。 3、在层次分析法中将准则层对目标层的权重进行综合计算，并计算出成对比较矩阵和权向量。 4、对于每一个成对比较阵计算最大特征根及对应特征向量，利用一致性指标，随机一致性指标和一致性比率做一致性检验。 4.3.3问题的综合分析与进一步研究 层次分析法从处理问题的类型来看，主要是决策、评价、分析、预测等。由于它在处理复杂的决策问题上的实用性和有效性，很快在世界范围内得到普遍的重视和广泛的应用，应用范围已遍及经济计划和管理、能源政策和分配、行为科学、军事指挥、运输、农业、教育、人才、医疗、环境等等领域。 从上面介绍的层次分析法的基本步骤看，建立层次结构模型是关键一步。假如给定空气污染的主要因素除了二氧化硫、氮氧化物和漂尘以外，还包括其他因素(比如一氧化碳、臭氧等)，那么，就必须重新给出相应的成对比较阵和全向量的计算模式，其他具体分类计算方法则不变。 建立层次模型的层次结构也不可能只有固定的二到三层，以管理信息系统为 例:当今任何部门每天都会接触到大量的信息，信息管理水平的高低直接关系着工作效率，甚至生存条件、财务、库存、销售、行政……各种各样的管理信息系统开发完成或准备推广时，通常要作全面的检查、测试和分析。这一类的信息管理系统指标层次结构，可以形成目标层、准则层、子准则层和方案层。 将问题包含的因素分层:最高层(解决问题的目的);中间层(实现总目标而采取的各种措施、必须考虑的准则等。也可称策略层、约束层、准则层等);最低层(用于解决问题的各种措施、方案等)。把各种所要考虑的因素放在适当的层次内。用层次结构图清晰地表达这些因素的关系。 此外，构造成对比较阵是整个工作的数据依据，当然是重要的，应当由经验和知识丰富、判断力强的专家给出，或者由专家给出数据作为凭据进行计算，还不妨采用群体判断的方式。 至于第三、第四步的计算工作，数学工作者是容易完成的。 另外，由于本次空气污染监测问题是依据国家标准，通过模型建立得出的分类排序计算方法，这种分类排序结果已不仅仅局限于本题目18个测试点，而是可以推广于其他城市与地区的、接近国家标准的综合分类标准。 4.4问题四 对于问题四，设计问卷表格，旨在根据表格量化评价公民对环境的关心程度。 在对问题2、3的分析研究过程中，发现对空气质量进行评价时，只考虑到二氧化硫、氮氧化物、漂尘这三个因素。由此可知，在模型应用过程中，因素偏少，降低了空气质量评价的标准。因此，在问卷调查中，应当增加必要的其他因素，提高空气质量评价结果的准确性。 设计问卷调查表如下: 以下设计的题中，1,16题为单选题，17,20题为多选题，请根据实际情况认真填写问卷，感谢您的积极支持和参与。 1、对于自然环境您如何看待 ? ( ) A、我们是自然之子，我热爱并敬畏大自然，非常关心 B、比较重视 C、我只关心我身边的环境，其他地方的环境我并不关心 D、完全不关心 2、对环境保护的迫切程度，您如何看待,( ) A、环境保护是一件很紧迫的事，人人应该参与进来 B、应先把生活水平提高到一定程度，再谈环境保护 C、我们的环境还没到非要刻意去保护的地步 D、确实很紧迫，但那是国家的事 3、认为现在哪种环境破坏对你的生活和工作影响最大,( ) A、空气污染 B、水源污染 C、噪声污染 D、固体垃圾污染 4、您是怎样看待空气污染的,( ) A、空气污染随时随处都有，是不可避免的，希望政府能够采取相应措施 B、自己积极主动考虑怎样治理污染问题，但不会落到实处 C、空气污染到处存在，只要不危及到自身生命安全就好 D、不关心，短期内不会有什么问题 5、您觉得空气质量如何,( ) A、存在一定的污染，但相信政府可以解决，自己也会参与其中 B、有污染，有机会支持并会适当的奉献自己的力量来解决问题 C、一般，勉强可以生活，所以自己不会去努力解决问题 D、很差，并且相信不会有相关部门来解决 6、在生活区空气受到污染的情况下，您最希望获得以下哪方面的帮助,( ) A、向相关行政部门申诉或向消协投诉 B、找媒体曝光 C、关乎到自己的切身利益时，会考虑向相关部门反映 D、无所谓，对政府部门不抱有希望 7、您最关心下面哪类内容,( ) A、环保类 B、时事政策类 C、财经类 D、娱乐类 8、在同等条件下，您选择居住地时最看重什么,( ) A、周围环境的绿化 B、周边设施的健全程度 C、周边交通的便利程度 D、其他 9、您做过对环境有害的事情吗,( ) A、完全没有 B、比较少 C、不清楚 D、做过，经常并且很随意 10、您会关注空气质量吗,( ) A、经常有意识的关注 B、有时间会关注 C、很少关注 D、不关注，关注也是做无用功 11、您对周边环境空气治理满意吗,( ) A、满意 B、一般 C、不满意 D、没注意 12、你认为政府当下应该花大量经费在环保措施上吗,( ) A、应该，并且应当尽快落实 B、应该，但对政府抱的希望不大 C、不应该，没有什么用处 D、随便，不关心 13、你觉得你所在地区酸雨严重吗,(如树叶枯黄、建筑腐蚀等)( ) A、很严重 B、比较严重 C、一般 D、不存在 14、你认为以下哪项对大气污染最严重,( ) A、煤炭燃烧 B、工业废气 C、交通尾气 D、自然因素 15、您认为国家对环境生态的重视程度怎样,( ) A、非常重视 B、一般 C、不重视 D、不知道 16、您对环境保护公益事业的态度如何,( ) A、很好，我很支持并想成为环境保护的志愿者 B、遇到这样的活动我会关注并支持，但不会成为志愿者 C、志愿者的奉献精神令人感动，但我认为这种方式不会产生实质效果 D、没事找事做，浪费时间和精力，不理解 17、您知道空气污染会带来哪些危害吗,(多选) ( ) A、臭氧层破坏 B、酸雨 C、各种呼吸病和皮肤病 D、影响全球气候 E、动植物生长受到危害 F、您认为还有, 18、您觉得影响空气质量较大的因素有,(多选)( ) 、二氧化硫 A B、氮氧化物 C、漂尘 D、一氧化碳 E、臭氧 F、您认为还有, 19、您觉得当今大气治理有哪些不足,(多选) ( ) A、法律体制不够完善 B、技术落后 C、资金投入不足 D、教育宣传匮乏 E、政策疲软 F、您认为还有, 20、您觉得你能为环保事业做些什么,(多选) ( ) A、不乱扔垃圾 B、节约用水用电 、开发新能源 C D、发明有效的治理污染的方法 E、制止他人污染环境 F、您认为还有, 单项选择题 多项选择题 选项 1 2 3 4 5 6 7 8 9 11111111112题目 0 1 2 3 4 5 6 7 8 9 0 A 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 3 3 3 3 B 4 4 4 3 4 4 3 3 4 4 4 3 4 4 3 3 3 3 3 3 C 2 2 4 2 2 2 3 3 2 2 2 2 4 4 2 2 3 3 3 3 D 1 1 4 1 1 1 3 1 1 1 1 1 4 4 1 1 3 3 3 3 E 3 3 3 3 六、模型的评价 从层次分析法的原理、步骤、应用等方面的讨论不难看出它有以下优点: (1) 系统性 层次分析把研究对象作为一个系统，按照分解、比较判 断、综合的思维方式进行分类决策，成为继机理分析、统计分 析之后发展起来的系统分析的重要工具。 (2) 实用性 层次分析把定性和定量方法结合起来，能处理许多用传 统的最优化技术无法着手的实际问题，应用范围很广。同时， 这种方法将决策者与决策分析者相互沟通，决策者甚至可以直 接应用，这就增加了决策的有效性。 (3) 简洁性 具有中等文化程度的人即可了解层次分析的基本原理 和掌握它的基本步骤，计算也非常简便，并且所得结果简单， 容易为决策者了解和掌握。 层次分析法的局限性可以用囿旧、主观等词来概括。就是说，第一，它只能从原有方案中选优或者分类，不能生成新方案;第二，从建立层次结构模型到给出成对比矩阵，人的主观因素的作用很大，这就使得决策结果可能难以为众人接受。 但是，值得一提的是，本文所建立的空气污染问题层次结构模型，有效地引入国家环境质量标准和国家空气污染指数作为计算判断的方法，在一定程度上克服了这个缺点。 此外，本表格设计的问题符合客观实际情况，能够围绕调查课题和研究假设设计必要的问题 体现出了自然环境与人类生存息息相关的事实，通过此表格能 够客观地评价公民对环境的关心程度。 附录1: 国家环境质量标准规定 1、二氧化硫(SO) 2 二氧化硫主要由燃煤及燃料油等含硫物质燃烧产生，其次是来自自然界，如火山爆发、森林起火等产生。二氧化硫对人体的结膜和上呼吸道粘膜有强烈刺激性，可损伤呼吸器管可致支气管炎、肺炎，甚至肺水肿呼吸麻痹。短期接触二氧化硫浓度为0.5毫克/立方米空气的老年或慢性病人死亡率增高，浓度高于0.25毫克/立方米，可使呼吸道疾病患者病情恶化。长期接触浓度为0.1毫克/立方米空气的人群呼吸系统病症增加。另外，二氧化硫对金属材料、房屋建筑、棉纺化纤织品、皮革纸张等制品容易引起腐蚀，剥落、褪色而损坏。还可使植物叶片变 /立方米。 黄甚至枯死。国家环境质量标准规定，居住区日平均浓度低于0.15毫克2、氮氧化物(NO) x 空气中含氮的氧化物有一氧化二氮(NO)、一氧化氮(NO)、二氧化氮(NO)、22三氧化二氮(NO)等，其中占主要成分的是一氧化氮和二氧化氮，以NOx(氮氧23 化物)表示。NO污染主要来源于生产、生活中所用的煤、石油等燃料燃烧的产x 物(包括汽车及一切内燃机燃烧排放的NOx);其次是来自生产或使用硝酸的工厂排放的尾气。当NO与碳氢化物共存于空气中时，经阳光紫外线照射，发生光x 化学反应，产生一种光化学烟雾，它是一种有毒性的二次污染物。NO比NO的2毒性高4倍，可引起肺损害，甚至造成肺水肿。慢性中毒可致气管、肺病变。吸入NO，可引起变性血红蛋白的形成并对中枢神经系统产生影响。NO对动物的x影响浓度大致为1.0毫克/立方米，对患者的影响浓度大致为0.2毫克/立方米。国家国家环境质量标准规定，居住区的平均浓度低于0.10毫克/立方米。 3、粒子状污染物 总悬浮颗粒物(TSP):总悬浮颗粒物是指粒径在100微米以下的颗粒物，简称TSP。其对人体的危害程度主要决定于自身的粒度大小及化学组成。TSP中粒径大于10微米的物质，几乎都可鼻腔和咽喉所捕集，不进入肺泡。对人体危害最大的是10微米以下的浮游状颗粒物，称为漂尘。可经过呼吸道沉积于肺泡。慢性呼吸道炎症、肺气肿、肺癌的发病与空气颗粒物的污染程度明显相关，当长年接触颗粒物浓度高于0.2毫克/立方米的空气时，其呼吸系统病症增加。国家环境质量标准规定居住区日平均浓度低于0.3毫克/立方米。 附录2: 空气污染指数(API) 空气质量周报就是根据国家《环境空气质量标准》中规定的几种常见污染物例行监测的结果，评价城市一周内的空气质量，并以空气污染指数的表征形式来向公众发布。 空气污染指数(AIR POLLUTION INDEX，简称API)是一种反映和评价空气质量的方法，就是将常规监测的几种空气污染物的浓度简化成为单一的概念性数值形式、并分级表征空气质量状况与空气污染的程度，其结果简明直观，使用方便，适用于表示城市的短期空气质量状况和变化趋势。 空气污染指数的确定原则:空气质量的好坏取决于各种污染物中危害最大的污染物的污染程度。空气污染指数是根据环境空气质量标准和各项污染物对人体健康和生态环境的影响来确定污染指数的分级及相应的污染物浓度限值。目前我国所用的空气指数的分级标准是:(1)空气污染指数(API)50点对应的污染物浓度为国家空气质量日均值一级标准;(2)API100点对应的污染物浓度为国家空气质量日均值二级标准;(3)API200点对应的污染物浓度为国家空气质量日均值三级标准;(4)API更高值段的分级对应于各种污染物对人体健康产生不同影响时的浓度限值，API500点对应于对人体产生严重危害时各项污染物的浓度。 根据我国空气污染的特点和污染防治工作的重点，目前计入空气污染指数的污染物项目暂定为:二氧化硫、氮氧化物和总悬浮颗粒物。随着环境保护工作的深入和监测技术水平的提高，再调整增加其它污染项目，以便更为客观地反应污染状况。 空气污染指数的计算与报告: 污染指数与各项污染物浓度的关系是分段线性函数(见表1和图1)，用内插法计算各污染物的分指数I(具n体计算方法请参见《环境监测简报》1997年第9期)，取各项污染物分指数中最大者代表该区域或城市的污染指数。即: API,max(I,I„I,„I) 12in 该指数所对应的污染物即为该区域或城市的首要污染物。当污染指数API值小于50时，不报告首要污染物。 表3污染指数API分级限值 污染指数 污染物浓度 ITSP SONOn 2 X 500 1.000 2.620 0.940 400 0.875 2.100 0.750 300 0.625 1.600 0.565 200 0.500 0.250 0.150 100 0.300 0.150 0.100 50 0.120 0.050 0.050 表4 空气污染指数及对应的空气质量级别 空气污染指数空气质量级别 空气质量状况 对健康的影响 对应空气(API) 质量的适 用范围 0,50 ? 优 可正常活动 自然保护 区、风景 名胜区和 其他需要 特殊保护 的地区 51,100 ? 良 可正常活动 城镇 中确定的 居住区、 商业交通 居民混合 区、文化 区、一般 工业区和 农村地区 101,200 ? 普遍(轻度污长期接触，易感特定工业 染) 人群症状有轻区 度加剧，健康人 群出现刺激症 状 201,300 ? 不佳(中度污一定时间接触 染) 后，心脏病和肺 病患者症状显 著加剧，运动耐 受力降低，健康 人群中普遍出 现症状 附录3 表5 随机一致性指标RI的数值 n 1 2 3 4 5 6 7 8 9 10 11 RI 0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45 1.49 1.51