VRS、CRS条件下非期望产出超效率SBM模型，以及普通SBM模型（可计算冗余度）

之前更新了一个帖子是VRS条件下非径向非导向条件下超效率SBM模型，现如今更新出CRS条件下（两者可配合使用），一般CRS条件下测出来的结果更好分析（没有效率值忽高忽低现象）；此外尽管现在众多期刊上的文章不在分析冗余度，鉴于大家需求，本期也更新了另外两个代码，分别是VRS、CRS条件下非导向SBM模型（可计算松弛变量，可进一步分析冗余度）

---------------------------------------------------------------------------------------------------------------------

部分代码如下：（该代码求解VRS条件下带非期望产出的SBM模型，可计算松弛变量）

X =[7.5455 13 16 18 9 7 6 14
    28.3636 70 80 63 34 26 86 79];%投入指标数据

Y = [12.6818 16 7 9 15 12 6 5
    0.5 0.7 0.65 0.5 0.62 0.42 0.3 0.34];%期望产出指标数据

Z = [4.5827 8 9 6 5.3 4.56 7 10 
    12 13 10 15 16 8 9 7];%非期望产出指标数据

[m,n]=size(X);

s=size(Y,1);

q=size(Z,1);

D=1./(m*X');

E=1./((s+q)*Y');

F=1./((s+q)*Z');

A=[];b=[];

LB=zeros(n+m+s+q+1,1);UB=[];

theta=zeros(n,1);

w=zeros(n+m+s+q+1,n);

for i=1:n

f=[zeros(1,n) -D(i,:) zeros(1,s+q) 1];