Aritalab:Lecture/NetworkBiology/Coupled Oscillator
From Metabolomics.JP
結合振動子
自然振動にばらつきのある素子が相互作用する系を考える。 同期のしやすさはネットワークの形状によるが、ここでは最も簡単な場合として2個の振動子を考える。
2つの振動子
振動子の振る舞いは以下の式で記述される。
は i に固有の周波数をあらわす。また の部分は振動子どうしの引き込み項に対応し、位相の大きさによって周波数に影響を与える。 周波数の差分 が のとき(つまり振動子2が1より前にいるとき)は振動子1が加速して振動子2が減速し、逆に差分が のときは振動子1が減速して振動子2が加速する。
k (>0) が十分に大きければ振動子は同期、つまり共通の振動数をとって が安定する。
もし であれば k = 0 でも同期するため、一般性を失わずに と仮定して、上式の差分を考えよう。
定常状態では だから
- 、つまり ()
値 付近の振動子の振る舞いをみる。
関数 は正のy切片を持ち の波形を描きながら において極小値 をとり、において正の値に戻る。よって k の値によって以下の場合がある。
- k が大きく極小値が負の場合: 解を2つ持つ。そのうち となるほうは振動子2が振動子1の前にある状態で、安定解となる。もう片方は不安定解になる。
- k が適切な値で極小値が0の場合: 解を1つ持つ。この値は不安定解になる。
- k が0に近く、極小値が正の場合: 解を持たない。つまり定常状態が存在しない。
結論として、 を満たす臨界値 が存在し、この値より k が大きい場合は安定解 が存在する。 それと等しいか、小さい場合は安定解が存在せず、二つの振動子が同期することはない。
N個の振動子
振動子がN個の場合は
結合する振動子どうしは それ以外は 0 とする。 全ての振動子が結合した完全グラフの場合は、 の極限において解けることが知られている(蔵本モデル)。
(ここから先、アイデアのみ) 振動子集団(平均場)の位相 を用いて
と簡略化することで解くことができる。