一、概述

目前,开发出一种技术,通过明确地考虑多种投入(即资源)的运用和多种产出(即服务)的产生,利用线性规划的方法,用来比较提供相似服务的多个服务单位之间的效率,这项技术被称为数据包络线分析(DEA)。

二、特点及应用

1.适合用于多输出-多输入的有效性综合评价问题,在处理所输出-多输入的有效性评价方面具有绝对优势
2.应用DEA方法建立模型前无需对数据进行量纲化处理
3.无需任何权重假设

比如,可以用于对各省发展状况进行评判、企业管理者评估一所快餐分销店、银行支行、健康诊所或初等学校的生产力等

三、算法详解

DEA是一个线形规划模型,表示为产出对投入的比率。通过对一个特定单位的效率和一组提供相同服务的类似单位的绩效的比较,它试图使服务单位的效率最大化。在这个过程中,获得100%效率的一些单位被称为相对有效率单位,而另外的效率评分低于100%的单位本称为无效率单位。

1) 定义变量

设Ek(k=1,2,……, K)为第k个单位的效率比率,这里K代表评估单位的总数。

设uj(j=1,2,……, M)为第j种产出的系数,这里M代表所考虑的产出种类的总数。变量uj用来衡量产出价值降低一个单位所带来的相对的效率下降。

设vI(I=1,2,……,N)为第I种投入的系数,这里N代表所考虑的投入种类的综合素。变量vI用来衡量投入价值降低一个单位带来的相对的效率下降。

设Ojk为一定时期内由第k个服务单位所创造的第j种产出的观察到的单位的数量。

设Iik为一定时期内由第k个服务单位所使用的第i种投入的实际的单位的数量。

2) 目标函数

目标是找出一组伴随每种产出的系数u和一组伴随每种投入的系数ν,从而给被评估的服务单位最高的可能效率。

式中,e是被评估单位的代码。 这个函数满足这样一个约束条件,当同一组投入和产出的系数(uj和vi)用于所有其他对比服务单位时,没有一个服务单位将超过100%的效率或超过1.0的比率。

3) 约束条件

k=1,2,……,K

式中所有系数值都是正的且非零。

为了用标准线性规划软件求解这个有分数的线性规划,需要进行变形。要注意,目标函数和所有约束条件都是比率而不是线性函数。通过把所评估单位的投入人为地调整为总和1.0,这样等式的目标函数可以重新表述为:

满足以下约束条件:

对于个服务单位,约束条件可类似转化为:

k=1,2,…,K

式中 uj≥0 j=1,2,…,M vi≥0 i=1,2,…,N

关于服务单位的样本数量问题是由在分析种比较所挑选的投入和产出变量的数量所决定的。下列关系式把分析中所使用的服务单位数量K和所考虑的投入种类数N与产出种类数M联系出来,它是基于实证发现和DEA实践的经验:

四、代码实现

例 利用DEA方法对天津市的可持续发展进行评价。在这里选取较具代表性的指标作为输入变量和输出变量,见表1。

表1 各决策单元输入、输出指标值

输入变量:政府财政收入占GDP的比例、环保投资占GDP的比例、每千人科技人员数。输出变量:经济发展(用人均GDP表示)、环境发展(用城市环境质量指数表示;计算过程中,城市环境指数的数值作了归一化处理)。

data.txt如下

14.40 0.65 31.30 3621.00 0.00
16.90 0.72 32.20 3943.00 0.09
15.53 0.72 31.87 4086.67 0.07
15.40 0.76 32.23 4904.67 0.13
14.17 0.76 32.40 6311.67 0.37
13.33 0.69 30.77 8173.33 0.59
12.83 0.61 29.23 10236.00 0.51
13.00 0.63 28.20 12094.33 0.44
13.40 0.75 28.80 13603.33 0.58
14.00 0.84 29.10 14841.00 1.00
代码:

clc,clearformat longload('data.txt');%把原始数据保存在纯文本文件data.txt中X=data(:,[1:3]);%X为输入变量,3为输入变量的个数X=X';Y=data(:,[4:5]);%Y为输出变量,5(3+2),2为输出变量的个数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:nf=[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);endtheta=diag(E)';fprintf('用DEA方法对此的相对评价结果为:\n');disp(theta);omega=w(1:m,:)mu=w(m+1:m+s,:)

参考链接

https://blog.csdn.net/qq_41686130/article/details/81866642

https://wiki.mbalib.com/wiki/DEA

美赛常用算法及代码——(1)DEA数据包络分析法相关推荐

  1. 数据包络分析法(DEA)在数学建模中的应用

    简介 数据包络分析方法(Data Envelopment Analysis,DEA)是运筹学.管理科学与数理经济学交叉研究的一个新领域.它是根据多项投入指标和多项产出指标,利用线性规划的方法,对具有可 ...

  2. 美赛算法及MATLAB代码,美赛常用算法及matlab代码——(3)熵权法

    熵权法真的用的太太太多啦(反正我们队是) 原理很简单,不再讲了(俺觉得就算不懂原理只要会用就OK) 举个栗子(以下数据全为我瞎编) 身高 气质 颜值 小王得分 5.5 3 2 小红得分 5.6 4 3 ...

  3. 2023美赛思路 | 2023美赛C题Matlab代码

    2023美赛思路 | 2023美赛C题Matlab代码 目录 2023美赛思路 | 2023美赛C题Matlab代码 基本介绍 程序设计 运行结果 参考资料 基本介绍 (1)问: 本文分两个小问,第一 ...

  4. 计算机图形学常见算法原理,计算机图形学常用算法及代码大全

    <计算机图形学常用算法及代码大全>由会员分享,可在线阅读,更多相关<计算机图形学常用算法及代码大全(41页珍藏版)>请在人人文库网上搜索. 1.2.1.1 生成直线的DDA算法 ...

  5. matlab包络分析函数,数据包络分析 (MATLAB代码).doc

    数据包络分析 (MATLAB代码) X=[] % 输入多指标输入矩阵X Y=[] % 输入多指标输出矩阵Y [n,m]=size(X) % n为输入的指标数,m为决策单元数 s=size(Y) % s ...

  6. 数据包络分析DEA(数学建模)

    一.模型背景 数据包络分析是线性规划模型的应用之一,常被用来衡量拥有相同目标的运营单位的相对效率. 数据包络分析是一种基于线性规划的用于评价同类型组织(或项目)工作绩效相对有效性的特殊工具手段.这类组 ...

  7. 数据包络分析--两阶段加性网络DEA--使用优化黄金分割比例算法(第一篇)

    数据包络分析-两阶段加性网络-使用黄金分割比例算法 文献介绍 乘子模型 基础信息 黄金分割比例算法 例子 文献介绍 如果阅读过我前面的文章<数据包络分析-加性网络DEA分解>,那么这篇文章 ...

  8. 仁慈型dea matlab程序,数据包络分析(DEA)方法..docx

    数据包络分析(DEA)方法. 二.数据包络分析(DEA)方法数据包络分析(data envelopment analysis, DEA)是由著名运筹学家Charnes, Cooper和Rhodes于1 ...

  9. DEA(数据包络分析)程序模板

    在人们的生产活动和社会活动中常常会遇到这样的问题:经过一段时间之后,需要对具有相同类型的部门或单位(称为决策单元)进行评价,其评价的依据是决策单元的"输入"数据和"输出& ...

  10. 2023美赛F题全部代码+数据+结果 数学建模

    2023年美赛F题全部思路 数据代码都已完成 全部内容见链接:https://www.jdmm.cc/file/2708700/ 1.根据文献选的GGDP的指标,发现GGDP与水资源等有关,由此可以筛 ...

最新文章

  1. 这本书,让我秒懂了微服务架构
  2. android+tv+自动切换,Android TV 重写GridView,实现焦点放大效果
  3. 分布式任务队列 Celery — 应用基础
  4. python 循环语句结果存储_Python条件语句和循环语句
  5. TypeScript constructor signature 类型的变量赋值方式
  6. 【bzoj5107】[CodePlus2017]找爸爸 dp
  7. NNIE-lite 为算法工程师而生
  8. 分布式事务处理--消息发送一致性的异常流程处理
  9. Java基础知识强化84:System类之exit()方法和currentTimeMillis()方法
  10. 2020-10-13 Comsol学习1
  11. JavaScript动画案例——筋斗云,手风琴,固定导航栏
  12. 人工神经网络基本构成有哪些具有什么特征
  13. ORA-1628: max # extents 32765 reached for rollback segment 错误处理
  14. python发微信提醒天气_基于Python实现定时自动给微信好友发送天气预报
  15. Md5工具包加密之MD5Util
  16. 【读书笔记】《数学之美》——一个好方法在形式上总是简单的
  17. i39 1.6 [TNT] 特别版 Mac iPhone 手机铃声制作软件
  18. 【MP4 QuickTime隐写】tcsteg、jpegEOFx、detect anomalies in video files
  19. 2_RP2040软件开发详解
  20. SQLServer 数据库 视图实验练习

热门文章

  1. 樽海鞘群算法原理详解
  2. 地图分幅组件的实现(七)——非标准制图过程的自动化
  3. hashcode值相同的字符串
  4. 国内外各大免费搜索引擎、导航网址提交入口
  5. 疯狂java 视频_疯狂Java讲义配书视频教程 下载
  6. driver接管当前浏览器
  7. 软件开发的流程和常用模式
  8. 【GlobalMapper精品教程】014:矢量线图层的创建及数字化操作
  9. LM2596降压DCDC芯片详解
  10. webpower携手Google、Paypal举跨境电商互动享会