1、1 小区开放对道路通行的影响 摘要 交通在城市的可持续发展中起着关键性的作用 ，老式封闭小区阻碍城市交通体系的发展，因此，研究小区开放对城市道路通行的影响具有重要意义。 对于问题一 ，本文从周边路网、路段服务水平及交叉口处交通状况三个大的方面来建立评价指标体系，用以分析小区开放对周边道路通行的影响。其中，周边路网情形主要从车辆通行能力及延误时间两个方面来衡量，路段服务水平选取了道路上行驶车辆的排队长度及饱和度两个指标，而在交叉口处则选择了车头间距和车辆平均行进速度两个指标。基于以上评价体系，为了客观、定量评价小区开 放前后对道路通行的影响，我们选取小区开放前后一天中九个时段的指标值建立投影寻踪
2、评价模型进行评价，运用模拟退火算法求得最优投影方向为 (0.4935,0.1349,0.4898,0.4962,0.1956,0.4624)， 最后求得各样本的投影评价值 ， 并将该小区开放前后的评价值做比较得出相应的结论 。 对于问题二 ，为研究小区开放对周边道路通行的影响，考虑到不同道路结构下交通流的差异，本文将道路结构细分为 直行车道、 T 字型车道以及入闸式车道三种情形分别建立车辆通行模型。在直行车道上，我们利用微分方程推导出了车速、车辆位 置、道路负荷系数等车辆行进特征的模型，类似于直行车道，我们得出了 T 型车道和入闸式车道的相应指标的模型。进一步地，我们使用元胞自动机进行仿真模拟
3、，将上述三种道路结构下的车辆行进特征作为元胞演变的规则，仿真得到不同道路结构下车辆通行模型的动态变化图。 对于问题三 ， 我们首先选取了正方形小区和三角形小区两种类型，分别对两种类型小区进行道路结构的设计，其中对于正方形小区设置三种周边道路结构，对于三角形小区设置一种周边道路结构。然后我们考虑了小区开放前后对高峰期车流量、平常期车流量的影响，利用第一问的评价模型得出开放前 后各类型小区不同时期的评价值，在高峰时期和平常时期方形小区的综合评价值分别为 1.2543、 0.3367，接着利用第二问的模型得出开放前后元胞所表示的车流量变化的平面演化图，结合仿真结果和评价值进行结果分析，得出在不同时期
4、各类型小区开放前后对道路通行的影响。通过观察，我们发现模型中方形小区 T型通道的开通并没有缓解交通拥堵，于是提出此种交通结构符合 Breass悖论的猜想，并利用仿真的数据进一步对我们的猜想进行了论证。最后我们设置元胞自动机中的变道概率和转弯概率在 20%范围左右波动，进行敏感性分析，发现评价值的波动在 5%的范围内，得出模型稳定性强的结论。 对于 问题四 ，我们根据以上研究得出的结论，从改善交通拥堵、解决安全问题等方面对有关部门提出了相应的建议。 最后 ， 我们对模型进行了中肯的评价和适当的推广 。 关键词 ： 投影寻踪 模拟退火算法 微分方程模型 元胞自动机 Breass 检验 2 一 问题
5、的重述 1.1 引言 关于进一步加强城市规划建设管理工作的若干意见文件的出台，引起了 社会各界广泛的关注和热议 ，小区开放究竟是利大于弊还是弊大于利，每个人都有自己独特的见解。一方面封闭式小区， 堵塞了城市“毛细血管 ”， 增加了交通的压力，阻碍着城市可持续发展；另一方面 小区开放后 ，也会遇到一系列的问题，比如 业主权利侵犯的问题，安全问题等。那么小区开放对道路通行会产生怎样的影响呢？ 1.2 问题的提出 （ 1）通过选取 构建恰当的评价指标体系，评价小区开放对其周围道路所产生的影响。 （ 2）通过建立车辆通行的模型，以此为基础分析小区开放对道路通行的影响。 （ 3）小区开放产生的效果，与诸
6、多因素相关，通过考虑不同类型的小区，在你们构建的模型的基础上，对各类型小区开放前后对交通产生的影响进行定量分析。 （ 4）依据研究结果，基于交通通行的 角度考虑，对城市规划和交通管理部门，提出你们的合理化建议。 二 问题的分析 2.1 问题一 要求我们通过建立评价指标体系，以此为基础来分析小区开放对其周围道路通行所产生的影响。小区开放会增加交通流，其对交通的影响是多方面多层次的，在建立评价指标体系时，采用定性定量分析相结合的原则，通过查阅文献和利用相关知识，我们选取了评价指标，同时对评价指标进行量化，构建评价指标体系，然后结合所选取的评价指标利用投影寻踪法进行综合评价。 2.2 问题二 要求我
7、们通过建立车辆通行的模型，以此讨论小区开放 对周围道路通行所产生的影响 。 首先我们把局部道路类型分为不同的情况进行考虑 ，即直行车道、 T 字型车道以及入闸式通道， 然后在道路类型的基础上研究正常路段以及交叉口路段的车辆通行 ， 最后在正常路段以及交叉口路段 ， 分别研究小区开放前后车辆位置 、 车速以及道路荷载量的变化 ，来分析小区开放对道路通行的影响。 2.3 问题三 通过 考虑不同类型的小区，利用所构建的模型，对小区开放前后对交通产生的影响进行定量分析。我们首先对小区类型进行划分，然后考虑小区周边道路不同的开路类型，由于不同时间段内车流量会存在不同，此时我们模拟每种类型下的小区开放前后
8、的车辆通 行状况，以此判断小区开放前后对道路通行的影响。 2.4 问题四 基于 交通通行的角度进行考虑，我们在前三问求解的基础上，通过对小区开放对其周围道路通行所产生的影响的研究，据此我们从不同角度向有关部门提出我们的建议。 三 问题的假设 3.1 假设小区周围 交通状况 处于理想状态 ， 没有交通意外事故的发生 。 3.2假设开放小区后 ， 其内部交通不会发生堵塞 。 3.3 假设车辆驶入公路时对于车道的选择是随机的 。 3.4假设小区周围公路均为双向四车道的主干道，小区内开放通道为次主干道 3 四 符号的说明 符号 说明 2PF 信号影响修正系数 PR 车辆到达率 车辆影响系数 iz 投影
9、特征值 R 局部密度的窗口半径 五 模型的建立与求解 5.1 问题一模型的建立与求解 要求我们通过建立评价指标体系，以此为基础来分析小区开放对其周围道路通行所产生的影响。小区开放会增加交通流，其对交通的 影响是多方面多层次的，在建立评价指标体系时，我们采用定性定量分析相结合的原则，通过查阅文献和利用相关知识，我们选取了评价指标，然后结合所选取的评价指标利用投影寻踪法进行综合评价。 5.1.1 评价指标的选取 为了使评价更具有科学性 、 全面性以及客观性 ，我们在选取交通影响度的评价指标时，采用定量指标与定性指标相结合，同时使指标之间具有可比性以及可行性的原则，从三个方面进行选取。 小区开放在一
10、定程度上会增加交 通流量 1，其作用于周边路网，进而对周边路段的服务水平以及交叉口路段产生影响，因此我们主要从这 3个方面进行 指标的选取。 结合我国城市发展现状 ， 同时综合考虑影响交通状况的主要因素 ， 我们最终选取了饱和度 、排队长度 、 通行能力 、 延误 、 车头间距以及 车流平均 速度 6个指标进行分析。 表 1 评价指标的选取 序号 指标 1 饱和度 2 排队长度 3 通行能力 4 延误 5 车头间距 6 平均速度 5.1.2 评价指标的量化计算 （ 1）饱和度 饱和度是反映路段服务水平的主要指标 以及衡量道路通行的重要参数。若饱和度大于 1，则表明此时交通处于拥堵状态；反之若小
11、于 1，则出现局部时间以及空间的过饱和现象。其计算公式如下： 其中 iv 表示实际交通量 ， 表示阻碍系数 ， ic 表示设计最大车流量 iiivx c4 （ 2）排队长度 排队长度主要是用来衡量交叉口状况，车辆在交叉口时，通常来说需要减速以及停车。当前的排队长度是分析时间段内滞留车辆的平均数和当前来车行程排队的车辆数目之和。其计算公式如下 : 22111111 m i n 1.0 ,LPLLPLLLvgRCsPFvgRCsgvCCQ P FgXC 其中 Lv 表示车道到达流量 ， Ls 表示车道饱和流率 ， C 表示信号周期 ， gC 表示绿信比 ， 2PF 表示信号影响修正系数 ， PR
12、表示车辆到达率 。 （ 3）通行能力 通行能力是指在单位时段内车辆通过断面的交 通流量 2， 其大小直接影响交通疏散及道路所承受压力的能力 。 我们首先计算建议值 得出理论上的通行能力 ， 然后 求出可能的通行能力，最后依据路段服务水平求解设计通行能力。其计算公式如下： 3600=/tppNhN N cN N v c 理 论理 论 其中 N理 论 /pcn h 表示理论通行能力 ， th s 表示 平均车头时距 ， pN / .lnpcn km表示可能通行能力 表示 车 辆 影响系数 ,c 表示交叉口影响系数 ， /vc 表示修正系数 。 修正系数在很大程度上有路段服务水平确定 ， 我们查阅相
13、关文献 ，得出建议值如下表： 表 2 道路修正系数 等级 快速路 主干道 次干道 支路 修正系数 0.75 0.8 0.85 0.9 （ 4）延误 延误是指车辆在行驶中 ，由于交通控制设施失误、受除自身以外其他车辆的影响等的阻碍而 造成 的时间损 失 ，主要有停车延误 、引道延误以及控制延误 。延误 可以用来反映 司机 的状态、油耗以及行驶时间的损失 。 我们 采用 交叉口 平均延误 进行求解，其计算公式如下： 5 0.5 11 m in 1,ggtTTdtxT 其中 T 表示信号周期长度 ， gt 表示有效绿灯时间 ， x 表示饱和度 （ 5）车头间距 车头间距是位于同一车道上处于行驶状态下
14、 ，对前后邻近两辆车辆的车头之间一种距离的度量。我们根 据车头时距来计算车头间距，其计算公式如下： 3.6ts hhV其中 shm 表示车头间距 ， ths表示车头时距 ， /V km h 表示车辆行驶速度 （ 6）车辆平均速度 车辆平均速度表示单位时间内车辆的行驶快慢程度和运动方向，其计算公式如下： sv t 5.1.3 评价指标体系的构建 小区开放 使交通量增加 ， 作用 于周边路网 ， 从而使周围路段的服务水平发生变化以及影响交叉口的交通状 况 3。三者之间存在关联性，进而各自所选取的指标之间也相互影响、相互制约，共同作用形成一个有机的整体，因此我们构建的小区交通影响度评价指标体系如下：
15、 图 1 小区交通影响指标体系 描述周边路网服务水平的指标有车流量 、车辆的行驶速度、 停车延误等，一方面由于在预测情况下，这些指标在获得准确的数据存在一定的难度，同时这些指标与通行能力、延误存在重复性，因此选择通行能力、延误来描述周边路网的影响。饱和度以及小区交通影响指标路段服务水平 饱和度 排队长度 周边路网 通行能力 延误时间 交叉口 车头间距 平均速度 6 排队长度能综合反映交通的拥堵情况，是 衡量服务水平的重要参数，因此选择它们来描述路段服务水平。研究新增加 的交通量 对 交叉口 产生 的影响， 实质 上 归结于交叉口服务水平 的改变 ， 而平均速度特别是高峰时段的车速和车头间距 ，
16、 最能反映交叉口的交通状况 ， 因此选择车头间距 、 平均速度来描述交叉口的交通 。 5.1.4 基于投影寻踪的评价模型的建立 我们打算通过比较小区开放前后交通状况的改变来判断小区开放的影响，由于我们选择了 6个评价指标，这些指标对小区开放前后的交通影响存在很大的差异，这样在进行评价时，容易出现仅根据个别指标，使评价带有主观性。为了更客观、科学的进行评价， 综合考虑各个指标的情况，我们选择投影寻踪进行评价。 投影寻踪法 3，通过某种组合使数据样本投影到低维子空间上 ,并通过极化投影指标 ,寻找出能反映原数据结构或特征的投影，即寻找出使投影指标函数达到最优的投影值 ,然后根据该投影值对样本集进行
17、相应的评价。这样得到的结果在一定程度上解决了仅根据单项指标评价而导致结果不相容的问题，从而可以提高综合评价各层次的分辨力和评价模型的精度。 投影寻踪综合评价模型建立与求解的具体步骤如下： Step1： 数据的处理 （ 1） 指标的正负向性分析 我们选取了 5个评价指标，在建立评价模型 之前，我们首先要对这些指标的正负向性进行分析，即分析该指标值的大小对道路通行的影响。 表 3各个指标的正负向性分析 指标 正负向性 饱和度 - 排队长度 - 通行能力 - 延误 - 车头间距 + 平均速度 + 定义 ： 正向性 ：指标值越大越好； 负向性 ：指标值越小越好。 （ 2） 归一化处理 由于所选取指标的
18、单位不统一，为了消除各评价指标量纲的影响，我们对正负向性不同的指标分别进行归一化处理。 对于正向性指标，归一化公式如下： m inm ax m injjjjjxxx xx 对于负向性指标，归一化 公式如下： m axm ax m injjjjjxxx xx 其中 max jx 为第 j 列指标的最大值， minjx 为第 j 列指标的最小值。 7 Step2:线性投影 从不同的方向或角度去观察数据，寻找最能充分挖掘数据特征的最优投影方向。 设 12 ma a a a 为 m 维单位向量 ， 即为各指标的投影方向向量 ， 则第 i 个样本在一维线性空间的投影特征值 iz 的表达式为 ： 1mi
19、j ijiz a x2101jmjjaa Step3:构造投影指标函数 在综合投影值时 ,要求投影值 1,.,6iZi 的散布特征满足局部投影点尽可能密集 ,最好凝聚 成若干个点团 ,而在整体上投影点团之间尽可能散开 2。为此 ,投影指标函数可构造为 : max aaQ a S D 其中 aS 为投影值 1,.,6iZi 的标准差 , aD 为投影值 1,.,6iZi 的局部密度 ,即 : 616711/1a i aia ij ij ij ijijS Z Z nD R r R r 式中 ,Z 为序列 1,.,iZ i n 的均值 ;R 为求局部密度的窗口半径 ,它的选取既要使包含在窗口内的投影
20、点的平均个数不太少 ,避免滑动平均偏差太大 ,又不 能使它随着 n的增大而增加太快 ;距离ij i jr Z Z; Uh为单位 阶跃函数， 综上所述，得到非线性优化模型 : 616711121m a x/1m a x / 2=101aaa i aia ij ij ij ijijklk l k lmi j ijimjjjQ a S DS Z Z nD R r R rR r ms t r z zz a xaa 8 Step4：优化投影方向 当各指标值的样本集给定时，投影指标函数 Qa只随着投影方向的投影寻踪法变化而变化 。 不同的投影方向反映不同的数据结构特征 ， 最佳投影方向就是最大可能体现高维
21、数据某类特征结构的投影方向 。 因此 ， 可通过求解投影指标函数最大化问题来寻找最佳投影方向 目标函数： max aaQ a S D 约束条件： 21. . 1njjst a 这是一个以 | 1 ja j n为优化变量的复杂非线性优化问题，运用 Lingo 等软件难以求出其解，考虑到模拟退火算法具有收敛速度快等特点，我们使用模拟退 火算法 4来求解最优投影方向。 求解最优投影方向的模拟退火算法模型可以表示为下式： 001 ( . 2 / ( . 2) ) b1 : ( ) ( ) 1 _ 0e xp( _ / ) _ 0*y sqr t b sum bfor k n b k randde lt
22、 a epde lt a e T de lt a eT q T 为 随 机 15 的 (0,1) 序 列产 生 新 解Step5: 求解投影评价值 利用求解所得的最佳投影方向 a 计算样本的投 影指标值： ,1ni j i jjz a xStep6： 求解结果及分析 我们 利用 matlab软件进行编程求解 ， 求解结果如下 ： 最优投影方向： (0.4935,0.1349,0.4898,0.4962,0.1956,0.4624) 表 4 18组样本的投影评价值 样本编号 投影评价值 样本编号 投影评价值 1 2 3 4 5 6 7 8 9 2.0159 1.5533 0.5033 0.993
23、5 1.1340 0.9731 0.4521 1.1493 1.7676 10 11 12 13 14 15 16 17 18 2.2560 1.9795 0.8626 1.5732 1.7559 1.6726 1.0097 1.5958 2.1200 其中 ， 1-9 组为该小区开放前各时段的投影评价值， 10-18 组为该小区开放后各时段的投影评价值。根据以上结果，小区在高峰期时的评价值普遍偏低，在平常期的评价9 值普遍偏高；且小区在开放后的周边道路的交通状况评价值明显提高，说明拥堵状况得到了一定程度上的改善。 5.2 问题二模型的建立与求解 要求我们通过建立车辆通行的模型，以此讨论小区开
24、放 对周围道路通行所产 生的影响 。 首先我们基于元胞自动机模型 ，通过规则的考虑， 把局部道路类型分为 3 种 情况进行考虑 ，即直行车道、 T 字型车道以及入闸式通道， 然后从正常路段以及交叉口路段两方面考虑 ， 分别研究小区开放前后车辆位置 、 车速以及道路负荷量的变化 ，来分析小区开放对道路通行的 影响 5。 初步解决思路 设 置 元 胞 自 动机 初 始 状 态 值映 射 到 二 维 元 胞 自 动 机上 并 显 示 初 始 值设 置 演 化 规 则进 行 演 化 模 拟循 环 次数 达 到输 出 图 像及 各 项 指 标Y E SN O选 择 基 本 的 道 路 结 构直 行 道
25、路 T 型 交 叉 口 入 匝 式 路 口通 过 微 分 方 程 推 导车 辆 各 类 通 行 特 征车 辆 速 度 车 辆 位 置道 路 负 荷系 数制 定 规 则图 2 求解流程图 5.2.1 基于直行路段的车辆通行模型的建立 Step1:车道变换的考虑 联系实际可知，左右车道都存在一定的概率相互变道，但由于靠边的车道是转向车道，所以车道 2 1,3 4的概率 mp 大于车道 1 2,4 3的概率 np 图 3小区外围道路车道示意图 10 Step2:车辆 位置与车流量变化 某一时刻 t车流量的变化为： 10t t mq q q q 其中 mq 表示该时刻末将要驶出该段路的车量数目 于是
26、远离交叉口处主干道上 车流量的变化公式如下： ,1 0 1 , 1 1 , 2t m n tttdM q M p M p qdt 其中 ,jtM 表示 t 时刻车道 j 的车流量， mp , np 表示两个车道内 车 辆 变道的概率。 Step3:道路负荷系数 道路负荷系数的求解公式为 ： tt Mx C 其中 C 表示公路的最大荷载量。 Step4:基于交通流的车速 交通流是车辆在道路上连续行驶形成的车流，宏观上交通流可以等于连续可压缩的流体。我们将物理学中的动量定理引入到交通流的概念中 ： 00mP Mv P Mv （ 1） 其中 P表示当前的交通压力， mv 是出口处的速度， 0v 是当
27、前驶入道路车辆的速度。 在元胞自动机的仿真中，我们假定在某一特定的时间段（如 高峰期）内驶入车辆的车速不变，驶入的车道变化为随机，则 000, 0MvSS （ 2） 对于 （ 1）式 我们对 车辆离入口处的距离 S 做微分 0()P M v MuS S S （ 3） 此时我们再引入车流密度 mk ，为车流量与速度的比值， m Mk v（ 4） 将 （ 4） 式带入 （ 3）式 得 0 0m( v )1- mmvkP v v vv v vk S t S k S S （ 5） 此时我们引入交通流中的连 续性方程 6 (k v)+0mmktS 利用分部积分法 ()mmmkk vkS t t （ 6）
28、 将 式 （ 6）带入（ 5）式得 ， 11 0( v )1 mvkdv P vvdt k S k S S 由于模型中车流密度和车流压力是常量，即 mk0, 0PSt 则 dv vvdt S （ 7） 对公式 （ 7） 进行离散化处理，我们就可以求得每辆车的速度公式： 11 ( t) ( t) nnnnvt v v v tt S S t 5.2.2 基于 “ T”型交叉口 通道 的车辆通行模型的建立 小 区 建 筑 群 A 小 区 建 筑 群 B图 4“ T”型交叉口 通道示意图 Step1:车量位置与车流量变化 当交叉口的交通灯开始转化为绿灯时，此时路口前存在一个排队长度为 Q 的车流，在我
29、们的仿真模型中，考虑到车辆行驶 的 安全性， 因此我们假定 前一辆车启动后，后一辆车在延迟时间 T 后才启动。 因此建立的 车 流 量位置模型如下： 2 m a x2m a xm a x 0L - ( n - 1 ) ( + D ) , 0 t0.01 for i=l:iter a1=suiji(n); y1=Target(x,a1); delta_e=y1-y; if delta_e0 y=y1; p=a1; else if exp(delta_e/temperature)rand() y=y1; p=a1; end end end l=l+1; temperature=temperatur
30、e*0.99; end w(k)=y; 23 e(k,:)=p; end toc disp(max(w); %求得各样本投影值 r a=e(find(w=max(w),:); for i=1:12 r(i)=sum(x(i,:).*a); end function a=suiji(n) for k=1:n b(k)=rand(); end temp=sum(b.2); a=sqrt(b.2/temp); end function y=Target(x,a) m,n=size(x); for i=1:m s1=0; for j=1:n s1=s1+a(j)*x(i,j); end z(i)=s1
31、; end Sz=std(z); R=0.1*Sz; s3=0; for i=1:m for j=1:m r=abs(z(i)-z(j); t=R-r; if t=0 u=1; else u=0; end s3=s3+t*u; end end Dz=s3; y=Sz*Dz; end 24 附录二：元胞自动机仿真程序（ matlab） clear all;close all;clc % %定义 button plotbutton=uicontrol(style,pushbutton,. string,Run,. fontsize,12,. position,100,400,50,20,. cal
32、lback,run=1;); erasebutton=uicontrol(style,pushbutton,. string,Stop,. fontsize,12,. position,300,400,50,20,. callback,freeze=1;); number=uicontrol(style,text,. string,1,. fontsize,12,. position,20,400,50,20); z=zeros(38,4); cells=z;%元胞矩阵 p1=0.4;%变道概率 p2=0.4;%转弯概率 t1=0;%时间变量 t2=0;%时间变量 %初始状态 cells(1,
33、1)=0; cells(1,2)=1; cells(38,3)=0; cells(38,4)=1; run=0; freeze=0; while 1 if run= %下行 for j=1:37 if cells(j,1)=1 cells(j+1,2)=1; else cells(j+1,1)=1; cells(j+1,2)=0; end end 25 if cells(j,1)=0 cells(j+1,2)=0; else cells(j+1,1)=0; cells(j+1,2)=1; end end if cells(j,1)=1 cells(j+1,2)=1; end if cells(
34、j,1)=0 cells(j+1,2)=0; end end %上行 for j=38:-1:2 if cells(j,3)=1 cells(j-1,4)=1; else cells(j-1,3)=1; cells(j-1,4)=0; end end if cells(j,3)=0 cells(j-1,4)=0; else cells(j-1,3)=0; cells(j-1,4)=1; end end if cells(j,3)=1 26 cells(j-1,4)=1; end if cells(j,3)=0 cells(j-1,4)=0; end end %显示图像 A,B=size(cel
35、ls); Area(1:A,1:B,1)=zeros(A,B); Area(1:A,1:B,2)=zeros(A,B); Area(1:A,1:B,3)=zeros(A,B); for i=1:A for j=1:B if cells(i,j)=1 Area(i,j,:)=255,222,0; elseif cells(i,j)=0 Area(i,j,:)=0,90,171; end end end Area=uint8(Area); Area=imagesc(Area); axis equal; axis tight; %计步 stepnumber=1+str2num(get(number,string); set(number,string,num2str(stepnumber); end if freeze=1 run=0; freeze=0; end drawnow end
