広田の方法

提供: testwiki
ナビゲーションに移動 検索に移動

テンプレート:Differential equations 広田の方法(ひろたのほうほう、テンプレート:Lang-en-short)は、ソリトン方程式のソリトン解を求めるための方法の一つで、簡便にして強力なことで知られる。広田良吾が考案した。双線形化法 (bilinearization method)、直接法 (direct method) とも呼ばれる。

テンプレート:Math微分などによる従属変数の変数変換により、非線形偏微差分方程式を双線形方程式に変換する。変換後の従属変数はしばしば テンプレート:Mvar 関数と呼ばれる。テンプレート:Mvar 関数は行列式またはパフィアン (Pfaffian) で、双線形方程式はPlucker関係式である。

ソリトン方程式の可積分性を保ったまま方程式の独立変数を離散化する際にも重要な役割を果たしている。

広田微分

定義

二つの関数の組 テンプレート:Math, テンプレート:Math に対して、

DxmDtnfg=(xx)m(tt)nf(x,t)g(x,t)|x=x,t=t

で定義される二項演算広田微分と呼ぶ。演算子 テンプレート:Mvar, テンプレート:Mvar広田のD演算子と呼ぶ。

実際の広田微分の計算例は次のようになる。

  • Dxfg=fxgfgx
  • Dx2fg=fxxg2fxgx+fgxx
  • Dx3fg=fxxxg3fxxgx+3fxgxxfgxxx
  • Dx4fg=fxxxxg4fxxxgx+6fxxgxx4fxgxxx+fgxxxx
  • DxDtfg=ftxgftgxfxgt+fgtx

双線形形式

二つの関数の組に、広田微分を作用させた場合、各項は二つの関数の導関数について、どちらも一次式の形になっており、これを双線形形式 (bilinear form) と呼ぶ。可積分系の非線形偏微分方程式は、適当な従属変数の変換の下、双線形形式の広田微分の方程式に変形できる。シンプルな形に表現された双線形形式の方程式に、広田微分の性質を組み合わせることで、見通しのよい計算で解を構成することが可能となる。

方程式 変数変換 双線形形式
KdV方程式: ut+6uux+uxxx=0  u=22x2logf  Dx(Dt+Dx3)ff=0
mKdV方程式: ut+6u2ux+uxxx=0  u=gf (Dt+Dx3)gf=0,Dx2ff=2g2
非線形Schrödinger方程式: iut+uxx+2|u|2u=0  u=gf (fは実数値関数、gは複素数値関数) (iDt+Dx2)gf=0,Dx2ff=2gg
サイン・ゴルドン方程式: utx=sinu  u=2ilogff (fは複素数値関数)  DxDtff=12(f2f2)
戸田格子: d2dt2rn=2ernern1ern+1  Vn=ern1,Vn=d2dt2logτn  12Dt2τnτn=τn+1τn1τn2
KP方程式: x(ut+6uux+uxxx)+uyy=0  u=22x2logf  (DxDt+Dy2+Dx4)ff=0

広田の方法

広田の方法では、可積分系の非線形偏微分方程式に対し、対数微分などの従属変数の変換を行った後、広田微分を用いて、双線形形式の微分方程式に帰着させる。さらに双線形形式の微分方程式を、べき級数の形式で展開し、各べき乗のオーダーを満たす関数形を定めていくことで解を構成する。逆散乱法では、非線形偏微分方程式をシュレディンガー方程式散乱問題に帰着させ、散乱データから元の非線形偏微分方程式の解に対応するポテンシャル関数を構成するという数学的技巧を要するが、広田の方法では直接的なアプローチで元の方程式を解くことができ、簡便性が高い。

KdV方程式の例

可積分系の代表的な例であるKdV方程式で、広田の方法を説明する。KdV方程式

ut+6uux+uxxx=0

において、

u=22x2logf

なる変数変換をすると、

Dx(Dt+Dx3)ff=0

なる双線形形式の方程式に帰着される。ここで テンプレート:Mvar

f=1+ϵf1+ϵ2f2+

テンプレート:Mvar によるべき級数で展開する。これを双線形形式の方程式に代入し、各べき テンプレート:Mvar のオーダー毎にまとめると、

ϵ:Dx(Dt+Dx3)(f11+1f1)=0
ϵ2:Dx(Dt+Dx3)(f21+f1f1+1f2)=0
ϵ3:Dx(Dt+Dx3)(f31+f2f1+f1f2+1f3)=0

となる。

1ソリトン解

1ソリトン解を構成するには次のような解の構成を行う。まず、

f1=e2(κxωt)

として、テンプレート:Math の項を考えると

ω=4κ3

の関係が満される必要があることがわかる。また、高次の テンプレート:Mvar の項については、特解として、

fn=0n2

をとることができる。よって、解 テンプレート:Mvar としては

u=22x2log(1+e2(κx4κ3t))

となる。

参考文献

  • R. Hirota, Phy. Rev. Lett., 27, p. 1192, 1971. テンプレート:Doi
  • 広田良吾, "直接法によるソリトンの数理", 岩波書店, 1992年, ISBN 978-4000056762

関連事項

テンプレート:Mathanalysis-stub