用matlab求解超效率DEA模型运行结果的辨认

6696

收藏 2016-02-20

输入变量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 的效率值在哪里找？求大神解决很急  谢谢