交叉效率DEA程序问题求教

chenxiongmi

5537

收藏 2015-12-13

如下做实证遇到困难，自己把参数代了几次都错了，公式见附件。

	x1	x2	x3	y1	y2
DMU1	30	25	130	35	60
DMU2	60	40	150	43	80
DMU3	55	70	120	76	53
DMU4	40	30	70	52	42
DMU5	70	90	180	63	71

求教各位大虾，这个交叉效率DEA参数怎样带进去。谢谢各位。

clear
X=[];
Y=[];
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=[];
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(:,i)']
      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 I]=sort(mean(E));
fliplr(I)

利用MATLAB进行DEA交叉评价分析_彭育威.pdf
大小:(261.34 KB)

马上下载

扫码加我拉你入群

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

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

全部回复

chenxiongmi

2015-12-15 09:56:25

唉，求助求助，搞到最后我自己弄出来了。
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(5,1);
LB=zeros(5,1);UB=[];
for i=1:5;
Aeq =[X(:,i)',zeros(1,2)];beq=1;
f=[zeros(1,3),-Y(:,i)'];
w(:,i)=linprog(f,A,b,Aeq,beq,LB,UB);
Eii=Y(:,i)'*w(4:5,i);
   for k=1:5;
      f=[zeros(1,3),Y(:,k)'];
      Aeq=[X(:,k)',zeros(1,2);Eii*X(:,i)',-Y(:,i)'];
      beq=[1;0];
      v=linprog(f,A,b,Aeq,beq,LB,UB);
      E(i,k)=(Y(:,k)'*v(4:5))/(X(:,k)'*v(1:3));
   end
end
E
mean(E)
[Y I]=sort(mean(E));
fliplr(I)
运行结果：
>> 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(5,1);
LB=zeros(5,1);UB=[];
for i=1:5;
Aeq =[X(:,i)',zeros(1,2)];beq=1;
f=[zeros(1,3),-Y(:,i)'];
w(:,i)=linprog(f,A,b,Aeq,beq,LB,UB);
Eii=Y(:,i)'*w(4:5,i);
   for k=1:5;
      f=[zeros(1,3),Y(:,k)'];
      Aeq=[X(:,k)',zeros(1,2);Eii*X(:,i)',-Y(:,i)'];
      beq=[1;0];
      v=linprog(f,A,b,Aeq,beq,LB,UB);
      E(i,k)=(Y(:,k)'*v(4:5))/(X(:,k)'*v(1:3));
   end
end
E
mean(E)
[Y I]=sort(mean(E));
fliplr(I)
Optimization terminated.
Optimization terminated.
Exiting: Maximum number of iterations exceeded; increase options.MaxIter.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Exiting: Maximum number of iterations exceeded; increase options.MaxIter.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Exiting: Maximum number of iterations exceeded; increase options.MaxIter.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Exiting: Maximum number of iterations exceeded; increase options.MaxIter.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Optimization terminated.
Optimization terminated.

E =

1.0000 0.6280 0.3155 0.5250 0.3287
0.9181 1.0000 0.5264 0.8700 0.5199
0.6520 0.4970 1.0000 0.9408 0.6183
0.3624 0.3859 0.6082 1.0000 0.4038
1.0000 1.0000 0.8013 1.0000 0.7497

ans =

0.7865 0.7022 0.6503 0.8672 0.5241

ans =

   4    1    2    3    5