X=xlsread('C:\MATLAB\yuxing.xlsx',2,'C2:JP4');

Y=xlsread('C:\MATLAB\yuxing.xlsx',2,'C5:JP7');

Z=xlsread('C:\MATLAB\yuxing.xlsx',2,'C8:JP10');

D=xlsread('C:\MATLAB\yuxing.xlsx',5,'H2:J275');

E=xlsread('C:\MATLAB\yuxing.xlsx',5,'K2:M275');

F=xlsread('C:\MATLAB\yuxing.xlsx',5,'S2:U275');

n=size(X',1);m=size(X,1);s=size(Y,1);q=size(Z,1);

A=zeros(1,n+m+s+q+1);b=0;

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

LB(n+m+s+q+1)=-Inf;

for i=1:n;

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

Aeq=[X eye(m) zeros(m,s+q) -X(:,i)

    Y zeros(s,m) -eye(s) zeros(s,q) -Y(:,i)

    Z zeros(q,m) zeros(q,s) eye(q) -Z(:,i)

    zeros(1,n+m) E(i,:) F(i,:) 1];

beq=[zeros(m,1)

    zeros(s,1)

    zeros(q,1)

     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_plusg=w(n+m+1:n+m+s,:)

s_plusb=w(n+m+s+1:n+m+s+q,:)

theta=w(n+m+s+q+1,:)

本人再写含有非期望产出的动态SBM模型的代码，希望大神指点。