保证区DEA模型
（1）CRS
[m,n]=size(X);s=size(Y,1);%X为m*n投入矩阵，m种投入，n个DMU；Y为产出矩阵，s种产出。
f=[1,zeros(1,n+1+2*m-2+2*s-2-1)];
P=zeros(m,2*m-2);%设置投入权重比例上下限，此处为无限制。
L12=0.2;U12=5;%设置产出权重比例上下限
Q=[L12,U12;-1,1;zeros(s-2,2)];
for i=1:169
xo=X(:,i);yo=Y(:,i);
A=[-xo,X,-P,zeros(m,2*s-2);...
  zeros(s,1),-Y,zeros(s,2*m-2),-Q];
b=[zeros(m,1);-yo];
Aeq=[];beq=[];lb=zeros(n+1+2*m-2+2*s-2,1);ub=[];
[x,fval]=linprog(f,A,b,Aeq,beq,lb,ub);
w(i,:)=x;eff(i,1)=fval;
end
（2）VRS
未完待续