超效率DEA-MATLAB代码

超效率DEAclearX=[];%用户输入多指标输入矩阵XY=[];%用户输入多指标输入矩阵Yn=size(X',1);m=size(X,1);s=size(Y,1);epsilon=10^-10;%定义非阿基米德无穷小量ε=10-10f=[zeros(1,n)-epsilon*ones(1,m+s)1];A=zeros(1,n+m+s+1);b=0LB=zeros(n+m+s+1,1);UB=[];LB(n+m+s+1)=-inf;fori=1:n;Aeq=[[X(:,1:i-1),zeros(m,1),X(:,i+1:n)]eye(m)zeros(m,s)-X(:,i)[Y(:,1:i-1),zeros(s,1),Y(:,i+1:n)]zeros(s,m)-eye(s)zeros(s,1)]beq=[zeros(m,1)Y(:,i)];w(:,i)=linprog(f,A,b,Aeq,beq,LB,UB);%解线性规划,得DMUi的最佳权向量wiendw%输出最佳权向量lambda=w(1:n,:)%输出λs_minus=w(n+1:n+m,:)%输出s-s_plus=w(n+m+1:n+m+s,:)%输出theta=w(n+m+s+1,:)%输出θclearX=[];%用户输入多指标输入矩阵XY=[];%用户输入多指标输入矩阵Yn=size(X',1);m=size(X,1);s=size(Y,1);