Matlab：数据包络分析（DEA）企业效率
输入变量：管理成本、营业成本、财务成本、销售成本、营业税金以及资产减值损失
输出变量：净利润
对数据进行min-max标准化（Min-maxnormalization）标准化处理
理解连享会示例代码：计划不通
data = xlsread('dea_data.xlsx',1); % 导入数据：第1列dum文本自动略过
X= data(:,1:2)'; % 投入指标数据：第1列-第2列
Y= data(:,3:4)'; % 产出指标数据：第3列-第4列（第3-第5列？）
n=size(X,2); % 决策单元数（虽然写的2，但识别出来是决策单元数31）
m=size(X,1); % 投入指标数（虽然写的2，但识别出来是投入指标数2）
q=size(Y,1); % 产出指标数（虽然写的2，但识别出来是录入的产出指标数2，实际3）
w = [];
fori = 1:n
f = [zeros(1,n) 1]; % 定义目标函数
Aeq = []; % 没有等式约束
beq = [];
LB = zeros(n+1,1); % 指定下界
UB = [];
A = [X -X(:,i);-Y zeros(q,1)]; % 设定不等式约束（位置2处的索引无效。数组索引必须为正整数或逻辑值）
b = [zeros(m,1);-Y(:,i)];
w(:,i) = linprog(f,A,b,Aeq,beq,LB,UB); % 模型求解end
theta_CCR_input = w(n+1,:)'; % 结果输出理解代码：结论与连享会CCR模型几乎一致
clc,clear
format long
data = xlsread('dea_data.xlsx',1); % 导入数据（借用连享会数据）
X= data(:,[1:2]); % X为输入变量，2为输入变量的个数
X=X';
Y=data(:,[3:5]); %Y为输出变量，5（2+3），3为输出变量的个数
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;
f=[zeros(1,m) -Y(:,i)'];
Aeq=[X(:,i)',zeros(1,s)];beq=1;
w(:,i)=linprog(f,A,b,Aeq,beq,LB,UB);
E(i,i)=Y(:,i)'*w(m+1:m+s,i);
end
theta=diag(E)';
fprintf('用DEA方法对此的相对评价结果为：\n');
disp(theta);
omega=w(1:m,:)
mu=w(m+1:m+s,:)逐年计算企业效率：计划通
clc,clear
format long
data = xlsread('data_dea_year.xlsx',1); % 导入数据
X= data(:,[2:7]); % X为输入变量，6为输入变量的个数
X=X';
Y=data(:,[8:8]); %Y 为输出变量，1为输出变量的个数
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;
f=[zeros(1,m) -Y(:,i)'];
Aeq=[X(:,i)',zeros(1,s)];beq=1;
w(:,i)=linprog(f,A,b,Aeq,beq,LB,UB);
E(i,i)=Y(:,i)'*w(m+1:m+s,i);
end
theta=diag(E)';参考资料
马苏,高良谋,赵光辉.基于Bootstrap-DEA模型的企业生命周期划分及其效率研究[J].中国软科学,2019(11):176-182.
几种常用数据标准化方法 - 简书
Matlab：数据包络分析 (DEA) 入门教程| 连享会主页 数学建模常用模型10 ：数据包络（DEA）分析法（投入产出法）_HaLosec_Wei-CSDN博客_dea模型