Legendre-polynoom

In de wiskunde is een legendre-polynoom een oplossing van de differentiaalvergelijking van Legendre. Soms echter bedoelt men de geassocieerde legendre-polynoom.

Differentiaalvergelijking

De vergelijking waarvan de polynomen een oplossing vormen luidt:

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

Beide zijn vernoemd naar de Franse wiskundige Adrien-Marie Legendre. De vergelijking komt regelmatig voor in de natuurkunde en de toegepaste wetenschappen, omdat de laplace-vergelijking in bolcoördinaten de gedaante van deze differentiaalvergelijking van Legendre heeft (althans voor rotatie-symmetrische gevallen, voor de θ-afhankelijkheid, waarbij θ de polaire hoek is)[1].

Een standaardmethode voor de oplossing van deze vergelijking is ontwikkeling van een machtreeks. De oplossing convergeert wanneer | x | < 1 {\displaystyle |x|<1} . Bovendien is de waarde eindig voor x = ± 1 {\displaystyle x=\pm 1} , vooropgesteld dat n {\displaystyle n} een niet-negatief geheel getal is. In dat geval vormen de oplossingen van de vergelijking een reeks van orthogonale polynomen, de legendre-polynomen

Iedere legendre-polynoom P n ( x ) {\displaystyle P_{n}(x)} is een veelterm van graad n {\displaystyle n} en kan uitgedrukt worden met de formule van Rodrigues:

P n ( x ) = 1 2 n n ! d n d x n [ ( x 2 1 ) n ] {\displaystyle P_{n}(x)={1 \over 2^{n}n!}{{\rm {d}}^{n} \over {\rm {d}}x^{n}}\left[(x^{2}-1)^{n}\right]}

Orthogonaliteit

Een belangrijke eigenschap van de legendre-polynomen is dat zij, zoals hierboven reeds vermeld, orthogonaal zijn met betrekking tot het L 2 {\displaystyle L^{2}} inproduct op het interval 1 x 1 {\displaystyle -1\leq x\leq 1} :

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

Daarin is δ m n {\displaystyle \delta _{mn}} de kronecker-delta die gelijk is aan 1 als m = n {\displaystyle m=n} en anders nul.

Een andere manier om de polynomen af te leiden is gebruik te maken van de gram-schmidtmethode op het inproduct van de polynomen { 1 , x , x 2 , } {\displaystyle \{1,x,x^{2},\ldots \}} De reden voor de orthogonaliteit van de polynomen is dat de legendre-vergelijking gezien kan worden als een sturm-liouvillevraagstuk

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

waar de eigenwaarde λ {\displaystyle \lambda } overeenkomt met n ( n + 1 ) {\displaystyle n(n+1)} .

Een andere belangrijke eigenschap van legendre-polynomen stelt dat:

1 1 f ( x ) P n ( x ) d x = 0 {\displaystyle \int \limits _{-1}^{1}f(x)P_{n}(x)\mathrm {d} x=0}

zodra f ( x ) {\displaystyle f(x)} een veelterm is van graad strikt kleiner dan n {\displaystyle n} .

Voorbeelden van legendre-polynomen

De eerste legendre-polynomen zijn:

n {\displaystyle n} P n ( x ) {\displaystyle P_{n}(x)} normalisatiefactor
0 1 {\displaystyle 1} 1 2 {\displaystyle {\sqrt {\tfrac {1}{2}}}}
1 x {\displaystyle x} 3 2 {\displaystyle {\sqrt {\tfrac {3}{2}}}}
2 1 2 ( 3 x 2 1 ) {\displaystyle {\begin{matrix}{\frac {1}{2}}\end{matrix}}(3x^{2}-1)} 5 2 {\displaystyle {\sqrt {\tfrac {5}{2}}}}
3 1 2 ( 5 x 3 3 x ) {\displaystyle {\begin{matrix}{\frac {1}{2}}\end{matrix}}(5x^{3}-3x)} 7 2 {\displaystyle {\sqrt {\tfrac {7}{2}}}}
4 1 8 ( 35 x 4 30 x 2 + 3 ) {\displaystyle {\begin{matrix}{\frac {1}{8}}\end{matrix}}(35x^{4}-30x^{2}+3)} 9 2 {\displaystyle {\sqrt {\tfrac {9}{2}}}}
5 1 8 ( 63 x 5 70 x 3 + 15 x ) {\displaystyle {\begin{matrix}{\frac {1}{8}}\end{matrix}}(63x^{5}-70x^{3}+15x)} 11 2 {\displaystyle {\sqrt {\tfrac {11}{2}}}}
6 1 16 ( 231 x 6 315 x 4 + 105 x 2 5 ) {\displaystyle {\begin{matrix}{\frac {1}{16}}\end{matrix}}(231x^{6}-315x^{4}+105x^{2}-5)} 13 2 {\displaystyle {\sqrt {\tfrac {13}{2}}}}

In de onderstaande figuur staan de grafieken van de eerste zes legendre-polynomen.

Toepassingen in de natuurkunde

Legendre-polynomen bewijzen hun nut bij de reeksontwikkeling van functies zoals

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

hierbij zijn r {\displaystyle r} en r {\displaystyle r'} respectievelijk de lengte van de vectoren x {\displaystyle \mathbf {x} } en x {\displaystyle \mathbf {x} ^{\prime }} en γ {\displaystyle \gamma } is de hoek tussen de beide vectoren. Deze ontwikkeling is geldig als r > r {\displaystyle r>r'} .

Deze uitdrukking wordt bijvoorbeeld gebruikt om de potentiaal van een puntlading te vinden zoals deze ondervonden wordt in het punt x {\displaystyle \mathbf {x} } wanneer de lading zich op punt x {\displaystyle \mathbf {x} '} bevindt. De ontwikkeling is vooral van nut wanneer men een integraal over een continue ladingsverdeling wil berekenen.

Legendre-polynomen komen voor als oplossingen van de laplace-vergelijking en de poissonvergelijking voor de potentiaal 2 Φ ( x ) {\displaystyle \nabla ^{2}\Phi (\mathbf {x} )} , indien bolcoördinaten worden gebruikt. Daarbij wordt gebruikgemaakt van de methode van scheiding van variabelen en wordt axiale symmetrie verondersteld, zodat er geen afhankelijkheid is van de azimuthale hoek. z ^ {\displaystyle {\widehat {\mathbf {z} }}} staat voor de symmetrieas en θ {\displaystyle \theta } is de hoek tussen de waarnemer en de as z ^ {\displaystyle {\widehat {\mathbf {z} }}} . De oplossing is dan:

Φ ( r , θ ) = k = 0 [ A k r k + B k r ( k + 1 ) ] P k ( cos θ ) {\displaystyle \Phi (r,\theta )=\sum _{k=0}^{\infty }\left[A_{k}r^{k}+B_{k}r^{-(k+1)}\right]P_{k}(\cos \theta )}

A k {\displaystyle A_{k}} and B k {\displaystyle B_{k}} moeten voor de randvoorwaarden van het betreffende probleem bepaald worden [2].

Legendre-polynomen in multipoolontwikkelingen

Figuur 2

De legendre-polynoom is ook van nut bij een reeksontwikkeling van de vorm

1 1 + η 2 2 η x = k = 0 η k P k ( x ) {\displaystyle {\frac {1}{\sqrt {1+\eta ^{2}-2\eta x}}}=\sum _{k=0}^{\infty }\eta ^{k}P_{k}(x)}

In feite is dit dezelfde ontwikkeling als hierboven in een wat andere vorm die voortvloeit uit de reeksontwikkeling van multipolen. De linkerzijde van de vergelijking genereert legendre-polynomen

Bijvoorbeeld, de elektrische potentiaal Φ ( r , θ ) {\displaystyle \Phi (r,\theta )} in bolcoördinaten die voortvloeit uit een puntlading op de z {\displaystyle z} -as in het punt z = a {\displaystyle z=a} (Fig. 2) is te schrijven als

Φ ( r , θ ) 1 R = 1 r 2 + a 2 2 a r cos θ {\displaystyle \Phi (r,\theta )\propto {\frac {1}{R}}={\frac {1}{\sqrt {r^{2}+a^{2}-2ar\cos \theta }}}}

Als de straal r {\displaystyle r} vanaf het punt van waarneming P {\displaystyle \mathbf {P} } veel groter is dan a {\displaystyle a} , is het mogelijk de potentiaal te ontwikkelen tot legendre-polynomen:

Φ ( r , θ ) 1 r k = 0 ( a r ) k P k ( cos θ ) {\displaystyle \Phi (r,\theta )\propto {\frac {1}{r}}\sum _{k=0}^{\infty }\left({\frac {a}{r}}\right)^{k}P_{k}(\cos \theta )}

hierbij is η = a / r < 1 {\displaystyle \eta =a/r<1} x = cos θ {\displaystyle x=\cos \theta } genomen.

In het omgekeerde geval r a {\displaystyle r\ll a} kunnen we ook ontwikkelen tot een legendre-polynoom, wanneer we de rol van r {\displaystyle r} en a {\displaystyle a} omkeren:

Overige eigenschappen

Legendre-polynomen zijn om en om symmetrisch en antisymmetrisch:

P k ( x ) = ( 1 ) k P k ( x ) {\displaystyle P_{k}(-x)=(-1)^{k}P_{k}(x)}

Aangezien de differentiaalvergelijking en de orthogonaliteit onafhankelijk zijn van de schaal zijn is de veelterm van nature gestandaardiseerd. Dit betekent dat P k ( 1 ) = 1 {\displaystyle P_{k}(1)=1} .

Soms noemt men dit genormaliseerd, maar dit is wat verwarrend omdat de norm niet een is, immers:

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

Een orthonormale versie van de polynoom kan verkregen worden door toevoeging van een normalisatiefactor n + 1 2 {\displaystyle {\sqrt {n+{\tfrac {1}{2}}}}}

P n , n o r m ( x ) = n + 1 2 2 n n ! d n d x n ( x 2 1 ) n {\displaystyle P_{n,norm}(x)={{\sqrt {n+{\tfrac {1}{2}}}} \over 2^{n}n!}{{\rm {d}}^{n} \over {\rm {d}}x^{n}}(x^{2}-1)^{n}}

De afgeleide aan de eindpunten wordt gegeven door:

P k ( 1 ) = k ( k + 1 ) 2 {\displaystyle P_{k}'(1)={\frac {k(k+1)}{2}}}

Legendre-polynomen kunnen opgebouwd worden door gebruik te maken van de relaties

( n + 1 ) P n + 1 = ( 2 n + 1 ) x P n n P n 1 {\displaystyle (n+1)P_{n+1}=(2n+1)xP_{n}-nP_{n-1}}

en

x 2 1 n d d x P n = x P n P n 1 {\displaystyle {x^{2}-1 \over n}{{\rm {d}} \over {\rm {d}}x}P_{n}=xP_{n}-P_{n-1}}

Een nuttige uitdrukking bij het integreren van de veelterm is:

( 2 n + 1 ) P n = d d x [ P n + 1 P n 1 ] {\displaystyle (2n+1)P_{n}={{\rm {d}} \over {\rm {d}}x}\left[P_{n+1}-P_{n-1}\right]}
Voetnoten
  1. Lorrain and Corson, Electromagnetic fields and waves, 2nd ed., Freeman, San Francisco 1970, p. 165
  2. Jackson, J.D. Classical Electrodynamics, 3rd edition, Wiley & Sons, 1999. page 103