全国大学生数学建模竞赛A题论文.docx
- 文档编号:9549131
- 上传时间:2023-02-05
- 格式:DOCX
- 页数:34
- 大小:342.84KB
全国大学生数学建模竞赛A题论文.docx
《全国大学生数学建模竞赛A题论文.docx》由会员分享,可在线阅读,更多相关《全国大学生数学建模竞赛A题论文.docx(34页珍藏版)》请在冰豆网上搜索。
全国大学生数学建模竞赛A题论文
2013高教社杯全国大学生数学建模竞赛
承诺书
我们仔细阅读了《全国大学生数学建模竞赛章程》和《全国大学生数学建模竞赛参赛规则》(以下简称为“竞赛章程和参赛规则”,可从全国大学生数学建模竞赛网站下载)。
我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。
我们知道,抄袭别人的成果是违反竞赛章程和参赛规则的,如果引用别人的成果或其他公开的资料(包括网上查到的资料),必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。
我们郑重承诺,严格遵守竞赛章程和参赛规则,以保证竞赛的公正、公平性。
如有违反竞赛章程和参赛规则的行为,我们将受到严肃处理。
我们授权全国大学生数学建模竞赛组委会,可将我们的论文以任何形式进行公开展示(包括进行网上公示,在书籍、期刊和其他媒体进行正式或非正式发表等)。
我们参赛选择的题号是(从A/B/C/D中选择一项填写):
A
我们的参赛报名号为(如果赛区设置报名号的话):
20009017
所属学校(请填写完整的全名):
湖南科技大学
参赛队员(打印并签名):
1.罗晓成
2.周燕霞
3.王志丽
指导教师或指导教师组负责人(打印并签名):
陈先伟
(论文纸质版与电子版中的以上信息必须一致,只是电子版中无需签名。
以上内容请仔细核对,提交后将不再允许做任何修改。
如填写错误,论文可能被取消评奖资格。
)
日期:
2013年9月16日
赛区评阅编号(由赛区组委会评阅前进行编号):
2013高教社杯全国大学生数学建模竞赛
编号专用页
赛区评阅编号(由赛区组委会评阅前进行编号):
赛区评阅记录(可供赛区评阅时使用):
评
阅
人
评
分
备
注
全国统一编号(由赛区组委会送交全国前编号):
全国评阅编号(由全国组委会评阅前进行编号):
车道被占用对城市道路通行能力的影响
【摘要】
随着汽车拥有量的不断增长,交通拥挤已经成为备受关注的世界性问题,几乎所有城市都不同程度地受到这一问题的困扰。
我国也不例外,自20世纪80年代以来,交通拥挤问题越来越严重,特别是大城市交通拥挤现象更为严重。
车辆排队是交通拥挤的一种典型的表现形式。
本文主要讨论交通事故引起车道被占时,在事故发生至撤离期间事故横断面的实际通行能力的变化过程以及车辆排队长度模型。
对于问题1,通过对视频1的观察,统计了交通事故发生至撤离期间,事故所处的横断面每30s内通过该横断面的车流量,根据所统计的数据应用三角函数进行拟合,建立了实际通行能力(TC)模型。
对拟合得到的实际通行能力数值与实际记录所得到的通行能力数值进行比较,发现存在误差,再对实际通行能力模型进行修改,得到了修正后的实际通行能力模型
(1)。
最后通过检验发现修正后模型拟合效果明显提高,与视频1中实际记录到数据基本相吻合。
所得到的模型
(1)及其模拟图形描述实际通行能力如何随时间的变化而变化。
对于问题2,对视频2进行观察,统计了每60s内通过事故所处横断面的车流量。
对此,我们采用与问题1相同的方法建立了实际通行能力模型。
同样进行了误差分析与修正,得到了比较理想的实际通行能力模型
(2)。
且对模型
(1)和模型
(2)进行比较分析,发现事故所占车道的不同对事故横断面的实际通行能力的影响也不同,二、三车道被占时事故横断面的实际通行能力比一、二车道被占时事故横断面的实际通行能力要小。
对于问题3,估算了视频1中交通事故所影响的路段车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量之间的表达式。
我们考虑到用传统的流量守恒方法求得事故路段排队长度所得到的结果与实际情况相差比较大,在此采用车流波动理论来建立车流波动模型,通过视频1中统计得到的数据,对车流波动模型进行验证,发现车辆排队长度模型具有较高精度,因此该模型具有较高的可靠信。
对于问题4,我们根据问题3中所建立的模型,将视频1中统计到的数据代入车流波动模型中,进行实际运算,得到了车辆从事故所处横断面排到上游路口所需要的时间。
关键词实际通行能力模型拟合分析波动模型排队长度
1、问题重述
车道被占用是指因交通事故、路边停车、占道施工等因素,导致车道或道路横断面通行能力在单位时间内降低的现象。
由于城市道路具有交通流密度大、联系性强等特点,一条车道被占用,也可能降低路段所有车道的通行能力,即使时间短,也可能引起车道排队,出现交通堵塞。
如处理不当,甚至出现区域性拥堵。
车道被占用的情况种类繁多、复杂、正确估算车道被占用对城市道路通行能力的影响程度,将为交通管理部门正确引导车辆行驶、审批占道施工、设计道路渠化方案、设置路边停车位和设置非港湾式公交车站等提供理论依据。
视频1(附件1)和视频2(附件2)中的两个交通事故处于同一路段的同一横截面。
且完全占用两条车道。
请研究一下问题:
1.根据视频1(附件1),描述视频中交通事故发生至撤离期间,事故所处横断面实际交通能力的变化过程。
2.根据问题1所得结论,结合视频2(附件2),分析说明同一横断面交通事故所占车道不同对该横截面实际通行能力影响的差异。
3.构建数学模型,分析视频1(附件1)中交通事故所影响的路段车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量间的关系。
4.假如视频1(附件1)中的交通事故所处横断面距离上游路口变为140米,路段下游方向需求不变,路段上游车流量为1500pcu/h,事故发生时车辆初始排队长度为零,且事故持续不撤离。
请估算,从事故发生开始,经过多长时间,车辆排队长度将达到上游路口。
2、问题分析
2.1对问题1的分析
问题1要求我们描述视频1中交通事故发生在二、三车道时,交通事故发生至撤离期间,事故所处横断面实际通行能力的变化过程。
主要是统计视频1中的车流量,通过这些数据分析事故所处横断面通行能力的变化。
利用软件拟合出实际通行能力的表达式,再根据拟合曲线与实际测得的值进行误差比较与误差分析,讨论模型是否须要修正,若须要修改,根据误差曲线再一次进行拟合得到误差曲线的表达式,然后对实际通行能力模型进行修正,最后对修正后的模型进行检验,且分析模型的合理性。
2.2对问题2的分析
视频2中发生事故的车辆是在车道一和车道二上,统计视频2中事故所处横断面的车流量,通过这些数据值建立实际通行能力模型,利用软件能够拟合得到实际通行能力表达式。
再结合问题1所得到的结论和附件3给出不同车道上所占总的车流量的所占的百分比,对交通事故发生在不同车道上该事故所处横断面的实际通行能力进行分析,从而可以知道车道不同对事故横断面实际通行能力之间的差异。
2.3对问题3的分析
从问题1所建实际通行能力模型来看,图中出现较多波峰与波谷,是因为上游信号灯和车道被占所引起的。
为了估算交通事故所影响的路段车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量之间的关系,我们采用车流波动理论来分析发生交通事故后,事故路段排队的形成,从而估算出事故路段车辆排队长度与事故横断面实际通行能力、事故持续时间、路段上游车流量之间的关系。
对此建立车流波动模型,最后根据实际统计到的数据对模型进行验证。
2.4对问题4的分析
由于问题4中给出了事故发生点到事故横断面上游的距离和交通事故上游的车流量的值,根据问题3中建立的车流波动模型和视频1中统计到的数据,依据问题3中所求的表达式进行合理运算。
3、模型假设
1、不考虑观察视频中统计到的数据人为的误差。
2、视频中占用的车道完全考虑为占用2条车道。
3、事故发生到事故撤离期间事故所处横断面的车流量相当于道路的通行能力。
4、假设事故横断面上游红绿灯处于正常工作状态。
5、假设事故发生后,未被占用的车道正常运行。
6、假设从小区出来的车对模型不作影响。
4、符号说明
N修正后的通行能力
修正前的通行能力
实际通行能力
D标准车当量
T事故持续时间
集散波的速度(i=1,2)
Q交通量
L事故路段车辆的排队长度
5、模型建立与求解
5.1问题1模型的建立
5.1.1模型的建立
通行能力指的是在一定的道路和交通条件下,道路上某一路段单位时间内通过某一断面的最大车辆数。
根据参考文献[1]的结论可知,事故发生至事故撤离这段时间内,事故点和事故点下游的各断面短时交通流量相当于事故点断面通行能力。
因此,在事故发生至事故撤离这段时间内,事故点和事故点下游的各断面短时的通行能力可以用交通流量表示。
通过对视频1的观察与统计,可以计算出某段时间内通过事故横断面车流量的值,并对问题1建立实际通行能力(TC)模型,应用模型描述实际通行能力从事故发开始至事故撤离这段时间内的变化过程。
考虑到四轮及以上机动车、电瓶车的交通流量,我们把它换算成标准车的当量数。
调查车型分类及车辆折算系数表
车型
折算
系数
荷载及功率
备注
机
动
车
汽
车
小客车
1.0
额定座位≤19座
大客车
1.5
额定座位>19座
小型货车
1.0
载质量≤2吨
中型货车
1.5
2吨<载质量≤7吨
包括吊车
大型货车
2.0
7吨<载质量≤14吨
特大型货车
3.0
载质量>14吨
拖挂车
3.0
包括半挂车、平板拖车
集装箱车
3.0
摩托车
1.0
包括轻骑、载货摩托车及载货(客)机动三轮车等
拖拉机
4.0
非机动车
人畜
力车
畜力车
4.0
人力车
1.0
包括人力三轮车、手推车
自行车
0.2
包括助动车
根据车辆折算系数表可知,电瓶车折算系数为1,小型车折算系数为1,大型车折算系数为2.对其进行求总和,就可以得到标准车当量数(D)。
标准车当量数的计算表达式为:
D=m1*1+m2*1+m3*2
其中,m1表示电瓶车的数量,
m2表示小型车的数量
m3表示大型车的数量
在视频1中,我们观察到发生交通事故后,事故横断面处车辆连续不断的驶来,把通过事故横断面处的车辆数记为这段时间内的车流量。
我们统计在事故点发生交通事故后,事故所处横断面每隔30s记录一次这一时段内通过的车辆数,可以得到30s内通过事故点横断面的车流量。
统计得到的数据(见附录1)。
可以得到公式:
N=
,
可由记录到的数据值进行三角函数拟合分析,得到以时间为变量实际通行能力(
)的拟合曲线,见图1:
图1.
根据实际通行能力散点分布图进行拟合,利用matlab软件可以得到实际通行能力(
)的表达式:
=a1*sin(b1*x+c1)+a1*sin(b1*x+c1)+a3*sin(b3*x+c3)+
a4*sin(b4*x+c4)+a2*sin(b2*x+c2)+a3*sin(b3*x+c3)+a4*sin(b4*x+c4)+a5*sin(b5*x+c5)+a6*sin(b6*x+c6)+a7*sin(b7*x+c7)+a8*sin(b8*x+c8)
参数取值范围
a1=62.74(-3.284e+008,3.284e+008)
b1=0.01281(-1.504e+004,1.504e+004)
c1=-1.293(-3.238e+006,3.238e+006)
a2=115.2(-1.987e+009,1.987e+009)
b2=0.02052(-5.272e+004,5.272e+004)
c2=0.1535(-1.181e+007,1.181e+007)
a3=81.24(-5.821e+008,5.821e+008)
b3=0.07535(-1337,1337)
c3=-0.7464(-1.602e+005,1.602e+005)
a4=74.47(-2.316e+009,2.316e+009)
b4=0.02354(-3.399e+004,3.399e+004)
c4=2.613(-7.751e+006,7.751e+006)
a5=3.317(-451.8,458.4)
b5=0.1676(-1.346,1.681)
c5=3.966(-372.6,380.5)
a6=3.155(-1115,1121)
b6=0.1418(-5.523,5.806)
c6=1.649(-1807,1810)
a7=77.44(-5.821e+008,5.821e+008)
b7=0.0757(-1227,1228)
c7=2.358(-1.12e+005,1.12e+005)
a8=2.352(-110.2,114.9)
b8=0.1209(-22.45,22.69)
c8=-1.957(-5452,5448)
Goodness of fit:
SSE:
45.1
R-square:
0.9052
Adjusted R-square:
-0.1855
RMSE:
4.749
通过matlab软件,可以知道这些数据值的标准差为4.749,标准差相对比较大,与我们估计的稳定程度不相符。
所以,根据拟合函数曲线上的值与对应实际测到的值进行比较,得到相对应的差值。
通过计算可以得到实际测得的值与曲线值相对应的差值为(1.7937,1.2742,1.1762,1.1085,0.2174,-2.1462,1.1248,-1.3291,1.5042,-0.6304,1.3237,-2.4043,-0.559,-0.9890,2.5147,-1.0371,1.0761,-0.1545,-1.6532,2.1854,-0.9791,0.5404,-0.4830,-0.2647,0.0268,1.0281),再根据这些差值点分布图进行拟合,利用matlab软件再次进行三角函数拟合,得到以时间为变量的通行能力(
)的误差拟合曲线,如图2所示:
图2.
通行能力(
)的误差关系式:
=a1*sin(b1*x+c1)+a2*sin(b2*x+c2)+a3*sin(b3*x+c3)
Coefficients(with95%confidencebounds):
参数取值范围
a1=0.8019(-0.0513,1.655)
b1=0.2077(0.1941,0.2213)
c1=-1.23(-4.371,1.912)
a2=0.708(-0.07101,1.487)
b2=0.05539(0.04433,0.06644)
c2=5.065(2.524,7.606)
a3=0.4332(-0.3942,1.261)
b3=0.1874(0.1571,0.2178)
c3=-2.772(-9.546,4.001)
Goodnessoffit:
SSE:
25.48
R-square:
0.3524
AdjustedR-square:
0.04765
RMSE:
1.224
此时的标准差为1.224,与事故点横断面的实际通行能力(
)的标准差相比较小,根据事故点横断面通行能力(
)的误差曲线和事故点横断面实际通行能力(
)的曲线进行修正,可以得到修正后的通行能力表达式:
N=
+
(1)
然后利用matlab进行模拟以上表达式,得到以时间为变量事故点横断面通行能力(N)的曲线,如图3所示:
图3.
5.1.2模型的求解分析
通过修正后事故点横断面的通行能力曲线波形曲线,并对其进行分析:
1.在事故刚发生时,事故所处横断面的通行能力基本没有受到影响,实际测到的值大致和拟合曲线上的值相吻合。
2.图形中出现的波峰与波谷主要是考虑到上游路口的信号灯,我们选取的时间间隔大致为30s,而交通灯信号周期为60s,当碰到信号绿灯时,事故点上游有大量车行驶至事故点,故事故点横断面的通行能力就相对增大。
反之,碰到信号红灯时,事故点上游行驶的车辆就少,事故点横断面的通行能力则减小。
3.根据波形曲线可以了解到,在视频1中修正后的事故点横断面的通行能力大致在20量/min上下波动,表明20量/min是道路在事故发生后正常的通行能力。
综上可得,我们得到修正后以时间为变量事故点横断面的通行能力的拟合曲线与实际记录到的事故点横断面的通行能力的值大致相符。
5.2问题2模型的建立
5.2.1模型的建立
对于问题2,我们利用问题1建立的模型方法,建立相同的实际通行能力(TC)模型。
对视频2进行统计过程中,取时间间隔为60s,记录通过事故点横断面的车辆数,可以得到每60s内通过事故点横断面的车流量。
统计得到数据(见附录2),根据实际通行能力散点分布图,应用matlab软件中的三角函数进行拟合,得到以时间为变量的实际通行能力(
)的拟合曲线,如图4所示:
图4.
视频2中实际通行能力(
)表达式:
=a1*sin(b1*x+c1)+a2*sin(b2*x+c2)+a3*sin(b3*x+c3)+
a4*sin(b4*x+c4)+a5*sin(b5*x+c5)+a6*sin(b6*x+c6)
Coefficients (with 95% confidence bounds):
参数取值范围
a1=33.67(-5694,5761)
b1=0.002258(-1.818,1.822)
c1=1.062(-302.8,304.9)
a2=2.466(-0.6079,5.54)
b2=0.382(-0.1521,0.916)
c2=-0.2386(-6.329,5.852)
a3=3.012(-0.1502,6.174)
b3=0.5844(0.3388,0.8301)
c3=1.586(-1.599,4.77)
a4=4.071(1.698,6.443)
b4=2.75(2.675,2.824)
c4=1.387(0.1944,2.581)
a5=1.66(-0.5355,3.855)
b5=2.451(2.236,2.667)
c5=5.083(1.77,8.397)
a6=1.896(-0.2983,4.09)
b6=1.464(1.308,1.619)
c6=2.797(0.2946,5.3)
Goodness of fit:
SSE:
147.9
R-square:
0.7725
Adjusted R-square:
0.421
RMSE:
3.667
同样,实际通行能力(
)标准差的值较大,为3.667,而实际测得到的值与曲线上的值存在一定的差值,根据实际记录到的值与曲线所对应点的值相差为(1.2648,-2.2770,-0.2859,1.3908,2.5745,-1.4306,0.0154,-3.4637,-0.6372,1.3576,7.0638,-3.4389,-3.1297,-0.7610,1.2877,-2.2006,4.4557,-0.0261,-1.4477,-0.7357,0.7348,-1.3340,0.4468,-0.5912,2.1186,0.9216,-2.3230,-0.5119,1.1056),再次用matlab软件三角函数拟合,得到以时间为变量通行能力(
)的误差拟合曲线,如图5所示:
图5.
视频2中事故横断面通行能力(
)的误差的表达式:
=a1*sin(b1*x+c1)+a2*sin(b2*x+c2)+a3*sin(b3*x+c3)+a4*sin(b4*x+c4)
+a5*sin(b5*x+c5)+a6*sin(b6*x+c6)+a7*sin(b7*x+c7)+a8*sin(b8*x+c8)
Coefficients (with 95% confidence bounds):
参数取值范围
a1=-4.264(-1.609e+004,1.608e+004)
b1=1.989(-23.74,27.71)
c1=-5.636(-938.2,927)
a2=1.257(0.8349,1.679)
b2=1.264(1.217,1.311)
c2=1.095(0.02663,2.164)
a3=0.9667(-0.1709,2.104)
b3=0.8423(0.6974,0.9873)
c3=-0.8194(-2.701,1.062)
a4=1.346(0.7019,1.991)
b4=1.06(0.8438,1.277)
c4=2.664(-0.1994,5.528)
a5=0.7074(-6.215,7.63)
b5=2.263(1.35,3.177)
c5=0.4878(-13.85,14.83)
a6=4.19(-1.599e+004,1.6e+004)
b6=2.001(-18.64,22.65)
c6=0.1446(-1011,1011)
a7=1.3(-7.042,9.642)
b7=1.711(0.9998,2.421)
c7=-3.869(-14.76,7.022)
a8=1.448(1.256,1.64)
b8=3.03(3.012,3.048)
c8=-0.2043(-0.6098,0.2012)
Goodness of fit:
SSE:
0.3036
R-square:
0.9979
Adjusted R-square:
0.9885
RMSE:
0.2464
经过这次拟合后,得到的视频2中事故横断面通行能力(
)误差的标准差的值相对比较小,且标准差为0.2464,符合我们的猜测。
根据分析,可以得到:
=
+
(2)
对其数据进行修正处理,得到以时间为变量修正后的通行能力(
)的拟合曲线,如图6所示:
图6.
5.2.2模型的分析
对曲线进行分析可得:
1.在事故刚发生时,通行能力有下降趋势,可能是受到上游信号灯的影响,行驶在事故点上游的车辆比较少,通过事故点横断面的实际通行能力就小。
2.视频2中的通行能力的值差不多在30量/min上下波动,表明这是车道一和车道二被占后,事故处横断面正常
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 全国大学生 数学 建模 竞赛 论文
![提示](https://static.bdocx.com/images/bang_tan.gif)