城市出租车的规划管理系统-数学建模

城市出租车的规划管理摘要本文通过数学建模的方法解决了城市交通管理中的部分出租车的规划问题。

在问题一的解决上，运用拟合和样条插值的方法预测出2005-2009，2011-2019+的城市市区人口规模。

经检验，样条插值较贴近实际。

在问题二的解决上，运用层次分析法计算出影响出租车数量因素的权重，建立该市出租车数量的动态数学模型。

一、问题重述城市中出租车的需求随着经济发展、城市规模扩大及居民生活方式改变而不断变化。

目前某城市中出租车行业管理存在一定的问题，城市居民普遍反映出租车价格偏高，另一方面，出租车司机却抱怨劳动强度大，收入相对来说偏低，整个出租车行业不景气，长此以往将影响社会稳定。

现为了配合该城市发展的战略目标，最大限度地满足城市中各类人口的出行需要，并协调市民、出租车司机和社会三者的关系，实现该城市交通规划可持续发展，需解决以下的问题：（1）从该城市当前经济发展、城市规模及总体人口规划情况出发，类比国城市情况，预测该城市居民的出行强度和出行总量，这里的居民指的是该城市的常住人口。

同时结合人口出行特征，进一步给出该城市当前与今后若干年乘坐出租车人口的预测模型。

（2）根据该城市的公共出行情况与出租车主要状况，建立出租车最佳数量预测模型。

（3）油价调整（3.87 元/升与4.30 元/升）会影响城市居民与出租车司机的双方的利益关系，给出能够使双方都满意的价格调节最优方案。

（4）针对当前的数据采集情况，提出更合理且实际可行的数据采集方案。

（5）从公用事业管理部门的角度考虑出租车规划的问题，写一篇短文介绍自己的方案。

二、模型假设1．由于第一类人口和第二类人口都对乘出租车产生重大影响，故只考虑人口的总规模。

2．由于城市地理状况和居民的生活习惯在短时期不易改变，所以在各交通小区之间采用的出行方式也相对固定，假定居民从A 地到B 地所习惯采用的出行方式在未来几年保持不变。

3．假设居民中出行人口占总人口数的比例不变。

4．假设对于出行人口而言，在出行方式选择方面的比例与出行人次的比例一样。

5．假设在未来几年，出租车固定营运成本不变。

6．由于每次一起打车的人数，与居民的生活习惯相关，所以假设出租车每趟载客人次不变，即不受出租车数目和收费方案的不同而改变。

7．基于题目给出的图表数据，假定出行与公交数据的统计口径只针对常住人口，不包括其他人口。

8．由于数据的采集统计等存在误差，本文假定所有计算数据在5%~10%误差围可以接受。

三、问题分析题目中要求考虑城市的发展战略目标，人民群众的出行需要，减少环境污染和资源消耗，并结合该城市经济和自身特点，类比国外城市情况，预测该城市居民出行强度和出行总量。

由于题目附录给出的历史数据几乎只有2004 年一年的数据，而做一次出行调查将耗费大量人力物力，所以对一个城市而言也无法得到太多出行特征的历史数据。

为了更好地预测该城市居民的出行强度，必须通过对我国其他城市特别是规模相近城市的居民出行特征的分析，总结出规律并以此来预测。

对于乘坐出租车人口的预测的问题，由于人们生活习惯相对固定，所以在各交通小区之间采用的出行方式也相对固定，又加上常住人口与流动人口都对乘出租车产生重大影响，故只考虑人口的总规模。

对于如何预测该城市出租车的最佳数量，除了考虑乘坐出租车人口数量外，还必须考虑与城区面积，消费水平三者之间的关系，从而确定出租车的最佳拥有量和实际运营数量。

四、符号定义五、模型建立1．城市居民的出行强度和出行总量的预测模型由于城市居民出行强度数据来源于对城市居民出行的大规模调查，缺乏同一城市连续多年的居民出行强度时间序列数据，因此，我们采用选取我国其他城市尤其是规模相近的城市居民出行强度数据，建立多元回归模型，并通过预测出行强度各解释变量的值来预测出行强度。

出行强度预测值与城市市区人口预测值相乘即得到城市居民的出行总量预测值。

出行强度影响因素模型根据我国各城市居民出行调查资料，初步认为城市的人口总量，经济发展水平，建成区面积影响这居民的出行强度，各影响因素与出行强度之间存在着线性相关系，可以建立多元线性回归方程。

选定城市市区人口规模，城市市区GDP ，城市建成区面积为解释变量，居民出行强度为被解释变量，建立多元线性回归方程为：0123D p GDP S ββββε=++++其中D 为出行强度，P 为城市市区人口总量，GDP 为城市市区GDP ，S 为建成区面积，0123,,,ββββ为待定系数，ε为随机变量。

利用SPSS 软件，通过回归方程与回归系数的各项统计检验后，确定出行强度影响因素模型，（1）城市市区人口规模预测：运用拟合和样条插值的方法预测出2005-2009，2011-2019的城市市区人口规模。

经检验，样条插值较贴近实际。

其数值如下表：时间人口（万）2011 329.7212019（2）市区建成面积线性趋势预测模型城市建成区面积指标时间序列在中短期时间存在着比较明显的时间线性变化趋势，对它们未来几年数值的预测属于短期预测，可以采用线性趋势预测模型对未来几年中城市建成区面积指标进行预测。

建立城市建成区面积线性趋势预测模型为：01t S a a t ε=++其中t S 为t 年城市建成区面积，t 为年份，01,a a 为待定系数，ε为随机变量。

利用历史时间系列数据，采用最小二乘法确定系数01,a a 的数值，通过检验后即可以根据灰色预测来预测未来几年市区建成区面积的数值。

（3）出行总量的预测模型根据（1）（2）中对影响出行强度各变量的预测，代入0123D p GDP S ββββε=++++中即可预测出未来T 年的出行强度为：0123t t t t D p GDP S ββββε=++++根据出行总量计算公式，得到未来T 年的出行总量为：t t t Q D P =⨯二．乘坐出租车人口的预测模型、 1.层次分析如果有一组物体，需要知道它们的重量，而又没有均衡器，那么我们就可以通过两两比较它们的相互重要，得出每对物体重量比的判断，从而构成判断矩阵，然后通过求解判断矩阵的最大特征值max λ和它所对应的特征向量，就可以得出这一组物体的相对重量，这样的求解方法称为层次分析法。

层次分析方法的基本过程，大体可以分为如下六个基本步骤 ：建立层次结构。

在这一个步骤中，要求将问题所含的因素进行分组，把每一组作为一个层次，按照最高层（目标层)、若干中间层（准则层)以 及最低层（措施层)的形式排列起来。

这种层次结构常用结构图来表示。

构造判断矩阵。

这一个步骤是层次分析法的一个关键步骤。

判断矩阵表 示针对上一层次中的某元素而言，评定该层次中各有关元素相对重要性的状况。

计算判断矩阵的最大特征值。

为了考察层次分析法得到的结果是否基本 合理，需要对判断矩阵进行一致性检验。

max 1nCI n λ-=-式中， n 代表比较层的因素个数，当 CI 0 时，判断矩阵具有完全一致性；反之，CI 愈大，则判断矩阵的一致性就愈差。

为了进一步检验判断矩阵是否具 有令人满意的一致性，需将 CI 与平均随机一致性指标 RI 比较。

判断矩阵的随机 一致性比例 CR 如下定义：CICR RI=若 CR 0.1 ，则判断矩阵具有令人满意的一致性；若0.1CR ≥，则调整判断矩阵，直到满意为止。

表 8:平均随机一致性指标表阶数 1 2 3 4 5 6 7 8 9 10 0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45 1.49 计算最大特征值对应的特征向量，并对此特征向量进行归一化，得到权向量。

将模型各层的权向量两两相乘得到组合权向量，然后将此权向量归一化，得到方案层对目标层的每一个权值。

若通过组合一致性检验，则将第六步得到的权值作为最终决策的根据。

a.建立乘坐出租车人口的层次结构 目标层：选择怎样的出行方式准则层：公交基本情况 出行习惯 居民消费能力 居民出行强度 方案层：公交 出租车 步行 自行车 摩托车 其它 b.构造成对比较阵在根据以前在市场中做问卷调查，统计后并结合专家给的意见，我们可给出的矩阵 （1）A 为准则层对目标层的正反矩阵1113351313555141111334A ⎛⎫ ⎪ ⎪ ⎪⎪= ⎪ ⎪ ⎪ ⎪⎝⎭运用和法可求出0.17150.21180.54080.0760w ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦4.7198λ= CI =max 1n n λ--= 0.239933333由层次分析给出的随机一致性指标RI 的数值表可知：当n=4 ，RI=0.90 因为0.2399333330.2666<10.90CI CR RI ===，因此所构造的成对比较阵满足一致性检验，下面计算权向量：(2)设(14)i B i =L 分别为方案层对准则层的正反矩阵111543533164331111115624311112144321143113311321153B ⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭21143322217554111111472351121133511311135111511243B ⎛⎫ ⎪ ⎪⎪ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎪ ⎪⎝⎭3111134435111112362111611432143211213121131512312B ⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪= ⎪⎪⎪ ⎪ ⎪⎪ ⎪ ⎪⎝⎭411412531131513411131142311111153114211231136315B ⎛⎫ ⎪ ⎪⎪ ⎪⎪⎪ ⎪=⎪ ⎪⎪⎪⎪⎪ ⎪⎪⎝⎭同理可求得：10.27110.37900.04140.06650.13290.1091w ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦，20.22430.41080.04480.11640.09460.1090w ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦ 30.13550.21420.16780.17210.17210.3104w ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦， 40.25320.18860.10310.08240.15300.2197w ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦ 1 6.0994λ= ，2 6.3319λ= ，39.1745λ= 48.7063λ= 分别对四个特征值作一致性检验可知：1234,,,w w w w 可作为权向量。

c.计算组合权向量i w 为列向量构成的矩阵：0.27110.22430.13550.25320.37900.41080.21420.18860.04140.04480.16780.01310.06650.11640.17210.08240.13290.09460.17210.15300.10910.10900.31040.2197i w ⎛⎫⎪⎪ ⎪=⎪ ⎪ ⎪⎪⎪⎝⎭方案层在目标中的组合权重应为它们相应项的两两乘积之和，即 q=wi*w=[]0.18650.28220.10830.13540.14750.2264对得到的组合权重向量q 归一化,我们就得到了方案层对目标层的组合权值，结果如下表示：表 ：方案层对目标层组合权值2.乘坐出租车人口的模型求解根据以上比例结合城市市区人口规模的预测，可得出2004-----2020年的乘坐出租车人口数量，例2004年的乘坐出租车人口为：240.15*0.2598=62.3910 同理其它年份的乘坐出租车的人口如下表：时间人口（万）20113 62019108.495三．城市出租车最佳数量预测模型（1）首先利用线性拟合预测出附表所缺少月份的数据，并算出02，03，04这三年的年消费支出的量，如下表所示：2002 2003 20041月545.3362 571.7787 452.74132月1000.3 1065.9 1041.23月1505.16 1644.24 1705.714月1986.33 2092.36 2243.845月2331.35 2520.65 2771.56月2738.88 2954.94 3313.47月3194.94 3455.41 3938.618月3657.52 3953.36 4505.789月4255.65 4590.63 5241.3510月4709.9 5103.6 5793.4211月5094.7 5512.9 6337.212月5549.6 6007 6925.7年消费量36569.67 39472.77 44270.45利用线性拟合可以预测出今后若干年的居民消费量，见下表年份 2002 2003 2004 2005 2006 2007 2008 消费量 36250 40100 43950 47810 51660 55510 59360 年份 2009 2010 2011 2012 2013 2014 2015 消费量 63210 67060 70910 74760 78610 8246086310年份 2016 2017 201820192020消费量901609401097860 101710 105560（2）假设城市出租车数量与居民消费水平 城建面积 乘坐出租车人口存在着一定的关系，对此我们对它们构造一个正反矩阵T 目标层：城市出租车的数量准则层：居民消费水平 城建面积 乘坐出租车人口 所构造的正反矩阵T ：111232131313T ⎛⎫⎪⎪= ⎪ ⎪ ⎪⎝⎭求出的0.17200.52420.3038w ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦ 3.2565λ=CI =max 1nn λ--= 0.1282由层次分析给出的随机一致性指标RI 的数值表可知：当n=3 ，RI=0.58 因为0.12820.2210<10.58CI CR RI ===，因此所构造的成对比较阵满足一致性检验， 所以w 可作为权向量。