dea的matlab程序

1、Matlab程序代碼程序代碼模型模型((PC2R)的MATLAB程序程序clearX=[]%用戶輸入多指標(biāo)輸入矩陣XY=[]%用戶輸入多指標(biāo)輸出矩陣Yn=size(X1)m=size(X1)s=size(Y1)A=[XY]b=zeros(n1)LB=zeros(ms1)UB=[]fi=1:nf=[zeros(1m)Y(:i)]Aeq=[X(:i)zeros(1s)]beq=1w(:i)=LINPROG(fAbAeqbeqLBUB)%解線

2、性規(guī)劃，得DMU的最佳權(quán)向量wE(ii)=Y(:i)w(m1:msi)%求出DMUi的相對效率值Eiiendw%輸出最佳權(quán)向量E%輸出相對效率值EiiOmega=w(1:m:)%輸出投入權(quán)向量。mu=w(m1:ms:)%輸出產(chǎn)出權(quán)向量。模型模型(D?C2R)的MATLAB程序程序clearX=[]%用戶輸入多指標(biāo)輸入矩陣XY=[]%用戶輸入多指標(biāo)輸出矩陣Yn=size(X1)m=size(X1)s=size(Y1)epsilon=10^

3、10%定義非阿基米德無窮小?=1010f=[zeros(1n)epsilonones(1ms)1]%目標(biāo)函數(shù)的系數(shù)矩陣：?的系數(shù)為0ss的系數(shù)為??e?的系數(shù)為1；A=zeros(1nms1)b=0%=約束；LB=zeros(nms11)UB=[]%變量約束；LB(nms1)=Inf%Inf表示下限為負(fù)無窮大。fi=1:nAeq=[Xeye(m)zeros(ms)X(:i)Yzeros(sm)eye(s)zeros(s1)]beq=[z

4、eros(m1)Y(:i)]w(:i)=LINPROG(fAbAeqbeqLBUB)%解線性規(guī)劃，得DMU的最佳權(quán)向量wendw%輸出最佳權(quán)向量lambda=w(1:n:)%輸出?s_minus=w(n1:nm:)%輸出ss_plus=w(nm1:nms:)%輸出stheta=w(nms1:)%輸出?Matlab程序代碼程序代碼模型模型((PC2R)的MATLAB程序程序clearX=[]%用戶輸入多指標(biāo)輸入矩陣XY=[]%用戶輸入多指

5、標(biāo)輸出矩陣Yn=size(X1)m=size(X1)s=size(Y1)A=[XY]b=zeros(n1)LB=zeros(ms1)UB=[]fi=1:nf=[zeros(1m)Y(:i)]Aeq=[X(:i)zeros(1s)]beq=1w(:i)=LINPROG(fAbAeqbeqLBUB)%解線性規(guī)劃，得DMU的最佳權(quán)向量wE(ii)=Y(:i)w(m1:msi)%求出DMUi的相對效率值Eiiendw%輸出最佳權(quán)向量E%輸出相對

6、效率值EiiOmega=w(1:m:)%輸出投入權(quán)向量。mu=w(m1:ms:)%輸出產(chǎn)出權(quán)向量。模型模型(D?C2R)的MATLAB程序程序clearX=[]%用戶輸入多指標(biāo)輸入矩陣XY=[]%用戶輸入多指標(biāo)輸出矩陣Yn=size(X1)m=size(X1)s=size(Y1)epsilon=10^10%定義非阿基米德無窮小?=1010f=[zeros(1n)epsilonones(1ms)1]%目標(biāo)函數(shù)的系數(shù)矩陣：?的系數(shù)為0ss的

7、系數(shù)為??e?的系數(shù)為1；A=zeros(1nms1)b=0%=約束；LB=zeros(nms11)UB=[]%變量約束；LB(nms1)=Inf%Inf表示下限為負(fù)無窮大。fi=1:nAeq=[Xeye(m)zeros(ms)X(:i)Yzeros(sm)eye(s)zeros(s1)]beq=[zeros(m1)Y(:i)]w(:i)=LINPROG(fAbAeqbeqLBUB)%解線性規(guī)劃，得DMU的最佳權(quán)向量wendw%輸出最佳