数理解析研究所講究録 1244 巻 2002 年 91-94

91
ODE の確率的切換えによって生ずるフフクタルな解集合と一般化次元
北大大学院工学研究科

郷原一寿 (Kazutoshi Gohara)

量子物理工学専攻

Department of Applied Physics,Hokkaido University,
Sapporo 060-8628, Japan
E-mail:gohara@nd-ap.eng.hokudai.ac.jp

1

問題設定

外界と相互作用するシステムを連続力学系 (ODE: 常微分方程式) の確率的切換えとしてモデル化することを考える。 以下の非自励力学系について考察する。
$\dot{x}=f(x, I(t))$$x,$ $I\in R^{n}$

ここで、$x,$の式は, 状態,(1).

は状態、 ペクトル場、 外部入力、 時間をそれぞれ表すものとする。 この時間発展が, 他のシステムから入力 によって影響を受けることを表現し$f,$ $I,t$$x$$I$ている。

2

切換え入力に対するダイナミクス

2.1

周期入力

周期入力 $I(t)=I(t+T)$ に対しては、角度変数連続力学系が定義できる。
$\dot{y}=fI(y)$$\theta\equiv\frac{2\pi}{T}t$(2),$y\equiv(x, \theta)\in R^{n}\cross S^{1}$

また、 $\theta=2\pi$ でポアンカレ断面$\Sigma$を導入することにより、以下の:$\mathcal{M}$.

を導入することにより、 以下の離散力学系が定義できる。
$x_{\tau+1}=g_{I}(x_{r})$$x_{r}\in R^{n}$

周期入力 I、ベクトル場書くことにする。

$f_{I}$、写像$g_{I}$:$\Sigma$,(3).の対応関係を強調するために模式的に以下のように$Iarrow f_{I}arrow g_{I}$.

(4)

92
2.2

切換え入力

次に、複数の入力が確率的に切換えられるダイナミクスを考察する. 入力は周期関数の-周期分と定義する。 例えば、周期関数をフーリエ級数で表せば、振幅 $A\in R^{m}$ と時間長 $T\in R^{1}$ をパラメータとして、 有限時間幅で時間的に変化する入力を、入力空間(5)$\mathrm{I}=\mathrm{I}(A,T)$の一点として、一般的に表すことができる。 この空間で $N$ 個の要素からなる入力集合.
を定義する。以下、集合の添え字は省略する。周期入力の場合と同様をに対して、 以下のようにベクトル場の集合 {f* }、写像の集合に、入力の集合$\{I_{*}$$=I_{i}(A:,T_{-})\}_{-=1}^{N}$$\{g:\}$$\cdot$$\{I_{i}\}$対応させる。
$\{I_{-}\}arrow\{f-\}arrow\{g_{i}\}$.

(6)

ここで、 写像の集合は Barnsley によって IFS と名付けられており、全ての写像 g- が縮小写像で、それらを等確率で切り換えた場合には、以下の式を満たす、吸引的で唯一の不変集合 $C$ に収束することが証明されている [1]。
$C=\cup=1Ng_{}(C)$.

(7)

従って、入力が等確率で切換わる場合には、方程式 (1) の解は円筒空間となる [2]。
変集合 を初期状態集合とする、 以下の不変軌道集合$\mathcal{M}$

で、$\Sigma$上の不$\Gamma(C)$$C$$\Gamma(C)=\bigcup_{=1}^{N}\gamma_{}(C)$.

(8)

に対する初期状態集合 $C$ から出発する軌道集合である。さらに gが縮小写像でない場合にも式 (7) (8) が満たされること、入力にノイズが重畳した場合のを包み込む閉包が存在するフラクタルな階層構造に対する効果、不変集合 $C$ およひことなどが数値的に示されている [2, 3, 4]。 また、 式 (7) (8) を検証するための実験が、電気回路 [7]、テニ
バネの振動リカレントニューラルネットワークモデル
スの運動 [8] などの異なる対象について行われており、 良い一致を見せている。
ここで、$\gamma_{i}(C)$は入力$I_{i}$$\Gamma(C)$$[6]_{\text{、}}\mathrm{L}\mathrm{C}\mathrm{R}$$[5]_{\text{、}}$

2.3

マルチフラクタル

式 (7) (8) のフラクタル集合を定量化することを考える。 : が縮小率 : で表せる相似
縮小 (similitude) な縮小写像で、かつ開集合条件 (open set condition) を満たし、切換え確率力 : であるとき、
$g$

$s$

$\mathrm{s}$

$p$,

(9)

$\tau(q)=(q-1)D_{q}$ ,(10)$\sum_{i=1}^{N}s_{-}^{-\tau(q)}p_{-}^{q}=1$$q\in R$,が成立する。 ここで、

93は一般化次元である。 また、ルジャンドル変換によってマルチフラクタルスペクトルが求まる $[9, 10]$問題の出発点は ODE なので、縮小率 を (1) 式から求めることを考える。例えば、

であり、$D_{q}$。

$f_{i}$$s_{i}$

が線形ならば も線形となり、 独 ( な線形空間での の縮小率は の固有値で表され、それによって、式 (9) を書き換え 6 ことが可能である。 ここでは、 が非線形でも適用可能な表式を導出したい。 円筒空間とポアンカレ断面 の体積変化率が等しいことおよび、 Liuville の定理より、.
(11),$\backslash \mathrm{I}$$g_{i}$$f_{i}$$g_{i}$$f_{i}$$\Sigma$$\mathcal{M}$$J(g:)=\mathrm{e}^{\int_{0}^{T}}.\cdot$が得られる。

は写像 のヤコビアンである。 この式は が未知でも、 のみによって、 のヤコビアンが求められることを意味している。 ここで、等方的な縮小性を仮定してみる。即ち、(12),$J(\mathrm{g}\text{ }\equiv J_{i}$$g_{i}$$g_{i}$$f.\cdot$$g_{i}$$s_{i}=J_{i}^{\frac{1}{n}}$である。 この仮定が成立する場合には、 式 (9) は、 以下のようになる。

.

(13)

$\sum_{i=1}^{N}\mathrm{e}^{-\int_{0}^{T}d\dot{\cdot}vf.dt}n..p_{i}^{q}=1\ovalbox{\tt\small REJECT}\tau$

この式は、 ODE の確率的切り換えによって構成されるフラクタルな解集合のマルチフラクタル表現を与えている。
ベクトル場 に条件を付加することで、 この式はさらに考察できるが、 ここでは例として以下の場合を考えて見る。
$f_{i}$.

$divf_{i}=\lambda=const$ $<0,$ $T_{i}=T=const.,\forall i$.

(14)

この条件では、$\ln\Sigma p_{i}^{q}N$$D_{q}= \frac{n}{q-1}\frac{i=1}{\lambda T}$となり、 一般化次元が陽に表現できる。

$\lambda Tarrow \mathrm{O}$,(15)の場合には、 開集合条件を満たさず、 成立しないが、 カオスを持つ Duffing 方程式を確率的に切換えた場合に生じるフラクタル集合についても、 良く近似していることが分かってきている [11]。 この式には外部入力の振(または、 マルチフラクタルスペクトが入っていない。 このことは、
幅パラメータ$A_{i}$$D_{q}$

ル)が時間長と切換え確率以外は外部入力の詳細に依らないことを示しており、興味深い。 この特徴は、 条件 (14) 以外のより一般的な条件でも成り立つ。

3

力学系集合
ODE の確率的切換えは、 もともと著者が脳機能を知るために考えた抽象的な数学モデルであるが、古典系と量子系の相互作用系、ブロッホ方程式の外部磁場による励起なども同様な枠組みに対応していることが分かってきた。 また、制御工学では離散時間と連続時間の混合システムはハイブリッドダイナミカルシステムと呼ばれ、研究が活発化してきて

94
いるが [12]、現在その理論的枠組みを模索している段階にあり、本研究は一つの重要な視点を与えていると考える。

最後に、複雑系研究と力学系集合について簡単に記す。 IFS によって構成される不変集
合およひそのマルチフラクタル表現については、 既に詳細な研究が進められている [10]。

これらは離散力学系集合に対する研究であるといえる。一方、連続力学系集合についてはこれまで殆ど研究がない。離散およひ連続力学系集合として、現象を整理して行くことは、現在複雑系で対象になっている、 または今後なり得る多体系全般を研究して行くための新しい重要な切り口であると考える。 ここでは力学系集合を用いて部分系の解析をしたが、 これを拡張して、部分から全体を構成する方法を力学系として確立することが今後の大きな課題である。

参考文献[1 M. Barnsley, fihctals Everywhere, Academic Press, Boston (1988).
[2] K.Gohara and AOkuyama, Nactals 7, 205 (1999).
[3] K.Gohara and AOkuyama, Ructals 7, 313 (1999).
[4] R.Wia and K.Gohara, Int. J. Bifurcation and Chaos 11, 755 (2001).
[5] S.Sato and K.Gohara, Int. J. Bifurcation and Chaos 11, 421 (2001).

[6] K.Gohara, H.Sakurai, and S.Sato, }$\mathrm{u}ctals8,67$ (2000).
$I$[7] J.Nishikawa and K.Gohara, Int. J. Bifurcation and Chaos, in press.

[8] Y.Yamamoto and K.Gohara, Hrrman Movement Science 19, 341 (2000).
[9]$\mathrm{T}.\mathrm{C}$. Halsey et al., Phys. Rev. A. 33,1141, 1986.

[10] K. J. Falconer, Technique in Rnctal Geometry, John Wiley&Sons (1997).
[11] K.Gohara and J.Nishikawa, submitted.
[12] AS. Matvee and AV. Savldn, Qualitative Theory of Hybrid Dynamical Systems,Birkhauser, Boston, 2000.