数値解析 第7回 (2) ルジャンドル多項式 $P_n(x)$
ルジャンドル多項式の定義
ガウスの積分公式はルジャンドル多項式の性質に基づいて設計されています。
Def.1 ルジャンドル多項式 $P_n(x)$ とは次式で決まる $n$ 次式である：
$\dps{P_n(x)=\frac{1}{2^n \times n!}\left(\frac{d}{dx}\right)^n(x^2-1)^n}$
$(\,n=0,1,\cdots\,)$
Ex.2 $n=0,1,2$ のときは
\begin{align}
P_0(x)&=\frac{1}{2^0 \times 0!}(x^2-1)^0=\frac{1}{1 \times 1}1=1 \\
P_1(x)&=\frac{1}{2^1 \times 1!}(x^2-1)'=\frac{1}{2 \times 1}(2x)=x \\
P_2(x)&=\frac{1}{2^2 \times 2!}\{(x^2-1)^2\}''=\frac{1}{4 \times 2}(x^4-2x^2+1)''=\frac{1}{2}(3x^2-1) \\
\end{align}
ルジャンドル多項式の性質
証明は省略しますが、基本的な定理を３つ。
Th.3
ルジャンドル多項式は次の漸化式から計算できる
$P_0(x)=1$, $P_1(x)=x$,
$\dps{P_n(x)=\frac{1}{n}\Big\{(2n-1)xP_{n-1}(x)-(n-1)P_{n-2}(x)\Big\}}$
$(\,n=2,3,\cdots\,)$
Th.4 $P_n(1)=1$, $P_n(-1)=(-1)^n$.
Th.5 内積
$\newcommand{\ip}[2]{\langle\,#1, #2\,\rangle}$ $\dps{\ip{f}{g}=\int_{-1}^1f\,g\,dx}$ に関して次が成り立つ:
- $m \neq n$ ならば $\ip{P_m}{P_n}=0$.
- $\dps{\ip{P_n}{P_n}=\frac{1}{2n+1}}$.
Th.5 の式は、関数をベクトルと考えるとき、
- $P_n$ たちは互いに直交していて、
- $P_n$ の長さは$\dps{\frac{1}{\sqrt{2n+1}}}$ である
ことを意味しています。(1) の意味で「 $P_n$ たちは直交関数系を成す」と言います。
言わば、こんなイメージです:
Cor.6 $k \lt m$ ならば $\ip{P_m}{x^k}=0$.
証明 $P_n$ は $n$ 次式なので、$x^k$ は
$P_0$, $P_1$, $P_2$, $\cdots$, $P_k$
の一次結合として書けることが帰納法で示せます:
$\dps{x^k=\sum_{n=0}^k c_n P_n}$,
$\exists \,c_k$
$\Big($ たとえば $\dps{x^2=\frac{1}{3}P_0+\frac{2}{3}P_2}$ $\Big)$
よって
Th.5 (1) と、内積の線形性より
$\dps{\ip{P_m}{x^k}=\ip{P_m}{\sum_{n=0}^k c_n P_n}=\sum_{n=0}^k c_n\ip{P_m}{P_n}=0}$.
(証明終)
ルジャンドル多項式の由来、他
Rem.7 電磁気学、熱力学、流体力学などで活躍する 3 次元のラプラス方程式を変数変換すると、
パラメータ $n$ の入ったルジャンドルの微分方程式というものが出てきます:
$\dps{(1-x^2)y''-2xy' + n(n+1)y=0}$
これは 2 階線形微分方程式なのでその解空間は 2 次元のベクトル空間になりますが、
その基底として、多項式の解と、無限級数の解を取ることができます。
その多項式解がルジャンドル多項式 $P_n$ です。
役に立つ関数は得てして、物理学などの、数学以外の分野に由来するという例です。
Rem.8 一次独立な $n$ 個のベクトル
$\newcommand{\vecv}{\boldsymbol{v}}$$\vecv_1$, $\vecv_2$, $\cdots$, $\vecv_n$
から、直交する $n$ 個のベクトルを作り出す「グラム・シュミットの直交化法」というアルゴリズムがあります。
関数列
$1$, $x$, $\cdots$, $x^n$
にこのアルゴリズムを適用するとルジャンドル多項式
$P_0$, $P_1$, $\cdots$, $P_n$
が得られます。