第15章数据包络分析法;;;2CCR模型的基本形式;;3DEA有效;;;;;;;15.1.3BCC模型;【例15-2】（续【例15-1】）试用带有非阿基米德无穷小量BCC模型12，判断决策单元DMU是否DEA有效。;epsilon=10^(-7)%定义非阿基米德无穷小量ε=10-7
f=[zeros(1,n),-epsilon*ones(1,m+s),1]
LB=zeros(n+m+s+1,1);UB=[];
A=[];b=[];
W=[];
forj=1:n
Aeq=[X,eye(m),zeros(m,s),-X(:,j);
Y,zeros(s,m),-eye(s),zeros(s,1);
ones(1,n),zeros(1,m+s+1)];
beq=[zeros(m,1);Y(:,j);1];
w=linprog(f,A,b,Aeq,beq,LB,UB);%解线性规划，得DMU的最佳权向量
W=[W,w];%输出最佳权向量W
end
lambda=W(1:n,:)%输出λ
s_minus=W(n+1:n+m,:)%输出s-
s_plus=W(n+m+1:n+m+s,:)%输出s+
theta=W(n+m+s+1,:)%输出θ;;超效率DEA原始模型：;【例15-3】（续【例15-1】）试用超效率DEA模型判断决策单元DMU是否DEA有效。;（1）超效率模型14
b=zeros(n-1,1);
LB=zeros(m+s,1);UB=[];
forj=1:n
Aeq=[X(:,j),zeros(1,s)];
beq=1;
f=[zeros(1,m),-Y(:,j)];
ifj==1
A=[-X(:,2:n),Y(:,2:n)];
elseifj==n
A=[-X(:,1:n-1),Y(:,1:n-1)];
else
A=[[-X(:,1:j-1),-X(:,j+1:n)],[Y(:,1:j-1),Y(:,j+1:n)]];
end
W(:,j)=linprog(f,A,b,Aeq,beq,LB,UB);%解线性规划，得DMUi的最佳权向量
e(1,j)=Y(:,j)*W(m+1:m+s,j);
end
E=e%效应值
V=W(1:m,:)%投入权向量
U=W(m+1:m+s,:)%产出权向量
[M,N]=sort(E,descend)%N从大到小排序
;;%（2）超效率模型15
epsilon=10^(-7)%定义非阿基米德无穷小量ε=10-7
f=[zeros(1,n),-epsilon*ones(1,m+s),1]
A=zeros(1,n+m+s+1)
b=0
LB=zeros(n+m+s+1,1)
UB=[]
LB(n+m+s+1)=-Inf
W=[];
forj=1:n
Aeq=[[X(:,1:j-1),zeros(m,1),X(:,j+1:n)],eye(m),zeros(m,s),-X(:,j);
[Y(:,1:j-1),zeros(s,1),Y(:,j+1:n)],zeros(s,m),-eye(s),zeros(s,1)]
beq=[zeros(m,1);Y(:,j)]
w=linprog(f,A,b,Aeq,beq,LB,UB)
W=[W,w];%最佳权向量
end
lambda=W(1:n,:)%输出λ
s_minus=W(n+1:n+m,:)%输出s-
s_plus=W(n+m+1:n+m+s,:)%输出s+
theta=W(n+m+s+1,:)%输出θ
[M,N]=sort(theta,descend)%N效率从大到小排序;;15.1.5???模效率和技术效率;;15.2案例分析;%每一指标数据按行输入
X=[