5

Legendre多項式の性質

821
0

Legendre多項式の性質

以前の記事の続きになります。

Legendre多項式はStrum=Liouville理論や超幾何関数の1本の枝葉として理解することができます。

一般論を展開する準備として、まずは具体について把握しておくと、抽象化の過程を10倍楽しむことができます。

Rodorigueの公式

Legendre多項式はRodorigue型多項式の1種です。
Hermiete多項式、Laguerre多項式、Gegenbauer多項式、Sonin多項式などもその一種です。

Rodorigueの公式

Pn(z)=12nn!dndzn(z21)n

モニックな2n次式をn階微分しているので、Pn(z)n次多項式を定義します。

母関数

Rodorigueの公式の右辺をマクローリン展開の展開係数とみなせば母関数表示ができます。閉じた形を作ろうと思った際には、複素関数の金字塔であるコーシーの積分定理を使うと感動します。

Legendre多項式の母関数表示

G(t,z)=112tz+t2=n=0Pn(z)tn

証明手法(任意)

コーシーの積分定理より、正則関数f(z)に対して、
12πiCf(s)szds=f(z)
より、12πiC(s21)nszds=(z21)n
両辺zn回微分すると、
Pn(z)=12nn!dndzn(z21)n=12nn!12πiCdndzn[(s21)nsz]ds
=12πiC[(s21)n2n(sz)n+1]ds
積分経路はz中心の円です。母関数G(t,z)は、
G(t,z)=n=0tnPn(z)=12πin=0tnC[(s21)n2n+1(sz)n]ds
=12πiCn=0[tn(s21)n2n(sz)n+1]ds
|t|<|2(sz)(s21)| となれば被積分関数は一様収束し、収束円内で正則関数を定義します。故に上記の無限和と積分の交換を実行可能になります。

G(t,z)=12πiCn=0[tn(s21)n2n(sz)n+1]ds
=12πiC1sz11t(s21)2(sz)dz=12πiC2t2z2sts2ds
sに関する二次方程式t2z2sts2=0の解s±=1t(1±12zt+t2)は、t0の極限でs+, s0となります。tが十分小さい時にはsのみが積分経路の内部に入るため、sの留数のみを計算すれば良いことがわかります。

G(t,z)=12πiC2t(ss+)(ss)ds=2t(ss+)=112zt+t2

正則関数列、一様収束性は非常に相性が良いです。正則関数列が(広義)一様収束すれば、積分と和を入れ替えられるのみならず、収束した関数まで正則関数になることがMoreraの定理より証明できます。
べき関数は収束円内で一様収束する正則関数列を定義することがM判定法からわかるので、収束円内で正則関数に収束することが自動的に理解できます。

母関数は係数数列から新しい情報をもたらしてくれます。それが以下に紹介する漸化式です。

漸化式

Legendre多項式の漸化式

(1):(2n+1)zPn(z)=(n+1)Pn+1(z)+nPn1(z)
(2):ddzPn+1(z)+ddzPn1(z)=2zddzPn(z)+Pn(z)

母関数を微分して係数を比較します。

(1)Pn(z)の微分を作りたくないので、zではなくtで微分します。
ddtG(t,z)=zt(12tz+t2)32=zt12tz+t2G(t,z)=n=1nPn(z)tn1
(12tz+t2)n=1nPn(z)tn1=(zt)G(t,z)=(zt)n=0Pn(z)tn
展開して整理すると、
n=1[(2n+1)zPn(z)(n+1)Pn+1(z)+nPn1(z)]tn+P1(z)zP0(z)=0

係数が恒等的に0になる条件は
(2n+1)zPn(z)=(n+1)Pn+1(z)+nPn1(z)
(2)Pn(z)の微分を作るため、zで微分します。
ddzG(t,z)=t(12tz+t2)32=t12tz+t2G(t,z)=n=0Pn(z)tn
両辺分母を払って同様に計算する。0次多項式P0(z)=0に注意すると、
P1(z)t+n=2[Pn(z)2zPn1(z)+Pn2(z)Pn1(z)]t2=0
係数ゼロの条件から添字をずらしてあげると、
ddzPn+1(z)+ddzPn1(z)=2zddzPn(z)+Pn(z)

直交性

Legendre多項式の直交性

11Pm(z)Pn(z)dz=22n+1δmn

定理4の証明

mnの場合を考える。被積分関数はm,nに関して対称なので、m>nとしても一般性を失わない。

11Pm(z)Pn(z)dz=12mm!11dmdzm(z21)mPn(z)dz
m回部分積分を施すと、表面項には(z21)kがかかるので、表面項は全て消える。

=(1)m2mm!11(z21)mdmdzm(Pn(z))dz
Pn(z)n次多項式なのでm回微分すると消える。故に積分値はゼロになるため、
11Pm(z)Pn(z)dz=0

m=nの場合、Pn(z)n次の係数は(2n)!/2nn!2なので、

11Pn(z)Pn(z)dz=(1)n2nn!11(z21)ndndzn(Pn(z))dz

=12nn!11(1z)n(z+1)n(2n)!2nn!2n!dz

ベータ関数の積分公式より、

11(1z)n(z+1)ndz=22n+101un(1u)ndu

=22n+1B(n+1,n+1)=22n+1n!n!(2n+1)!

11Pn(z)Pn(z)dz=12nn!(2n)!2nn!22n+1n!n!(2n+1)!

=22n+1

従って、直交性が示された。

Legendreの微分方程式

母関数から逆に、Pn(z)の満たす微分方程式を導出することもできます。

Legendreの微分方程式

Legendre多項式Pn(z)は以下の微分方程式を満たす。
(1z2)d2Pndz22zdPndz+n(n+1)Pn=0

母関数にz微分、t微分を適用する。

計算は煩雑だが、G(t,z)に微分作用素(1z2)ddz2zddz を施すと、
(1z2)dGdz2zdGdz=t(3t(1z2)2z(z22tz+1))(z22tz+1)52
n(n+1)を出すために、G(t,z)tをかけてから2回微分すると、
d2dt2(tG(t,z))=t(3t(1z2)2z(z22tz+1))(z22tz+1)52
となる。従って
(1z2)dGdz2zdGdz+d2dt2(tG(t,z))=0
級数で表示して、
n=0[(1z2)d2Pndz22zdPndz+n(n+1)Pn]tn=0
係数がゼロになる条件から、微分方程式が導かれた。

投稿日:2023109
更新日:2024114
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

noho1024
noho1024
22
9599
専攻は物理学です。 有機化学、物性理論、確率数理、場の量子論の勉強を経て、科学の面白さを世に広める活動をしていきたいと思っています。

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. Legendre多項式の性質
  2. Rodorigueの公式
  3. 母関数
  4. 漸化式
  5. 直交性
  6. Legendreの微分方程式