The science behind the simulator — why plain IDM misses the interesting
physics, and two Bayesian extensions that fix it:
MA-IDM (Gaussian-process driver
noise) and dynamic-regression IDM
(AR(p) driver noise).
The colour of the noise matters.
After fitting IDM to a real driver, the acceleration residual
\(\eta_n(t)\) is not white — its memory effect extends
over the past several seconds of driving. Feeding the right kind
of correlated noise back into a simulator is what turns a clean
deterministic trace into the persistent stop-and-go waves you see
in real highway data.
This page explains the two Bayesian models the simulator is built
on, and why they produce qualitatively different jam dynamics from
the white-noise baseline.
White — B-IDM (i.i.d.)AR(p) — DR-IDMGP — MA-IDM
Sample realisations of \(\eta(t)\) for each noise model — same
marginal variance, very different temporal structure.
01
Motivation — why stochastic car-following?
Deterministic car-following models are elegant, parsimonious, and capture
the first-order dynamics of traffic flow very well. But real trajectory
data — HighD, NGSIM, any modern drone survey — looks nothing like a
deterministic simulation. Two facts drive this work:
① Residuals are large and persistent
After fitting IDM to a driver's trajectory, the acceleration residual
\(\eta_n(t) = \ddot x_n^{\mathrm{obs}}(t) - f_{\mathrm{IDM}}(\cdot)\)
has a marginal standard deviation on the order of
0.2 m/s² and a memory horizon spanning several seconds:
the MA-IDM GP lengthscale fitted on HighD is
\(\ell\approx 1.4\) s, and the residual still carries useful
information from the past ~5 s of driving (MA-IDM
analysis) — pushed out to ~10 s by the AR(p)
calibration. It is emphatically not white noise.
② Noise colour changes the macroscopic flow
White-noise residuals produce high-frequency jitter with little
effect on flow. Correlated residuals with the same marginal
variance produce persistent stop-and-go waves and measurable
fundamental-diagram hysteresis. The
compare page shows this directly.
Parameters: desired speed \(v_0\), safe time headway \(T\), minimum gap
\(s_0\), max acceleration \(a\), comfortable braking \(b\), acceleration
exponent \(\delta\). The stochastic extensions all replace the
deterministic acceleration with
Key idea. Model the residual \(\eta_n(t)\) as a
zero-mean Gaussian process with a stationary covariance
kernel \(k(t,t')\). The GP endows the residual with memory:
how strongly \(\eta_n(t)\) depends on \(\eta_n(t-\tau)\) is
controlled entirely by the kernel.
The simulator exposes four kernels (RBF and Matérn 5/2, 3/2,
1/2), all with the same marginal variance \(\sigma_k^{2}\) and
lengthscale \(\ell\). Smoother kernels (RBF) produce smoother
acceleration noise; rougher kernels (Matérn 1/2, the
Ornstein–Uhlenbeck limit) produce visibly jittery traces with the
same long-run variance.
Simulation. Exact GP sampling is \(O(N^3)\) in the
trajectory length, so we use a random-Fourier-feature approximation
\(\eta(t)\approx \sigma\sqrt{2/M}\sum_{m=1}^{M}\cos(\omega_m t + b_m)\)
with \(\omega_m\) drawn from the kernel's spectral density — this
brings sampling down to \(O(M)\) per step while preserving the
target covariance.
04
Dynamic-regression IDM — AR(p) noise
AR(p) · DR-IDM
Calibrating Car-Following Models via Bayesian Dynamic Regression
The coefficients \(\{\rho_1,\ldots,\rho_p\}\) are inferred jointly
with the IDM parameters in a fully Bayesian framework. The paper
compares orders via short-horizon RMSE / CRPS on HighD and
recommends \(p\approx 4\!-\!6\) (using AR(5) for its
demonstrations); the simulator exposes \(p=1,\ldots,7\) and uses
the posterior-mean \(\boldsymbol\rho\) reported in Table 1 of
the paper (HighD, 5 Hz sampling).
Why AR? AR(p) is the discrete-time analogue of a
low-order linear SDE. It is cheap to simulate (\(O(p)\) per step,
no matrix factorisations), exactly stationary for any admissible
\(\boldsymbol\rho\), and the coefficients themselves are
interpretable — \(\rho_1\) is the lag-1 autocorrelation, and the
dominant root of the characteristic polynomial sets the memory
timescale.
On the benchmark data, AR(p) matches the GP's predictive
log-likelihood while being roughly an order of magnitude faster to
fit — attractive for real-time or online calibration.
05
Baseline — B-IDM (white noise)
White · B-IDM
I.i.d. Gaussian residual (naive baseline)
The classical Bayesian IDM likelihood.
What falls out of a naive least-squares or Gaussian-likelihood
calibration: residuals are assumed independent across time.
Useful as a baseline precisely because its failure modes — washed-out
jam waves, unrealistic short-timescale oscillation — are exactly what
the two correlated models above are designed to fix.
Zhang, C., & Sun, L. (2024).
Bayesian Calibration of the Intelligent Driver Model.
IEEE Transactions on Intelligent Transportation Systems.
doi:10.1109/TITS.2024.3354102
[arXiv:2210.03571]
Zhang, C., Wang, W., & Sun, L. (2024).
Calibrating Car-Following Models via Bayesian Dynamic Regression.
Transportation Research Part C: Emerging Technologies, 104719.
doi:10.1016/j.trc.2024.104719
[arXiv:2307.03340]
BibTeX
@article{zhang2024maidm,
title = {Bayesian Calibration of the Intelligent Driver Model},
author = {Zhang, Chengyuan and Sun, Lijun},
journal = {IEEE Transactions on Intelligent Transportation Systems},
year = {2024},
doi = {10.1109/TITS.2024.3354102}
}
@article{zhang2024dynamicidm,
title = {Calibrating Car-Following Models via Bayesian Dynamic Regression},
author = {Zhang, Chengyuan and Wang, Wenshuo and Sun, Lijun},
journal = {Transportation Research Part C: Emerging Technologies},
year = {2024},
pages = {104719},
doi = {10.1016/j.trc.2024.104719}
}