[论文解读] Conjectures about Traffic Light Queues

Steven R. Finch, Guy Louchard|arXiv (Cornell University)|Oct 9, 2018

Stochastic processes and statistical mechanics参考文献 3被引用 4

一句话总结

本文提出了离散时间红绿灯系统中最大队列长度的猜想渐近分布，该系统具有周期性的红绿相位，周期长度 ℓ ≥ 1。通过大规模行列式方程的数值求解与高精度计算下的代数识别，作者推导出最大队列大小 Mn 尾部衰减速率常数 χₗ(p) 的闭式表达式，将 ℓ=1 时已知结果扩展至 ℓ=2 和 ℓ=3 的情形。主要贡献在于提出了在周期性信号控制下最坏情况拥塞的极端值分布猜想。

ABSTRACT

In discrete time, $\ell$-blocks of red lights are separated by $\ell$-blocks of green lights. Cars arrive at random. The maximum line length of idle cars is fully understood for $\ell = 1$, but only partially for $2 \leq \ell \leq 3$.

研究动机与目标

• 理解具有 ℓ 段红灯与绿灯相位的周期性红绿灯系统中，最坏情况下的交通拥塞（即最大队列长度）。
• 将 ℓ=1 时已知的精确渐近结果推广至 ℓ≥2 的情形，其中目前尚无法进行严格处理。
• 猜想控制最大队列大小分布尾部衰减的速率常数 χₗ(p) 的闭式表达式。
• 基于高精度行列式方程的数值求解，开发一种计算框架，以推断 χₗ(p) 的代数公式。

提出的方法

• 定义随机过程 Sj = max(Sj−1 + Xj, 0)，其中 Xj 取决于交通灯周期模 2ℓ 的相位。
• 将车辆到达建模为独立同分布的 Bernoulli(p)，绿灯相位时以确定性 -1 速率离开，红灯相位时无离开。
• 构建 (k+1)×(k+1) 矩阵 Uk 和 Vk，以表示 ℓ 段红灯与 ℓ 段绿灯相位的转移动力学。
• 数值求解行列式方程 det[I - Uₖ^ℓ Vₖ^ℓ z] = 0，寻找最接近 1 的实数 z > 1，得到 zₖ。
• 定义 χₗ(p) = limₖ→∞ (zₖ - 1)/(p/q)^{2k}，从渐近行为中提取速率常数。
• 利用有理数 p 的高精度数值结果 χₗ(p)，通过多项式回归与根式识别方法推断代数表达式。

实验结果

研究问题

• RQ1 在周期长度为 2ℓ 且包含 ℓ 个红灯/绿灯相位的周期性红绿灯系统中，最大队列长度 Mn 的渐近分布是什么？
• RQ2 在尾部衰减 P(Mn ≤ log(q²/p²)n + h) ∼ exp[−χₗ(p)/(2ℓ)(q²/p²)^{-h}] 中，速率常数 χₗ(p) 如何依赖于到达概率 p 和周期长度 ℓ？
• RQ3 能否开发一种系统性的计算方法，以对 ℓ ≥ 2 的情形猜想 χₗ(p) 的闭式代数表达式？
• RQ4 在不同策略下（包括确定性周期信号与随机信号相位），最坏情况拥塞表现如何比较？

主要发现

• 对于 ℓ=2，猜想的尾部分布为 P(Mn ≤ log(q²/p²)n + h) ∼ exp[−χ₂(p)/4 ⋅ (q²/p²)^{-h}]，其中 χ₂(p) = (q−p)²/(4q⁶)[(1−8p²+16p³−8p⁴)+(q−p)√(1+4pq)]。
• 对于 ℓ=3，猜想的尾部分布为 P(Mn ≤ log(q²/p²)n + h) ∼ exp[−χ₃(p)/6 ⋅ (q²/p²)^{-h}]，其中 χ₃(p) 为包含根式与 p 的多项式的复杂代数表达式。
• 通过行列式方程的高精度数值计算得到 zₖ，成功识别出 χₗ(p) 的代数形式，经验证 p=1/3, 1/5, 1/17, 1/19 时可得到精确的根式表达式。
• n=10¹⁰ 时最大队列长度的经验直方图与理论预测高度吻合，证实了 ℓ=2 和 ℓ=3 的猜想分布。
• 期望最大队列长度 Eₗ(n,p) 在 ℓ=1 时最小，ℓ=2 和 ℓ=3 时仅略差，而随机信号情形（ℓ=0）表现显著更差。
• 该方法成功识别出 χₗ(p) 中的代数结构，包括当 x=1/p 为奇数时分母为 12(−1+x)^9 的有理形式，以及通过回归得到的多项式系数。