2013全国数学建模大赛a题优秀论文.pdf
车道被占用对城市道路通行能力的影响车道被占用对城市道路通行能力的影响摘要随着城市化进程加快，城市车辆数的增加，致使道路的占用现象日益严重，同时也导致了更多交通事故的发生。而交通事故发生过程中，路边停车、占道施工、交通流密增大等因素直接导致车道被占用，进而影响了城市道路的通行能力。本文在视频提供的背景下通过数据采集，利用数据插值拟合、差异对比、车流波动理论等对这一影响进行了分析，具体如下：针对问题一，首先根据视频 1 中交通事故前后道路通行情况的变化过程运用物理观察测量类比法、数学控制变量法提取描述变量（如事故横断面处的车流量、车流速度以及车流密度）的数据，从而通过研究各变量的变化，来分析其对通行能力的影响。而视频 1 中有一些时间断层，我们可根据现有的数据先用统计回归对各变量数据插值后再进行拟合，拟合过程中利用残差计算值的大小来选择较好的模型来反应各变量与事故持续时间的关系，进而更好地说明事故发生至撤离期间，事故所处横断面实际通行能力的变化过程。针对问题二：沿用问题一中的方法，对视频2 中影响通行能力的各个变量进行数据采集，同样使用 matlab 对时间断层处进行插值拟合处理，再将所得到的的变化图像与题一中各变量的变化趋势进行对比分析，其中考虑到两视频的时间段与两视频的事故时长不同，从而采用多种对比方式（如以事故发生前、中、后三时段比较差值、以事故相同持续时间进行对比、以整个事故时间段按比例分配时间进行对比）来更好地说明这一差异。由于小区口的位置不同、时间段是否处于车流高峰期以及1、2、3 道车流比例不同等因素的影响，采用不同的数据采集方式使采集的变量数据的实用性更强，从而最后得到视频 1 中的道路被占用影响程度高于视频 2 中的影响程度，再者从差异图像的变化波动中得到验证，使其合理性更强。针对问题三：运用问题 1、2 中三个变量与持续时间的关系作为纽带，再根据附件 5中的信号相位确定出车流量的测量周期为一分钟，测量出上游车流量随时间的变化情况，而事故横断面实际通行能力与持续时间的关系已在 1、2 问中由拟合得到，所以再根据波动理论预测道路异常下车辆长度模型的结论，结合采集数据得到的函数关系建立数学模型，最后得出事故发生后，车辆排队长度与事故横断面实际通行能力、事故持续时间以及路段上游车流量这三者之间的关系式。针对问题四：在问题3 建立的模型下，利用问题4 中提供的变量数据推导出其它相关变量值，然后代入模型，估算出时间长度，以此检验模型的操作性及可靠性。关键词：通行能力车流波动理论车流量车流速度车流密度1一、 问题重述车道被占用是指因交通事故、路边停车、占道施工等因素，导致车道或道路横断面通行能力在单位时间内降低的现象。 由于城市道路具有交通流密度大、 连续性强等特点，一条车道被占用，也可能降低路段所有车道的通行能力，即使时间短，也可能引起车辆排队，出现交通阻塞。如处理不当，甚至出现区域性拥堵。车道被占用的情况种类繁多、复杂，正确估算车道被占用对城市道路通行能力的影响程度，将为交通管理部门正确引导车辆行驶、审批占道施工、设计道路渠化方案、设置路边停车位和设置非港湾式公交车站等提供理论依据。视频 1（附件 1）和视频 2（附件 2）中的两个交通事故处于同一路段的同一横断面，且完全占用两条车道。请研究以下问题：1.1. 根据视频 1（附件 1） ，描述视频中交通事故发生至撤离期间,事故所处横断面实际通行能力的变化过程。2.2. 根据问题 1 所得结论，结合视频 2（附件 2） ，分析说明同一横断面交通事故所占车道不同对该横断面实际通行能力影响的差异。3.3. 构建数学模型，分析视频 1（附件 1）中交通事故所影响的路段车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间的关系。4.4. 假如视频 1（附件 1）中的交通事故所处横断面距离上游路口变为 140 米，路段下游方向需求不变，路段上游车流量为 1500pcu/h,事故发生时车辆初始排队长度为零，且事故持续不撤离。请估算，从事故发生开始，经过多长时间，车辆排队长度将到达上游路口。注：只考虑四轮及以上机动车、电瓶车的交通流量，且换算成标准车当量数。二、问题分析（一）题1、2 中都是比较从交通事故发生到撤离期间，事故所处横断面实际通行能力的变化过程，因而（1）对问题1 首先设定在外部因素不变的情况下，2,3 车道被占用，实际上可视为局部缩短了道路宽度。而车道通行能力为每秒通过横断面的最大车流量。故为能体现出这一变化，我们选取：横断面车流量，道路的车流速，横断面前的车辆密度来描述。 （而这些变量都是在事故这一段时间内发生变化，所以可与事故时间段建立函数关系） 。（2） 对问题 2 也是沿用第一题的做法对事故横断面处的道路通行能力进行估计运算（其中车道占用变为 1,2 号车道） 。交通事故对道路行车造成的影响，不仅跟事故的严重性本身有关，而且跟事故发生的时间和地点有密切的关系。此问主要为研究占用不同车道对城市道路通行能力的影响。（二） 对问题 3： 从视频一可以直接反映出车辆的排队长度与事故持续时间的变化情况，而事故处横断面的通行能力与事故持续时间已在问题1 中得到解答（即附件1 中的拟合函数） ，从而通过时间建立通行能力与排队长度的关系。并且视频 1 中可观测出上游车流量与时间的关系，同理能够得出上游车流量与排队长度的关系。综上所述，我们再利用车流波动理论推导出车辆排队长度与道路通行能力、事故持续时间、上游车流量三者间的关系式。（三）对问题 4 当交通事故所处横断面距离上游路口变为 140 米时，下游方向需求不变，上游路段车流量为1500pcu/h,发生事故时车辆排队长度为 0，代入问题3 的模型中可预估到车队长度。2三、模型假设1、假设路口不发生事故，且所有司机都遵守交通规则。2、假设路面状况良好。3、信号灯转为绿灯时，车发动时间忽略不计。4、交叉口无人行道。5、只考虑红绿灯对车辆的影响。6、上游路口到下游到路口为纵向道路段。7、假设一定时间内不同车型的速度及性能是一致的。8、假设视频 1 与视频 2 发生在不同时段，不影响该道路的实际流量。四、符号说明1、pcu标准车当量数2、v车辆流速3、k车辆密度 4、kj阻塞密度值为110pcu/km 5、ks1事故处的车密 6、k1路段上游的车密 7、s1为事故处的通行能力 8、q1上游的车流量 9、t1为持续时间五、模型的建立与求解（一）流程图模型分解通过视频 1、2，得到通行能力的 3 个变量与时间关系车辆排队长度再分别作出车流量、速度、车密度与时间的关系图事故持续时间由这些关系式总结出实际道路通行能力与时间的关系，而得出表达式得到车辆排队长度与的关系上游车流量可由视频测量，得出关系表达式3（二）问题的求解分析1、问题一的求解分析通过观察视频 1 我们发现：每分钟的00s-30s 时车辆进行纵向通行，其他时间车辆进行横向通行，这便于我们更好地确定观测周期为 1 分钟。在外部因素（车道宽度、侧向净空、沿线状况）不变的情况下， 2,3 车道被占用，实际上可想为局部缩短道路宽度。而车道通行能力为每秒通过横断面的最大车流量。故为能体现出这一变化，我们选取：横断面车流量，道路的车流速，横断面前的车辆密度来描述。 （而这些变量都是在事故这一段时间内发生变化，所以可与事故时间段内建立函数关系）车流量：通过视频 1 观测出事故发生处横断面的车流量随时间的变化，见附件 1.根据 excel 表格视频 1描绘出的各时间段车流量图如下：时间段内通过事故横断面的车流量2520151050车流量(pcu)车总流量16:39:00-16:41:00-16:43:00-16:45:00-16:47:00-16:49:00-16:51:00-16:53:00-16:55:00-16:57:00-16:59:00-17:01:30-时间17:03:00-再利用 matlab 插值作出图像如图 a。图 a由图 a 可得:事故前车流量较为平缓，之后波动较为明显，说明车流量受影响比较严重，然后撤离后又渐渐恢复平缓。4再又绘出总流通量随时间的变化曲线如下：图（1）对上图总流通量的拟合函数进行求导得到事故横断面每分的流通量曲线图如下：图（2）由图（2）可得： 在事故发生的这一短时间内，车流量的变化忽高忽低，不同于正常时的规律性，当车流量低时乃是因为车道被占用所造成的堵车拥挤时车速减缓，在那一时间段内通行受阻而进行的测量，得到的值就会很小，而当此受阻结束后，一旦通行畅通一点后再次测量时得到的值就会很大，由于之前的高密度车密使其在后面时间内的通行量大增使其通行能力上升，而这种反复的降低上升现象中的时间段数也会因事故持续时间的延长而变得紧凑；因此事故其间的通行能力的变化规律为：先是下降一段时间后就会有小段时间回升，可能超过正常能力值；上升后又会持续循环上述规律且总体上看下降的时间段会越来越频繁随着事故持续时间的增多。5最后利用几何画板对上图的每一时间段选取最大流通量作为估计此段时间的道路通行能力的评定。分别做出它们的图如下：图（3）进一步说明道路通行能力的变化情况。道路车流速：测出事故前后与事故中的各时间段车流速的变化情况：横断面前240米作为测速距离得到：6时间段时间间隔（s） 速度（km/h）16:40:02-16:40:232141.1428571416:40:02-16:40:323028.816:40:10-16:40:403028.832.5661538516:40:12-16:40:402830.8571428616:40:16-16:40:422633.23076923事故前16:41:01-16:41:181750.8235294116:41:02-16:41:222043.216:41:04-16:41:292534.5641.7167058816:41:07-16:41:25184816:41:09-16:41:36273216:42:37-16:43:153822.7368421122.716:43:15-16:44:105515.7090909121.3545454516:44:07-16:44:39322716:44:53-16:45:192633.2307692333.2316:45:19-16:46:145515.7090909115.716:46:16-16:46:52362425.0909090916:46:52-16:47:253326.18181818事故中16:47:29-16:48:275814.8965517214.916:48:37-16:49:406313.7142857113.716:50:06-16:50:554917.6326530617.616:51:01-16:52:511107.8545454557.916:53:01-16:54:28879.9310344839.916:55:03-16:56:036014.414.416:58:04-16:59:187411.6756756811.716:59:18-17:01:377910.9367088610.9撤离后17:01:29-17:03:221147.5789473687.6视频一中车速的变化情况6050辆（pcu）403020100系列116:40:1616-:411:6:401602:4-1:421:6:1607-41:2:412:6:421637-1:416254:40167-3:1:415:6:51619-44:3:416:6:4961652:4-1:148:6:1637-47:516251:409161:5-1:405:6:1603-52:5:519:6:51186-1:037:01:37时间图 b7车速在事故发生时出现一段时间下降，后又有小幅度上升，总体上一直呈现下降趋势随着事故时间的增加。根据视频 1 中选取 240 米作为实验，再测量时间，根据速度时间公式得到各个时间段的速度，见附件 1；横断面前的车密度： （车密度=车流通量/车流速）的变化情况图如下：图 c在事故发生前车流密度的幅度变化较小，而在事故期间，在图b 中速度降低的影响下，车流的堆积，使其在事故后密度持续增加，中途有小段时间由于车流的暂时疏通使其车密有小段减缓。2、问题二的求解分析对第二题分析也是沿袭第一题的做法对事故横断面处的道路通行能力进行估计运算（其中车道占用变为 1,2 号车道） ：车流量：起始时间终止时间小车大车Pcu（min）事故前17:29:0017:30:001801817:30:0017:31:002012217:31:0017:32:0017:32:0017:33:0017:33:0017:34:0017:34:0017:35:0017:35:0017:36:002322717:36:0017:37:001622017:37:0017:38:0019325事故中17:38:0017:39:002112317:39:0017:40:001532117:40:0017:41:001622017:41:0017:42:001912117:42:0017:43:002312517:43:0017:44:001932517:44:0017:45:0014116817:45:0017:46:0017:47:0017:48:0017:49:0017:50:0017:51:0017:52:0017:53:0017:54:0017:55:0017:56:0017:57:0017:58:0017:59:0018:00:0018:01:0018:02:0018:03:0017:46:0017:47:0017:48:0017:49:0017:50:0017:51:0017:52:0017:53:0017:54:0017:55:0017:56:0017:57:0017:58:0017:59:0018:00:0018:01:0018:02:0018:03:0018:04:0020161417222017191620171517141317182225022201212124223212420201821222221212022212321181921202633图 5视频2中观测的事故横断面的车流量变化（一分钟的时间段内）35302520151050车流量（pcu）pcu17:2179:0:30171:0:30173:0:30175:0:30177:0:30179:0:40171:0:40173:0:40175:0:40177:0:40179:0:50171:0:50173:0:50175:0:50177:0:50189:0:00181:0:003:00时间段图 d道路车流速：起始时间终止时间17:29:0017:29:2817:29:0317:29:25事故前17:29:1017:29:4017:30:0017:30:2417:30:0817:30:31时间 （s） 速度（km/h） 每分内的流速2822302423930.8571439.2727328.83637.5652234.4527732.9766217:30:1417:34:2217:35:0417:35:2717:36:1417:36:3017:36:5017:37:1717:38:0717:38:2417:39:1317:39:1417:39:2717:40:20事故中17:40:0817:41:0317:41:1917:42:0317:42:0517:42:3017:30:4317:35:0817:35:2917:36:1117:36:3417:36:5117:37:1917:37:4017:38:2617:38:5017:39:2717:39:3017:39:5817:40:4017:40:2217:41:1917:41:4117:42:3317:42:3117:43:0729462544202129231926141631201416223026371525221928141029.793118.7826134.5619.6363643.241.1428629.793137.5652245.4736833.2307761.714295427.8709743.261.714295439.2727328.833.2307723.3513557.634.5639.2727345.4736830.8571461.7142918.7827.0981842.1714333.6791639.3522347.8617552.4571446.6363628.4607117:43:0917:43:2617:43:2517:43:5017:43:5317:44:1517:44:1717:44:3617:44:3817:45:0617:45:0717:45:2146.0838.5345244.3571417:45:2417:45:5617:46:1517:46:3417:46:4317:47:0217:47:1417:47:2817:47:2917:47:5817:48:1317:48:4517:49:0017:49:4017:49:4017:50:2832191914293240482745.4736845.4736861.7142929.79312721.61815.157891216.3018913.512.8955211.6756814.46.8571437.57894713.2923115.7090920.0930218.7826122.7368431.6856719.815.151216.313.512.911.714.46.867.5813.315.720.118.822.745.753692745.4736817:50:3017:51:235717:51:2317:52:357217:52:3517:53:285317:53:2917:54:336417:54:3417:55:416717:55:4118:56:557417:57:0919:58:096017:58:3918:00:4512618:01:1518:03:0811418:02:5518:04:006518:03:0718:04:025518:03:2018:04:0343撤离后18:03:2518:04:114618:03:2818:04:0638平均速度表 a11706050403020100km/h17:29:0017:29:1017:30:0817:34:2217:35:2717:36:30706050403020100横断面前的车密度：时间视频2事故时间段内速度变化视频2整个时间段速度（km/h）图 f图 e1217:3417:2:32617:1:34177:1:37179:1:43170:2:40171:1:49172:3:40173:5:43175:0:47176:4:43178:1:53017:3:50317:2:597:1809:02:55速度（km/h）系列1对 1,2 题的各变量进行比较得到：1.在车流量（图a 与图 d 比较得到）上：视频一从16:42:32 到 17:01:21 这段时间上事故堵车造成车流量与之前形成的变化高于视频二 1,2 道被占用所形成的变化，速度（图 b 与图 d 比较得到）上：在视频一中占据了 2,3 车道比视频二占距 1、2 车道对事故前车流速造成的影响大（其中考虑道两视频的时间段也不同，但总体上还是可以评估的，因为在一车道上还有两小区路口及三个车道流量的比列也不同可以相互调节，再者又是在城市小道路， 所以各种因素相互调节下， 采集的数据可以算作衡量数据使用）车密度（图 c 与图 f）的变化：视频一在前期时间段上车密的变化幅度明显高于视频二；但随着视频二的时间进一步加长，密度的堆积有一些升高，但就比较上还是低于视频二中的影响变化。所以综上所述可以得到：对于车道占用的不同影响到道路通行能力，关键是在外侧车道的车辆流通量较大，故而影响较为明显，而内侧车道的车流量较小对占用时在影响上与正常情况下通行能力的影响较小。3、问题三的求解分析对第三题进行分析：从视频一可以看出车辆的排队长度与事故持续时间的的变化情况，事故的横断面的通行能力关系可有1 题中的时间关系与排队长度的关系。而从交通信号相位中选取一分钟作为周期进行采集数据得出上游车流量与时间的关系，故而能够得出路段车辆长度与各变量的关系如下：上游车流量与持续时间的表格图像如下： （大车=2pcu，小车=1pcu）起始时间终止时间小车数大车数Pcu/min16:39:0016:40:00911116:40:0016:41:001511716:41:0016:42:001632216:42:0016:43:0016:43:0016:44:001411616:44:0016:45:001701716:45:0016:46:00151171316:46:0016:47:0016:48:0016:49:0016:50:0016:51:0016:52:0016:53:0016:54:0016:55:0016:47:0016:48:0016:49:0016:50:0016:51:0016:52:0016:53:0016:54:0016:55:0016:56:0011211819212316181413010101220011231821212520221413图 9车辆长度=a 道路通行能力+b 事故持续时间+c 上游车流量；道路通行能力 s 与持续时间 t 的关系：145138289223877yx x x x31.4单位为（pcu/min）（a）50050050002500上游车流量 q 与持续时间 t 的关系为：y 0.1x53.6x461.2x3453.4x21421.3x 401.7本 文再采 用 G R E E N S H IE L D 流一密模 型 , 如 图 4所 示 , 并规定 需求流量 q l 属 于 高速低 密 的 畅流 态 而s1属 于低 速高密的拥挤态。假设当交通事故发生后, 本车道上游的需求流量下降为q1, 对应的密度记为k1, 瓶颈点的通行能力下降为s1, 车流密度相应地上升为ks1 , 事故持续时间为t1, 故障排除后, 排队车辆以饱和流率s驶出, 对应密度记为ks。14y (t )事故发生点ks1s1k1q1k1q1事故影响区段0 t1 A Dt R Cy车流阻塞-消散过程的波形时-距图一般异常事件持续时间的定义是指从交通异常产生到交通流状态恢复正常所需的时间。它由4 个阶段构成, 第1 阶段是交通异常事件产生到AID系统检测并确认事件;第2 阶段是响应阶段 , 即从确认事件到救援车辆到达事发现场 ; 第3 阶段是清除时间 ,即从救援车辆到达到离开现场; 第4 阶段是交通流恢复阶段,即从事件清除到排队完全消散, 交通流恢复正常。这里的事故持续时间是指前3 个阶段的总时间,也可称为事故清除时间。流密关系曲线图其中得到结论：15波速：OBs1q1ks1k1，由参考文献 2 可以得到以下结论：k kOB uf(1s11)k12kdqh(ks1) |k ks1 uf(1s1)dkkjk1ks1从而得出车辆长度的与三变量的关系如下：t1(kjks1k1)(t t1)2(s1q1)(t t1)（2 s1q1）y= (1)(kj ks1k1)k1ks1kjks1k1tB(kj2ks1)t1其中s1为事故处的通行能力，q1为上游的车流量，kj为阻塞密度值为 110pcu/km；ks1为事故处的车密，k1路段上游的车密，t1为持续时间。(根据上游车流量与车速的关系可知具体车流量下的车密值，因为车流量与车速在事故发生期间都是事故持续时间的函数)4、问题四的分析求解当交通事故所处横断面距离上游路口变为 140 米时，而下游方向需求不变，上路段车流量为 1500pcu/h,发生事故时车辆排队长度为 0；车辆长度 y 的最大值实际上为 140 米排队到上游路口所以对问题三中的（）式对时间 t求导 dy/dt=0 得到车辆排队最大长度与事故横断面实际通行能力、事故持续时间、路段上游车游量的关系为：ymaxh(ks1) h(k1)2(tBt1)/4h(k1) (2)2kdqh(ks1) |k ks1 uf(1s1)dkkjtB(kj2ks1)t1k1ks1s1q1kj(s1q1)ks1k1得到：ufks1 k1(ks1k1)(kjks1k1)1kjuf(1ymaxt (k 2ks1)2ks12k)uf(11)2(1j)kjkjk1ks1 (3)2k14uf(1)kj题中给出上游车流量 q1 为 1500pcu/h车流量 q1 与车密(k1)间的关系式：q1=1/10*t5-18/5*t4+306/5*t3-2267/5*t2+14213/10*t-4017/10速度与时间的关系：v=-3/2500*t5+13/250*t4-83/100*t3+61/10*t2-23*t+60)16k1=(1/10*t5-18/5*t4+306/5*t3-2267/5*t2+14213/10*t-4017/10)/(-3/2500*t5+13/250*t4-83/100*t3+61/10*t2-23*t+60)q1 上游流通量为 1500pcu/h 时对应的车密 k1 的值为为 13,57，再计算时考虑到函数关系试为实测数据拟合而得，所以在计算上采用先将上游车流量与持续时间的关系图中进行预测得到其时间大致在 7.337min 左右，而后将时间代入到（3）中进行磨合与车辆长度的 140 米真实值进行比较，最后选取得到其值大致在6min 到 8min 中，这与经过三式直接得到的 6.45min 相接近。实际道路通行能力 s 与事故持续时间 t 的关系为公式（a）事故持续时间与 ks 的关系：ks=(1/500*t4-51/500*t3+8289/5000*t2-23877/2500*t+8831615164268885/281474976710656)/(-3/2500*t5+13/250*t4-83/100*t3+61/10*t2-23*t+60)将以上数据代入到（3）式得到时间 t=6.45min。六、模型的评价与推广该模型根据现实视频采集的实际数据， 所以在反映真实性上可靠性强， 实际操作性上也更有用，而且该模型的建立过程中采用大量的数据采集与分析方法使其数据实用性更强，并且在影响变量上从一开始就考虑较为全面，拥有多种变量来体现通行能力的变化，各变量之间又相互印证使其准确性更高。且在建立模型时也是在原有模型的基础上进行复合加工使其在理论上更实用；但在采集数据时虽说运用多种方法采集了多种数据，但难免由于人工采集和方法的先进性不足使其数据存在一定的误差导致计算结果不够准确，同时选择原始模型作为基础难免会有一些不适用之处使得模型存在一些误差。但该模型在对城市的车道占用问题上还算比较实用，故可考虑推广到对城市各个路段进行预测进而推广到对整个城市的道路进行多路预测。17七、参考文献【1】臧华,彭国雄.高速道路异常状况下车辆排队长度的预测模型J.交通与计算机【2】姜启源,谢金星,叶俊.数学模型M（第四版),高等教育出版社【3】胡守信,李柏年.基于 matlab 的数学实验M,科学出版社【4】王能超,计算方法算法设计及其 matlab 实现M,高等教育出版社【5】郭冠英,邹智军.道路阻塞时的车辆排队长度计算法J,中国公路学报【6】样少辉.固定式配时车辆延误与排队长度计算模型研究D,吉林大学附录1：1、车流量叠加后拟合的函数值：x1=1:3;x2=5:17;x=x1,x2,23;y=13 32 52 73 92 106 125 141 161 179 198 214 232 251 272 292 303;clfx3=1:25;y3=interp1(x,y,x3,cubic);plot(x,y,*,x3,y3,linewidth,5)a1=polyfit(x3,y3,5);z1=polyval(a1,x3);plot(x3,y3,x3,z1,g,linewidth,5)xlabel(time(min)ylabel(pcu/min)MSEa1=sum(y3-z1).2)2、上游车密度：y1=660 1020 1320 960 1020 1020 660 1380 1080 1260 1260 1500 1200 1320 840 780;y2=32.56615385 41.71670588 22.7 21.35454545 33.23 15.7 25.09090909 14.9 13.71817.6 7.9 9.9 14.4 11.7 10.9 7.6;x=1:16;y=y1./y2;x3=1:0.1:16;y3=interp1(x,y,x3,cubic);plot(x,y,*,x3,y3,linewidth,5)title(上游车密度)3、事故后速度的变化：y1=660 1020 1320 960 1020 1020 660 1380 1080 1260 1260 1500 1200 1320 840 780;y2=32.56615385 41.71670588 22.7 21.35454545 33.23 15.7 25.09090909 14.9 13.717.6 7.9 9.9 14.4 11.7 10.9 7.6;x=1:16;y=y1./y2;x3=1:0.1:16;y3=interp1(x,y,x3,cubic);plot(x,y,*,x3,y3,linewidth,5)title(上游车密度)4、上游车流量变化x1=1:3;x2=5:17;x=x1,x2,23y=13 19 20 21 19 14 19 16 20 18 19 16 18 19 21 20 21;clfx3=1:25;y3=interp1(x,y,x3,cubic);plot(x3,y3,*,linewidth,5)xlabel(time)ylabel(pcu)title(bianhua)5、视频1，拟合后的函数：x=1:0.1:25;y=0.0004*x.5-0.0255*x.4+0.5526*x.3-4.7754*x.2+31.3762*x-14.6281;plot(x,y,linewidth,5)title(y=0.0004*x5-0.0255*x4+0.5526*x3-4.7754*x2+31.3762*x-14.6281)xlabel(time(min)ylabel(pcu)y1=diff(y);x=1:0.1:25;y1=1/500*x.4-51/500*x.3+8289/5000*x.2-23877/2500*x+8831615164268885/281474976710656;plot(x,y1,linewidth,5)title(y1=1/500*x4-51/500*x3+8289/5000*x2-23877/2500*x+8831615164268885/281474976710656)xlabel(time(min)19ylabel(pcu/s)6、道路通行能力离散点：x=1:25;y=30.00,22.36,16.58,13.57,12.29,12.87,14.19,15.65,17.03,18.30,19.07,19.33,19.33,19.05,18.14,16.70,14.80,12.51,10.26,7.91,6.08,4.85,5.91,9.28,14.76;plot(x,y,*,linewidth,5)xlabel(time(min)ylabel(pcu/min)title(道路通行能力的离散点)7、速度与时间的拟合函数y=32.56615385 41.71670588 22.7 21.35454545 33.23 15.7 25.09090909 14.9 13.717.6 7.9 9.9 14.4 11.7 10.9 7.6;x=1:16;plot(x,y)x1=1:0.1:16;y1=interp1(x,y,x1,cubic);plot(x,y,*,x1,y1,linewidth,5)a1=polyfit(x1,y1,7);z1=polyval(a1,x1);plot(x1,y1,r-,x1,z1,linewidth,5)MSEa1=sum(y1-z1).2)8、车密度变化图；y1=660 1020 1320 960 1020 1020 660 1380 1080 1260 1260 15001200 1320 840 780;y2=32.56615385 41.71670588 22.7 21.35454545 33.23 15.7 25.09090909 14.9 13.717.6 7.9 9.9 14.4 11.7 10.9 7.6;x=1:16;y=y1./y2;x3=1:0.1:16;y3=interp1(x,y,x3,cubic);plot(x,y,*,x3,y3,linewidth,5)title(车密度的变化图)xlabel(time)ylabel(pcu/km)9、整个过程车流量变化值：x1=1:3;x2=5:17;x=x1,x2,23y=780 1140 1200 1260 1140 840 1140 960 1200 1080 1140 960 1080 1140 1260 12001260;plot(x,y)x3=1:0.1:25;y3=interp1(x,y,x3,cubic);plot(x,y,r,x3,y3,linewidth,5)20title(整个过程车流量值的变化)xlabel(time)ylabel(pcu/h)视频1车流量注：小车=1pcu，大车=2pcu编号起始时间终止时间小车大车pcu/minpcu/hpcu/min 叠加值116:39:0016:40:001111378013事故前216:40:0016:41:0017119114032316:41:0016:42:0014320120052416:42:0016:43:0052516:43:0016:44:0017221126073616:44:0016:45:0017119114092716:45:0016:46:0014014840106816:46:0016:47:00171191140125916:47:0016:48:00141169601411016:48:0016:49:002002012001611116:49:0016:50:001421810801791216:50:0016:51:00171191140198事故中1316:51:0016:52:00160169602141416:52:0016:53:001611810802321516:53:0016:54:001521911402511616:54:0016:55:001912112602721716:55:0016:56:001812012002921816:56:0016:57:001916:57:0016:58:002016:58:0016:59:002116:59:0017:00:002217:00:0017:01:002317:01:3017:02:00172211260303撤离后2417:02:0017:02:302517:02:3017:03:00视频1数据240米速度编号时间段时间间隔（s） 速度（km/h）116:40:02-16:40:23 2141.14285714216:40:02-16:40:32 3028.832.566153316:40:10-16:40:40 3028.885416:40:12-16:40:40 2830.85714286事故前516:40:16-16:40:42 2633.23076923616:41:01-16:41:18 1750.823529
