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,1)];%(n+1)*1
A2=[-input(i,:)' input';zeros(m,1) -output'];
b2=[zeros(k,1);-output(i,:)'];
Aeq2=[0,ones(1,n)];
beq2=1;
lb2=zeros(1+n,1);
ub2=ones(1+n,1)
[x3,fval]=linprog(f2,A2,b2,[],[],lb2,[]);
[x4,fvalvrs]=linprog(f2,A2,b2,Aeq2,beq2,lb2,[]);
Di_ccr(i,1)=x3(1,1);
Di_bcc(i,1)=x4(1,1);
end;
sbm_eff=[];
sminus=[];
splus=[];
这是我的非期望产出的SBM模型，做出来的结果与实际有很大差距，请问错在哪里呢？
output=[output,unoutput];
for i=1:n
c=ones(k,1)./input(i,:)';%k*1,k是输入个数
d=ones(m+m1,1)./output(i,:)';%(m+m1)*1，m是期望输出个数
f3=[1;(-1/k)*c;zeros(m+m1,1);zeros(n,1)];%(1+k+m+n)*1
Aeq3=[1,zeros(1,k),(1/(m+m1))*d',zeros(1,n);-input(i,:)',eye(k,k),zeros(k,m),zeros(k,m1),input'
-output(i,:)',zeros(m+m1,k),-eye(m+m1,m+m1),output'];
%Aeq3=[1,(1/m)*d',zeros(1,k),zeros(1,n);-input(i,:)',eye(k,k),zeros(k,m),input'
%-output(i,:)',zeros(m,k),-eye(m,m),output'];
beq3=[1;zeros(k,1);zeros(m+m1,1)];
lb3=zeros(1+k+m+m1+n,1);
options=optimset;
[x5,fvalsbm,exitflag]=linprog(f3,[],[],Aeq3,beq3,lb3,[],[],options);
sbm_eff(i,:)=fvalsbm;
sminus(i,1:k)=x5(2:k+1,1)';
sg(i,1:m)=x5(k+2:k+m+1,1)';
sb(i,1:m1)=x5(k+2+m:k+m+1+m1,1)';
end;
save D:\sbm_05.mat Do_ccr Do_bcc Di_ccr Di_bcc sbm_eff sminus splus
