DEA的Matlab程序(數(shù)據(jù)包絡(luò)分析).doc

DEA的Matlab程序(數(shù)據(jù)包絡(luò)分析).doc

ID:52564232

大小:29.50 KB

頁數(shù):1頁

時間:2020-03-28

DEA的Matlab程序(數(shù)據(jù)包絡(luò)分析).doc_第1頁
資源描述:

《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

當(dāng)前文檔最多預(yù)覽五頁,下載文檔查看全文

此文檔下載收益歸作者所有

當(dāng)前文檔最多預(yù)覽五頁,下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學(xué)公式或PPT動畫的文件,查看預(yù)覽時可能會顯示錯亂或異常,文件下載后無此問題,請放心下載。
2. 本文檔由用戶上傳,版權(quán)歸屬用戶,天天文庫負責(zé)整理代發(fā)布。如果您對本文檔版權(quán)有爭議請及時聯(lián)系客服。
3. 下載前請仔細閱讀文檔內(nèi)容,確認(rèn)文檔內(nèi)容符合您的需求后進行下載,若出現(xiàn)內(nèi)容與標(biāo)題不符可向本站投訴處理。
4. 下載文檔時可能由于網(wǎng)絡(luò)波動等原因無法下載或下載錯誤,付費完成后未能成功下載的用戶請聯(lián)系客服處理。