输入变量X

DMU1 1963.55         2208.81         2518.71         3069.07         3205.69         3392.79         3228.64         3418.75         3873.98         4325.83

DMU2 270.3         280.7         283.9         316.7         327.4         367.8         380.4         407.7         457.8         468.0

DMU3 0.479         0.448         0.437         0.469         0.480         0.484         0.500         0.513         0.519         0.524

输出变量Y为

1807.986597        2185.066684        2566.301791        2981.353881        3378.236821        3832.468925        4152.755332        4623.733404        5278.30992        6145.52

程序为

X=[1963.55         2208.81         2518.71         3069.07         3205.69         3392.79         3228.64         3418.75         3873.98         4325.83

;270.3         280.7         283.9         316.7         327.4         367.8         380.4         407.7         457.8         468.0

;0.479         0.448         0.437         0.469         0.480         0.484         0.500         0.513         0.519         0.524

];Y=[1807.986597        2185.066684        2566.301791        2981.353881        3378.236821        3832.468925        4152.755332        4623.733404        5278.30992        6145.52

];

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=[];

for i=1:n

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)];

beq=[zeros(m,1);Y(:,i)];

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,:)';

运行结果为

Optimization terminated.

Optimization terminated.

Optimization terminated.

Optimization terminated.

Optimization terminated.

Optimization terminated.

Optimization terminated.

Optimization terminated.

Optimization terminated.

Optimization terminated.

w =

0    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

0.0000         0    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

0.0000    0.0000         0    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

0.0000    0.0000    0.0000         0    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

0.0000    0.0000    0.0000    0.0000         0    0.0000    0.0000    0.0000    0.0000    0.0000

0.0000    0.0000    0.0000    0.0000    0.0000         0    0.0000    0.0000    0.0000    0.0000

0.0000    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.0000    0.0000    0.0000

0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000         0    0.0000    0.0000

0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000         0    1.1643

0.2942    0.3556    0.4176    0.4851    0.5497    0.6236    0.6757    0.7524    0.8589         0

0.0000    0.0000    0.0000  101.6121  141.0135    0.0000    0.0000    0.0000    0.0000  478.0289

37.5067   29.0613    8.1815    0.0000    0.0000    0.5909   28.1590   36.0184   37.1009    6.6765

0.1563    0.1256    0.0946    0.0820    0.0891    0.0581    0.0986    0.0941    0.0477    0.0000

0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000

0.6481    0.6963    0.7172    0.7169    0.7858    0.7951    0.9054    0.9520    0.9591    1.1532

问题是

DMU1 的效率值在哪里找？求大神解决 很急  谢谢