《2011数学建模大赛A题完整论文.pdf》由会员分享，可在线阅读，更多相关《2011数学建模大赛A题完整论文.pdf（19页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、20112011高教社杯全国大学生数学建模竞赛高教社杯全国大学生数学建模竞赛承承诺诺书书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料） ，必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从 A/B/C/D 中选择一项
2、填写） ：A我们的参赛报名号为（如果赛区设置报名号的话） ：所属学校（请填写完整的全名） ：电子科技大学参赛队员 (打印并签名) ：1.吴建辉2.赖邱袭3.谭仲伟指导教师或指导教师组负责人(打印并签名)：毛勇日期： 2011年 9月 12日赛区评阅编号（由赛区组委会评阅前进行编号）：20112011高教社杯全国大学生数学建模竞赛高教社杯全国大学生数学建模竞赛编编 号号 专专 用用 页页评阅人评分备注赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：基于城市表层土壤重金属污染分析
3、基于城市表层土壤重金属污染分析摘要摘要关键词：关键词：1一、问题重述一、问题重述随着城市经济的快速发展和城市人口的不断增加， 人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价， 研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为 1 类区、2 类区、5 类区，不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进行调查。为此， 将所考察的城区划分为间距1 公里左右的网格子区域，按照每平方公里 1 个采样点对表
4、层土（010 厘米深度）进行取样、编号，并用 GPS 记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面， 按照 2 公里的间距在那些远离人群及工业活动的自然区取样， 将其作为该城区表层土壤中元素的背景值。附件 1 列出了采样点的位置、海拔高度及其所属功能区等信息，附件 2 列出了 8 种主要重金属元素在采样点处的浓度， 附件 3 列出了 8 种主要重金属元素的背景值。现要求你们通过数学建模来完成以下任务：(1) 给出 8 种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。(2) 通过数据分析，说明重金属污染的主要原因。(3)
5、 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二、问题分析二、问题分析随着城市经济的快速发展和城市人口的不断增加， 人类活动对城市环境质量的影响日显突出。对于重金属在环境中的影响更为明显，因此研究城市表层土壤重金属污染是迫在眉下的事。对于问题一，经过对数据和题目的分析，直接使用matlab（附录一）使用二次插值法可以画出8种主要重金属元素在该城区的空间分布图。对于第二小问，我们首先根据所给的数据将已区分好的各个区域归在一起， 求出各种重金属元素在该区域的平
6、均值，建立综合污染指数评价法模型；分析各种重金属元素在各个区域的污染指数来分析该城区内不同区域重金属的污染程度。对于问题二，由所给的数据求出每个区域的各种金属的平均值，采用主成分分析法对其对进行综合评价分析。 通过主成分的主要元素因子和土壤中的重金属2元素的载荷矩阵和各个因子的累积贡献率以及其百分比来具体分析各种重金属元素对环境污染的主要原因。对于问题三，根据题目和有关资料对重金属污染物的传播特征的分析，可将8种重金属污染物分为两类。一类是主要在大气中传播，一类是主要在土壤中传播。 对于在大气中传播的重金属污染物，我们建立重金属污染物在气体中扩散模型， 根据所在的空间任意位置土壤表面的重金属污
7、染物浓度的多少来确立污染源的位置，函数的最大值即为污染源的位置；同理，对于另一类，我们也可以建立重金属污染物在土壤中扩散模型来确立污染源的位置，只是所考虑的方面不一样，从而所涉及的式子有所差别而已。三、模型假设三、模型假设（1）假设所给的采样点的数据都是可靠准确的；（2）假设每个样本点都能很好的反映该平方公里的实际情况；（3）污染物的排放瞬时完成，且排放速率恒定；（4）风向为水平方向，且风速和风向不随时间的变化而变化；（5）气体的传播服从扩散定律，即单位时间通过单位扩展面积的流量与它的浓度梯度成正比；（6）污染物的沉降速率恒定；四、符号说明四、符号说明Pi：土壤重金属i的污染指数；Ci： 土壤
8、重金属i的实测值；Ri：土壤重金属i的背景值；xi：第i种指标变量；xij：第i个评价对象的第j个指标的取值；rij：第i个评价对象的第j个指标的相关系数；bj：主成分yj的信息贡献率；ap：主成分yp的累积贡献率；s：污染物的扩散速度；k：风速，单位 m/s；H:污染源O距有效地面的高度；t：任意扩散时刻；3Cx, y,z,t：空间任意一点扩散污染物的浓度；vs：空气向土壤的沉降系数；wd：地面干沉结率系数；：冲洗系数；kii x, y,z：空间任意一点扩散污染物的扩散系数；Q0:从污染源放出的污染物总量；Q1:在t,t t内通过某空间域的流量；Q2:在某空间域污染物的增量；x,y,z：分别
9、为用浓度标准差表示的x,y,z轴上的扩散参数；p0：为单位体积内污染物排放的速率；五、模型建立及其求解五、模型建立及其求解5.1 模型一经过对数据和题目的分析，直接使用matlab（附录一）使用二次插值法可以画出8种主要重金属元素在该城区的空间分布图， 由图形来反映各种的分布情况，具体图形如下；4注：红 o 代表生活区，绿*代表工业区，蓝*代表山区，红*代表交通道路区，蓝 o 代表公园区由图形可以得出： 某种重金属元素对某个区域具有非常显著的影响，因此我们再次利用软件matlab（附录二） 求出各种重金属在各个区域最大的浓度和所在的区域，具体的信息可以在下表体现：表一 元素分布元素最大值编号所
10、在区域(g / g)AsCd30.131619.858495交通区（4）交通区（4）Cr920.8422交通区（4）Cu2528.58工业区（2）Hg160009交通区（4）Ni142.522交通区（4）Pb472.4816生活区（1）Zn3760.861交通区（4）由表格可以知道，Pb主要分布在生活区，Cu元素主要分布在工业区，其他元素主要分布在交通区。对于第二小问，我们首先根据所给的数据将已区分好的各个区域归在一起，求出各种重金属元素在该区域的平均值，建立综合污染指数评价法模型；分析各种重金属元素在各个区域的污染指数来分析该城区内不同区域重金属的污染程度。综合污染指数采用国内相关研究常用的
11、内梅罗（Nemerow）污染指数来评价土壤重金属的复合污染状况， 综合污染指数评价可用来评价土壤中某一重金属元素富集和污染程度，而污染程度一般用污染指数来表示，计算见公式：PiCi/Ri（1）综合污染评价公式：2 1nP max P inii1P （2）22单因子指数评价标准为：Pi1，无富集（无污染） ；1 Pi 2，轻度富集（轻度污染） ；2 P；3 Pi，过度富集（严重污染）。i3，中富集（中度污染）由matlab（附录三、四）软件可以得到各种重金属元素在各个区域的污染指数。表二 污染指数区域生活区工业区山区主干道路区公园绿地区元素PiAsCdCrCuHgNiPbZnP1.74182.0
12、1431.12341.58561.73992.23053.02391.17172.76932.15802.22641.72291.25681.87271.40763.74279.66181.31194.71332.28732.658318.35301.170212.76643.28551.49121.61071.25641.43231.24312.22923.00131.17922.04951.95833.43494.02791.06223.51962.23543.170613.53301.25319.4262.7344通过以上表格分析可知，生活区为严重污染，工业区为严重污染，山区为轻6度污染
13、接近无污染，主干道路区为严重污染，公园绿地区为中度污染。5.2 模型二由所给的数据求出每个区域的各种金属的平均值， 采用主成分分析法对其对进行综合评价分析。5.2.1 模型的建立1）对原数据进行标准化处理假设进行主成分分析的指标变量为8种重金属元素：x1,x2,xp，评价对象可以为五个区域，第i个评价对象的第j个指标的取值为xij。将各个指标xij转换成标准化指标xij，xij xjxij，(i 1,2, p)（3）sj1n1n其中xjxij,sj即xj,sj为第j个指标的样本xij xj,j 1,2, p，ni1n1i12均值和样本标准差。对应的，称xij xjx ij，(i 1,2, p)
14、sj为标准化指标变量。2）计算相关系数矩阵R相关系数矩阵R rijpp,i, j 1,2, p（4）n1式中rij1,rij rji,rij是第i个指标与第j个指标的相关系数。rijk1xnkikjx3）计算特征值和特征向量计算相关系数矩阵R的特征值12p 0，及对应的特征向量u1,u2,up，其中uju1j,u2 j,unp,由特征向量组成m个新的指标变量T1u xy1u11x212un1xn1u xy2u12x222un2xn（5）y ux1m1u2mx2unmxnm式中y1是第一主成分，y2是第二主成分，ym是第m主成分。74）选择p(p m)个主成分，计算综合评价值a、计算特征值jj
15、1,2,m的信息贡献率和累积贡献率。称bjjk1mj 1,2,m（6）k为主成分yj的信息贡献率；pk1k1mpk（7）k为主成分y1, y2, yp的累积贡献率，当p接近于1（p 0.85,0.90,0.95）时，则选择前p个指标变量y1, y2, yp作为p个主成分，代替原来m个指标变量，从而可对p个主成分进行综合分析。b、计算综合得分：Z bjyj（8）j1p其中bj为第j个主成分的信息贡献率，根据综合得分就可进行评价。5.2.2 模型的求解经过对土壤单点样重金属元素含量的数据标准化处理后， 利用matlab软件进行主成分分析，可得出如下结果；表三 变量相关矩阵指标 AsCdCrCuHg
16、NiPbZnAsCdCrCuHgNiPbZn1.000.739-.296.473.731.319.762.943.100.225.595.423.892.830.809.656.442.590.646.242.292.205.501.585.7391.000-.546.473.319.100.892-.296-.5461.000-.780-.842-.015-.523.731-.7801.000.762-.842.225-.015.830-.523.595.809.9431.000.656.4231.000.4421.000.590.646.242.292.205.501.5851.000再利
17、用SPSS统计软件可得到前几个因子的累积贡献率和初始特征值， 如表四所示；8表四 解释的总方差初始特征值因子合计123456784.8281.6761.1090.3683.913E-82.643E-92.194E-9累 积 贡 献方差百分比% 值%合计60.35520.95613.8624.8284.891E-73.304E-82.743E-860.35581.31095.172100.000100.000100.000100.0004.8281.6761.109提取平方和载入方 差 百 分 累积贡 献比%值%60.35520.95613.86260.35581.31095.1726.468E
18、-188.085E-17100.000由表四可见，前三个主成分分别累计提取了总方差的95.172%，表明所提取的主成分能够较好地代表源数据所蕴涵的信息。将表四的前三个主成分的主要元素因子和土壤中的重金属元素的载荷矩阵反映在表五中；表五 载荷矩阵因子元素AsCdCrCuHgNiPb1.730.909-.660.928.861.495.9422.368.104.736-.272-.422.274.1763-.508-.245.146.245.187.806-.144Zn.549.794.070从表五中可见，Cd、Cu和Pb在主成分1中具有非常高的载荷，Hg和As在主成分1中的载荷尽次于前三种。研究
19、证实Pb和Cd在来源上关系较密切，并且Pb主要来源于汽车尾气排放而汽车轮胎磨损可以释放Cd， 再加上Cu可以来源于汽油和车体的磨损，从而可以推断出因子1可能主要代表汽车交通运输活动的影响。其次Hg污染不仅与燃煤污染有关，在很大的程度上还与汽车尾气的排放有关；As和S一样在空气中以AsO2为主，也可以是由汽车燃料燃烧所释放的，进一步的加大了推断的可能性。由表格可以看出，第一主成分解释的原变量方差最多（60.355%） ， 而其后的主成分所解释的方差逐渐在减少，因此可以认为第一主成分是污染的主要来源，鉴于此， 可认为主干道路区上的尾气排放和轮胎磨损是重金属Cd、Cu、Pb、Hg和As污染的主要原因
20、。在表五中体现了，Cr和Zn在主成分2中具有较高的载荷。有学者研究表明，9厂矿企业的三废排放中含有大量的Zn和Cr，因此可以断定因子2为工业活动的影响。与此我们可以认为工业区的工业废气排放是重金属Zn和Cr的主要原因。在表五中也体现了，Ni单独在主成分3中具有较高的载荷，且与其他元素关系不密切。由资料显示，废水中含有较多的Ni，因此也可以认为生活区的生活废水、污水是重金属Ni的主要原因。综合以上所述， 可知Ni使生活区的污染变严重，Cr和Zn使工业区的污染变严重，Cd、Cu、Pb、Hg和As使主干道路区的污染变严重。与第一问的结果一致，验证了两个模型的准确性。5.3 模型三根据题目和有关资料对
21、重金属污染物的传播特征的分析，可将8种重金属污染物分为两类。一类是主要在大气中传播，一类是主要在土壤中传播。对于在大气中传播的重金属污染物， 我们建立重金属污染物在气体中扩散模型来确立污染源的位置； 而对于另一类，我们则建立重金属污染物在土壤中扩散模型来确立污染源的位置。1）重金属污染物在气体中的扩散模型记污染物开始扩散的时刻为t 0， 并且以污染源作为扩散中心， 则可将t 时刻无穷空间中任一点（x,y,z）的污染物浓度记为 C(x,y,z,t)。假设单位时间内通过单位法向量面积的流量与浓度梯度成正比，则有 q kigradC其中kii x, y,z是扩散系数，grad表示浓度梯度，其中的负号
22、代表重金属粉尘的浓度是由高到低的地方扩散。假设空间域的体积为 V，包围空间域的曲面为一规则的球面，设其表xy面面积为S表,外法线向量为 n ,1，zz则在t,t t内流通过空间域的流量可表示为：Q1ttt qnddt（9）s空间域所包围的区域内污染的重金属的增量可表示为：Q2Cx, y,z,t tCx, y,z,tdV（10）V由污染源排放出的该重金属的总量为：Q0由质量守恒可得出：Q0Q1Q2,即，tttp dVdt（11）010Cx, y,z,t tCx, y,z,tdV Vttt ttqnddt （12）p0dVdst 根据曲面积分Gauss公式得：qnddivqdV（13）sV则式子（
23、12）可以转换成ttttCx, y,z,t tCx, y,z,t（14）tdV divqdVdt p0dVdttttVVCx, y,z,t tCx, y,zClimlimt由于t0tt0tttkdivgradCdtt（15）C故式子（14）即可转换成dVtdivqdVtp0dVt（16）tVVC化简得：dV divqdV P0 (17)tVV根据A.Fick扩散微分方程式中：CCCC2V2V2VUXUYUZ kXkYkZ（18）222tXYZXYZ其中：C为气体浓度；t为时间；ux,uy,uy为x,y,z方向风速；kx,ky,kz为x,y,z方向上的扩散系数。假设重金属污染物在无风条件下扩散，
24、此时有，C2V2V2V kx2ky2kz2 (19)txyz结合式子（18）可解出（17）的结果为：Cx, y,z,tz Hyxexp (20 )0.51.54kzt4tkxkykz4kxt4kytp0由此得到地面土壤表层的浓度：Cx, y,z,tvswdCx,y,z,tx2 y2 z2 s2t2（21）由（21）式可知，土壤表层的浓度是一个向四周逐渐减少的函数。2）重金属污染物在土壤中的扩散模型记污染物开始扩散的时刻为t 0， 并且以污染源作为扩散中心， 则可将t 时刻无穷平面中任一点（x,y）的污染物浓度记为 C(x,y,t)。假设单位时间内通过单位路程的流量与浓度梯度成正比，则有 q k
25、igradC11其中kii x, y是扩散系数，grad表示浓度梯度，其中的负号代表重金属污染物的浓度是由高到低的地方扩散。x假设在平面S上，设其传播路径为l,法线向量为n ,1，y则在t,t t内流通过平面S的流量可表示为：Q1ttt qndsdt（9）l平面S所包围的区域面积内污染的重金属的增量可表示为：Q2Cx, y,t tCx, y,td（10）s在有水气共存的土壤中,即土壤非饱和带中,污染物随水运移时,水中的污染物会与气体及土壤颗粒发生物质交换。 水气两相间的物质交换平衡可用亨利定律来表示:Ca HCw (22)式中，Ca为气体中的污染物浓度mg cm3,Cw为水中的污染物浓度mg
26、cm3,H为亨利系数。不同物质的H值可在文献中查到。对于水中污染物与固体颗粒间的物质交换,经常用到以下这一简单的线性方程:qs KdCw (23)式中，qs为固体颗粒表面的污染物浓度mg g1,Kd为污染物在水与固体颗粒间的分配系数cm3g1。对于土壤中的水-固物质交换，可以使用Freundlich方程。它是一个经常被使用的非线性吸附等温线方程,其表达式为：1/nqs KFCw(24)KF为Freundlich分配系数mg11/ng1cm3/n，1/ n表示水-固吸附非线性指数。假设气相为静止,且不考虑污染物在气相中的扩散,污染物在非饱和带中随水运移可由下列一维对流扩散方程来表示:CwqsCa
27、Cw2Cwa wvxwDww (25)2tttxx式中，w，a分别为土壤含水率与含气率，为土壤干容重g cm3，vx为孔隙水流速度cmmin1，x为平面距离cm，Dw为扩散系数cm2min1。12综上式子可得：RtCwCw2Cw (26) vx Dw2txx式中，Rt为延滞系数，即Rt1KdW (27)aHaH延滞系数Rt代表污染物运移时由于发生了水气及水固体颗粒之间的物质交换而产生的延滞作用。则土壤中污染物的总浓度：Q2CwCaqs由污染源排放出的该重金属的总量为：Q0由质量守恒可得出：Q0Q1Q2,即，Cx, y,t tCx, y,tdstttttp dsdt（11）0lt ttqndsd
28、t p0dsdt（12）ltl 根据曲面积分Gauss公式得：qnds divqd（13）ls则式子（12）可以转换成ttttCx, y,t tCx, y,tt ddivqdp0dsdt（14）tttsslCx, y,t tCx, yClimlimt由于t0tt0tttkdivgradCdtt（15）C故式子（14）即可转换成dt divqdt p0dst（16）tsslC化简得：ddivqd P0 (17)tss根据A.Fick扩散微分方程式中：CCCC2V2V2VUXUYUZ kXkYkZ（18）tXYZX2Y2Z2其中：C为气体浓度；t为时间；ux,uy,uy为x,y,z方向风速；kx,
29、ky,kz为x,y,z方向上的扩散系数。假设重金属污染物在无风条件下扩散，此时有，C2V2V2V kx2ky2kz2 (19)txyz13结合式子（18）可解出（17）的结果为：z HyxCx, y,z,texp (20 )0.51.54k t4k t4k txyz4tkxkykzp0由此得到地面土壤表层的浓度：Cx, y,z,tvswdCx,y,z,tx2 y2 z2 s2t2（21）由（21）式可知，土壤表层的浓度是一个向四周逐渐减少的函数。六、模型的评价及改进六、模型的评价及改进七、参考文献七、参考文献1 杨苏才，曾静静，王胜利，南忠仁，兰州市表层土壤 Cu、 Zn、 Pb 污染评价及成
30、因分析，市场周刊理论研究第，2004.11。2 吴邵华，周生路，潘贤章，赵其国，城市扩建过程对土壤重金属积累影响的定量分析，土壤学报，2011.5。3 刘丽琼，魏世江，江韬，三峡库区消落带土壤重金属分析特征及潜在风险评价，中国环境科学，2011，31（17） ：1204-1211。4 彭胜，陈家军，王红旗，挥发性有机污染物在土壤中的运移机制与模型，土壤学报，第 38 卷第 3 期 2001 年 8 月5 王雄军， 赖健青等，基于因子分析法研究太原市土壤重金属污染的主要来源2008.17（2）:671-676。6 放射气体模型的预估模型7 姜启源，谢金星，叶俊，数学模型（第三版） ，高等教育出版
31、社，2003,1-12。8 王沫然，MATLAB 与科学计算（第二版） ，电子工业出版社，2004,89-90。9 王冬琳，数学建模及实验，国防工业出版社，2004,197-199八、附录附录一：plot(xm(1),ym(1),ow) holdfor i=1:319switch gongneng(i)case 1 plot(xm(i),ym(i), or)case 2 plot(xm(i),ym(i), *g)case 3 plot(xm(i),ym(i), *b)14case 4 plot(xm(i),ym(i), *r)case 5 plot(xm(i),ym(i), ob)endend
32、load data.txtload data2.txtx=data(:,1);y=data(:,2);z=data2(:,8);xi=0:100:max(x);yi=min(y):100:max(y);zi=griddata(x,y,z,xi,yi,cubic);contour(xi,yi,zi)title(Zn的空间分布图)附录二：for i=1:8; a,b=max(data(:i);End附录三：%划分区域for i=1:size(qu,1)if qu(i)=1 shenghuo=shenghuo;i;endif qu(i)=2 gongye=gongye;i;endif qu(i)=3
33、 shanqu=shanqu;i;endif qu(i)=4 daolu=daolu;i;endif qu(i)=5 gongyuan=gongyuan;i;endendshenghuo,gongye,shanqu,daolu,gongyuan,15结果：1）生活区为以下采样点：12 13 16 18 20 21 23 25 33 36 39 4263 67 68 72 94 106 107 152 154 156 157 158176 180 183 184 186 187 195 198 199 243 254 258259 267 268 269 272 273 276 3062）工业区
34、为以下采样点：4 6 8 29 30 31 103 153 177 185 197 221223 224 225 226 228 230 231 232 233 237 238 239246 247 248 249 255 262 270 271 274 275 277 2783）山区为以下采样点：62 76 80 81 82 92 100 101 104 108 109110 111 112 114 115 116 117 119 120 121 122 123124 125 126 127 128 129 131 132 133 134 135 136137 138 189 194 207
35、 209 210 211 217 218 219 220251 280 281 282 286 287 288 290 291 295 296 297298 299 300 301 302 304 3054）主干道路区为以下采样点：1 2 3 5 7 9 10 14 15 17 19 2226 27 28 32 34 35 37 38 40 41 43 4445 46 47 48 49 50 51 52 53 54 55 5657 58 59 60 61 64 65 66 71 73 75 7778 79 83 84 85 86 87 88 89 90 91 9395 96 97 98 99
36、102 105 113 118 130 150 151155 159 160 161 162 163 164 165 166 167 168 169170 171 172 173 174 175 178 179 181 182 188 191192 196 201 203 204 205 206 208 212 214 215 216222 227 229 234 235 236 240 241 242 244 245 250252 256 257 260 261 263 264 265 266 279 283 284285 289 292 293 303 3185）公园绿地区为以下采样点：1
37、1 24 69 70 74 139 140 141 142 143 144 145146 147 148 149 190 193 200 202 213 253 294 307308 309 310 311 312 313 314 315 316 317 319附录四：load data.txtqu=data(:,4);shenghuo=;gongye=;shanqu=;daolu=;gongyuan=;for i=1:size(qu,1)if qu(i)=1 shenghuo=shenghuo;i;end16if qu(i)=2 gongye=gongye;i;endif qu(i)=3 s
38、hanqu=shanqu;i;endif qu(i)=4 daolu=daolu;i;endif qu(i)=5 gongyuan=gongyuan;i;endendshenghuo=shenghuo,gongye=gongye,shanqu=shanqu,daolu=daolu,gongyuan=gongyuan,load data2.txtload data3.txtpingjunzhi=data3(:,1);fangcha=data3(:,2);shenghuo=data2(shenghuo,:);a1=std(shenghuo)shenghuo=mean(shenghuo)./ping
39、junzhi;gongye=data2(gongye,:);a2=std(gongye)gongye=mean(gongye)./pingjunzhi;shanqu=data2(shanqu,:);a3=std(shanqu)shanqu=mean(shanqu)./pingjunzhi;daolu=data2(daolu,:);a4=std(daolu)daolu=mean(daolu)./pingjunzhi;gongyuan=data2(gongyuan,:);a5=std(gongyuan)gongyuan=mean(gongyuan)./pingjunzhi;pi=shenghuo,gongye,shanqu,daolu,gongyuanp=sqrt(mean(pi).2+(max(pi).2)./sqrt(2)17