ヤコビ楕円関数 cnの2 乗 による Korteweg–De Vries 方程式のクノイド波 解(パラメータm の値は0.9 )。KdV方程式u t + u u x + δ 2 u x x x = 0 ( δ = 0.022 ) の数値解。初期条件はu ( x , 0) = cos(π x ) である。時間発展はザブスキー・クラスカル法によって行われた。初期の余弦波は孤立波の列へと発展する。 KdV方程式の2つのソリトン解が相互作用し(紫色)、それらが互いを通過する際に生じる位相シフトを強調している。赤と青の解は、一方のソリトンが存在しない状態での個々のソリトンの運動を示している[ 2 ] 。 数学 において、コルテヴェク・ド・フリース(KdV)方程式 は、浅い水面の波の数学的モデル として機能する偏微分方程式 (PDE)である。この方程式は積分可能なPDE の典型的な例として特に注目されており、PDEを一般的に扱いにくくする非線形性にもかかわらず、多数の明示的解(特にソリトン 解)や無限数の保存量などの典型的な挙動を示す。KdVは 逆散乱法 (ISM)によって解くことができる。 実際、クリフォード・ガードナー 、ジョン・M・グリーン、 マーティン・クラスカル 、ロバート・ミウラは KdV方程式を解くために古典的な逆散乱法を開発した。
KdV方程式は、 ジョセフ・ヴァレンティン・ブシネスク (1877年 、360ページの脚注)によって初めて導入され、 1895年にディーデリック・コルテウェグ とグスタフ・デ・フリース によって再発見され、最も単純な解である1ソリトン解が発見されました。方程式と解の挙動の理解は、1965年のノーマン・ザブスキー およびクラスカルのコンピュータシミュレーションと、1967年の逆散乱変換の開発によって大きく進歩しました。
1972年にT.Kawaharaは、KdV方程式の係数が非常に小さくなるかゼロになる場合の分散波を記述するKawahara方程式 として知られる5次のKdV型方程式を提案した。[ 6 ]
意味 KdV方程式は、(空間的に)1次元の非線形 分散 非散逸波をモデル化する 偏微分方程式 であり、次式に従う関数によって記述される: ϕ ( × 、 t ) {\displaystyle \phi (x,t)}
∂ t ϕ + ∂ × 3 ϕ − 6 ϕ ∂ × ϕ = 0 × ∈ R 、 t ≥ 0 、 {\displaystyle \partial _{t}\phi +\partial _{x}^{3}\phi -6\,\phi \,\partial _{x}\phi =0\,\quad x\in \mathbb {R} ,\;t\geq 0,} ここで は分散を表し、非線形要素は移流 項です。 ∂ × 3 ϕ {\displaystyle \partial _{x}^{3}\phi } ϕ ∂ × ϕ {\displaystyle \phi \partial _{x}\phi }
浅瀬の波をモデル化する場合、水面の高さの平衡高さからの変位です。 ϕ {\displaystyle \phi }
最後の項の前の定数は慣例的なものですが、大きな意味はありません。 、、およびを定数で乗算すると、3 つの項のいずれかの係数を任意のゼロ以外の定数に等しくすることができます。 6 {\displaystyle 6} t {\displaystyle t} × {\displaystyle x} ϕ {\displaystyle \phi }
ソリトンソリューション
1ソリトン解 で与えられる固定波形が、 位相速度 で右方向に移動する際にその形状を維持する解を考える。このような解は で与えられる。これをKdV方程式に代入すると、常微分方程式が得られる。 f ( X ) {\displaystyle f(X)} c {\displaystyle c} φ ( × 、 t ) = f ( × − c t − 1つの ) = f ( X ) {\displaystyle \varphi (x,t)=f(x-ct-a)=f(X)}
− c d f d X + d 3 f d X 3 − 6 f d f d X = 0 、 {\displaystyle -c{\frac {df}{dX}}+{\frac {d^{3}f}{dX^{3}}}-6f{\frac {df}{dX}}=0,} または、 について積分すると、 X {\displaystyle X}
− c f + d 2 f d X 2 − 3 f 2 = あ {\displaystyle -cf+{\frac {d^{2}f}{dX^{2}}}-3f^{2}=A} ここで、積分定数 である。上記の独立変数を仮想時間変数として解釈すると、これは ニュートンの立方体ポテンシャルにおける単位質量の粒子の 運動方程式を満たすことを意味する。 A {\displaystyle A} X {\displaystyle X} f {\displaystyle f}
V ( f ) = − ( f 3 + 1 2 c f 2 + A f ) {\displaystyle V(f)=-\left(f^{3}+{\frac {1}{2}}cf^{2}+Af\right)} 。もし
A = 0 , c > 0 {\displaystyle A=0,\,c>0} とすると、ポテンシャル関数はで極大値 をとります。 には、 の解が存在し、 は「仮想時間」 にこの点から始まり 、最終的に極小値 まで滑り落ち、その後反対側へ上昇して同じ高さに達し、その後方向を反転して の時点で再び極大値 に達します。つまり、は に近づくにつれて に近づきます。これが孤立波 解の特徴的な形状です。 V ( f ) {\displaystyle V(f)} f = 0 {\displaystyle f=0} f ( X ) {\displaystyle f(X)} − ∞ {\displaystyle -\infty } ∞ {\displaystyle \infty } f ( X ) {\displaystyle f(X)} 0 {\displaystyle 0} X → − ∞ {\displaystyle X\to -\infty }
より正確に言えば、解決策は
ϕ ( x , t ) = − 1 2 c sech 2 [ c 2 ( x − c t − a ) ] {\displaystyle \phi (x,t)=-{\frac {1}{2}}\,c\,\operatorname {sech} ^{2}\left[{{\sqrt {c}} \over 2}(x-c\,t-a)\right]} ここで は双曲正割 を表し、は任意の定数である。これは速度 で右方向に動くソリトン を表す。 sech {\displaystyle \operatorname {sech} } a {\displaystyle a} c {\displaystyle c}
N ソリトン解解を表す式が知られている。これは-ソリトン解であり、後期には別々の単一ソリトンに分解される。この解は、減少する正のパラメータの集合と非ゼロのパラメータの集合に依存する。この解は、 行列の成分が N {\displaystyle N} N {\displaystyle N} χ 1 > ⋯ > χ N > 0 {\displaystyle \chi _{1}>\cdots >\chi _{N}>0} β 1 , ⋯ , β N {\displaystyle \beta _{1},\cdots ,\beta _{N}} ϕ ( x , t ) = − 2 ∂ 2 ∂ x 2 l o g [ d e t A ( x , t ) ] {\displaystyle \phi (x,t)=-2{\frac {\partial ^{2}}{\partial x^{2}}}\mathrm {log} [\mathrm {det} A(x,t)]} A ( x , t ) {\displaystyle A(x,t)} A n m ( x , t ) = δ n m + β n e 8 χ n 3 t e − ( χ n + χ m ) x χ n + χ m . {\displaystyle A_{nm}(x,t)=\delta _{nm}+{\frac {\beta _{n}e^{8\chi _{n}^{3}t}e^{-(\chi _{n}+\chi _{m})x}}{\chi _{n}+\chi _{m}}}.}
これは逆散乱法を使用して導出されます。
運動の積分 KdV方程式には、時間とともに変化しない解の関数である運動の積分が無限に存在する。 これらは次のように明示的に表すことができる。 ϕ ( t ) {\displaystyle \phi (t)}
∫ − ∞ + ∞ P 2 n − 1 ( ϕ , ∂ x ϕ , ∂ x 2 ϕ , … ) d x {\displaystyle \int _{-\infty }^{+\infty }P_{2n-1}(\phi ,\,\partial _{x}\phi ,\,\partial _{x}^{2}\phi ,\,\ldots )\,{\text{d}}x\,} ここで多項式は再帰的に次のように定義される。 P n {\displaystyle P_{n}}
P 1 = ϕ , P n = − d P n − 1 d x + ∑ i = 1 n − 2 P i P n − 1 − i for n ≥ 2. {\displaystyle {\begin{aligned}P_{1}&=\phi ,\\P_{n}&=-{\frac {dP_{n-1}}{dx}}+\sum _{i=1}^{n-2}\,P_{i}\,P_{n-1-i}\quad {\text{ for }}n\geq 2.\end{aligned}}} 運動の最初の積分は次のとおりです。
質量∫ ϕ d x , {\displaystyle \int \phi \,\mathrm {d} x,} 勢い∫ ϕ 2 d x , {\displaystyle \int \phi ^{2}\,\mathrm {d} x,} エネルギー。∫ [ 2 ϕ 3 − ( ∂ x ϕ ) 2 ] d x {\displaystyle \int \left[2\phi ^{3}-\left(\partial _{x}\phi \right)^{2}\right]\,\mathrm {d} x} 奇数項のみが非自明な(つまりゼロでない)運動積分をもたらす。P 2 n + 1 {\displaystyle P_{2n+1}}
ラックスペア KdV方程式
∂ t ϕ = 6 ϕ ∂ x ϕ − ∂ x 3 ϕ {\displaystyle \partial _{t}\phi =6\,\phi \,\partial _{x}\phi -\partial _{x}^{3}\phi } はラックス方程式 として再定式化できる。
L t = [ L , A ] ≡ L A − A L {\displaystyle L_{t}=[L,A]\equiv LA-AL\,} シュトゥルム・リウヴィル演算子 を用いると: L {\displaystyle L}
L = − ∂ x 2 + ϕ , A = 4 ∂ x 3 − 6 ϕ ∂ x − 3 [ ∂ x , ϕ ] {\displaystyle {\begin{aligned}L&=-\partial _{x}^{2}+\phi ,\\A&=4\partial _{x}^{3}-6\phi \,\partial _{x}-3[\partial _{x},\phi ]\end{aligned}}} ここで はとなる交換子 である。 Lax 対はKdV 方程式の最初の積分の無限数を説明する。 [ ∂ x , ϕ ] {\displaystyle [\partial _{x},\phi ]} [ ∂ x , ϕ ] f = f ∂ x ϕ {\displaystyle [\partial _{x},\phi ]f=f\partial _{x}\phi }
実際、は(定数を無視すれば)時間に依存しないシュレーディンガー作用素 であり、そのポテンシャルはである。このラックス定式化により、固有値は実際には に依存しないことがわかる。L {\displaystyle L} ϕ ( x , t ) {\displaystyle \phi (x,t)} t {\displaystyle t}
ゼロ曲率表現 Lax接続 の成分をKdV方程式に 設定することは 、Lax接続のゼロ曲率方程式と同等である。 L x = ( 0 1 ϕ − λ 0 ) , L t = ( − ϕ x 2 ϕ + 4 λ 2 ϕ 2 − ϕ x x + 2 ϕ λ − 4 λ 2 ϕ x ) , {\displaystyle L_{x}={\begin{pmatrix}0&1\\\phi -\lambda &0\end{pmatrix}},L_{t}={\begin{pmatrix}-\phi _{x}&2\phi +4\lambda \\2\phi ^{2}-\phi _{xx}+2\phi \lambda -4\lambda ^{2}&\phi _{x}\end{pmatrix}},} ∂ t L x − ∂ x L t + [ L x , L t ] = 0. {\displaystyle \partial _{t}L_{x}-\partial _{x}L_{t}+[L_{x},L_{t}]=0.}
最小作用原理 Korteweg-De Vries 方程式
∂ t ϕ + 6 ϕ ∂ x ϕ + ∂ x 3 ϕ = 0 , {\displaystyle \partial _{t}\phi +6\phi \,\partial _{x}\phi +\partial _{x}^{3}\phi =0,} はラグランジュ密度から導かれる オイラー・ラグランジュの 運動方程式であり、L {\displaystyle {\mathcal {L}}\,}
L := 1 2 ∂ x ψ ∂ t ψ + ( ∂ x ψ ) 3 − 1 2 ( ∂ x 2 ψ ) 2 {\displaystyle {\mathcal {L}}:={\frac {1}{2}}\partial _{x}\psi \,\partial _{t}\psi +\left(\partial _{x}\psi \right)^{3}-{\frac {1}{2}}\left(\partial _{x}^{2}\psi \right)^{2}} 1
定義 ϕ {\displaystyle \phi }
ϕ := ∂ ψ ∂ x . {\displaystyle \phi :={\frac {\partial \psi }{\partial x}}.}
長期漸近解析 十分に速く減衰する滑らかな解は、最終的には右方向に移動するソリトンの有限の重ね合わせと、左方向に減衰する分散成分に分裂することが示される。これはZabuskyとKruskal (1965)によって初めて観察され、振動型 リーマン・ヒルベルト問題 に対する非線形最急降下法 解析を用いて厳密に証明できる。
歴史 KdV 方程式の歴史は、1834 年のジョン・スコット・ラッセルの実験から始まり、その後 1870 年頃に レイリー卿 とジョセフ・ブシネスク による理論的研究が行われ、最終的に 1895 年にコルテウェグとド・フリースによって研究されました。
KdV方程式はその後、ZabuskyとKruskal (1965) が 数値的にその解が長時間にわたって「ソリトン」、すなわち十分に離れた孤立波の集合に分解するように見えることを発見するまで、あまり研究されていませんでした。さらに、ソリトンは互いに通過しても形状にほとんど影響を及ぼさないようです(ただし、これにより位置が変化する可能性はあります)。彼らはまた、 KdV方程式がFPUT 系の連続極限 であることを示し、フェルミ、パスタ、ウラム、ツィンゴウによる以前の数値実験との関連を示しました。 逆散乱変換 を用いた解析解の開発は、1967年にガードナー、グリーン、Kruskal、三浦によって行われました。
KdV方程式は現在ではホイヘンスの原理 と密接に関係していると考えられています。
アプリケーションと接続 KdV方程式は物理問題といくつかの関連があります。連続体極限におけるフェルミ-パスタ-ウラム-ツィンゴウ問題 における弦の支配方程式であることに加えて、KdV方程式は、以下のような多くの物理的設定における長い1次元波の発展を近似的に記述します。
KdV 方程式は、非線形シュレーディンガー方程式 に適用されるような逆散乱変換 を使用して解くこともできます。
KdV方程式とグロス・ピタエフスキー方程式フォームの簡略化されたソリューションを考慮すると
ϕ ( x , t ) = ϕ ( x ± t ) {\displaystyle \phi (x,t)=\phi (x\pm t)} KdV方程式は次のようになる。
± ∂ x ϕ + ∂ x 3 ϕ + 6 ϕ ∂ x ϕ = 0 {\displaystyle \pm \partial _{x}\phi +\partial _{x}^{3}\phi +6\,\phi \,\partial _{x}\phi =0\,} または
± ∂ x ϕ + ∂ x ( ∂ x 2 ϕ + 3 ϕ 2 ) = 0 {\displaystyle \pm \partial _{x}\phi +\partial _{x}(\partial _{x}^{2}\phi +3\phi ^{2})=0\,} 積分し、積分定数がゼロとなる特殊なケースをとると、次のようになります。
− ∂ x 2 ϕ − 3 ϕ 2 = ± ϕ {\displaystyle -\partial _{x}^{2}\phi -3\phi ^{2}=\pm \phi \,} これは一般化定常グロス・ピタエフスキー方程式 (GPE) の特別な場合である。λ = 1 {\displaystyle \lambda =1}
− ∂ x 2 ϕ − 3 ϕ λ ϕ = ± ϕ {\displaystyle -\partial _{x}^{2}\phi -3\phi ^{\lambda }\phi =\pm \phi \,} したがって、一般化GPEの特定の解のクラス(真の1次元凝縮体の場合と、 1次元で3次元方程式を用いた場合)では、2つの方程式が1つになる。さらに、負の符号と実数の場合を例にとると、明るいソリトン を生み出す引力的な自己相互作用が得られる。 λ = 4 {\displaystyle \lambda =4} λ = 2 {\displaystyle \lambda =2} λ = 3 {\displaystyle \lambda =3} ϕ {\displaystyle \phi }
バリエーション KdV方程式には様々なバリエーションが研究されてきました。いくつかを次の表に示します。
名前 方程式 コルテヴェク・デ・フリース(KdV) ∂ t u + ∂ x 3 u + 6 u ∂ x u = 0 {\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u+6u\partial _{x}u=0} KdV(円筒形) ∂ t u + ∂ x 3 u − 6 u ∂ x u + 1 2 t u = 0 {\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u-6u\partial _{x}u+{\tfrac {1}{2t}}u=0} KdV(変形) ∂ t u + ∂ x ( ∂ x 2 u − 2 η u 3 − 3 u ( ∂ x u ) 2 2 ( η + u 2 ) ) = 0 {\displaystyle \displaystyle \partial _{t}u+\partial _{x}\left({\frac {\partial _{x}^{2}u-2\eta u^{3}-3u(\partial _{x}u)^{2}}{2(\eta +u^{2})}}\right)=0} KdV(一般化) ∂ t u + ∂ x 3 u = ∂ x 5 u {\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u=\partial _{x}^{5}u} KdV(一般化) ∂ t u + ∂ x 3 u + ∂ x f ( u ) = 0 {\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u+\partial _{x}f(u)=0} KdV(修正版) ∂ t u + ∂ x 3 u ± 6 u 2 ∂ x u = 0 {\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u\pm 6u^{2}\partial _{x}u=0} ガードナー方程式 ∂ t u + ∂ x 3 u − ( 6 ε 2 u 2 + 6 u ) ∂ x u = 0 {\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u-(6\varepsilon ^{2}u^{2}+6u)\partial _{x}u=0} KdV(修正版) ∂ t u + ∂ x 3 u − 1 8 ( ∂ x u ) 3 + ( ∂ x u ) ( A e a u + B + C e − a u ) = 0 {\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u-{\tfrac {1}{8}}(\partial _{x}u)^{3}+(\partial _{x}u)(Ae^{au}+B+Ce^{-au})=0} KdV(球状) ∂ t u + ∂ x 3 u − 6 u ∂ x u + 1 t u = 0 {\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u-6u\partial _{x}u+{\tfrac {1}{t}}u=0} 広田・薩摩方程式 { u t − 1 2 u x x x + 3 u u x − 3 ( v w ) x = 0 v t + v x x x − 3 u v x = 0 w t + w x x x − 3 u w x = 0 {\displaystyle \displaystyle {\begin{cases}u_{t}-{\frac {1}{2}}u_{xxx}+3uu_{x}-3(vw)_{x}=0\\v_{t}+v_{xxx}-3uv_{x}=0\\w_{t}+w_{xxx}-3uw_{x}=0\end{cases}}} KdV(スーパー) { ∂ t u = 6 u ∂ x u − ∂ x 3 u + 3 w ∂ x 2 w ∂ t w = 3 ( ∂ x u ) w + 6 u ∂ x w − 4 ∂ x 3 w {\displaystyle \displaystyle {\begin{cases}\partial _{t}u=6u\partial _{x}u-\partial _{x}^{3}u+3w\partial _{x}^{2}w\\\partial _{t}w=3(\partial _{x}u)w+6u\partial _{x}w-4\partial _{x}^{3}w\end{cases}}} KdV(移行型) ∂ t u + ∂ x 3 u − 6 f ( t ) u ∂ x u = 0 {\displaystyle \displaystyle \partial _{t}u+\partial _{x}^{3}u-6f(t)u\partial _{x}u=0} KdV(変数係数) ∂ t u + β t n ∂ x 3 u + α t n u ∂ x u = 0 {\displaystyle \displaystyle \partial _{t}u+\beta t^{n}\partial _{x}^{3}u+\alpha t^{n}u\partial _{x}u=0} KdV-バーガース方程式 ∂ t u + μ ∂ x 3 u + u ∂ x u − ν ∂ x 2 u = 0 {\displaystyle \displaystyle \partial _{t}u+\mu \partial _{x}^{3}u+u\partial _{x}u-\nu \partial _{x}^{2}u=0} KdV-バーガース-フィッシャー方程式(コチャック方程式) ∂ t u + μ ∂ x 3 u + ϵ u ∂ x u − ν ∂ x 2 u = r u ( 1 − u ) {\displaystyle \displaystyle \partial _{t}u+\mu \partial _{x}^{3}u+\epsilon u\partial _{x}u-\nu \partial _{x}^{2}u=ru(1-u)} 川原方程式 ∂ t u + α u ∂ x u + β ∂ x 3 u − γ ∂ x 5 u = 0 , {\displaystyle \partial _{t}u+\alpha u\partial _{x}u+\beta \partial _{x}^{3}u-\gamma \,\partial _{x}^{5}u=0,} 非均質KdV ∂ t u + α u + β ∂ x u + γ ∂ x 2 u = A i ( x ) , u ( x , 0 ) = f ( x ) {\displaystyle \partial _{t}u+\alpha u+\beta \partial _{x}u+\gamma \partial _{x}^{2}u=Ai(x),\quad u(x,0)=f(x)}
参照
注記
参考文献 ベレスト、ユーリ・Y.ロウツェンコ、イーゴリ M. (1997)。 「ミンコフスキー空間におけるホイヘンスの原理とコルテヴェーグ・ド・フリース方程式のソリトン解」。数理物理学におけるコミュニケーション 。190 ( 1 ): 113–132。arXiv : solv-int/ 9704012 。ビブコード : 1997CMaPh.190..113B 。土井 :10.1007/s002200050235 。ISSN 0010-3616 。 Boussinesq, J. (1877)、Essai sur la thetheorie des eaux courantes 、Memoires presentes par divers savants ` l'Acad.科学。研究所ナット。フランス、 XXIII、 1–680 ページ Chalub, Fabio ACC; Zubelli, Jorge P. (2006). 「双曲型作用素と可積分階層に対するホイヘンスの原理」 (PDF) . Physica D: 非線形現象 . 213 (2): 231– 245. Bibcode : 2006PhyD..213..231C . doi : 10.1016/j.physd.2005.11.008 . ダリゴル、オリヴィエ(2005年)『フローの世界 』オックスフォード、ニューヨーク:オックスフォード大学出版局。ISBN 978-0-19-856843-8 。 ドーソワ, ティエリー; ペラール, ミシェル (2006). 『ソリトンの物理学』 ケンブリッジ, イギリス; ニューヨーク: ケンブリッジ大学出版局. ISBN 0-521-85421-0 . OCLC 61757137 . ディンゲマンス, MW (1997).不均一な底部における水波伝播 . ニュージャージー州リバーエッジ: ワールドサイエンティフィック. ISBN 981-02-0427-2 。 ドゥナイスキー、マチェイ (2009)。ソリトン、インスタントン、ツイスター 。オックスフォード;ニューヨーク: OUP オックスフォード。ISBN 978-0-19-857063-9 . OCLC 320199531 . ガードナー, クリフォード・S.; グリーン, ジョン・M.; クラスカル, マーティン・D.; ミウラ, ロバート・M. (1967). 「コルテウェグ-デ・フリース方程式の解法」.フィジカル・レビュー・レターズ . 19 (19): 1095– 1097.書誌コード : 1967PhRvL..19.1095G . doi : 10.1103/PhysRevLett.19.1095 . ISSN 0031-9007 . Grunert, Katrin; Teschl, Gerald (2009), "Long-Time Asymptotics for the Korteweg–De Vries Equation via Nonlinear Steepest Descent", Math. Phys. Anal. Geom. , vol. 12, no. 3, pp. 287– 324, arXiv : 0807.5041 , Bibcode : 2009MPAG...12..287G , doi : 10.1007/s11040-009-9062-2 , S2CID 8740754 Korteweg, DJ; de Vries, G. (1895). 「XLI. 長方形の運河を進む長波の形態変化と新しいタイプの長定常波について」.ロンドン、エディンバラ、ダブリン哲学雑誌・科学ジャーナル . 39 (240): 422– 443. doi : 10.1080/14786449508620739 . ISSN 1941-5982 . Lax, Peter D. (1968). 「発展の非線形方程式と孤立波の積分」. Communications on Pure and Applied Mathematics . 21 (5): 467– 490. Bibcode : 1968CPAM...21..467L . doi : 10.1002/cpa.3160210503 . ISSN 0010-3640 . OSTI 4522657 . 三浦, ロバート・M.; ガードナー, クリフォード・S.;クラスカル, マーティン・D. (1968), 「コルテウェグ–ド・フリース方程式とその一般化。II. 保存則と運動定数の存在」, J. Math. Phys. , 9 (8): 1204– 1209, Bibcode : 1968JMP.....9.1204M , doi : 10.1063/1.1664701 , MR 0252826 ポリアニン, アンドレイ・D.; ザイツェフ, ヴァレンティン・F. (2003). 『非線形偏微分方程式ハンドブック 』 フロリダ州ボカラトン: チャップマン・アンド・ホール/CRC. ISBN 978-1-58488-355-5 。 Vakakis, Alexander F. (2002). 『非線形システムにおける正規モードと局在』 ドルドレヒト; ボストン: Springer Science & Business Media. ISBN 978-0-7923-7010-9 。 Zabusky, NJ; Kruskal, MD (1965). 「無衝突プラズマにおける『ソリトン』の相互作用と初期状態の再帰」. Physical Review Letters . 15 (6): 240– 243. Bibcode : 1965PhRvL..15..240Z . doi : 10.1103/PhysRevLett.15.240 . ISSN 0031-9007 .
外部リンク