這是代碼
clear
X=[ ]; %用戶鍵入輸入矩陣X
Y=[ ]; %用戶鍵入輸入矩陣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);
E(i,i)=Y(:,i)'*w(m+1:m+s,i);
for k=1:n
f=[zeros(1,m) Y(:,k)'];
Aeq=[X(:,k)' zeros(1,s);E(i,i)*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)); %得到交叉效率值Eik
end
end
轉載請註明出處，本文鏈接：https://www.uj5u.com/shujuku/17929.html
標籤：疑難問題