投入導(dǎo)向超效率DEA程序總是運(yùn)行出錯(cuò)。但是對(duì)程序檢查了很多遍，也不知道究竟是什么錯(cuò)誤。求各位大神指導(dǎo)！謝謝！
Warning: Input arguments must be scalar.
??? Error using ==> vertcat
All rows in the bracketed expression must have the same
number of columns.
程序如下，其中data1為33*11的矩陣，前四列為輸出量，后七列為投入量。
a=data1(:,5:11);
b=data1(:,1:4);
x=a';
y=b';
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;
Aeq=[zeros(m,1) x(:,(i+1):n) eye(m) zeros(m,s) -x(:,i);zeros(s,1) y(:,(i+1):n) zeros(s,m) -eye(s) zeros(s,1);zeros(1,1) ones(1,(i+1):n) zeros(1,m) zeros(1,s) zeros(1,1)];
beq=[zeros(m,1);y(:,i);ones(1)];
w(:,i)=LINPROG(f,A,b,Aeq,beq,LB,UB);
end
for i=2:n-1;
Aeq=[x(:,1:(i-1)) zeros(m,1) x(:,(i+1):n) eye(m) zeros(m,s) -x(:,i);y(:,1:(i-1)) zeros(s,1) y(:,(i+1):n) zeros(s,m) -eye(s) zeros(s,1);ones(1,1:(i-1)) zeros(1,1) ones(1,(i+1):n) zeros(1,m) zeros(1,s) zeros(1,1)];
beq=[zeros(m,1);y(:,i);1];
w(:,i)=LINPROG(f,A,b,Aeq,beq,LB,UB);
end
for i=n;
Aeq=[x(:,1:i-1) zeros(m,1) eye(m) zeros(m,s) -x(:,i);y(:,1:i-1) zeros(s,1) zeros(s,m) -eye(s) zeros(s,1);ones(1,1:i-1) zeros(1,1) zeros(1,m) zeros(1,s) zeros(1,1)];
beq=[zeros(m,1);y(:,i);1];
w(:,i)=LINPROG(f,A,b,Aeq,beq,LB,UB);
end
w
lambda=w(1:n,:)
s_minus=w(n+1:n+m,:)
s_plus=w(n+m+1:n+m+s,:)
theta=w(n+m+s+1,:)