matlab物元可拓模型程序,基于物元可拓模型的水资源承载力分析

当前,在水资源短缺日益严重的情况下,进行水资源承载力研究具有重要意义。水资源承载能力是指某一地区的水资源,在一定社会历史和科学发展阶段,在不破坏社会和生态系统时,最大可承载的农业、工业、城市规模和人口的能力,是一个随着社会、经济、科学技术发展而变化的综合目标[1]。目前,研究水资源承载力的方法可分为两类,一类是基于模拟和优化的多目标决策法可分为两类,一类是综合评价法,如模糊综合评判、主成分分析、投影寻踪法等。优化方法模型求解极其困难,而综合评价法计算简单,具有很强的操作性。但水资源承载力研究较为复杂,有时单项指标间相互影响、相互制约,评价结果往往不相容。物元分析理论主要研究事物的可拓性,可以解决矛盾问题[2],其原理是以物元为基元建立物元模型,以物元可拓为依据,应用物元变换化矛盾问题为相容问题。应用熵值法确定权重系数可以避免主观性问题[3]。因此,笔者将模糊物元与熵权结合起来,构建物元可拓模型,对焦作市水资源承载力进行了分析。1物元可拓模型给定事物名称M,它关于特征c的量值为v,把有序三元组R=(M,c,v)作为描述事物的基本元,简称物元。如果事物M有n个特征c1、c2、…、cn,相应的量值为v1、v2、…、vn,称R为n维物元[4],简记为R=(M,c,v),即R=Mc1v1c2v2cnvn(1)1.1经典域与节域研究水资源承载力的量化指标为n项,记为c1、c2、…、cn,量化指标的等级为m项,则经典域物元可表示为Rj=Mjc1v1jc2v2jcnvnj=Mjc1(a1j,b1j)c2(a2j,b2j)cn(anj,bnj)(2)式中:Mj为事物的第j个评价类别;vij为M关于ci的取值范围,vij=(aij,bij)(i=1,2,…,n);(aij,bij)即为经典域。评价指标的节域物元指各指标的所有取值范围,即Rp=Mc1v1pc2v2pcnvnp=Mc1(a1p,b1p)c2(a2p,b2p)cn(anp,bnp)(3)式中:M为事物的全体;vnp为所有等级指标的最大取值范围,vnp=(anp,bnp);(anp,bnp)即为节域,且(aij,bij)(anp,bnp)。1.2关联度及等级的确定(1)对经典域Rj中的元素进行标准化处理,即将经典域的量值除以vp最右端点的取值[5]:Rj=(M,c,v)=Mc1(a1j/b1p,b1j/b1p)c2(a2j/b2p,b2j/b2p)cn(anj/bnp,bnj/bnp)(4)(2)将待评物元事物Rq=qc1v1c2v2cnvn也进行标准化处(5)-a理,得到新的待评物元体,即Rq=qc1(v1/b1p)c2(v2/b2p)cn(vn/bnp)(3)定义关联度:+aDj(vi)=(vi,vij)=|vi-b2|-b2(6)Dj(q)=1-ni=1WiDj(vi)(7)式中:Dj(vi)为待评物元体与各个经典域的距;vi为标准化后的实际点值;(a,b)为标准化后经典域的取值范围;Dj(q)为待评事物q关于各类别等级j的关联度,即待评对象的等级贴近度[5],若Dj0=max[Dj(q)],则评价对象属于等级j0。1.3权重的确定m根据信息熵的定义[6],Hj=-ki=1jfilnjfi,说明信息熵是系统无序程度的度量,信息是系统有序程度的度量。可以根据各指标值的变异程度,利用信息熵计算出各指标的权重。(1)将各指标进行标准化处理后,计算第i个样本下第j项指标的比重:jfi=rijmi=1rij(8)(2)计算第j项指标的熵值:mHj=-ki=1jfilnfji(9)假定fji=0时,jfilnjf