2次方程式で摂動論とは何かを理解する

摂動論とは、対象としている問題Aを、厳密に解析解が求められる問題Bに小さな「ずれ」(摂動項) が加えられた問題と見なすことで、Aの近似解を求める方法です。本稿ではその分かりやすい導入として、2次方程式を例に摂動論を理解することを目指します。

 摂動論の考え方

多重振り子や天体の運動などの力学や、シュレディンガー方程式などの量子力学に関する問題を扱う際に、厳密解を求められない問題に遭遇することが多々あります。そのようなときに活躍するのが「摂動論」という考え方です（この方法を指す場合は特に「摂動法」と呼ばれます）。普通は次のように説明されることが多いです。

摂動法では基本的にこれだけの操作しか要求されないのですが、この説明だけでは摂動論が具体的に何をする方法なのか分かりませんし、全くありがたみが伝わってきませんよね。

そこで2次方程式を例に、摂動論が具体的にどんな方法で、どんなご利益があるものなのかを直感的に理解していきましょう！

 2次方程式を摂動法で解く

例えば、次の2次方程式を解くのは容易です。$$(★)x^2+x-6=0$$これは$$(x-2)(x+3)=0$$ $$\therefore x=2,-3$$と厳密に解くことができます。

では、次の2次方程式はどうでしょうか？$$(\ast )x^2+x-6=0.01$$皆さんは「2次方程式の解の公式」という便利な道具を知っていますので、方程式$(\ast )$の厳密解を求めるのは然程難しくはないでしょう。しかしこれが3次方程式や4次方程式になってくると厳密解を求めるのは困難です。そういう場合に有効なのが摂動法なのです。

では早速、2次方程式$(\ast )$を摂動法で解いてみます。

ⅰ. 摂動項$\epsilon$の導入

より一般的に考えるために、$0.01$の部分を摂動と見なして$\epsilon$と置きます。因みに摂動項としてイプシロンの文字が用いられるのは、英語の “error”（誤差）に由来しています。

さて、$$(\ast )x^2+x-6=\epsilon$$の解$x$をイプシロン$\epsilon$の冪級数として$$x(\epsilon )=x_0+\epsilon x_1+{\epsilon }^2{x}_2+\cdots$$のように「$\epsilon$の関数」として展開できると考えます（これについては、そういうものなのね、と思って下さい）。ここで$x_0$は2次方程式$(★)$の厳密解を表します。$x_1$と$x_2$は摂動にくっついているオマケの定数です。

ここでは「2次の摂動法」を用いることにします。いま、$\epsilon$は微小なので、3次以降の項はもっと微小な数値になっています。つまり、無視してもそれほど問題が無い大きさだということです。

${\epsilon }^2$という2次の項まで考えて、近似解$x$を$$x(\epsilon )=x_0+\epsilon x_1+{\epsilon }^2{x}_2$$という式で$\epsilon$の関数として表すことにします。これを2次方程式$(\ast )$に代入します。

ⅱ. $\epsilon$の式を代入して整理

$(\ast )$の式に $x(\epsilon )=x_0+\epsilon x_1+{\epsilon }^2{x}_2$ を代入して$\epsilon$について整理すると、$$\begin{array}{rl}& \{x_0^2+2\epsilon x_0{x}_1+{\epsilon }^2(x_1^2+2x_0{x}_2)\}\\ & +(x_0+\epsilon x_1+{\epsilon }^2{x}_2)-6=\epsilon \end{array}$$ $$\therefore (x_1^2+2x_0{x}_2+x_2){\epsilon }^2+(2x_0{x}_1+x_1)\epsilon +(x_0^2+x_0-6)=\epsilon$$となります。

得られた方程式の両辺の係数を比較することで、以下の連立方程式が得られます。$$\{\begin{array}{lr}{x}_0^2+x_0-6=0& \cdots \mathrm{①}\\ 2x_0{x}_1+x_1=1& \cdots \mathrm{②}\\ x_1^2+2x_0{x}_2+x_2=0& \cdots \mathrm{③}\end{array}$$次はこの連立方程式を解いていきます。

ⅲ. 連立方程式を解く

①は厳密解を与える方程式で、ここから $x_0=2$ もしくは $x_0=-3$ となります。

$x_0=2$ のとき、②より$$5x_1=1\therefore x_1={\displaystyle \frac{1}{5}}$$となり、③より$${\displaystyle \frac{1}{25}}+4x_2+x_2=0$$ $$\therefore x_2=-{\displaystyle \frac{1}{125}}$$となります。よってこのとき近似解は$$x(\epsilon )=2+{\displaystyle \frac{1}{5}}\epsilon -{\displaystyle \frac{1}{125}}{\epsilon }^2$$と求められます。

同様に $x_0=-3$ のとき、②より$$x_1=-{\displaystyle \frac{1}{5}}$$となり、③より$${\displaystyle \frac{1}{25}}-6x_2+x_2=0$$ $$\therefore x_2=\frac{1}{125}$$となるので、近似解は$$x(\epsilon )=-3-{\displaystyle \frac{1}{5}}\epsilon +{\displaystyle \frac{1}{125}}{\epsilon }^2$$と求められます。

ⅳ. 近似解の$\epsilon$に値を代入する

以上より摂動を2次まで考慮した近似解$$x(\epsilon )=\{\begin{array}{l}2+{\displaystyle \frac{1}{5}}\epsilon -{\displaystyle \frac{1}{125}}{\epsilon }^2\\ -3-{\displaystyle \frac{1}{5}}\epsilon +{\displaystyle \frac{1}{125}}{\epsilon }^2\end{array}$$が得られました。

もとの2次方程式$(\ast )$は$$(\ast )x^2+x-6=0.01$$だったので、これは $\epsilon ={\displaystyle \frac{1}{100}}$ の場合に相当します。このとき解は$$\begin{array}{rl}x& =\{\begin{array}{r}2+{\displaystyle \frac{1}{500}}-{\displaystyle \frac{1}{1250000}}\\ -3–{\displaystyle \frac{1}{500}}+{\displaystyle \frac{1}{1250000}}\end{array}\\ & =\{\begin{array}{r}2.0019992\\ -3.0019992\end{array}\end{array}$$となります。

厳密解は$$x={\displaystyle \frac{\sqrt{625}-5}{100}},-{\displaystyle \frac{\sqrt{625}+5}{100}}$$であり、それぞれの数値は 2.001999200639、-3.001999200639 なので、求められた近似解は小数点以下8桁もの精度で厳密解に一致していることが確かめられます。これは結構凄い精度ですね！

摂動法の流れは、これまで見てきたように、

ⅰ 対象としている方程式に摂動項$\epsilon$を導入する
ⅱ $\epsilon$の冪級数を$x$に代入して整理する
ⅲ 係数比較により得られる連立方程式を解く
ⅳ 近似解の$\epsilon$に求めた定数を代入する

という手順になっています。求めた近似解が厳密解とかけ離れていたらどこかで計算ミスをしています。検算も忘れずに！

 誤差の範囲と摂動の大きさ

ここで気になるのが、摂動の大きさがどのくらいの範囲までなら許容できる誤差範囲内に収められるのか、という点です。何万回も繰り返し処理するような計算では小さな誤差でも、蓄積していくと膨大なズレに繋がります。誤差の程度は数値シミュレーションにおいて死活問題と言えます。

というわけで、近似解と厳密解の相対誤差を調べてみましょう。

2次方程式$$x^2+x-6=\epsilon$$の厳密解は$$x={\displaystyle \frac{\pm \sqrt{25+4\epsilon }-1}{2}}$$であり、このうち $x=2$ 近傍の解は $x={\displaystyle \frac{\sqrt{25+4\epsilon }-1}{2}}$ です。よって、$x=2$ 近傍の厳密解と近似解の相対誤差は$$E=|1-{\displaystyle \frac{2+{\displaystyle \frac{1}{5}}\epsilon -{\displaystyle \frac{1}{125}}{\epsilon }^2}{{\displaystyle \frac{\sqrt{25+4\epsilon }-1}{2}}}}|$$により評価できます。

厳密解と近似解の相対誤差が1%以内に収まるのは $E<{\displaystyle \frac{1}{100}}$ のときですから、これを解くと$$-2.54047<\epsilon <3.85867$$という範囲が得られます。同様にして $x=-3$ 近傍の場合は$$E^{\prime }=|1-{\displaystyle \frac{-3-{\displaystyle \frac{1}{5}}\epsilon +{\displaystyle \frac{1}{125}}{\epsilon }^2}{{\displaystyle \frac{-\sqrt{25+4\epsilon }-1}{2}}}}|$$で相対誤差が評価できるので $E^{\prime }<{\displaystyle \frac{1}{100}}$ を解いて$$-2.92956<\epsilon <4.36904$$という範囲を得ます。誤差1%以内なら、摂動$\epsilon$を割と大きめにしても良いようです。

誤差の範囲を0.001%以内と厳しくすると、$x=2$ 近傍の解については $$-0.308416<\epsilon <0.321644$$ $x=-3$ 近傍の解については $$-0.353384<\epsilon <0.367828$$という範囲が得られます。誤差1%のときよりも$\epsilon$の範囲が結構狭くなりましたね。

実際、$\epsilon ={\displaystyle \frac{1}{10}}$ のときの近似解は小数点以下5桁の精度で厳密解に一致し、$\epsilon =1$ の場合でも小数点以下2桁の精度で一致します。

一方で $\epsilon =10$ と大きめの値にすると、近似解は $x={\displaystyle \frac{16}{5}}$、$-{\displaystyle \frac{21}{5}}$（3.2…, -4.2…）となり、厳密解は $x={\displaystyle \frac{\sqrt{65}-1}{2}}$、$-{\displaystyle \frac{\sqrt{65}+1}{2}}$（3.5311…, -4.5311…）となるので、さすがに無視できない大きさのズレになります。