高次元自律系におけるNeimark-Sacker 分岐の一計算法

自律系にみられるリミットサイクルにおいて，パラメータの変化により特性乗数が複素単位円を横切る場合，Neimark-Sacker分岐を生じる．本論文ではこの分岐パラメータ値を求める新しい計算方法を提案する．固定点条件，Poincaré断面条件に加えて，共役複素数の特性乗数を代入した特性方程式を実部と虚部に分離し，それぞれを独立した条件として採用する．これら条件式を連立させ，ニュートン法により解く．新たな独立変数として，共役複素な特性乗数の偏角を用いた．ニュートン法のヤコビ行列の各要素は，複素数演算，数値微分を用いることなく，変分方程式の解と行列操作によって求めることができるため，ニュートン法の収束性もよい．例として8次元自律系のリミットサイクルのNeimark-Sacker分岐の計算結果をあげる．