给你提供一个程序,按提示添加X和Y矩阵就行了
%数据包络分析程序
clear
X=[ … ]; %用户输入多指标输入矩阵X
Y=[ … ]; %用户输入多指标输出矩阵Y
n=size(X',1);m=size(X,1);s=size(Y,1);
epsilon=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 eye(m) zeros(m,s) -X(:,i)
Y 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,:) %输出l *
s_minus=w(n+1:n+m,:) %输出s*-
s_plus=w(n+m+1:n+m+s,:) %输出s*+
theta=w(n+m+s+1,:) %输出q *
|