Agentic discovery of interpretable car-following laws from multi-source trajectory data
Organizations: College of Transportation Engineering & Key Laboratory of Road and Traffic Engineering of the Ministry of Education, Tongji University, Shanghai, 201804, China
Abstract
Car-following (CF) models is a fundamental component of microscopic traffic simulation and underpin the analysis of collective traffic dynamics, while its modeling becomes more challenging due to diverse human driving behaviors and the increasing presence of automated vehicle. Analytical CF models offer interpretable representations of longitudinal vehicle interactions, but their prescribed functional structures may limit the representation of heterogeneous driving behavior. This study proposes an agent-centered, physics-aware framework for the empirical discovery of CF laws from multi-source trajectory data. The framework iteratively combines agent-guided hypothesis generation, informed by traffic flow knowledge and state-dependent reconstruction errors, with evolutionary symbolic search and physics-aware multi-objective evaluation. Using NGSIM US-101, highD, and exiD, the framework discovers the Adaptive Equilibrium Response Model (AERM), which couples state-dependent relative-speed response with spacing-dependent equilibrium-speed restoration. AERM admits analytical equilibrium speed-spacing and flow-density relationships; its stability conditions are derived for the zero-delay formulation, and platoon simulations examine disturbance propagation. Evaluated against seven classical CF models on Shanghai naturalistic driving data, NGSIM I-80, and Waymo trajectories, AERM attains the lowest spacing and speed reconstruction errors under trajectory-specific calibration, reducing spacing RMSE by 9.11%, 7.59%, and 21.73% relative to the strongest baselines, and remains competitive under dataset-level shared parameters. These results demonstrate the fitting flexibility and analytical tractability of the discovered CF law, providing an interpretable basis for microscopic traffic simulation.
Figures & tables
| Research direction | Representative studies | Main approach |
| I. Analytical Car-Following Models | ||
| Classical CF modeling | GHR, Gipps, Newell, OVM, FVDM, IDM | Explicit mathematical formulations based on predefined behavioral assumptions and traffic flow theory |
| Trajectory-based calibration | Kesting and Treiber (2008) | Numerical optimization of model parameters within prescribed analytical structures |
| II. Data-Driven and Symbolic Discovery of CF Models | ||
| Data-driven CF modeling | Locally weighted regression, LSTM | Learning nonlinear vehicle-response relationships directly from observed trajectories |
| Neural-symbolic discovery | SciNet-CFM | Combining neural representations and symbolic regression to discover interpretable CF relationships |
| Dataset | Traffic environment | Role |
|---|---|---|
| NGSIM US-101 | Freeway traffic | Discovery |
| highD | Highway traffic | Discovery |
| exiD | Highway entry and exit areas | Discovery |
| Shanghai NDS | Naturalistic driving | External |
| NGSIM I-80 | Freeway traffic | External |
| Waymo | Autonomous-vehicle-recorded traffic | External |
| Component | Configuration |
|---|---|
| Sampling interval | s |
| Subtree crossover probability | |
| Coefficient optimization | Bounded multistart Nelder–Mead |
| Acceleration clipping bounds | m/s 2 |
| Minimum safeguarded spacing | m |
| Directional-test tolerance | m/s 2 |
| Parameter | Configuration |
|---|---|
| Encoding | Real-valued |
| Population size | 64 |
| Maximum generations | 10,000 |
| Elite individuals | 4 |
| Selection | Tournament (size 3) |
| Crossover probability | 0.9 |
| Parameter | Unit | Physical interpretation |
|---|---|---|
| – | Relative-speed response intensity | |
| s -1 | Equilibrium-speed restoration rate | |
| m/s | Desired free-flow speed | |
| m | Standstill equilibrium spacing | |
| s | Desired time-headway scale |
| Configuration | Discovered CF equation |
|---|---|
| Full Framework (AERM) | |
| w/o Agent-Guided Pathway | |
| w/o Evolutionary Pathway | |
| w/o Physical Constraints |
| Configuration | Gap RMSE (m) | Speed RMSE (m/s) | Collision rate (%) |
|---|---|---|---|
| Full Framework | 1.2683 | 0.8454 | 0.0407 |
| w/o Agent-Guided Pathway | 4.4877 | 1.0800 | 0.0407 |
| w/o Evolutionary Pathway | 4.5314 | 1.0812 | 0.0407 |
| w/o Physical Constraints | 8.7521 | 1.4986 | 0.0407 |
| Shanghai NDS | NGSIM I-80 | Waymo | ||||
|---|---|---|---|---|---|---|
| Model | Gap | Speed | Gap | Speed | Gap | Speed |
| GHR | 10.0970 | 0.9948 | 5.6453 | 1.2050 | 5.5636 | 0.8506 |
| Helly | 9.3016 | 1.0749 | 5.3515 | 1.2242 | 5.9596 | 0.9303 |
| Gipps | 9.6389 | 1.0731 | 5.1731 | 1.1270 | 5.1695 | 0.8484 |
| Newell | 11.0970 | 1.7086 | 6.1859 | 1.3207 | 7.7446 | 1.4549 |
| OVM | 12.7097 | 1.9593 | 6.3365 | 1.4418 | 22.1557 | 2.4185 |
| Shanghai NDS | NGSIM I-80 | Waymo | ||||
|---|---|---|---|---|---|---|
| Model | Gap | Speed | Gap | Speed | Gap | Speed |
| GHR | 4.4980 | 0.7551 | 2.8763 | 1.0125 | 2.0739 | 0.5697 |
| Helly | 2.0142 | 0.5962 | 1.5750 | 0.8954 | 0.9933 | 0.4322 |
| Gipps | 1.7206 | 0.5289 | 1.3725 | 0.8658 | 0.6683 | 0.3551 |
| Newell | 2.0394 | 0.6002 | 1.7771 | 0.9102 | 1.0450 | 0.4665 |
| OVM | 2.1060 | 0.6162 | 1.8906 | 0.9504 | 1.3486 | 0.5488 |