专利名称：驾驶员对路径信息板动态响应行为的动态模拟方法
技术领域：
本发明涉及驾驶员出行行为的计算机模拟技术，特别涉及一种驾驶员对 路径信息板动态响应行为的动态模拟方法。 '
背景技术：
目前，应用计算机模拟技术，对驾驶员在提供实时交通信息条件下的信 息响应行为(例如，路径选择、出发时刻选择)进行模拟，收集驾驶员行为 数据，是智能交通运输研究领域的前沿和热点，相关文献报道已经不少。研 究成果用途很多，包括为驾驶员行为建模提供基础数据，更有效地设计交 通信息服务产品，制定更科学的信息发布策略，对交通信息诱导系统的效益 进行建模和定量分析。以往研究中，学者们主要结合所在国家的智能交通技 术特点，针对某种或某几种交通信息发布手段(如可变信息标志、交通广 播、车载导航设备)，进行研究。
近年来，亚洲和欧洲一些国家开始在快速路上应用或实验以地图(发光 元件组成的光带)形式显示快速路交通状态的图状路径信息板(英文称 GRIP, graphical route information panel),诱导交通流。也有文献称 GRIP为图形式可变情报板。GRIP这一项新技术的应用对交通科学界提出了 一个新挑战"如何正确理解GRIP影响下的驾驶员动态路径选由行为，并对 其进行建模？"目前，关于GRIP影响下驾驶员动态路径选择行为的研究主 要还停留在理论研究初级阶段。纵观现有文献，尚未见有文献提出针对GRIP 的驾驶员路径选择行为动态模拟方法。理论研究的滞后，使得工程实践中无 法科学地回答的"该怎样进行GRIP路网优化布局？ GRIP的交通诱导效何量化评价？"等关键问题。
发明内容
本发明是针对上述现有技术存在的问题，提出了一种驾驶员对路径信息 板动态响应行为的动态模拟方法。该方法利用人机交互的计算机模拟系统，
模拟驾驶员在提供GRIP信息条件下的动态出行决策过程，出行决策包括出 发时的路径选择和出行途中的改道行为。模拟系统为驾驶员提供较为逼真的 快速路出行情境和满足人际功效学原理的人机交互界面，驾驶员在虚拟出行 情境中出行，模拟系统自动记录驾驶员在整个出行过程中的决策行为，收集 丰富的行为数据，为研究人员探索GRIP信息与驾驶员行为之间的动态复杂 关系以及建立行为模型提供基础数据。
本发明的技术方案为 一种驾驶员对路径信息板动态响应行为的动态模 拟方法，方法包括如下步骤
1)利用计算机作为界面，应用交通流METANET模型，建立快速路网宏观交 通流仿真数学模型路网由节点和有向路段组成，节点为存在出(入)口匝道 或道路属性发生变化的位置；路段是属性，如车道数、坡度，相同的一段高 速公路；路段分为若干等长小段，路段^分为^个长度为A^^的小段，小段 车道数记为、，路段历的第i小段在A时刻的交通状态用密度、速度和流量 表示，分别记作A^(的、vW和^(A:)，给出METANET的基本数学公式 A:描述小段交通状态变化的公式如下 <formula>formula see original document page 6</formula><formula>formula see original document page 7</formula>
公式中，/V,，,W表示A^("中终点为J的分密度；7为模型取样周期； ;V,,,~,y//Vy，称为交通组成比例；r，v，K是常量参数，^(.)为稳态速度-密 度关系； 为路段/z 的自由流速度，p，是路段历的临界密度，^是与路段
"有关的参数；
B:节点处的交通流运行特征由以下公式描述
>7历
公式中，/"为流入节点n的路段集合，O"为从节点n流出的路段集合；人为
经路段m能够到达的终点集合；a，,ot)为从节点出发驶往终点J的流量；
为a,^)中选择路段m的流量所占比例，称作分流比例；u()t)为路段m入口 处流量，？V。，/W为^。")中流向终点J'的比例；
C:路网起点，包括主线起始点和入口匝道的车辆排队现象，由以下公式表
小
<formula>formula see original document page 7</formula>
公式中，M^.为起点0处车辆数W。中驶往目的地J'的数量，W。=;W。，；；《是O点交通需求；^是0D比例，是《中驶往目的地J的车流所占比例。《。是
实际流入路网的流量；y。，d是9。中驶往目的地J'的车流所占比例；r。"0，l]为 匝道调节率；G。为起点o的通行能力，义。为起点o的车道数，^一是阻塞密 度。A,为起点o下游紧邻路段s的第一小段的密度；
2 )已知路网初始时刻(k=0)的Pm,,,; (0)、 (0)、 (0)，已知所有时刻的 《W、《"Q、 A^(A:)、 r。(W，根据建立的快速路网宏观交通流仿真数学模 型；驾驶员从路网的起点进入路网，驾驶员所驾驶的车辆就根据建立的快速 路网宏观交通流仿真数学模型计算出的平均车速，自动向前行驶；
3)在快速路网宏观交通流仿真数学模型基础上，建立功能模块，同时在计 算机界面上对应功能视窗-
A:交通流运行状况视窗视窗中运动的小方格的数量根据小段密度确定， 小方格稠密程度表征了道路拥堵程度，驾驶员的车辆在运动的小方格中，驾 驶员通过操作避开运动的小方格；
B : GRIP信息视窗GRIP上的快速路分割为若干信息单元段，信息单元段 就是驾驶员能够辨识出的像素，其显示状态有畅通为绿色、拥挤为黄色、阻 塞为红色三种；
C:基本出行信息视窗出发时间、当前时间、已行驶距离、路网地图，路 网地图驾驶员可根据需要开启或关闭；
D:车辆方位信息视窗显示驾驶员所在的局部路网和车辆在该路网中的位 置，当车辆临近快速道路转换点，如立交、出口匝道时视窗中会提供前方相 关道路的信息，模拟现实世界中的指路标志和预告标志； E:其他辅助信息视窗标示速度的汽车仪表盘；4)驾驶员在计算机界面营造的出行情境视窗中出行，出行途中，接收到分 流点处的GRIP信息后做出路径选择，通过鼠标和键盘把所作的选择输入给 计算机，驾驶员的车辆继续前进，直至到达目的地，系统自动记录整个出行 过程中驾驶员作出的动态响应行为，包括出行前的路径选择和出行途中的路 径改换，行为数据存储于数据库。
本发明的有益效果在于本发明驾驶员对路径信息板动态响应行为的动 态模拟方法，能够为驾驶员提供符合网络交通流动态演化特征的较逼真的快 速路网虚拟出行情境，具有高效的人机界面设计，能够采集丰富的驾驶员信 息响应行为数据；较之高逼真的模拟驾车器，系统开发成本低得多；模拟系 统收集数据简便易操作，可以在台式机上进行，也可以用便携式笔记本进 行，试验场地不受限制；试验过程快速，仅需几分钟；试验成本较低，克服 了传统的问巻调査"出行情景设计不宜太复杂、问题不宜过多、易引起驾驶 员疲劳"等不足。收集的驾驶员信息响应行为数据为驾驶员行为的分析与建 模提供基础数据。
图1为本发明驾驶员对路径信息板动态响应行为的动态模拟方法中功能模块 示意图2为本发明驾驶员对路径信息板动态响应行为的动态模拟系统中出行情境
人机界面示意具体实施例方式
如图1所示驾驶员对图状路径信息板动态响应行为的计算机模拟方法中 功能模块示意图，驾驶员和计算机之间互动，驾驶员可从快速路网虚拟出行情境的计算机界面得到需要信息，也可将实际的信息或需要得到的信息输入
快速路网虚拟出行情境的计算机界面，如图2所示示意图，快速路网虚拟出
行情境的计算机界面由以下几个功能模块支持基本出行信息生成模块、交
通流运行状态生成模块、车辆方位信息生成模块、驾驶员任务模块、GRIP信
息生成模块、辅助信息模块，上述功能模块的信息通过快速路网宏观交通流 仿真模型从数据库中调入，同时功能模块将信息通过快速路网宏观交通流仿 真模型存入数据库。 .
下面为驾驶员对图状路径信息板动态响应行为的计算机模拟方法的实现
1、建立快速路网宏观交通流仿真数学模型(简称网络流仿真模型)。 应用交通流理论、软件工程等知识，建立快速路网宏观交通流仿真数学
模型。本发明运用经典的METANET模型(由Papageorgiou等学者提出)，构 建快速路网宏观交通流仿真数学模型。
METANET中，路网由节点和有向路段组成。节点为存在出(入)口匝道或 道路属性发生变化的位置。路段是属性(如车道数、坡度)相同的一段高速公 路。路段分为若干等长小段。例如，路段^分为^个长度为Am的小段，小 段车道数记为A。路段辺的第i小段在A时刻的交通状态用密度、速度和流 量表示，分别记作/v,询、^W和^,询。以下给出METANET的基本数学公 式。
描述小段交通状态变化的公式如下。
&,」+1) = (A:) + (77 ;i加Am (A:)gm，M (A) — W^,i W] ' (1)
^,,询=&,(^ #)^ (2) ("1)= ("+—『五( W) _ V, ("] + T' v加,'(".[V"— v加,'("]
T m (3)
v.r pm，,+1("-/y,("
K(/V,仰=v,,m exp(—丄(^i^户) (4)
以上公式中，A^W表示/v"W中终点为J的分密度；7为模型取样周 期；^ZA^/zv"称为交通组成比例；r,v,K是常量参数，&(.)为稳态速 度-密度关系； 为路段仿的自由流速度，户，是路段历的临界密度，&是 与路段a有关的参数。 ' 路网节点处的交通流运行特征的描述如下
a,#)= s w"y w v(",力 (5)
^,0(Q= i: /^;(A:)a,/^ VmeO" (6)
Vo,j (" = d(" m,。 W Vm ee厶 (7)
以上公式中，/ 为流入节点n的路段集合，"为从节点n流出的路段集 合；人为经路段m能够到达的终点集合；么，,.(W为从节点出发驶往终点j的
流量；/C为K"中选择路段m的流量所占比例，称作分流比例；UW为 路段m入口处流量，乙，。，,(fc)为^，。("中流向终点j的比例。
路网起点(包括主线起始点和入口匝道)的车辆排队现象，由以下公式
w。 ,乂 (A: +1) = w。,y W + 7[《,)(" -W《。("] (8) ^W = r。Ot).minK^) + ，,a .义。.min{l,P,—(9)
以上公式中，W。j为起点o处车辆数w。中驶往目的地乂的数量，w。=2>w;夂是O点交通需求；《w是0D比例，是《中驶往目的地J的车
乂
流所占比例。《。是实际流入路网的流量；7。，,是《。中驶往目的地'J的车流所 占比例；r。e
为匝道调节率；a为起点O的通行能力，A。为起点0的车道 数，p，是阻塞密度。ZV为起点O下游紧邻路段S的第一小段的密度。
针对某 一 个具体的快速路网，只要己知初始时刻(k=0)的&。. (0)、 、"0)、 w。,(0)，已知所有时刻的J。W、《jW、《,(W、 r。(A)，就可以模拟 出任意时刻路网中任意位置(小段)的平均车速、流量、密度。
一旦驾驶员从路网的起点进入路网，驾驶员所驾驶的车辆就根据网络流 仿真模型计算出的平均车速自动向前行驶。 '
此外，网络流仿真模型还能够模拟事故的发生。例如，在某个时段内， 通过改变路网中某小段的车道数，就能近似模拟该小段发生事故的情形。
2、 在网络流仿真模型基础上，开发交通流运行状况生成模块。交通流 运行状况生成模块的基本工作原理是车辆在路网中行驶过程中，会遇到不 同的交通状况，即遇到不同的车流密度，该车流密度由网络流仿真模型计算 得到；(注密度的物理意义是驾驶员当前所在位置(小段)上每公里的车辆 数。)假设交通流运行状况视窗中呈现的道路长度是X公里，那么视窗中运 动的小方格的数量即为密度值与道路长度X的乘积。如图2示意图中的交通 流运行状况视窗③。
3、 建立驾驶员驾车任务模块。
交通流运行状况视窗中的运动小方格(驾驶员所驾车辆周围的车辆)由 计算机自动产生，匀速地向前运动。驾驶员所驾车辆(用黑色小方格表示) 被周边车辆所包围，需要驾驶员通过键盘对其进行操控。驾驶员通过键盘上的指定键位确保车辆保持前进状态，并避免与周边车辆相撞。如图2所示示 意图中的驾驶员所操控"车辆"⑥。
4、 建立GRIP信息生成模块。
GRIP信息生成模块运行原理是模拟现实世界中的GRIP信息生成方 法，GRIP上的快速路分割为若干信息单元段，信息单元段就是驾驶员能够辨 识出的像素，其显示状态有畅通(绿色)、拥挤(黄色)、阻塞(红色)三 种。信息单元段的划分与小段的划分一致。该模块以网络流仿真模型计算的 小段速度为依据，根据一定规则确定信息单元段的显示状态，并每隔一定间 隔(GRIP更新周期)更新显示状态。
确定信息单元段显示状态的规则为假设GRIP更新周期为Tvnis， Tvms 包含P个模型取样周期T，设当前时刻是k时刻，那么我们把k时刻前p个 T时间内的小段速度进行平均，将平均值对照预先设定的"速度-显示状态映 射关系"，确定信息单元段的显示状态。例如，速度-显示状态映射关系为 速度在20公里/小时以下发布红色，速度在20-40公里/小时情况下发布黄 色，速度大于40公里/小时情况下发布绿色。举例的规则与上海快速路上的 GRIP基本一致。根据实际需要，显示状态确定规则也可以采用第它方法。如 图2所示示意图中的GRIP信息视窗④。
5、 建立基本出行信息生成模块。
基本出行信息模块生成出发时间、当前时间、已行驶距离、路网地图等 信息。出发时间就是虚拟出行开始的时刻，由用户预先设定(例如，假定出 发时间是早上6点50分)。当前时间随着虚拟出行的推进，由计算机自动计 算。己行驶距离也是由计算机自动计算。需要注意的事，虚拟出行中时间推进快慢(即时钟快慢)与现实世界的时钟不同，用户可以自行设定虚拟出行
的时钟快慢，例如，假设现实世界中的10秒钟对应虚拟出行中的1秒钟。
这样虚拟出行的时间相对于实际出行，可以大大縮短。这也是本发明优于真 车模拟驾驶(类似玩真车驾驶游戏)的一大特点。
路网地图在计算机界面中可以有开启和关闭两种模式。开启模式用于模
拟车辆中有其他乘客，其他乘客可以帮忙看地图。如图2示意图所示基本出 行信息视窗①。
该模块生成的信息通过计算机界面提供给驾驶员。
6、 建立车辆方位信息生成模块。
该模块用于生成车辆方位信息，并通过计算机界面向驾驶员发布。信息 包括驾驶员所在的局部快速路网和车辆在该局部路网中的位置、以及相关注 释信息。
这些信息显示在位于计算机屏幕左上方的视窗中。视窗中通常显示驾驶 员所在的局部路网和车辆在该路网中的位置，当车辆临近快速道路转换点 (如立交、出口匝道)时视窗中还会提供前方相关道路的信息，以模拟现实 世界中的指路标志和预告标志等。如图2所示示意图中的车辆方位信息视窗 ②；
7、 建立其他辅助信息生成模块。
该模块用于生成其他辅助信息。信息包括标示速度的汽车仪表盘、发动 机轰鸣声，预告前方GRIP的语音提示等。仪表盘显示的速度近似为车辆当 前所处位置上的平均车速，由网络流仿真模型计算得到。发动机轰鸣声的响 度与车速有关系，该关系通过对车辆发动机声音的实际采样和编辑后得到。预告前方GRIP的提示音模拟驾车过程中驾驶员在可视范围内发现GRIP,提 示音在临近GRIP的位置响起。如图2所示示意图中的标示速度的汽车仪表 盘⑤。
8、 进行出行情境人机界面的开发。
采用软件工程、图形学等知识，应用面向对象的程序语言Visual 0++进 行人机界面开发，设计出符合人机工效学原理的出行情境人机界面，确保以 上各模块生成的信息在计算机屏幕中有合理的布局。
9、 构建驾驶员对图状路径信息板动态响应行为的计算机模拟系统。 以上述介绍的网络流仿真模型和各个模块为主要程序模块，应用面向对
象的程序语言Visual C++，构筑驾驶员对图状路径信息板动态响应行为的计 算机模拟系统。驾驶员在模拟系统的计算机界面营造的出行情境中出行，出 行途中，接收到分流点处的GRIP信息后做出路径选择，通过鼠标和键盘把 所作的选择输入给计算机，驾驶员的车辆继续前进，直至到达目的地。系统 自动记录整个出行过程中驾驶员作出的动态响应行为(包括出行前的路径选 择和出行途中的路径改换)，行为数据存储于数据库。该系统在模拟出行前 会询问和记录下驾驶员信息(年龄、性别、驾龄、收入、教育程度等)；在模 拟出行过程中会自动记录驾驶员的动态响应行为(包括出发前的路径选择和 出行途中看到GRIP信息后的路径改换等)，自动记录下驾驶员阅读GRIP信 息的时间；在模拟出行后记录驾驶员对GRIP的有用性和可理解性的评分， 并且存储其它出行相关数据，如，路径行程时间，道路信息、出行目的、出 行起讫点、事故信息、GRIP信息图像等。
权利要求
1、一种驾驶员对路径信息板动态响应行为的动态模拟方法，其特征在于，方法包括如下步骤1)利用计算机作为界面，应用交通流METANET模型，建立快速路网宏观交通流仿真数学模型路网由节点和有向路段组成，节点为存在出(入)口匝道或道路属性发生变化的位置；路段是属性，如车道数、坡度，相同的一段高速公路；路段分为若干等长小段，路段m分为Nm个长度为Δm的小段，小段车道数记为λm，路段m的第i小段在k时刻的交通状态用密度、速度和流量表示，分别记作ρm，i(k)、vm，i(k)和qm，i(k)，给出METANET的基本数学公式A描述小段交通状态变化的公式如下ρm，i，j(k+1)＝ρm，i，j(k)+(T/λmΔm)[γm，i-1，j(k)qm，i-1(k)-γm，i，j(k)qm，i(k)]qm，i(k)＝ρm，i(k)vm，i(k)λm<maths id="math0001" num="0001" ><math><![CDATA[ <mrow><msub> <mi>v</mi> <mrow><mi>m</mi><mo>,</mo><mi>i</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo></mrow><mo>=</mo><msub> <mi>v</mi> <mrow><mi>m</mi><mo>,</mo><mi>i</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><mfrac> <mi>T</mi> <mi>τ</mi></mfrac><mo>[</mo><msub> <mi>V</mi> <mi>E</mi></msub><mrow> <mo>(</mo> <msub><mi>ρ</mi><mrow> <mi>m</mi> <mo>,</mo> <mi>i</mi></mrow> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>-</mo><msub> <mi>v</mi> <mrow><mi>m</mi><mo>,</mo><mi>i</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>]</mo><mo>+</mo><mfrac> <mi>T</mi> <msub><mi>Δ</mi><mi>m</mi> </msub></mfrac><mo>·</mo><msub> <mi>v</mi> <mrow><mi>m</mi><mo>,</mo><mi>i</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>·</mo><mo>[</mo><msub> <mi>v</mi> <mrow><mi>m</mi><mo>,</mo><mi>i</mi><mo>-</mo><mn>1</mn> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>v</mi> <mrow><mi>m</mi><mo>,</mo><mi>i</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>]</mo> </mrow>]]></math></maths><maths id="math0002" num="0002" ><math><![CDATA[ <mrow><mo>-</mo><mfrac> <mrow><mi>v</mi><mo>·</mo><mi>T</mi> </mrow> <mrow><mi>τ</mi><mo>·</mo><msub> <mi>Δ</mi> <mi>m</mi></msub> </mrow></mfrac><mo>·</mo><mfrac> <mrow><msub> <mi>ρ</mi> <mrow><mi>m</mi><mo>,</mo><mi>i</mi><mo>+</mo><mn>1</mn> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>-</mo><msub> <mi>ρ</mi> <mrow><mi>m</mi><mo>,</mo><mi>i</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow> <mrow><msub> <mi>ρ</mi> <mrow><mi>m</mi><mo>,</mo><mi>i</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><mi>κ</mi> </mrow></mfrac> </mrow>]]></math></maths><maths id="math0003" num="0003" ><math><![CDATA[ <mrow><msub> <mi>V</mi> <mi>E</mi></msub><mrow> <mo>(</mo> <msub><mi>ρ</mi><mrow> <mi>m</mi> <mo>,</mo> <mi>i</mi></mrow> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow> <mo>)</mo></mrow><mo>=</mo><msub> <mi>v</mi> <mrow><mi>f</mi><mo>,</mo><mi>m</mi> </mrow></msub><mi>exp</mi><mrow> <mo>(</mo> <mo>-</mo> <mfrac><mn>1</mn><msub> <mi>a</mi> <mi>m</mi></msub> </mfrac> <msup><mrow> <mo>(</mo> <mfrac><mrow> <msub><mi>ρ</mi><mrow> <mi>m</mi> <mo>,</mo> <mi>i</mi></mrow> </msub> <mrow><mo>(</mo><mi>k</mi><mo>)</mo> </mrow></mrow><msub> <mi>ρ</mi> <mrow><mi>cr</mi><mo>,</mo><mi>m</mi> </mrow></msub> </mfrac> <mo>)</mo></mrow><msub> <mi>a</mi> <mi>m</mi></msub> </msup> <mo>)</mo></mrow> </mrow>]]></math></maths>公式中，ρm，i，j(k)表示ρm，i，j(k)中终点为j的分密度；T为模型取样周期；γm，i，j＝ρm，i，j/ρm，j，称为交通组成比例；τ，v，κ是常量参数，VE(·)为稳态速度-密度关系；vf，m为路段m的自由流速度，ρcr，m是路段m的临界密度，am是与路段m有关的参数；B节点处的交通流运行特征由以下公式描述<maths id="math0004" num="0004" ><math><![CDATA[ <mrow><msub> <mi>Q</mi> <mrow><mi>n</mi><mo>,</mo><mi>j</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><munder> <mi>Σ</mi> <mrow><mi>μ</mi><mo>∈</mo><msub> <mi>I</mi> <mi>n</mi></msub> </mrow></munder><msub> <mi>q</mi> <mrow><mi>μ</mi><mo>,</mo><msub> <mi>N</mi> <mi>μ</mi></msub> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><msub> <mi>γ</mi> <mrow><mi>μ</mi><mo>,</mo><msub> <mi>N</mi> <mi>μ</mi></msub><mo>,</mo><mi>j</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow>]]></math> id="icf0004" file="A2009100514620002C4.tif" wi="52" he="7" top= "253" left = "25" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths><maths id="math0005" num="0005" ><math><![CDATA[ <mrow><mo>∀</mo><mrow> <mo>(</mo> <mi>n</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo></mrow> </mrow>]]></math> id="icf0005" file="A2009100514620002C5.tif" wi="11" he="4" top= "253" left = "84" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths><maths id="math0006" num="0006" ><math><![CDATA[ <mrow><msub> <mi>q</mi> <mrow><mi>m</mi><mo>,</mo><mn>0</mn> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><munder> <mi>Σ</mi> <mrow><mi>j</mi><mo>∈</mo><msub> <mi>J</mi> <mi>m</mi></msub> </mrow></munder><msubsup> <mi>β</mi> <mrow><mi>n</mi><mo>,</mo><mi>j</mi> </mrow> <mi>m</mi></msubsup><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><msub> <mi>Q</mi> <mrow><mi>n</mi><mo>,</mo><mi>j</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow>]]></math> id="icf0006" file="A2009100514620002C6.tif" wi="46" he="7" top= "264" left = "25" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths><maths id="math0007" num="0007" ><math><![CDATA[ <mrow><mo>∀</mo><mi>m</mi><mo>∈</mo><msub> <mi>O</mi> <mi>n</mi></msub> </mrow>]]></math> id="icf0007" file="A2009100514620002C7.tif" wi="14" he="4" top= "264" left = "88" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths><maths id="math0008" num="0008" ><math><![CDATA[ <mrow><msub> <mi>γ</mi> <mrow><mi>m</mi><mo>,</mo><mn>0</mn><mo>,</mo><mi>j</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><msubsup> <mi>β</mi> <mrow><mi>n</mi><mo>,</mo><mi>j</mi> </mrow> <mi>m</mi></msubsup><msub> <mi>Q</mi> <mrow><mi>n</mi><mo>,</mo><mi>j</mi> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>/</mo><msub> <mi>q</mi> <mrow><mi>m</mi><mo>,</mo><mn>0</mn> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow>]]></math> id="icf0008" file="A2009100514620003C1.tif" wi="50" he="4" top= "30" left = "26" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths><maths id="math0009" num="0009" ><math><![CDATA[ <mrow><mo>∀</mo><mi>m</mi><mo>∈</mo><msub> <mi>O</mi> <mi>n</mi></msub><mo>,</mo><mo>∀</mo><mi>j</mi><mo>∈</mo><msub> <mi>J</mi> <mi>m</mi></msub> </mrow>]]></math> id="icf0009" file="A2009100514620003C2.tif" wi="28" he="4" top= "30" left = "91" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>公式中，In为流入节点n的路段集合，On为从节点n流出的路段集合；Jm为经路段m能够到达的终点集合；Qn，j(k)为从节点出发驶往终点j的流量；βn，jm为Qn，j(k)中选择路段m的流量所占比例，称作分流比例；qm，0(k)为路段m入口处流量，γm，0，j(k)为qm，0(k)中流向终点j的比例；C路网起点，包括主线起始点和入口匝道的车辆排队现象，由以下公式表示wo，j(k+1)＝wo，j(k)+T[θo，j(k)do(k)-γo，j(k)qo(k)]<maths id="math0010" num="0010" ><math><![CDATA[ <mrow><msub> <mi>q</mi> <mi>o</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>=</mo><msub> <mi>r</mi> <mi>o</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>·</mo><mi>min</mi><mo>{</mo><msub> <mi>d</mi> <mi>o</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow><mo>+</mo><mfrac> <mrow><msub> <mi>w</mi> <mi>o</mi></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow> <mi>T</mi></mfrac><mo>,</mo><msub> <mi>Q</mi> <mi>o</mi></msub><mo>·</mo><msub> <mi>λ</mi> <mi>o</mi></msub><mo>·</mo><mi>min</mi><mo>{</mo><mn>1</mn><mo>,</mo><mfrac> <mrow><msub> <mi>ρ</mi> <mi>jam</mi></msub><mo>-</mo><msub> <mi>ρ</mi> <mrow><mi>s</mi><mo>,</mo><mn>1</mn> </mrow></msub><mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo></mrow> </mrow> <mrow><msub> <mi>ρ</mi> <mi>jam</mi></msub><mo>-</mo><msub> <mi>ρ</mi> <mrow><mi>cr</mi><mo>,</mo><mi>s</mi> </mrow></msub> </mrow></mfrac><mo>}</mo><mo>}</mo> </mrow>]]></math></maths>公式中，wo，j为起点o处车辆数wo中驶往目的地j的数量，<maths id="math0011" num="0011" ><math><![CDATA[ <mrow><msub> <mi>w</mi> <mi>o</mi></msub><mo>=</mo><munder> <mi>Σ</mi> <mi>j</mi></munder><msub> <mi>w</mi> <mrow><mi>o</mi><mo>,</mo><mi>j</mi> </mrow></msub><mo>;</mo> </mrow>]]></math> id="icf0011" file="A2009100514620003C4.tif" wi="23" he="8" top= "133" left = "150" img-content="drawing" img-format="tif" orientation="portrait" inline="yes"/></maths>do是o点交通需求；θo，d是OD比例，是do中驶往目的地j的车流所占比例。qo是实际流入路网的流量；γo，d是qo中驶往目的地j的车流所占比例；ro∈
为匝道调节率；Qo为起点o的通行能力，λo为起点o的车道数，ρjam是阻塞密度。ρs，1为起点o下游紧邻路段s的第一小段的密度；2)已知路网初始时刻(k＝0)的ρm，i，j(0)、vm，i(0)、wo，j(0)，已知所有时刻的do(k)、θo，j(k)、βn，jm(k)、ro(k)，根据建立的快速路网宏观交通流仿真数学模型；驾驶员从路网的起点进入路网，驾驶员所驾驶的车辆就根据建立的快速路网宏观交通流仿真数学模型计算出的平均车速，自动向前行驶；3)在快速路网宏观交通流仿真数学模型基础上，建立功能模块，同时在计算机界面上对应功能视窗A交通流运行状况视窗视窗中运动的小方格的数量根据小段密度确定，小方格稠密程度表征了道路拥堵程度，驾驶员的车辆在运动的小方格中，驾驶员通过操作避开运动的小方格；BGRIP信息视窗GRIP上的快速路分割为若干信息单元段，信息单元段就是驾驶员能够辨识出的像素，其显示状态有畅通为绿色、拥挤为黄色、阻塞为红色三种；C基本出行信息视窗出发时间、当前时间、已行驶距离、路网地图，路网地图驾驶员可根据需要开启或关闭；D车辆方位信息视窗显示驾驶员所在的局部路网和车辆在该路网中的位置，当车辆临近快速道路转换点，如立交、出口匝道时视窗中会提供前方相关道路的信息，模拟现实世界中的指路标志和预告标志；E其他辅助信息视窗标示速度的汽车仪表盘；4)驾驶员在计算机界面营造的出行情境视窗中出行，出行途中，接收到分流点处的GRIP信息后做出路径选择，通过鼠标和键盘把所作的选择输入给计算机，驾驶员的车辆继续前进，直至到达目的地，系统自动记录整个出行过程中驾驶员作出的动态响应行为，包括出行前的路径选择和出行途中的路径改换，行为数据存储于数据库。
全文摘要
本发明涉及一种驾驶员对图状路径信息板动态响应行为的计算机动态模拟方法，模拟系统为驾驶员提供较为逼真的快速路出行情境和满足人际功效学原理的人机交互界面，出行情境信息包括基本出行信息、交通流运行状况、GRIP信息、车辆方位信息、其它辅助信息。驾驶员出行过程中经历的交通状况和GRIP信息都有快速路网宏观交通流仿真数学模型仿真而得，使得驾驶员出行过程中经历的交通状况更真实，增强收集的驾驶员响应行为数据的可靠性。驾驶员在模拟出行中通过鼠标和键盘把所作的路径选择输入给计算机。模拟方法具有开发成本低、模拟试验简便且易操作、试验场地不受限制、试验过程快速、试验成本较低等特点，可以为驾驶员动态路径选择行为的分析与建模提供基础数据。
文档编号G06F17/50GK101556625SQ20091005146
公开日2009年10月14日 申请日期2009年5月19日 优先权日2009年5月19日
发明者干宏程, 晴 汪 申请人:上海理工大学