Matlab DEA 编程 请高手赐教

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

模型见图片文件，以下代码不能运行，请高手赐教：

clear
X=[3696.50 290.90 610.17 63.86 823.02 851.40 354.75 404.42];  
Y=[25374.85 1025.83 1552.10 358.26 5359.02 5507.94 2566.74 1668.75];  
U=[12554.03 1092.84 2232.17 252.84 3307.42 2257.41 844.61 743.12];  
N=size(X,1);
K=size(X',1);
M=size(Y,1);
J=size(U,1);

lb=zeros(K+N+M+1+1,1);
lb(K+N+M+1:K+N+M+1+1)=-inf;
ub=[];
ub(K+N+M+1)=1;
A=zeros(1,K+N+M+1+1);
b=0;

epsilon=10^-10;
f=[zeros(1,K) -epsilon*ones(1,N+M) 0 1];
for i=1:K
Aeq=[X eye(N) zeros(N,M) -X(:,i) zeros(N,1)
     Y zeros(M,N) -eye(M) zeros(M,1) zeros(M,1)
     U zeros(J,N) zeros(J,M) zeros(J,1) -U(:,i)
     ones(1,K) zeros(1,N) zeros(1,M) -1 0];
beq=[zeros(N,1)
     Y(:,i)
     zeros(J,1)
     0];
w(:,i)=linprog(f,A,b,Aeq,beq,lb,ub);   
end

w
z=w(1:K,:)
s_minus=w(K+1:K+N,:)
s_plus=w(K+N+1:K+N+M,:)
beta=w(K+N+M+1,:)
lamda=w(K+N+M+1+1,:)

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：MATLAB matla atlab Lab Atl MATLAB DEA 高手 赐教

楼主，对于方程个数较多的，建议使用lingo，比较好循环

ywh19860616 发表于 2010-11-1 13:01
楼主，对于方程个数较多的，建议使用lingo，比较好循环

谢谢经济学家matlab斑竹的建议，这个程序有错误，斑竹能指出一下？希望有高手出来指正～～～～

数据设置错了 做成一个矩阵

liuxin9023 发表于 2010-11-2 16:52
数据设置错了 做成一个矩阵

版主，请仔细讲解一下，好吗？

我也在此处给困惑了，求解

好东西 谢谢楼主~~~~~~~~