我也在找，加你好友，找到的话互相说下啊，谢谢

matlab code

帮同学写的.... X, Y   列数要相同

unction [w, lambda, s_minus, s_plus, theta]=eda(X,Y)
[m,n]=size(X);
[s,l]=size(Y);
if n~=l
    disp('the Matrix X, Y must have the same num of col');
else
    epsilon=10^-6;    %定义非阿基米德无穷小
    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:n
        Aeq=[X  eye(m)   zeros(m,s)   -X(:,i);Y  zeros(s,m)   -eye(s)  zeros(s,1) ];
        beq=[zeros(m,1); Y(:,i)];
        w(:,i)=LINPROG(f,[],[],Aeq,beq,LB,UB); %解线性规划的DMU的最佳权向量
    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,:);
end

转载来的~非原创~有懂的人可以帮忙看看么~

http://blog.csdn.net/aris_zzy/article/details/1336087

同问，不知道解决了没？？？

这个用dea-xp1可以求θ和λ吗？

有非阿基米德无穷小量的话，求解的时候不是要用两阶段算法么？应该在编程的时候也要分两步来吧。。。我是这么想的，但其实直接给无穷小量赋一个比较小的值也能求出效率值，就是不知道结果对不对啦。还是慎重一些，多看些文献吧

MATLAB可以

MAXDEA可以解决