カオスについて2
daruma3940.hatenablog.com
院試も終わったし前回の続きをかくのじぇ
前回2次元調和振動子は摂動を加えられなければトーラスになり、
摂動が弱い場合でも元のトーラスは変形して新しいトーラスができ、摂動が強くなるとトーラスが壊れてカオスが発生するらしいことについて書いたのじぇ。
ここでポアンカレマップについて説明するのじぇ。
ポアンカレマップとはトーラスの断面を切り取ったものなのじぇ。
2次元調和振動子の場合は変数はx,y,px,pyの4つであるけれど、
ハミルトニアンが一定という条件
から変数(例えばpy)を一つ消せ、
y=0の点となった時に点をプロットするということをすれば
変数はx,pxだけになるので2次元平面に点をプロットできるようになるというわけのじぇ。
↓こんな感じ
ここで非加積分系について考えようじぇ。
周期軌道から少しずらしたものがポアンカレマップのどこに来るのかを考えると
となるのじぇ。
このMは安定性行列と呼ばれ、保存系ではdet(M)=1となるのらしいのじぇ(よくわかってない)
これはリウビルの定理と呼ばれるものなのじぇ(リウビルの定理って正準変換のところで出てきたけど...)
とすると
となり、
↓
面積S=ab-bc=det(M)=1
になるのでdet(M)は面積が何倍されたかを表すのじぇ。
ここで行列Mの固有値について考えると
ここでもし固有値が実数で重解ではない場合、どちらか片方は固有値が1より大きく片方が固有値が1より小さくなるのじぇ
ということは片方は写像を繰り返すたびに指数関数的な拡大をし、
片方は写像を繰り返すたびに指数関数的に小さくなっていくのじぇ。
安定方向の直線上にない、つまり安定成分と不安定成分を持っているものは不安定性の指数関数的増加に引っ張られて不安定方向に流されていくのじぇ
もし固有値が複素数の場合、
λ1はλ2の複素共役かつλ1*λ2=|λ1|^2=1なので
になりこのような行列Mは
であり、これは回転行列なのじぇ。
これは指数関数的に拡大しないため安定なのじぇ。
もし重解を持った場合、
λ1=λ2=1であり、このような固有値を持つ行列は
であり
となるのじぇ。
図に書くとこんな感じ。
これも指数関数的に大きくなっていくわけではないので安定(中立安定)なのじぇ(この図あんまり理解できていない)
もし上のハミルトニアンでk=0のときこれは加積分系であり、周期軌道から少し初期値をずらしても別の周期軌道に乗るだけなのでこれは指数関数的に大きくなっていかず、中立安定的なのじぇ
しかしkが大きくなっていくとトーラスが壊れ、もしそれが複素数のほうへ壊れたのなら、安定固定点が発生し、実数のほうへ壊れたのなら不安定固定点が発生するのじぇ
このようすをkicked rotatormodelで描いてみるとこんな感じになったのじぇ
kicked rotatormodelとは
周期的に撃力を受ける回転子のことで....って書くのはめんどくさいので
http://a-phys.eng.osaka-cu.ac.jp/suri-g/phys8.html
このページを見てくれだじぇ
これがK(キックの強さ)=0の時
キックのじゃまを受けないのでひたすらθは大きくなっていくのじぇ
こんな感じ まさしく中立安定
K=0.01のとき
K=1.2の時
島みたいなのが現れているのじぇ。これが複素数の方向に壊れたやつなのじぇ。
島と島の間には実数の方向に壊れたやつが存在しているけどカオス的になってしまっているので構造はみてとれないのじぇ
ちなみにこの島を拡大すると島の中にも島が存在するのじぇ。これはいわゆるフラクタル的構造なのじぇ
さらにキックを強くするとこんな感じ
これらの図は範囲が違ったりするので注意なのじぇ
Kが大きくなると完全にカオス的になるのはすべて実数方向にずれるから??
安定固定点と不安定固定点が交互に出てくるのはなぜ???
さっきの
についてもやってみたけどそんなにきれいなのは見られなかったのじぇ
でもまあ一応張っておくか
なんかうまくいかないなーー
なんでだろうだじぇ?
う~~んう~~ん......
分かったのじぇ
ハミルトニアンが
で表されるので運動方程式は
これは線形なので非線形性がないためカオスは現れないのじぇ。
そしてk=0の時トーラスっぽいものが現れていないのは
これは単純な調和振動子なので一度ポアンカレマップに乗った後運動をして、もう一度ポアンカレマップに乗ったときに完全に同じ点に乗ってしまう、
つまりトーラス上の一部しか運動をしないためトーラスが現れていないのじぇ
それに変数変換をすれば、2つの2次元調和振動子に分離もできるのじぇ
A'は体格化されているのでx',y'のように変数変換をすれば2つの調和振動子に分けることができるのじぇ。
つまりこれはkの値がなんであれトーラスでカオス性は持たないのじぇ
でもkが0.6の時にはトーラスが出てきてないよね???これはなんで???
そしてトーラスの中心が2つに分かれているのはなんでなんだじぇ???
~~~追記~~~
k=0.6のときトーラスが現れていないのは新しい座標では振動数がそれぞれsqrt(1+k), sqrt(1-k) になっているのでk=0.6のとき周期の比が2になるので
片方が1周する間にもう片方がちょうど2周するので同じ点に帰ってきてしまい、トーラスが見えにくくなっているのじぇ
トーラスの中心は片方の振幅がゼロになる箇所に対応するのでどの方向の振幅がゼロになるかによって2種類あるので中心は2つあるのじぇ