勒壤得多項式

數學上，勒壤得函數指以下勒壤得微分方程式的解：

${\displaystyle (1-x^{2}){\frac {\mathrm {d} ^{2}P(x)}{\mathrm {d} x^{2}}}-2x{\frac {\mathrm {d} P(x)}{\mathrm {d} x}}+n(n+1)P(x)=0.}$

為求解方便一般也寫成如下史特姆-萊歐維爾形式:

${\displaystyle {\mathrm {d} \over \mathrm {d} x}\left[(1-x^{2}){\mathrm {d} \over \mathrm {d} x}P(x)\right]+n(n+1)P(x)=0.}$

上述方程式及其解函數因法國數學家阿德里安-馬里·勒壤得而得名。勒壤得方程式是物理學和其他技術領域常常遇到的一類常微分方程式。當試圖在球坐標中求解三維拉普拉斯方程式（或相關的其他偏微分方程式）時，問題便會歸結為勒讓德方程式的求解。

勒壤得方程式的解可寫成標準的冪級數形式。當方程式滿足 ${\displaystyle |x|<1}$ 時，可得到有界解（即解級數收斂）。並且當n 為非負整數，即${\displaystyle n=0,1,2,\ldots }$.

正交性

勒壤得多項式的一個重要性質是其在區間 ${\displaystyle -1\leq x\leq 1}$ 關於L²內積滿足正交性，即：

${\displaystyle \int _{-1}^{1}P_{m}(x)P_{n}(x)\,\mathrm {d} x={2 \over {2n+1}}\delta _{mn}}$

其中 ${\displaystyle \delta _{mn}}$ 為克羅內克δ記號，當${\displaystyle m=n}$ 時為1，否則為0。 事實上，推導勒壤得多項式的另一種方法便是關於前述內積空間對多項式 ${\displaystyle {1,x,x^{2},\ldots }}$ 進行格拉姆-施密特正交化。之所以具有此正交性是因為如前所述，勒壤得微分方程式可化為標準的Sturm-Liouville問題：

${\displaystyle {\mathrm {d} \over \mathrm {d} x}\left[(1-x^{2}){\mathrm {d} \over \mathrm {d} x}P(x)\right]=-\lambda P(x),}$

其中本徵值 ${\displaystyle \lambda }$ 對應於原方程式中的 ${\displaystyle n(n+1)}$。

部分實例

下表列出了前11階（n 從0到10）勒壤得多項式的表達式：

n	${\displaystyle P_{n}(x)\,}$
0	${\displaystyle 1\,}$
1	${\displaystyle x\,}$
2	${\displaystyle {\begin{matrix}{\frac {1}{2}}\end{matrix}}(3x^{2}-1)\,}$
3	${\displaystyle {\begin{matrix}{\frac {1}{2}}\end{matrix}}(5x^{3}-3x)\,}$
4	${\displaystyle {\begin{matrix}{\frac {1}{8}}\end{matrix}}(35x^{4}-30x^{2}+3)\,}$
5	${\displaystyle {\begin{matrix}{\frac {1}{8}}\end{matrix}}(63x^{5}-70x^{3}+15x)\,}$
6	${\displaystyle {\begin{matrix}{\frac {1}{16}}\end{matrix}}(231x^{6}-315x^{4}+105x^{2}-5)\,}$
7	${\displaystyle {\begin{matrix}{\frac {1}{16}}\end{matrix}}(429x^{7}-693x^{5}+315x^{3}-35x)\,}$
8	${\displaystyle {\begin{matrix}{\frac {1}{128}}\end{matrix}}(6435x^{8}-12012x^{6}+6930x^{4}-1260x^{2}+35)\,}$
9	${\displaystyle {\begin{matrix}{\frac {1}{128}}\end{matrix}}(12155x^{9}-25740x^{7}+18018x^{5}-4620x^{3}+315x)\,}$
10	${\displaystyle {\begin{matrix}{\frac {1}{256}}\end{matrix}}(46189x^{10}-109395x^{8}+90090x^{6}-30030x^{4}+3465x^{2}-63)\,}$

前6階（n 從0到5）勒壤得多項式的曲線如下圖所示：

在物理學中的應用

在求解三維空間中的球對稱問題，譬如計算點電荷在空間中激發的電勢時，常常要用到勒壤得多項式作如下形式的級數展開：

${\displaystyle {\frac {1}{\left|\mathbf {x} -\mathbf {x} ^{\prime }\right|}}={\frac {1}{\sqrt {r^{2}+r^{\prime 2}-2rr'\cos \gamma }}}=\sum _{\ell =0}^{\infty }{\frac {r^{\prime \ell }}{r^{\ell +1}}}P_{\ell }(\cos \gamma )}$

其中${\displaystyle r}$和${\displaystyle r'}$分別為位置向量${\displaystyle \mathbf {x} }$和${\displaystyle \mathbf {x} ^{\prime }}$ 的長度(其中${\displaystyle r}$和${\displaystyle r'}$分別為對位置向量${\displaystyle \mathbf {x} }$和${\displaystyle \mathbf {x} ^{\prime }}$ 的長度進行測量的結果)，${\displaystyle \gamma }$為兩向量的夾角(${\displaystyle \gamma }$為對兩向量的夾角展開估計的結果)。當${\displaystyle r>r'}$時上式成立。該式計算了在${\displaystyle \mathbf {x} '}$處的點電荷激發的電場在${\displaystyle \mathbf {x} }$點引起的電勢大小。在對空間中連續分佈的電荷引起的電勢大小進行計算時(當計算由連續分佈之電荷所產生的電勢時)，將涉及對上式進行積分(需積分上式中間項)。這時，上式右邊的勒壤得多項式展開將對此積分的計算帶來很大的方便(逐項積分上式右邊的展開式可得一級數解，此級數之第一項叫做電單極矩，第二項叫做電偶極矩，第三項叫做電四極矩)。

靜電場中具有軸對稱邊界條件的問題可以歸結為在球坐標系中用分離變量法求解關於電勢函數的拉普拉斯方程式${\displaystyle \nabla ^{2}\Phi (\mathbf {x} )=0}$（與和對稱軸的夾角無關）。若設${\displaystyle {\widehat {\mathbf {z} }}}$為對稱軸，${\displaystyle \theta }$為觀測者位置向量和${\displaystyle {\widehat {\mathbf {z} }}}$軸的夾角，則勢函數的解可表示為：

${\displaystyle \Phi (r,\theta )=\sum _{\ell =0}^{\infty }\left[A_{\ell }r^{\ell }+B_{\ell }r^{-(\ell +1)}\right]P_{\ell }(\cos \theta ).}$

其中${\displaystyle A_{\ell }}$和${\displaystyle B_{\ell }}$由具體邊界條件確定^[1]。

其他性質

勒壤得多項式的奇偶性由其階數確定。當階數k為偶數時，${\displaystyle P_{k}(x)}$為偶函數；當階數k為奇數時，${\displaystyle P_{k}(x)}$為奇函數，即：

${\displaystyle P_{k}(-x)=(-1)^{k}P_{k}(x).\,}$

遞歸關係

相鄰的三個勒壤得多項式具有三項遞歸關係：

${\displaystyle (n+1)P_{n+1}=(2n+1)xP_{n}-nP_{n-1}\,}$

另外，考慮微分後還有以下遞歸關係：

${\displaystyle {x^{2}-1 \over n}{\mathrm {d} \over \mathrm {d} x}P_{n}=xP_{n}-P_{n-1}.}$

${\displaystyle (2n+1)P_{n}={\mathrm {d} \over \mathrm {d} x}\left[P_{n+1}-P_{n-1}\right].}$

其中最後一個式子在計算勒壤得多項式的積分中較為有用。

使多項式的值：

#include <iostream>
using namespace std;

int main()
{
	float n,x;
	float polyaendl;

	return 0;
}

float polya(float n, float x)
{
	if (n == 0) return 1.0;
	eurn x;
	else return ((2.0 * n - 1.0) * x * polya(n - 1.0, x) - (n - 1.0) * polya(n - 2.0, x)) / n;
}

移位勒壤得多項式

移位勒壤得多項式${\displaystyle {\tilde {P_{n}}}(x)}$的正交區間定義在${\displaystyle [0,1]}$上，即：

${\displaystyle \int _{0}^{1}{\tilde {P_{m}}}(x){\tilde {P_{n}}}(x)\,\mathrm {d} x={1 \over {2n+1}}\delta _{mn}.}$

其顯式表達式為：

${\displaystyle {\tilde {P_{n}}}(x)=(-1)^{n}\sum _{k=0}^{n}{n \choose k}{n+k \choose k}(-x)^{k}.}$

相應的羅德里格公式為：

${\displaystyle {\tilde {P_{n}}}(x)=(n!)^{-1}{\mathrm {d} ^{n} \over \mathrm {d} x^{n}}\left[(x^{2}-x)^{n}\right].\,}$

下表列出了前4階移位勒壤得多項式：

n	${\displaystyle {\tilde {P_{n}}}(x)}$
0	1
1	${\displaystyle 2x-1}$
2	${\displaystyle 6x^{2}-6x+1}$
3	${\displaystyle 20x^{3}-30x^{2}+12x-1}$

分數階勒壤得多項式

分數階勒壤得多項式通過將分數階微分和通過Γ函數定義的非整數階乘代入羅德里格公式中來定義。

極限關係

大Q勒壤得多項式→勒壤得多項式

令大q雅可比多項式中的${\displaystyle c=0}$,即勒壤得多項式

令連續q勒壤得多項式 q->1得勒壤得多項式

${\displaystyle \lim _{q\to 1}P_{n}(x|q)=P_{n}(x)}$

小q勒壤得多項式→勒壤得多項式

${\displaystyle \lim _{q\to 1}p_{n}(x|q)=P_{n}(1-2x)}$

參見

• 高斯求積
• 伴隨勒壤得多項式
• 勒壤得有理函數

外部連結

• 沃爾夫勒姆（Wolfram）數學世界對勒壤得多項式的介紹（英文） （頁面存檔備份，存於互聯網檔案館）

參考文獻

1. 嚴鎮軍編，《數學物理方程式》，第二版，中國科學技術大學出版社，合肥，2002，ISBN 7-312-00799-6，第140頁

• 2. Milton Abramowitz and Irene A. Stegun, eds. (1965). Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. New York: Dover. ISBN 0-486-61272-4（參見 第8章 （頁面存檔備份，存於互聯網檔案館）和第22章 （頁面存檔備份，存於互聯網檔案館））