超效率DEA MATLAB代码
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
超效率DEA
clear
X= [];%用户输入多指标输入矩阵X
Y= [];%用户输入多指标输入矩阵Y
n=size(X', 1);
m=size(X,1);
s=size(Y,1);
epsilon=10^-10;%定义非阿基米德无穷小量ε=10-10
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;
for i=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 的最佳权向量wi end
w % 输出最佳权向量
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,:) %输出θclear
X= [];%用户输入多指标输入矩阵X
Y= [];%用户输入多指标输入矩阵Y
n=size(X', 1);
m=size(X,1);
s=size(Y,1);