工業数学シリーズ第8回：級数解とordinary point

これまでは特性方程式や標準公式を利用して比較的きれいに解ける問題を見てきました。しかし実際にはそのような問題ばかりではありません。ここで登場する重要なアイデアが解をべき級数で仮定する方法です。

なぜ級数解が必要なのか

ある微分方程式には初等関数できれいに表現できる解が存在しません。だからといって解を諦めるわけではありません。代わりに解を

$y(x) = \sum_{n=0}^{\infty} a_n x^n$

のように無限級数で仮定し、係数を見つけていくことができます。

この方式は後にベッセル関数、ルジャンドル関数などの特殊関数につながる出発点です。

ordinary pointとは何か

2次線形方程式を

$y'' + P(x)y' + Q(x)y = 0$

の形で書いたとき、展開したい点 $x_0$ の近くで $P(x)$ と $Q(x)$ が解析的であれば、その点をordinary pointと考えることができます。入門段階では「その点の近くで係数関数が十分に滑らかで、最高次微分項の係数が0でないため安全に展開できる点」くらいに理解すれば良いでしょう。

基本手順

べき級数解法の流れは通常以下の通りです。

解を $\sum a_n x^n$ で仮定する
微分して $y'$ 、 $y''$ も級数で書く
元の式に代入する
同じ次数の項をまとめる
係数比較で漸化式を得る

つまり計算の核心は「無限級数を微分方程式の中に押し込んで係数を比較すること」です。

手で見る例題

次の式を見てみましょう。

$y'' + y = 0$

解を

$y = \sum_{n=0}^{\infty} a_n x^n$

とすると

$y'' = \sum_{n=0}^{\infty} (n+2)(n+1)a_{n+2}x^n$

となります。したがって

$\sum_{n=0}^{\infty} \left[(n+2)(n+1)a_{n+2} + a_n\right]x^n = 0$

で、各次数の係数はすべて0でなければならないので

$a_{n+2} = -\frac{a_n}{(n+2)(n+1)}$

を得ます。

この漸化式は偶数係数と奇数係数がそれぞれ別々に決まるという意味で、実際に展開するとコサインとサインの級数が現れます。

この方法が与える意味

この例題はすでに解がわかっている方程式ですが、級数解がどのように機能するかを見るのに良い例です。重要なのは「解を一つの関数で一発で当てられなくても、係数を順番に作っていける」という点です。

工学応用

特殊関数の出発点

円筒座標や球面座標の問題を解いていくとベッセル関数、ルジャンドル関数が登場します。これらの関数の多くは級数解から定義されたり理解されたりします。

近似計算

ある点の近くで解を素早く近似する必要があるとき、べき級数は非常に有用です。

アルゴリズム実装

開発者の観点では、漸化式はそのままアルゴリズムです。係数生成規則をコードに変換すれば数値計算ライブラリの一部になります。

よくある間違い

指数のインデックスを整理しない

級数解はインデックスの移動をきれいにしないとすぐに混乱します。 $n$ 、 $n+1$ 、 $n+2$ を正確に合わせる習慣が重要です。

ordinary pointの概念を重く考えすぎる

入門段階では「展開したい点の近くで係数が問題ないか」を確認する程度で始めても十分です。

漸化式まで書いて意味を見失う

漸化式は単純な計算装置ではなく、解の係数構造と自由度の数を示す情報です。

一行まとめ

級数解は閉形式の解をすぐに見つけにくい微分方程式で、解を係数の規則として復元する強力な方法です。

次回予告

次の記事ではODEパートを一旦整理しながら、問題を積分問題に変えてくれる強力なツールであるラプラス変換に移ります。

参考資料

Erwin Kreyszig, Advanced Engineering Mathematics, 10th Edition
George B. Arfken, Mathematical Methods for Physicists
James Ward Brown, Ruel V. Churchill, Complex Variables and Applications

工業数学シリーズ 第8回：級数解とordinary point