matlab如何求解sbm模型,非期望产出的SBM模型代码问题

clc,clear all;

load('OUTPUT.mat')

output = OUTPUT

load('INPUT.mat')

input = INPUT

load('UNOUTPUT.mat')

unoutput = UNOUTPUT

[n ,m]=size(output)%n是单元数。m是期望产出数目

[n,k]=size(input)%n是单元数。k是输入数目

[n,m1]=size(unoutput)%n是单元数。m1是非期望产出数目

Do_ccr=[];

Do_bcc=[];

Di_ccr=[];

Di_bcc=[];

for i=1:n

f1=[-1;zeros(n,1)];

A1=[output(i,:)' -output';zeros(k,1) input'];

b1=[zeros(m,1);input(i,:)'];

Aeq1=[0,ones(1,n)];

beq1=1;

lb1=[1;zeros(n,1)];

[x1,fval]=linprog(f1,A1,b1,[],[],lb1,[]);

[x2,fvalvrs]=linprog(f1,A1,b1,Aeq1,beq1,lb1,[]);

Do_ccr(i,1)=1/x1(1,1);

Do_bcc(i,1)=1/x2(1,1);

end;

for i=1:n

f2=[1;zeros(n