以前書いた Pythonでカオス・フラクタルを見よう! と似たような数理モデル一発ネタであります. Optimal Velocity Model 坂東昌子さんという素粒子畑の物理屋さんがつくった渋滞モデルです. これは自然渋滞を表現するモデルとして有名です. 「事故なんかなくても車の数がある程度ありゃあ渋滞は起きるよ!」ということです. このモデルの骨子は, 前の車との車間距離に応じて目標速度が変化することにあります. 車間距離が大きければ法定速度で走ろうとするし, 小さければ速度を下げようとします. これを微分方程式で表現すると以下のようなものになります: \frac{d}{dt}x_n = v_n, \hspace{1cm}\frac{d}{dt}v_n = a\left[V(x_{n+1} - x_n) - v_n\right] $c$ は適当な定数. $V(h)$ はこんな概形です
![Pythonで自然発生する交通渋滞を見よう! - Qiita](https://cdn-ak-scissors.b.st-hatena.com/image/square/3827cc8bffdd3f0932538703817afbd093362f39/height=288;version=1;width=512/https%3A%2F%2Fqiita-user-contents.imgix.net%2Fhttps%253A%252F%252Fcdn.qiita.com%252Fassets%252Fpublic%252Farticle-ogp-background-412672c5f0600ab9a64263b751f1bc81.png%3Fixlib%3Drb-4.0.0%26w%3D1200%26mark64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZ3PTk3MiZoPTM3OCZ0eHQ9UHl0aG9uJUUzJTgxJUE3JUU4JTg3JUFBJUU3JTg0JUI2JUU3JTk5JUJBJUU3JTk0JTlGJUUzJTgxJTk5JUUzJTgyJThCJUU0JUJBJUE0JUU5JTgwJTlBJUU2JUI4JThCJUU2JUJCJTlFJUUzJTgyJTkyJUU4JUE2JThCJUUzJTgyJTg4JUUzJTgxJTg2JUVGJUJDJTgxJnR4dC1hbGlnbj1sZWZ0JTJDdG9wJnR4dC1jb2xvcj0lMjMyMTIxMjEmdHh0LWZvbnQ9SGlyYWdpbm8lMjBTYW5zJTIwVzYmdHh0LXNpemU9NTYmcz1mYzdkMjQ4ZTJiYjNiYWQxY2I4OTNkNjVmYjM5OWU5Zg%26mark-x%3D142%26mark-y%3D57%26blend64%3DaHR0cHM6Ly9xaWl0YS11c2VyLWNvbnRlbnRzLmltZ2l4Lm5ldC9-dGV4dD9peGxpYj1yYi00LjAuMCZoPTc2Jnc9NzcwJnR4dD0lNDBqYWJiZXJ3b2NreTAxMzkmdHh0LWNvbG9yPSUyMzIxMjEyMSZ0eHQtZm9udD1IaXJhZ2lubyUyMFNhbnMlMjBXNiZ0eHQtc2l6ZT0zNiZ0eHQtYWxpZ249bGVmdCUyQ3RvcCZzPTBmNDFkNjNkZGM4ZmQ0MTUxN2EyNzVkOGQ0MTlmYzY4%26blend-x%3D142%26blend-y%3D486%26blend-mode%3Dnormal%26s%3D9f4ce3859445c50ca530d77281c70c3e)