勒让德多项式

维基百科,自由的百科全书
跳转至: 导航搜索

数学上,勒让德函数指以下勒让德微分方程的解:

(1 - x^2 )\frac{{d^2 P(x)}}{{dx^2 }} - 2x\frac{{dP(x)}}{{dx}} + n(n + 1)P(x) = 0.

为求解方便一般也写成如下施图姆-刘维尔形式Sturm-Liouville form):

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

上述方程及其解函数因法国数学家阿德里安-马里·勒让德而得名。勒让德方程是物理学和其他技术领域常常遇到的一类常微分方程。当试图在球坐标中求解三维拉普拉斯方程(或相关的其他偏微分方程)时,问题便会归结为勒让德方程的求解。

勒让德方程的解可写成标准的幂级数形式。当方程满足 |x| < 1 时,可得到有界解(即解级数收敛)。并且当n 为非负整数,即n = 0, 1, 2,... 时,在x = ± 1 点亦有有界解。这种情况下,随n 值变化方程的解相应变化,构成一组由正交多项式组成的多项式序列,这组多项式称为勒让德多项式(Legendre polynomials)。

勒让德多项式Pn(x)是n 阶多项式,可用罗德里格公式表示为:

P_n(x) = {1 \over 2^n n!} {d^n \over dx^n } \left[ (x^2 -1)^n \right].

正交性[编辑]

勒让德多项式的一个重要性质是其在区间 −1 ≤ x ≤ 1 关于L2内积满足正交性,即:

\int_{-1}^{1} P_m(x) P_n(x)\,dx = {2 \over {2n + 1}} \delta_{mn}

其中 δmn克罗内克δ记号,当m = n 时为1,否则为0。 事实上,推导勒让德多项式的另一种方法便是关于前述内积空间对多项式{1, x, x2, ...}进行格拉姆-施密特正交化。之所以具有此正交性是因为如前所述,勒让德微分方程可化为标准的strum-liouville问题

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

其中本征值 λ 对应于原方程中的 n(n+1)。

部分实例[编辑]

下表列出了头11阶(n 从0到10)勒让德多项式的表达式:

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

头6阶(n 从0到5)勒让德多项式的曲线如下图所示:

Legendre poly.svg

在物理学中的应用[编辑]

在求解三维空间中的球对称问题,譬如计算点电荷在空间中激发的电势时,常常要用到勒让德多项式作如下形式的级数展开:

\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)

其中rr'分别为位置向量\mathbf{x}\mathbf{x}^\prime 的长度,\gamma为两向量的夹角。当r>r'时上式成立。该式计算了在\mathbf{x}'处的点电荷激发的电场\mathbf{x}点引起的电势大小。在对空间中连续分布的电荷引起的电势大小进行计算时,将涉及对上式进行积分。这时,上式右边的勒让德多项式展开将对此积分的计算带来很大的方便。

静电场中具有轴对称边界条件的问题可以归结为在球坐标系中用分离变量法求解关于电势函数的拉普拉斯方程\nabla^2 \Phi(\mathbf{x})=0(与和对称轴的夹角无关)。若设\widehat{\mathbf{z}}为对称轴,\theta为观测者位置向量和\widehat{\mathbf{z}}轴的夹角,则势函数的解可表示为:


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

其中A_\ellB_\ell由具体边界条件确定[1]


其他性质[编辑]

奇偶性[编辑]

当阶数k偶数时,P_k(x)偶函数;当阶数k奇数时,P_k(x)奇函数,即:

P_k(-x) = (-1)^k P_k(x). \,

递推关系[编辑]

相邻的三个勒让德多项式具有三项递推关系式:

 (n+1) P_{n+1} = (2n+1) x P_n - n P_{n-1}\,

另外,考虑微分后还有以下递推关系:

{x^2-1 \over n} {d \over dx} P_n = xP_n - P_{n-1}.
(2n+1) P_n = {d \over dx} \left[ P_{n+1} - P_{n-1} \right].

其中最后一个式子在计算勒让德多项式的积分中较为有用。

使用C++语言,利用递归方法求n阶勒让德多项式的值:

#include<iostream>

using namespace std;

int main()

{

float n,x;


float polya(float ,float );

cout<<"please input x and n:";

cin>>x>>n; //输入时加上小数位 如: n为1 输入:1.0


cout<<polya(n,x)<<endl;

return 0;

}


float polya(float n,float x)

{

if(n==0) return 1.0;

else if(1==n) return x;

else return ( (2.0*n-1.0)*x*polya(n-1.0,x) - (n-1.0)*polya(n-2.0,x) )/n;

}

移位勒让德多项式[编辑]

移位勒让德多项式\tilde{P_n}(x)的正交区间定义在[0,1]上,即:

\int_{0}^{1} \tilde{P_m}(x) \tilde{P_n}(x)\,dx = {1 \over {2n + 1}} \delta_{mn}.

其显式表达式为:

\tilde{P_n}(x)=(-1)^n \sum_{k=0}^n {n \choose k} {n+k \choose k} (-x)^k.

相应的罗德里格公式为:

\tilde{P_n}(x) = ( n!)^{-1} {d^n \over dx^n } \left[ (x^2 -x)^n \right].\,

下表列出了头4阶移位勒让德多项式:

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

分数阶勒让德多项式[编辑]

分数阶勒让德多项式通过将分数阶微分(定义参见分数微积分理论)和通过Γ函数定义的非整数阶乘代入罗德里格公式中来定义。

参见[编辑]

外部链接[编辑]

参考文献[编辑]

  1. ^ 严镇军编,《数学物理方程》,第二版,中国科学技术大学出版社,合肥,2002,ISBN 7-312-00799-6/O·177,第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章