数学建模B题全国一等奖论文.docx

1、数学建模B题全国一等奖论文基于供求匹配率的出租车资源配置模型摘 要 本文针对城市出租车资源配置问题，采用定性与定量相结合的研究方法，建立衡量出租车供求匹配程度的指标，分析打车软件各种补贴方案对所建指标的影响，在充分考虑各方利益的前提下，得到打车软件的最优补贴方案，对城市出租车行业资源优化配置、持续良性发展具有一定的参考意义。 为分析不同时空出租车资源的供求匹配程度，引入出租车资源供求匹配率这一指标，指标的定义为城市中实际运行的出租车辆数与居民出行需要的出租车辆数之比，反映城市中实际运行的出租车辆数与居民出行需要的出租车辆数之间的差异。计算得出2013年出租车供求匹配率为0.7766，表示供不应

2、求。居民出行需要的出租车辆数与居民人均日出行次数、城市总人口数量、居民出行选择乘坐出租车的比例有关，也与每辆出租车日均载客次数、每单载客人数和车辆满载率有关。对于居民人均日出行次数，利用十五个国大中城市的数据，将十二个城市经济指标聚类分析选出每类指标中典型的经济指标，建立居民人均日出行次数与这些典型经济指标间的多元线性回归方程，而与居民出行需要的出租车辆数相关的其他指标可查阅文献或年鉴获得。分析市每天6:00-8:30，11:00-12:30,13:30-14:30,17:00-18:30四个时间段得供求匹配率分别为0.4111，0.5678,0.6062,0.5631，结果显示供不应求。得到

3、、八座城市的出租车资源供求匹配率分别为1.0936、0.8827、0.9430、0.7040、0.7049、0.7666、0.6583、0.5252，表明只有的出租车资源是供大于求，而其余七座城市为供小于求。 为了分析各公司的出租车补贴方案对缓解打车难是否有帮助，定性分析出租车日均载客次数、出租车满载率随打车软件对出租车司机每单补贴金额的变化趋势，分别建立阻滞增长模型，进而分析打车软件对出租车司机每单补贴金额的变化对所建指标的影响。得到的结论为：对于使用打车软件的乘客来说，出租车补贴方案能够缓解打车难的问题；而对于不使用打车软件的乘客来说，出租车补贴方案则不能缓解打车难的问题。 针对打车软件服

4、务平台的最优补贴问题，综合考虑乘客、出租车司机和打车软件公司三方的满意度，利用熵值法确定这三方各自满意度的权重，将三方满意度加权之和作为综合满意度，进而以综合满意度为目标函数，以打车软件对出租车司机每单补贴金额为控制变量，以补贴金额设置的围为约束条件建立优化模型。遍历所有可能的方案得到最优补贴方案为对出租车司机每单补贴9元,综合满意度为0.5710。关键词：聚类分析；回归分析；灰色预测；阻滞增长模型；熵值法；最优化一、问题重述随着经济的发展，近年来，人们对出行的要求不断提高，城市出租车以其方便、快捷、舒适和私密性的特点成为越来越多人的出行选择。但是，国各大城市交通问题日趋严重，“打车难”也是人

5、们关注的一个社会热点问题。数据显示，包括、等众多大城市，出租车非高峰期的空驶率始终在30%上下徘徊，而高峰期却打不到车。这与众多市民反映的打车难背后所隐藏的强烈需求看似形成了一个矛盾。究其原因，最主要的莫过于司机与乘客需求信息不对称，缺乏及时沟通交流的平台。随着“互联网+”时代的到来，有多家公司依托移动互联网建立了打车软件服务平台，实现了乘客与出租车司机之间的信息互通，同时推出了多种出租车的补贴方案，吸引了越来越多的司机和乘客使用打车软件。然而，打车软件同时也导致出租车行业乱象丛生，存在马路扬招成功率降低、乘客怕司机接到大单拒载、司机分心忙于抢单影响行车安全等问题。请你们搜集相关数据，建立数学

6、模型研究如下问题： 问题一：试建立合理的指标，并分析不同时空出租车资源的“供求匹配”程度。 问题二：分析各公司的出租车补贴方案是否对“缓解打车难”有帮助？ 问题三：如果要创建一个新的打车软件服务平台，你们将设计什么样的补贴方案，并论证其合理性。二、问题分析“打车难”是人们关注的一个社会热点问题。但是，在也不是无论何时何地都难打到车。打车难往往出现在特殊时间和地点：上班高峰的住宅区，下班高峰的商务区，凌晨和深夜的郊区或偏僻地点，遇到雨雪天气 “互联网+”时代的出租车资源配置是一个十分复杂的社会问题。要想准确得出合理的资源配置方案难以实现，同时也难以准确收集大量出租车的各项数据如出租车的每天跑单数

7、，收费，拒载情况等，以及不同城市不同城区不同时间居民的出行行为特征数据。为了建立合理的指标，分析不同时空出租车资源的供求匹配程度，首先从城市居民出行对出租车的需求量入手，分析与需求量有关的主要指标，如城市居民出行量。为分析城市居民出行量与城市经济指标的相关性，先将这些指标进行聚类分析，继而得出每类最具代表性的经济指标，再将最具代表性的经济指标与居民出行总量进行回归分析，得到多元线性回归模型，从而预测居民的出行总量。通过查阅文献可以确定居民出行选择出租车作为出行方式的比例从而，计算得出城市的出租车运输量的需求量。然后根据供需平衡法预测出城市出租车需求量。将城市实际出租车数量与城市出租车需求数量作

8、比，得到衡量出租车资源的供求匹配程度的指标即供求匹配率。对未来城市的出租需求量进行灰色关联预测，得到未来城市的出租需求量，通过计算不同城市的出租车需求量，进行不同时空的出租车资源供求匹配的分析。对于各公司的出租车补贴方案是否对“缓解打车难”有帮助问题，由于难以得到各公司不同时间的补贴方案对居民打车难度的实际影响效果数据，我们从公司对每单的补贴金额入手，分析每单补贴金额围为015元，认为补贴金额再高对公司利益有较大损失。司机的每日跑单数为平均每台出租车每天接单数，有效载客率等效为里程利用率。以每单补贴金额为自变量，司机的每日跑单数和有效载客率为因变量，建立阻滞增长模型，基于前面的供求匹配率关系式

9、算出不同每单补贴金额对应的供求匹配率这一指标。为一个新的打车软件服务平台设计合理的补贴方案，我们基于已建立的模型计算出租车单车日均载客次数、有效载客率与打车软件平台对每单补贴金额的关系式。综合考虑每单补贴金额对出租车司机、乘客、打车软件平台三方各自的满意度，利用熵值法求出三方满意度的权重，继而将三方的标准化数据乘以权重得出综合满意度，比较得出最优的补贴方案。三、模型假设1假设城市中的黑车现象对居民出行没有造成影响；2假设所研究的城市没有发生严重的自然灾害和社会动荡；3假设所研究的城市政府对出租车行业的政策基本不变；4假设司机和乘客都是为自身利益考虑，即经济人假设；5假设参考文献中的数据来源可靠

10、，真实可信。四、符号说明符号含义，单位P出租车供求匹配率M市民出行需要的出租车辆数，辆N城市实际运行的出租车辆数，辆Y市民日均出行次数，次/日W城市总人口数量，人市民选择打车出行的比例s出租车单车日均载客次数，次/车日出租车单车日均每次载客人数，人/车次出租车满载率vi城市的第i个经济指标RPearson相关系数w打车软件每单对司机的补贴钱数，元x每辆出租车每天接单的数量，单y 出租车每天的满载率五、模型建立与求解5.1建模前的准备由参考文献2可得到现有指标体系见表1。表1 衡量出租车供求的三大指标指标里程利用率意义指营业里程与行驶里程之比，一般以一辆车为单位，公式为：里程利用率=营业里程（公

11、里）/行驶里程（公里）100%评价这一指标反映车辆载客效率，如果比例高，说明车辆行驶中载客比例高，空驶比较低，对与要车的乘客来说可供租用的车辆不多，乘客等待时间增加，说明供求关系比例紧。如果比例低，则车辆空驶比例高，乘客租用比较方便，但经营者的经济效益下降。国外一些城市把这一指标作为出租汽车发展规划的主要数据。如日本的东京、横滨等城市，把出租汽车里程利用率控制在52%左右，以方便乘客租车，如果里程利用率高于52%，则发展出租汽车，使之降到52%左右。指标车辆满载率意义通过在客流集散较为集中的地点选取几个长期观测点，对单位时间通过道路的载有乘客的出租汽车数量占总通过出租汽车数量的比，公式为：车辆

12、满载率=载客车数（辆）/总通过车辆（辆）100%评价经验认为，车辆满载率应达70%。在实际操作过程中，通过控制出租车的满载率实现运力与运量的适当平衡。在中心城市，当出租车汽车载客率低于70%时，限制出租汽车运力增加；高于70%,增加出租汽车运力，这样对于提高服务质量，满足高峰时运力需求具有重要作用。指标万人拥有量意义是人均设备普指标，用来描述一定城市规模车辆的占有量。公式为：万人拥有量=车辆（辆）/人口规模（万人）评价目前对城市出租车拥有量的控制标准中并没有上限规定，现行的城市道路交通规划设计规（GB50220-95）仅给出了出租车拥有量的下限，即大城市不少于每万人20辆，小城市不少于每千人5

13、辆，中等城市可在其间取值。5.2问题一：建立评价出租车资源供求匹配程度的指标为评价出租车资源的“供求匹配”程度，引入出租车资源的供求匹配率这一指标，指标的定义为城市中实际运行的出租车辆数与市民出行需要的出租车辆数之比，即 (1)其中，引入P表示出租车资源的供求匹配率，M表示市民出行需要的出租车辆数，N表示城市中实际运行的出租车辆数。市民出行需要的出租车辆M的意义是指这些出租车辆能够恰好满足市民打车出行的需求，即城市出租车资源供求平衡时的车辆数。供求匹配率P反映了城市中实际运行的出租车辆数与市民出行需要的出租车辆数之间的差异。供求匹配率P1为出租车资源供求平衡状态，供求匹配率P越接近1，则说明城

14、市出租车资源供求匹配程度越高，出租车数量配置越合理；当供求匹配率P大于1时，表明城市中现有的出租车数量超过市民出行需要的数量，会增加出租车的空驶率，造成出租车司机的收益降低；当供求匹配率P小于1时，表明城市中现有的出租车数量少于市民出行需要的数量，需要增加出租车的数量来缓解打车难的情况。5.2.1建立市民出行需要的出租车辆数M的预测模型 市民出行需要的出租车辆数M与市民人均日出行次数、城市总人口数量、市民出行选择乘坐出租车的比例有关，也与每辆出租车日均载客次数、每单载客人数和车辆满载率有关，具体关系式为：符号表达式为： (2)其中，M表示市民出行需要的出租车辆数（辆），Y表示市民人均日出行次数

15、（单位：次/人日），W为城市总人口数量（人），表示市民选择打车出行的比例，s表示出租车单车日均载客次数（单位：次/车日），表示出租车单车日均每次载客人数（人/车次），为出租车满载率。根据参考文献3，选取出租车单车日均载客次数s=35（次/车日），出租车单车日均每次载客人数=2.0（人/车次），出租车满载率=65%，居民选择打车出行的比例为6%。 接下来建立市民人均日出行次数与城市经济指标关联的量化模型 市民人均日出行次数是居民出行强度的最直接反映，其与城市人口总数量的乘积即为市民的出行总量，而市民人均日出行次数与城市经济指标有着极大联系。通常情况下，市民人均日出行次数的多少与出行目的、城市布局

16、、交通设施、城市环境质量等因素有关。对于某一城市来说，影响居民人均出行次数的因素又间接地反映在该城市的相关经济指标上。因此，多种因素与市民人均日出行次数的在关联可以转化为多种经济指标与市民人均日出行次数的在关联。STEP1：各经济指标的聚类分析；STEP2：典型指标的选取；STEP3：回归模型的建立STEP4：模型的检验。聚类分析是根据事物本身的特性来定量研究分类问题的一种多元统计分析。其基本思想是按照距离的远近将数据分为若干个类别，以使类别数据的“差异”尽可能小，类别间“差异”尽可能大。所用的变量可以被大致分成两类：对样本个体进行聚类通常称为Q型聚类，对研究变量进行聚类称为R型聚类。选用欧几

17、里得距离（欧式距离）来度量指标之间接近的程度。欧式距离就是空间中两点之间的直线距离，其中各特征参数是等权的，记dij表示指标vi和vj之间的距离，则有计算公式如下： 聚类分析具体过程如下：（1）首先将各聚类单位各自作为一类（这时有p类），按照所选取的距离计算各数据点之间的距离，形成一个距离阵。（2）将距离最近的两个单位并为一个类别，形成n-1个类别，计算新产生的类别与其他各类别之间的距离，形成新的距离阵。（3）按照和第二步相同的原则，再将距离最接近的两个类别合并，这时如果类别个数仍然大于1，则继续重复这一步骤，直到所有的数据都被合并为成为一个类别为止。STEP1：选取、十五个大中城市为研究对象

18、，分析各城市人均日出行次数和十二个经济指标之间的关联。十五个大中城市2001年人均日出行次数和各经济指标见表2。表2 大中城市居民人均出行次数和经济指标城市居民人均出行次数/（次/人天）市区土地面积/(km2)市区总人口/(万人)第三产业值/亿元人均GDP/元工业产品销售收入/（亿元）社会消费品零售总额/（亿元）经济指标符号v1v2v3v4v5v62.431650209.5669303842643.7391.542.442599371.9544.5205971672.5465.832.8628089.539.119704944.7222.611.9553001262.42509.8373827

19、2131861.32.641041676.81660.9253003006.91593.42.447417748856.6199862983.3832.72.073068379.5613250741828.3458.822.721043153.843618034699.7386.281.8637195771452.6380002811.31248.281.5939283.2879.6433442971.6609.263.04163375.9148.648931634.2128.441.882415270.7550223481077.7534.22.543583298423.514274920.

20、9358.32.848494758.2667.917882980.4685.822.561418341.5682.214665726.4627.52城市居民人均出行次数（次/人.天）房地产开发资/（亿元）城乡居民储蓄存款/（亿元）居民消费价格总指数/（%）在岗+职工平均工资/元市区居民人均可支配收入/元农民纯收入/元经济指标符号v7v8v9v10v11v122.4368.6936.999.51367010512.0057902.44111.0716.199.9165758848.0043112.8632.2412.7100.1131089406.0047191.95630.73001.9100

21、.02178112883.0058502.64783.83536.3103.11915511578.0050992.44161.31285.0101.2143088959.0048252.07140.9941.899.51831910896.0048962.7289.5550.198.7127609053.0040201.86387.02600.498.92277214694.0064461.59302.61373.497.82594123544.0098693.0434.1248.998.61704015870.0048001.88115.5839.099.5134937418.003900

22、2.5448.5608.1102.3110906339.0028752.84115.3802.099.5113147305.0031002.56170.8995.5100.8124938128.003111 注：数据来源：参考文献1。以12个经济指标为聚类单位，指标与指标间的距离选用欧式距离，采用组间平均联接法，进行聚类分析。利用SPSS19.0软件进行聚类分析，得到各经济指标间的相关系数矩阵如表3，聚类过程中的运算结果参数见表4所示，聚类分析的谱系图如图1。表3 相关系数矩阵Rv1v2v3v4v5v6v7v8v9v10v11v12v11.00.711.281-.186.262.307.027

23、.112.137-.173-.326-.281v2.7111.00.830-.014.743.850.651.693.345.208-.144-.064v3.281.8301.00.305.900.972.899.924.207.591.254.342v4-.186-.014.3051.00.420.198.291.279-.537.766.861.734v5.262.743.900.4201.00.815.729.753.062.624.356.502v6.307.850.972.198.8151.00.935.959.329.519.164.247v7.027.651.899.291.72

24、9.9351.00.973.343.625.329.370v8.112.693.924.279.753.959.9731.00.346.611.281.365v9.137.345.207-.537.062.329.343.3461.00-.332-.527-.469v10-.173.208.591.766.624.519.625.611-.3321.00.866.864v11-.326-.144.254.861.356.164.329.281-.527.8661.00.918v12-.281-.064.342.734.502.247.370.365-.469.864.9181.00表4 聚类分

25、析参数聚类表阶群集组合系数首次出现阶群集下一阶群集 1群集 2群集 1群集 2136831.753004279945.2510033271490.7300244232535.5703155284117.3344066257813.49050771212482.022068811216855.22370109101122008.99200101011049884.8798911111495789.8111000图1 聚类分析的谱系图根据聚类分析的谱系图可以看出，这十二个经济指标可分为三大类：第一大类包括人均GDP（v4）；第二大类包括在岗职工平均工资（v10）和市区居民人均可支配收入（v11）；

26、第三大类包括市区面积（v1）、市区总人口数（v2）、第三产业值（v3）、工业产品销售收入（v5）、社会消费品零售总额（v6）、房地产开发投资（v7）、城乡居民储蓄存款（v8）、居民消费价格总指数（v9）、农民人均纯收入（v12）。聚类分析的结果见表5表5 聚类分析结果指标代码v1v2v3v4v5v6v7v8v9v10v11v12分类类别333133333223STEP2：将十二个指标进行聚类分析得到三类指标后，出来第一类指标外，其他两类都包含多个指标，为了选取各类指标中的典型指标，分别计算类中每一变量与其余变量的Pearson相关系数，继而求得类中每一变量与其余变量的Pearson相关系数的平

27、方和的平均值，最后把该值最大的变量作为典型指标。 Pearson相关系数的定义为 （3）相关系数R的数值围是介于-1与+1之间: 如果|R|0，表明两个变量没有线性相关关系。 如果|R|1，则表示两个变量完全直线相关。线性相关的方向通过相关系数的符号来表示，“+”号表示正相关，“-”表示负相关。 （4）式中， 表示变量 与其余变量 的Pearson相关系数R的平方和的平均值，m为变量的个数。 将表3中的值代入公式(3)、(4)，计算得到第二类和第三类指标中的各变量的，具体见表6所示。表6 同一类中每一变量与其余变量的相关系数平方和的均值指标第二类指标第三类指标v10v11v1v2v3v5v6v7v8v9v120.4580.3760.1850.4150.5320.4670.5200.4930.5060.2210.376 从表6中可以看出，第二类指标中v10对应的值最大，第三类指标v3中对应的值最大。所以，选择v10、v3为典型数据。而第一类指标只有一个指标，故可选取v4、v10、 v3分别作为第一、二、三类指标的典型指标。STEP3：根据聚类分析得到的典型变量，利用表2中的数据，建立2001年人均日出行次数Y与第三产业值（v3）、人均GDP（v4）、在岗+职工平均工资（v10）的多元线性回归方程。