From 6ce0bb9a8f9fd7d169cbb414a9537d68c5290aae Mon Sep 17 00:00:00 2001 From: Prefetch Date: Fri, 14 Oct 2022 23:25:28 +0200 Subject: Initial commit after migration from Hugo --- source/know/concept/legendre-polynomials/index.md | 119 ++++++++++++++++++++++ 1 file changed, 119 insertions(+) create mode 100644 source/know/concept/legendre-polynomials/index.md (limited to 'source/know/concept/legendre-polynomials') diff --git a/source/know/concept/legendre-polynomials/index.md b/source/know/concept/legendre-polynomials/index.md new file mode 100644 index 0000000..338b23f --- /dev/null +++ b/source/know/concept/legendre-polynomials/index.md @@ -0,0 +1,119 @@ +--- +title: "Legendre polynomials" +date: 2021-09-08 +categories: +- Mathematics +layout: "concept" +--- + +The **Legendre polynomials** are a set of functions that sometimes arise in physics. +They are the eigenfunctions $u(x)$ of **Legendre's differential equation**, +which is a ([Sturm-Liouville](/know/concept/sturm-liouville-theory/)) +eigenvalue problem for $\ell (\ell + 1)$, +where $\ell$ turns out to be a non-negative integer: + +$$\begin{aligned} + \boxed{ + (1 - x^2) u'' - 2 x u' + \ell (\ell + 1) u = 0 + } +\end{aligned}$$ + +The $\ell$th-degree Legendre polynomial $P_\ell(x)$ +is given in the form of a *Rodrigues' formula* by: + +$$\begin{aligned} + P_\ell(x) + &= \frac{1}{2^\ell \ell!} \dvn{\ell}{}{x}(x^2 - 1)^\ell +\end{aligned}$$ + +The first handful of Legendre polynomials $P_\ell(x)$ are therefore as follows: + +$$\begin{gathered} + P_0(x) = 1 + \qquad \quad + P_1(x) = x + \qquad \quad + P_2(x) = \frac{1}{2} (3 x^2 - 1) + \\ + P_3(x) = \frac{1}{2} (5 x^3 - 3 x) + \qquad \quad + P_4(x) = \frac{1}{8} (35 x^4 - 30 x^2 + 3) +\end{gathered}$$ + +And then more $P_\ell$ can be computed quickly +using **Bonnet's recursion formula**: + +$$\begin{aligned} + \boxed{ + (\ell + 1) P_{\ell + 1}(x) = (2 \ell + 1) x P_\ell(x) - \ell P_{\ell - 1}(x) + } +\end{aligned}$$ + +The derivative of a given $P_\ell$ can be calculated recursively +using the following relation: + +$$\begin{aligned} + \boxed{ + \dv{}{x}P_{\ell + 1} + = (\ell + 1) P_\ell(x) + x \dv{}{x}P_\ell(x) + } +\end{aligned}$$ + +Noteworthy is that the Legendre polynomials +are mutually orthogonal for $x \in [-1, 1]$: + +$$\begin{aligned} + \boxed{ + \Inprod{P_m}{P_n} + = \int_{-1}^{1} P_m(x) \: P_n(x) \dd{x} + = \frac{2}{2 n + 1} \delta_{nm} + } +\end{aligned}$$ + +As was to be expected from Sturm-Liouville theory. +Likewise, they form a complete basis in the +[Hilbert space](/know/concept/hilbert-space/) +of piecewise continuous functions $f(x)$ on $x \in [-1, 1]$, +meaning: + +$$\begin{aligned} + \boxed{ + f(x) + = \sum_{\ell = 0}^\infty a_\ell P_\ell(x) + = \sum_{\ell = 0}^\infty \frac{\Inprod{P_\ell}{f}}{\Inprod{P_\ell}{P_\ell}} P_\ell(x) + } +\end{aligned}$$ + +Each Legendre polynomial $P_\ell$ comes with +a set of **associated Legendre polynomials** $P_\ell^m(x)$ +of order $m$ and degree $\ell$. +These are the non-singular solutions of the **general Legendre equation**, +where $m$ and $\ell$ are integers satisfying $-\ell \le m \le \ell$: + +$$\begin{aligned} + \boxed{ + (1 - x^2) u'' - 2 x u' + \Big( \ell (\ell + 1) - \frac{m^2}{1 - x^2} \Big) u = 0 + } +\end{aligned}$$ + +The $\ell$th-degree $m$th-order associated Legendre polynomial $P_\ell^m$ +is as follows for $m \ge 0$: + +$$\begin{aligned} + P_\ell^m(x) + = (-1)^m (1 - x^2)^{m/2} \dvn{m}{}{x}P_\ell(x) +\end{aligned}$$ + +Here, the $(-1)^m$ in front is called the **Condon-Shortley phase**, +and is omitted by some authors. +For negative orders $m$, +an additional constant factor is necessary: + +$$\begin{aligned} + P_\ell^{-m}(x) = (-1)^m \frac{(\ell - m)!}{(\ell + m)!} P_\ell^m(x) +\end{aligned}$$ + +Beware, the name is misleading: +if $m$ is odd, then $P_\ell^m$ is actually not a polynomial. +Moreover, not all $P_\ell^m$ are mutually orthogonal +(but some are). -- cgit v1.2.3