ExcelでHaken-Kelso-Bunzモデルを。

HKBモデルというのがあって、メトロノームに合わせて左右の指を往復させるときのモデルらしいんだけれど、左右の指の動きの位相差をφとして、

dφ/dt = -a sinφ - 2b sin2φ

なるそうだ。ある程度速くなると急に指の運動パターンが別のものになるということらしい。いつものようにRunge-Kutta8次のDormand-Prince法で計算。

a=1として、上からb=0,0.2,0.5。bが大きくなると位相差が反転してしまうみたいだな。

速いと確かについていけなくなって指がずれていきそう。

参考にしたのはこの本。