帮忙纠错！matlab求非阿基米德无穷小的DEA.运行不出结果。

1699

收藏 2013-11-07

clear
x=[60 70 85 106 35;156 200 157 263 105;50 180 100 86 30];
y=[80 60 90 96 30;12 13 20 17 8;27 25 15 28 3;4 2 5 5 1];
n=size(x',1);m=size(x,1);s=size(y,1);
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);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,l)];
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,:)
大家帮看看，哪里错了~ 感谢~~

扫码加我拉你入群

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

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

全部回复

sdrain

2013-11-7 10:32:59

clear
x=[30 60 55 40 70;25 40 70 30 90;130 150 120 70 180];
y=[35 43 76 52 63;60 80 53 42 71];
n=size(x',1);m=size(x,1);s=size(y,1);
A=[-x' y'];
b=zeros(n,1);
LB=zeros(m+s,1);UB=[];
options=optimset('maxiter',1e4,'tolfun',1e-6);
for i=1:n
Aeq=[x(:,i)',zeros(1,s)];beq=1;
f=[zeros(1,m)-y(:,i)'];
w(:,i)=linprog(f,A,b,Aeq,beq,LB,UB);
Eii=y(:,i)'*w(m+1:m+s,i);
for k=1:n;
f=[zeros(1,m),y(:,k)'];
Aeq=[x(:,k)',zeros(1,s);Eii*x(:,i)',-y(:,j)'];
beq=[1;0];
v=linprog(f,A,b,Aeq,beq,LB,UB);
E(i,k)=(y(:,k)'*v(m+1:m+s))/(x(:,k)'*v(1:m));
end
end
E
mean(E)
[Y 1]=sort(mean(E));
fliplr(I)
这是交叉效率，也出错了。请帮忙纠错。谢谢

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群