水资源短缺风险综合评价 (修复的)

1、水资源短缺风险综合评价摘要：本文通过建立模型来判定北京市水资源短缺风险的主要因子对北京市水资源短缺风险进行综合评价，进而提出调控办法。对于问题一一，影响响水资源源短缺的的因子很很多，主主要有四四方面：第一，农农业用水水；第二二，工业业用水；第三，人人口规模模；第四四，气候候条件与与水利工工程设施施。以上上四方面面分别对对应附表表中农业业用水量量，工业业用水量量，第三三产业及及生活等等其他用用水量与与水资源源总量。对于主要因子，本文采用关联分析对关联度进行计算量化处理。首先对数据进行了预处理，以缺水量（总用水量-水资源总量）作为参考数列，把农业用水量，工业用水量，第三产业及生活等其他用水量与水资

2、源总量作为参考数列，然后对个数列进行初始化处理，利用matlab分别计算出以上四方面对缺水量（总用水量-水资源总量）的相关性。得出总体相关性大小排序如下：0.64777 00.63327 00.59971 0.558444即：水资源源总量第三产产业与生生活等其其他用水水量农农业用水水量工工业用水水量为检验该模模型的合合理性，本文采用matlab作出以上四个量以及缺水量（总用水量-水资源总量）对时间的关系图，从图中可以直观显示农业与缺水量的相关性较大，与该模型结果吻合，模型具有较好的准确性。对于问题二二，本文文建立了了合适的的模型对对北京市市水资源源短缺风风险进行行综合评评价， 作出风风险等级级

3、划分。本本文将改改革开放放以来的的三十年年分成六六个阶段段，每个个阶段分分为五个个点。采采用熵值值确定农农业用水水量，工工业用水水量，第第三产业业及生活活等其他他用水量量三方面面对水资资源短缺缺影响的的权重，得得出水资资源短缺缺的综合合测评指指数Q，再再利用六六个阶段段的Q值值与实际际数据对对比的结结果，定定义出反反映水资资源短缺缺程度的的程度系系数e。由于水水资源总总量相对对于其他他三方面面因子的的特殊性性，本文文决定分分两个阶阶段（分分别是11979920000和和20001220088）拟合合出（水水资源总总量/总总用水量量）的比比值相对对于时间间的一次次函数，根根据函数数走势对对e进行

4、行修正，再对程度系数进行区间划分，作为风险等级的指标。最后计算出：197920000：ee=0.96775 风险等等级为高200120008：ee=0.80113 风险等等级为较高本文采用模模糊集对对模型，以以集对分分析为基基础，重重视信息息的相对对性和模模糊性，综综合的评评价了三三十年来来北京地地区水资资源短缺缺风险情况况。 本本文所采采用的方方法主要要优点是是注重信信息的关关联性，模糊性和相对性，并且能够对结果进行模型检验和对结果进行修正处理。本文结尾还对模型进了评价与推广。关键词：北北京地区区 水水资源短短缺 关联分分析 熵 程度系系数 模糊集集对模型型一 问题重重述水资源，是是指可供供

5、人类直直接利用用，能够够不断更更新的天天然水体体。主要要包括陆陆地上的的地表水水和地下下水。风风险，是是指某一一特定危危险情况况发生的的可能性性和后果果的组合合。水资资源短缺缺风险，泛泛指在特特定的时时空环境境条件下下，由于于来水和和用水两两方面存存在不确确定性，使使区域水水资源系系统发生生供水短短缺的可可能性以以及由此此产生的的损失。近年来，我我国、特特别是北北方地区区水资源源短缺问问题日趋趋严重，水水资源成成为焦点点话题。以以北京市市为例，北北京是世世界上水水资源严严重缺乏乏的大都都市之一一，其人人均水资资源占有有量不足足3000m3，为全全国人均均的1/8，世世界人均均的1/30，属属重

6、度缺缺水地区区，附表表中所列列的数据据给出了了19779年至至20000年北北京市水水资源短短缺的状状况。北北京市水水资源短短缺已经经成为影影响和制制约首都都社会和和经济发发展的主主要因素素。政府府采取了了一系列列措施, 如南南水北调调工程建建设, 建立污污水处理理厂,产产业结构构调整等等。但是是，气候候变化和和经济社社会不断断发展，水水资源短短缺风险险始终存存在。如如何对水水资源风风险的主主要因子子进行识识别，对对风险造造成的危危害等级级进行划划分，对对不同风风险因子子采取相相应的有有效措施施规避风风险或减减少其造造成的危危害，这这对社会会经济的的稳定、可可持续发发展战略略的实施施具有重重要

7、的意意义。北京20009统统计年鉴鉴及市市政统计计资料提提供了北北京市水水资源的的有关信信息。利利用这些些资料和和你自己己可获得得的其他他资料，讨讨论以下下问题：1评价判定定北京市市水资源源短缺风风险的主主要风险险因子是是什么？影响水资源源的因素素很多,例如：气候条条件、水水利工程程设施、工工业污染染、农业业用水、管管理制度度，人口口规模等等。2建立一个个数学模模型对北北京市水水资源短短缺风险险进行综综合评价价， 作作出风险险等级划划分并陈陈述理由由。对主主要风险险因子,如何进进行调控控，使得得风险降降低？二 问题分分析问题一：根据附表中中信息，要要求对北北京市水水资源短短缺风险险的主要要风险

8、因因子进行行判定，因因数分析析的基本本方法有有回归分分析和关关联分析析，但回回归分析析有很多多欠缺，要要求数据据量大，计计算量大大及可能能产生反反常现象象等，故故本文采采用关联联分析。影影响水资资源短缺缺的因子子很多，主主要有四四方面：第一，农农业用水水；第二二，工业业用水；第三，人人口规模模；第四四，气候候条件与与水利工工程设施施。以上上四方面面分别对对应附表表中农业业用水量量，工业业用水量量，第三三产业及及生活等等其他用用水量与与水资源源总量。由由于是在在评价风风险，故故本文在在进行关关联分析析时，将将缺水量量作为参参考数列列，将以以上四方方面因素素作为比比较数列列。先将将各数列列进行初初

9、始化处处理，由由于水资资源总量量越多，风风险越小小，故初初始化处处理时应应区别对对待。通通过建立立关联模模型计算算出各个个数列的的关联度度，从而而进行量量化比较较，容易易得出结结论。为为了保证证结果的的可靠性性，本文文还要做出了了图形进进行模型型检验。问题二：本题主要是是对北京京市水资资源短缺缺风险进进行综合合评价，作作出风险险等级划划分和提提出调控控建议。我我们认为为对北京京市水资资源短缺缺风险的的综合评评价应该该从农业业用水量量，工业业用水量量，第三三产业及及生活等等其他用用水量与与水资源源总量四四方面对对水资源源短缺的的影响出出发，从从客观的的实际数数据出发发，给定定各类影影响因素素各一

10、个个权值反反映其影影响大小小，利用用相对比比较的原原理进行行评测。对对三十年年来的水水资源短短缺的综综合分析析，本文文采用分分为六个个小时间间段，将将每个时时间段划划分为五五个小点点，对个个点的水水资源短短缺情况况进行评评测，再再采用平平均检测测值的方方法判断断一个时时期内的的风险情情况。采采用熵值值确定农农业用水水量，工工业用水水量，第第三产业业及生活活等其他他用水量量三方面面对水资资源短缺缺影响的的权重，得得出水资资源短缺缺的综合合测评指指数Q，再再利用六六个阶段段的Q值值与实际际数据对对比的结结果，定定义出反反映水资资源短缺缺程度的的程度系系数e，由由于水资资源总量量相对于于其他三三方面

11、因因子的特特殊性，本本文决定定分两个个阶段（分分别是11979920000和和20001220088）拟合合出（水水资源总总量/总总用水量量）相对对于时间间的一次次函数，根根据函数数走势对对e进行行修正，再再对程度度系数进进行区间间划分，作作为风险险等级的的指标。最最后计算算出e ，按照照风险等等级划分分的指标标，得出出结论。三 模型假假设假设附件中中统计数据据正确无无误假设计算中中的微小小误差忽忽略不计计假设所提供供的数据据，默认认成涵盖盖各种可可能性，覆覆盖范围围大四 符号说说明 比较较数列XX对参考考数列在在k时刻刻的关联联系数。ri 第i类类用水量量与风险险的关联联度Cij： 归一化化

12、矩阵元元素Xij 实际置置矩阵元元素Hj 第j种种评价指指标的熵熵 W 熵权矩矩阵j 熵权矩矩阵元素素Dj 第j类类污染物物相对于于标准值值的相对对差求和和Q 综合测测评指数数Qq 第q个个月的综综合测评评指数五 模型建建立与求求解51问题题一51.11模型的的建立与与求解首先本文对对题目附附表中数数据进行行预处理理，得出出缺水量量（即总总用水量量-总水水资源量量）的数数据（见见附表）。在进行关联度分析之前，为了使无量钢化和所有数列具有共同点，需要对附表中各个数列进行初始化处理：定义数列xx=（xx（1），xx(2),x（nn），称称为原始数列列X的初初始化数数列。利用以上公公式对表表中农业业

13、用水量量，工业业用水量量，第三三产业及及生活等等其他用用水量，水水资源总总量和缺缺水量进进行初始始化处理理，但由由于水资资源总量量增大意意味着风风险减小小，初始始化处理理时，采采用以下下公式：依照问题的的要求，本本文选取取缺水量量数列作作为参考考数列进进行关联联性分析析。参考数列：比较数列（mm个）则称 为比较数列列X对参参考数列列在k时时刻的关关联系数数。根据关联度度公式求出关联度度。得出：r1=0.59771 r2=0.558444 rr3=00.63327 r44=0.64777（求解程序序见附录录）即：农业用用水量的的关联度度为0.59771，工工业用水水量的关关联度为为0.55844

14、4，第三三产业与与生活等等其他用用水量的的关联度度为0.63227，水水资源总总量的关关联度为为0.664777。5.1.22模型检检验：根据附表中中数据，分分别画出出农业用用水量，工工业用水水量，第第三产业业与生活活等其他他用水量量，水资资源总量量和缺水水量对时时间的曲曲线走势势，如下下：备注：蓝色色为总用用水量；紫色为水资资源总量量；绿色为农业业用水量量；红色为工业业用水量量；淡蓝色为第第三产业业及生活活用水量量；黄色为缺水水量。从表中可以以看出：水资源总量量与缺水水量吻合合度较高高，第三三产业与与生活等等其他用用水量次次之，与与模型结结果一致致，证明明模型具具有较高高的科学学性与合合理性

15、。故根据对应应关系可可以得出出结论：北京市市水资源源短缺风风险的主主要风险险因子是是气候条条件、水水利工程程设施，次次主要因因子为管管理制度度.人口口规模，工工业用水水量对风风险影响响最小。5.2问题题二5.2.11模型的建立与求求解 首先本本文将330年分分成六个个时期，每每个时期期中包含含五个点点。用55个点对对农业用用水量，工工业用水水量，第第三产业业及生活活等其他他用水量量3个方面面用水量量值构成成决策矩矩阵。具体分期如如下表：第一时期第二时期第三时期第四时期第五时期第六时期第一点197919841989199419992004第二点198019851990199520002005第三

16、点198119861991199620012006第四点198219871992199720022007第五点198319881993199820032008 这里类比比信息论论中的熵熵概念。信信息论中中。信息息熵反映映了信息息无序化化的程度度熵越小小信息作作用越大大；熵越越大，信信息作用用越小。通通过度量量评价指指标的效效用大小小，从而而获得对对水质影影响的权权重。而而权重源源于数据据本身（附附表），因因此可以以避免人人为主观观因素判判断而形形成偏差差，从而而可以客客观全面面地从各各点数据据中的得得到对北北京地区区水资源源短缺的的综合评评测指标标。（以以下用nn来描述述各个点点，m来来描述农

17、农业用水水量，工工业用水水量，第第三产业业及生活活等其他他用水量量与水资资源总量量，n=5，mm=4）根根据n个个点，mm个水量量，建立立判断矩矩阵：（Xij）mn（ii=1,2,33n，jj=1,2，3m）根据评价指指标的属属性差异异，可将将评价指指标分为为以下两两种：（1）递增增型（随随评价指指标的递递增，水水资源风风险越大大，如农业业用水量量，工业业用水量量，第三三产业及及生活等等其他用用水量）（2）递减减型（随随评价指指标的递递剪减，水水资源风风险越大大）在计算各指指标的权权重之前前，有必必要每一一个采样样样本进进行归一一化处理理，具体体操作如如下：将判断矩阵阵做归一一化处理理，得到到

18、归一化化后的判判断矩阵阵：Cij=（XXij-XXminn）/(Xmaax-XXminn)式中，在同同种评价价指标下下：Xmmax表表示五个个点中最差的的，即最最有可能能加大风风险者。 Xmmin表表示五个个点中最最好的，即即最有利利于减小小风险者者。例如：当考考虑评价价递增型型指标农农业用水水量时，XXmaxx为各个个样品中中农业用用水量的的的最大大值，XXminn为各个个各个样样品中农农业用水水量的最最小值。 当考考虑评价价递减性性指标水水资源总总量时，XXmaxx为各个个样品中中的最小小值，XXminn为各个个样品中中的最大大值。根据熵的定定义，nn个点mm个评价价指标，可可以确定定评价

19、指指标的熵熵为：Hj=-11/lnn(ijlnnij)在上式中，有有ij=Ciij/(Cij)为使lniij有意意义，一一般需要要假定当当ij=00时，ijlnnij=00。但当当ij=11时ijjlnijj也等于于零，显显然不符符合实际际，与熵熵的含义义相悖，故需要对ij加以修正，将其定义为：=（1+CCij）/ij)矩阵元素的的计算公公式为：j=（1-Hj）/(j)评价指标熵熵权W矩矩阵公式式如下：W=（j）1*mm熵权W具有有如下性性质：j=1下面定义标标准值jj：考虑到北京京地区气气候干旱旱，人口口众多，通过网络查询得出：第一产业，第二产业和第三产业用水量的科学比例大约为3.5:2.5

20、:4，首先算出北京地区三十年年平均水资源总量，然后分别用公式：1=*0.352=*0.253=*0.4计算三个产产业的用用水量标标准值。定义相对差差dijj：dij= Xijj-Ij上式用来描描述农业业用水量量等这类递递增型指指标显然值越大大意味着着i点的的j类评评价指标标对山水水资源短短缺风险险影响越越大，风风险越高高。对相对差求求和：Dj=ijjDj指j类类指标在在各个点点的相对对差之和和，Djj可以整整体反映映j类指指标对风风险的影影响。定义风险的的综合评评测指数数Q：Q=jjQ是从各类类指标所所占权重重与求和和相对差差的角度度全面考考虑每个个时期的的水资源源短缺风风险，通通过进一一步计

21、算算，可以以求出三三十年来来北京地地区水资资源短缺缺风险状状况。定义衡量标标准Qss：由于近几年年来，北北京加大大产业结结构调整整和南水水北调等等水利工工程设施施建设力力度，从从附表中中数据可可以得出出，北京京20001220088年水资资源短缺缺风险已已趋于好好转，因因此我们们算出农农业用水水量，工工业用水水量，第第三产业业与生活活等其他他用水量量八年中中平均比比例为：4:2:44，作为为较为合合理的次次标准比例例，K1=I*0.44K2=I*0.22K3I*00.4而我们可以以算出三三十年年年平均水水资源总总量约为为30亿亿方，故故定义衡衡量标准准为：Qs=*nn*j衡量标准合合理性分分析

22、：我我们希望望的水资资源利用用状况是是按照利利用量不不超过水水资源总总量而且且利用量量比例为为3.55:2.5:44，而上上文中次次标准用用水比例例虽不够够合理，但但利用量量未超过过水资源源总量，作作为此标标准，故故本文采采用他们们的差值值作为衡衡量标准准具有合合理性。定义程度系系数：e=q-ss/q-s上式中，ee表示三三十年内内北京地地区水资资源短缺缺总体风风险程度度； K表示时时期数（本本题k=6）； q表示第第q个时时期的综综合评测测指数；分子的意义义：各个个时期实实际综合合评测指指数与衡衡量标准准的差求求和，因因为每个个时期与与衡量标标准的差差有正有有负，求求和反映映了k个个时期风风

23、险有与与无在量量上的相相对差。分母的意义义：各个个时期实实际综合合评测指指数与衡衡量标准准的差绝绝对值求求和，值值为正，反反映了风风险有无无在量上上的总和和。比值的意义义：水资资源短缺缺风险有有无的相相对优势势。从表达式中中可以看看出：e-1,1,按照我我们的模模型用ee的均匀匀分配把把风险综综合评价价等级定定义为以以下几类类比较合合理：e-1,0根根据上述述模型，基基本上不不存在。e0,11/6定义为为低e1/66,1/3定定义为偏低e1/33,1/2定定义为较低e1/22,2/3定定义为中等e2/33,5/6定定义为较高e5/66,1定义为为高（求解e的的程序见见附录）结果：横坐标代表表第

24、一到到第六时时期；纵坐标表示示Q（即即综合评评价指标标）值，绿绿线为QQs值。e= 0.96775下面分两个个阶段对对e值进进行修正正：根据附表中中数据，利利用maatlaab分别别画出11979920000和和20001220088两个阶阶段的（水水资源总总量/总总用水量量）的值值对时间间的散点点图与拟拟合曲线线，如下下：197920000散点点图：200120008散点点图和拟拟合曲线线：图注：斜率率K= 0.0060888197920000阶段段：散点点图过于于分散，拟合误差将会很大，故e只在此阶段保持不变，即19792000年，北京水资源短缺风险等级为高200120008阶段段：将散散

25、点图拟拟合成一一条直线线，斜率率K= 0.0060888，由由此可知知近年来来水资源源总量占占总用水水量比例例逐渐增增大，意意味着风风险一点点程度上上有所下下降，故故根据下下式：e=e*（1-k）3对e进行修修正，得得e= 0.80113 。即22001120008年年，北京京水资源源短缺风风险等级级为较高。5.2.22主要因因子调控控建议：考虑到北京京水资源源短缺风风险较高高，而北北京市水水资源短短缺风险险的主要要风险因因子是气气候条件件、水利利工程设设施，次次主要因因子为管管理制度度.人口口规模，工工业用水水量对风风险影响响最小。为为使风险险降低，本本文提出出如下调调控建议议：1合理利用用

26、本地水水资源，提提高城市市供水安安全保证证程度水资源短缺缺是北京京市经济济发展与与社会发发展的主主要制约约因素。北北京市水水资源已已远不能能满足未未来用水水需求，为为了提高高城市供供水安全全的保证证程度，本本着“节流优优先，治治污为本本，多渠渠道开源源”的城市市水资源源可持续续利用的的新战略略，鉴于于连续干干旱时地地表水供供水保证证程度降降低，为为保证城城市供水水安全，应应科学地地适度增增加地下下水开采采量，合合理开发发利用。对对已确定定的应急急供水水水源地应应尽快投投入勘探探和开发发工作，对对其它地地区继续续开展调调查工作作，寻找找新的后后备应急急水源。节约用水是是当务之之急也是是长远发发展

27、战略略方针，在在优先保保证城市市生活和和重点工工业供水水的前提提下，在在无法满满足需水水时，适适度压缩缩农业用用水。加加强工业业、农业业节水力力度，调调整产业业结构，大大力发展展节水型型工业、农农业。继续开展污污水资源源化、雨雨洪利用用的研究究和应用用。把城城市污水水排放规规划管理理、污水水处理厂厂建设、再再生污水水利用三三个环节节综合起起来，全全面规划划考虑，实实现污水水资源化化。收集集和利用用城市雨雨洪，既既可防治治雨洪灾灾害，缓缓解城市市雨洪压压力；同同时又增增加了可可用水资资源，并并可通过过回灌补补给蓄养养地下水水。建议统筹安安排官厅厅水库上上游地区区的地表表水开发发利用，加加大官厅厅

28、水库上上游地区区的污水水治理力力度，维维持一定定的入库库水量，并并逐步恢恢复官厅厅水库的的饮用源源水功能能。2应该坚持持水资源源可持续续利用，支支持城市市可持续续发展“南水北调调”中线工工程，是是解决北北京乃至至海河流流域缺水水问题的的战略设设施。南水北北调中线线工程的的实施后后，应建建立外来来水源、本本地水源源相互协协调的供供水网络络，实现现本地地地表水源源与外来来水源的的联合调调蓄、地地下水与与地表水水的联合合调蓄，提提高北京京城市供供水安全全保证程程度，支支持城市市可持续续发展。充分发挥南南水北调调外来水水源的供供水功能能，调整整地下水水开采布布局，减减少地下下水开采采，全面面恢复地地下

29、水超超采区的的生态环环境，养养蓄地下下水，使使地下水水资源可可持续利利用。在永定河、潮潮白河冲冲洪积扇扇的中上上部地区区，利用用洪水、水水库弃水水和地表表水进行行人工回回灌，恢恢复地下下水环境境。3、加强地地下水环环境保护护，建立立完善的的地下水水动态监监测系统统在水资源短短缺的今今天，地地下水作作为北京京的主要要供水水水源，其其水质的的好坏直直接影响响到城市市的发展展，因此此，为从从区域上上保护地地下水，防防止水质质恶化，应应从源头头上即地地下水补补给区进进行保护护。建立立一套完完善的地地下水监监测体系系，逐步步做到自自动化监监测和地地下水水水质、水水位的统统一监测测，资料料要及时时、准确确

30、反映环环境的现现状。动动态监测测工作在在满足向向社会发发布公益益性信息息的同时时，还要要考虑能能够反映映地下水水的动态态变化特特征，为为今后深深入研究究地下水水环境系系统的变变化、演演变规律律提供基基础资料料和数据据。六 模型评评价与推推广问题一中，本本文采用用了关联联分析的的方法讨讨论了农农业用水水量，工工业用水水量，第第三产业业与生活活等其他他用水量量，水资资源总量量和缺水水量的关关联度，摒摒弃了回回归分析析要求数数据量大大，计算算量大及及可能产产生反常常现象等等欠缺，考考察全面面，并进进行了模模型检验验，所建建立的模模型具有有较高的的准确性性与科学学性。问题二中，本本文采用用模糊集集对模

31、型型对北京京地区水水资源短短缺风险险进行综综合评价价，采用用熵的概概念，分分阶段对对e值进进行修正正。是模模型具有有较高的的综合性性和可靠靠性。当然，由于于计算是是忽略了了一些误误差，对对模型的的准确性性有一定定的影响响，仍有有改进的的余地。此模型具有有较多优优点，可可以推广广到许多多其他领领域，对对于事物物的相互互关联性性及对事事物的综综合评价价都有很很多参考考与应用用价值。参考文献1王家家文、王王皓、刘刘海，MMATLLAB77.0编编程基础础 机机械工业业出版社社，20007.72姜启启源，数数学模型型（第第三版），北北京：高高等教育育出版社社，20003.83徐瑞瑞，黄兆兆东，阎阎凤玉

32、MMATLLAB220077科学计计算与工工程分析析，科科学出版版社，220088.94同济济大学数数学系高等数数学（第六六版下册册），高等教教育出版版社，20007.66附录：数据预处理理表格：年份总用水量(亿立方方米)农业用水(亿立方方米)工业用水(亿立方方米)第三产业及及生活等等其它用用水(亿亿立方米米)水资源总量量（亿方方）缺水量（亿亿方）197942.92224.18814.3774.3738.2334.69198050.54431.83313.7774.942624.544198148.11131.612.2114.32424.111198247.22228.81113.8994.

33、5236.610.622198347.56631.611.2444.7234.712.866198440.05521.84414.37764.017739.3110.74198531.71110.12217.24.3938-7.711198636.55519.4669.917.1827.0339.52198730.9559.6814.0117.2638.666-8.299198842.43321.99914.0446.439.1883.25198944.64424.42213.7776.4521.55523.099199142.03322.711.97.4342.299-0.266199246

34、.43319.94415.51110.98822.44423.999199345.22220.35515.2889.5919.67725.555199445.87720.93314.57710.37745.4220.45199544.88819.33313.78811.77730.34414.544199640.01118.95511.7669.345.877-5.866199740.32218.12211.111.122.25518.077199840.43317.39910.84412.237.72.73199941.71118.45510.56612.714.22227.49920004

35、0.416.49910.52213.39916.86623.544200138.917.49.212.319.219.7200234.615.57.511.616.118.5200335.813.88.413.618.417.4200434.613.57.713.421.413.2200534.513.26.814.523.211.3200634.312.86.215.324.59.8200734.812.45.816.623.811200835.112.05.217.934.20.9数据初始化化与关联联性分析析程序load x.ttxt %把原始始数据存存放在纯纯文本文文件x.txtt 中f

36、or ii=1:4x(i,:)=xx(i,:)/x(ii,1); %标准化化数据endfor ii=5x(i,:)=xx(i,1)./x(i,:); %标准准化数据据enddata=x;n=sizze(ddataa,1);ck=daata(1,:);mm1=ssizee(ckk,1);bj=daata(2:nn,:);m22=siize(bj,1);for ii=1:m1for jj=1:m2t(j,:)=bbj(jj,:)-ckk(i,:); enndjc1=mmin(minn(abbs(tt);jjc2=maxx(maax(aabs(t);rho=00.5;ksi=(jc11+rhho*j

37、jc2)./(abss(t)+rhho*jjc2);rt=suum(kksi)/ssizee(kssi,22);r(i,:)=rrt;endrrs,rrindd=ssortt(r,deesceend) %对对关联度度进行排排序r = 00.59971 0.558444 0.63227 00.64477rs =0.64777 00.63327 0.559711 0.58444x.txtt:4.69 24.5424.1110.6212.860.774-6.299.552-7.713.22523.095.226-0.2623.9925.550.44514.54-5.8618.072.77327.49

38、23.54 119.77 188.5 17.4 113.22 111.3 9.88 111 0.9924.18831.8331.628.8131.621.8410.1219.469.66821.9924.4221.7422.719.9420.3520.9319.3318.9518.1217.3918.4516.49 17.415.513.813.513.212.812.412.014.37713.7712.2113.8911.2414.376617.29.99114.0114.0413.7712.3411.915.5115.2814.5713.7811.7611.110.8410.5610.5

39、2 9.227.558.447.776.886.225.885.224.374.9944.334.5524.7724.00174.3397.1187.2266.446.4457.0047.44310.989.55910.3711.779.3311.112.212.713.39 112.0010.813.012.813.413.713.914.738.233262436.634.739.313827.0338.6639.1821.5535.8642.2922.4419.6745.4230.3445.8722.2537.714.2216.8619.216.118.421.423.224.523.8

40、34.2作用水量月月缺水量量对时间间的图形形走势程程序：x1=119799 19980 19881 119822 19983 19884 119855 19986 19887 119888 19989 19990 119911 19992 19993 119944 19995 19996 119977 19998 19999 220000 20001 20002 220033 20004 20005 220066 20007 20008x3=224.11831.8331.628.8131.621.8410.1219.469.66821.9924.4221.7422.719.9420.3520.

41、9319.3318.9518.1217.3918.4516.49 17.415.513.813.513.212.812.412.0 x4=114.337 113.777 112.221 113.889 111.224 114.3376 17.2 99.911 144.011 144.044 133.777 122.344 111.9 15.51 15.28 14.57 13.78 11.76 11.1 110.884 110.556 110.552 99.27.558.447.776.886.225.885.22x2=442.992 550.554 448.111 447.222 447.55

42、6 440.005 331.771 336.555 330.995 442.443 444.664 441.112 442.003 446.443 445.222 445.887 444.888 440.001 440.332 440.443 441.771 440.44 388.934.635.834.634.534.334.835.1x5=44.377 4.94 4.33 4.52 4.772 44.0117 44.399 7.18 7.226 66.4 6.445 77.044 7.43 10.98 9.559 110.337 111.777 99.3 11.1 112.22 122.7

43、 13.39 12.010.813.012.813.413.713.914.7x6=338.223 226 224 336.66 344.7 39.31 38 27.03 38.66 39.18 21.55 35.86 42.29 22.44 19.67 45.42 30.34 45.87 22.25 37.7 114.222 116.886 119.2216.118.421.423.224.523.834.2x7=x22-x66plot(x1,x2,x1,x3,x1,x4,x1,x5,x1,x6,x1,x7)问题二中标标准计算算程序：x6=338.223 226 224 336.66 344

44、.7 39.31 38 27.03 38.66 39.18 21.55 35.86 42.29 22.44 19.67 45.42 30.34 45.87 22.25 37.7 114.222 116.886 119.22 166.118.421.423.224.523.834.2averaage=meaan(xx6)I1=avveraage*0.335;I2=avveraage*0.225;I3=avveraage*0.445;求Q和e的的相关程程序：functtionn c=juddge(x)%输入参数数x为监监测值三三维矩阵阵for kk=1:6 foor ii=1:3 ffor j=11:5 cc(i,j,kk)=(x(ii,j,k)-minn(x(i,:,k)./(mmax(x(ii,:,k)-miin(xx(i,:,kk); eend enndendfuncttionn w=powwersshann(h)%输入为熵熵值矩阵阵，输出出为各评评价指标标熵权矩矩阵for kk=1:6 foor i