西安市经开区公共自行车服务系统设计.doc

西安市经开区公共自行车服务系统设计摘要本文以西安市经济开发区公共自行车服务系统为背景的车辆分配调度和选址问题。建立快速、便捷的城市公共交通体系是道路拥堵和空气污染问题的有效手段之一，而公共自行车租赁服务系统的纳入使公共交通服务网络趋于更加完善。本文从居民出行需求和交通设施供给角度出发，分析了目前公共自行车的使用特征与问题，建立模型进行求解，对题中三个基本问题进行了全面综合的回答。在现有自行车租赁点信息中，首先根据车辆需求数据建立了车辆分配和调度模型，接着结合西安市的实际数据，采用一种改进的遗传模拟退火算法来求解公共自行车分配和调度问题。为了扩大自行车租赁规模，为广大市民提供便捷的租赁平台，在待选点中确定扩建租赁点数目和位置。本文构建分层评价体系，按人体行为、建设费用、运营协调三个准则量化评价指标，基于TOPSIS选址评价模型，建立指标评价体系进行分析确定网点的具体位置并分配车辆。最后，对第以上问题进一步研究，根据需求平衡确定车辆在限定时间内的调度方案，做到将自行车合理分配。通过实例对模型进行验证结果表明：以上模型能够有效解决城市公共自行车租赁点的布局问题，使公共自行车租赁系统更加有效地运行，达到资源最大化的利用以及最大限度的满足消费者需求的目的. 关键字：公共自行车，交通系统，遗传退火算法，TOPSIS模型，优化目录一、问题重述21.1 问题背景21.2 目标任务3二、问题假设3三、符号说明3四、模型建立与求解44.1 问题一44.1.1车辆分配模型44.1.2.车辆调度模型64.1.3模型算法设计84.1.3.1遗传模拟退火算法的结构流程84.1.3.2 适应度函数84.1.3.3 选择、交叉和变异操作94.1.3.4 模拟退火操作94.1.3.5模型计算104.2问题二124.2.1三层评价体系建立问题的简化124.2.2租赁点方案评价体系建立134.2.3 TOPSIS 模型 选址评价方案154.2.4 模型求解174.3问题三214.3.1车辆调度模型修正214.3.2模型求解22五、模型的评价22参考文献23附录231. 数据图表232.程序代码252.1个体适应度计算252.2比例操作计算262.3交叉变异26一、问题重述1.1 问题背景随着经济的不断发展，我国各级城市的机动车保有量都进入了持续高速增长时期，交通拥堵问题、能源问题、环境问题日益突出，引起了政府以及百姓的极大关注。众所周知，建立快速、便捷的城市公共交通体系是解决这一问题的有效手段之一。然而，居民居住地和交通站点通常都有一段距离，这段不远的距离以及现实存在的公共交通拥挤现象则使居民乘坐公共交通的意愿降低。于是，自行车这种“绿色”交通工具重新得到人们的重视，公共自行车服务系统已被证明能够从一定程度上缓解这一现象。公共自行车租赁服务系统纳入城市公共交通体系，有助于解决公交出行“最后一公里”问题，使公共交通服务网络趋于更加完善。由于其公用性、利用率高、易于管理、中短距离出行成本低、投资成本低的特点,各地政府将其纳入城市公共交通体系并进行大力推广。目前，北京、上海、深圳、济南、郑州、武汉、无锡、佛山、西安等全国30多个大中城市正在逐步建设公共自行车租赁服务系统，它是国内新兴起的一个行业。西安市经开区公共自行车服务系统于2011年4月开始建设，到目前为止，已建成租赁点30个，自行车总量达到850辆。目前正在筹备第三期建设。发展慢行交通,建立公共自行车系统,鼓励更多的出行者采用非机动交通工具,引导居民形成公共自行车+公共交通的出行模式，有助于提高西安城市交通运行效率，有利于减少环境污染。公共自行车系统效益的有效发挥不仅仅与运营模式、租赁点的布局、租赁点车辆配置有关,更与车辆调配密切相关。车辆调配直接影响到公共自行车系统运营效果,因此研究公共自行车实际运营中的车辆调配问题具有很高的研究价值和实际意义。1.2 目标任务 根据西安市经济开发区公共自行车租赁点的设置、需求及位置限制、运营的成本等信息，完成以下问题：问题一：根据目前经开区网点自行车需求情况等信息，若要求调度平均耗时尽量少，请针对已有的30个租赁点设计最优车辆分配方案、调度方案，并给出完成调度所耗费的时间。 问题二：假设经开区公共自行车服务系统三期建设准备投入建设经费200万元，据此建立数学模型，确定新增租赁点数目、位置以及合适的放置车辆数目。问题三：针对问题二，进一步研究，如果要求在150min内完成调度，是否需要增加调度车辆（购置调度车辆费用由其它项目经费解决，不包含在三期建设提供的200万元经费中间）？并给出该情形下的自行车调度方案。0二、问题假设1）调度车可以在任意自行车站点停放，且可以随时出发完成调度任务。2）每一天各个站点需求量基本相同，一天内需求变化规律也不变。3）路网图中描线部分为城市道路，调运车安该路网行驶，其余部分无道路分布。三、符号说明表 1 符号说明符号意义符号意义时间成本（消耗时间）运输车辆数目租赁点数目二进制变量租赁点i的需求量租赁点i到j的最短距离调度车服务完i后服务j时拥有自行车量调运车所能调运的最大车辆数A待选租赁点数效益指标四、模型建立与求解4.1 车辆调度模型首先根据车辆需求数据建立了车辆分配和调度模型，接着结合西安市的实际数据，采用一种改进的遗传模拟退火算法来求解公共自行车分配和调度问题。4.1.1车辆分配模型依据西安市30个租赁点车辆需求数据，采取动态分配车辆模型，首先将7:008:30车辆需求数作为该点初始分配数，共709辆。由于每个站点日变化程度不同，故算得各个站点的需求变化标准差如表2： 表2 各个网点车辆需求变化标准差编号12345678910站点位置经发大厦可口可乐北门经发国际会馆昆仑银行赛高街区西安中学西门运动公园东门运动公园南门管委会出口加工区广场标准差6.135.444.196.604.509.9310.1410.203.307.93编号11121314151617181920站点位置鼎新花园西安外国语学校雅荷花园御道华城天地时代广场市图书馆文景观园长庆电视台移动公司凤城五路标准差3.746.659.673.0911.436.688.066.163.771.25编号21222324252627282930站点位置凤城六路中登家园北门万华园粤华凤城家园市人大市委政务大厅首创国际城中登广场运动公园北门标准差8.7313.449.8811.564.991.252.361.253.688.22表2中反映了各个站点在一天内车辆需求变化的程度，其中中登家园北门变化程度最大，粤华凤城家园次之。一个站点变化程度越大则该站点车辆被调动的可能越大。当需求车辆数从时间段7:008:30到时间段11:0012:30由小变大时称为正变化，相反称为负变化。故将剩余141辆车分配到正变化需求数中标准差变化大的站点中，如表3 给出了待分配站点标准差及变化车辆数。表3 正变化待分配网点变化信息编号站点位置标准差变化量1经发大厦6.1372可口可乐北门5.4415赛高街区4.5067运动公园东门10.142210出口加工区广场7.931311鼎新花园3.74613雅荷花园9.672014御道华城3.09716市图书馆6.68517文景观园8.061818长庆电视台6.161521凤城六路8.731823万华园9.88528首创国际城1.25329中登广场3.68930运动公园北门8.2212依据表3选择标准差大的站点优先分配，分配车辆的数量不能超过表2中的变化量，各站点最终分配车辆不能超过40辆。例如表:3标准差最大为运动公园东门，优先分配22辆，且总量为35辆，没有超过40辆，因此分配成功，其他站点以此类推。最终各个站点分配的车辆如表4： 表4 各网点最终分配车辆编号站点位置初始车辆分配编号站点位置初始车辆分配1经发大厦2216市图书馆222可口可乐北门2417文景观园393经发国际会馆3818长庆电视台384昆仑银行3819移动公司285赛高街区2220凤城五路236西安中学西门3221凤城六路337运动公园东门3522中登家园北门358运动公园南门4023万华园209管委会3924粤华凤城家园3410出口加工区广场1825市人大2111鼎新花园1826市委2312西安外国语学校3527政务大厅3513雅荷花园2728首创国际城1814御道华城1229中登广场1315天地时代广场3830运动公园北门30表4给出了各个站点分配的车辆数，该车辆数是一天中的初始车值，即在早上7:00开始，各个站点达到需求平衡，当有人借车或者存在还车时，该平衡会被打破，此时需要调度车完成各站点之间调度以重新达到平衡。4.1.2.车辆调度模型本次调运系统有2辆调运车，每辆调运车拥有负荷数为，当有租赁点达到上下限时（小于20%或大于90%）,调运车从最近的停车站点出发，负责对各租赁点进行自行车的需求调度服务。完成调度服务后就近回到停车站点，各个租赁点之间的距离以及各自需求量已经确定（需求量见表 4， 各租赁点距离见图 1）。设为所有租赁点的集合，为租赁点数目(n=30)；，m为运输车辆的数目；C为固定时间成本，即每辆自行车装卸平均耗时，为车辆的最大载重数（）；如果车辆被使用，则二进制变量。租赁点，即将服务的车辆的当前拥有车辆数为。对于两个不同的租赁点表示两者之间的最短距离。如果车辆k在服务i后再服务j，则。图1 各租赁点的位置及道路情况图1 中租赁点位置在图中用带圆圈的数字所示，圆圈中数字代表租赁点序号。字代表路线长度（单位：米）。已知运输车速度为，模型的目标函数即运输时间成本，运输时间成本（记为Z）的数学模型如下： （1） （2） （3） （4） （5）式（1）是目标函数，表示最小运输时间成本；式（2）规定了从调度车出发时车上的自行车数量不超过m；式（3）和（4）规定了每个租赁点都服务一次且只服务一次；式（5）规定了每次服务都能完成并且不超过车辆最大载车数。4.1.3模型算法设计在智能优化算法中，遗传算法（GeneticAlgorithm，简称 GA）具有收敛速度快的优点，但是具有局部搜索能力较差并容易早熟收敛的致命弱点。相反，模拟退火算法（Simulated Annealing，简称 SA）能通过概率突跳方式避免陷入局部最小并最终趋于全局最优，但是收敛速度比较慢。基于 GA 和 SA 具有很强的互补性，将 SA 和 GA 有机结合则能增强算法的全局搜索能力和效率。本文采用一种改进的遗传模拟退火算法来求解公共自行车调度问题。4.1.3.1遗传模拟退火算法的结构流程本文提出的遗传模拟退火算法思想是以遗传算法运算流程为主体流程，融入模拟退火机制来调整优化群体。其流程如图 2 所示图 2 遗传模拟退火算法流程4.1.3.2 适应度函数由于要求的是最少的运输成本，是一个最小值问题，因此在设计适应度函数时要把原始目标值转换为适应度值，以确保优秀个体具有大的适应值。通过下式的尺度变换可以将目标值转换为适应度值，即：，其中， I 为当前种群的第i 个染色体， Fitness (I)为适应度函数值，Dmax 为当前种群的最小目标值，Dmin为当前种群的最小目标值， 为需要转换的目标值，本文中取值0.5。使用 可以防止上式被整除，还可以将选择行为从适应度值比例选择调整为纯随机数选择。如果染色体间适应度值的差距较大，则采用适应度值比例选择；如果区间相对较小，则选择趋向于在相互竞争的染色体中进行随机选择。4.1.3.3 选择、交叉和变异操作采用轮盘赌操作对适应度值进行选择。首先生成随机数 （0 1），然后再按照下式进行选择。 （6）其中， 为群体中的第 j 个个体， 为第 j 个个体的适应度值，n为自行车租赁点数量，pop-size为群体大小，通过该操作可以选择出需要繁殖的父代群体。采用单点交叉和均匀变异算子，交叉、变异概率采用自适应的 和，计算表达式如下： （7）式中， 是群体中最大的适应值，是每代群体的平均适应度值，是要交叉的二个个体中较大的适应 度 值 ， 是 要 变 异 个 体 的 适 应 度 值 ， 且。4.1.3.4 模拟退火操作首先选取一个足够大的初始温度 ，因为要使得算法在合理的时间内搜索尽可能大的解空间，只有足够大的 才能满足这个要求；接着设定一个合理的退火率，温度控制参数 的下降函数为=T+1，其中衰减参数是一个略小于 1 的系数；最后终止温度应该设置为足够小。在使用智能优化算法求解问题时，参数的控制十分重要。对以上遗传模拟退火混合算法进行多次测试后，最终选择算法的参数如下表（表 6）。表6 算法参数种群大小（pop-size）60迭代系数（g）100初始温度（）1000降温速度（）0.95交叉概率（）0.30变异概率（）0.40初始接收概率（）0.9994.1.3.5模型计算本模型的目的是寻求目标函数最小，其n个城市之间的距离实质构成了一个的矩阵，同时遗传算法的诸多算子（如选择、交叉、变异），都是针对所谓染色体的，而染色体实质上是一个向量，可以看成的矩阵，因此这些算子实质上是一些矩阵的运算。（1）种群初始化公式（6）已经给出了适应度函数，采用二进制编码法确定变量的编码，同时要对种群大小、最大迭代次数、交叉概率、变异概率等赋值，随机生成种群大小。染色体长度（租赁点数目）为num的MATLAB程序为：Popm=zeros(M,num);For i=i:MPopm(i,:)=randperm(num);end部分适应度计算结果如下：表7 部分路段适应度Fitness 编号两点间距离适应度函数值Fitness (I)编号两点间距离适应度函数值Fitness (I)14-13156.5861160.99687360926-27250.01964990.94853296313-111091.6953850.51306673527-30814.97284710.65623745411-11231.897140.44052914530-9405.40091680.8681417971-15150.543380519-8763.78370310.68272169615-16167.24612250.9913583358-102082.8585590.0002586916-29420.59368890.8602813610-241790.014674018628.4275910.75275224724-231095.3628680.51116925318-19542.68615910.79711315223-22555.36863280.79055149419-5179.72776480.98490058222-6376.72827010.8829764525-20168.47929120.9907203186-171228.9252890.44206672220-21411.21375330.86513435117-3160.4884390.99485462521-251106.2075350.5055584393-2539.71581140.79864995125-26292.51617340.926546112-14615.22098130.759585083（2） 比例选择操作具体执行过程分三步，第一步计算所有个体适应度个体总和；第二步计算每个个体在选择操作中别选中的概率（即使用度大小）；第三步模拟各个个体被选中的次数得到中间群体（程序见附录）(3) 交叉变异操作交叉概率从中间群体中随机的选出需要进行交叉的个体，对这些个体随机的两两配对。在，这两个数表示交叉点的位置。接着对已经配对的两个个体，相互对应的交换（变异操作类似）。（4） 计算结果由以上操作过程可以寻求公式1中的最优解，但两辆车的出发点和结束点不同同样会影响计算结果，具体可以分为三种情况（如下图）。上图红色圆圈表示调度车的起点或者终点，蓝色箭头表示调度车运行的方向，图3 表示两调度车的运行方向一致，手尾相连；图4表示两调度车有共同的起点，最终在同一终点相遇；图5表示两调度车起点不同，终点相同；对表3中各租赁点的需求量按照图3、图4、图5 三种方式调度分配自行车，得到各个方案的最小运行成本和运行轨迹如表8.表8 模型计算结果运输行车线路最小时间成本（min）运行轨迹同方向手尾相连（如图3）62min14-13-11-1-15-16-29-18-19-5-20-21-25-26-27-30-9-8-10-24-23-22-6-17-3-2-14同起点同终点（如图4）71min线路一：21-20-5-19-18-29-16-15-1-11-13-14线路二：21-25-26-27-30-9-8-10-24-23-22-6-17-3-2-14同终点不同起点（图5）79min线路1：21-20-5-19-18-29-16-15-1-11-26线路2：27-30-9-8-10-24-23-22-6-17-3-2-14-25从表8可以看出采用同方向手尾相连运行线路所用时间成本最低，约为52分钟运行轨迹为：14-13-11-1-15-16-29-18-19-5-20-21-25-26-27-30-9-8-10-24-23-22-6-17-3-2-14。其他两种情况调度时间均偏大，主要原因是一条当一个调度车完成调度任务后另一个调度车可能没有完成任务或者有时间间歇，导致资源浪费。而采用同方向手尾相连线路调度能够最大程度的弥补资源浪费的现象，当一个车完成调度时可以协助另一辆车尽早完成调度。4.2选址模型本文构建分层评价体系，按人体行为、建设费用、运营协调三个准则量化评价指标，基于TOPSIS选址评价模型，建立指标评价体系进行分析确定网点的具体位置并分配车辆。4.2.1三层评价体系建立问题的简化确定网点数目及公共自行车的数量，管理者会面对一个问题，就是在增加网点数目和自行车数量的同时投资金额也会同时增加。在本题中，市政能提供的资金为200万元，建设一个网点需要的金额为50000元，投入一辆自行车的成本为1000元。在投资金额一定的前提下，网点数目与自行车数目是此消彼长的，因此在两者之间必须寻求一个平衡，在尽可能满足站点之间距离适当的情况下同时站点能够提供足够使用的自行车。（1）建设的网点数需要不多于筛选出的网点数（2）总资金为200万元，资金约束为由于自行车和网点的数目都为整数，两者又满足资金的约束条件，因此可以通过式求得所有满足条件的网点数以及对应的自行车数总量。这样网点数目和自行车总数都为已知量，为进一步简化计算，还可以对过少的网点数和过多的网点数加以剔除，对这种情况不做考虑。在资金为200万的前提下，最多可设立的网点数为40个，再根据每个网点至多可停放40辆自行车，可得网点数下限为22个。取能够有效减少计算量。根据题意推理，当租赁点停车率小于20%或大于90%，容易发生租赁点的自行车短缺或堆积现象，而每个网点的分配车辆数不能大于40辆，进而可以推出租赁点停车上下限为836辆，折中取安全值为22辆，作为预设网点平均配车数，在此推理下通过上述约束条件（2）可以求得预设网点数目为28个。在已有70个待选网点中在进行28个租赁点的选址决策，我们采用TOPSIS 的选址评价模型，建立指标评价体系。4.2.2租赁点方案评价体系建立1）租赁点选址方案评价体系建立原则方案评价体系的结构和单项评价指面性有着直接的影响。评价标的优劣，对城市公共自行车租赁点选址的科学合理性以及全体系的建立要求能全面、精确的反映目标的本质，并且要具有操作实用性。在建立城市公共自行车租赁点选址方案评价指标体系时，应遵循如下原则：租赁点选址方案评价体系建立原则系统性实用性科学性独立性可比性定性与定量结合图6 评价体系建立原则2）选址方案评价步骤城市公共自行车租赁点选址方案评价体系的建立包含的步骤有：收集资料；分析城市公共自行车租赁点选址决策目标；收集、分析、筛选指标；确定准则层和指标层；指标体系的建立；指标值的确定；选择租赁点选址方案的评价方法；确定指标权重；综合评价和决策 9 个步骤。如图所示：图7 公共自行车租赁点选址指标评价步骤3）选址方案评价体系的建立v 方案评价体系的构建城市公共自行车的选址受多项因素的影响，在决策之前首先要明确站点选址问题的目标。根据前面的指标评价步骤，我们可以将评价体系分为三部分：目标层，准则层和指标层。目标层：充分协调好租赁点规划与实施过程中实际问题之间的矛盾，优化租赁点选址。准则层：选择“以人为本”、“建设费用”、“功能协调”三项指标作为准则层指标。指标层：准确的对该项指标进行分析、量化并进行分解。v 评价指标分析及量化表9 评价体系的构建目标层西安市公共自行车租赁点选址准则层人体行为建设费用运营协调指标层停车步行距离可换乘便捷度服务网点建设自行车购买及养护车辆需求数调度的时空成本停车步行距离公共自行车停车后或者由别的交通方式转乘公共自行车时，出行者步行的距离。步行距离影响公共自行车的利用率，既要满足出行者的换乘要求，又不能过于密集造成资源的浪费。停车步行距离可以利用实际距离来进行量化。换乘便利性租赁点应分散在城市的多处设置以方便租借，可以利用到达某个自行车租赁点所用时间的平均值来进行量化。服务网点建设一个租赁服务网点需要50000元自行车购买及养护在使用周期内，购买、养护一辆自行车需要1000元。车辆需求数每个租赁点能够放置的车辆数目有限，不能超过40辆；为了更好满足居民对车辆的租赁要求、简化调度、提高车辆使用率，通常车辆总数至少应超出需求量的10%；调度的时空成本目前用于运送公共自行车的调度车有2辆，尽量不新增调度车辆，每辆每次可运50辆自行车，调度车平均时速30km/h，每辆自行车装（或卸）平均耗时1min；4.2.3 TOPSIS 模型 选址评价方案TOPSIS法1981年由wangC.L.H和Yoon.K.S首次提出的，这是一种逼近于理想解的排序法。TOPSIS法根据有限个评价对象与理想化目标的接近程度进行排序，评价现有的对象中的相对优劣。“正理想解”和“负理想解”是TOPSIS法的两个基本概念。“正理想解”即一设想的最优解（或方案），它的各个属性值都达到各备选方案中的最好值；而“负理想解”是一设想的最劣解（或方案），它的各个属性值都达到各备选方案中的最坏值。方案排序的规则：比较可行解与“正理想解”和“负理想解”，若其中发现可行解接近正理想解，同时又远离负理想解，那么该可行解为密集的满意解，反之则为最差。 建立 TOPSIS 模型 （1）设有m个待选的公共自行车租赁点A=，影响指标有n个，C=,M=,N=, ,影响指标的权重为,表示待选租赁点对影响指标集的决策矩阵。 （8）Q 为初始决策矩阵。考虑到评价指标的含义和计算方法不同，量纲各异，应先对其进行标准化处理。本文采用线性变换对决策矩阵进行标准化处理得到Q ，该矩阵标准化过程是将指标统一为效益型指标的过程。 对于效益型指标，指标的优越度表示在同类指标中距离最小指标的相对距离，最大指标对最小指标值的优越度为 1；对于成本型指标，指标优越度指在同类指标中距离最大指标的相对距离，最小指标对最大指标的优越度为 1。故可令： （9）其中为效益型指标，为成本型指标。则Q 可表示为 （10）（3）确定影响指标的权重 对于各属性指标权重的确定有多重方法，主要集中在德尔菲（DELPHI）法和层次分析法（AHP）法，通过计算得到各属性的权重，并形成权重向量如下： （11）式中，表示第j中影响指标的权重。（4）形成加权判断矩阵 将归一化的决策矩阵与决策指标的权重系数相结合，构造加权判断矩阵 Z ，其中， （12）（5）确定正负理想方案（14）（13）（6）计算各个备选方案与理想方案之间的距离：待选公共自行车租赁点到正理解（方案）的距离为; 到负理想解（方案）的距离为。（15） , 7）计算各备选到理想方案的贴近度（16）当接近0时，愈接近 0，待选公共自行车租赁点愈靠近负理想方案，该方案可行性越低；当接近1时， 愈接近 0，方案愈靠近理想方案，该备选租赁点的可行性越高。将贴近度进行由小到大排序，贴近度最大的待选租赁点就是最优租赁点的位置。4.2.4 模型求解通过收集资料，对70个待选点进行分析，根据本文提供的6个决策指标，分别为“停车步行距离”，“可换乘便捷度”，“服务网点建设”，“自行车购买及养护，“车辆需求数”，“调度的时空成本”。6个指标的数据如下表所示：表10 待选点决策指标待选点停车步行距离/m可换乘便捷度/min服务网点建设/元自行车购买及平均养护/元平均车辆需求数/辆调度的时空成本/元31600.75000042025 10322002.45000050021 33331201.45000046023 20348015000062022 1335180230 1003203.8500003801253（1）.形成决策矩阵 Q 并标准化得到（2）运用次分析法（AHP）法，通过计算得到各属性的权重，并形成权重向量： (3) 计算加权判断矩阵（4）确定正负理想方案 决策矩阵标准化的过程是将指标统一为效益性指标的过程，因此，以下的计算均以效益型指标计算。 正理想解： 负理想解：（5）计算各个备选方案与理想方案之间的距离自行车选址到正理想方案的距离为到负理想方案的距离为 得到16 23（6）计算各备选到理想方案的贴近度贴近度按大小排列结果为，根据贴近度的大小可得31号待选点是最优的，可以在此置设置公共自行车租赁点；100号待选点贴近度最差，不宜选作租赁点。同理可以求得31100号待选点共70个备选方案的贴近度，依次按大小排列，取前28个点作为最终决策点。最终决策租赁点的位置如下图中五角心点所示。如下图所示，黄色星形点为三期自行车网点的最终的确定位置，三期租赁点最终决策后，经开区公共自行车网点共58个。图8 第三期建设网点决策位置通过TOPSIS 的选址评价模型完成租赁点的最终决策后，各租赁点的相应自行车的配备量可通过问题一所述方法确定，以达到降低调运，最大限度满足居民用车需求的目的。4.3模型三4.3.1车辆调度模型修正为了保证自行车调度系统在150min内完成调度，设调度车辆数为，则采用同方向调度方式时，所消耗的时间由调运时间最长的调运车决定。为了合理配置调配资源，避免资源浪费每辆调运车服务的租赁点个数尽量保持均匀，则对问题一中的调运时间模型修正如下： （17） （18） （19）由于遗传算法适合求解最小值函数，故采用式（18）、（19）的变换。式（17）的约束条件和问题一中式（2）（3）（4）（5）相同，这里不再赘述。4.3.2模型求解采用问题一中遗传算法求解式（17）得出每辆调运车在150min内服务的站点数可以达到20个。因此要完成58个站点的调度工作需要约三个调运车，故增加一辆调运车。五、模型的评价本次将遗传算法应用到车辆调度问题中，由于该算法的整体搜索策略和优化搜索方法在计算方面不依赖梯度信息或其他辅助知识，而只需要影响搜索方向的目标函数和相应的适度函数，所以遗传算法可以应用到复杂系统求最优值中，具体优缺点如下：（1）遗传算法优点：1. 与问题领域无关切快速随机的搜索能力。2. 搜索从群体出发，具有潜在的并行性，可以进行多个个体的同时比较。3. 搜索使用评价函数启发，过程简单。4. 使用概率机制进行迭代，具有随机性。5. 具有可扩展性，容易与其他算法结合。遗传算法的缺点： 1、遗传算法的编程实现比较复杂,首先需要对问题进行编码,找到最优解之后还需要对问题进行解码。 2、另外三个算子的实现也有许多参数,如交叉率和变异率,并且这些参数的选择严重影响解的品质,而目前这些参数的选择大部分是依靠经验。 3、没有能够及时利用网络的反馈信息,故算法的搜索速度比较慢，要得要较精确的解需要较多的训练时间。 4、算法对初始种群的选择有一定的依赖性，能够结合一些启发算法进行改进。（2）Topsis模型优点：topsis模型是一种逼近于理想解的排序法，根据有限个评价对象与理想化目标的接近程度进行排序的方法，是对现有的对象中进行相对优劣的评价，本办法先确定各项指标的正负理想方案，通过计算各方案与理想方案的贴近度来确定最优解，该方法通过对原数据进行同趋势和归一化的处理后，消除了不同指标量纲的影响，并能充分利用原始数据的信息，所以能充分反映各方案之间的差距、客观真实的反映实际情况，具有真实、直观、可靠的优点，而且其对样本资料无特殊要求，故应用日趋广泛。缺点：确定评价指标的权重指标时一般采用专家意见调查法或层次分析法(AHP)等方法，这些方法存在着较大的主观因素，不同人对各个指标的重要程度有不同的评价。参考文献1何流,李旭宏,陈大伟,卢静,吴圆圆. 公共自行车动态调度系统需求预测模型研究J. 武汉理工大学学报(交通科学与工程版),2013,02:278-282. 2李婷婷. 城市公共自行车租赁点选址规划研究D.北京交通大学,2010. 3叶丽霞. 城市公共自行车调度系统研究D.南京理工大学,2013. 4鲍娜. 城市公共自行车租赁点选址决策及调度模型研究D.长安大学,2012 5张建国. 城市公共自行车车辆调配问题研究D.西南交通大学,2013.附录1. 数据图表附表1 经开区已有自行车网点经纬度编号地点经度纬度1经发大厦108.95295434.3248282可口可乐北门108.94856234.3237623经发国际会馆108.94319934.3263414昆仑银行108.94338734.3332655赛高街区108.95316134.3348746西安中学西门108.94463634.3387877运动公园东门108.95253234.3508658运动公园南门108.95023234.3475729管委会108.94552534.35336910出口加工区广场108.93607534.36264411鼎新花园108.96689634.32186212西安外国语学校108.95851534.32028913雅荷花园108.95413134.32076614御道华城108.95323334.31953615天地时代广场108.95464334.32522316市图书馆108.9550234.32669917文景观园108.9433634.32777918长庆电视台108.95974534.33291419移动公司108.95343934.33327220凤城五路108.95401434.33623821凤城六路108.95417634.33993422中登家园北门108.94028834.33911423万华园108.94459134.34282524粤华凤城家园108.93298534.3467325市人大108.96238634.34755726市委108.95901834.34751227政务大厅108.9561734.3478428首创国际城108.94543534.36229429中登大厦108.95328734.33023930运动公园北门108.9501734.353466附表2 部分已有网点位置距离计算表第一点经度1纬度1第二点经度2纬度2距离c1108.95295434.3248282108.94856234.323762392.53010540.9999937691108.95295434.3248283108.94319934.326341848.8766660.999970861108.95295434.3248284108.94338734.3332651244.9226520.9999373271108.95295434.3248285108.95316134.3348741117.2054110.9999495271108.95295434.3248286108.94463634.3387871708.1550330.999882011108.95295434.3248287108.95253234.3508652895.4123060.9996610031108.95295434.3248288108.95023234.3475722539.8400470.9997391481108.95295434.3248289108.94552534.3533693237.6408030.9995761361108.95295434.32482810108.93607534.3626444451.5393650.9991987611108.95295434.32482811108.96689634.3218621231.897140.9999386321108.95295434.32482812108.95851534.320289691.78979430.9999806471108.95295434.32482813108.95413134.320766462.64530720.9999913441108.95295434.32482814108.95323334.319536588.92214370.9999859751108.95295434.32482815108.95464334.325223150.54338050.9999990841108.95295434.32482816108.9550234.326699272.66219050.9999969941108.95295434.32482817108.9433634.327779882.10111330.9999685351108.95295434.32482818108.95974534.3329141070.4043670.9999536671108.95295434.32482819108.95343934.333272939.84661990.999964281108.95295434.32482820108.95401434.3362381271.9786560.9999345731108.95295434.32482821108.95417634.3399341682.9787020.9998854621108.95295434.32482822108.94028834.3391141924.3002890.9998502611108.95295434.32482823108.94459134.3428252126.1420750.9998172021108.95295434.32482824108.93298534.346732980.1192710.9996408791108.95295434.32482825108.96238634.3475572654.445650.9997150781108.95295434.32482826108.95901834.3475122575.7416590.9997317221108.95295434.32482827108.9561734.347842573.7401970.9997321391108.95295434.32482828108.94543534.3622944216.5499150.9992811111108.95295434.32482829108.95328734.330239602.34830390.9999853281108.95295434.