1.一种基于方向距离与超效率模型的图纸设计效率评估方法，其特征在于，包括如下步骤：(1)对所有工程数据集C1中的不良产出进行处理，建立工程效率计算数据集合PD；

(2)以设计单位为标准，将其承接过的工程按工程类型绑定到设计单位各工程类型下，建立设计单位工程集合DP；

(3)根据方向距离函数的超效率模型计算工程效率计算数据集合PD，得到每个类型下的工程效率，即计算结果集MC，具体步骤如下：步骤3.1：定义一个类型下的所有工程编号集合PIb＝{pi1,pi2,…,piN}，所有工程类型集合PS＝{PI1,PI2,…,PI13}，其中N≥0，b∈[1,13]，执行步骤3.2；

步骤3.2：定义循环变量n1，用来遍历PS，其中，n1∈[1,13]，n1赋初值为1，执行步骤

3.3；

步骤3.3：定义循环变量n2，用来遍历PIb，其中，n2∈[1,N]，n2赋初值为1，执行步骤3.4；

步骤3.4：定义mcg为集合PS的计算结果，idg为工程编号，typeg为工程类型，typeg＝{1,

2,…,13}，MC为计算结果集，MC＝{{id1,mc1,type1},{id2,mc2,type2},…,{idG,mcG,typeG}}，其中，g∈[1,G]，MC赋初值为空，执行步骤3.5；

步骤3.5：定义tempn1n2为一个工程类型下的一个工程输入产出数据，Temp为输入产出数据集，其中，Temp赋初值为空，执行步骤3.6；

步骤3.6：如果n1≤13，执行步骤3.7，否则执行步骤3.15；

步骤3.7：如果n2≤N，执行步骤3.8，否则执行步骤3.10；

步骤3.8：根据PTn1n2找到IOn1n2，tempn1n2＝IOn1n2，Temp＝Temp∪tempn1n2，执行步骤3.9；

步骤3.9：n2＝n2+1，执行步骤3.7；

步骤3.10：定义循环变量t1，用来遍历Temp，其中，t1∈[1,N]，t1赋初值为1，执行步骤

3.11；

步骤3.11：如果t1≤N，执行步骤3.12，否则执行步骤3.14；

maxξt1

s.t.

步骤3.12： 其中，mcg＝1‑ξt1，将IOn2中的idn2赋值给idg，MC＝MC∪{idg,mcg,n1}，执行步骤3.13；

步骤3.13：t1＝t1+1，执行步骤3.11；

步骤3.14：n1＝n1+1，n2＝1，执行步骤3.6；

步骤3.15：得到工程效率计算结果集MC；

(4)根据设计单位工程集合DP得到设计单位在每个工程类型下的所有工程效率，并评估设计单位在每个工程类型下的综合效率，得到设计单位工程类型效率集合PE；

(5)根据设计单位的所有工程类型效率评估设计单位综合设计效率，得到设计单位综合设计效率集合DE。

2.根据权利要求1所述的一种基于方向距离与超效率模型的图纸设计效率评估方法，其特征在于，所述步骤(1)中建立工程效率计算数据集合PD的具体步骤如下：(1.1)定义工程的输入产出数据集IOg＝{idg,x1,x2,...,x6,...,x17}，工程计算数据集C1＝{IO1,IO2,…,IOG}，其中，idg表示第g个第工程编号，x1,x2,…,x6是输入数据，x7,…,x17是产出数据，x15,x16,x17是不良产出，xi>0，i∈[1,17]，g∈[1,G]，G为C1工程输入与产出的总个数，执行步骤(1.2)；

(1.2)定义循环变量m，用来遍历C1，pdm为处理IOm中不良产出得到的输入产出数据集，PD为所有工程输入产出数据的集合，其中，m∈[1,G]，m赋初值为1，PD赋初值为空，执行步骤(1.3)；

(1.3)如果m≤G，执行步骤(1.4)，否则执行步骤(1.7)；

(1.4)pdm＝IOm，执行步骤(1.5)；

(1.5) 其中pdmi表示处理后的不良产出，c1mi表示处理前的不良产出，PD＝PD∪pdm，执行步骤(1.6)；

(1.6)m＝m+1，执行步骤(1.3)；

(1.7)得到PD＝{pd1,pd2,…,pdG}。

3.根据权利要求1所述的一种基于方向距离与超效率模型的图纸设计效率评估方法，其特征在于，所述步骤(2)中建立设计单位工程集合DP的具体步骤如下：(2.1)定义一个工程类型下的工程编号集合PTh＝{p1,p2,…,pN}，设计单位各工程类型的工程编号集合为DTr1＝{PT1,PT2,…,PT13}，设计单位工程集DP＝{{corpname1,DT1},{corpname2,DT2},…,{corpnameR1,DTR1}}，其中，corpnamer1为第r1个设计单位名称，DTr1为不同工程类型下的工程编号集合，r1∈[1,R1]，R1为r1设计单位工程类型总数，PT赋初值为空，DT赋初值为空，DP赋初值为空，N≥0，h∈[1,13]；

(2.2)定义设计单位名称集合DC＝{dname1,dname2,…,dnameK}，工程类型信息集PR＝{{taskId1,dname1,type1},{taskId2,dname2,type2},…,{taskIdV,dnameV,typeV}}，其中，dnamek为第k个设计单位名称，taskIdv为第v个工程编号，typev为第v个工程类型，typev＝{1,2,…,13}，k∈[1,K]，v∈[1,V]，K为设计单位名称总数，V为工程类型总数；

(2.3)定义循环变量k1，用来遍历DC，其中，k1∈[1,K]，k1赋初值为1；

(2.4)定义循环变量v1，用来遍历PR，其中，v1∈[1,V]，v1赋初值为1；

(2.5)定义循环变量t2，其中，t2∈[1,13]，t2赋初值为1；

(2.6)如果k1≤K，执行步骤(2.7)，否则执行步骤(2.16)；

(2.7)如果t2≤13，执行步骤(2.8)，否则执行步骤(2.14)；

(2.8)如果v1≤V，执行步骤(2.9)，否则执行步骤(2.12)；

(2.9)如果dnamek1＝dnamev1且t2＝typev1，执行步骤(2.10)，否则执行步骤(2.11)；

(2.10)PTt2＝PTt2∪taskIdv1，执行步骤2.11；

(2.12)DTr1＝DTr1∪PTt2，执行步骤2.13；

(2.13)t2＝t2+1，v1＝1，PT赋值为空，执行步骤2.7；

(2.14)DP＝DP∪{dnamek1,DTr1}，执行步骤2.15；

(2.15)k1＝k1+1，t2＝1，DTr1赋值为空，执行步骤2.6；

(2.16)得到DP＝{{corpname1,DT1},{corpname2,DT2},…,{corpnameR1,DTR1}}。

4.根据权利要求1所述的一种基于方向距离与超效率模型的图纸设计效率评估方法，其特征在于，所述步骤(4)中得到设计单位工程类型效率集合PE的具体步骤如下：(4.1)定义设计单位各工程类型的评估效率集为PEi＝{{corpName1,efficiency1,type1},{corpName2,efficiency2,type2},…,{corpNameI,efficiencyI,typeI}}，其中，corpNamei为第i个设计单位名称，efficiencyi为第i个设计单位在typei工程类型下的设计效率，PE赋初值为空；

(4.2)将MC中的效率分为四个区间，分别是(‑∞,0.5),[0.5,1),[1,1.5),[1.5,+∞)，统计各个区间上效率值的个数，选出效率值个数最多的区间，并求此区间上效率值的平均数，定义pe1，其值为求得的平均数；

(4.3)定义变量bool，用来表示设计单位在一个工程类型下是否有工程，值为0时表示没有工程，值为1时表示有工程；

(4.4)定义循环变量d1，遍历DP，d1∈[1,R1]，d1赋初值1，其中DP为设计单位工程集合，R1为设计单位工程类型总数；

(4.5)定义循环变量d2，遍历DTd1，d2∈[1,13]，d2赋初值1；

(4.6)定义pe为设计单位一个工程类型的效率，赋初值0；

(4.7)如果d1≤R1，执行步骤(4.8)，否则执行步骤(4.18)；

(4.8)如果d2≤13，执行步骤(4.9)，否则执行步骤(4.17)；

(4.9)如果DTd1d2表示没有工程，bool＝0，否则bool＝1，执行步骤(4.10)；

(4.10)如果bool＝0，执行步骤(4.11)，否则执行步骤(4.13)；

(4.11)求该设计单位所有工程效率平均值，并赋值给pe2；

(4.12)pe＝pe1+pe2，执行步骤(4.15)；

(4.13)求该设计单位在DTd1d2工程类型下的所有工程效率平均值，并赋值给pe3；

(4.14)pe＝pe1+pe3；

(4.15)将DPd1的corpName赋值给name，PE＝PE∪{name,pe,d2}；

(4.16)d2＝d2+1，执行步骤(4.8)；

(4.17)d1＝d1+1，d2＝1，执行步骤(4.7)；

(4.18)得到设计单位各工程类型的效率评估集PE。

5.根据权利要求1所述的一种基于方向距离与超效率模型的图纸设计效率评估方法，其特征在于，所述步骤(5)中得到设计单位综合设计效率集合DE的具体步骤如下：(5.1)定义设计单位综合效率评估集为DE，DE＝{{corp1,effic1},{corp2,effic2},…,{corpf,efficf}}，其中，corpf为第f个设计单位，efficf为第f个设计单位的综合效率，f>0，DE赋初值为空；

(5.2)定义循环变量z，用来遍历DP，赋初值为1；

(5.3)如果z≤f，执行步骤(5.4)，否则执行步骤(5.7)；

(5.4)根据corpz选出PE中第z个设计单位的所有工程类型效率，并求其平均数，定义cpe，其值为求得的平均数赋值；

(5.5)DE＝DE∪{corpnamez,cpe}；

(5.6)z＝z+1，执行步骤(5.3)；

(5.7)得到设计单位综合效率评估集DE。