3段階アルゴリズム:初期解を初期散乱データに変換し、初期散乱データを進化させ、進化した散乱データを進化した解に変換する 数学 において、逆散乱変換 (または非線形フーリエ変換)は、 波の散乱 に関連する数学的手法を用いて非線形 偏微分方程式 の初期値問題 を解く手法である。 : 4960 直接散乱変換は、関数 が どのように波を散乱するか、または束縛状態 を生成するかを記述する。 : 39–43 逆散乱変換は 、波の散乱データを用いて波の散乱を担う関数を構築する。 : 66–67 直接散乱変換と逆散乱変換は、線形 偏微分方程式を解くために使用される直接フーリエ変換と逆フーリエ変換 に類似している。 : 66–67
一対の微分演算子 を用いた3段階アルゴリズムで非線形微分方程式 を解くことができる。すなわち、初期解を散乱データに変換し(直接散乱変換)、散乱データを時間的に前進させ(時間発展)、散乱データを時間的に前進させて解を再構築する(逆散乱変換)。:66–67
このアルゴリズムは、非線形偏微分方程式を解くことを2つの線形常微分方程式 と1つの常積分方程式を解くことに簡略化し、最終的には、解くのが困難な多くの非線形偏微分方程式の 解析解 を導く方法である。:72
逆散乱問題は、少なくとも1次元空間の方程式の場合、リーマン・ヒルベルト分解問題と等価である。 この定式化は、2以上の次数の微分演算子や周期問題にも一般化できる。 より高い空間次元では、代わりに「非局所的」リーマン・ヒルベルト分解問題(乗算の代わりに畳み込みを使用)またはdバー問題が存在する。
歴史 逆散乱変換は孤立波の研究から生まれました。J.S .ラッセルは 浅瀬で発生する「並進波」または「孤立波」について説明しました。最初にJ.V.ブシネスクが 、後にD.コルテウェグ とG.デフリースが 、これらの波を記述する非線形偏微分方程式であるコルテウェグ・デフリース(KdV)方程式 を発見しました。 その後、N.ザブスキーとM.クラスカルは、フェルミ・パスタ・ウラム・ツィンゴウ問題を 数値解析的に調べ、孤立波は衝突粒子の弾性特性を持つことを発見しました。つまり、波の初期および最終的な振幅と速度は、波の衝突後も変化しません。これらの粒子のような波はソリトン と呼ばれ、分散効果と非線形効果の間の弱いバランスのために非線形方程式で発生します。
ガードナー、グリーン、クラスカル、三浦は、コルテヴェク・ド・フリース方程式 を解くための逆散乱変換を導入した。 ラックス、アブロウィッツ、カウプ、ニューウェル、セギュールはこのアプローチを一般化し、非線形シュレーディンガー方程式 、サインゴードン方程式 、修正コルテヴェク・ド・フリース方程式 、カドムツェフ・ペトビアシビリ方程式 、石森方程式 、戸田格子方程式、 ダイム方程式 などの他の非線形方程式を解くことに成功した。 このアプローチは、微分差分方程式、偏差分方程式、多次元方程式、分数積分非線形システムなど、さまざまな種類の非線形方程式にも適用されている。
説明
非線形偏微分方程式 独立変数は空間変数と時間変数である。添え字または微分演算子()は微分を表す。関数 は、初期条件(値) を持つ非線形偏微分方程式 の解である。 : 72 × {\displaystyle x} t {\displaystyle t} ∂ × 、 ∂ t {\textstyle \partial _{x},\partial _{t}} あなた ( × 、 t ) {\displaystyle u(x,t)} あなた t + 北 ( あなた ) = 0 {\textstyle u_{t}+N(u)=0} あなた ( × 、 0 ) {\textstyle u(x,0)}
要件 微分方程式の解は積分可能性とファデーエフ条件を満たす: : 40
積分可能性条件:∫ − ∞ ∞ | あなた ( × ) | d × < ∞ {\displaystyle \int _{-\infty }^{\infty }\ |u(x)|\ dx\ <\infty } ファデエフの状態:∫ − ∞ ∞ ( 1 + | × | ) | あなた ( × ) | d × < ∞ {\displaystyle \int _{-\infty }^{\infty }\ (1+|x|)|u(x)|\ dx\ <\infty }
微分演算子ペア Lax微分作用素 、 、は、係数に関数またはその導関数を含む可能性のある線型常微分作用素である。自己随伴作用素は 時間微分を持ち、固有関数 と時定数固有値 (スペクトルパラメータ )を持つ固有値(スペクトル)方程式 を生成する。 : 4963 : 98 L {\textstyle L} M {\textstyle M} あなた ( × 、 t ) {\textstyle u(x,t)} L {\textstyle L} L t {\textstyle L_{t}} ψ {\textstyle \psi } λ {\textstyle \lambda }
L ( ψ ) = λ ψ 、 {\displaystyle L(\psi )=\lambda \psi ,\ } そして L t ( ψ ) = d e f ( L ( ψ ) ) t − L ( ψ t ) {\textstyle \ L_{t}(\psi ){\overset {def}{=}}(L(\psi ))_{t}-L(\psi _{t})} 演算子は、固有関数が時間の経過とともにどのように進化するかを記述し、の固有関数から の新しい固有関数演算子を生成します。 : 4963 M {\textstyle M} ψ 〜 {\textstyle {\widetilde {\psi }}} L {\textstyle L} ψ {\textstyle \psi } L {\textstyle L}
ψ ~ = ψ t − M ( ψ ) {\displaystyle {\widetilde {\psi }}=\psi _{t}-M(\psi )\ } Lax演算子を組み合わせると、固有関数の微分演算子ではなく乗法演算子が形成されます。:4963 ψ {\textstyle \psi }
( L t + L M − M L ) ψ = 0 {\displaystyle (L_{t}+LM-ML)\psi =0} Lax演算子は乗法演算子を非線形微分方程式と等しくするために選択される。:4963
L t + L M − M L = u t + N ( u ) = 0 {\displaystyle L_{t}+LM-ML=u_{t}+N(u)=0} アブロウィッツ、カウプ、ニューウェル、セガーによって開発されたAKNS微分演算子は、ラックス微分演算子の代替であり、同様の結果を達成します。[ 4964
直接散乱変換は初期散乱データを生成する。これには、この微分方程式の固有関数解の反射係数、透過係数、固有値データ、および正規化定数が含まれる場合がある。:39–48
L ( ψ ) = λ ψ {\displaystyle L(\psi )=\lambda \psi }
散乱データの時間発展 散乱データの時間的変化を記述する方程式は、時間に関する1階線形常微分方程式の解として現れる。様々なアプローチを用いることで、この1階線形微分方程式は、線形微分演算子(Laxペア、AKNSペア)、線形微分演算子と非線形微分方程式の組み合わせ、あるいは追加の代入、積分、または微分演算によって得られる。空間漸近方程式()は、これらの微分方程式の解を簡素化する。:4967–4968 :68–72 x → ± ∞ {\textstyle x\to \pm \infty }
マルチェンコ方程式は 散乱データを線形フレドホルム積分方程式 に統合する。この積分方程式の解は非線形微分方程式の解u(x,t)となる。 :48–57
例: Korteweg-De Vries 方程式非線形微分 Korteweg-De Vries 方程式は : 4 です。
u t − 6 u u x + u x x x = 0 {\displaystyle u_{t}-6uu_{x}+u_{xxx}=0}
緩い演算子 Lax演算子は以下のとおりである: : 97–102
L = − ∂ x 2 + u ( x , t ) {\displaystyle L=-\partial _{x}^{2}+u(x,t)\ } そして M = − 4 ∂ x 3 + 6 u ∂ x + 3 u x {\textstyle \ M=-4\partial _{x}^{3}+6u\partial _{x}+3u_{x}} 乗算演算子は次のとおりです。
L t + L M − M L = u t − 6 u u x + u x x x = 0 {\displaystyle L_{t}+LM-ML=u_{t}-6uu_{x}+u_{xxx}=0}
この微分方程式の解は
L ( ψ ) = − ψ x x + u ( x , 0 ) ψ = λ ψ {\textstyle L(\psi )=-\psi _{xx}+u(x,0)\psi =\lambda \psi } 連続的な固有値範囲を持つ散乱解( 連続スペクトル )と離散的な固有値を持つ境界状態 解(離散スペクトル )が含まれる場合がある。散乱データには、透過係数、左反射係数、右反射係数、離散固有値、および左および右境界状態正規化定数(ノルム定数)が 含まれる。:4960 T ( k , 0 ) {\textstyle T(k,0)} R L ( k , 0 ) {\textstyle R_{L}(k,0)} R R ( k , 0 ) {\textstyle R_{R}(k,0)} − κ 1 2 , … , − κ N 2 {\textstyle -\kappa _{1}^{2},\ldots ,-\kappa _{N}^{2}}
c ( 0 ) L j = ( ∫ − ∞ ∞ ψ L 2 ( i k j , x , 0 ) d x ) − 1 / 2 j = 1 , … , N {\displaystyle c(0)_{Lj}=\left(\int _{-\infty }^{\infty }\ \psi _{L}^{2}(ik_{j},x,0)\ dx\right)^{-1/2}\ j=1,\dots ,N} c ( 0 ) R j = ( ∫ − ∞ ∞ ψ R 2 ( i k j , x , 0 ) d x ) − 1 / 2 j = 1 , … , N {\displaystyle c(0)_{Rj}=\left(\int _{-\infty }^{\infty }\ \psi _{R}^{2}(ik_{j},x,0)\ dx\right)^{-1/2}\ j=1,\dots ,N}
散乱データの時間発展 空間的に漸近的な左および右のジョスト関数は このステップを簡素化する。:4965–4966 ψ L ( k , x , t ) {\textstyle \psi _{L}(k,x,t)} ψ R ( k , x , t ) {\textstyle \psi _{R}(k,x,t)}
ψ L ( x , k , t ) = e i k x + o ( 1 ) , x → + ∞ ψ L ( x , k , t ) = e i k x T ( k , t ) + R L ( k , t ) e − i k x T ( k , t ) + o ( 1 ) , x → − ∞ ψ R ( x , k , t ) = e − i k x T ( k , t ) + R R ( k , t ) e i k x T ( k , t ) + o ( 1 ) , x → + ∞ ψ R ( x , k , t ) = e − i k x + o ( 1 ) , x → − ∞ {\displaystyle {\begin{aligned}\psi _{L}(x,k,t)&=e^{ikx}+o(1),\ x\to +\infty \\\psi _{L}(x,k,t)&={\frac {e^{ikx}}{T(k,t)}}+{\frac {R_{L}(k,t)e^{-ikx}}{T(k,t)}}+o(1),\ x\to -\infty \\\psi _{R}(x,k,t)&={\frac {e^{-ikx}}{T(k,t)}}+{\frac {R_{R}(k,t)e^{ikx}}{T(k,t)}}+o(1),\ x\to +\infty \\\psi _{R}(x,k,t)&=e^{-ikx}+o(1),\ x\to -\infty \\\end{aligned}}} 従属定数は 右と左のジョスト関数と右と左の正規化定数を関連付けます。:4965–4966 γ j ( t ) {\textstyle \gamma _{j}(t)}
γ j ( t ) = ψ L ( x , i κ j , t ) ψ R ( x , i κ j , t ) = ( − 1 ) N − j c R j ( t ) c L j ( t ) {\displaystyle \gamma _{j}(t)={\frac {\psi _{L}(x,i\kappa _{j},t)}{\psi _{R}(x,i\kappa _{j},t)}}=(-1)^{N-j}{\frac {c_{Rj}(t)}{c_{Lj}(t)}}} Lax微分演算子は、他の固有関数の時間依存線形結合として表現できる固有関数を生成する。:4967 M {\textstyle M}
∂ t ψ L ( k , x , t ) − M ψ L ( x , k , t ) = a L ( k , t ) ψ L ( x , k , t ) + b L ( k , t ) ψ R ( x , k , t ) {\displaystyle \partial _{t}\psi _{L}(k,x,t)-M\psi _{L}(x,k,t)=a_{L}(k,t)\psi _{L}(x,k,t)+b_{L}(k,t)\psi _{R}(x,k,t)} ∂ t ψ R ( k , x , t ) − M ψ R ( x , k , t ) = a R ( k , t ) ψ L ( x , k , t ) + b R ( k , t ) ψ R ( x , k , t ) {\displaystyle \partial _{t}\psi _{R}(k,x,t)-M\psi _{R}(x,k,t)=a_{R}(k,t)\psi _{L}(x,k,t)+b_{R}(k,t)\psi _{R}(x,k,t)} これらの微分方程式の解は散乱と境界状態の空間漸近的ジョスト関数を用いて決定され、透過係数は時間一定であるが、反射係数と正規化係数は時間依存であることを示す。:4967–4968 T ( k , t ) {\textstyle T(k,t)}
R L ( k , t ) = R L ( k , 0 ) e − i 8 k 3 t R R ( k , t ) = R R ( k , 0 ) e + i 8 k 3 t c L j ( t ) = c L j ( 0 ) e + 4 κ j 3 t , j = 1 , … , N c R j ( t ) = c R j ( 0 ) e − 4 κ j 3 t , j = 1 , … , N {\displaystyle {\begin{aligned}R_{L}(k,t)&=R_{L}(k,0)e^{-i8k^{3}t}\\R_{R}(k,t)&=R_{R}(k,0)e^{+i8k^{3}t}\\c_{Lj}(t)&=c_{Lj}(0)e^{+4\kappa _{j}^{3}t},\ j=1,\ldots ,N\\c_{Rj}(t)&=c_{Rj}(0)e^{-4\kappa _{j}^{3}t},\ j=1,\ldots ,N\end{aligned}}}
マルチェンコ核は :4968–4969で ある。F ( x , t ) {\textstyle F(x,t)}
F ( x , t ) = d e f 1 2 π ∫ − ∞ ∞ R R ( k , t ) e i k x d k + ∑ j = 1 N c ( t ) L j 2 e − κ j x {\displaystyle F(x,t){\overset {def}{=}}{\frac {1}{2\pi }}\int _{-\infty }^{\infty }R_{R}(k,t)e^{ikx}\ dk+\sum _{j=1}^{N}c(t)_{Lj}^{2}e^{-\kappa _{j}x}} マルチェンコ積分方程式 は、について解かれる線形積分方程式である。:4968–4969 K ( x , y , t ) {\textstyle K(x,y,t)}
K ( x , z , t ) + F ( x + z , t ) + ∫ x ∞ K ( x , y , t ) F ( y + z , t ) d y = 0 {\displaystyle K(x,z,t)+F(x+z,t)+\int _{x}^{\infty }K(x,y,t)F(y+z,t)\ dy=0} マルチェンコ方程式の解は、非線形偏微分方程式の解を生成する。 : 4969 K ( x , y , t ) {\textstyle K(x,y,t)} u ( x , t ) {\textstyle u(x,t)}
u ( x , t ) = − 2 ∂ K ( x , x , t ) ∂ x {\displaystyle u(x,t)=-2{\frac {\partial K(x,x,t)}{\partial x}}}
積分可能な方程式の例
参照
引用
参考文献 Konopelchenko, BG; Dubrowsky, VG (1991). 「石森方程式の局所ソリトン」. Sattinger, David H.; Tracy, CA; Venakides, Stephanos (編). 『逆散乱とその応用 』アメリカ数学会 pp. 77– 90. ISBN 978-0-8218-5129-6 。 Oono, H. (1996). 「逆散乱法によるハリー・ダイム方程式のNソリトン解」. Alfinito, E.; Boiti, M.; Martina, L. (編).非線形物理学:理論と実験 . World Scientific Publishing Company Pte Limited. pp. 241– 248. ISBN 978-981-02-2559-9 。
さらに読む
外部リンク