完全直交多項式系
最初の6つのルジャンドル多項式
数学 において 、 ルジャンドル多項式は、 アドリアン・マリー・ルジャンドル (1782)にちなんで名付けられた 、完全な 直交多項式 の体系であり、さまざまな数学的特性と多数の応用を備えています。ルジャンドル多項式はさまざまな方法で定義でき、さまざまな定義によってさまざまな側面が強調されるとともに、さまざまな数学的構造や物理的および数値的応用への一般化や関連性が示唆されます。
ルジャンドル多項式に密接に関連しているものとして、 ルジャンドル陪関数 、 ルジャンドル関数 、第 2 種ルジャンドル関数、 大 q-ルジャンドル多項式 、 ルジャンドル陪関数 などがあります。
定義と表現
直交系としての構成による定義
このアプローチでは、多項式は区間 上の 重み関数に関する直交系として定義されます 。つまり、 は次数 の多項式であり 、
わ
(
x
)
=
1
{\displaystyle w(x)=1}
[
−
1
、
1
]
{\displaystyle [-1,1]}
ポ
ん
(
x
)
{\displaystyle P_{n}(x)}
ん
{\displaystyle n}
∫
−
1
1
ポ
メートル
(
x
)
ポ
ん
(
x
)
d
x
=
0
もし
ん
≠
メートル
。
{\displaystyle \int _{-1}^{1}P_{m}(x)P_{n}(x)\,dx=0\quad {\text{if }}n\neq m.}
追加の標準化条件 により 、すべての多項式を一意に決定できます。次に、構築プロセスを開始します。 は、 次数 0 の唯一の正しく標準化された多項式です。は に直交する必要があり、 につながります。 また、 は、 および への直交を要求することで決定され 、以下同様に続きます。は、 で すべて に直交することを要求することで固定されます。これにより、 標準化とともにの すべての 係数を固定する条件 が得られます 。 の作業により、すべての多項式のすべての係数を体系的に決定でき、以下に示す の累乗の明示的な表現につながります 。
ポ
ん
(
1
)
=
1
{\displaystyle P_{n}(1)=1}
ポ
0
(
x
)
=
1
{\displaystyle P_{0}(x)=1}
ポ
1
(
x
)
{\displaystyle P_{1}(x)}
ポ
0
{\displaystyle P_{0}}
ポ
1
(
x
)
=
x
{\displaystyle P_{1}(x)=x}
ポ
2
(
x
)
{\displaystyle P_{2}(x)}
ポ
0
{\displaystyle P_{0}}
ポ
1
{\displaystyle P_{1}}
ポ
ん
{\displaystyle P_{n}}
ポ
メートル
{\displaystyle P_{m}}
メートル
<
ん
{\displaystyle m<n}
ん
{\displaystyle n}
ポ
ん
(
1
)
=
1
{\displaystyle P_{n}(1)=1}
ん
+
1
{\displaystyle n+1}
ポ
ん
(
x
)
{\displaystyle P_{n}(x)}
x
{\displaystyle x}
のこの定義は 最も単純なものです。微分方程式の理論には関係ありません。第 2 に、多項式の完全性は、べき乗 1, の完全性から直接導かれます。最後に、 上の ルベーグ測度 に関する直交性を介してそれらを定義することにより 、ルジャンドル多項式を 3 つの 古典的な直交多項式システム の 1 つとして設定します。他の 2 つは、 重み を持つ 半直線上で直交する ラゲール多項式 と、 重み を持つ 全直線上で直交する エルミート多項式 です。
ポ
ん
{\displaystyle P_{n}}
x
、
x
2
、
x
3
、
…
{\displaystyle x,x^{2},x^{3},\ldots }
[
−
1
、
1
]
{\displaystyle [-1,1]}
[
0
、
∞
)
{\displaystyle [0,\infty )}
e
−
x
{\displaystyle e^{-x}}
(
−
∞
、
∞
)
{\displaystyle (-\infty ,\infty )}
e
−
x
2
{\displaystyle e^{-x^{2}}}
生成関数による定義
ルジャンドル多項式は、生成関数 の べき乗の形式的展開における係数として定義することもできる [1]
t
{\displaystyle t}
の係数は、 の 次数 の多項式です 。 まで展開すると、 が得られ
ます。
より高い次数への展開はますます面倒になりますが、体系的に行うことは可能であり、再び以下に示す明示的な形式のいずれかにつながります。
t
ん
{\displaystyle t^{n}}
x
{\displaystyle x}
ん
{\displaystyle n}
|
x
|
≤
1
{\displaystyle |x|\leq 1}
t
1
{\displaystyle t^{1}}
ポ
0
(
x
)
=
1
、
ポ
1
(
x
)
=
x
。
{\displaystyle P_{0}(x)=1\,,\quad P_{1}(x)=x.}
しかし、テイラー級数 を直接展開しなくても、 より高い を得ることは可能です 。 式 2 は両辺を
t について微分し、次のように整理されます。 平方根の商を式 2
の定義に置き換え 、 結果として得られる展開で t の累乗の 係数を等しくすると、 ボネットの再帰式が得られます。 この関係と最初の 2 つの多項式 P 0 と P 1
により 、残りはすべて再帰的に生成できます。
ポ
ん
{\displaystyle P_{n}}
x
−
t
1
−
2
x
t
+
t
2
=
(
1
−
2
x
t
+
t
2
)
∑
ん
=
1
∞
ん
ポ
ん
(
x
)
t
ん
−
1
。
{\displaystyle {\frac {xt}{\sqrt {1-2xt+t^{2}}}}=\left(1-2xt+t^{2}\right)\sum _{n=1}^{\infty }nP_{n}(x)t^{n-1}\,.}
(
ん
+
1
)
ポ
ん
+
1
(
x
)
=
(
2
ん
+
1
)
x
ポ
ん
(
x
)
−
ん
ポ
ん
−
1
(
x
)
。
{\displaystyle (n+1)P_{n+1}(x)=(2n+1)xP_{n}(x)-nP_{n-1}(x)\,.}
生成関数アプローチは、以下で説明するように、静電気学における 多極展開 に直接関連しており、1782 年にルジャンドルによって多項式が初めて定義された方法です。
微分方程式による定義
3番目の定義はルジャンドルの微分方程式 の解に基づいている 。
この 微分方程式は x = ±1 で 正則特異点 を持つ ため、標準的な フロベニウス 法または べき級数 法を使用して解を求める場合、一般に 原点の周りの級数は | x | < 1でのみ収束します。 nが整数のとき、 x = 1 で正則な 解 P n ( x )は x = −1 でも正則であり 、この解の級数は終了します(つまり、多項式です)。 これらの解の直交性と完全性は、 スチュム・リウヴィル理論 の観点から最もよくわかります。の代わりに
固有値を使用して 、微分方程式を固有値問題として書き直します
。 で解が正則であることを要求すると
、 左側の 微分演算子は エルミートになります。 固有値は n ( n + 1) の形式であることがわかり
、 であり 、固有関数は です 。 この解の集合の直交性と完全性は、スチュム・リウヴィル理論のより大きな枠組みからすぐにわかります。
d
d
x
(
(
1
−
x
2
)
d
d
x
)
ポ
(
x
)
=
−
λ
ポ
(
x
)
、
{\displaystyle {\frac {d}{dx}}\left(\left(1-x^{2}\right){\frac {d}{dx}}\right)P(x)=-\lambda P(x)\,,}
λ
{\displaystyle \lambda}
ん
(
ん
+
1
)
{\displaystyle n(n+1)}
x
=
±
1
{\displaystyle x=\pm 1}
ん
=
0
、
1
、
2
、
…
{\displaystyle n=0,1,2,\ldots }
ポ
ん
(
x
)
{\displaystyle P_{n}(x)}
微分方程式には、別の非多項式解、 第2種のルジャンドル関数があります。(式 1 ) の2パラメータの一般化は ルジャンドルの 一般 微分方程式と呼ばれ、 ルジャンドルの準多項式 で解かれます。 ルジャンドル関数は、ルジャンドルの微分方程式(一般化または非一般)の 非整数 パラメータによる解です 。
質問
ん
{\displaystyle Q_{n}}
物理的な設定では、 ラプラス方程式 (および関連する 偏微分方程式)を 球面座標 で変数分離して解くときはいつでも、ルジャンドルの微分方程式が自然に生じます 。この観点から、ラプラス演算子の角度部分の固有関数は 球面調和関数 であり、そのルジャンドル多項式は(乗法定数を除いて)極軸の周りの回転によって不変のままになる部分集合です。多項式は のように表され、 は 極角です。ルジャンドル多項式へのこのアプローチは、回転対称性との深いつながりを提供します。解析手法を通じて苦労して見つけられるルジャンドル多項式の特性の多くは、たとえば加法定理など、対称性および群論の方法を使用するとより簡単に見つけられ、深い物理的および幾何学的意味を持ちます。
ポ
ん
(
コス
θ
)
{\displaystyle P_{n}(\cos \theta )}
θ
{\displaystyle \theta}
ルジャンドル多項式の特に簡潔な表現は ロドリゲスの公式 で与えられる:
ポ
ん
(
x
)
=
1
2
ん
ん
!
d
ん
d
x
ん
(
x
2
−
1
)
ん
。
{\displaystyle P_{n}(x)={\frac {1}{2^{n}n!}}{\frac {d^{n}}{dx^{n}}}(x^{2}-1)^{n}\,.}
この公式により、の多数の特性を導出することができる 。これらの中には、次のような明示的な表現がある。
ポ
ん
{\displaystyle P_{n}}
ポ
ん
(
x
)
=
[
t
ん
]
(
(
t
+
x
)
2
−
1
)
ん
2
ん
=
[
t
ん
]
(
t
+
x
+
1
)
ん
(
t
+
x
−
1
)
ん
2
ん
、
ポ
ん
(
x
)
=
1
2
ん
∑
け
=
0
ん
(
ん
け
)
2
(
x
−
1
)
ん
−
け
(
x
+
1
)
け
、
ポ
ん
(
x
)
=
∑
け
=
0
ん
(
ん
け
)
(
ん
+
け
け
)
(
x
−
1
2
)
け
、
ポ
ん
(
x
)
=
1
2
ん
∑
け
=
0
⌊
ん
/
2
⌋
(
−
1
)
け
(
ん
け
)
(
2
ん
−
2
け
ん
)
x
ん
−
2
け
、
ポ
ん
(
x
)
=
2
ん
∑
け
=
0
ん
x
け
(
ん
け
)
(
ん
+
け
−
1
2
ん
)
、
ポ
ん
(
x
)
=
{
1
π
∫
0
π
(
x
+
x
2
−
1
⋅
コス
(
t
)
)
ん
d
t
もし
|
x
|
>
1
、
x
ん
もし
|
x
|
=
1
、
2
π
⋅
x
ん
⋅
|
x
|
⋅
∫
|
x
|
1
t
−
ん
−
1
t
2
−
x
2
⋅
コス
(
ん
⋅
アルコス
(
t
)
)
罪
(
アルコス
(
t
)
)
d
t
もし
0
<
|
x
|
<
1
、
(
−
1
)
ん
/
2
⋅
2
−
ん
⋅
(
ん
ん
/
2
)
もし
x
=
0
そして
ん
平
、
0
もし
x
=
0
そして
ん
奇数
。
{\displaystyle {\begin{aligned}P_{n}(x)&=[t^{n}]{\frac {\left((t+x)^{2}-1\right)^{n}}{2^{n}}}=[t^{n}]{\frac {\left(t+x+1\right)^{n}\left(t+x-1\right)^{n}}{2^{n}}},\\[1ex]P_{n}(x)&={\frac {1}{2^{n}}}\sum _{k=0}^{n}{\binom {n}{k}}^{\!2}(x-1)^{n-k}(x+1)^{k},\\[1ex]P_{n}(x)&=\sum _{k=0}^{n}{\binom {n}{k}}{\binom {n+k}{k}}\left({\frac {x-1}{2}}\right)^{\!k},\\[1ex]P_{n}(x)&={\frac {1}{2^{n}}}\sum _{k=0}^{\left\lfloor n/2\right\rfloor }\left(-1\right)^{k}{\binom {n}{k}}{\binom {2n-2k}{n}}x^{n-2k},\\[1ex]P_{n}(x)&=2^{n}\sum _{k=0}^{n}x^{k}{\binom {n}{k}}{\binom {\frac {n+k-1}{2}}{n}},\\[1ex]P_{n}(x)&={\begin{cases}\displaystyle {\frac {1}{\pi }}\int _{0}^{\pi }{\left(x+{\sqrt {x^{2}-1}}\cdot \cos(t)\right)}^{n}\,dt&{\text{if }}|x|>1,\\x^{n}&{\text{if }}|x|=1,\\\displaystyle {\frac {2}{\pi }}\cdot x^{n}\cdot |x|\cdot \int _{|x|}^{1}{\frac {t^{-n-1}}{\sqrt {t^{2}-x^{2}}}}\cdot {\frac {\cos \left(n\cdot \arccos(t)\right)}{\sin \left(\arccos(t)\right)}}\,dt&{\text{if }}0<|x|<1,\\\displaystyle (-1)^{n/2}\cdot 2^{-n}\cdot {\binom {n}{n/2}}&{\text{if }}x=0{\text{ and }}n{\text{ even}},\\0&{\text{if }}x=0{\text{ and }}n{\text{ odd}}.\end{cases}}\end{aligned}}}
多項式をべき級数 として表すと 、 のべき乗の係数 も一般的な公式を使用して計算できます。ルジャンドル多項式は 、 が奇数の 場合 、が偶数の 場合 の 2つの定数 と に使用される値によって決まります 。 [2]
P
n
(
x
)
=
∑
a
k
x
k
{\textstyle P_{n}(x)=\sum a_{k}x^{k}}
x
{\displaystyle x}
a
k
+
2
=
−
(
n
−
k
)
(
n
+
k
+
1
)
(
k
+
2
)
(
k
+
1
)
a
k
.
{\displaystyle a_{k+2}=-{\frac {(n-k)(n+k+1)}{(k+2)(k+1)}}a_{k}.}
a
0
{\textstyle a_{0}}
a
1
{\textstyle a_{1}}
a
0
=
0
{\textstyle a_{0}=0}
n
{\displaystyle n}
a
1
=
0
{\textstyle a_{1}=0}
n
{\displaystyle n}
4 番目の表現では、は 以下の最大の整数 を表します 。 最後の表現も再帰式から直接得られ、ルジャンドル多項式を単純な単項式で表し、 二項係数 の一般化された形式 を含みます。
⌊
n
/
2
⌋
{\displaystyle \lfloor n/2\rfloor }
n
/
2
{\displaystyle n/2}
最初のいくつかのルジャンドル多項式は次のとおりです。
これらの多項式のグラフ( n = 5 まで)を以下に示します。
最初の 6 つのルジャンドル多項式のプロット。
主な特性
直交性
標準化は、 ルジャンドル多項式の正規化を固定します( 区間 −1 ≤ x ≤ 1の L 2 ノルム に関して)。これらは 同じノルムに関して 直交しているので、2 つのステートメント [ 明確化が必要 ] を 1 つの方程式に組み合わせることができます
(ここで δ mn は クロネッカーのデルタ を表し、 m = n の場合は 1、それ以外の場合は 0 になります)。この正規化は、以下に示す ロドリゲスの公式 を使用することによって最も簡単に見つけることができます 。
P
n
(
1
)
=
1
{\displaystyle P_{n}(1)=1}
∫
−
1
1
P
m
(
x
)
P
n
(
x
)
d
x
=
2
2
n
+
1
δ
m
n
,
{\displaystyle \int _{-1}^{1}P_{m}(x)P_{n}(x)\,dx={\frac {2}{2n+1}}\delta _{mn},}
完全
多項式が完全であるということは、次のことを意味する。 区間 [−1, 1] に有限個の不連続点を持つ任意の区分連続関数が与えられたとき、和の列は
平均で に収束する 。 ただし、
f
(
x
)
{\displaystyle f(x)}
f
n
(
x
)
=
∑
ℓ
=
0
n
a
ℓ
P
ℓ
(
x
)
{\displaystyle f_{n}(x)=\sum _{\ell =0}^{n}a_{\ell }P_{\ell }(x)}
f
(
x
)
{\displaystyle f(x)}
n
→
∞
{\displaystyle n\to \infty }
a
ℓ
=
2
ℓ
+
1
2
∫
−
1
1
f
(
x
)
P
ℓ
(
x
)
d
x
.
{\displaystyle a_{\ell }={\frac {2\ell +1}{2}}\int _{-1}^{1}f(x)P_{\ell }(x)\,dx.}
この完全性特性は、この記事で説明するすべての展開の基礎となっており、多くの場合、
−1 ≤ x ≤ 1 および −1 ≤ y ≤ 1
の形式で表現されます 。
∑
ℓ
=
0
∞
2
ℓ
+
1
2
P
ℓ
(
x
)
P
ℓ
(
y
)
=
δ
(
x
−
y
)
,
{\displaystyle \sum _{\ell =0}^{\infty }{\frac {2\ell +1}{2}}P_{\ell }(x)P_{\ell }(y)=\delta (x-y),}
アプリケーション
逆距離ポテンシャルの拡大
ルジャンドル多項式は、1782年に アドリアン・マリー・ルジャンドル [3]によって、 ニュートンポテンシャル
の展開における係数として初めて導入されました
。ここで、 r と r ′は それぞれベクトル x と x ′ の長さ、 γは これら2つのベクトル間の角度です。この級数は、 r > r ′ のときに収束します。この式は、 質点 に関連付けられた 重力ポテンシャル 、または 点電荷 に関連付けられた クーロンポテンシャル を与えます。ルジャンドル多項式を使用した展開は、たとえば、連続的な質量または電荷分布にわたってこの式を積分する場合に役立ちます。
1
|
x
−
x
′
|
=
1
r
2
+
r
′
2
−
2
r
r
′
cos
γ
=
∑
ℓ
=
0
∞
r
′
ℓ
r
ℓ
+
1
P
ℓ
(
cos
γ
)
,
{\displaystyle {\frac {1}{\left|\mathbf {x} -\mathbf {x} '\right|}}={\frac {1}{\sqrt {r^{2}+{r'}^{2}-2r{r'}\cos \gamma }}}=\sum _{\ell =0}^{\infty }{\frac {{r'}^{\ell }}{r^{\ell +1}}}P_{\ell }(\cos \gamma ),}
ルジャンドル多項式は、 電荷のない空間領域における 静的 ポテンシャル ∇ 2 Φ( x ) = 0の ラプラス方程式 の解に、 変数分離 法を用いて現れる。この場合、 境界条件は軸対称性を持つ( 方位角 に依存しない )。ここで、 ẑ は対称軸、 θ は観測者の位置と ẑ 軸(天頂角)の間の角度であり 、ポテンシャルの解は次のようになる。
Φ
(
r
,
θ
)
=
∑
ℓ
=
0
∞
(
A
ℓ
r
ℓ
+
B
ℓ
r
−
(
ℓ
+
1
)
)
P
ℓ
(
cos
θ
)
.
{\displaystyle \Phi (r,\theta )=\sum _{\ell =0}^{\infty }\left(A_{\ell }r^{\ell }+B_{\ell }r^{-(\ell +1)}\right)P_{\ell }(\cos \theta )\,.}
A l と B l は 各問題の境界条件に応じて決定される。 [4]
これらは、中心力について 3 次元の
シュレーディンガー方程式 を解くときにも現れます。
多極展開では
電位の多極展開の図。
ルジャンドル多項式は、次の形式の関数を展開する場合にも役立ちます (これは前と同じですが、少し書き方が異なります)。
これは、 多重極展開
で自然に発生します 。方程式の左側は、 ルジャンドル多項式の
生成関数です。
1
1
+
η
2
−
2
η
x
=
∑
k
=
0
∞
η
k
P
k
(
x
)
,
{\displaystyle {\frac {1}{\sqrt {1+\eta ^{2}-2\eta x}}}=\sum _{k=0}^{\infty }\eta ^{k}P_{k}(x),}
例えば、 z 軸上の z = a に位置する 点電荷 による 電位 Φ( r , θ ) ( 球座標 )は、次のように変化する(右図参照)。
Φ
(
r
,
θ
)
∝
1
R
=
1
r
2
+
a
2
−
2
a
r
cos
θ
.
{\displaystyle \Phi (r,\theta )\propto {\frac {1}{R}}={\frac {1}{\sqrt {r^{2}+a^{2}-2ar\cos \theta }}}.}
観測点 P の半径 rが a より大きい場合 、ポテンシャルはルジャンドル多項式で展開される可能性があり
、 η =
Φ
(
r
,
θ
)
∝
1
r
∑
k
=
0
∞
(
a
r
)
k
P
k
(
cos
θ
)
,
{\displaystyle \Phi (r,\theta )\propto {\frac {1}{r}}\sum _{k=0}^{\infty }\left({\frac {a}{r}}\right)^{k}P_{k}(\cos \theta ),}
1つの / r < 1 かつ x = cos θ 。この展開は、通常の 多重極展開を 展開するために使用されます。
逆に、観測点 P の半径 rが a より小さい場合 、ポテンシャルは上記と同様にルジャンドル多項式で展開されますが、 a と r は交換されます。この展開は内部多重極展開の基礎となります。
三角法では
三角関数 cos nθ は 、チェビシェフ多項式 T n (cos θ ) ≡ cos nθ とも表記され 、ルジャンドル多項式 P n (cos θ ) によって多重極展開することもできます。最初のいくつかの次数は次のとおりです。
T
0
(
cos
θ
)
=
1
=
P
0
(
cos
θ
)
,
T
1
(
cos
θ
)
=
cos
θ
=
P
1
(
cos
θ
)
,
T
2
(
cos
θ
)
=
cos
2
θ
=
1
3
(
4
P
2
(
cos
θ
)
−
P
0
(
cos
θ
)
)
,
T
3
(
cos
θ
)
=
cos
3
θ
=
1
5
(
8
P
3
(
cos
θ
)
−
3
P
1
(
cos
θ
)
)
,
T
4
(
cos
θ
)
=
cos
4
θ
=
1
105
(
192
P
4
(
cos
θ
)
−
80
P
2
(
cos
θ
)
−
7
P
0
(
cos
θ
)
)
,
T
5
(
cos
θ
)
=
cos
5
θ
=
1
63
(
128
P
5
(
cos
θ
)
−
56
P
3
(
cos
θ
)
−
9
P
1
(
cos
θ
)
)
,
T
6
(
cos
θ
)
=
cos
6
θ
=
1
1155
(
2560
P
6
(
cos
θ
)
−
1152
P
4
(
cos
θ
)
−
220
P
2
(
cos
θ
)
−
33
P
0
(
cos
θ
)
)
.
{\displaystyle {\begin{alignedat}{2}T_{0}(\cos \theta )&=1&&=P_{0}(\cos \theta ),\\[4pt]T_{1}(\cos \theta )&=\cos \theta &&=P_{1}(\cos \theta ),\\[4pt]T_{2}(\cos \theta )&=\cos 2\theta &&={\tfrac {1}{3}}{\bigl (}4P_{2}(\cos \theta )-P_{0}(\cos \theta ){\bigr )},\\[4pt]T_{3}(\cos \theta )&=\cos 3\theta &&={\tfrac {1}{5}}{\bigl (}8P_{3}(\cos \theta )-3P_{1}(\cos \theta ){\bigr )},\\[4pt]T_{4}(\cos \theta )&=\cos 4\theta &&={\tfrac {1}{105}}{\bigl (}192P_{4}(\cos \theta )-80P_{2}(\cos \theta )-7P_{0}(\cos \theta ){\bigr )},\\[4pt]T_{5}(\cos \theta )&=\cos 5\theta &&={\tfrac {1}{63}}{\bigl (}128P_{5}(\cos \theta )-56P_{3}(\cos \theta )-9P_{1}(\cos \theta ){\bigr )},\\[4pt]T_{6}(\cos \theta )&=\cos 6\theta &&={\tfrac {1}{1155}}{\bigl (}2560P_{6}(\cos \theta )-1152P_{4}(\cos \theta )-220P_{2}(\cos \theta )-33P_{0}(\cos \theta ){\bigr )}.\end{alignedat}}}
もう一つの性質は、 sin ( n + 1) θ の式であり 、
sin
(
n
+
1
)
θ
sin
θ
=
∑
ℓ
=
0
n
P
ℓ
(
cos
θ
)
P
n
−
ℓ
(
cos
θ
)
.
{\displaystyle {\frac {\sin(n+1)\theta }{\sin \theta }}=\sum _{\ell =0}^{n}P_{\ell }(\cos \theta )P_{n-\ell }(\cos \theta ).}
リカレントニューラルネットワークでは
d 次元のメモリベクトル を含む リカレントニューラルネットワークは 、そのニューラル活動が 次の 状態空間表現によって与えられる 線形時間不変システム に従うように最適化できます 。
m
∈
R
d
{\displaystyle \mathbf {m} \in \mathbb {R} ^{d}}
θ
m
˙
(
t
)
=
A
m
(
t
)
+
B
u
(
t
)
,
{\displaystyle \theta {\dot {\mathbf {m} }}(t)=A\mathbf {m} (t)+Bu(t),}
A
=
[
a
]
i
j
∈
R
d
×
d
,
a
i
j
=
(
2
i
+
1
)
{
−
1
i
<
j
(
−
1
)
i
−
j
+
1
i
≥
j
,
B
=
[
b
]
i
∈
R
d
×
1
,
b
i
=
(
2
i
+
1
)
(
−
1
)
i
.
{\displaystyle {\begin{aligned}A&=\left[a\right]_{ij}\in \mathbb {R} ^{d\times d}{\text{,}}\quad &&a_{ij}=\left(2i+1\right){\begin{cases}-1&i<j\\(-1)^{i-j+1}&i\geq j\end{cases}},\\B&=\left[b\right]_{i}\in \mathbb {R} ^{d\times 1}{\text{,}}\quad &&b_{i}=(2i+1)(-1)^{i}.\end{aligned}}}
この場合、過去の 時間単位にわたる のスライディングウィンドウは、 の時点 における の要素で重み付けされた 最初のシフトされたルジャンドル多項式の線形結合によって 最もよく近似されます 。
u
{\displaystyle u}
θ
{\displaystyle \theta }
d
{\displaystyle d}
m
{\displaystyle \mathbf {m} }
t
{\displaystyle t}
u
(
t
−
θ
′
)
≈
∑
ℓ
=
0
d
−
1
P
~
ℓ
(
θ
′
θ
)
m
ℓ
(
t
)
,
0
≤
θ
′
≤
θ
.
{\displaystyle u(t-\theta ')\approx \sum _{\ell =0}^{d-1}{\widetilde {P}}_{\ell }\left({\frac {\theta '}{\theta }}\right)\,m_{\ell }(t),\quad 0\leq \theta '\leq \theta .}
これらのネットワークは、ディープラーニング 手法と組み合わせることで、より少ない計算リソースで、 長期短期記憶 ユニットや関連アーキテクチャよりも優れた パフォーマンスを発揮するようにトレーニングすることができます。 [5]
追加のプロパティ
ルジャンドル多項式は明確な偶奇性を持つ。つまり、それらは 偶数か奇数か であり、 [6] によれば、
P
n
(
−
x
)
=
(
−
1
)
n
P
n
(
x
)
.
{\displaystyle P_{n}(-x)=(-1)^{n}P_{n}(x)\,.}
もう1つの有用な特性は、
との直交関係を考慮すると となることです 。これは、関数または実験データを近似するためにルジャンドル級数 を使用するときに便利です。 区間 [−1, 1]での級数の 平均は 、単に主要な展開係数 で与えられます 。
∫
−
1
1
P
n
(
x
)
d
x
=
0
for
n
≥
1
,
{\displaystyle \int _{-1}^{1}P_{n}(x)\,dx=0{\text{ for }}n\geq 1,}
P
0
(
x
)
=
1
{\displaystyle P_{0}(x)=1}
∑
i
a
i
P
i
{\textstyle \sum _{i}a_{i}P_{i}}
a
0
{\displaystyle a_{0}}
微分方程式と直交性はスケーリングに依存しないため、ルジャンドル多項式の定義は「標準化」(「正規化」と呼ばれることもありますが、実際のノルムは1ではありません)され、
P
n
(
1
)
=
1
.
{\displaystyle P_{n}(1)=1\,.}
終点における導関数は次のように与えられる。
P
n
′
(
1
)
=
n
(
n
+
1
)
2
.
{\displaystyle P_{n}'(1)={\frac {n(n+1)}{2}}\,.}
ルジャンドル多項式に対するアスケイ・ガスパー不等式は次のように
なる 。
∑
j
=
0
n
P
j
(
x
)
≥
0
for
x
≥
−
1
.
{\displaystyle \sum _{j=0}^{n}P_{j}(x)\geq 0\quad {\text{for }}\quad x\geq -1\,.}
単位ベクトル の スカラー積 のルジャンドル多項式は 、 を用いて
球面調和関数 で展開することができます。
ここで、単位ベクトル r と r ′ は、それぞれ 球面座標 ( θ , φ ) と ( θ ′, φ ′) を持ちます 。
P
ℓ
(
r
⋅
r
′
)
=
4
π
2
ℓ
+
1
∑
m
=
−
ℓ
ℓ
Y
ℓ
m
(
θ
,
φ
)
Y
ℓ
m
∗
(
θ
′
,
φ
′
)
,
{\displaystyle P_{\ell }\left(r\cdot r'\right)={\frac {4\pi }{2\ell +1}}\sum _{m=-\ell }^{\ell }Y_{\ell m}(\theta ,\varphi )Y_{\ell m}^{*}(\theta ',\varphi ')\,,}
2つのルジャンドル多項式[7]
の積。
ここでは 第一種完全楕円積分 である 。
∑
p
=
0
∞
t
p
P
p
(
cos
θ
1
)
P
p
(
cos
θ
2
)
=
2
π
K
(
2
t
sin
θ
1
sin
θ
2
t
2
−
2
t
cos
(
θ
1
+
θ
2
)
+
1
)
t
2
−
2
t
cos
(
θ
1
+
θ
2
)
+
1
,
{\displaystyle \sum _{p=0}^{\infty }t^{p}P_{p}(\cos \theta _{1})P_{p}(\cos \theta _{2})={\frac {2}{\pi }}{\frac {\mathbf {K} \left(2{\sqrt {\frac {t\sin \theta _{1}\sin \theta _{2}}{t^{2}-2t\cos \left(\theta _{1}+\theta _{2}\right)+1}}}\right)}{\sqrt {t^{2}-2t\cos \left(\theta _{1}+\theta _{2}\right)+1}}}\,,}
K
(
⋅
)
{\displaystyle K(\cdot )}
再帰関係
上で述べたように、ルジャンドル多項式は、ボネの再帰式として知られる3項再帰関係に従う。
これは、
および
または、別の表現で、端点でも成り立つ。
(
n
+
1
)
P
n
+
1
(
x
)
=
(
2
n
+
1
)
x
P
n
(
x
)
−
n
P
n
−
1
(
x
)
{\displaystyle (n+1)P_{n+1}(x)=(2n+1)xP_{n}(x)-nP_{n-1}(x)}
x
2
−
1
n
d
d
x
P
n
(
x
)
=
x
P
n
(
x
)
−
P
n
−
1
(
x
)
{\displaystyle {\frac {x^{2}-1}{n}}{\frac {d}{dx}}P_{n}(x)=xP_{n}(x)-P_{n-1}(x)}
d
d
x
P
n
+
1
(
x
)
=
(
n
+
1
)
P
n
(
x
)
+
x
d
d
x
P
n
(
x
)
.
{\displaystyle {\frac {d}{dx}}P_{n+1}(x)=(n+1)P_{n}(x)+x{\frac {d}{dx}}P_{n}(x)\,.}
ルジャンドル多項式の積分に便利なのは
(
2
n
+
1
)
P
n
(
x
)
=
d
d
x
(
P
n
+
1
(
x
)
−
P
n
−
1
(
x
)
)
.
{\displaystyle (2n+1)P_{n}(x)={\frac {d}{dx}}{\bigl (}P_{n+1}(x)-P_{n-1}(x){\bigr )}\,.}
上記から、または
同等に、
‖ P n ‖は区間 −1 ≤ x ≤ 1 上のノルムである
ことがわかる。
d
d
x
P
n
+
1
(
x
)
=
(
2
n
+
1
)
P
n
(
x
)
+
(
2
(
n
−
2
)
+
1
)
P
n
−
2
(
x
)
+
(
2
(
n
−
4
)
+
1
)
P
n
−
4
(
x
)
+
⋯
{\displaystyle {\frac {d}{dx}}P_{n+1}(x)=(2n+1)P_{n}(x)+{\bigl (}2(n-2)+1{\bigr )}P_{n-2}(x)+{\bigl (}2(n-4)+1{\bigr )}P_{n-4}(x)+\cdots }
d
d
x
P
n
+
1
(
x
)
=
2
P
n
(
x
)
‖
P
n
‖
2
+
2
P
n
−
2
(
x
)
‖
P
n
−
2
‖
2
+
⋯
{\displaystyle {\frac {d}{dx}}P_{n+1}(x)={\frac {2P_{n}(x)}{\left\|P_{n}\right\|^{2}}}+{\frac {2P_{n-2}(x)}{\left\|P_{n-2}\right\|^{2}}}+\cdots }
‖
P
n
‖
=
∫
−
1
1
(
P
n
(
x
)
)
2
d
x
=
2
2
n
+
1
.
{\displaystyle \|P_{n}\|={\sqrt {\int _{-1}^{1}{\bigl (}P_{n}(x){\bigr )}^{2}\,dx}}={\sqrt {\frac {2}{2n+1}}}\,.}
漸近解析
漸近的には、に対して 、ルジャンドル多項式は [8]および 1
より大きい引数に対して [9] と表すことができます。
ここで、 J0 、 J1 、 および I0 は ベッセル関数 です 。
ℓ
→
∞
{\displaystyle \ell \to \infty }
P
ℓ
(
cos
θ
)
=
θ
sin
(
θ
)
{
J
0
[
(
ℓ
+
1
2
)
θ
]
−
(
1
θ
−
cot
θ
)
8
(
ℓ
+
1
2
)
J
1
[
(
ℓ
+
1
2
)
θ
]
}
+
O
(
ℓ
−
2
)
=
2
π
ℓ
sin
(
θ
)
cos
[
(
ℓ
+
1
2
)
θ
−
π
4
]
+
O
(
ℓ
−
3
/
2
)
,
θ
∈
(
0
,
π
)
,
{\displaystyle {\begin{aligned}P_{\ell }(\cos \theta )&={\sqrt {\frac {\theta }{\sin \left(\theta \right)}}}\left\{J_{0}{\left[\left(\ell +{\tfrac {1}{2}}\right)\theta \right]}-{\frac {\left({\frac {1}{\theta }}-\cot \theta \right)}{8(\ell +{\frac {1}{2}})}}J_{1}{\left[\left(\ell +{\tfrac {1}{2}}\right)\theta \right]}\right\}+{\mathcal {O}}\left(\ell ^{-2}\right)\\[1ex]&={\sqrt {\frac {2}{\pi \ell \sin \left(\theta \right)}}}\cos \left[\left(\ell +{\tfrac {1}{2}}\right)\theta -{\tfrac {\pi }{4}}\right]+{\mathcal {O}}\left(\ell ^{-3/2}\right),\quad \theta \in (0,\pi ),\end{aligned}}}
P
ℓ
(
cosh
ξ
)
=
ξ
sinh
ξ
I
0
(
(
ℓ
+
1
2
)
ξ
)
(
1
+
O
(
ℓ
−
1
)
)
,
P
ℓ
(
1
1
−
e
2
)
=
1
2
π
ℓ
e
(
1
+
e
)
ℓ
+
1
2
(
1
−
e
)
ℓ
2
+
O
(
ℓ
−
1
)
{\displaystyle {\begin{aligned}P_{\ell }\left(\cosh \xi \right)&={\sqrt {\frac {\xi }{\sinh \xi }}}I_{0}\left(\left(\ell +{\frac {1}{2}}\right)\xi \right)\left(1+{\mathcal {O}}\left(\ell ^{-1}\right)\right)\,,\\P_{\ell }\left({\frac {1}{\sqrt {1-e^{2}}}}\right)&={\frac {1}{\sqrt {2\pi \ell e}}}{\frac {(1+e)^{\frac {\ell +1}{2}}}{(1-e)^{\frac {\ell }{2}}}}+{\mathcal {O}}\left(\ell ^{-1}\right)\end{aligned}}}
ゼロ
のすべての 零点 は実数で、互いに異なり、区間 内にあります 。さらに、これらを区間を部分区間 に分割するものと見なすと 、各部分区間には の零点が正確に 1 つ含まれます。これはインターレース特性として知られています。パリティ特性により、 が の零点であれば も の零点 である ことは明らかです。これらの零点は 、ガウス積分法 に基づく数値積分で重要な役割を果たします 。 に基づく特定の積分法は、 ガウス・ルジャンドル積分法 として知られています 。
n
{\displaystyle n}
P
n
(
x
)
{\displaystyle P_{n}(x)}
(
−
1
,
1
)
{\displaystyle (-1,1)}
[
−
1
,
1
]
{\displaystyle [-1,1]}
n
+
1
{\displaystyle n+1}
P
n
+
1
{\displaystyle P_{n+1}}
x
k
{\displaystyle x_{k}}
P
n
(
x
)
{\displaystyle P_{n}(x)}
−
x
k
{\displaystyle -x_{k}}
P
n
{\displaystyle P_{n}}
この性質と という事実から、 は において極小値と極大値を 持つ ことがわかります 。同様に、 は においてゼロを 持ちます 。
P
n
(
±
1
)
≠
0
{\displaystyle P_{n}(\pm 1)\neq 0}
P
n
(
x
)
{\displaystyle P_{n}(x)}
n
−
1
{\displaystyle n-1}
(
−
1
,
1
)
{\displaystyle (-1,1)}
d
P
n
(
x
)
/
d
x
{\displaystyle dP_{n}(x)/dx}
n
−
1
{\displaystyle n-1}
(
−
1
,
1
)
{\displaystyle (-1,1)}
ポイントごとの評価
パリティと正規化により、境界における値は 次のように
表される。
原点においては 、値は次のように表される。
x
=
±
1
{\displaystyle x=\pm 1}
P
n
(
1
)
=
1
,
P
n
(
−
1
)
=
(
−
1
)
n
{\displaystyle P_{n}(1)=1\,,\quad P_{n}(-1)=(-1)^{n}}
x
=
0
{\displaystyle x=0}
P
2
n
(
0
)
=
(
−
1
)
n
4
n
(
2
n
n
)
=
(
−
1
)
n
2
2
n
(
2
n
)
!
(
n
!
)
2
=
(
−
1
)
n
(
2
n
−
1
)
!
!
(
2
n
)
!
!
{\displaystyle P_{2n}(0)={\frac {(-1)^{n}}{4^{n}}}{\binom {2n}{n}}={\frac {(-1)^{n}}{2^{2n}}}{\frac {(2n)!}{\left(n!\right)^{2}}}=(-1)^{n}{\frac {(2n-1)!!}{(2n)!!}}}
P
2
n
+
1
(
0
)
=
0
{\displaystyle P_{2n+1}(0)=0}
シフトルジャンドル多項式
シフトされたルジャンドル多項式は 次 のように定義されます。
ここで、「シフト」関数 x ↦ 2 x − 1 は、 区間 [0, 1] を区間 [−1, 1]に 全単射にマッピングする アフィン変換 であり 、多項式 P̃ n ( x ) が [0, 1] 上で直交することを意味します。
P
~
n
(
x
)
=
P
n
(
2
x
−
1
)
.
{\displaystyle {\widetilde {P}}_{n}(x)=P_{n}(2x-1)\,.}
∫
0
1
P
~
m
(
x
)
P
~
n
(
x
)
d
x
=
1
2
n
+
1
δ
m
n
.
{\displaystyle \int _{0}^{1}{\widetilde {P}}_{m}(x){\widetilde {P}}_{n}(x)\,dx={\frac {1}{2n+1}}\delta _{mn}\,.}
シフトされたルジャンドル多項式の明示的な表現は次のように与えられる。
P
~
n
(
x
)
=
(
−
1
)
n
∑
k
=
0
n
(
n
k
)
(
n
+
k
k
)
(
−
x
)
k
.
{\displaystyle {\widetilde {P}}_{n}(x)=(-1)^{n}\sum _{k=0}^{n}{\binom {n}{k}}{\binom {n+k}{k}}(-x)^{k}\,.}
シフトされたルジャンドル多項式に対する
ロドリゲスの公式 の類似は
P
~
n
(
x
)
=
1
n
!
d
n
d
x
n
(
x
2
−
x
)
n
.
{\displaystyle {\widetilde {P}}_{n}(x)={\frac {1}{n!}}{\frac {d^{n}}{dx^{n}}}\left(x^{2}-x\right)^{n}\,.}
最初のいくつかのシフトされたルジャンドル多項式は次のとおりです。
ルジャンドル有理関数
ルジャンドル 有理関数は、[0, ∞) 上の 直交関数 の列です。ルジャンドル有理関数は ケーリー変換 とルジャンドル多項式を組み合わせることで得られます 。
n 次の有理ルジャンドル関数は 次のように定義されます。
R
n
(
x
)
=
2
x
+
1
P
n
(
x
−
1
x
+
1
)
.
{\displaystyle R_{n}(x)={\frac {\sqrt {2}}{x+1}}\,P_{n}\left({\frac {x-1}{x+1}}\right)\,.}
これらは、 特異 シュトゥルム・リウヴィル問題の 固有関数 であり、
固有値は
(
x
+
1
)
d
d
x
(
x
d
d
x
[
(
x
+
1
)
v
(
x
)
]
)
+
λ
v
(
x
)
=
0
{\displaystyle \left(x+1\right){\frac {d}{dx}}\left(x{\frac {d}{dx}}\left[\left(x+1\right)v(x)\right]\right)+\lambda v(x)=0}
λ
n
=
n
(
n
+
1
)
.
{\displaystyle \lambda _{n}=n(n+1)\,.}
参照
注記
^ アルフケン&ウェーバー 2005、p.743
^ ボアス、メアリー L. (2006)。 物理科学における数学的手法 (第 3 版)。ホーボーケン、ニュージャージー: ワイリー 。ISBN 978-0-471-19826-0 。
^ ルジャンドル、A.-M. (1785) [1782]。 「球状物の均質性の魅力に関する研究」 (PDF) 。 Mémoires de Mathématiques et de Physique、présentés à l'Académie Royale des Sciences、par divers savans、et lus dans ses Assemblées (フランス語)。 Vol. X.パリ。 411–435ページ。 2009 年 9 月 20 日の オリジナル (PDF) からアーカイブされました。
^ ジャクソン、JD (1999)。 古典電気力学 (第3版)。ワイリー&サンズ。p. 103。ISBN 978-0-471-30932-1 。 {{cite book}}: CS1 maint: location missing publisher (link)
^ Voelker, Aaron R.; Kajić, Ivana; Eliasmith, Chris (2019). Legendre メモリ単位: リカレントニューラルネットワークにおける連続時間表現 (PDF) 。ニューラル情報処理システムの進歩。
^ アルフケン&ウェーバー 2005、p.753
^ レナード・C・マキシモン (1957)。 「2 つのルジャンドル多項式の積の生成関数」。 Norske Videnskabers Selskab Forhandlinger 。 29 :82~86。
^ Szegő, Gábor (1975). 直交多項式 (第4版). プロビデンス:アメリカ数学会. pp. 194(定理8.21.2) 。ISBN 0821810235 . OCLC 1683237.
^ 「DLMF: 14.15 一様漸近近似」。
参考文献
アブラモウィッツ、ミルトン 、 ステグン、アイリーン ・アン編 (1983) [1964 年 6 月]。「第 8 章」。 数式、グラフ、および数学表付き数学関数ハンドブック 。応用数学シリーズ。第 55 巻 (第 10 刷の訂正を加えた第 9 刷 (1972 年 12 月)、初版)。ワシントン DC、ニューヨーク: 米国商務省、国立標準局、ドーバー出版。332、773 ページ 。ISBN 978-0-486-61272-0 LCCN 64-60036。MR 0167642。LCCN 65-12253 。 第22章も参照してください。
Arfken, George B. ; Weber, Hans J. (2005). 物理学者のための数学的手法 . Elsevier Academic Press. ISBN 0-12-059876-0 。
Bayin, SS (2006). 科学と工学における数学的手法 . Wiley. 第 2 章. ISBN 978-0-470-04142-0 。
Belousov, SL (1962)。 正規化準ルジャンドル多項式の表 。数学表。第 18 巻。Pergamon Press。ISBN 978-0-08-009723-7 。
クーラント、リチャード ; ヒルベルト、デイヴィッド (1953)。 数理物理学の方法 。第 1 巻。ニューヨーク、NY: インターサイエンス 。ISBN 978-0-471-50447-4 。
Dunster, TM (2010)、「Legendre と関連関数」、 Olver, Frank WJ 、Lozier, Daniel M.、Boisvert, Ronald F.、Clark, Charles W. (編)、 NIST Handbook of Mathematical Functions 、Cambridge University Press、 ISBN 978-0-521-19225-5 、 MR 2723248 。
El Attar, Refaat (2009)。Legendre 多項式 と関数 。CreateSpace。ISBN 978-1-4414-9012-4 。
Koornwinder, Tom H. ; Wong, Roderick SC; Koekoek, Roelof; Swarttouw, René F. (2010)、「直交多項式」、 Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.)、 NIST Handbook of Mathematical Functions 、Cambridge University Press、 ISBN 978-0-521-19225-5 、 MR 2723248 。
外部リンク
ウィキメディア・コモンズには、ルジャンドル多項式 に関連するメディアがあります 。
水素の量子力学の文脈におけるルジャンドル多項式の簡単な非公式な導出
「ルジャンドル多項式」、 数学百科事典 、 EMS Press 、2001 [1994]
Wolfram MathWorld のルジャンドル多項式に関するエントリ
ジェームズ・B・カルバート博士の個人的数学コレクションからのルジャンドル多項式に関する論文
カーライル・E・ムーア著『ルジャンドル多項式』
超物理学からのルジャンドル多項式