数学建模论文——城市公共系统规划自行车租赁服务.doc

装订线目 录一、问题的重述21.1 问题背景21.2 问题提出2二、问题的分析32.1 问题一的分析32.2 问题二的分析32.3 问题三的分析3三、模型的假设3四、符号的说明4五、模型的建立、求解与结果分析55.1 问题一的模型5模型建立5模型求解7结果分析75.2 问题二的模型13模型建立13模型求解13结果分析145.3 问题三的模型15模型建立15模型求解17结果分析18六、模型的评价196.1 模型的优点196.2 模型的不足196.3 模型的改进19七、参考文献19八、附录20装订线一、问题的重述1.1 问题背景近年来，随着经济的发展，我国各级城市的机动车保有量都进入了持续高速增长时期，但由此所引发的道路拥堵、空气污染也引起了政府以及百姓的极大关注众所周知，建立快速、便捷的城市公共交通体系是解决这一问题的有效手段之一然而，居民居住地和交通站点通常都有一段距离，这段不远的距离以及现实存在的公共交通拥挤现象则使居民乘坐公共交通的意愿降低，公共自行车服务系统已被证明能够从一定程度上缓解这一现象1.2 问题提出公共自行车服务系统是指在*个区域内，隔一定距离规划出一些停放自行车的租赁点（如地铁出口、城市中心等人员密集的地方），一个租赁点放置一定数量的自行车，很多的自行车租赁点共同组成一个网络以形成一个服务系统，居民可以在任意租赁点租、还车辆，费用全免（*些城市收取少量的超时费用，但目的只是用来提高自行车的利用率，不以盈利为目的），根据租赁点自行车的使用频率，避免部分租赁点的自行车短缺或堆积现象发生，将通过调度专用车进行合理调度，以最大程度地满足居民对车辆需求，提高车辆利用率（系统有自动报警功能。

公共自行车租赁服务系统纳入城市公共交通体系，有助于解决公交出行“最后一公里”问题，使公共交通服务网络趋于更加完善西安市经开区公共自行车服务系统于2011年4月开始建设，到目前为止，已建成租赁点30个（1），自行车总量达到850辆目前正在筹备第三期建设根据题目中给出的条件和附录中的信息解决如下几个问题：（1）根据目前经开区网点自行车需求情况等信息，若要求调度平均耗时尽量少，请针对已有的30个租赁点设计最优车辆分配方案、调度方案，并给出完成调度所耗费的时间2）假设经开区公共自行车服务系统三期建设准备投入建设经费200万元，据此建立数学模型，确定新增租赁点数目、位置以及合适的放置车辆数目3）针对问题（2），进一步研究，如果要求在150min内完成调度，是否需要增加调度车辆（购置调度车辆费用由其它项目经费解决，不包含在三期建设提供的200万元经费中间）.并给出该情形下的自行车调度方案装订线二、问题的分析2.1 问题一的分析问题一是一个TSP问题（即货郎担问题），没有提及任何的费用问题，只需要考虑时间以及路程，而在分析过程中路程可以用时间来衡量，装卸自行车也使用时间来衡量，因此这可以看成是一个单目标规划问题，利用lingo编程用穷举的的方法对最短路程进行求解，进而转化成时间，目标是使调度车行驶的时间最少。

2.2 问题二的分析问题二中，经开区公共自行车服务系统三期建设准备投入建设经费200万元，由此建立线性规划数学模型，确定新增租赁点数目、位置以及合适的放置车辆数目我们在此首先考虑居民需求，也就是在满足租赁点自行车数量上限的前提下，使得公共自行车数量达到最大值，确定租赁点相应数目对应地点后，利用几何关系，分析比较新增租赁点分布合理性，从而解出新增租赁点的数目、位置和各个租赁点分配的自行车数量2.3 问题三的分析问题三中，针对问题（2），进一步研究，题目要求在150min内完成调度，若总用时低于150min则不需要增加调度车辆，若总用时大于150min则需要增加调度车辆该问相当于是问题一的拓广，我们在此利用问题一的模型，只需要考虑时间以及路程，而在分析过程中路程可以用时间来衡量，由此转化成关于时间的单目标函数求解，即求解一辆调度车时的最短路径，两辆调度车时的调度总用时，三辆车时的调度总用时……当车辆数目满足调度用时小于150min时，确定此时的调度路线，即为所求解装订线三、模型的假设1. 租赁点之间的距离使用经纬度计算，即假设到达任意租赁点所走的路程就是两租赁点之间的距离，不考虑街道等因素对本题模型的影响。

2. 为简化模型，不考虑调度车启动和停止的时间，假设调度车运输过程匀速行驶，不会受到交通事故、红绿灯等外界因素的影响3. 假设各个租赁点每天的自行车需求量不变，且自行车的需求时间刚好处在二中的时间段4. 假设居民的骑行距离不超过2km，在*个租赁点还车的概率与租车点和还车点的距离成反比，骑行距离超过2km的情况一定不会发生5. 求解过程中，自行车数量出现小数时采用四舍五入取整四、符号的说明第i个租赁点在第s个时段车辆需求量（i=1,2,…,30;s=1,2,3）第i个租赁点在第s个时段车辆实际量（i=1,2,…,30;s=1,2,3）第j个租赁点在第i个租赁点还车的概率（i,j=1,2,…,30）第j个租赁点在第i个租赁点还车的比例系数（i,j=1,2,…,30）第i个租赁点装卸车用时（i=1,2,…,30）第i个租赁点到第j个租赁点调度车用时Z调度总用时Z1调度车行驶用时Z2调度车装卸车辆用时W新增租赁点及车辆总花费K新增租赁点个数编号为l的*租赁点日平均车辆需求量（l=1,2,…,70）m新增车辆的数目编号为l的*租赁点实际车辆数装订线五、模型的建立、求解与结果分析5.1 问题一的模型模型建立在问题一中，我们考虑编号为1到30的租点。

建立以调度车行驶时间最少为目标的模型，我们结合调度车行驶的方向性，运用lingo编程找出任意两个租赁点的可行最短路径，也就是运用穷举的方法筛选出最短的行驶路径，从而求出从一个点到另一个点花费的时间设租赁点之间的距离用矩阵来表示，表示租赁点与租赁点之间的距离而租赁点之间的距离本模型按照题目要求用经纬度进行计算，此运算过程用Matlab编程实现，设0-1矩阵用来表示经过的个城市之间的路线设：考虑每个租赁点后只有一个租赁点，则考虑每个租赁点前只有一个租赁点，则但仅有以上约束条件不能避免在一次遍历中产生多于一个互不连通回路为此我们引入额外变量，附加一下充分约束条件，即该约束的解释为：①与不会构成回路，若构成回路有=1，=1，则从而有,导致矛盾；②，与不会构成回路，若构成回路有=1，=1，=1，则-1，，从而有，导致矛盾其他情况以此类推于是我们可以得到如下的模型：然后运用lingo程序进行实现，求得在一辆调度车的行驶的最短路径，即本全局最优解本文中有两辆调度车，可以将租赁点合理地分为两组，得出两组的最佳行驶路线，而行驶速率一定，则时间达到最短lingo程序代码如附录二所示租赁点车辆的分配数，需求量（二中网点格式简单的需求量）满足关系第个租赁点下一时间段的分配数表示点在时段的自行车需求量，与相区别。

由于居民可以在任意一个租赁点还车，在*个租赁点还车的概率与租车点和还车点的距离成反比，且假设居民的骑行距离不超过2km；在此我们设在*租赁点还车的概率，依据概率论相关知识必然事件概率和为1，则由Matlab编程可以求得三十个方程的结果，得到反比例系数的三十个数值Matlab程序如附录三所示又由此可以确定，共有900个数值其中包括租赁点的距离大于2km情况，但其概率为0，为避免数据量大导致误差，我们在此并未将概率算出，而是直接将概率使用到了求解自行车分配数的求解程序中，既简化了模型求解步骤，又保证了数据的准确性从其余租赁点还到租赁点的自行车数，=经过Matlab编程可以求得的一系列数值，Matlab程序如附录四所示也就可以得到问题一调度方案中需要调度的最少自行车数第个租赁点装卸车用时=，从而求得目标函数：完成调度所耗费的时间=模型求解经过lingo编程能够得到最短路线经过的租赁点序号：*(1,2)=1, *(2,3)=1, *(3,17)=1, *(4,22)=1, *(5,19)=1, *(6,23)=1, *(7,8)=1, *(8,27)=1, *(9,30)=1, *(10,28)=1, *(11,12)=1, *(12,14)=1, *(13,15)=1, *(14,13)=1, *(15,16)=1, *(16,1)=1, *(17,4)=1, *(18,11)=1, *(19,29)=1, *(20,5)=1, *(21,20)=1, *(22,6)=1, *(23,24)=1, *(24,10)=1, *(25,21)=1, *(26,25)=1, *(27,26)=1, *(28,9)=1, *(29,18)=1, *(30,7)=1而其余全为0。

由以上数据可以得到一辆调度车时的最佳行车回路为：→1→2→3→17→4→22→6→23→24→10→28→9→30→7→8→27→26→25→21→20→5→19→29→18→11→12→14→13→15→16→1→总共用时27.87min本文问题一中有两辆调度车进行运输，考虑到调度问题的复杂性，结合实际情况把租赁点进行分区管理，使调度过程简单、效率高，本模型将只有一辆调度车时最佳行驶路线经过的租赁点划分为两部分，两辆调度车分别经过各自租赁点的时间基本相等分组如下第一组：→11→12→14→13→15→16→1→2→3→17→4→22→6→23→24→第二组：→10→28→9→30→7→8→27→26→25→21→20→5→19→29→18→假设每辆调度车负责的区域一定，则每组的租赁点调度时应该构成行驶路线回路，使用附录二中的程序分别对两组租赁点的最佳行驶路线进行求解，得到最佳行驶路线：第一组：→7→30→28→10→9→8→21→20→5→19→18→25→26→27→7→所花费时间18.75min第二组：→1→15→11→12→13→14→2→3→17→4→22→24→23→6→29→16→1→所花费时间20.02min如下图所示由附录四程序可以得到每个租赁点在三个时间段所还车辆，所需车辆对比如下图：7:00-8:3011:00-12:3017:30-19:00编号所还车辆所需车辆编号所还车辆所需车辆编号所还车辆所需车辆120 15 122 30 120 15 222 23 224 35 222 23 323 38 327 31 323 38 421 38 422 22 421 38 522 17 523 28 522 17 623 32 623 10 623 32 725 13 722 34 725 13 826 40 824 16 826 40 926 26 924 19 926 26 1023 18 1022 37 10。