多次元非線形マップ
📂動力学 多次元非線形マップ 非線形マップの定義 マップ f : R m → R m \mathbf{f} : \mathbb{R}^{m} \to \mathbb{R}^{m} f : R m → R m が線形 でない場合、非線形 と言われる。 ビルドアップ あるマップが線形であることを示すのは難しいが、非線形であることを示すのは簡単だし、線形問題は簡単だが、非線形問題は難しい。この宇宙のほとんど全てが非線形であり、非線形は難しいからこそ、人間は最初に非線形を線形に変換しようとする。
上の図では、示された曲線は確かに曲がっているが、非常に小さい範囲で見ると、ほぼ直線と同じであることが分かる。少しの誤差を無視すれば、十分に小さい範囲内で非線形はほぼ線形と同じだ。この方法を線形マップ に適用すれば、非線形マップも事実上線形マップと見なせるだろう。
テイラー近似 によると、f ( x + h ) ≈ f ( x ) + h f ′ ( x ) f(x+h) \approx f(x) + h f '(x) f ( x + h ) ≈ f ( x ) + h f ′ ( x ) だが、x x x の近傍、B ( x ; h ) B(x; h) B ( x ; h ) では、f ( x + h ) f(x+h) f ( x + h ) は h h h に関して一次式で表される。これを多変数関数 に拡張するには、ヤコビアン を導入すれば良い。このようなアイデアのために解析学の発展は絶対に必要だったと言える。
f \mathbf{f} f の固定点 p \mathbf{p} p の近傍 B ( p ; ∣ h ∣ ) B \left( \mathbf{p} ; | \mathbf{h} | \right) B ( p ; ∣ h ∣ ) にて、f \mathbf{f} f が線形性を持つだけの小さいベクトル h \mathbf{h} h が与えられていた場合、
f ( p + h ) ≈ f ( p ) + D f ( p ) ⋅ h ⟹ f ( p + h ) − f ( p ) ≈ D f ( p ) ⋅ h
\mathbf{f} ( \mathbf{p} + \mathbf{h} ) \approx \mathbf{f} ( \mathbf{p} ) + D \mathbf{f} ( \mathbb{ p} ) \cdot \mathbf{h}
\\ \implies \mathbf{f} ( \mathbf{p} + \mathbf{h} ) - \mathbf{f} ( \mathbf{p} )\approx D \mathbf{f} ( \mathbb{ p} ) \cdot \mathbf{h}
f ( p + h ) ≈ f ( p ) + D f ( p ) ⋅ h ⟹ f ( p + h ) − f ( p ) ≈ D f ( p ) ⋅ h
f \mathbf{f} f が線形性を持つならば、f ( p + h ) ≈ f ( p ) + f ( h ) \displaystyle \mathbf{f} ( \mathbf{p} + \mathbf{h} ) \approx \mathbf{f} ( \mathbf{p} ) + \mathbf{f} ( \mathbf{h} ) f ( p + h ) ≈ f ( p ) + f ( h ) であるから、
f ( h ) = D f ( p ) ⋅ h
\mathbf{f} ( \mathbf{h} ) = D \mathbf{f} ( \mathbb{ p} ) \cdot \mathbf{h}
f ( h ) = D f ( p ) ⋅ h
つまり、f \mathbf{f} f は p \mathbf{p} p の近傍の h ∈ R m \mathbf{h} \in \mathbb{R}^{m} h ∈ R m を D f ( p ) ⋅ h ∈ R m D \mathbf{f} ( \mathbb{ p} ) \cdot \mathbf{h} \in \mathbb{R}^{m} D f ( p ) ⋅ h ∈ R m へと送る線形マップ T D f ( p ) : R m → R m T_{D \mathbf{f} ( \mathbb{ p} )} : \mathbb{R}^{m} \to \mathbb{R}^{m} T D f ( p ) : R m → R m であることを意味している。
このように、線形マップ T A T_{A} T A に対応する行列 A A A のアイゲンバリューを使って動的システムについて話したように、非線形マップ f \mathbf{f} f の固定点 p \mathbf{p} p についても、同じ議論を
D f ( p ) = [ ∂ f 1 ∂ x 1 ( p ) ⋯ ∂ f 1 ∂ x n ( p ) ⋮ ⋱ ⋮ ∂ f m ∂ x 1 ( p ) ⋯ ∂ f m ∂ x n ( p ) ]
D \mathbf{f} ( \mathbb{ p} ) = \begin{bmatrix}
{{\partial f_{1} } \over {\partial x_{1} }} ( \mathbb{ p} ) & \cdots & {{\partial f_{1} } \over {\partial x_{n} }} ( \mathbb{ p} )
\\ \vdots & \ddots & \vdots
\\ {{\partial f_{m} } \over {\partial x_{1} }} ( \mathbb{ p} ) & \cdots & {{\partial f_{m} } \over {\partial x_{n} }} ( \mathbb{ p} )
\end{bmatrix}
D f ( p ) = ∂ x 1 ∂ f 1 ( p ) ⋮ ∂ x 1 ∂ f m ( p ) ⋯ ⋱ ⋯ ∂ x n ∂ f 1 ( p ) ⋮ ∂ x n ∂ f m ( p )
で続けていけば良い。
非線形マップにおけるハイパーボリックとサドルの定義 D f ( p ) D \mathbf{f} ( \mathbf{p}) D f ( p ) のアイゲンバリュー を λ 1 , ⋯ , λ m \lambda_{1} , \cdots , \lambda_{m} λ 1 , ⋯ , λ m としよう。
2. ∣ λ 1 ∣ ≠ 1 , ⋯ , ∣ λ m ∣ ≠ 1 | \lambda_{1} | \ne 1, \cdots , | \lambda_{m} | \ne 1 ∣ λ 1 ∣ = 1 , ⋯ , ∣ λ m ∣ = 1 ならば、p \mathbf{p} p は ハイパーボリック hyperbolic と言われる。
3. ハイパーボリックの p \mathbf{p} p に対して、{ ∣ λ i ∣ > 1 ∣ λ j ∣ < 1 \begin{cases} | \lambda_{i} | >1
\\ | \lambda_{j} | <1 \end{cases} { ∣ λ i ∣ > 1 ∣ λ j ∣ < 1 を満たす i , j i,j i , j が少なくとも一つ存在すれば、p \mathbf{p} p は サドル saddle と言われる。
サドルでない時の判定法 p \mathbf{p} p がサドルでない場合、次の定理に従って、シンクまたはソース かを判断できる。
[1]: ∣ λ 1 ∣ < 1 , ⋯ , ∣ λ m ∣ < 1 | \lambda_{1} | < 1, \cdots , | \lambda_{m} | < 1 ∣ λ 1 ∣ < 1 , ⋯ , ∣ λ m ∣ < 1 ならば、p \mathbf{p} p はシンクである。 [2]: ∣ λ 1 ∣ > 1 , ⋯ , ∣ λ m ∣ > 1 | \lambda_{1} | > 1, \cdots , | \lambda_{m} | > 1 ∣ λ 1 ∣ > 1 , ⋯ , ∣ λ m ∣ > 1 ならば、p \mathbf{p} p はソースである。