2013年数学建模A题-优秀论文

2013高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》（以下简称为”竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：A我们的参赛报名号为（如果赛区设置报名号的话）：A0114所属学校（请填写完整的全名）：参赛队员（打印并签名）：1.2.3.指导教师或指导教师组负责人（打印并签名）：（论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。）日期：2013年9月15日赛区评阅编号（由赛区组委会评阅前进行编号）：2013高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：车道被占用对城市道路通行能力的影响摘要车道被占用是指因交通事故、路边停车、占道施工等因素，导致车道或道路横断面通行能力在单位时间内降低的现象。车道被占用的情况种类繁多、复杂，正确估算车道被占用对城市道路通行能力的影响程度，将为交通管理部门正确引导车辆行驶、审批占道施工、设计道路渠化方案、设置路边停车位和设置非港湾式公交车站等提供理论依据。针对问题一，通过对视频1中交通事故发生至撤离期间,各种数据的采集，确定了事故横断面实际通行能力，运用Gibbs抽样仿真方法，通过excel软件解决数据缺失问题，并用eviews进行ARMA模型拟合，发现实际通行能力的变化过程为yt=−2.0232+ϵt−0.7243ϵt−1−0.2757ϵt−2。针对问题二，我们结合视频1，采用与问题一同样的方法对视频2进行数据采集和分析，得到了横断面的实际通行能力，并将视频1和视频2仔细对比，通过视频1和视频2中对交通事故发生至撤离期间,事故所处横断面实际通行能力，运用SPSS软件进行两独立样本的曼―惠特尼U检验，再根据路段附近交通设置，车辆流向比例，司机心理，周围地形等因素，分析出产生差异的原因是主要是在不同车道车辆流量的比例。针对问题三，通过分析视频1塞车情况，分别建立非稳态排队论模型和分段差分方程模型，并运用matlab软件编程绘制图像，解释出视频1中路段车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间的关系为：实际通行能力越大交通事故所影响的路段车辆排队长度越小，或者说路段车辆排队长度增长越慢；路段上游车流量越大交通事故所影响的路段车辆排队长度越大，或者说路段车辆排队长度增长越快；事故持续时间越长，交通事故所影响的路段车辆排队长度越大。针对问题四，利用问题三所建立的模型和题目所给数据，如果采用非稳态排队论模型，运用matlab软件求解，得出从事故发生开始，经过6.656分钟，车辆排队长度将到达上游路口；如果通过分段差分方程模型，事故发生时上游路口刚好是绿灯时需要经过6.42分钟，事故发生时上游路口刚好是红灯时需要经过7.97分钟，与上一模型所得结果相吻合。说明两模型有很高的合理性和实用性。【关键词】通行能力Gibbs抽样仿真ARMA模型拟合曼―惠特尼U检验非稳态排队论分段差分方程11问题重述车道被占用是指因交通事故、路边停车、占道施工等因素，导致车道或道路横断面通行能力在单位时间内降低的现象。一条车道被占用，也可能降低路段所有车道的通行能力。附件中视频1和视频2的两个交通事故处于同一路段的同一横断面，且完全占用两条车道。问题1：根据视频1（附件1），描述视频中交通事故发生至撤离期间,事故所处横断面实际通行能力的变化过程。问题2：根据问题1所得结论，结合视频2（附件2），分析说明同一横断面交通事故所占车道不同对该横断面实际通行能力影响的差异。问题3：构建数学模型，分析视频1（附件1）中交通事故所影响的路段车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间的关系。问题4：假如视频1（附件1）中的交通事故所处横断面距离上游路口变为140米，路段下游方向需求不变，路段上游车流量为1500pcu/h,事故发生时车辆初始排队长度为零，且事故持续不撤离。请估算，从事故发生开始，经过多长时间，车辆排队长度将到达上游路口。2模型假设1.视频中所统计数据真实可靠。2.排队所占车道车辆数与对应车道行驶方向车辆数成正比，即：车道一车辆数：车道二车辆数：车道三车辆数=0.21:0.44:0.35。3.除事故车辆外的其他车辆严格遵守交通规则，红灯停，绿灯行。4.车辆到达率与正在排队车辆数量无关，无论有多少车在排队，车辆到达率不变。5.车辆来源是无限的。6.堵车期间该路段没有其他交通事故发生。7.在堵车状况下相邻两辆车车头之间间距为7米。3符号说明•yt.............................................................t时刻视频1实际通行能力•ϵt...........................................................................t时刻的误差•I..............................................................................单位矩阵•i............................................................................虚数i2=−12•t..........................................................................堵车持续时间•λj................................................................j状态下的车辆到达率•µj................................................................j状态下的车辆离开率•Pj(t)...............................................................t时刻为状态j的概率•A...............................................................常微分方程组系数矩阵•xk.................................................................矩阵A的第k个特征值•e.....................................................................自然指数e≈2.718•ck.................................................................常微分方程通解系数•m(t)................................................................t时刻的平均车辆数•s(t)............................................................t时刻路段车辆排队长度•a1.............................................................小区1进出的车净到达率•a2.............................................................小区2进出的车净到达率•λ′1............................................................红灯时的上游路口到达率•λ′2............................................................绿灯时的上游路口到达率•P′j(t)........................................................t时刻状态j下的概率变化率•l0..........................................................平均车间距与平均车长之和•rk.........................................................第k个特征值对应的特征向量•Pi..........................................车道i所含排队车辆数占总排队车辆数比例•P{N(t+∆t)|N(t)}......t时刻状态为N(t)时，t+∆t时刻状态为N(t+∆t)的条件概率4模型的建立与求解4.1问题一4.1.1问题一的分析为分析事故所处横断面实际通行能力的变化，首先对视频1中交通事故发生至撤离期间的录像进行处理，数出以30s为时间间隔，红绿灯交替1min为一个周期时单位时间内横断面通过的车量个数，计算实际通行能力。由于视频中有多处出现片段的中断，影响计数，为了得到完整的数据，同时考虑到经过横断面的车辆数为随机数这一情况，我们选用Gibbs抽样仿真方法，对缺失的数据进行填补。用填补后的数据绘制折线图，描述交通事故发生至撤离期间,事故所处横断面实际通行能力的变化过程。34.1.2问题一的解答4.1.2.1数据的处理与分析：①对视频1中事故所处横断面通过车辆进行计数：对经过事故所处横断面的车辆进行计数，按照规定，只考虑四轮及以上机动车、电瓶车的交通流量，且换算成标准车当量数，根据《公路工程技术标准》（JTGB01-2003）规定的换算标准[1]，对视频中符合要求的车辆数据进行换算。具体换算规则如下：表1：四轮以上车辆折算规则车型载荷及功率折算系数小客车额定座位≤19座1.0大客车额定座位≥19座1.5小型货车载质量≤2吨1.0大型货车载质量≥2吨1.5从视频可知事故发生时间为16:42:32附近，所以我们以30s为时间间隔，以16:42:30为时间起点，开始计数。数据结果如下：表2：视频1车辆计数及当量转换时间小车大车当量转换通行能力时间小车大车当量转换通行能力42:307210120043:00819.5114043:309110.5126044:00819.5114044:3080896045:0080896045:30718.5102046:00718.5102046:30909108047:0070784047:30718.5102048:0011011132048:3010010120049:00909108049:30909108050:00718.5102050:3010010120051:00909108051:3080896052:009110.5126052:30819.5114053:0080896053:30909108054:00909108054:30819.51140