数理解析研究所講究録
第 806 巻 1992 年 157-179

157

ベクトル場のホモクリニック分岐と Lorenz-type attractor
京都大学理学部数学教室

国府寛司

目次
1. Introduction
2. Guckenheimer-Williams の geomehic Loren-z attractor と Lorenz 血 ifurcafion
3. Robinson による invariant foliation の構成
4. Keller-Robinson による Lorenz map の不変測度の存在とその性質
5. Rychlik のホモクリニック分岐定理
6. Robinson のホモクリニック分岐定理$b$\S 1. Introduction
Massachusetts 工科大学で気象学を研究していた E. N. Lorenz は 1963 年に発表した論文[Lorl] で, 大気循環による気候変動を説明するために温度勾配を持つ 3 次元流体層における熱対流現象をあらわす偏微分方程式を Galerkin 近似して得られた Saltzman のモデルからさらに reducfion を行って次のような常微分方程式系を導いた.

(L)$(= \frac{d}{d\iota})$$\{\begin{array}{l}\dot{x}=\sigma(y-x)\dot{y}=rx-y-xz-\dot{z}=-bz+\eta\end{array}$はこの方程式のパラメータで,は Prandfl 数,

は Rayleigh 数と呼ばれる
流体力学的定数を表わし, また は対流現象を考えている領域の幾何学的形状によって決まる定数である. 詳しくは丹羽 [Niw; 付録 2] を参照せよ. その後, Lorenz [Lor2] は
W. Malkus が実験室中での water wheel の運動方程式からこの方程式を直接導いたことを報告している. これについては Sparrow [Spa;Appendix Bl を参照せよ.
Lorenz がこの方程式で特にここで$\sigma,$$r,$ $b$$\sigma$$r$$b$$\sigma=10$,$r=28$ ,$b= \frac{8}{3}$としてコンピュータによる数値実験を試みたところ, 次の図のような解を発見した.
この図 1 は次のようにして得られた 原点 (0,0,0) の近く (例えば $(x,y,z)=(0,0.01,0)$ ):に初期値をとり, その点から方程式 (L) を数値積分によって解いて解 $(x(t),y(t),z(t))$ の描

く近似曲線を求め, その x-z 平面への射影を図示する. Lorenz の数値実験によればこの

158

図1

$xarrow$

曲線は時間刻みを $\Delta t=0.O1$ として約 6000 ステップまでの計算をしても, 単純な解に収束することなく極めて不規則な仕方で 2 つの平衡点のまわりをまわり続けるようにみえる.
コンピュータによる数値計算は誤差を必然的に含んでおりまたこの図形は初期値やパラメータの値を微小に変化させても依然として観察されることから, Lorenz はこのような不規則な解が方程式 (L) において何らかの意味で安定に存在しており, それが実は大気循環の不規則性, ひいては気象現象の予測の困難さを示していると考えたようである.
このような不規則性を持った解は今日ではカオス (chaos) $-$ 渾沌一現象を示すものとして数学的考察の対象となっているが, Lorenz の発見した解はそのような研究の端緒を開いたものの一つであり, 特にこの図にみられる奇妙な解は Lorenz attractor と呼ばれている.

Lorenz 方程式の簡単な性質(4) 解の大域的存在と一意性 :Lorenz 方程式 (L) は完備である. すなわち任意の初期値に対しそこを通る解が時間大域的かつ一意的に存在する. 詳しくは Coomes $[Coo;Lemma$3.1] を参照せよ.
i)$(\Pi)$

不変領域

: 関数

$V(x,y,z)=rx^{2}+oy^{2}+\sigma(z-2r)^{2}$ は方程式 (L) に沿って$\frac{dV}{dt}=-2\sigma(rx^{2}+y^{2}+bz^{2}-2brz)$

159
をみたす. そこで領域 D を$D= \{(x,y,z)|\frac{dV}{dt}=-2\sigma(rx^{2}\cdot+y^{2}+bz^{2}-2brz)\geq 0\}$と定めるとこれは当な正の数$\epsilon$$R^{3}$

内の有界閉領域であり,

$D$

上での関数

$V$

の最大値を $M$ とするとき適に対して$E=\{(x,y,z)|V(x,y,z)\leq M+\epsilon\}$

ととると,

$E$

の境界上では$\frac{dV}{dt}\leq 0$

上至るところ内向きである.

であるから方程式 (L) の定めるベクトル場は $E$ の境界

さらに $E$ の外の任意の点から出る軌道は有限時間のうちにE に到達する. 従って Lorenz 方程式にみられる漸近的挙動はすべて楕円体領域 E の内部に閉じ込められている.
負の発散 :Lorenz 方程式 (L) の定めるベクトル場の divergence は $-(\sigma+b+1)$ で至るところ負の一定値をとる. 従って (D) で述べたこととあわせると, 図 1 で観察され$(\nearrow\backslash )$る Lorenz attractor を含む方程式 (L) のいわゆる極大コンパクト不変集合は体積

$0$

の連結な

閉集合に含まれることがわかる.

(二) 対称性: 方程式 (L) は involution$\iota:(x,y,z)arrow(-x,-y,z)$で不変である. またこの対称性から, 特に z-軸は方程式 (L) に関して不変である.

:

平衡点とその分岐 原点 $O=(0,0,0)$ は Lorenz 方程式の平衡点であり, $r>1$ のときはさらに 2 つの平衡点 $E_{\pm}=(\pm\sqrt{b(r-1)},\pm\sqrt{b(r-1)},r-1)$ を持つ. これらおよび関連するLorenz 方程式の分岐について簡単にまとめておく. 以下ではパラメータの具体的な数値はとくに断らないかぎり $\sigma=10,$ $b=$ 男のときのものである.
$(’+\backslash )$[1] $0<r<1--$ 原点 $O$ が唯一の attractor である.

実際このときは Lyapunov 関数として$V(x,y,z)=x^{2}+\sigma(y^{2}+z^{2})$がとれる.
[21 r=l—— 原点 O が pitchfork 分岐を起こし,のときの原点での固有値は[3] $1<r<^{470}/19$

-

原点

$O$

$0,$

2 つの安定な平衡点 E\pm が出現する.

$-b,$ $-(1+\sigma)$ である.

は不安定化し saddle となる. その固有値は$\lambda_{1}=\frac{1}{2}\{-(\sigma+1)+\sqrt{(\sigma-1)^{2}+4\alpha}\},$ $\lambda_{2}=\frac{1}{2}\{-(\sigma+1)-\sqrt{(\sigma-1)^{2}+4or}\},$ $\lambda_{3}=-b$であり,$r>1+ \frac{b(\sigma+1+b)}{\sigma}=4.644$のとき$-\lambda_{2}>\lambda_{1}>-\lambda_{3}>0$,,と書くこともある.
となる. そこでこれらの固有値をはまたどの固有値も$r>$ 1.346... では複素共役な固有値を持ち,
一方平衡点$\lambda_{1}=\lambda_{u}$$E_{\pm}$$\lambda_{3}=\lambda_{s}$$\lambda_{2}=\lambda_{ss}$

こ

160

$r< \frac{\sigma(\sigma+b+3)}{\sigma-b-1}=\frac{470}{19}=24.74$

はこの範囲で安定である.
数値実験によれば 1<r<2406... ではこの平衡点 4 が大域的にも安定であって, すべての軌道はこれに漸近するようである. しかし r=13.926... で原点のホモクリニック軌
道が形成され, $r>$ 13.926... ではそこから不安定な周期軌道が分岐する. この周期軌道
は が大きくなるにつれて次第に平衡点に近づき, r=24.74... で盈につぶれてしまう.
周期軌道が十分に平衡点に近づくと, r<24.06... であってもほとんどの軌道は平衡点に収束する前に過渡的ではあるが極めて長い時間にわたって非周期的で複雑な動きをみせるようになる. このような現象は preturbulence (Kaplan-Yorke $[KaYo]$ ) または metastablechaos (Yorke-Yorke $[YoYo]$ ) と呼ばれている. さらに $r=24.06\ldots$ をこえると, 安定な平衡点以外に strange attractor が出現する. (McLaughlin-Martin[McMa] は strange attractor が出現す
るのは r=2474... の Hopf 分岐の後であると述べたが, それは誤りである. )の範囲では負の実部を持つので平衡点$E_{\pm}$

$r$

[4]

$E_{\pm}$$r=^{470_{19}}/$–

このとき平衡点$E_{\pm}$

は subcrifical Hopf 分岐を起こして不安定化し,時に不安定な周期軌道が消滅する. このときの$-(\sigma+b+1)$ ,$E_{\pm}$

それと同

の固有値は

$\pm i\sqrt{\frac{2\sigma(\sigma+1)}{\sigma-b-1}}$である. 詳しい計算については Marsden-McCracken $[MaMc]$ を参照せよ.
[5] $470_{19}/<r--$ この範囲の$r$

では $r<30$ 程度までは strange attractor が存在し, それ以外のattractor は数値実験では存在しないようである. 一方十分大きい については様々な
auractor や分岐現象が調べられている. 以下 Sparrow $[Spa;Chap.4]$ にしたがってその概略をまとめておく.
Franceschini [Fra] は 99.524 <r<100.795... のときにみられる非対称周期軌道の周期倍分岐について調べた.
Manneville-Pomeau $[MaPo1,2]$ は $145$ <r<166.07... のときにみられる対称周期軌道の周期倍分岐および $r=166$ .07... の近くでみられる intermittent chaos について調べた.
Robbins [Robb], Lorenz [Lor3] は r>214.364... のときにみられる対称周期軌道の周期倍分岐を調べた. また $r=\infty$ では Lorenz 方程式が積分可能であることに注目して, これらの分岐をからの摂動として調べた.
さらに Shimizu-Morioka $[ShMo]$ , Shimada-Nagashima $[ShNa]$ は $\sigma=16,$ $b=4$ の場合に同様の周期倍分岐を調べた.
$r$$\ldots$’$r=\infty$(へ) Lorenz lot と 1 次元写像 . 図 1 の Lorenz attractor 上の点から出る軌道 $(x(t),y(t),z(t))$) をプロットするとに対し, その -座標 $z(t)$ の 番目の極大値をとするとき, 点 (zn’) はすべてあ図 2 のような極めて 1 次元性の高い図形が得られる. このことは点 (zn’
る単一の 1 次元写像のグラフ上の点としてよく近似されることを示している. この図形
を Lorenz plot と呼ぶ. 図 1 の Lorenz attractor 上でこのような極大値を達成する点の集合を考えるとこれはすべて $-bz+\eta=0$ で与えられる 2 次元曲面 の上にあって, そこから出発して再び 上に戻る点を対応させる写像は非常に強い 1 次元性を持ち, その 1 次
元部分が図 2 のような写像に支配されることを示唆する.
この図 2 の写像は Lorenz 方程式の対称性を良く反映したものでないので, 図 3 のような断面に取りかえてそれについて上と同様の操作を行うと, 数値実験ではやはり図
$n$

$z$

$z_{n+1}$$z_{n}$$z_{n+1}$$Z_{n}$$S$$S$$\Sigma$161
4 のような 1 次元的図形が得られる. これを表現する 1 次元写像を Lorenz map と呼ぶ.
断面が対称性 $\iota:(x,y,z)arrow(-x,-y,z)$ で不変であるので, この Lorenz map は奇関数の対称性を持つ.
$\Sigma$$-\wedge+$$\backslash _{\vee}\backslash$$\in\alpha x$

図2

図3

162

$\uparrow$

$Y$

図4

$- 10$$-v^{\ulcorner}$

$0$

10

5

(ト) Lorenz map の分岐項で述べた Lorenz 方程式のパラメータ を変えたときの分岐は図 5 $(a)-(e)$ にみられるように Lorenz map に良く反映する. 実際, 図 5(a) は熊手が発生した状況を示し, 図 5(b) では原点において型分岐によって安定な平衡点
saddle-loop 分岐が起こり, そこから発生した不安定周期軌道が Lorenz map の不安定な不動点に対応する. それが図 5(c) の状況に近づくにつれて過渡的な非周期軌道がみられる
ようになることがわかる. さらに図 5(d) になると安定平衡点 以外に原点を含む閉区間が不変になり, その上では Lorenz map は拡大的になる. これが Lorenz attractor の Lorenz
map における対応物であり, さらに逆向きの Hopf 分岐を経て図 5 (e) に至ると, Lorenz
として残る.
attractor だけが唯一の$(’+\backslash )$$r$$E_{\pm}$$E_{\pm}$$a\mathfrak{n}raCtor$

ii)

本稿の目的
このノートは, 最近 Lorenz attractor の数学的理論に関して得られた著しい成果である
Rychlik [Rycl と Robinson $[Rob3,4]$ のホモクリニック分岐定理の紹介を中心にして, これまで Lorenz 方程式とその attractor に関して得られた成果を概観することを目的とする.
第 2 節では Guckenheimer-Williams $[GuWi]$ による geometric Lorenz attractor の理論と, その構造不安定性について述べる. 第 3 節では Robinson [Robl] による geomeffic brenz
amactor における strong stable foliation の摂動に対する存続性について紹介する. このstrong stable foliation により geometric Lorenz attractor の研究は 1 次元の Lorenz map の研究に帰着されるが, 第 4 節ではその Lorenz map の性質についての Keller [Kel] と Robinson[Rob2] の結果を中心に述べる. 第 5 節と第 6 節ではそれぞれ Rychlik [Rycl と Robinson[ROb3,4] の定理の概略を述べる.

163

図5

164
\S 2. $Guckenheimer- Williams$ の $geometric$ Lorenz attractor Lorenz bifurcation$\succeq$

Lorenz attractor の性質を数学的に解明するために Guckenheimer と Williams は論文[GuCl, [Wil], $[GuWi]$ においてもとの Lorenz 方程式で数値的に観察される strange attractorの性質を抽象して幾何学的に定式化された model flow を提出し, さらにその model flowの持つ auractor の構造についてほぼ完全な記述を与えた. これを geometric Lorenz modelと呼び, その attractor を Lorenz 方程式の attractor と区別して geometric Lorenz attractor と呼ぶ. ここではこの geometric Lorenz attractor の性質, 特にそれが Lorenz map と呼ばれる
1 次元写像によってほぼ分類され, そのことから geometric Lorenz attractor の構造不安定性が導かれることを示す. この最後の結果は Guckenheimer-Williams $[GuWi]$ によるもの
で, generic なベクトル場のほとんどすべての軌道は構造安定な attractor に漸近するであろうという Thom の予想 $([Tho;p.59])$ を否定的に解いたものである.
i)geometric Lorenz model と geometric Lorenz a 廿 ractor$R^{3}$

上の$C^{\infty}$

級ベクトル場 $X$ が次の (FP) をみたすものとする.

(FP) [fundamental properties of Lorenz attractor](a) $X(O)=0$ かつ DX(O) の固有値

(b)

$X$

$\lambda_{u},$$\lambda_{s},$$\lambda_{ss}$es involution

は実で$\lambda_{u}>0,\lambda_{s},\lambda_{ss}<0$をみたす ;$\iota:(x,y,z)arrow(-x,-y,z)$で不変である ;の固有空間である ;(C) -軸は$\lambda_{s}$

$z$

(d)

$|\lambda_{ss}|\geq|\lambda_{u}|+|\lambda_{s}|\geq 2|\lambda_{s}|$.

Lorenz 方程式はこれらの性質を確かにみたしており, さらに適当な定数 $C,$
まる長方形$\delta$について決$\Sigma=\{(x,y,z)||x|+|y|<C, z=\delta\}$が Lorenz 方程式の定める flow の局所断面を与える. 以下では座標を規格化して$\Sigma=I\cross I=[-1,1]^{2}$とする. 図 1 のような数値実験によればこの局所断面
ことはもっともらしい.

$\Sigma$について次の性質を仮定する(GLO) [confinemem of Poincar\’e map]

上の点は flow により再び 上に戻る. これにより定まる Poincar\’e 写像を
は $x=0$ でのみ不連続である.
であらわす.
(GL1) [existence of strong stable foliation]
上に で不変な曲線の族 $F$ が存在する. $W^{s}(O)\cap\Sigma$ は $F$ に属する. この $F$を 上の strong stable foliafion という.
$\Pi$$\Sigma$$\Sigma$

$\Pi$

$\Pi$

$\Sigma$$\Sigma$(GL2) [hyperbolicity]

Poincar\’e 写像$\Pi$

は $F$ の方向に縮小し,$F$に横断的な方向に拡大的である.

165
これは解析的には次のように定式化される.
(GLO)

Poincar\’e 写像 \Pi は\Sigma からそれ自身への C2 級写像で x=o 以外では l 対 l,でのみ不連続である.
上に $x=const$ . が で不変な曲線族 $F$ を定めるような適当な座標系 (X, ) が
(GL1)
とれる.
は $x=0$ に対応する. Poincar\’e 写像
は適当な関数 $f,$
を用いて
$x=0$
$\Sigma$

$\Pi$

$y$

$W^{s}(O)\cap\Sigma$$\Pi$$\Pi(x,y)=(f(x),g(x,y))$$g$

for $x\neq 0$と表現される. また対称性より$\Pi(-x,-y)=-\Pi(x,y)$がなりたつ.
は(GL2) 関数 $f,$$g$$f’(x)>^{\exists} \kappa>1,0<\frac{\partial g}{\phi}(x,y)<^{\exists}c<1$$f’(x)arrow\infty,$$\frac{\partial g}{\phi}(x,y)arrow 0$

for $x\neq 0$

as $xarrow 0$をみたす.
この flow の様子は図 6 のようなものである.

図6

逆にこれらの性質を持つ 内の flowの flow について正の方向に不変であり,$R^{3}$$\phi^{t}$が与えられれば, 図 6 に示される集合 はこ$\Lambda=\bigcap_{r\geq 0}\psi^{t}(V)$

$V$

166
は attractor を定める. 上の条件 $(FP),$ $(GLO)-(GL2)$ をみたす flow を \check geometric Lorenz modelを geometric Lorenz attractor と呼ぶ.と呼び, それから上式によって与えられる attractor
もとの Lorenz 方程式 (L) が条件 $(GLO)-(GL2)$ をみたすかどうかは数学的に確かめられていないので図 1 の attractor が geomettic Lorenz attractor であることはいまのところ証明されて$\Lambda$

いない.

geometric Lorenz attraCtor の逆極限による表現

ii)

Williams は geometric Lorenz model に対し 2 次元 branched surface およびその上の semiflow を定め, geometric Lorenz attractor はその逆極限として記述されることを示した. それによれば, geometric Lorenz model が与えられればその strong stable foliation の leaf に沿った射影による商空間をとることにより 2 次元 branched surface が得られ, 逆に 2 次元
branched surface とその上の semi-flow からその逆極限をとることにより geometric Lorenz
model と位相同値な flow が得られる. さらにこの逆極限で与えられる flow の性質からgeomehic Lorenz attractor の構造が数学的に調べられる. ここではその詳細に立ち入らな
い. 詳しくは Williams [Wil] を参照せよ.

Lorenz map の位相的性質について

iii)

ここでは幾何的 Lorenz attractor が対応する Lorenz map によってほぼ完全に分類され,そのことから幾何的 Lorenz attractor の構造不安定性が導かれることをを示す.

q 匣質 . geometric Lorenz attractor の定義 (GL1), (GL2) で与えられる 1次元写像 $f(x)$ を Lorenz map と呼ぶ. 第 1,4 節も参照せよ. この写像を規格化して区間$I=[-1,1]$ からそれ自身への写像と考えたとき, これは次の性質をみたす.
(イ) Lorenzm

$c$

(LM1)

f(X) は x=0 で不連続である.

(LM2)$\lim_{\uparrow X0}f(x)=+1$(LM3)$f\in C^{1}$(LM4)$\lim_{xarrow\pm 0}$on,$\lim_{\downarrow X0}f(x)=-1$$[-1,1]\backslash \{0\}$$\frac{f}{fx}(x)=+\infty$,$f(0)=0$.

, and $\exists\kappa>1s.t$ . $\frac{f}{fx}(x)>\kappa(\forall x\neq 0)$.

.

このうち, 条件 (LM3), (LM4) は第 i) 項の条件 (GL2) に対応している. Lorenz map はこれ以外にも対称性などの性質を持つが, 以下この項では主としてこれらの性質を持つ写像を扱う.
Lorenz map のカオスと拡大性 : 一般に区間 $I=[-1,1]$ 上の写像 $f$ が初期値に敏感に依存する (sensitive dependence on initial conditions) または chaofic であるとは, ある定数が存在して区間 の任意の点 と の任意の近傍 に対し$(\mathfrak{o})$$\delta>0$

(2.1)

$I$

$x$

$x$

$|f^{n}(x)-f^{n}(y)|>\delta$

$U$

167
となるような 内の点 と自然数 が必ず存在することをいう. すなわちどの点のどん
なに小さな近傍にも $f$ を何回か施すことによりやがてはある一定の距離以上離れてしまう点が見つかるということである. 写像 $f$ が拡大的 (expanding) であるとはこれよりさらに強く, 点 の近傍 のすべての点 について (2.1) が成り立つような自然数 が存在す
ることをいう.
Lorenz map については次が成り立つ.
$U$

$U$

$x$

命

V- 1

(1)
(2)

$n$

$y$

Lorenz map は拡大的であり, したがって特にカオス的である.
Lorenz map の周期点は区間 $I=[-1,1]$ の中で稠密である.

(ノ ) kneading invariant$\backslash$

$n$

$y$

:

$x\in I$に対し
$k_{0}(x)=\{\begin{array}{l}1(x>0)0(x=0)-1(x<0)\end{array}$とし, Lorenz map $f(x)$ に対し$k_{i}(x)=k_{0}(f^{i}(x))$

と定める.

$t$

についての形式べき級数 $k(x)$ を$k(x)= \sum_{i=0}^{\infty}k_{i}(x)t^{i}$と定義し, これを $x\in I$ の kneading series と呼ぶ.
順序, すなわち$\sum\theta_{j}t^{i}<\sum\theta_{i}’t^{i}\Leftrightarrow def\sum(\theta_{i}’-\theta_{i})t^{i}$

1 変数形式べき級数環 $Z[[t]]$ に辞書式の最初の零でない係数が正によって順序を定めると, 写像 $x\vdash k(x)$ は狭義単調増加である. 従って$f^{n}(x)=x\Leftrightarrow k(x)$: $periodic\Leftrightarrow k(x)=P_{n-1}(t)(1-t^{n})^{-1}$

がなりたつ. ここで $P_{n-1}(t)$ は $(n-1)$ 次以下の多項式である.

1 変数形式べき級数環 $Z[[t]]$ 上の距離を$d( \sum\theta_{i}t^{j},\sum\theta_{i’}t^{i})=\sum_{i=0}^{\infty}|\theta_{i}-\theta_{i}1^{2^{-i}}$で定義し,$k(x+):= \lim_{y^{\downarrow X}}k(y)$$k(x-);= \lim_{y^{\uparrow X}}k(y)$と定める. 特に $k_{\pm}(f)=k_{\pm};=k(0\pm)$ とするとき, 組う.
定理 2.2(Rand, Williams)2 つの Lorenz map $f,$

すなわち

$g$

$(k_{+},k_{-})$

を $f$ の kneading invanant といが同じ kneading invariant を持つ,

168$(k_{+}(f),k_{-}(f))=(k_{+}(g),k_{-}(g))$

であるならば,

$f$

と

$g$

は向きを保って位相共役である. また、$(k_{+}(f),k_{-}(f))=(-k_{-}(g),-k_{+}(g))$

であるならば、 $f$ と$g$

は向きを逆にする同相写像で位相共役である.

と Lorenz attractor の関係 . geometric Lorenz model においては Lorenz map
Lorenzが geometric Lorenz attractor の位相型を決定する. したがって上の定理よりつぎが得られ$(_{-}^{-})$$ma_{-}o$る.

をみたすならば,(LM3) の定数 がLorenz a 廿 ractor の位相型の完全不変量を与える.
系 2.3$\kappa$$\kappa>\sqrt{2}$$\{k_{+}, k_{-}, -k_{+}, -k_{-}\}$

はgeometric

実際, Williams [Wil] は $(LM3)$ の定数 が\kappa 1\kappa
をみたすならば Lorenz map は locally
eventually onto (LEO) であり, したがってその invariant set は 1 つの区間でその上でLorenz map は位相推移的であることを示した. よって上の定理 22 より Lorenz map が位相共役であることと geomehic Lorenz attractor が位相同値であることとが 1 対 1 に対応す$>\sqrt{2}$$\kappa$

る.

第 4 節も参照せよ.

特に対称性を考慮して奇関数のみ考えるならば,$\{k_{+}, k_{-}\}$

が完全不変量となる.

geometric Lorenz attractor の構造不安定性
Gucke 油 eimer-Williams-[GuWi] および Robinson [Robl] によれば, geometric Lorenz modelの条件は -open である. したがって geometric Lorenz model flow の 近傍に属する flowはやはり geometric Lorenz model であり, それには Lorenz map が対応する. この摂動された Lorenz map はもとの Lorenz map の近傍にある. ところが Williams [Wil], Rand [Ran]iv)$C^{2}$$C^{2}$$C^{0}$

により, つぎの結果が示されている.

ズヒ理 2.

$L$

2 つの Lorenz map $f,$$g$が$f(x)>g(x)$ for$\forall x\neq 0$を満たすならば, その kneading invariant は$k_{+}(f)>k_{+}(g)$,$k_{-}(f)>k_{-}(g)$

したがって特に任意の Lorenz map のを持つ Lorenz map が存在する.

である.

$C^{0}$

近傍には非可算無限個の異なる位相型

以上の結果をあわせると, つぎの geometric Lorenz model の構造不安定性定理が証明される.

定理 2.5 (Guckenheimer-Williams $\lceil GuWil$ ) geometric Lorenz model の 近傍に非可算無限個の異なる位相型を持っ geometric Lorenz attractor を持つ flow が存在する.
$C^{2}$

169\S 3. Robinson による invariant foliation の存続性
geometric Lorenz model を与える条件が open であることを示すためにはその摂動において strong stable foliation が存続し, Poincaree 写像から 1 次元の Lorenz map が得られることを示さなければならない. この strong stable foliation の存続性に関して Robinson [Robl]
はつぎの 2 つの定理を与えた.

geometric Lorenz model の原点における固有値が (FP-d) を満たすならば,
その geometric Lorenz model に 位相で近い flow は geometric Lorenz a 廿 ractor の近傍
で 級の strong stable foliation を持つ.
$\ovalbox{\tt\small REJECT}^{-\neq}31$$C^{2}$$C^{1}$

この定理の証明の概略は次のとおりである. まず geometric Lorenz model は自明な
strong stable foliation を持つ. その存在と微分可能性は flow の positive invariance と
hyperbolicity esnmate から保証される. 実際, strong stable 方向の縮小率はそれにになるにはさらに最大の拡
transversal な方向の縮小率よりも大きい. また foliation が大率と最大の縮小率の積も 1 より小さいことが必要である. これらの条件はすべてopen なので, 摂動された flow についてもやはり同様の条件が成り立ち, 従って strong
stable foliation が存続する.
$C^{1}$

笈月 B2 正方形領域$\Sigma$

上の写像

$\Pi$

が (GLO)-(GL2) に加えてつぎの性質を満たすと

する.

(a)

$| \frac{\partial g}{\partial x}(x,y)|\leq 1$(b)$g|_{\{x\}\cross J}$(C)$xarrow\pm 0$.

\rightarrow constam function$b^{\pm},$$\frac{\partial g}{\phi}(x,y)arrow 0$

as $xarrow\pm 0$ ;とするとき次はすべて一様有界である.

$\frac{\partial^{2}g}{\phi^{2}}(x,y),$$\frac{\partial^{2}g}{\partial x\phi}(x,y),$$\frac{\partial^{2}g}{\partial x^{2}}(x,y)\frac{\partial g}{b}(x,y)($$\frac{f}{fx}(x,y))^{-1}$$\frac{\partial^{2}f}{\partial x^{2}}(x,y)\frac{\partial g}{\phi}(x,y)(\frac{\#}{\partial x}(x,y))^{-1}$にで位相で近いものとと同じ性質を満たしこのとき写像はと定めると,
写像に 位相で近いものに対し $\Phi=J$
foliation の近くに 級の strong stable foliation を持つ.
$\Pi’$

$\Pi$

$\Pi$

。$\Pi’$$C^{1}$$C^{2}$$\Phi$

級写像$C^{2}$

$\Pi$

の

$J$

で恒等

strong stable$C^{1}$

この定理の証明にはいわゆるグラフ変換の方法を用いる. すなわち求める invariant
foliation を適当な関数空間の上のグラフ変換の不動点としてとらえ, 縮小写像定理を用いてその存在を証明する. その foliafion の微分可能性はグラフ変換の微分で定義される写像に対しファイバー縮小写像定理を適用することによりわかる. いずれにおいてもLipschitz 定数の評価にはもとの写像の hyperbolicity が重要な役割を果たす. ただし写像は不連続点を持つので求められた invariant foliation がそのまわりで 級になっていることを確かめることは注意する必要がある.
$\Pi$$C^{1}$

170
定理 3.1, 3.2 の証明の詳細は [Robl], $[HiPu]$ , [HPS] を参照せよ.

定理 32 から定理 3.1 を証明することができる. 実際 geometric Lorenz model を摂動した
flow に対しその Poincar\’e 写像を, 原点近傍において Cl 級微分同相写像で線形化することにより $\Phi=J\circ\Pi’$ の形の写像にとることができる. したがって定理 32 より, その写像は再び geomettic Lorenz model を与える. しかし, 摂動した flow が原点の近傍で 級の線形化が可能であるためにはもとの geometric Lorenz model の原点での固有値についてのある種の非退化条件が必要である. 実際の Lorenz 方程式の原点での固有値は$C^{1}$$\lambda_{u}=11.83,$ $\lambda_{s}=-2.67,$ $\lambda_{ss}=-22.83$

級にとらなければなであり, この値から 級の線形化のためには摂動が少なくともが 級でなければらないことがしたがう. さらに定理 32 を適用するためには写像であることが必要である. ただしこの条件は少し改良
ならず, そのためには摂動はすることもできる. 例えば Rychlik [$Ryc$ ;Appendix Il をみよ.
$C^{14}$$C^{1}$$\Pi’$$C^{20}$$C^{2}$

171\S 4. Keller-Robinson による Lorenz map の不変測度の存在とその性質
第 2 節で見たように geometric Lorenz model から得られる Poincar\’e 写像は$\Pi(x,y)=(f(x),g(x,y))$の形をしていた. これから導かれる 1 次元写像 $f(x)$ が geometric Lorenz model の性質をほぼ決定するが, とくにこの写像が Lebesgue 測度に同値な不変測度を持ち, それに関してこの写像が ergodic であるならば, もとの geometric Lorenz model は strange attractor を持つといえるであろう. この節では Lorenz map を含む 1 次元写像のあるクラスについて
このような性質が得られるための条件を与える.

Lorenz map に対する ergodic な不変測度の存在
Loren map を含む区分的に単調な 1 次元写像に対する不変測度の存在は Keller [Kel] によって与えられた. ここでは [Kel] の論文における定理 3.3 と定理 35 をまとめたつぎのi)形で述べておく.

笈月 AJ
(a)
(b)
(C)

区間

$I$

上の区分的に単調な写像 f がつぎの条件を満たすとする.

は各単調区間上で微分可能 ;の H\"older 連続関数である ;$/1f$ ’ は指数$f$$0<\theta\leq 1$$\exists\kappa>1\exists n\in Ns.t$.

$|(f^{n})’|\geq\kappa$.

このとき LebeSgue 測度に関して絶対連続な f.4–\変確率測度が存在し, つぎが成り立つ.

区間

$I$

の分割

$\{L_{k.l}\}_{\iota\leq k\leq r.1\leq\iota\leq}p_{k}$が取れて$f(L_{k.l})=L_{k.l+1}(mod p_{k})$

かつ

はすべての

$f^{p_{k}}|_{L_{k1}}$

$k,$

$l$

について $weak|y$ mixing,

したがって ergodic である.

実際, geometric Lorenz model から導かれる Lorenz map については上の 3 つの条件は満たされると考えてよい.
ii)

Lorenz map$0$) locally eventually onto property

上の Keller の定理の結論に加えてつぎの locally eventually onto property があれば得られた不変測度は f の定義域上の Lebesgue 測度と同値であり, その上で写像が位相推移的であることが保証される.
定義 4.2
部分区間

$J$

$f$

: $[-1,1]arrow[-1,1]$ が locallv eventuallv onto (LEO) とは区間 $[-1,1]$ の任意の開に対し自然数$k$が存在して $f^{k}(J)=[-1,1]$ となることをいう.

定湛 AU $f:[-1,1]arrow[-1,1]$ が Es-什 H\sim 府 擬 (topologically transitive) とは区間 $[-1,1]$ 内の点でその -極限集合が $[-1,1]$ 全体になるものが存在することをいう.
$\omega$

写像 $f$ が LEO ならば位相推移的である. より強く, 区間 $[-1,1]$ 内の residual set が存在してに属する点の -極限集合はすべて $[-1,1]$ 全体になる. 実際, 区間 $[-1,1]$ の有
$R$

$R$

$\omega$

172

限開被覆の列$O_{n}$

で任意の元 $U\in O_{n}$ はその長さが$R=$

寡$U\text{寡_{}\hslash}$火$/^{\iota_{n}}$以下であるものをとると,$f^{-k}(U)$で与えられる.

第 2 節でみたように Lorenz map $f(x)$ はつぎの性質を持っていた.
(LM1)$f(x)$ は区間 $[-1,1]$ 内のある点 $x=0$ で不連続である(LM2)$f(O-)=1,$ $f(O+)=-1,$ $f(-1)<0<f(1)$ ;(LM3)$f\in C^{1}$(LM4)$\lim_{xarrow\pm 0}f’(x)=+\infty$

on$[-1,1]\backslash \{O\}$か $\cdot 0\exists\kappa>1s.t.f’(x)>\kappa$ for;$\forall x\neq 0$;.

この写像 $f$ が LEO になるための条件として, Williams [Wil] は次を示した.

t 理生 4 条件 $(LM1)-(LM3)$ のもとで $f’(x)>\sqrt{2}$$(\forall x\neq C)$

ならば $f$ は LEO である.

この定理は写像 f が LEO になるための簡潔な条件を与えているが, 数値実験によればLorenz- 方程式については $f’(x)$ は端点の近くで 1 に近く, したがって上の定理の仮定が成り立つことは期待できない. そこで Robinson [Rob2] は $f$ が LEO になるためのより弱い条件を与えた.

笈月 A.5 条件 (LMI)-(LM3) およびっぎの条件のもとで $f$ は LEO である.
2 以上の自然数 について $f^{j}(-1)\leq 0$ , and $f^{j}(1)\geq 0$ for $0\leq j\leq k$ であり,$k$$f(B^{\pm})=0$

かつ $-1\leq B^{-}\leq 0\leq B^{+}\leq 1$$\exists\Lambda>0s.t$

かつ.

$M(f’)(x)\geq\Lambda$$\Lambda\lambda^{k}>2$

for

を満たす

$B^{-}\leq\forall x\leq B^{+}$,$B^{-},$ $B^{+}$

に対し

for

$f’(x)\geq\kappa>1$$\forall x\neq 0$が成り立つ. ここで $M(f’)(x)$ は $f’$ の から への平均値
$0$

$x$

$M(f’)(x)= \frac{1}{x}\int_{0^{X}}f’(t)dt=\frac{f(x)-f(0)}{x}$である.

ここで上の条件は必要である. 実際 (LMl)-(LM3) だけでは f が LEO にならない例がある.
詳しくは [Spa] または [Rob2; p.607] を参照せよ.

geometric Lorenz model ) transitive attractor
Robinson [Rob2] は第 1 項の Keller [Kel] の結果を用いて, geometric Lorenz model から導かれる Lorenz map は Lebesgue 測度に同値な不変測度を持ち, それに関して ergodic かつ位相推移的になることを示した.
iii)

カユ里\subset 6$\sigma$geometric Lorenz model の原点

$O$

での固有値

$\lambda_{u}$,$\lambda_{s}$,$\lambda_{ss}$が適当な正の数$\alpha$173
に対し$0<-\lambda_{s}<\lambda_{u}<-\lambda_{ss},$ $0<-\lambda_{ss}+\lambda_{s}-(1+\alpha)\lambda_{u},$ $0<\lambda_{u}+(1+\alpha)\lambda_{s}$を満たし, かつ の近傍で 線型化可能 (であるような固有値の非退化条件が満たされる) ならば, geometric Lorenz model の -近傍が存在してその中のすべての flow は
やはり geometric Lorenz model となり, それから導かれる 1 次元写像 $f$ は次の性質を満たす.
(a) $f$ は Lebesgue 測度 に絶対連続な不変測度 $\mu=hdm$ を持つ.
$C^{2}$$O$$C^{2}$$m$の support は $f$ の定義域全体に一致する.

(b)

密度関数

(c)

$f$

は不変測度 $\mu=hdm$ に関して ergodic である.

(d)

$f$

は LEO 従って位相推移的である.

$h\in L^{1}$

証明の概略はつぎのとおり. まず, 第 3 節の Robinson の結果から geometric Lorenz
model に 位相で近い flow はやはり strong stable foliation をもち, 第 2 節第 i) 項の条件
(GLO)-(GL2) を満たす. 従って, それから導かれる 1 次元写像も第 2 節第 iii) 項で挙げた
条件 (LMI)-(LM4) を満たす. 一方上の定理 45 の条件は自然数 を固定すれば open な条
件であるから 位相で近い flow に対しても成り立ち, その 1 次元写像は LEO になる.
Keller による第 1 項の定理 4.1 を適用するためにはその条件 (b) を確かめなければならない
が, これは $\exists\alpha>0s.t.f^{-1}$ は区間 $J_{-}=[f(-1),1],$ $J_{+}=[-1,f(1)]$ 上でそれぞれ Cl+\alpha 級関数に拡張することと同値であり, このことは strong stable foliation の微分可能性から導かれる.
これにより定理 4.1 が適用できて不変測度の存在が導かれるが, さらに Lorenz map はた
だ 1 つの不連続点しか持たないので, Li-Yorke $[LiYo]$ の定理により不変測度の support
は disjoint な区間であってそれぞれの区間は少なくとも 1 つの不連続点を含まねばならないことから, support は実は 1 つの区間をなすことがわかり, しかもその上の不変測度の一意性も導かれる. この support は, 写像 $f$ が LEO であることから $f$ の定義域全体となり, 従って不変測度は Lebesgue 測度と同値である. また $f$ が LEO であることから
位相推移的であることもしたがう. 詳しくは Robinson [Rob2] を参照せよ.
$C^{2}$$k$$C^{2}$

これに関連する結果として, Bunimovich と Sinai $[BuSi]$ は,

2 次元 Poincar\’e 写像がある種の双曲性条件を満たすならばいわゆる Bowen-Ruelle-Sinai measure が存在することを示した. しかしこれから上に述べた性質がただちに導かれるかどうかは明らかではない.
また Parry [Par] は区分線形写像で Lorenz- map と位相共役なものについては Lebesgue 測度に絶対連続な不変測度が存在することを示したが, それがもとの Lorenz map の不変測度になるためには位相共役写像に適当な微分可能性の条件が必要である.

174\S 5. Rychlik のホモクリニック分岐定理
1985 年頃 Rychlik は geomeffic Lorenz attractor がある種のホモクリニック軌道から分岐して出現することを証明したと発表した. この証明には誤りがあったが, 1989 年頃にはそれは修正された ([Rycl). ほぼ同じころ Robinson [ROb3] も Rychlik の結果にヒントを得て同様の分岐に関する結果を証明した. ホモクリニック軌道からの Lorenz-type attractor の分岐を考えるのは非常に自然なことに思われる (例えば Afraimovich-Bykov-Shilnikov[ABS] にもこのようなアイディアがのべられている) が, 実際にこれを証明しようとすると generic なホモクリニック軌道からは Lorenz attractor は分岐せず preturbulence においてみられるような双曲型不変集合しか出現しない. そこで attractor が分岐するようにもう つパラメータを導入して余次元 2 の分岐として考えようというのが Rychlik やRobinson の方法である. この第 2 の余次元のとりかたに複数の可能性があり, Rychlikは安定多様体の幾何学的配置に関する退化の条件を, また Robinson は原点における固有値についての共鳴の条件を課した. この節ではまず Rychlik によるホモクリニック分
岐定理を述べ, つぎの節で Robinson の定理を述べる.
$1$i)geometric Lorenz attractor from critically twisted double homoclinic loop
上のが次の条件 (FP) をみたすものとする.
級ベクトル場$R^{3}$$X$$C^{\infty}$(FP) [fundamental properties](a) $X(O)=0$ かつ $DX(O)$ の固有値
す;$\lambda_{u},$$\lambda_{s},$$\lambda_{ss}$は実で$\lambda_{u}>0,\lambda_{s},\lambda_{ss}<0$

をみた(b) $Xlh$ involution$\iota:(x,y,z)arrow(-x,-y,z)$で不変である ;$\Xi^{\infty}$

を$R^{3}$

上の-軸はの固有空間である ;

(C)

$z$

(d)

$|\lambda_{ss}|\geq|\lambda_{u}|+|\lambda_{s}|\geq 2|\lambda_{s}|$$C^{\infty}$$\lambda_{s}$.

級ベクトル場で, involution$\iota:(x,y,z)arrow(-x,-y,z)$を(との中で上の (FP) をみたすベクトル場
で不変なもの全体とし,同値なもの) 全体とする. 町に属するベクトル場に対しては次の 4 種類の不変多様体が存在する.
$\Xi_{0}^{\infty}$$X$$\Xi^{\infty}$$C^{\infty}$$X$

2 次元安定多様体 $W^{s}(O)$ および 1 次元不安定多様体 $W^{u}(O)$ ;1 次次元元不不変変多多様様体体$W^{ss}(O)=\{X\in R^{3}|-\iota\Vert_{\psi^{t^{t}}(x)||e^{-\lambda r}\leq Cfor\forall\leq^{\geq}0^{0}\}^{\};}}^{\phi(x)||e_{-\lambda_{s}^{ss,}}\leq Cfor\forall_{t^{t}}}$.

2 次元不変多様体 $W^{u}(O)=\{x\in R^{3}$
またこれから次のような 2 次元平面場も定義される.

安定多様体 $W^{s}(O)$ の接平面場
Pu=不変多様体の接平面場$P^{s}=$$W^{\iota}’(O)$175
不変多様体 $W^{us}(O)$ は一意的ではないが平面場は一意に存在する. 詳しくは HirschPugh-Shub 田 PS] を参照せよ.
を次の条件をみたすベクトル場の全体とする.
$P^{u}$$\Xi_{1}^{\infty}\subset\Xi_{0}^{\infty}$(HB1) [existence of homoclinic loop]$O$に関する double homoclinic loop が存在する, すなわち $W^{u}(O)\subseteq W^{s}(O)$ ;(HB2) [principal homoclinic loop]$W^{u}(O)\neq W^{ss}(O)$;(HB3) [critical twist condition]on$P^{u}=P^{s}$$W^{u}(O)\backslash \{O\}$;(HB4) [non-resonance condition]$|\lambda_{u}|\neq|\lambda_{s}|$(HB5)$\lambda_{u}<2|\lambda_{s}|$;.

ここで条件 (HB2) と (HB4) は (HBI) のもとで generic な条件であるが (HB1) と (HB3) ははそれぞれ余次元 1 の条件であり, 従っての余次元 2 の部分多様体になる. また条件 (HB5) は strange attractor が存在するために必要な open condition である. 以上の条
件をみたす double homoclinic loop をここでは crifically twisted double homoclinic loop と呼ぶ.
$\Xi_{0}^{\infty}$$\Xi_{1}^{\infty}$$R^{3}$

上の$C^{\infty}$

級ベクトル場の S-パラメータ族 $X_{\mu}(\mu\in U\subset R^{3})$ が次の意味で非退化であるとする.
(PD) [non-degeneracy of the parameter dependence]$\forall\mu\in U$,$X_{\mu}\in\Xi_{0}^{\infty}$か$\cdot\supset\exists\mu_{0}\in U$,$X_{\mu_{0}}\in\Xi_{1}^{\infty}$

and$(\mu\vdash\Rightarrow X_{\mu})$ns$\Xi_{1}^{\infty}$

このとき critically twisted double homoclinic loop から geometric Lorenz attractor が分岐する.

すなわち定理 5.1$\mu_{0}$$R^{3}$

上の$c^{\infty}$にいくらでも近い$\mu$が条件 (PD) をみたすならば,が geometric Lorenz a 廿 ractor を持つものが存在する.

級ベクトル場の S-パラメータ族
で$X_{\mu}$$X_{\mu}$ii)

Poincar\’e 写像の構成
Rychlik は Poincar\’e 写像を構成するために Roussarie の結果に基づいた smoothlinearization の方法を用いているが, それには原点における固有値についての非共鳴条件が必要である. ここではそれを用いないで, Shilnikov-Deng [Shi1-4, Den1-3] の方法を使って Poincar\’e 写像を表現する.
をとる. また
Poincar\’e 写像を考えるためにまず局所断面$\Sigma_{\pm}^{u}=\{x=\pm\delta\},$$v= \frac{-\lambda_{s}}{\lambda_{u}}$,$\eta=v+\omega$,$\Sigma^{s}=\{z=\delta\}$$\omega<\min(\frac{-\lambda_{s}}{\lambda_{u}},1,\frac{\lambda_{s}-\lambda_{ss}}{\lambda_{u}})$とおく. 定数\omega を上の条件を満たす範囲で十分大きくとるならば仮定 (FP) および (HB5)
より $\}_{2}’<v<1<\eta<2$ となる.

176
適当なパラメータの変換が存在して,
Poincar\’e 写像は次のような形で与えられる.
$\ovalbox{\tt\small REJECT}^{\wedge-}p52$$\mu\vdash\div\epsilon=(\epsilon_{1},\epsilon_{2})$$(x,y)\in\Sigma^{s}$に対し$\Pi;\Sigma^{s}arrow\Sigma^{s}$$\Pi(x,y)=(\begin{array}{l}\epsilon_{1}p_{2}(\epsilon)\end{array})+(\begin{array}{l}\epsilon_{2}r_{2}(\epsilon)\end{array})h(y,\epsilon)x^{\nu(\epsilon)}+k(x,y,\epsilon)x^{\eta(\epsilon)}$さらに$p_{2}(0)=0,$$r_{2}(0)=(\begin{array}{l}0*\end{array})h(0,0)=1$が成り立っ.

パラメータ空間の$0$の近傍をスケール集合$M_{\sigma}(c)= \{(\epsilon_{1},\epsilon_{2})|\epsilon_{2}=-\sigma\frac{c}{h(0,\epsilon)}|\epsilon_{1}|^{1- v(\epsilon)}$,$\sigma=sgn(\epsilon_{1})\}$

によって分割する.
命題 5.3
Poincar\’e 写像$\epsilon=(\epsilon_{1},\epsilon_{2})\in M_{\sigma}(c)$$\Pi$ならば変数 $(x,y)$ の適当な変換によって命題 5.2 のは正方形領域 $Q=(-1,1)^{2}$ 上の次のような形に帰着される.
$\Pi(x,y)=\Pi_{0}(x,y)+H(x,y)$

ここで、$\Pi_{0}$および $H$ は$\Pi_{0}(x,y)=\sigma sgn(x)\{(\begin{array}{l}1p_{2}(\epsilon)\end{array})-(\begin{array}{l}c0\end{array})h(sgn(x)y,\epsilon)|x|^{v(\epsilon)}\}$$H(x,y)=\sigma sgn(x)k(|x|,sgn(x)y,\epsilon)|x|^{\eta(\epsilon)}$であり,さらに適当な正の数\mbox{\boldmath $\zeta$} に対して次が成り立つ.
$p_{2}(\epsilon)=O(|\epsilon_{1}|^{\zeta})$

iii)

,

$h(0,\epsilon)=1$,$\frac{\partial h}{\phi}(y,\epsilon)=O(|\epsilon_{1}|^{\zeta})$,$\Vert\frac{\partial^{k+l}H}{\partial x^{k}\phi^{l}}\Vert\leq const|\epsilon_{1}|^{\zeta}|x|^{\eta- k}$

不変葉層の存在

以下では
像 $\Pi;Qarrow Q$ の不変葉層は$\epsilon=(\epsilon_{1},\epsilon_{2})\in M_{\sigma}(c)$$\Pi$が $/^{1}2<y_{c’}<v<1$ をみたす場合について考える. Poincar\’e 写で不変に保たれる

求めるために, 次のようなグラフ変換を考える.

$Q$上の曲線族として与えられる. それを$\Gamma(u)(x,y)=\frac{d(x,y)u(\Pi(x,y))-c(x,y)}{a(x,y)-b(x,y)u(\Pi(x,y))}$

177
参考文献
[ABS]

$[AfSi]$[Bun]$[BuSi]$

[Coo]
[Fra]
[Guc]

$[GuHo]$$[GuWi]$$[HiPu]$[HPS]$[KaYo]$[Kel]$[LiYo]$[Lorl][Lor2][Lor3]$[MaMc]$V. S. Afraimovich, V. V. Bykov, &L. P.Shilnikov, On structurally unstable attracting
limit sets of Lorenz attractor type, Trans. Moscow Math. Soc. (1983),2, pp.153-216.
V. S. Afraimovich&L. P. Shilnikov, Strange attractors and quasiattractors, inNonlinear Dynamics and Turbulence (Eds. G. I. Barenblatt, G. Iooss, &D. D.
Joseph), 1983, Pitman, pp. 1-34.
L. A. Bunimovich, Statistical properties of Lorenz attractors, in Nonlinear Dynamics
and Turbulence (Eds. G. I. Barenblatt, G. Iooss, &D. D. Joseph), 1983, Pitman, pp.

71-92.
L. A. Bunimovich&Ja. G. Sinai, Stochasticity of the attractor in the Lorenz model, in
Nonlinear Waves (Ed. A. V. Gaponov-Grekhov), 1979, Nauka, pp. 212-226.
B. A. Coomes, The Lorenz system does not have a polynomial flow, J. Diff. Eq.
(1989), 82, pp. 386-407.
V. Franceschini, A Feigenbaum sequence of bifurcations in the Lorenz model, J. Stat.
Phys. (1980), 22, pp. 397-406.
J. Guckenheimer, A strange strange attractor, in The Hopf Bifurcation and Its
Applications, Appl. Math. Sci., Vol. 19, (1976), Springer-Verlag, pp. 368-381.
J. Guckenheimer&P. J. Holmes, Nonlinear Oscillations, Dynamical Systems, andBifurcations of Vector Fields, Appl. Math. Sci., Vol. 42, (1983), Springer-Verlag.
J. Guckenheimer&R. F. Williams, Stmctural stability of Lorenz attractors, Publ.
Math. IHES (1979), 50, pp. 59-72.

M. W. Hirsch&C. C. Pugh, Stable manifolds and hyperbolic sets, in Proc. Symp.
Pure Math., Vol. 14, (1970), pp. 133- 163.
M. W. Hirsch, C. C. Pugh, &M. Shub, Invariant Manifolds, Lecture Notes in
Math., Vol. 583, (1977), Springer-Verlag.

J. L. Kaplan&J. A. Yorke, Preturbulence: aregime obserbed in a fluid flow model ofLorenz, Comm. Math. Phys. (1979), 67, pp. 93- 108.
G. Keller, Generalized bounded variation and applications to piecewise linear transformations, Z. Wahrscheinlichkeitstheorie verw. Gebiete (1985), 69, pp. 461-478.
T.-Y. Li&J. A. Yorke, Ergodic transformations from an interval into itself, Trans.
Amer. Math. Soc. (1978), 235, pp. 183- 192.
E. N. Lorenz, Deterministic nonperiodic flow, J. Atom. Sci. (1963), 20, pp. 130- 141.
E. N. Lorenz, On the prevalence of aperiodicity in simple systems, In Global Analysis(Eds. M. Grmela&J. E. Marsden), Lecture Notes in Math., Vol. 755, (1979),
Springer-Verlag, pp. 53-75.
E. N. Lorenz, Noisy periodicity and reverse bifurcation, in Nonlinear Dynamics (Eds.
R. H. G. Helleman), New York Acad. Sci., 357, (1980), pp. 1685-1699.
J. Marsden&M. McCracken, The Hopf Bifurcation and Its Applications, Appl.

178
Math. Sci., Vol. 19, (1976), Springer-Verlag.
[MaPol] P. Manneville&Y. Pomeau, Intermittency and the Lorenz model, Phys. Lett. (1979),$75A$ , p. 1.
$[MaPo2]$ P. Manneville&Y. Pomeau, Different ways to turbulence in dissipative dynamicalsystems, Physica 1D (1980), pp. 216 226.
$[McMa]$ J. B. McLaughlin&P. C. Martin, Transitions to turbulence in a statistically stressed
fluid system, Phys. Rev. (1975), A12, pp. 186-203.
丹羽敏雄, 微分方程式と力学系の理論入門, 遊星社, 1988.
[Niw]
-

[Par]

W. Parry, The Lorenz attractor and a related population model, in Ergodic Theory
Conference Oberwalfach (Eds. M. Denker&K. Jacobs), Lecture Notes in Math., Vol.
729, Springer-Verlag, pp. 169 -187.

[Ran]

D. Rand, The topological classification of Lorenz attractors, Math. Proc. Camb. Phil.
Soc. (1978), 83, pp.451 $- 460$ .

[Robb]

K. Robbins, Periodic solutions and bifurcation structure at high in the Lorenz system,
SIAM J. Appl. Math. (1979), 36, pp. 457-472.

[Robl]

C. Robinson, Differentiability of the stable foliation for the model Lorenz equations, In
Dynamical Systems and Turbulence (Eds. D. Rand&L.-S. Young), Lecture Notes in
Math. Vol. 898 (1981), Springer-Verlag, pp. 302-315.
C. Robinson, Transitivity and invariant measures for the geometric model of the Lorenz
equations, Ergod. Th. &Dynam. Sys. (1984), 4, pp.605- 611; Errata, ibid (1986), 6,
p.323.
C. Robinson, Homoclinic bifurcation to a transitive attractor of Lorenz type,
Nonlinearity (1989), 2, pp.495 518.
C. Robinson, Homoclinic bifurcation to a transitive attractor of Lorenz type, II, IMA
Preprint Series #719, 1990.
D. Ruelle, The Lorenz attractor and the problem of turbulence, in Turbulence andNavier-Stokes equations (Ed. R. Temam), Lecture Notes in Math., Vol. 565, (1976),
Springer-Verlag, pp. 146-155.
M. Rychlik, Lorenz attractors through Sil’nikov-type bifurcation, Part I, Ergod. Th. &Dynam. Sys. (1989), 10, pp. 793-821.
I. Shimada&T. Nagashima, The iterative transition phenomenon between periodic andturbulent states in a dissipative dynamuical system, Prog. Theor. Phys. (1979), 59, pp.

[Rob2][Rob3]$r$-[Rob4][Rue][Ryc]$[ShNa]$1605-1616.
$[ShMo]$T. Shimizu&N. Morioka, Chaos and limit cycles in the Lorenz equations, Phys. Lett.
(1978), $66A$ , p. 182.

$[SiVu]$Ja. G. Sinai&E. B. Vul, Hyperbolicity conditions for the Lorenz model, Physica $2D$
(1980), pp. 3-7.

[Spa]

C. Sparrow, The Lorenz Equations:Bifurcations, Chaos, and Strange Attractors,
Applied Mathematical Sciences Vol.41, (1982), Springer-Verlag.

179$[Tho;p.59]$$[YoYo]$E. D. Yorke&J. A. Yorke, Metastable chaos: transition to sustained chaotic behavior in
the Loren model, J. Stat. Phys. (1979), 21, pp. 263-277.

[Wil]

R. F. Williams, Structure of Lorenz attractors, Publ. Math. IHES (1979), 50, pp. 73-

100.
[Wonl][Won2]

S. Wong, Some metric properties of piecewise monotonic mappings of the unit interval,
Trans. Amer. Math. Soc. (1978), 246, pp. 493-500.
S. Wong, H\"older continuous derivatives and ergodic theory, London Math. Soc.
(1980), 22, pp. 506-520.