基于直方图优选的车速分布模型研究_曹瑾鑫
- 格式:pdf
- 大小:396.03 KB
- 文档页数:7
2012年1月内蒙古大学学报(自然科学版)Jan.2012第43卷第1期Journal of Inner Mongolia University(Natural Science Edition)Vol.43No.1 文章编号:1000-1638(2012)01-0089-07基于直方图优选的车速分布模型研究*曹瑾鑫1,肖崇紫1,2,武 钧1(1.内蒙古大学交通学院,呼和浩特010070;2.北京交通大学交通运输学院,北京100044)摘要:通过对国内外已有的研究分析,发现车速的分布规律变化不一,很难用单一的某一种分布模型进行拟合.因此,通过对呼和浩特市昭乌达路的速度数据采集,利用直方图优选的方法对速度分布直方图进行改进,并利用正态分布、威布尔分布、伽玛分布和Logistic分布模型对改进后的直方图进行了拟合分析.通过对拟合分析得出的相关系数的比较,给出了这四种分布在拟合速度分布时的适用范围.关键词:速度分布;直方图;相关系数;拟合函数中图分类号:U491.1 文献标志码:A 车辆的行车速度是一个随机变量,随着驾驶员在行驶过程中所驾车辆、驾驶习惯及行车路段的不同而不同.在不同流量下,车速是在一定范围内随机出现的,且其分布符合一定的数学规律.国内外的学者大都认为车速应符合正态分布、威布尔分布、伽玛分布或Logistic分布,但并没有人给出在拟合速度分布时各类分布模型的应用条件.因此,有必要研究这几种分布的拟合特性.1 对已有研究的分析当前,不少学者已对车速的分布做出了相应的研究,并得出了各自的结论.孙逢春等人对北京市的公共汽车车速做了相应的研究,并提出了用对数正态分布去拟合车速分布〔1〕;周宏敏〔2〕,李江〔3〕等人认为高速公路运行车速服从正态分布;阎莹等人通过研究认为车速的分布应视情况而定〔4〕;SunCarlos等人绘出了不同时段的速度分布图,并对速度的分布特性进行了一系列的总结〔5〕.综合已有的研究成果,可以看出,由于车速受到多方面因素的影响,其分布特性会因时间、地点、交通参与者、道路环境以及车况的不同而有所差异.车速的分布是比较灵活多变的.在实际中,试图寻找一种通用的分布模型是不太实际的,且大多数的情况下,车速的分布都存在着一定的偏度.既存在左偏,又存在右偏,这就要求研究者针对不同的车速分布类型找出适用于各自特性的分布函数.2 数据采集针对本研究,项目组选取呼和浩特市主干道昭乌达路作为研究对象,对该路段进行长达6小时的视频拍摄,并利用VirtualDub软件对拍摄的视频进行标准化处理,采集路段见图1.由于市区干道的交通流是以间断流的形式存在的,故以一分钟为时间步长来截取每一个信号周期内的车流峰值视频,而后对视频进行合并和格式转换.利用先进的视频处理设备Traficon VIP/D来提取相应的交通流数据,共得到五组速度分布数据(如图2到图6所示).*收稿日期:2011-06-30;修回日期:2011-012-08基金项目:国家自然科学基金项目(71061010);内蒙古自治区自然科学基金项目(2010MS0705);内蒙古自治区高等学校科学研究计划项目(NJ10009);内蒙古大学高层次人才引进启动基金项目(210221)作者简介:曹瑾鑫(1982-),男,内蒙古呼和浩特市人,副教授,博士.E-mail:caojinxin@gmail.com.图1 数据采集路段Fig.1 Road of data acquisition图2 流量小于100veh/h(第1组)Fig.2 Volume is less than 100veh/h(the first group of data)图3 流量在1440-1800veh/h间(第2组)Fig.3 Volume is between 1440and 1800veh/h(the second group of data)图4 流量在1980-2340veh/h间(第3组)Fig.4 Volume is between 1980and 2340veh/h(the third group of data)图5 流量在2520-2880veh/h间(第4组)Fig.5 The volume is between 2520and 2880veh/h(the fourth group of data)图6 流量大于3060veh/h(第5组)Fig.6 The volume is greater than 3060veh/h(the fifth group of data)3 几种拟合模型的特性国外对车速数据进行正态分布、威布尔分布、伽玛分布以及Logistic分布的研究表明,在乡村公路和高速公路上,车速一般呈正态分布;在城市道路或高速公路匝道入口处,车速比较集中,一般呈偏态分布〔6〕.我国市区干道速度分布具有离散性较大的特点,其特性还有待进一步研究.09内蒙古大学学报(自然科学版)2012年正态分布又被称为高斯分布〔7〕,是1809年德国著名数学家在研究天文学的观测误差时推导出来的.若随机变量x的概率密度函数为:f(x)=1/(2槡πδ)exp{-(x-μ)2/(2δ2)},-∞!x!+∞(1) 则称x服从均值为μ,方差为的δ2正态分布,记作N(μ,δ2),式中μ为均值,δ为标准差.威布尔分布〔7〕是瑞典科学家威布尔(Weibull)于1937年提出的一种分布函数,是为失效分析而开发的.其概率密度函数为:f(x)=(β/λ)(x/λ)β-1e-(x/λ)β,x>0(2) 式中λ称为刻度参数,决定该分布的分散度;β称为形状参数,决定该分布的形态.伽玛分布和威布尔分布一样,其基本形式为一种双参数分布〔7〕,其概率密度函数为:f(x)=λγxγ-1e-λx/(∫∞0tγ-1e-tdt),x>0(3) 式中λ称为刻度参数;γ称为形状参数.Logistic函数是由比利时数学家P.F.Veerhulist于1844年创建〔7〕,最早来自于生长曲线的需要,后来用于经济和人口统计学中,从十九世纪末以来已逐渐普及,近代又用于回归和判别之中.其概率密度函数为:f(x)=exp-(x-α)/{}β/(β1+exp-(x-α)/{}[]β2);-∞<x,α<+∞,0<β(4) 式中α称为位置参数;β称为尺度参数.通过对以上四种拟合模型图像的分析发现,正态分布函数和Logistic函数具有强对称性,而伽玛函数具有左偏的特性.威布尔函数最为灵活,根据参数的不同,它既可以变换成具有左偏特性的函数,又可以变换成右偏的函数,并且它还可以是对称的.4 直方图优选对于一定量的数据,直方图能很直观的描述出速度的分布特点.数据量越大,直方图越能够精确的表述出这些数据的实际分布特性.但是在实际中往往不可能取到足够多的数据,甚至数据总量常常会小于200.在这种情况下,如果所作的某一直方图的组数和组距都为一个定值,而第一组的起点值不一样时,得到的直方图会有很大的差异.以图5中的数据为例,取直方图的组数为8,组距为1,如果第一组的起点值分别为7.2、7.6和7.8,那么得到的三组频率直方图如图7(a)到图7(c)所示.(a)起点为7.2(b)起点为7.6(c)起点为7.8图7 不同起点时的直方图Fig.7 Histogram with the different starting points由图7可知,在数据总量较小时,不同的取值起点对直方图的分布形态会产生极大的影响.实际中,作直方图时的取值起点往往是随机选取的,倾向于取某一个整数.为了使直方图更加平滑,有必要去确定一种评选形状最优直方图取值起点的方法.假设直方图的组数为n,组距为d,且完全服从某一种分布g(x),其关系如图8所示:当n趋向于无穷大时,那么d也就无限的趋近于零,对于其中的任意m(1<m<n),有以下关系成立[g(m-1)+g(m+1)2-g(m)]2→0(5)19第1期曹瑾鑫等 基于直方图优选的车速分布模型研究 故有∑n-1m=2[g(m-1)+g(m+1)2-g(m)]2→0(6) 依据这一关系,可以建立一个直方图优劣的判断指标#.令χ=∑n-1m=2[g(m-1)+g(m+1)2-g(m)]2(7) 当#越小时,可以认为直方图的平滑性越好.因此,最优直方图的选取算法如下:Z=min(χt)=min{∑n-1m=2[gt(m-1)+gt(m+1)2-gt(m)]2}(8) 上式中下标t表示第一组数据的取值起点,Z表示最优目标值.直方图的形态由组距、组数以及第一组的取值起点三个参数决定.在组距和组数确定的条件下,可以通过Z值来确定最佳的取值起点t,进而得出形态最好的直方图.运用这一算法对图2的数据进行处理,得到不同取值起点t和对应的目标函数Z的关系值,见表1.图8 直方图和分布曲线之间的关系Fig.8 Relationship between histogram and distribuing line表1 t和Z的数值关系Table 1 Relationship between t and Zt 4 4.1 4.2 4.3 4.4 4.5 4.6 4.7 4.8 4.9Z 494 1077 1235 561 485 554 1485 1521 535 361 由于事先并不知道直方图是服从那一种分布,故用借助曲线拟合软件1thopt,用威布尔,Logis-tic、正态和对数正态分布等多种曲线拟合出不同取值起点下所对应的直方图,选取出最优的拟合系数R2,得到Z和R2的相互关系如表2所列.表2 Z和R2的数值关系Table 2 The relationship between Zand R2Z 494 1077 1235 561 485 554 1485 1521 535 361R2 0.78 0.69 0.71 0.81 0.84 0.78 0.67 0.72 0.79 0.85 利用一次线性方程对这些数据进行拟合,得出的拟合方程为y=-0.0001x+0.8684,复相关系数为0.8055.可见,目标值Z和直方图的拟合相关系数之间有着很好的线性关系,故可以认为上述方法适用于选取直方图.5 速度拟合为了对图2到图6中的数据进行拟合,需先对其相应的直方图进行优化.利用前述优化算法,取29内蒙古大学学报(自然科学版)2012年组距d=1,取值起点根据各自的数值范围取一合适的整数,然后求出目标值Z.之后,再将取值起点的值增加0.1,组距保持不变,然后再求对应的目标值.这样循环计算,直到找出最小的目标值Z,并画出相应的直方图.利用MATLAB编程可实现这一运算过程,并绘制出最佳的直方图.其中图2至图6的数据所对应的最优直方图如图10至图14所示.在找出最优直方图之后,再利用前述四种概率密度曲线对优化后的直方图进行拟合,即可得到的拟合参数,如表3所列.表3 拟合参数Table 3 Matching parameters6 拟合模型分析表3仅给出了拟合参数,并没有说明哪种拟合方法是最优的.为了找出最优的拟合曲线,需要求出曲线的拟合度.复相关系数R2是一种常见的拟合度的判断方法,常用于判断线性曲线的拟合度,但将这种线性相关系数用于非线性拟合曲线是不合理的〔8-9〕.为此,可引入一种通用的相关系数Rc的计算方法〔10〕.其公式如下:Rc=∑ni=1(yi-珔y)(zi-珔z)/(∑ni=1(yi-珔y)2∑ni=1(zi-珔z)[]212)(9) 式中yi是实际值,zi为预测值,珔y=∑ni=1yi/n,珔z=∑ni=1zi/n.运用式(9)对所有的拟合模型进行检验,可得出相关系数Rc,其数值关系可见表4和图9. 从图9中可以发现,第四组的四种拟合模型拟合效果最差,主要原因是其数据量不大,分布不规则所引起的.第一组,第二组和第五组数据的拟合效果均比较好,而第三和第四组的数据拟合效果稍差.表4 相关系数值Table 4 Related coefficients组别1 2 3 4 5均值分布正态分布相关系数0.929 0.945 0.768 0.669 0.880 0.838分布威布尔分布相关系数0.923 0.987 0.866 0.715 0.938 0.886分布伽玛分布相关系数0.920 0.939 0.915 0.671 0.985 0.886分布Logistic分布相关系数0.922 0.948 0.812 0.654 0.972 0.862 通过对图像的观察和拟合相关系数的分析发现,图10和图14具有比较强的对称性,而正态分布在拟合它们时并不比其他的模型占据优势,甚至在图14的拟合中正态分布模型的拟合效果要远差于其他模型的拟合效果.对于峰值比较突出,离散程度较小,偏度不明显的图12,正态函数拟合的相关系数仅为0.768,远小于伽玛分布的0.915和另外两种拟合的相关系数.图11和图13有明显的左偏,在拟合它们时,威布尔函数所得到的相关系数要比其他拟合都要高.但是,第四组数据的整体拟合效果要明显差于另外几组数据的拟合效果.这主要是由于速度数据左偏比较严重,且在速度为14m/s附近时出现了较大的跳跃,因而使得四种模型的拟合效果较差.从相关系数的均值看,正态分布的相关系数均值最小,威布尔分布和伽玛分布的相关系数值相等且高于其他均值.由此,可以得出以下结论:1)正态分布在拟合速度分布时并没有明显的优势,且当速度分布不对称时所得到的拟合结果不39第1期曹瑾鑫等 基于直方图优选的车速分布模型研究理想.2)在拟合离散性小、峰值突出的速度分布时,伽玛函数和Logistic函数具有明显的优势.3)当速度分布存在左偏时,威布尔函数能很好的实现拟合.对表4中的四组相关系数的均值作比较,可见威布尔分布和伽玛分布的值相等且大于另外两组的均值,但由于伽玛函数比威布尔函数形式复杂,且仅存在右偏,该模型不能很好的拟合具有左偏特性的速度数据.由此,可以认为威布尔函数最适宜用来描述速度的分布特性.图9 相关系数比较图Fig.9 Comparing of the related coefficients图10 第一组数据拟合图Fig.10 Fitting of the first group of data图11 第二组数据拟合图Fig.11 Fitting of the second group of data图12 第三组数据拟合图Fig.12 Fitting of the third group of data图13 第四组数据拟合图Fig.13 Fitting of the fourth group of data图14 第五组数据拟合图Fig.14 Fitting of the fifth group of data49内蒙古大学学报(自然科学版)2012年7 结语本文通过对国内外已有速度分布研究结果的分析,发现速度分布的形式多种多样,很难单一的用某一种模型去拟合.因此,有必要找出各种常用拟合模型的适用特性.根据数据统计分析,当原始的统计数据总量小于200时,所绘制的直方图会因为计数起点的不同而存在较大的差异.因此,本文提出了一种直方图的优选算法,用以选取同一组数据下外观最为平滑的直方图.研究表明,随着直方图平滑程度的增大,其拟合效果也会随之变好.将此算法用于呼和浩特市主干路的实测数据优化处理,并用四种最常用的速度分布模型分析这些数据的拟合效果,从而得出拟合参数.最终通过相关系数的比较,找到了这几种分布的适用特性.参考文献:[1] 孙逢春,王震坡,王军.北京市公共汽车平均车速统计分析[J].汽车工程,2003,25(3):219-222.[2] 周宏敏,马玉成,王君,等.高速公路断面运行车速分布的研究[J].华东交通大学学报,2008,25(5):32-35.[3] 李江,孙立军.一种基于自由流车速的交通数据检验方法[J].交通运输系统工程与信息,2008,8(4):30-35.[4] 阎莹,王晓飞,张宇辉,等.高速公路断面运行车速分布特征研究[J].中国安全科学学报,2008,18(7):171-176.[5] Sun C,Ritchie,Stephen G.Individual Vehicle Speed Estimation Using Single Loop Inductive Waveforms[R].California:California Partners for Advanced Transit and Highways(PATH),1999.[6] Fitzpatrick K,Elefteriadou L,Douglas W.Speed Prediction for Two Lane Rural Highways[R].Washington,CD.:Federal Highway Administration,2000.[7] 茆诗松,濮晓龙,刘忠.寿命数据中的统计模型与方法[M].北京:中国统计出版社,1998.[8] 张世强.曲线回归的拟合优度指标的探讨[J].中国卫生统计,2002,19(1):9-11.[9] 张世强.非线性生物模型回归参数计算的一个新方法及应用[J].重庆医科大学学报,2003,28(6):754-757.[10] 张世强,吕杰能,蒋峥,等.关于相关系数的探讨[J].数学的实践与认识,2009,39(19):102-107.(责任编委 柴金义)An Analysis of Curve Fitting Functions for vehickSpeed Distribution Based on Histogram Improving MethodCAO Jin-xin1,XIAO Chong-zi 1,2,WU Jun1(1.Faculty of Transportation,Inner Mongolia University,Hohhot 010070,China;2.School of Traffic and Transportation,Beijing Jiaotong University,Beijing100044,China) Abstract:According to the findings of domestic and foreign studies,it is found that the range ofvehicle speed distribution is enormous,and it’s very difficult to use a single model to describe thespeed distribution.So,the related data was collected in the Zhaowuda Road of Hohhot by the projectteam,then a histogram improving method is proposed to deal with bar charts of speed distributingand four kinds of curve fitting functions(Normal,Weibull,Gamma,Logistic)are adopted to describethe distribution.After the comparison of the correlation coefficient,the applicability of each functionis presented. Key words:speed distribution;histogram;correlation coefficient;curve fitting function59第1期曹瑾鑫等 基于直方图优选的车速分布模型研究。