经济预测与决策技术及MATLAB实现 第15章 数据包络分析法.pptx

第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=[