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

【例15-2】（续【例15-1】）试用带有非阿基米德无穷小量BCC模型12，判断决策单元DMU是否DEA有效。clear%投入矩阵，每一指标数据按行输入X=[35 43.5;2.54.63.8 3.1;23.52.22.5];%产出矩阵，每一指标数据按行输入Y=[2.83.62.53.3;4.16.43.62.6];[m,n]=size(X);%m投入指标个数，n决策单元个数s=size(Y,1);%s产出指标个数epsilon=10^(-7)%定义非阿基米德无穷小量ε=10-7f=[zeros(1,n),-epsilon*ones(1,m+s),1]LB=zeros(n+m+s+1,1);UB=[];A=[];b=[];W=[];forj=1:nAeq=[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];%输出最佳权向量Wendlambda=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,:)%输出θs_minus=0.00000.00000.63640.00000.00000.00000.95450.00000.00000.00000.00000.0000s_plus=0.00000.00000.30000.00000.00000.00000.50000.0000theta=1.00001.00000.90911.0000超效率DEA原始模型：线性规划问题模型：对偶线性规划模型：模型1515.1.4超效率DEA评价模型模型14【例15-3】（续【例15-1】）试用超效率DEA模型判断决策单元DMU是否DEA有效。clear%投入矩阵，每一指标数据按行输入X=[35 43.5;2.54.63.8 3.1;23.52.22.5];%产出矩阵，每一指标数据按行输入Y=[2.83.62.53.3;4.16.43.62.6];[m,n]=size(X);%m投入指标个数，n决策单元个数s=size(Y,1);%s产出指标个数（1）超效率模型14b=zeros(n-1,1);LB=zeros(m+s,1);UB=[];forj=1:nAeq=[X(:,j),zeros(1,s)];beq=1;f=[zeros(1,m),-Y(:,j)];ifj==1A=[-X(:,2:n),Y(:,2:n)];elseifj==nA=[-X(:,1:n-1),Y(:,1:n-1)];elseA=[[-X(:,1:j-1),-X(:,j+1:n)],[Y(:,1:j-1),Y(:,j+1:n)]];endW(:,j)=linprog(f,A,b,Aeq,beq,LB,UB);%解线性规划，得DMUi的最佳权向量e(1,j)=Y(:,j)*W(m+1:m+s,j);endE=e%效应值V=W(1:m,:)%投入权向量U=W(m+1:m+s,:)