資源描述:
《DEA的Matlab程序(數(shù)據(jù)包絡(luò)分析).doc》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、模型((PC2R)的MATLAB程序clearX=[];%用戶輸入多指標(biāo)輸入矩陣XY=[];%用戶輸入多指標(biāo)輸出矩陣Yn=size(X',1);m=size(X,1);s=size(Y,1);A=[-X'Y'];b=zeros(n,1);LB=zeros(m+s,1);UB=[];fori=1:n;f=[zeros(1,m)-Y(:,i)'];Aeq=[X(:,i)'zeros(1,s)];beq=1;w(:,i)=LINPROG(f,A,b,Aeq,beq,LB,UB);%解線性規(guī)劃,得DMU;的最佳權(quán)向量w;E(i,i)=Y(:,
2、i)'*w(m+1:m+s,i);%求出DMUi的相對效率值Eiiendw%輸出最佳權(quán)向量E%輸出相對效率值EiiOmega=w(1:m,:)%輸出投入權(quán)向量。mu=w(m+1:m+s,:)%輸出產(chǎn)出權(quán)向量。模型(DeC2R)的MATLAB程序clearX=[];%用戶輸入多指標(biāo)輸入矩陣XY=[];%用戶輸入多指標(biāo)輸出矩陣Yn=size(X',1);m=size(X,1);s=size(Y,1);epsilon=10^-10;%定義非阿基米德無窮小e=10-10f=[zeros(1,n)-epsilon*ones(1,m+s)1];%
3、目標(biāo)函數(shù)的系數(shù)矩陣:l的系數(shù)為0,s-,s+的系數(shù)為-e×e,q的系數(shù)為1;A=zeros(1,n+m+s+1);b=0;%<=約束;LB=zeros(n+m+s+1,1);UB=[];%變量約束;LB(n+m+s+1)=-Inf;%-Inf表示下限為負無窮大。fori=1:n;Aeq=[Xeye(m)zeros(m,s)-X(:,i)Yzeros(s,m)-eye(s)zeros(s,1)];beq=[zeros(m,1)Y(:,i)];w(:,i)=LINPROG(f,A,b,Aeq,beq,LB,UB);%解線性規(guī)劃,得DMU的
4、最佳權(quán)向量w;endw%輸出最佳權(quán)向量lambda=w(1:n,:)%輸出ls_minus=w(n+1:n+m,:)%輸出s-s_plus=w(n+m+1:n+m+s,:)%輸出s+theta=w(n+m+s+1,:)%輸出q