2014第七届优秀第二阶段一等奖

编号专用页竞赛统一编号（由竞赛送至评委团前编号竞赛评阅编号（由竞赛评委团评阅前进行编号2 关键 土地储备风险；数据包络分析；模糊聚类分析 要本文针对土地储备风险评价问题，首先，对五部门专家评分进行相关性分析，并剔个指标对风险的隶属度作为输出，通过数据包络分析方法，运用加以实现，得72个项目分类并求所选题目:C3（ 必须部分，选填可加分，加分不超 总分的Aimingattheproblemofriskevaluationoflandreserve,first,onfivesectorexpertsscorewereyzed,andremovesomeerrordata;then,onthebasisofcorrelationysis,theenvelopeysismethodforprocessingthefirststagemodel,establishesthefuzzyevaluationmodelbasedonpacketnetworkysis,inordertoenhancethecredibilityoftheevaluationresults;finally,bythemeansoffuzzyclusteringysis,riskstructureoflandreservesintheprojectadifferentclassification,accordingtothedifferentcategoriesofitems,withtheenvelopeysismethod,todeterminetheriskindexintheircredibilitywithinrange.Aimingattheproblemofareference,expertevaluation,thefuzzyevaluationofthefirststageofthemodelismodified,themodeltoachieveahigherdegreeofconfidence.Expertysisofratingdata,andestablishesafuzzyriskassessmentmodelbasedondata ysismethod.First,selectionofthefivesectorexpertsscoreasinput,capitalstructure,capitalreturnrate,threeindexdynamicrecoverycyclemembershipofriskdegreeastheoutput,throughthemethodofdataenvelopment ysis,usingtoachieve,thefuzzyriskassessmentofthecredibilityandtheinputandoutputvectors,andthecausesofthelowreliabilityofsomeproject;then,inthetargetconditionundertheumreliability,theduallinear ysis,getthecorrectresultsofthreeoutputs;finally,theoutputoftheoutputweightvectorsandcorrectionhasbeenobtained,getthebiggestmodelcredibilityriskForquestiontwo,envelopeysismethodtoestablishthemodeloffuzzyclusteringysisbasedontheriskassessment,giveninsomeindexestousethemethodofthecolumnsinthepermissiblevariation.Firstofall,thecapitalstructure,capitalreturnrateandthedynamicrecoveryperiodofthreeriskindexasthebasisforclassification,usingfuzzyclusteringysisofoptimalclassificationmethod,theclassificationof72projectsandobtainedtheclusteringcenter;then,ysisoftheeffectivenessofclassification;finally,theenvelopeysisriskassessmentprojectcredibility,sothatintherangeofcredibilityriskInthispaper,combinedwithdataenvelopmentysismethodandfuzzymathematicstheory,areliabletheoreticalbasis;comprehensiveconsiderationofannexestwoandthreedata,comprehensiveysis;atthesametime,thetativemathematicalsoftwareanddataysis,theresultistrueandreliable;theresultsmostlyfromthecharts,clearand:Landreserverisk;Dataenvelopmentysis;Fuzzyclusteringysis;4§1问题的重述《土地储备管理办法》（资发号）第三条规定：土地储备机构应为市、县批准成立，具有独人资格，隶属于资源管理部门，统一承担门设立的土地储备机构（部门）和出资设立的企业性质的土地储备机构两种。 §2问题的分析3阶段一中的模糊综合评价指标体系模型：Fwifif1f2f3金结构、率、动态回收周期对风险集合隶属度，wi为各指标对应权重。输p输输p输f1,f2,ff1,f2,f输输 §3模型的假设 的的数额比例关系1F2345e率67c结构=额/自有§5模型的建立与求解的有效性越高。因此，针对第一阶段的模糊综合评价指标体系模型，可将专家评分作为投入，将结构、率及动态回收周期的风险隶属度作为产出，利f1,f2,pf1,f2,ppfDf： 0

0 0

0 0

6 0

0 6

0 1个部门打分相对较高，而另外一个部门打分可能较低。这是由于虽然有可能同数据包络分析，简称DEA（DataEnvelopmentysis），是以相对效率概念为为决策单元（DMU））进行相对有效性或效益评价的法。y1iy2i,ysi（xjiyji0）1……n…………y2…… …i…n…………记X[x1,x2,,xn Y[y1,y2,,ynv(v,v,v u(u,u,,u1 IvxvxV xTv和OuyuyuyyT 1 2 1 2 s 与总输入之比的大小来衡量DMUi的有效性。令Eii

xiiEii称为DMUi的效率评价指数。在上式中，权向量u和v都是待定的，他们的每一个DEA的C2R模型（P）：DMU，解以下极大化问题：i yT max

xi

（Pjs.t.yju1(1j1),u0,vj xT maxyT s.t.yjTxT(1jn),xT1,0, 有效（C2R）的；⑵若线性规划（P）的解中存在解i*0,i*0EiiyiTI*1，DMUiDEA有效（C2R）的。变量且有非无穷小）： min((eTseTs s.t.xsx,ys （Dj j

j

j 0,s0,s其中，s(sss是m项输入的松弛变量；ssss是s 的松弛变量；(,,)是n个DMU的组合系数；eT 2Te

；是12很n的正数（一般取10）

i⑴若*1DMUDEA有效（C2R）i⑵若*1s*0,s*0DMUiDEA有效（C2R）回收周期对风险的隶属度（f1,f2,f2）作为产出。通过 得到各决策项目可信度（见表1）及投入权向量（见附表1）和产出的权向量（见附表2）11235E11119E11DMUE11DMUDMUE11E11E11E1E11EEE1EDMU1DMU2DMU31 11DMU40.40046DMU20.42795DMU290056554DMU460.061141DMU {5.0271 6.5433 0.65416可信度达到最大

构（c、率（e）和动态回收周期（pt。ce以五部门专家评分作为输入，三个风险指标作为输出，即结构（c、报ce定义模糊CXX1X2,XnRpRp表示P维实数向量空间。令uik表示第k 第i0uik1uik1,0uikN,1kN,1icvi为第ii k J(U,V)Nc(u) k1idik

xk

k个序列到第i

(,)JUN(uik(l))mNv(l)k

(uik(l)修正U

k2uik(l1) 2d ik)mj1d

i,对给定的0，实际计算时应对取定的初始值进行迭代计算直至

t

x第j类。

）得到各指标的可信度见下表2：表21110.9091510.921090.965770.80717110.866860.988490.956550.751810.918770.9758110.829510.80739 0.95124 0.80523110.778250.9715410.931380.666330.967940.89588 0.74772 0.659250.805640.797710.70870.700970.872690.762130.77227 0.85950.83897 0.742370.657190.955260.813380.913920.969690.839840.89583注：数据顺序是剔除了原始数据第26和69使风险评估的可信度达到最高的结构c、e和动态回收周期pt可以30%12、图3：3210

图 结构允许波动范100321

项目序号图 率允许波动范0 项目序号3首先，基于模糊聚类分析最优分类方法，对72个项目进行分类。运用实结果见图4和表3：12 结1率率0 结构4-1

1 结构4-21 率4-3项E389426和69 对A、B、C、D和E五个项目类别进行聚类的有效性分析，见下 1 Column5聚类有效性——A1 Column6聚类有效性——B 1 Column7聚类有效性——C1 Column8聚类有效性——D10 Column9聚类有效性——E得出各个类别的聚类中心，见表4。4率ABCDE26695各类别的风险评价指标可信度AE10.955860.62152 0.95481 0.63623 0.56593 0.838350.685720.53220.878610.669011110.835740.790350.8605810.72010.86950.725480.711861 0.909090.5796210.838420.469180.498780.79318 0.685530.522970.573540.573540.773440.799140.533443信度内，A类项目的风险评价指标变动幅度必须在40%以内。信度内，B类项目的风险评价指标变动幅度必须在30%以内。信度内，C类项目的风险评价指标变动幅度必须在20%以内。信度内，D类项目的风险评价指标变动幅度必须在30%以内。信度内，E类项目的风险评价指标变动幅度必须在40%以内。6ABCDE32102103动态回收周期动态回收周期10

类图10结构（C）允许变动范 类图11率（e）允许变动范 类12动态回收周期（pt）§6误差分析数，这显然与实际情况是不相符合的。这是数据包络分析法所应注意的。分析法对其进行改善时，所得结果仍然不可能反映各个指标的冗余量或者详细信息，

§7模型的评价DEADEADEA§8模型的改进与推广在聚类分析时，由于各指标的离散程度不同，所以在对72个项目分类时，要对聚类的三个指标赋予不同的权重，使得分类更加客观有效。模糊C均值聚类可以表示为如下的规划 minJ(U,V,c)

2)m dik2k10s.t.k

1ic,1k 1kN,1i1mk

其中wjjj[1].,经济数学与建模[M],：[2].立,数学建模[M],合肥：中国科技大学[3]，数据包络分析——理论和应用[M].：科学普及，2006:5-[4]马占新.任慧龙.戴仰山，DEA方法在多综合评价中的应用研究[J]系统工程与电2002（08）:7-11；[5]. 风险机制[M].:经济管理,2002-10，120-[6].，数据包络分析法(DEA)的研究与应用，中国工程科学,2003年06期[7] /view/30b4379af：2011-5-29,50-51附scores_1=scores([1:25,27:68,70:742669X=scores_1';%用户输入多指标输入矩阵Xfori=1:n;f=[zeros(1,m)-Aeq=[X(:,i)'w(:,i)=LINPROG(f,A,b,Aeq,beq,LB,UB);%DMUiwiE(i,i)=Y(:,i)'*w(m+1m+s,i);%求出DMUi的相对效率值EiiwE%输出相对效率值Eiiomega=w(1:m,:)输出投入权向量ωmu=w(m+1:m+s,:)输出产出权向量scores_1=scores([1:25,27:68,70:742669X=scores_1';%用户输入多指标输入矩阵Xepsilon=10^­10;%定义非无穷小ε=10^(－10)f=[zeros(1,n)-epsilon*ones(1,m+s)1];A=zeros(1,n+m+s+1);b=0;forAeq=[Xeye(m)zeros(m,s)-w(:,i)=LINPROGf,A,b,Aeq,beq,LB,UB);%DMUiwiw%输出最佳权向量globallambdas_minus=w(n+1:n+m,:)%输出s*－s_plus=w(n+m+1:n+m+s,:)%s*＋theta=w(n+m+s+1,:)%输出θX=scores_1';%用户输入多指标输入矩阵XY=[c_1';1./e_1';pt_1Yfori=1:n;f=[zeros(1,m)-Aeq=[X(:,i)'w(:,i)=LINPROG(f,A,b,Aeq,beq,LB,UB);%DMUiwiE(i,i)=Y(