数据包络分析法（matlab）.zip

%本程序用于计算DEA:C^2R-improved,计算时去掉本身 %修改程序时间2010.8.15 %清空内存变量与工作空间，以防干扰 clear; clc; %从Excel文件中读取决策单元的投入产出数据 %将原始数据按照规定格式输入到data_dea.xls中 %将数据文件data_dea.xls存放到MATLAB的当前工作目录下 disp('请在弹出的Excel文件data_dea.xls中，选择各决策单元的投入数据') %tmp_input存放各决策单元的投入数据 tmp_input=xlsread('data_dea.xls',-1); disp('请在弹出的Excel文件data_dea.xls中，选择各决策单元的产出数据') %tmp_output存放各决策单元的产出数据 tmp_output=xlsread('data_dea.xls',-1); %从投入产出数据中获取决策单元个数、投入变量个数及产出变量个数信息 %m表示投入的变量个数 %n表示决策单元个数 [m,n]=size(tmp_input); %s表示产出的变量个数 [s,n]=size(tmp_output); %确定等式约束系数矩阵及常数项向量的维数 %mA表示系数矩阵的行数 mA=m+s; %nA表示系数矩阵的列数 nA=n+mA+1; %mb表示常数项向量的行数 mb=mA; %给系数矩阵A预分配存储空间 A=zeros(mA,nA); %构造系数矩阵 %将投入的数据放到系数矩阵的左上区域 A(1:m,1:n)=tmp_input; %将产出的数据放到系数矩阵的左下区域 A(m+1:mA,1:n)=tmp_output; %投入的松弛变量的系数矩阵，为单位对角阵 A(1:m,n+1:n+m)=eye(m); %产出的剩余变量的系数矩阵，为负单位对角阵 A(m+1:mA,n+1+m:nA-1)=-eye(s); %构造目标函数，用向量形式表示 %目标函数系数，依次分别是lamda的系数（=0）,松弛变量和剩余变量及theta的系数 %非阿基米德无穷小量取(1e-6) F=[zeros(n,1);-(1e-6)*ones(mA,1);1]; %变量的下限约束lb，保证变量大于等于0，无上限约束 lb=zeros(nA-1,1); %预分配常数项向量的存储空间 b=zeros(mA,1); %预分配输出变量的存储空间 outputx=zeros(nA-1,n); %预分配目标函数值的存储空间 outputfval=zeros(1,n); %调用linprog函数，计算各决策单元的投入产出效率 for i=1:n %theta前的系数，移到等式左边，为负 A(1:m,nA)=-A(1:m,i); %b向量，取自待评价单元的对应数据 b(m+1:mA,1)=A(m+1:mA,i); %为了防止A在若干次循环后维数（列）不断缩减，引入临时变量AA，保证A的完整性 AA=A; %改进的算法中，不考虑本身那一列 AA(:,i)=[]; %为了防止F在若干次循环后维数（列）不断缩减，引入临时变量FF，保证F的完整性 FF=F; %改进的算法中，不考虑本身那一列，待评价列的系数也要同时删除 FF(i)=[]; [x,fval]=linprog(FF,[],[],AA,b,lb); %将结果保存，包含每个单元前的系数，松弛变量，剩余变量，theta outputx(1:nA-1,i)=x; %目标函数的值 outputfval(1,i)=fval; end %运行结果的格式化输出 %输出文件为output.xls %output.xls的存放路径为MATLAB的当前工作目录 %提示输出lamda值 xlswrite('output.xls','lamda值',1,'A1'); %列出lamda的值 xlswrite('output.xls',outputx(1:n-1,:),1,'A2'); %数据输出时单元格位置定位 k1=n+2; string_1=['A',num2str(k1)]; %提示输出松弛变量值 xlswrite('output.xls','松弛变量值',1,string_1); %列出松弛变量 string_2=['A',num2str(n+3)]; xlswrite('output.xls',outputx(n:n+m-1,:),1,string_2); %提示输出剩余变量 string_3=['A',num2str(n+m+4)]; xlswrite('output.xls','剩余变量值',1,string_3); %列出剩余变量值 string_4=['A',num2str(n+m+6)]; xlswrite('output.xls',outputx(n+m:nA-2,:),1,string_4); %提示输出theta string_5=['A',num2str(n+m+s+7)]; xlswrite('output.xls','theta值',1,string_5); %输出theta值 string_6=['A',num2str(n+m+s+8)]; xlswrite('output.xls',outputx(nA-1,:),1,string_6); %提示输出目标函数值 string_7=['A',num2str(n+m+s+10)]; xlswrite('output.xls','目标函数值',1,string_7); %输出目标函数的值 string_8=['A',num2str(n+m+s+11)]; xlswrite('output.xls',outputfval,1,string_8);