• 罗德里格斯公式
• 正交性
• 归一化
• 应用

罗德里格斯公式

勒让德多项式的一个重要物理背景就是球坐标下的静电场满足的拉普拉斯方程， ∇ 2 V = 0 \nabla^2V=0 ∇2V=0，即
r − 2 ∂ r 2 ( r 2 V ) + r − 2 ( s i n ( θ ) ) − 1 ∂ θ 2 V + r − 2 ( s i n ( θ ) ) − 2 ∂ ϕ 2 V = 0 r^{-2}\partial_r^2(r^2V)+r^{-2}(sin(\theta))^{-1}\partial_\theta^2V+r^{-2}(sin(\theta))^{-2}\partial_\phi^2V=0 r−2∂r2​(r2V)+r−2(sin(θ))−1∂θ2​V+r−2(sin(θ))−2∂ϕ2​V=0
分离变量： V = R ( r ) Θ ( θ ) Φ ( ϕ ) V=R(r)\Theta(\theta)\Phi(\phi) V=R(r)Θ(θ)Φ(ϕ)，得 Θ \Theta Θ满足得常微分方程为：
1 s i n θ d d θ ( s i n θ d Θ d θ ) + l ( l + 1 ) Θ − m 2 s i n 2 θ Θ = 0 \frac{1}{sin\theta}\frac{d}{d\theta}\left(sin\theta\frac{d\Theta}{d\theta}\right)+l(l+1)\Theta-\frac{m^2}{sin^2\theta}\Theta=0 sinθ1​dθd​(sinθdθdΘ​)+l(l+1)Θ−sin2θm2​Θ=0.
进行变量代换： x = c o s ( θ ) , d x = − s i n θ d θ x=cos(\theta), dx=-sin\theta d\theta x=cos(θ),dx=−sinθdθ，并假定 m = 0 m=0 m=0，即得到勒让德方程：
d d x ( ( 1 − u 2 ) d Θ d x ) + l ( l + 1 ) Θ = 0 , x ∈ [ − 1 , 1 ] \frac{d}{dx}((1-u^2)\frac{d\Theta}{dx})+l(l+1)\Theta=0,x\in[-1,1] dxd​((1−u2)dxdΘ​)+l(l+1)Θ=0,x∈[−1,1]
该方程可通过级数法求解，所得级数与由罗德里格斯公式：
P l ( x ) = 1 2 l l ! u l ( l ) ( x ) , u l ( x ) = ( x 2 − 1 ) l P_l(x)=\frac{1}{2^ll!}u_l^{(l)}(x),u_l(x)=(x^2-1)^l Pl​(x)=2ll!1​ul(l)​(x),ul​(x)=(x2−1)l
展开后相同（参考数学物理方法教材）。
由于 u l ( x ) u_l(x) ul​(x)为2n次多项式，因此求 l l l次导后所得 P l ( x ) P_l(x) Pl​(x)为 l l l次多项式，且前几阶勒让德多项式为：
P 0 ( x ) = 1 , P 1 ( x ) = x , P 2 ( x ) = 1 2 ( 3 x 2 − 1 ) P_0(x)=1,P_1(x)=x,P_2(x)=\frac{1}{2}(3x^2-1) P0​(x)=1,P1​(x)=x,P2​(x)=21​(3x2−1)

正交性

正交性即证明 P m , P n = ∫ − 1 1 P m ( x ) P n ( x ) d x = 0 , m ≠ n \langle P_m,P_n\rangle=\int_{-1}^1P_m(x)P_n(x)dx=0,m\neq n Pm​,Pn​=∫−11​Pm​(x)Pn​(x)dx=0,m​=n.
设 m > n ≥ 0 m>n\geq0 m>n≥0，则 P n ( x ) P_n(x) Pn​(x)为n次多项式，因此只需证明对任意 k ≤ n k\leq n k≤n，均有 ∫ − 1 1 P m ( x ) x k d x = 0 \int_{-1}^1P_m(x)x^kdx=0 ∫−11​Pm​(x)xkdx=0.

1. 对于 k = 0 k=0 k=0，有：
   ∫ − 1 1 P m ( x ) d x = 1 2 m m ! u m ( m − 1 ) ∣ − 1 1 \int_{-1}^1P_m(x)dx=\frac{1}{2^mm!}u_m^{(m-1)}\big|^{1}_{-1} ∫−11​Pm​(x)dx=2mm!1​um(m−1)​∣∣​−11​
   而对任意 m ′ < m m'<m m′<m，由莱布尼兹公式：
   u m ( m ′ ) = ∑ m ′ ′ = 0 m ′ [ ( x + 1 ) m ] ( m ′ ′ ) [ ( x − 1 ) m ] ( m ′ − m ′ ′ ) u_m^{(m')}=\sum_{m''=0}^{m'}[(x+1)^m]^{(m'')}[(x-1)^m]^{(m'-m'')} um(m′)​=m′′=0∑m′​[(x+1)m](m′′)[(x−1)m](m′−m′′)
   由于 m ′ ′ , m ′ − m ′ ′ < m m'',m'-m''<m m′′,m′−m′′<m
   u m ( m ′ ) ( 1 ) = u m ( m ′ ) ( − 1 ) = 0 u_m^{(m')}(1)=u_m^{(m')}(-1)=0 um(m′)​(1)=um(m′)​(−1)=0
   （注意当 m ′ = 0 m'=0 m′=0时上式同样成立）
   因此， ∫ − 1 1 P m ( x ) d x = 0 \int_{-1}^1P_m(x)dx=0 ∫−11​Pm​(x)dx=0
2. 对于 k ≥ 1 k\geq1 k≥1，使用分部积分：
   ∫ − 1 1 P m ( x ) x k d x = 1 2 m m ! ( x k u m ( m − 1 ) ∣ − 1 1 − k ∫ − 1 1 u m ( m − 1 ) ( x ) x k − 1 d x ) \int_{-1}^1P_m(x)x^kdx=\frac{1}{2^mm!}\left(x^ku_m^{(m-1)}\big|^1_{-1}-k\int_{-1}^1u_m^{(m-1)}(x)x^{k-1}dx\right) ∫−11​Pm​(x)xkdx=2mm!1​(xkum(m−1)​∣∣​−11​−k∫−11​um(m−1)​(x)xk−1dx)
   其中第一项为0，重复使用分部积分，最终得到：
   ∫ − 1 1 P m ( x ) x k d x = ( − 1 ) k k ! 2 m m ! ∫ − 1 1 u m ( m − k ) ( x ) d x = u m ( m − k − 1 ) ∣ − 1 1 = 0 \int_{-1}^1P_m(x)x^kdx=\frac{(-1)^kk!}{2^mm!}\int_{-1}^1u_m^{(m-k)}(x)dx\\ \quad \\ =u_m^{(m-k-1)}\big|^1_{-1}=0 ∫−11​Pm​(x)xkdx=2mm!(−1)kk!​∫−11​um(m−k)​(x)dx=um(m−k−1)​∣∣​−11​=0

正交性得证；

归一化

勒让德多项式的模为：
∥ P n ∥ L 2 = ∫ − 1 1 P n ( x ) P n ( x ) d x \|P_n\|_{L^2}=\sqrt{\int_{-1}^1P_n(x)P_n(x)dx} ∥Pn​∥L2​=∫−11​Pn​(x)Pn​(x)dx ​
当 n = 0 n=0 n=0，有 ∥ P 0 ∥ L 2 = 2 \|P_0\|_{L^2}=\sqrt{2} ∥P0​∥L2​=2 ​
当 n ≥ 1 n\geq1 n≥1，有：
∫ − 1 1 P n ( x ) P n ( x ) d x = 1 ( 2 n n ! ) 2 ∫ − 1 1 u n ( n ) ( x ) u n ( n ) ( x ) d x = 1 ( 2 n n ! ) 2 ( u n ( n ) u n ( n − 1 ) ∣ − 1 1 − ∫ − 1 1 u n ( n − 1 ) ( x ) u n ( n + 1 ) ( x ) d x ) = ( − 1 ) n ( 2 n n ! ) 2 ∫ − 1 1 u n ( x ) u n ( 2 n ) ( x ) d x \int_{-1}^1P_n(x)P_n(x)dx=\frac{1}{(2^nn!)^2}\int_{-1}^1u_n^{(n)}(x)u_n^{(n)}(x)dx\\ \quad \\ =\frac{1}{(2^nn!)^2}\left(u_n^{(n)}u_n^{(n-1)}\big|^1_{-1}-\int_{-1}^1u_n^{(n-1)}(x)u_n^{(n+1)}(x)dx\right) \\\quad\\=\frac{(-1)^n}{(2^nn!)^2}\int_{-1}^1u_n(x)u_n^{(2n)}(x)dx ∫−11​Pn​(x)Pn​(x)dx=(2nn!)21​∫−11​un(n)​(x)un(n)​(x)dx=(2nn!)21​(un(n)​un(n−1)​∣∣​−11​−∫−11​un(n−1)​(x)un(n+1)​(x)dx)=(2nn!)2(−1)n​∫−11​un​(x)un(2n)​(x)dx
由于 u n u_n un​最高次项为 x 2 n x^{2n} x2n， u n ( 2 n ) ( x ) = ( 2 n ) ! u_n^{(2n)}(x)=(2n)! un(2n)​(x)=(2n)!。令 I n = ∫ − 1 1 u n ( x ) d x I_n=\int_{-1}^1u_n(x)dx In​=∫−11​un​(x)dx，使用分部积分：
I n = x u n ∣ − 1 1 − 2 n ∫ − 1 1 x 2 ( x 2 − 1 ) n − 1 d x = − 2 n I n − 2 n I n − 1 I_n=xu_n\big|_{-1}^1-2n\int_{-1}^1x^2(x^2-1)^{n-1}dx=-2nI_n-2nI_{n-1} In​=xun​∣∣​−11​−2n∫−11​x2(x2−1)n−1dx=−2nIn​−2nIn−1​
移项，并重复这一步骤，可得：
I n = ( − 1 ) n ∏ k = 1 n ( 2 k ) ∏ k = 1 n ( 2 k + 1 ) I 0 I_n=\frac{(-1)^n\prod_{k=1}^n(2k)}{\prod_{k=1}^n(2k+1)}I_0 In​=∏k=1n​(2k+1)(−1)n∏k=1n​(2k)​I0​
而 I 0 = 2 I_0=2 I0​=2,于是，勒让德多项式模为：
∥ P 0 ∥ L 2 = ( − 1 ) n ( 2 n ) ! ( 2 n n ! ) 2 2 ( − 1 ) n ∏ k = 1 n ( 2 k ) ∏ k = 1 n ( 2 k + 1 ) = 2 2 n + 1 \|P_0\|_{L^2}=\sqrt{\frac{(-1)^n(2n)!}{(2^nn!)^2}\frac{2(-1)^n\prod_{k=1}^n(2k)}{\prod_{k=1}^n(2k+1)}}=\sqrt{\frac{2}{2n+1}} ∥P0​∥L2​=(2nn!)2(−1)n(2n)!​∏k=1n​(2k+1)2(−1)n∏k=1n​(2k)​ ​=2n+12​ ​
归一化的勒让德多项式即为：
Q n = 2 n + 1 2 P n Q_n=\sqrt{\frac{2n+1}{2}}P_n Qn​=22n+1​ ​Pn​

应用

由维尔斯特拉斯定理，任意闭区间上连续函数 f f f可用某一多项式逼近，而这一多项式就可表示为勒让德多项式之和：
p ( x ) = ∑ n = 0 ∞ Q n , f Q n p(x)=\sum_{n=0}^{\infty}\langle Q_n,f\rangle Q_n p(x)=n=0∑∞​Qn​,fQn​

与泰勒展开不同，这里未假设连续函数 f f f可导。
作为例子，设
f ( x ) = t r i ( 2 x − 1 ) , x ∈ [ − 1 , 1 ] f(x)=tri(2x-1),x\in[-1,1] f(x)=tri(2x−1),x∈[−1,1]
将该函数展开到不同阶，如图所示

代码如下：

x=-1:.01:1;
len=length(x);
f_0=@f0;
coef(f_0,0)
subplot(3,1,1)
p_n=zeros(1,len);
for n=0:3
   p_n=p_n+coef(f_0,n)*Q(x,n,0);
end
plot(x,f0(x),x,p_n),axis([-1,1,-0.2,1.1])
ylabel('p_3(x)')

subplot(3,1,2)
for n=4:8
   p_n=p_n+coef(f_0,n)*Q(x,n,0);
end
plot(x,f0(x),x,p_n),axis([-1,1,-0.2,1.1])
ylabel('p_{8}(x)')

subplot(3,1,3)
for n=9:40
   p_n=p_n+coef(f_0,n)*Q(x,n,0);
end
plot(x,f0(x),x,p_n),axis([-1,1,-0.2,1.1])
ylabel('p_{40}(x)')
xlabel('x')

function f=f0(x)
n=length(x);
temp=zeros(1,n);
for i=1:n
if x(i)>0&&x(i)<1
   temp(i)=1-abs(2*(x(i)-1/2));
end
end
f=temp;
end

function f=coef(f0,n)
p_f0=@(x)(f0(x).*Q(x,n,0));
f=integral(p_f0,-1,1);
end

function f=Q(x,n,m)
temp=legendre(n,x);
f=sqrt((2*n+1)/2)*temp(m+1,:);
end

在求解电磁波散射问题时，就可以在球坐标下将入射场展开为勒让德多项式，并将散射场展开而其系数待定，再通过边界条件求得散射场展开系数。

参考：汪德新，数学物理方法