CN110266528A - 基于机器学习的车联网通信的流量预测方法 - Google Patents
基于机器学习的车联网通信的流量预测方法 Download PDFInfo
- Publication number
- CN110266528A CN110266528A CN201910504819.4A CN201910504819A CN110266528A CN 110266528 A CN110266528 A CN 110266528A CN 201910504819 A CN201910504819 A CN 201910504819A CN 110266528 A CN110266528 A CN 110266528A
- Authority
- CN
- China
- Prior art keywords
- data
- traffic
- communication
- vehicle
- index
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004891 communication Methods 0.000 title claims abstract description 110
- 230000006855 networking Effects 0.000 title claims abstract description 26
- 238000010801 machine learning Methods 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000004088 simulation Methods 0.000 claims abstract description 34
- 238000007637 random forest analysis Methods 0.000 claims abstract description 24
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 238000004458 analytical method Methods 0.000 claims abstract description 12
- 238000002156 mixing Methods 0.000 claims abstract description 3
- 230000006870 function Effects 0.000 claims description 38
- 238000012549 training Methods 0.000 claims description 34
- 238000012360 testing method Methods 0.000 claims description 23
- 230000029305 taxis Effects 0.000 claims description 16
- 238000012546 transfer Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 8
- 238000010606 normalization Methods 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 5
- 238000000926 separation method Methods 0.000 claims description 5
- 238000010219 correlation analysis Methods 0.000 claims description 3
- 238000007405 data analysis Methods 0.000 claims description 3
- 238000005538 encapsulation Methods 0.000 claims description 3
- 238000006116 polymerization reaction Methods 0.000 claims description 3
- 241001269238 Data Species 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000002790 cross-validation Methods 0.000 claims description 2
- 235000013399 edible fruits Nutrition 0.000 claims description 2
- 238000013277 forecasting method Methods 0.000 claims description 2
- 238000013139 quantization Methods 0.000 claims description 2
- 230000006978 adaptation Effects 0.000 claims 1
- 230000002708 enhancing effect Effects 0.000 abstract description 2
- 238000003066 decision tree Methods 0.000 description 14
- 238000005516 engineering process Methods 0.000 description 5
- 238000007477 logistic regression Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 238000012706 support-vector machine Methods 0.000 description 3
- 206010039203 Road traffic accident Diseases 0.000 description 2
- 230000002411 adverse Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 206010036590 Premature baby Diseases 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
- H04L41/14—Network analysis or design
- H04L41/145—Network analysis or design involving simulating, designing, planning or modelling of a network
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L41/00—Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
- H04L41/14—Network analysis or design
- H04L41/147—Network analysis or design for predicting network behaviour
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Traffic Control Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于机器学习的车联网通信的流量预测方法,利用交通数据平台发布的流量速度数据库,选取八类数据指标,通过优化参数后的随机森林算法完成全天候的车流量预测;利用openstreetmap导出某地城市道路车辆交通场景,获取交通数据,再配置通信仿真文件,获取通信数据,两种数据混合,分析车流量与通信流量二者关系;利用openstreetmap导出在交通数据平台所选的路段,配置通信仿真文件,获取通信数据,从交通数据平台发布的流量速度数据数据与通信数据中选取九类相关指标,通过Bagging模型进行通信流量预测。本发明的方法泛化性能好,准确度较高,能为后期利用经济高效的数据分发提供可靠的车载通信分析方法,增强车辆用户驾驶的安全性。
Description
技术领域
本发明涉及城市道路车辆交通场景下车辆的流量预测技术领域,具体是一种利用机器学习算法对现有的交通数据进行预测,并结合通信仿真完成车联网通信流量预测的方法。
背景技术
车载自组织网络是依托计算机网络、现代无线通信和云计算等新一代信息技术的革命性发展,它的开发是为了通过经济高效的数据分发提供可靠的车载通信。车辆通信可用于减少交通事故,交通拥堵,行驶时间,燃料消耗等。车载通信允许道路使用者通过交换一些信息,了解可能发生在他们身上的危急和危险情况时的周围环境。因此对车联网通信流量的研究能进一步促进交通智能化,能有效地提升交通的效率、减少甚至避免交通事故的发生。
在预测交通流量和网络流量的非机器学习方法中,大部分都只能针对特定区域的短时交通流进行预测。如两步优化选择法是对时间序列进行预测的统计方法,但它只能基于单条时间序列进行检测和统计,而对多条时间序列曲线则相对困难些。另一种结合波动理论分析和频谱分析的组合方法,是根据频谱分析将流量数据分成三类成分,而不同的流量成分采用相应的模型预测,但它在处理庞大的交通流量和网络通信流量数据时,不仅耗时耗力,而且数据指标的分布特征都没有很好地得到展示。而近年来陆续被使用在交通流量和网络流量的机器学习方法,都能综合地考虑到交通场景中的时空性数据,且应对海量的流量数据时,能高效率地统计交通流量分布规律,预测未来流量情况。
发明内容
本发明的目的在于提供一种基于机器学习的车联网通信的流量预测方法,提高预测模型的泛化性能和全天候流量预测的准确性。。
实现本发明目的的技术解决方案为:一种基于机器学习的车联网通信的流量预测方法,步骤如下:
步骤1、车流量预测:利用交通数据平台发布的流量速度数据,选取8类指标进行数据处理,完成全天候的车流量预测;
步骤2、分析两种流量关系:利用openstreetmap导出某地城市道路车辆交通场景,配置车流探测文件,获取交通数据,再配置车辆通信仿真文件,获取通信数据,两类数据混合,分析车流量与通信流量二者关系;
步骤3、通信流量预测:在交通数据平台选择某些路段,利用openstreetmap导出这些路段的城市道路车辆交通场景后,配置车辆通信仿真文件,获取通信数据,依据已得的两种流量关系,从交通数据平台发布的流量速度数据和获取的仿真车辆通信数据中,选取9类指标进行数据处理,完成车载网络的通信流量预测。
本发明与现有技术相比,其显著优点:1)本发明采用机器学习算法,将海量经验数据加载给学习算法训练出模型,泛化性能好,能直观看出各指标下高低流量分布情况,针对城市道路车辆交通场景下的数据都有很好的预测性能;2)不局限于其他只针对仿真场景下的车流量短时预测的研究,本发明有效地对真实城市车辆交通场景的未来全天候车流量进行了很好的预测;3)本发明更进一步研究了车流量等交通信息与车联网联实时的通信质量之间的关系,为车联网的技术发展甚至智能交通的构建提供较大的便利。
附图说明
图1是本发明据实测数据预测车流量的方案流程图。
图2是本发明对关联度大且多值指标的处理流程图。
图3是本发明据车联网仿真预测通信流量的方案流程图。
图4是本发明预测通信流量预测的方案流程图。
具体实施方式
下面结合说明书附图和实施例对本发明作进一步说明。
本发明基于机器学习的车联网通信的流量预测,包括如下步骤:
步骤1、车流量预测:利用交通数据平台发布的流量速度数据,选取8类指标进行数据处理,如图1所示完成全天候的车流量预测的具体方法为:
获取交通数据平台发布的某些路段一周全天候的数据,利用这些数据,预测这些路段未来某周某天全天候的车流量,即形成7天的训练数据集train.csv,1天的测试数据集test.csv;8类数据集指标包括6类数字型指标和2类类目型指标,数字型指标:车流量Flow、时间Time、速度Speed、交通量Volume、交通指数Index、路段Place,类目型指标:即星期Week、天气Weather;
利用isna函数判断各指标是否有缺失值,并用len函数对训练集数据train.csv进行统计,包括:各数字型指标的数量count、均值mean、标准差std、最小值min、最大值max;
为了直观看出各指标下高低车流量分布情况,以车流量的均值mean作为高低流量的分界点,用sns.kdeplot函数画内核密度图查看车流量的分布状况:1)类目型指标week、weather上的高流量与低流量分布与量化,2)Week各天的高流量与低流量分布情况,3)交通指数Index上高流量与低流量分布情况,4)Time全天的高流量与低流量分布情况,5)各位置Place的高流量与低流量分布情况,6)volume、speed上的高流量与低流量分布情况;
利用sns.heatmap画图展示各指标间的关联度,关联度值的绝对值越大,两指标间的相互影响关系越深,用groupby函数得到车流量Flow与其它指标的关联度排序。
关联度排序在最后的两种指标,对最终预测结果的准确性程度贡献小,为了提高运算效率,课将这两种指标合并为新指标;关联度排序在前三的三种数值连续性指标,依次设定指标数值分成8、6、4段取值范围,利用pd.cut函数自动寻找各段交界点,以各范围边界点来划分数值连续性指标,获得类目型指标;再使用get_dummies对类目型指标进行数据因子化,将指标平展开后,取值均为0、1,即当某展开指标的取值为1时,它会确定Flow的一个模糊的取值范围,再结合其它指标的贡献,逐渐收敛,一步步将Flow的取值范围缩小,最后确定预测值。由于特征值的取值只有0和1,这会大大加快机器学习运算效率。对关联度大的指标进行取值范围划分,对处理海量的机器学习算法来说,不仅大大提高了运算效率,还在数据收敛方向上,提供了机器学习算法自动划分的的分界点与聚集中心;
如果连续性数据指标取值范围太大,将对逻辑回归的收敛性造成不利影响,利用preprocessing.StandardScaler()函数,对取值超过100种的多样性数据指标进行数据标准化,将数据指标分布调整成标准正太分布,使得其均值为0,方差为1;
将处理后的训练数据集train.csv,输入scikit-learn工具包下的六种常用经典机器学习算法进行模型训练:随机森林模型Random Forest、逻辑回归模型LogisticRegression、K近邻模型KNN、支持向量机模型Support Vector Machines、朴素贝叶斯模型Naive Bayes 、决策树模型Decision Tree等。由于scikit-learn工具包为模型对象提供predict()接口,经过训练的模型,可以用这个接口来进行预测,同时也提供了score()接口来评价一个模型的好坏,根据置信度得分越高则模型预测性能越好,来进行预测模型的选择和验证。
根据给出的置信度得分排序,初步选择置信度得分排在前两位的随机森林模型和决策树模型,对置信度得分最高的随机森林模型和决策树模型进行优化:利用不同的机器学习算法的可调参数以及它们对机器学习算法的性能以及准确率的影响,通过GridSearchCV()、StratifiedKFold()和StratifiedShuffleSplit()函数,寻找最优参数下的决策树模型Decision_best和最优参数下的随机森林模型RandomForest_best。
将处理后的训练数据集train.csv,输入最优参数下的决策树模型Decision_best和最优参数下的随机森林模型RandomForest_best进行模型训练,并计算置信度得分,用置信度得分最高的优化参数下的随机森林模型RandomForest_bes对测试集test.csv进行预测,并提交车流量预测结果fianl_submission1.csv
为了一目了然获得训练数据集的样本大小是如何影响所选模型的置信度得分,避免欠拟合和过拟合情况的发生,需要利用优化参数后的随机森林模型画学习曲线:横坐标为训练集train.csv所训练样本个数,纵坐标为训练集train.csv拟合的准确性得分和交叉验证集预测的准确性得分,画出学习曲线,验证当前训练样本的数量大小下,随机森林模型不仅能很好地拟合训练数据集,而且对测试数据集的预测准确性很高。
本发明采用最优参数下的随机森林对全天候的车流量进行预测,在处理海量的流量速度数据时,能很好地避免出现过拟合问题,且模型泛化性能高。
目前车联网相关的基础设施搭建还未成熟,相关的车载网络通信数据获取成本高且不能保障其有效性,故本发明的通信流量预测工作是基于仿真平台获取的通信数据,利用交通仿真平台SUMO的城市道路车辆交通场景与网络仿真平台NS-3的车载通信网络的连接与搭建,获得仿真车流量与通信流量,分析车流量与通信流量二者关系,进而可以选取对通信流量的预测准确性贡献较大的指标,进而完成由真实车流量来预测未来车联网网络通信状况。
步骤2、分析两种流量关系:利用openstreetmap导出某地城市道路车辆交通场景,配置车流探测文件,获取交通数据,再配置车辆通信仿真文件,获取通信数据,两类数据混合,分析车流量与通信流量二者关系的具体方法为:
基于开源平台openstreetmap导出某实地城市道路交通路口,作为用来实验的城市道路车辆交通场景进行建模:配置车流文件rou.xml、道路文件net.xml、探测器loopsDetectors文件等,然后用交通仿真软件SUMO进行交通仿真,获取交通数据:交通量nVehContrib(采集周期内完整通过传感器的车辆数量)、车流量Flow(每小时通过传感器的车辆数量)、时间占有率occupancy(车辆通过的累计值与时间周期的比值)、平均速度speed(采集周期内的平均车速)、通过探测器的车辆平均长度length、接触到传感器车辆的数量nVehEntered。
基于ubuntu16.04操作系统,完成通信仿真软件NS-3与交通仿真软件SUMO的连接,使NS-3能够遵循SUMO实验城市道路车辆交通场景模式输出通信数据,方法为:将车流文件.rou.xml和道路文件.net.xml合并转化为.sumo.tr文件;随后利用java -jar命令下的traceExporter.jar,将道路文件.net.xml和.sumo.tr文件合并转化为NS-3需要的.mobility.tcl文件,最后通过ns2-mobility-trace连接NS-3与SUMO。
配置车辆通信仿真文件dsdv-njust1.cc,要完成路由协议dsdv、车速Speed、通信结点个数Nodes、仿真时间Time、输出文本类型等配置,主要添加std::string CSVfileName= "njust-dsdv1.csv"语句和std::string tr_name = "Dsdv_Manet_" + t_nodes + "Nodes_" + sTotalTime + "SimTime";std::cout << "Trace file generated is " <<tr_name << ".tr\n"语句,运行车辆通信仿真文件dsdv-njust1.cc,获取通信数据:路由协议RoutingProtocol、接收到的数据包数目PacketsReceived、数据包接受率ReceiveRate和可供网络链路数据分析的网络抓包.pcap文件。
使用Wireshark的统计功能,选用网络配适滤波器对网络抓包.pcap文件进行分组封装,对通信结点链路和信道的数据包进行数据撷取,获取另一部分通信数据:数据包传送速率PacketsV和平均带宽Band。
对获取到的车流量和通信流量进行相关性分析:使用sns.catplot查看仿真时间内,车流量Flow与数据包传送速率PacketsV、平均带宽Band的相关度,验证交通车流量和车联网通信流量之间的相关性:交通车流量和车联网通信流量存在一定程度的正相关,证明可选取城市道路车辆交通场景与车载通信网络下的数据指标,进行通信流量预测。
步骤3、通信流量预测:在交通数据平台选择某些路段,利用openstreetmap导出这些路段的城市道路车辆交通场景后,配置车辆通信仿真文件,获取通信数据,依据已得的两种流量关系,从交通数据平台发布的流量速度数据和获取的仿真车辆通信数据中,选取9类指标进行数据处理,完成车载网络的通信流量预测的具体方法为:
在交通数据平台选择某些路段,利用openstreetmap导出这些路段的城市道路车辆交通场景后,配置车辆通信仿真文件dsdv-njust2.cc:配置路由协议dsdv、车速Speed、通信结点个数Nodes、仿真时间Time、输出文本类型等;添加std::string CSVfileName = "njust-dsdv2.csv"语句,使得运行通信文件dsdv-njust2.cc后,获得njust-dsdv.csv文件,文件包括一部分通信数据:接收到的数据包数目PacketsReceived、路由协议RoutingProtocol、数据包接受率ReceiveRate;添加std::string tr_name = "Dsdv_Manet_" + t_nodes + "Nodes_" + sTotalTime + "SimTime";std::cout << "Trace file generated is " <<tr_name << ".tr\n"语句,使得运行车辆通信仿真文件dsdv-njust2.cc后,获得网络抓包.pcap文件;
使用Wireshark的统计功能,选用网络配适滤波器对网络抓包.pcap文件进行分组封装,对通信结点链路和信道的数据包进行数据撷取,获得另一部分通信数据:数据包传送速率PacketsV和平均带宽Band;
依据已得的交通车流量和车联网通信流量之间的相关性,从交通数据平台发布的这些路段流量速度数据和仿真车辆通信数据中,选取9类指标:平均带宽Band、数据包传送速率PacketsV、车辆数目Vehicles、交叉口数目Junctions、道路数目Lane、速度Speed、数据包接受率ReceiveRate、接收到的数据包数目PacketsReceived、路由协议RoutingProtocol;
利用所选路段在不同车速Speed和不同通信结点个数Nodes下获取的9类数据指标,预测这些路段未来某周某天的通信流量,即形成训练数据集train.csv和测试数据集test.csv;
对9类指标进行统计:利用isna函数判断各指标是否有缺失值,并用len函数对训练集数据train.csv进行统计,包括:各指标的数量count、均值mean、标准差std、最小值min、最大值max ;以平均带宽Band的均值mean作为高低通信流量的分界点,用sns.kdeplot函数画内核密度图查看通信流量的分布状况:1)速度Speed上的高流量与低流量分布,2)Lane各道路的高流量与低流量分布情况,3)车辆数目Vehicles上高流量与低流量分布情况,4)数据包传送速率PacketsV上的高流量与低流量分布情况,5)交叉口数目Junctions上的高流量与低流量分布情况,6)接收到的数据包数目PacketsReceived、数据包接受率ReceiveRate上的高流量与低流量分布情况;
对9类指标进行处理:利用sns.heatmap画图展示各指标间的关联度,用groupby函数得到平均带宽Band与其它指标的关联度,并进行关联度排序;对关联度排序在最后的三种指标,合并为新指标;对关联度排序在前四的连续性指标,依次设定指标数值分成8、6、4、2段取值范围,利用pd.cut函数自动寻找各段交界点,以各范围边界划分连续性指标,获得类目型指标;再使用get_dummies对类目型指标进行数据因子化,将指标平展开;
利用preprocessing.StandardScaler()函数,对取值超过50种的多样性数据指标进行数据标准化,将数据指标分布调整成标准正太分布,使得其均值为0,方差为1
利用scikit-learn工具包下的经典机器学习算法自助聚合BaggingClassifier模型对处理后的9类指标进行训练,然后对测试集test.csv进行预测,并提交结果通信流量预测结果fianl_submission2.csv。
实施例1
为了验证本发明的有效性,利用交通数据平台发布的流量速度数据先进行车流量预测,具体如下:
采用上海市大数据联合创新实验室(交通领域)平台发布的2018年9月1日至9月7日的上海延安高架12个路段全天候的数据,来预测这些路段在9月8日全天候的车流量,即形成7天的训练数据集共60480组(train.csv),1天的测试数据集共8640组(test.csv)。数据集指标包括8类:车流量 Flow、星期 Week、天气Weather、时间Time、速度Speed、交通量Volume、交通指数Index、路段Place。
利用isna函数判断各指标是否有缺失值,并用len函数对训练集数据train.csv进行统计,包括:各数字型指标的数量count、均值mean、标准差std、最小值min、最大值max;以车流量的均值107作为高低流量的分界点,用sns.kdeplot函数画内核密度图查看车流量的分布状况:1)类目型指标week、weather上的高流量与低流量分布与量化,2)Week各天的高流量与低流量分布情况,3)交通指数Index上高流量与低流量分布情况,4)Time全天的高流量与低流量分布情况,5)各位置Place的高流量与低流量分布情况,6)volume、speed上的高流量与低流量分布情况;
利用sns.heatmap画图展示各指标间的关联度,关联度值的绝对值越大,两指标间的相互影响关系越深,用groupby函数得到与车流量Flow的其它指标关联度大小排序:交通量Volume(0.872605)、时间Time(0.533182)、速度Speed(0.328509)、路段Place(0.215406)、星期Week(0.064199)、交通指数Index(0.050326)、天气Weather(0.031915)。
可以得到Volume、Time、Speed这三个多值特征,与Flow关联度很大,可以做相同的数据处理,而Place、Week这两个离散值,与Flow关联度较低,也可以做相同的一类数据处理,排在最后的Index(一天交通指数的均值)、Weather这两项指标,在一天之内没有变化的特征值,同样可以做相同的一类数据处理。
根据关联度排序,将关联度小的week、weather合并为新指标地区region,如图2所示,将关联度高的连续性指标速度Speed、时间Time、交通量Volume分别划分为取值范围是4段、6段、8段,得到类目型指标Speed_group, Time_group, Volume_group。
使用get_dummies进行数据因子化:将类目型指标Place, Speed_group, Time_group, Volume_group平展开,取值均为0、1,比如当一个展开的Place指标取值为1时,它会确定Flow的一个模糊的取值范围,再结合其它指标的贡献性,逐渐收敛,一步步将Flow的取值范围缩小,最后确定预测值。最后Place平展为12个新指标, Speed_group平展为4个新指标, Time_group平展为6个新指标, Volume_group平展为8个新指标,由于特征值的取值只有0和1时,这会大大加快运算效率。完成数据因子化后,再通过drop函数将原先的group,Weather, Index, Speed这四个指标从训练数据集中去掉。
使用preprocessing.StandardScaler()进行数据标准化:注意Time和Volume的数据取值范围太大,这将对逻辑回归的收敛造成不利的影响。进行数据标准化,把这两个数据指标分布调整成标准正太分布,使得其均值为0,方差为1。
由于scikit-learn工具包为模型对象提供predict()接口,经过训练得出模型,可以用这个接口来进行预测,同时也提供了score()接口来评价一个模型的好坏,得分越高则模型预测性能越好,进行模型的选择和验证:将处理后的训练集的8种指标,输入scikit-learn工具包下六种常用机器学习算法进行模型训练,并得出各模型置信度:随机森林模型Random Forest(86.375000)、决策树模型Decision Tree(85.061000)、支持向量机模型Support Vector Machines(82.366000)、K近邻模型KNN(73.132000)、朴素贝叶斯模型Naive Bayes(58.350000)、逻辑回归模型Logistic Regression(35.813000)。
随机森林模型和决策树模型置信度最高,利用不同的机器学习算法的可调参数以及它们对机器学习算法的性能以及准确率的影响,下面对这两种模型进行优化:通过GridSearchCV()、StratifiedKFold()和StratifiedShuffleSplit()函数寻找最优参数下的决策树模型Decision_bes和最优参数下的随机森林模型RandomForest_best。
将处理后的训练集的8种指标,输入最优参数下的决策树模型Decision_bes和最优参数下的随机森林模型RandomForest_best,计算模型置信度得分,用置信度得分最高的最优参数下的随机森林模型RandomForest_best对测试集test.csv进行预测,最后以87%的准确性提交车流量预测结果fianl_submission1.csv。
为了一目了然获得训练数据集的样本大小是如何影响所选模型的置信度得分,避免欠拟合和过拟合情况的发生,需要利用优化参数后的随机森林模型画学习曲线:横坐标为训练集train.csv所训练样本的个数,纵坐标为训练集train.csv拟合的准确性得分和交叉验证集预测的准确性得分,画出学习曲线,验证当前训练样本的数量大小下,随机森林模型不仅能很好地拟合训练数据集,而且对测试数据集的预测准确性很高,即在保证准确度的同时也提升了模型防止过拟合的能力。
以南京理工大学7号门外丁字路口为实验城市道路车辆交通场景进行建模:配置车流文件rou.xml、道路文件net.xml、探测器loops Detectors等,通过改变车辆交通场景中的车辆数Vehicles,车速Speed,车道Lane和交叉口Junctions的设定,即是选择不同的主干路来构建道路交通模式,用交通仿真软件SUMO进行交通建模仿真,获取交通数据:交通量nVehContrib(采集周期内完整通过传感器的车辆数量)、车流量Flow(每小时通过传感器的车辆数量)、时间占有率occupancy(车辆通过的累计值与时间周期的比值)、平均速度speed、通过车辆的平均长度length、接触到传感器车辆的数量nVehEntered。
基于ubuntu16.04操作系统,完成通信仿真软件NS-3与交通仿真软件SUMO的连接,使NS-3能够遵循南京理工大学7号门外丁字路口的道路车辆交通场景输出通信数据,方法为:将车流文件.rou.xml和道路文件.net.xml合并转化为.sumo.tr文件;随后利用java -jar命令下的traceExporter.jar,将道路文件.net.xml和.sumo.tr文件合并转化为NS-3需要的.mobility.tcl文件,最后通过ns2-mobility-trace连接NS-3与SUMO。
配置车辆通信仿真文件dsdv-njust1.cc,要完成路由协议dsdv、车速Speed、通信结点个数Nodes、仿真时间Time、输出文本类型等配置,主要添加std::string CSVfileName= "njust-dsdv1.csv"语句和std::string tr_name = "Dsdv_Manet_" + t_nodes + "Nodes_" + sTotalTime + "SimTime";std::cout << "Trace file generated is " <<tr_name << ".tr\n"语句,运行车辆通信仿真文件dsdv-njust1.cc,获取通信数据:路由协议RoutingProtocol、接收到的数据包数目PacketsReceived、数据包接受率ReceiveRate和可供网络链路数据分析的网络抓包.pcap文件。
使用Wireshark的统计功能,选用Npcap Loopback Adapter滤波器对网络抓包.pcap文件进行分组封装,对通信结点链路和信道的数据包进行数据撷取,获取另一部分通信数据:数据包传送速率PacketsV和平均带宽Band。
对获取的交通流量和通信流量进行相关性分析:使用sns.catplot查看仿真时间内车流量Flow与数据包传送速率PacketsV、平均带宽Band的相关性,结果表明交通车流量和车联网通信流量存在一定程度的正相关。当车流量增大时,整体上通信网络流量随着增大,且数据传送质量要求越来越高,车载通信网络容易发生数据拥塞。
利用openstreetmap导出在交通数据平台所选的延安高架路段的城市道路车辆交通场景,配置车辆通信仿真文件dsdv-njust2.cc,获取通信数据。依据已得的交通车流量和车联网通信流量之间的相关性,从交通数据平台发布的这些路段流量速度数据和仿真通信数据中,选取9类指标:平均带宽Band、数据包传送速率PacketsV、车辆数目Vehicles、交叉口数目Junctions、道路数目Lane、速度Speed、数据包接受率ReceiveRate、接收到的数据包数目PacketsReceived、路由协议RoutingProtocol。
用groupby函数得到各指标与平均带宽Band的关联度排序:平均带宽Band(1.000000)、车辆数目Vehicles(0.895059)、接收到的数据包数目PacketsReceived(0.877838)、数据包接受率ReceiveRate(0.871777)、数据包传送速率PacketsV(0.655794)、交叉口数目Junctions(0.030816)、道路数目Lane(0.030816)、速度Speed(0.030816)、路由协议RoutingProtocol(NaN),其中NaN表示数据集中该指标取值唯一,无变化。
将9类指标经过图3所示处理:对关联度排序在最后的三种指标Speed、LaneJunctions,合并为新指标;对关联度排序在前四的连续性指标Vehicles、PacketsReceived、ReceiveRate、PacketsV,依次设定指标数值分成8、6、4、2段取值范围,利用pd.cut函数自动寻找各段交界点,以各范围边界划分连续性指标,获得类目型指标;再使用get_dummies对类目型指标进行数据因子化,将指标平展开;利用preprocessing.StandardScaler()函数,对取值超过50种的多样性数据指标Vehicles、PacketsV进行数据标准化,将数据指标分布调整成标准正太分布,使得其均值为0,方差为1。
将处理后的训练数据集的9类指标,输入scikit-learn工具包下十个经典机器学习算法进行模型训练,并计算各模型置信度:自助聚合模型Bagging(92.683)、决策树模型Decision Tree(90.244)、随机森林模型Random Forest(87.805)、朴素贝叶斯模型NaiveBayes(86.585)、梯度提升决策树模型GradientBoosting(84.146)、高斯过程模型GaussianProcess(81.707)、K近邻模型KNN(75.610)、逻辑回归模型Logistic Regression(48.780)、AdaBoostClassifier模型(48.780)、支持向量机模型Support Vector Machines(39.024)。用置信度最高的Bagging模型进行通信流量预测,最后以92%的准确性提交通信流量预测结果fianl_submission2.csv。
从本实施例可以看出,本发明能够从机器学习的常用的经典算法出发,不仅有效的对真实交通场景的未来全天候车流量进行了很好的预测,还研究了车流量等交通信息与车联网实时的通信质量之间的关系,模型泛化性能好,准确度较高,能为后期利用经济高效的数据分发提供可靠的车载通信分析方法,增强车辆用户驾驶的安全性。
Claims (4)
1.一种基于机器学习的车联网通信的流量预测方法,其特征在于步骤如下:
步骤1、车流量预测:利用交通数据平台发布的流量速度数据,选取8类指标进行数据处理,完成全天候的车流量预测;
步骤2、分析两种流量关系:利用openstreetmap导出某地城市道路车辆交通场景,配置车流探测文件,获取交通数据,再配置车辆通信仿真文件,获取通信数据,两类数据混合,分析车流量与通信流量二者关系;
步骤3、通信流量预测:在交通数据平台选择某些路段,利用openstreetmap导出这些路段的城市道路车辆交通场景后,配置车辆通信仿真文件,获取通信数据,依据已得的两种流量关系,从交通数据平台发布的流量速度数据和获取的仿真车辆通信数据中,选取9类指标进行数据处理,完成车载网络的通信流量预测。
2.根据权利要求1所述的基于机器学习的车联网通信的流量预测方法,其特征在于,步骤1所述选取8类指标进行数据处理,完成全天候的车流量预测的方法为:
1.1获取交通数据平台发布的某些路段一周全天候的数据,利用这些数据,预测这些路段未来某周某天全天候的车流量,即形成7天的训练数据集train.csv,1天的测试数据集test.csv;8类数据集指标包括6类数字型指标和2类类目型指标,数字型指标:车流量Flow、时间Time、速度Speed、交通量Volume、交通指数Index、路段Place,类目型指标:星期Week、天气Weather;
1.2利用isna函数判断各指标是否有缺失值,并用len函数对训练集数据train.csv进行统计,包括:各数字型指标的数量count、均值mean、标准差std、最小值min、最大值max ;
1.3以车流量的均值mean作为高低流量的分界点,用sns.kdeplot函数画内核密度图查看车流量的分布状况:1)类目型指标week、weather上的高流量与低流量分布与量化,2)Week各天的高流量与低流量分布情况,3)交通指数Index上高流量与低流量分布情况,4)Time全天的高流量与低流量分布情况,5)各位置Place上的高流量与低流量分布情况,6)volume、speed上的高流量与低流量分布情况;
1.4利用sns.heatmap画图展示各指标间的关联度,用groupby函数得到车流量Flow与其它指标的关联度排序;
1.5将关联度排序在最后的两种指标,合并为新指标;对关联度排序在前三的三种数值连续性指标,依次设定指标数值分成8、6、4段取值范围,利用pd.cut函数自动寻找各段交界点,以各范围边界点来划分数值连续性指标,获得类目型指标;再使用get_dummies对类目型指标进行数据因子化,将指标平展开;
1.6利用preprocessing.StandardScaler()函数,对取值超过100种的多样性数据指标进行数据标准化,将数据指标分布调整成标准正太分布,使得其均值为0,方差为1;
1.7利用scikit-learn工具包下的经典机器学习算法随机森林模型,对处理后的8类指标进行训练和预测,并给出置信度得分;再通过随机森林模型的可调参数对模型进行优化:通过GridSearchCV()、StratifiedKFold()和StratifiedShuffleSplit()函数寻找最优参数下模型;最后对测试集test.csv进行预测,并提交结果fianl_submission1.csv;
1.8利用优化参数后的随机森林模型画学习曲线:横坐标为训练集train.csv所训练样本的个数,纵坐标为训练集train.csv拟合的准确性得分和交叉验证集预测的准确性得分,画出学习曲线。
3.根据权利要求1所述的基于机器学习的车联网通信的流量预测方法,其特征在于,步骤2分析两种流量关系的方法为:
2.1基于开源平台openstreetmap导出某实地城市道路路口,作为实验城市道路车辆交通场景进行建模:配置车流文件rou.xml、道路文件net.xml、探测器loops Detectors文件;然后用交通仿真软件SUMO进行交通仿真,获取交通数据:交通量nVehContrib即采集周期内完整通过传感器的车辆数量、车流量Flow即每小时通过传感器的车辆数量、时间占有率occupancy即车辆通过的累计值与时间周期的比值、平均速度speed即采集周期内的平均车速、通过探测器的车辆平均长度length、接触到传感器车辆的数量nVehEntered;
2.2基于ubuntu16.04操作系统,完成通信仿真软件NS-3与交通仿真软件SUMO的连接,使NS-3能够遵循SUMO实验城市道路车辆交通场景模式输出通信数据,方法为:将车流文件.rou.xml和道路文件.net.xml合并转化为.sumo.tr文件;随后利用java -jar命令下的traceExporter.jar,将道路文件.net.xml和.sumo.tr文件合并转化为NS-3需要的.mobility.tcl文件,最后通过ns2-mobility-trace连接NS-3与SUMO;
2.3配置车辆通信仿真文件dsdv-njust1.cc,完成路由协议dsdv、车速Speed、通信结点个数Nodes、仿真时间Time、输出文本类型配置,添加std::string CSVfileName = "njust-dsdv1.csv"语句和std::string tr_name = "Dsdv_Manet_" + t_nodes + "Nodes_" +sTotalTime + "SimTime";std::cout << "Trace file generated is " << tr_name <<".tr\n"语句,运行车辆通信仿真文件,获取通信数据:路由协议RoutingProtocol、接收到的数据包数目PacketsReceived、数据包接受率ReceiveRate和可供网络链路数据分析的网络抓包.pcap文件;
2.4使用Wireshark的统计功能,选用网络适配滤波器对网络抓包.pcap文件进行分组封装,对通信结点链路和信道的数据包进行数据撷取,获取另一部分通信数据:数据包传送速率PacketsV和平均带宽Band;
2.5对获取到的车流量和通信流量进行相关性分析:使用sns.catplot查看仿真时间内,车流量Flow与数据包传送速率PacketsV、平均带宽Band的相关度,验证交通车流量和车联网通信流量之间的相关性。
4.根据权利要求1所述的基于机器学习的车联网通信的流量预测方法,其特征在于,步骤3完成通信流量预测的具体方法为:
3.1在交通数据平台选择某些路段,利用openstreetmap导出这些路段的城市道路车辆交通场景后,配置车辆通信仿真文件dsdv-njust2.cc:配置路由协议dsdv、车速Speed、通信结点个数Nodes、仿真时间Time、输出文本类型等,添加std::string CSVfileName = "njust-dsdv2.csv"语句和std::string tr_name = "Dsdv_Manet_" + t_nodes + "Nodes_" + sTotalTime + "SimTime";std::cout << "Trace file generated is " <<tr_name << ".tr\n"语句,最后运行车辆通信仿真文件dsdv-njust2.cc,获取仿真车辆通信数据;
3.2依据已得的交通车流量和车联网通信流量之间的相关性,从交通数据平台发布的这些路段流量速度数据和仿真车辆通信数据中,选取9类指标:平均带宽Band、数据包传送速率PacketsV、车辆数目Vehicles、交叉口数目Junctions、道路数目Lane、速度Speed、数据包接受率ReceiveRate、接收到的数据包数目PacketsReceived、路由协议RoutingProtocol;
3.3利用所选路段在不同车速Speed和不同通信结点个数Nodes下获取的9类数据指标,预测这些路段未来某周某天的通信流量,即形成训练数据集train.csv和测试数据集test.csv;
3.2对9类指标进行统计:利用isna函数判断各指标是否有缺失值,并用len函数对训练集数据train.csv进行统计,包括:各指标的数量count、均值mean、标准差std、最小值min、最大值max ;以平均带宽Band的均值mean作为高低通信流量的分界点,用sns.kdeplot函数画内核密度图查看通信流量的分布状况:1)速度Speed上的高流量与低流量分布,2)Lane各道路的高流量与低流量分布情况,3)车辆数目Vehicles上高流量与低流量分布情况,4)数据包传送速率PacketsV上的高流量与低流量分布情况,5)交叉口数目Junctions上的高流量与低流量分布情况,6)接收到的数据包数目PacketsReceived、数据包接受率ReceiveRate上的高流量与低流量分布情况;
3.3对9类指标进行处理:利用sns.heatmap画图展示各指标间的关联度,用groupby函数得到平均带宽Band与其它指标的关联度排序;将关联度排序在最后的三种指标,合并为新指标;对关联度排序在前四的数值连续性指标,依次设定指标数值分成8、6、4、2段取值范围,利用pd.cut函数自动寻找各段交界点,以各范围边界划分数值连续性指标,获得类目型指标;再使用get_dummies对类目型指标进行数据因子化,将指标平展开;
3.4利用preprocessing.StandardScaler()函数,对取值超过50种的多样性数据指标进行数据标准化,将数据指标分布调整成标准正太分布,使得其均值为0,方差为1
3.5利用scikit-learn工具包下的经典机器学习算法自助聚合BaggingClassifier模型对处理后的9类指标进行训练,然后对测试集test.csv进行预测,并提交通信流量预测结果fianl_submission2.csv。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910504819.4A CN110266528B (zh) | 2019-06-12 | 2019-06-12 | 基于机器学习的车联网通信的流量预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910504819.4A CN110266528B (zh) | 2019-06-12 | 2019-06-12 | 基于机器学习的车联网通信的流量预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110266528A true CN110266528A (zh) | 2019-09-20 |
CN110266528B CN110266528B (zh) | 2022-04-08 |
Family
ID=67917728
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910504819.4A Expired - Fee Related CN110266528B (zh) | 2019-06-12 | 2019-06-12 | 基于机器学习的车联网通信的流量预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110266528B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111696369A (zh) * | 2020-04-10 | 2020-09-22 | 北京数城未来科技有限公司 | 一种基于多源地理空间大数据的全市道路分时分车型交通流预测方法 |
CN112216104A (zh) * | 2020-09-17 | 2021-01-12 | 广东新时空科技股份有限公司 | 一种基于多源数据融合的城市路口交通流预测方法 |
CN113313937A (zh) * | 2021-04-14 | 2021-08-27 | 交通运输部规划研究院 | 一种基于土地利用信息的路网交通流量动态模拟方法 |
CN113487105A (zh) * | 2021-07-27 | 2021-10-08 | 成都理工大学 | 一种滑坡地质灾害的空间预测方法、装置和计算机设备 |
CN114125932A (zh) * | 2020-09-01 | 2022-03-01 | 中国移动通信集团湖南有限公司 | 数据分流的方法、装置及网络设备 |
CN114580086A (zh) * | 2022-05-05 | 2022-06-03 | 中汽研汽车检验中心(天津)有限公司 | 一种基于监督式机器学习的车辆部件建模方法 |
CN117455553A (zh) * | 2023-12-08 | 2024-01-26 | 青岛理工大学 | 一种地铁站客流量预测的方法 |
CN118277452A (zh) * | 2024-04-15 | 2024-07-02 | 中国科学院地理科学与资源研究所 | 基于要素流动大数据的边界效应量化方法及系统 |
CN118314727A (zh) * | 2024-04-22 | 2024-07-09 | 北京博派通达科技有限公司 | 基于机器学习的车辆交通事故实时检测方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110084854A1 (en) * | 2009-10-09 | 2011-04-14 | David Randal Johnson | Monitoring Management and Presentation of Preemption Control Data of Centrally Managed Traffic Signals |
CN102496284A (zh) * | 2011-12-24 | 2012-06-13 | 重庆大学 | 一种道路交通流量采集及预测方法 |
US9171041B1 (en) * | 2011-09-29 | 2015-10-27 | Pivotal Software, Inc. | RLE-aware optimization of SQL queries |
CN108364490A (zh) * | 2018-01-17 | 2018-08-03 | 佛山杰致信息科技有限公司 | 城市公路交通系统车辆运行调控方法 |
CN108900346A (zh) * | 2018-07-06 | 2018-11-27 | 西安电子科技大学 | 基于lstm网络的无线网络流量预测方法 |
-
2019
- 2019-06-12 CN CN201910504819.4A patent/CN110266528B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110084854A1 (en) * | 2009-10-09 | 2011-04-14 | David Randal Johnson | Monitoring Management and Presentation of Preemption Control Data of Centrally Managed Traffic Signals |
US9171041B1 (en) * | 2011-09-29 | 2015-10-27 | Pivotal Software, Inc. | RLE-aware optimization of SQL queries |
CN102496284A (zh) * | 2011-12-24 | 2012-06-13 | 重庆大学 | 一种道路交通流量采集及预测方法 |
CN108364490A (zh) * | 2018-01-17 | 2018-08-03 | 佛山杰致信息科技有限公司 | 城市公路交通系统车辆运行调控方法 |
CN108900346A (zh) * | 2018-07-06 | 2018-11-27 | 西安电子科技大学 | 基于lstm网络的无线网络流量预测方法 |
Non-Patent Citations (2)
Title |
---|
CHRISTIAN BACKFRIEDER等: "Increased Traffic Flow Through Node-Based Bottleneck Prediction and V2X Communication", 《IEEE TRANSACTIONS ON INTELLIGENT TRANSPORTATION SYSTEMS》 * |
赵海涛等: "基于连通概率感知的车联网资源优化技术研究", 《仪器仪表学报》 * |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111696369A (zh) * | 2020-04-10 | 2020-09-22 | 北京数城未来科技有限公司 | 一种基于多源地理空间大数据的全市道路分时分车型交通流预测方法 |
CN114125932B (zh) * | 2020-09-01 | 2023-08-15 | 中国移动通信集团湖南有限公司 | 数据分流的方法、装置及网络设备 |
CN114125932A (zh) * | 2020-09-01 | 2022-03-01 | 中国移动通信集团湖南有限公司 | 数据分流的方法、装置及网络设备 |
CN112216104A (zh) * | 2020-09-17 | 2021-01-12 | 广东新时空科技股份有限公司 | 一种基于多源数据融合的城市路口交通流预测方法 |
CN113313937A (zh) * | 2021-04-14 | 2021-08-27 | 交通运输部规划研究院 | 一种基于土地利用信息的路网交通流量动态模拟方法 |
CN113313937B (zh) * | 2021-04-14 | 2022-05-06 | 交通运输部规划研究院 | 一种基于土地利用信息的路网交通流量动态模拟方法 |
CN113487105B (zh) * | 2021-07-27 | 2023-04-07 | 成都理工大学 | 一种滑坡地质灾害的空间预测方法、装置和计算机设备 |
CN113487105A (zh) * | 2021-07-27 | 2021-10-08 | 成都理工大学 | 一种滑坡地质灾害的空间预测方法、装置和计算机设备 |
CN114580086B (zh) * | 2022-05-05 | 2022-08-09 | 中汽研汽车检验中心(天津)有限公司 | 一种基于监督式机器学习的车辆部件建模方法 |
CN114580086A (zh) * | 2022-05-05 | 2022-06-03 | 中汽研汽车检验中心(天津)有限公司 | 一种基于监督式机器学习的车辆部件建模方法 |
CN117455553A (zh) * | 2023-12-08 | 2024-01-26 | 青岛理工大学 | 一种地铁站客流量预测的方法 |
CN117455553B (zh) * | 2023-12-08 | 2024-04-30 | 青岛理工大学 | 一种地铁站客流量预测的方法 |
CN118277452A (zh) * | 2024-04-15 | 2024-07-02 | 中国科学院地理科学与资源研究所 | 基于要素流动大数据的边界效应量化方法及系统 |
CN118277452B (zh) * | 2024-04-15 | 2024-09-24 | 中国科学院地理科学与资源研究所 | 基于要素流动大数据的边界效应量化方法及系统 |
CN118314727A (zh) * | 2024-04-22 | 2024-07-09 | 北京博派通达科技有限公司 | 基于机器学习的车辆交通事故实时检测方法及系统 |
CN118314727B (zh) * | 2024-04-22 | 2024-09-24 | 北京博派通达科技有限公司 | 基于机器学习的车辆交通事故实时检测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110266528B (zh) | 2022-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110266528A (zh) | 基于机器学习的车联网通信的流量预测方法 | |
CN104464321B (zh) | 一种基于交通运行指数发展态势的智能交通诱导方法 | |
CN104778834B (zh) | 一种基于车辆gps数据的城市道路交通拥堵判别方法 | |
CN110634287B (zh) | 一种基于边缘计算的城市交通状态精细化判别方法 | |
CN108629974B (zh) | 顾及城市道路交通网络特征的交通运行指数建立方法 | |
CN111341095B (zh) | 一种基于边缘侧在线计算的交通信号控制系统及方法 | |
CN105679025B (zh) | 一种基于可变权重混合分布的城市干道行程时间估计方法 | |
WO2019047905A1 (zh) | 一种交通路况分析系统、方法以及装置 | |
CN107248283A (zh) | 一种考虑路段关键度的城市区域路网运行状态评价方法 | |
CN104021671B (zh) | 一种svm与模糊判断相结合的道路实时路况的判断方法 | |
CN111210621B (zh) | 基于实时路况的信号绿波协调路线优化控制方法及系统 | |
CN108072381A (zh) | 一种路径规划的方法及装置 | |
CN109615851B (zh) | 一种在群智感知系统中基于关键路段的感知节点选取方法 | |
Chen et al. | Discrimination and prediction of traffic congestion states of urban road network based on spatio-temporal correlation | |
CN105743783A (zh) | 基于bs-ts和自编码网络的车联网网络节点筛选方法及其通达性路由机制 | |
CN114708746A (zh) | 智慧城市的交通信号提示方法及系统 | |
Zhang et al. | A hybrid method of traffic congestion prediction and control | |
CN109544967A (zh) | 一种基于低频avl数据的公交线网运行状态监控方法 | |
CN108091131A (zh) | 一种交通事件识别方法及装置 | |
CN107644533A (zh) | 基于移动网络数据的高速公路虚拟断面车流量监测方法 | |
Mohammadi et al. | User throughput optimization for signalized intersection in a connected vehicle environment | |
CN108171367A (zh) | 一种公交车服务水平可靠性预测方法 | |
CN115662121B (zh) | 混行信控交叉口以排队方式估算通行能力上下限的方法 | |
CN105654415A (zh) | 一种面向交通管理者的路网通行效率变化率计算方法 | |
CN105243439A (zh) | 基于路网结构的城市群城际铁路线网规模测算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220408 |