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

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

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,:)
大家帮看看，哪里错了~ 感谢~~

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：MATLAB matla atlab Mat Lab 阿基米德 matlab

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)
这是交叉效率，也出错了。请帮忙纠错。谢谢