数学建模论文-深圳关内外交通拥堵探究与治理.doc

2014年全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： B 我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：洛阳师范学院 参赛队员 (打印并签名) ：1.王佳佳 2.张清娟 3.李艳 指导教师或指导教师组负责人 (打印并签名)：数学建模组 日期： 2014 年 8 月 25 日赛区评阅编号（由赛区组委会评阅前进行编号）：2014年全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：深圳关内外交通拥堵探究与治理摘要交通拥堵不仅是当今许多大城市日益严重的社会问题，同时也给国家和城市造成巨大的经济损失。对于不同的城市和地区其深层原因存在差异。本文以深圳市为例研究了深圳市关内外交通拥堵的深层原因并建立数学模型讨论如何治理深圳交通拥堵这个难题。针对问题一，以梅林关为例定义关口连接道路的拥堵指数，分析关口广场区域高峰期拥堵的直接原因及深层原因，并对进一步研究关口广场拥堵问题所需交通数据的采集侧重内容提出建议。道路拥堵评价指标有很多，为了能使评价客观，全面和科学，我们采用单位时间每千米内单一车道上运行的车辆数来判断道路的拥堵指数。考虑到信息的不完备性，我们有选择地选取了数据较完善的几个关口来计算拥堵指数，将拥堵的强度划分为5个等级，分别为：畅通、基本畅通、缓行、轻度拥堵、拥堵。得出各关口早晚高峰期的拥堵指数后，利用以时间为横坐标，交通拥堵指数为纵坐标绘图从中找出了道路拥堵的直接原因，通过进一步的讨论，最终得到造成拥堵的深层原因。最后，我们根据建立的模型参数，确定了进一步研究拥堵问题应侧重采集的数据。针对问题二，在不增加关内外通道情况下，能否通过调整城市分区功能，改变关口区域功能架构以及改善交通管控措施来缓解梅林、布吉等关口的交通拥堵。考虑到不同产业对从业人员的吸引力程度不同，为了定量分析这些数据，本文建立了吸引力模型。吸引力模型是考虑到总量中第一产业、第二产业、第三产业等因素对分区内从业人员数量的作用，可以认为更高的能吸引更多的从业人员。分别给第一、二、三产业赋权重得到了某区域总量这一数学模型，从而将不同分区的不同产业与人数建立起联系，根据不同吸引力指数的大小，提出了缓解关口拥堵问题的解决方案。 针对问题三，如果可以增加关内通道，应选在哪些地方增加通道。通过对前两个问题的分析，可知深圳市有些关口道路存在不同程度的拥堵，而有些关口道路却处于畅通状态，故本文利用分流疏通的方法，通过建设两条关内新通道以达到将拥堵路段的车流引到畅通路段的目的。增加通道的选址，需要考虑到城市的分区，关口拥堵指数和关内外人口分布等因素。关键词：交通拥堵指数 吸引力模型 吸引力指数 分流疏通 一、问题重述交通拥堵是目前中国各大城市面临的共同难题，但拥堵的成因各不相同，因而需要在摸清规律的基础上有针对性地提出解决方案。由于历史的原因，深圳由关内关外两个区域组成。关外由宝安、龙岗两个行政区和光明新区、龙华新区、坪山新区、大鹏新区四个功能区组成；关内含罗湖、福田、南山、盐田四个行政区。关外与关内由自然山丘隔开，由于有相当的一部分人口在关外居住，在关内上班，导致在上下班高峰期各关口进出通道经常成为交通最拥堵的地方，尤其以布吉关、梅林关等处为甚，在高峰期发生道路交通事故更会严重影响到广大市民的工作和生活。 为了解决这一长期困扰深圳发展的问题，政府在道路建设上投入了大量的资源。目前，主要关口道路的互联互通程度越来越高，直接增加了关口交通管控工作的复杂度。与此同时，大规模的基础设施建设也对交通信息采集设备的完好性和可靠性造成了不良影响，从而使关口交通管控和事故应急处理决策愈加困难。 因此，使用数学方法对不完整的交通信息进行建模分析，就成为定量分析关口交通特性及构成要素的重要手段。如果能在不断修正、调整的基础上取得较可靠的分析结果，将对制定有效、合理的交通管控及事故应对方案提供有益的帮助。 本题附件给出了交警部门记录的各主要关口进出通道瓶颈断面代表时段的交通流量、对应车速数据和行车道数，对拥堵严重的梅林关还以样本抽取方法给出了部分与关口广场连接道路对应时段的相关参考数据。请根据这些数据以及收集到的深圳城市功能分区规划、以及实际城市发展等方面的相关资料分析讨论以下问题： 问题一：分析造成各关口拥堵的深层原因。以梅林关为例，考虑信息不完备的影响因素构建关口交通模型，分析造成关口区域高峰期拥堵的直接原因，对关口广场各连接道路进行分类或定出拥堵指数，根据你的模型参数，给出今后进一步研究关口广场拥堵问题所需交通数据的采集侧重内容建议。 问题二：在不增加关内外通道的情况下，能否通过调整城市分区功能，改变关口区域功能架构以及改善交通管控措施来缓解梅林、布吉等关口的交通拥堵。 问题三：如果可以增加关内通道，试问应选在哪些地方增加通道（不考虑建设成本）。二、问题分析 对于问题一，本文用单位时间每千米内单一车道上的车辆数（假定每辆车的平均长度为4米）以及车速来描述交通拥堵指数。公式如下：其中为单位时间内的车流量，为车速，为车道数，为车长此处设为4米通过对交通拥堵指数的定义，对各个关口的各个时间段的拥堵指数进行算，利用绘图，分析图形进一步探讨造成关口拥堵直接原因和深层原因，再通过建立模型的参数给出了今后进一步研究关口广场拥堵问题所需交通数据的采集侧重内容建议。对于问题二，要在不改变交通道路布局的情况下，缓解交通压力，可以调整城市的功能分区，改变关口的区域功能构架，进一步改善关口与关外的交通管控措施来缓解关口的交通压力。本文通过查询地图，了解深圳的道路网布局及深圳市的未来规划，对深圳市的拥堵问题提出了些改善的建议。通过建立吸引力模型，考虑各分区的总量及深圳市关内各区的第一，第二，第三产业的，总量越大，越能吸引更多的从业人员，故可用各个区的总量代表每个区吸引的人数。，可以将总量的函数表示为：通过这个函数和本文搜集到的深圳市关内各区的本地生产总值及第一，第二，第三产业的总值，得到深圳市关内各区吸引的人数，结合道路网布局提出合理的解决方案。对于问题三，如果可以增加关内通道，应选在哪些地方增加通道，即增加关内通道的选址问题，鉴于只可增加关内通道，本文结合地图采用分流疏通的方法将拥堵路段的车流引到畅通路段，缓解关口路段的拥堵。三、模型假设（1）假设关口检验站不对车辆进行检查；（2）假设关口无红绿灯；（3）假设题目所给出的有效数据真实可靠；四、符号说明：单位小时内的车流量；：车速；：交通拥堵指数；：车道数；：车辆的平均长度；：车辆的平均长度；：总量；：吸引力系数；五、模型的建立与求解5.1问题一模型的建立与求解5.1.1交通拥堵指数模型的建立交通拥堵指数（）：单位时间每千米内单一车道上运行的车辆数，运算公式如下其中为单位时间内的车流量，为车速，为车道数，为车长此处设为4米。交通拥堵指数是一个连续的变量，定义取值范围为0到1，值越大则评价时段内的道路运行状态越差，拥堵强度越大；反之，道路运行状态越好，拥堵强度越轻。本文通过交通拥堵指数公式，结合附件所给出的数据计算出各个时间段的拥堵指数，再将拥堵指数进行划分，总共分为5个等级分别为：畅通、基本畅通、缓行、轻度拥堵、拥堵，表示如下：5.1.2问题一模型的求解（1）分析造成关口拥堵的直接原因本文根据表1中的数据随机抽取其中两个关口（丹砂路沙湾、松柏路）在给出的数据的7天内各小时的交通拥堵指数变化情况如图1、图2。图1 松柏路1月7号至13号路况图2 丹沙路1月7号至13号路况通过分析图形可以看出，拥堵指数变化大致呈,交通拥堵主要集中在每天7点左右和18点左右，其余时间较少。而7点左右和18点左右正是早高峰期和晚高峰期，由此可知造成高峰期关口拥堵的直接原因可能有：关外人员给道路带来了较大的车流量，超过了道路负荷极限；机动车辆使用增加，道路建设相对缓慢，导致有限的城市道路无法满足日益增长的交通需求；现行以功能分区为主的城市规划理念为拥堵埋下隐患；人们生活水平日益提升，导致城市私人汽车的数量迅速增加；人们对出行时间及路口的选择。（2）分析造成关口拥堵的深层原因根据表1中数据我们任选出其中两天多个关口的交通数据，分析其交通拥堵指数，得到的图形如图3，图4。图3 1月7号各路路况图4 1月11号各路路况通过分析图形可以看出，同一时间段不同关口的拥堵指数情况相差较大，尤其图中选出的几个关口中，107国道南头关口的交通拥堵指数明显大于其他关口。从这三个图中还能看出：在同一时间不同道路上的车流量不同，拥堵程度也不同即有些关口处于拥堵状态，而有些关口却非常畅通。所以造成高峰期拥堵的直接原因还有：城市布局、功能分区造成城市交通时间上、空间上的不均衡，即交通运行在个别路口上时间及空间上的过度集中。综上两点，可以推断造成各个关口的深层原因：关口道路建设失去平衡；关口道路网的布局形式缺陷；公交车道较少。（3）关口交通数据采集的侧重点分析通过对关口交通拥堵指数模型的建立，得出了深圳市交通拥堵的深层原因和直接原因。对这些原因的分析我们给出以下建议：1）加强对私家车数量增长的数据的采集。人们生活水平的提高导致人们的购买力也相应的提高，所以私家车的数量迅速增加，给交通运输带来了很大的压力，所以应当加强对私家车数量增长的数据的采集。2）增加对车流量以及车速的采集次数与采集区域的广度，车流量以及车速是进行交通拥堵指数模型分析的重点，如果能得到更加准确的和大量的关于它们的有效数据，那么得到的模型参数会更加的准确，这样就能根据这些参数提出相应的解决措施与方法。3）不同程度地去采集成品油价格的相关数据。油是整个交通系统的动力所在，如果成品油的价格增长过快、过高就会造成交通拥挤的短暂性地减缓，然而经过一段时间后，交通拥挤问题会又凸显出来。4）调查公共交通的客运输量，进一步研究公共交通在城市交通系统的作用，制定相应的措施。5）早晚拥堵时段的数据对研究车流量高峰很重要，应当作为重点数据采集的采集重点。5.2问题二模型的建立与求解5.2.1吸引力模型的建立（1）模型介绍 引力/空间互动模型最早由美国的提出，随后经过多方的修正与完善。该模型的创新主要体现在两个方面：第一，考虑了商店与顾客之间的互动。该模型认为顾客选择商店主要基于两个指标：商店的吸引力和商店的可达性。前者主要通过停车方便、价格等因素来衡量；后者则根据距离、交通等因素衡量。模型假设在同等吸引力情况下，顾客优先选择可达性更好的商店；同样，可达性相同时，顾客优选选择吸引力高的商店；但同时又考虑到了顾客行为的不确定性，认为优先选择并不一定是决定性的。第二，该模型考虑了商店选址对企业连锁网络上的其他商店以及竞争对手的影响，而不是只进行孤立的分析。通过分析商店潜在业绩的构成和来源，发现有多少是以其他商店的业绩损失作为代价的，有多少是从竞争对手那里赢得的，有多少是通过开发新市场获得的。通过这种分析，使得企业的选址决策更为合理。（2）吸引力模型的建立建立吸引力模型，需求出对各区人口数量有影响的吸引力指数，每生产一亿元需要的劳动力的人数是不一样的，可以将表示为总量的函数 （1） （2） （3） （4）根据对吸引力模型的分析，本文考虑到深圳市2012年的总量及关内各区的第一，第二，第三产业的总量等因素可能影响各区的从业人员数量，可以认为更高的能吸引更多的从业人员，进而就有更多的车流量通过该区，所以本文通过对深圳市关各区的三个产业对从业人员的吸引力，可以得到哪个区的人口数量较多，哪个区的人口数量较少，然后再通过调整各区的产业规划来适当地减少或增加人口数量，达到缓解交通压力的目的。通过查找资料得到深圳市关内各区的三个产业如下表3表3 各区三个产业的福田区罗湖区盐田区南山区全市第一产业（万元）899711012631100863018第二产业（万元）186485410720047498131667936357376373第三产业（万元）218511711250890529076971161586072061210生产总值（万元）2372502213582010365777328306231129500601利用将表三中数据绘制成条形图如下图5所示图5 深圳市关内各区生产总值分析图5可知关内各区经济以第二、三产业为主。福田以第三产业为主且第三产业是关内发展最好的，南山以第二产业为主且它的第二产业是关内所有区中最好的，罗湖盐田都以第三产业为主，但发展比较缓慢。 表4 2012年深圳市三个产业从业人数及总量统一表第一产业第二产业第三产业生产值（万元）630185737637372061210人口数（万人）0.1311373.2419397.824吸引力系数k（10(-6)）2.086.5050.552由表4可知深圳市从事第二、三产业的工作人员很多，第二、三产业创造的在深圳市总中也占很大比例。表5 深圳市关内吸引力统计表第一产业（万元）第二产业（万元）第三产业（万元）吸引人数（万人）福田区899718648542185117124.2131罗湖区110110720041250890513.8816盐田区26374981329076976.4834南山区110081667936311615860114.9374 由表5可知南山区吸引人数最多，福田区次之，说明南山区，福田区提供的工作条件和机会，更能吸引人。5.2.2问题二模型的求解（1）功能分区解决从表3可以看出，南山区具有最高的吸引力，而盐田区却有最小的吸引力，这说明各个分区的功能是有密切关系的。由于福田区位于深圳市特区的中心，在该区的三个产业中，第三产业占有较大的比重。位于东部的罗湖区生产总值也侧重于第三产业，可见罗湖区的功能重在第三产业。但最靠近东部的盐田区的经济发展则远远落后与其他三区。南山区中的三个产业中第二产业的生产总值较为显著，而且该区具有独特的地理位置：离珠江进，港口大，对外贸易进出方便，使该区吸引了大量从事第二产业的人。中部的福田区由于聚集了大量从事第三产业的群体，导致该区早、晚高峰通过梅林关的车辆多，因此该区应该削减第三产业的继续发展，可以增加东部罗湖区和盐田区的第三产业生产总值，增加该区的第三产业比重，可以有效地吸引从梅林关到福田区的车流，去东部的车流可以从清坪快速路入关，进而缓解布吉关的交通压力。由于西部的南山区独特的地理位置，并且是重要的工业集中地，所以应侧重发展该区的第一产业，积极鼓励第二产业的开发，从而吸引中部通过梅林关的车辆从107国道、广深高速、沙河路和白芒路、福龙隧道等处入关，这样可以缓解梅林关的压力。经济发展速度较慢的盐田区，第一，第二，第三产业都应加快其发展速度，吸引其他三个区的从业人口数量。（2）关口区域功能解决位于布吉关口的清水河是深圳市最大的物流中心，其中的清水河国际汽车物流产业园是深圳市最大的汽车物流中心，因此，清水河物流中心每天都会有大量的车辆出入，从而导致布吉关口因有大量车辆通行并造成高峰期的拥堵。于是，可以把清水河的物流产业适当地移到东部的盐田区和西部的南山区，利用东部快捷的交通以及西部南山区发达的工业，可以减少物流的成本，并可以有效的缓解布吉关口的布吉路、清水河联检站的交通压力。（3）交通管控措施上解决从关外到关内必须要通过梅林关和布吉关等关口才能入关，为了缓解梅关和布吉关的交通压力，本文列出一下措施：第一政府可以适当提高梅林关和布吉关的通关收费并减少其它关口的收费，利用价格调控作用，让一部分本来可以从梅林、布吉关出入关的车辆，绕道从其它收费低的关口入关。第二从题目中我们看到，公交车道很少，政府可提倡改变人们的出行方式或适当提高油价，减少私家汽车数量缓解交通拥堵。第三位于布吉关口的清水河每天都有大量的车辆出入，因此交警部门通过建立一些规定等，限制或减少清水河物流公司的车辆在早晚高峰时禁止出行。第四还应该提高城市道路突发状况的快速处理能力，当城市的某条道路上遇到车祸等一些突发事件，如果交通警察能够及时、快速地到达现场进行处理，将事故拖延的时间降到最低，进而就可以减少交通拥堵的时间。第五要建立和健全交通监控设施。对特别拥堵的路段可以利用移动通讯方式发给行驶中的出租车、大货车和私家车主。5.3问题三模型的建立与求解5.3.1分流模型的建立通过对前两个问题的分析，我们得到了关口拥堵的直接原因：车流量超过关口所能承受的容量。城市分区构架的不合理导致了车辆的分配不均匀，从而导致了有些道路存在不同程度的拥堵，而另一些关口道路却处于畅通状态。鉴于只能在关内增加通道的限制，我们利用分流疏通的方法，通过建设关内新通道以达到将拥堵道路段的车流引到畅通道路段的目的。 把深圳市关内分为东、中、西三个地区，对各个地区的1月7号到1月13号各道路早、晚高峰期出入关的车流量进行统计，统计如下： 表6 东部出入关的统计表早高峰车流量晚高峰车流量高峰车流量合计清坪快速清水河联检站出口北行-南-战略205716743731清坪快速清水河联检站入口南行-北-战略148510982583保洁路清水河检查站出-南-战略66813091977西环路清水河检查站入-北-战略443348791布吉路-深惠路南向北行驶350837867294深惠路-布吉路北向南行驶868.59321800.5丹沙路沙湾检查站入-北-战略8329061738沙湾路沙湾检查站出-南-战略315454769车流量合计10176.51050720683.5 从表6中我们可以看出布吉路进出关车流量较大，清平快速清水河次之，进一步我们可以推断出东部地区，布吉路、清平快速、丹砂路在高峰期可能会比较拥堵。表7中部出入关统计表早高峰车流量晚高峰车流量高峰车流量合计福龙隧道-南向北行驶432655519877福龙隧道-北向南行驶416142208381梅观公路南坪立交桥下北行-南-战略298735436530梅观公路普滨加油站南行-北-战略5639471010349车流量合计171131802435137 从表7中可看出，中部地区整体上车流量都很大，在中部地区交通拥堵尤为严重。表8 西部出入关统计表 早高峰车流量晚高峰车流量高峰车流量合计107国道(广深公路）南头检查站南行-北-战略5533459910132深南大道南头检查站出-南-战略445742998756G4(广深速)(2260km+969m)南行-北-战略274833656113广深高速同乐检查站出-东-战略296715974564沙河西路白芒关检查站出-南-战略91110651976松白路白芒检查站入-北-战略118112522433车流量合计177971617733974由表8可知，西部地区交通拥堵程度相差很大，其中107国道车流量最大。5.3.2模型的求解从表4、表5、表6可以看出东、中、西三个地区的车流量之间存在这巨大的差别，而造成这种现象的主要原因就是关内车道连通率不高，结合对第二题的分析，我们提出可以在关内修建两条路，如下图6示图6关内增加两通道 第一条路修建在北环路和南坪快速路之间，西起107国道，向东依次连接广深高速公路、福龙隧道、梅观路、209省道、清水河检查站和布吉关。优点有下两点： 可以增加关内道路的车通率，使得行驶道路多样化。原来从关口到关内市区只能走南北路，如果要走其它关口道路那么只能在市区里面选路，并且还导致大部分车辆没有余地去选择出关关口，通过修建此路可以将市区内几个主要交通线路连为一体。此路线西部可以连接深圳宝安国际机场、南山区工业园，东部可连接福田区行政中心和梅林关、布吉关等关口，大大增加了车辆出入关口的道路选择。 缓解了北环路等市区道路的通行压力。由于北环路非常靠近深圳市中心，只能疏导位于福田区和罗湖区的部分车流。原来从梅林关入关后只能走梅观路到市区，然后就只能走北环路到南山区或罗湖区，这样就大大增加了梅林关口、梅观路和北环路的交通压力，此路修建后，可以缓解北环路的交通压力。第二条路可以修建在深圳市东南部地区，此路起点在福田区的广深高速起点，东部可连接到沿河北路、丹平快速和360省道的交叉口，并且在罗湖区内可以高架桥为连接方式连接各道路与之联系。此路优点有以下三个： 缓解福田区和罗湖区市区内的交通压力，起到疏散区域交通的作用。由于大量车流通过梅林关进入关内，使得梅林关附近区域早晚高峰时的交通量非常大，修建了这条路后，福田区内交通线路更加优化，布局更加合理，因而能有效缓解和疏散区域交通压力。 与第一条形成环绕深圳市关内的大环路，这两条路将深圳市关内包围起来，形成了环区公路，有机的将关内交通线路连接起来，大大地方便了关内出行选择道路，并且形成环路网有利于使深圳市中心与各分区、郊区、市区外围相邻各区之间的交通联系。 增加靠近香港和海边沿岸部分的出关道路选择性。广深高速靠近深圳市最南边，并且从福田区市中心穿过，而这条路只能通向西边。通过修建该路后，延伸了，扩大了广深高速的东向通行能力。可以有效地缓解西区交通压力，还可以扩大东边罗湖区的通行能力，起到了很好的分流作用。六、模型的评价模型优点：（1）问题一模型通过定义拥堵指数，将附件中大量复杂数据转化成简单易懂的数据再将数据绘图，可以更直观的反应数据特点，便于分析交通拥堵的原因。（2）问题二吸引力模型，巧妙巧妙地将GDP与从业人数联系起来，采用了数理统计的方法，制定出了吸引力指数，为改变城市功能分区，调整产业结构做出合理的解释。（3）问题三分流疏通方法，综合考虑城市的分区，关口拥堵指数和关内外人口分布等因素，选取的通道较为合理。模型缺点：为了简化模型，在建立模型的时候建立了一些假设，这使得论文描述的情况与实际情况有一定的误差。参考文献1魏雪萍.Excel2003中文版从入门到精通M.北京：人民邮电出版社，2009.2刘卫国.MATLAB程序设计教程M.北京：中国水利水电出版社，2010.3李学文，李炳照，王宏洲.数学建模优秀论文精选与点评M.北京：清华大学出版社，2011.4雷功炎.数学模型讲义M.北京：北京大学出版社，2004.5朱道元.数学建模案例精选M.北京：科学出版社，2003.6姜启源，谢金星，叶俊.数学建模（第三版）M.北京：高等教育出版社，2003.7姜启源，谢金星，叶俊.数学建模（第四版）M.北京：高等教育出版社，2011.8叶其孝.大学生数学建模竞赛辅导教材M.长沙：湖南教育出版社，2000.9深圳产业情况分析报告，/p-315524042.html，2014,8,25.10市统计局，深圳市2012年统计年鉴，/xxgk/tjsj/tjnj/201304/t20130424_2130913.htm，2014,8,25.11深圳市港务管理局， /，2014,8,25.附录附录1 松白路和丹沙路在1月7号到13号的路况（1）q=27118714912014714329495612381349137912591055124815651591170417241677111096974153739829818917212712815130310001379144614351924107312761553140214231371650999856728531393291195138124127120274952128394413271337104013051607159116501727165611299787155313462611751241301151573209041366130813751267108612371527163816271784167011109277265803622591961221191081603319481426138513981322100912671610162117451719164411039158326354413362431781431231712527561044123211611113103610811308132313921347120293078772950739831725015312212717122354671496110499808631003120411861166117510057727075784603545464244043062282043528448418127998386928439221038111011311046842831747697583573479402330299297449896806911809113073987591188996191366583183275472564454241436422824327142082083752569085176382291195810361059109986084281866958658449838131129834444991078681782177165779789796410891143112283984375874061958643539730229432144892286680271084569584283899210921227110391787679569072568252953041434033145161979881883881972983689993010491074106296292183878470564063246043139035739155552876577275565977781789898010551003869851791728659;v= 58.4156.1954.6652.6548.7751.2651.7655.4458.8156.4455.4453.8456.1356.8557.1755.8457.4756.654.7651.3554.355.2759.0756.4456.7653.2752.3852.5754.2448.3551.2953.3756.945553.1556.1355.3155.5856.9955.9556.1956.7155.5454.555.8156.157.8158.7358.0955.9854.2951.4948.6852.1752.3455.558.5354.5354.1353.9656.2555.6256.8254.7954.8757.4653.9553.8252.354.5555.7358.5156.9156.2654.3550.9549.447.6947.2853.1458.3856.254.9353.4654.0254.7456.356.5554.9156.4853.452.1453.752.255.3858.2856.358.9160.7353.449.4250.951.0451.2756.2255.8454.3655.6154.6755.9855.6455.5656.2557.6454.4552.0852.953.7356.762.7958.6757.7456.4551.8153.7150.1347.0151.8555.9154.125556.3856.7756.9756.9856.2654.3956.7354.1452.4853.757.8758.762.2360.6362.5462.3963.5751.9953.9950.7951.354.8256.958.1657.8759.2359.6558.8259.3358.4259.0557.4156.2456.358.3461.6862.0348.9543.6938.8435.0728.7324.7421.3523.5333.0433.7335.5432.6237.4235.2637.433.2327.7930.3934.8436.4634.5837.9441.6347.444.2441.0137.9227.0725.4421.5618.8418.5328.1232.0435.534.1136.3931.8635.9132.9529.8729.4139.4240.9641.0441.1544.9148.7649.2344.0340.6834.6429.0725.7121.4423.5633.2330.5433.7133.1936.4732.2332.7134.629.3430.4835.0337.1735.2838.9543.4645.6247.2337.5836.5831.6925.322318.5620.4228.931.3234.232.9838.1129.2134.1432.9130.2627.8235.1637.7434.7340.0238.6841.4747.0138.6735.2831.0325.8921.7821.0120.3328.2830.3334.5531.2536.9831.2335.7536.3927.0232.0234.6736.6338.3241.7141.9142.9744.8639.3934.6531.9225.0921.1921.1824.6826.8631.6333.0333.2937.8334.8533.7532.1230.1833.2635.2233.336.6137.5939.4742.6839.4941.6436.2431.9626.522.8220.7823.7226.2335.8438.7541.442.842.1340.3139.5237.7536.6939.7140.5640.9942.3444.4747.44;n=44444442222222;l=4;for i=1:7 for j=1:24t(i,j)=(q(i,j)/(v(i,j)*n(i,1)/(1000/l); endendfor i=8:14 for j=1:24t(i,j)=(q(i,j)/(v(i,j)*n(i,1)/(1000/l); endendt（2）x=1:24;y1=t(1,:);y2=t(2,:);y3=t(3,:);y4=t(4,:);y5=t(5,:);y6=t(6,:);y7=t(7,:);plot(x,y1,-.,x,y2,-o,x,y3,-x,x,y4,-+,x,y5,-*,x,y6,-s,x,y7,-p)legend(1月7号,1月8号,1月9号,1月10号,1月11号,1月12号,1月13号)title(松白路路况)xlabel(一天24小时)ylabel(交通拥堵指数)（3）x=1:24;z1=t(8,:);z2=t(9,:);z3=t(10,:);z4=t(11,:);z5=t(12,:);z6=t(13,:);z7=t(14,:);plot(x,z1,-.,x,z2,-o,x,z3,-x,x,z4,-+,x,z5,-*,x,z6,-s,x,z7,-p)legend(1月7号,1月8号,1月9号,1月10号,1月11号,1月12号,1月13号)title(丹沙路路况)xlabel(一天24小时)ylabel(交通拥堵指数) 附录2 1月7号和1月9号各路段路况(1)q=19311341104879475786623985043607054865749476947485601602453775805595649284578429141393095218914326062922101591736991847317432243439307529503469387839784212439738663478275123511868180227118714912014714329495612381349137912591055124815651591170417241677111096974153739811075493235591811537151714001343107093110401265134414481671162798768053437124337226017415110711821355364680694795687983187710461084125114511463101490576450415468744450531494344634305344053474404444524614853953143353222541945464244043062282043528448418127998386928439221038111011311046842831747697583188115311132951881106021334051494244454204402238994556469743665699576651834985445440633638293313318693312121762057001831320634363775391938284215458549454349419943834083399434202691206825919612211910816033194814261385139813221009126716101621174517191644110391583263544114210367505363174124812951031913792660777975100810901335133276857443036728438330919315414216925958768881281195388682798510311116127313771278118395993868916514172617066200430505448391405331398431447430459380317300344283224586435397302294321448922866802710845695842838992109212271103917876795690725;v=52.7352.8553.0550.0447.1949.1249.7941.2437.1947.7343.7746.1544.9943.7740.0945.8141.6340.9937.8935.5938.5537.9541.6647.1282.1779.4773.686865.362.3270.6880.1380.4278.5276.9276.0275.4976.4174.773.0771.8971.0765.1567.2469.7572.1473.0776.6658.4156.1954.6652.6548.7751.2651.7655.4458.8156.4455.4453.8456.1356.8557.1755.8457.4756.654.7651.3554.355.2759.0756.4454.551.4147.0646.0535.6334.439.7859.2556.9353.9853.2952.7550.2153.6655.0455.0653.9354.453.1649.0248.351.2854.347.3350.5751.6150.9150.6849.0249.241.3640.3240.8138.8238.3239.339.5639.0540.4739.0640.2938.529.2630.1739.8841.9843.3745.1434.4733.0136.6935.8136.0435.5125.8823.5427.0126.4623.4526.4625.5923.5524.0724.0920.8923.8619.8323.2724.9126.1133.8935.8448.9543.6938.8435.0728.7324.7421.3