勒让德多项式（Legendre Polynomials）推导

推导Legendre Polynomials（勒让德多项式）

问题：

证明Legendre Polylnomails（勒让德多项式）是Legendre Differential Equation（勒让德微分方程）的解
$$勒让德微分方程:(1-x^2)\frac{d^{2}y}{d{x}^2}-2x\frac{dy}{dx}+k(k+1)y=0$$
有解
$$勒让德多项式：P_k(x)=\sum _{m=0}^{\frac{k}{2}|\frac{k-1}{2}}\frac{(-1{)}^m(2k-2m)!}{m!(2^k)(k-m)!(k-2m)!}$$

推导出递推关系（Recurrence Relation）

证明：首先我们变换一下微分方程。检查$\frac{2x}{(1-x^2)}$，和$\frac{k(k+1)}{(1-x^2)}$ 满足是用幂级数（Power Series）的方法。

$$\frac{d^{2}y}{d{x}^2}-\frac{2x}{(1-x^2)}\frac{dy}{dx}+\frac{k(k+1)}{(1-x^2)}y=0$$
我们可以套用幂级数（Power Series）的方法求解微分方程。
$$y(x)=\sum _{n=0}^{\infty }{a}_n{x}^n$$一阶微分$$y^{{}^{\prime }}(x)=\sum _{n=0}^{\infty }n{a}_n{x}^{n-1}$$二阶微分$$y^{{}^{\prime \prime }}(x)=\sum _{n=0}^{\infty }n(n-1)a_n{x}^{n-2}$$
把$y^{}(x)$，$y^{{}^{\prime }}(x)$和$y^{{}^{\prime \prime }}(x)$带入微分方程中，化简可得。
$$\begin{array}{rl}(1-x^2)\sum _{n=0}^{\infty }n(n-1)a_n{x}^{n-2}-2x\sum _{n=0}^{\infty }n{a}_n{x}^{n-1}+k(k+1)\sum _{n=0}^{\infty }{a}_n{x}^n& =0\\ \sum _{n=0}^{\infty }n(n-1)a_n{x}^{n-2}-x^2\sum _{n=0}^{\infty }n(n-1)a_n{x}^{n-2}-2x\sum _{n=0}^{\infty }n{a}_n{x}^{n-1}+k(k+1)\sum _{n=0}^{\infty }{a}_n{x}^n& =0\\ \sum _{n=0}^{\infty }n(n-1)a_n{x}^{n-2}-\sum _{n=0}^{\infty }n(n-1)a_n{x}^n-2\sum _{n=0}^{\infty }n{a}_n{x}^n+k(k+1)\sum _{n=0}^{\infty }{a}_n{x}^n& =0\end{array}$$
针对方程中的第一项${\sum }_{n=0}^{\infty }n(n-1)a_n{x}^{n-2}$中的角标变量$n$用虚拟变量$m=n-2$来替换。当$n=0$的时候，$m=-2$并且$n=m+2$。所以第一项可以整理得：
$$\sum _{n=0}^{\infty }n(n-1)a_n{x}^{n-2}=\sum _{m=-2}^{\infty }(m+2)(m+1)a_{m+2}{x}^m$$
我们展开第一项，当$m=-1$和$m=-2$的时候均等于$0$。所以角标变量可以从$m=0$开始。因为$m$是虚拟变量所以我们用$n$来替换回$m$。原微分方程整理可得：

$$\begin{array}{rl}\sum _{n=0}^{\infty }(n+2)(n+1)a_{n+2}{x}^n-\sum _{n=0}^{\infty }n(n-1)a_n{x}^n-2\sum _{n=0}^{\infty }n{a}_n{x}^n+k(k+1)\sum _{n=0}^{\infty }{a}_n{x}^n& =0\\ \sum _{n=0}^{\infty }\left[(n+2)(n+1)a_{n+2}-n(n-1)a_n-2n{a}_n+k(k+1)a_n\right]x^n& =0\end{array}$$

所以系数等于$0$。
$$(n+2)(n+1)a_{n+2}-n(n-1)a_n-2n{a}_n+k(k+1)a_n=0$$
于是我们找到了一个递推关系（recurrence relation）的等式，化简可得。
$$\begin{array}{rl}{a}_{n+2}& =a_n\frac{n(n-1)+2n-k(k+1)}{(n+2)(n+1)}\\ & =a_n\frac{(n-k)(n+k+1)}{(n+2)(n+1)}\end{array}$$

推导勒让德多项式

假设$a_0$ 是勒让德多项式有偶数自由度的$k$的首项系数，$a_1$ 是勒让德多项式有基数自由度的$k$的首项系数。
假设$a_k$ 是勒让德多项式最后一项系数。为了数学上的方便，我们让$a_k$赋予一个特定常数值：
$$a_k=\frac{(2k)!}{2^k(k!{)}^2}$$
假设为了数学上的方便，我们让$P_k(x=1)=1$。
由出递推关系（Recurrence Relation）得：
$$\begin{array}{rl}{a}_n& =\frac{(n+2)(n+1)}{(n-k)(n+k+1)}{a}_{n+2}\\ a_{k-2}& =\frac{k(k-1)}{(-2)(2k-1)}{a}_k\\ & =\frac{k(k-1)}{(-2)(2k-1)}\frac{(2k)!}{2^k(k!{)}^2}\\ & =\frac{-1(k)(k-1)}{2(2k-1)}\frac{2k(2k-1)(2k-2)!}{2^k{k}^2(k-1{)}^2((k-2)!{)}^2}\\ & =\frac{-1(2k-2)!}{2^k(k-1)!(k-2)!}\\ & ⋮\\ a_{k-2m}& =\frac{(-1{)}^m(2k-2m)!}{m!(2^k)(k-m)!(k-2m)!}\end{array}$$其中$m\in [0,k/2]$ 勒让德多项式有偶数自由度的$k$；$m\in [0,\frac{k-1}{2}]$ 勒让德多项式有基数自由度的$k$。
由此证明的：
$$P_k(x)=\sum _{m=0}^{\frac{k}{2}|\frac{k-1}{2}}\frac{(-1{)}^m(2k-2m)!}{m!(2^k)(k-m)!(k-2m)!}{x}^{k-2m}$$

推导 Rodrigues Formular （罗德里格公式）

问题

求证恒等式左边等于右边
$$\frac{1}{2^{k}k!}\frac{{\mathrm{d}}^k}{\mathrm{d}{x}^k}\left[(x^2-1{)}^k\right]=\sum _{m=0}^{\frac{k}{2}|\frac{k-1}{2}}\frac{(-1{)}^m(2k-2m)!x^{k-2m}}{m!(2^k)(k-m)!(k-2m)!}$$

应用二项式定理

解：
我们在$(x^2-1{)}^k$上应用二项式定理。由二项式定理可知：
$$\begin{array}{rl}(x+y{)}^k& =\sum _{i=0}^n\frac{n!}{i!(n-i)!}{x}^{n-i}{y}^i\\ (x^2-1{)}^k& =\sum _{i=0}^k\frac{k!}{i!(k-i)!}(x^2{)}^{k-i}(-1{)}^i\end{array}$$

公式证明

我们把有二项式定理分解的$(x^2-1{)}^k$带入罗德里格公式部分得：
$$\frac{1}{2^{k}k!}\frac{{\mathrm{d}}^k}{\mathrm{d}{x}^k}\left[\sum _{i=0}^k\frac{k!}{i!(k-i)!}(x^2{)}^{k-i}(-1{)}^i\right]$$
因为微分是线性运算符所以可知$\frac{\mathrm{d}}{\mathrm{d}x}\sum =\sum \frac{\mathrm{d}}{\mathrm{d}x}$，由此化简可得：
$$\begin{array}{rl}& =\frac{1}{2^{k}k!}\sum _{i=0}^k\frac{k!}{i!(k-i)!}(-1{)}^i\frac{{\mathrm{d}}^k}{\mathrm{d}{x}^k}\left[(x^2{)}^{k-i}\right]\\ & =\frac{1}{2^{k}k!}\sum _{i=0}^k\frac{k!}{i!(k-i)!}(-1{)}^i\frac{(2k-2i)!}{(k-2i)!}{x}^{k-2i}\\ & =\sum _{i=0}^{\frac{k}{2}|\frac{k-1}{2}}\frac{(-1{)}^i(2k-2i)!}{2^{k}i!(k-i)!(k-2i)!}{x}^{k-2i}\end{array}$$
其中
$$\begin{array}{rl}\frac{{\mathrm{d}}^n}{\mathrm{d}{x}^n}{x}^{r}forn=1& :r{x}^{r-1}=\frac{r!}{(r-1)!}{x}^{r-1}\\ n=2& :r(r-1)x^{r-2}=\frac{r!}{(r-2)!}{x}^{r-2}\\ & ⋮\\ n=n& :r(r-1)\cdots (r-n)x^{r-n}=\frac{r!}{(r-n)!}{x}^{r-n}\end{array}$$
由此可得：
$$\begin{array}{rl}\frac{{\mathrm{d}}^k}{\mathrm{d}{x}^k}\left[(x^2{)}^{k-i}\right]& =\frac{{\mathrm{d}}^k}{\mathrm{d}{x}^k}\left(x^{2k-2i}\right)\\ & =\frac{(2k-2i)!}{(2k-2i-k)!}{x}^{2k-2i-k}\\ & =\frac{(2k-2i)!}{(k-2i)!}{x}^{k-2i}\end{array}$$
由此证明结论：
$$P_k(x)=\frac{1}{2^{k}k!}\frac{{\mathrm{d}}^k}{\mathrm{d}{x}^k}\left[(x^2-1{)}^k\right]$$