Lagrangesche Inversionsformel

Die Lagrangesche Inversionsformel in der Mathematik entwickelt zu einer gegebenen analytischen Funktion die Potenzreihe der Umkehrfunktion.

Aussage

Gegeben sei eine Gleichung

$z=f(w)$

mit einer am Punkt $a$ analytischen Funktion $f$ und Vorlage:Nowrap Dann ist es möglich, $f$ zu invertieren, also die Gleichung nach $w$ in Form einer formalen Potenzreihe $w=g(z)$ aufzulösen:

$g(z)=a+\sum _{n=1}^{\mathrm{\infty }}{g}_n\frac{(z-f(a){)}^n}{n!}$

mit

$g_n=\underset{w\to a}{\lim }\left[\frac{d^{n-1}}{d{w}^{n-1}}{\left(\frac{w-a}{f(w)-f(a)}\right)}^n\right].$

Die Potenzreihe $g(z)$ hat einen von 0 verschiedenen Konvergenzradius, d. h., sie ist eine analytische Funktion in einer Umgebung des Punktes Vorlage:Nowrap Die Formel invertiert $f$ als formale Potenzreihe in Vorlage:Nowrap Sie kann zu einer Formel für $H(g(z))$ mit einer beliebigen formalen Potenzreihe $H$ erweitert und in vielen Fällen mit $f^{\prime }(a)=0$ (dann eine „mehrwertige“ Funktion) verallgemeinert werden.

Der Satz wurde von Lagrange bewiesen und von Hans Heinrich Bürmann verallgemeinert, beides im späten 18. Jahrhundert. Es gibt Weiterentwicklungen in Richtung komplexe Analysis und Kurvenintegrale.

Taylorreihe

Die obige Formel gibt für eine formale Potenzreihe $f$ nicht direkt die Koeffizienten der formalen Umkehrfunktion $g$ ausgedrückt in den Koeffizienten von Vorlage:Nowrap Kann man die Funktionen $f$ und $g$ als formale Potenzreihe

$f(w)=\sum _{k=0}^{\mathrm{\infty }}{b}_k\frac{w^k}{k!}\mathrm{u}\mathrm{n}\mathrm{d}g(z)=\sum _{k=0}^{\mathrm{\infty }}{c}_k\frac{z^k}{k!}$

mit $b_0=0$ und $b_1\ne 0$ ausdrücken, dann können die Koeffizienten der Inversen $g$ mithilfe von Bell-Polynomen angegeben werden:

$c_n=\frac{1}{b_1^n}\sum _{k=1}^{n-1}(-1{)}^k{n}^{(k)}{B}_{n-1,k}(a_1,a_2,\dots ,a_{n-k})(n\ge 2)$ ,

mit $a_k=\frac{b_{k+1}}{(k+1)b_1},$    $c_1=\frac{1}{b_1}$   und   $n^{(k)}=n(n+1)\cdots (n+k-1)$  als steigender Faktorielle.

Explizite Formel

Die folgende explizite Formel gilt nicht nur für analytische Funktionen (über $ℝ$ oder $ℂ$), sondern für alle formalen Potenzreihen über einem Ring $R$ mit Eins. Ist nämlich

$A(X)=\sum _{k=1}^{\mathrm{\infty }}{a}_k{X}^k\in R[[X]]$

eine formale Potenzreihe, dann hat $A$ genau dann eine (formale) Umkehrfunktion (ein formales kompositionelles Inverses)

$B(X)=\sum _{n=1}^{\mathrm{\infty }}{b}_n{X}^n\in R[[X]]$ ,

wenn der Koeffizient $a_1=A^{\prime }(0)$ invertierbar (eine Einheit) in $R$ ist.

Der einfacheren Rechnung halber substituieren wir $X$ durch $a_1^{-1}Y$ und schreiben

$C(Y):=A(a_1^{-1}Y)=:Y+\sum _{k=2}^{\mathrm{\infty }}{c}_k{Y}^k$

mit $c_k:=a_1^{-k}{a}_k$ für Vorlage:Nowrap

Die zugehörige formale Umkehrfunktion sei

$D(Y):=Y+\sum _{n=2}^{\mathrm{\infty }}{d}_n{Y}^n$,

so dass $C(D(Y))=D(Y)+\sum _{k=2}^{\mathrm{\infty }}{c}_k(D(Y){)}^k=Y$ ist. Die Koeffizienten von $D$ lassen sich durch Koeffizientenvergleich in der Gleichung

$D(Y)=Y-\sum _{k=2}^{\mathrm{\infty }}{c}_k(D(Y){)}^k$

für $n\ge 2$ sofort zu

$\begin{array}{ccc}{d}_n& =-\left[Y^n\right]& \sum _{k=2}^{\mathrm{\infty }}{c}_k(D(Y){)}^k\end{array}$

$\begin{array}{ccccccc}& =-\left[Y^n\right]& (& c_2(1Y^1& +d_2{Y}^2& +\cdots +& d_{n-1}{Y}^{n-1}{)}^2\\ & & +& c_3& (1Y^1& +\cdots +& d_{n-2}{Y}^{n-2}{)}^3\\ & & ⋮\\ & & +& c_{n-1}& & (1Y^1+& d_2{Y}^2{)}^{n-1}\\ & & +& c_n& & & (1Y^1{)}^n\\ & & )\end{array}\}\text{ (RF)}$

ausrechnen mit dem Operator $\left[Y^n\right]$ für Koeffizientenextraktion. Da die Formel $\text{(RF)}$ auf ihrer rechten Seite nur Koeffizienten $d_j$ mit Indizes $j<n$ enthält, stellt sie eine rekursive Spezifikation der $d_n(n\ge 2)$ dar.

Bemerkung

Da die Formel $\text{(RF)}$ nur Ringoperationen (nur Additionen und Multiplikationen und keine Division) enthält, sind die Koeffizienten $d_n$ ganzzahlige Polynome in den $c_k$; das hat zur Folge, dass $\text{(RF)}$ über allen kommutativen unitären Ringen unabhängig von der Charakteristik und also gewissermaßen universell gültig ist.

Eine Herleitung der expliziten Auflösung

$d_n$	$=$	$\sum (-1{)}^{i_2+i_3+\cdots }$	$\frac{(n-1+i_2+i_3+\cdots )!}{n!i_2!i_3!\cdots }$	$c_2^{i_2}{c}_3^{i_3}\cdots$
	$=$	$\sum (-1{)}^{\sum _{\nu =2}^{\mathrm{\infty }}{i}_{\nu }}$	$\frac{(n-1+\sum _{\nu =2}^{\mathrm{\infty }}{i}_{\nu })!}{n!\prod _{\nu =2}^{\mathrm{\infty }}{i}_{\nu }!}$	$\prod _{\nu =2}^{\mathrm{\infty }}{c}_{\nu }^{i_{\nu }}$ ,

bei der über alle Kombinationen $i_2,i_3,\dots \in {\mathbb{N}}_0$ mit   $1i_2+2i_3+\cdots ={\sum }_{\nu =2}^{\mathrm{\infty }}(\nu -1)i_{\nu }=n-1$   zu summieren ist, findet sich bei Morse und Feshbach.

Die ersten paar Koeffizienten von $D$ sind:

$d_1$	$=1$
$d_2$	$=-c_2$
$d_3$	$=-c_3$	$+2c_2^2$			$(=-\frac{(2+1+0)!}{3!1!0!}{c}_3^1{c}_2^0+\frac{(2+0+2)!}{3!0!2!}{c}_3^0{c}_2^2)$
$d_4$	$=-c_4$	$+5c_3{c}_2$	$-5c_2^3$		$(=-\frac{(3+1+0+0)!}{4!1!0!0!}{c}_4^1{c}_3^0{c}_2^0+\frac{(3+0+1+1)!}{4!0!1!1!}{c}_4^0{c}_3^1{c}_2^1-\frac{(3+0+0+3)!}{4!0!0!3!}{c}_4^0{c}_3^0{c}_2^3)$
$d_5$	$=-c_5$	$+6c_4{c}_2$	$+3c_3^2$	$-21c_3{c}_2^2$	$+14c_2^4$
$d_6$	$=-c_6$	$+7c_5{c}_2$	$+7c_4{c}_3$	$-28c_4{c}_2^2$	$-28c_3^2{c}_2$	$+84c_3{c}_2^3$	$-42c_2^5$
$d_7$	$=-c_7$	$+8c_6{c}_2$	$+8c_5{c}_3$	$-36c_5{c}_2^2$	$+4c_4^2$	$-72c_4{c}_3{c}_2$	$+120c_4{c}_2^3$	$-12c_3^3$	$+180c_3^2{c}_2^2$	$-330c_3{c}_2^4$	$+132c_2^6$

Die Monome sind hier in den Zeilen lexikographisch absteigend geordnet, d. h., $c_3^2=c_3{c}_3$ kommt vor $c_3{c}_2$ kommt vor $c_3$ kommt vor Vorlage:Nowrap Die (ganzzahligen) Koeffizienten dieser Polynome sind in dieser Anordnung zusammengestellt in der Vorlage:OEIS. Die Vorlage:OEIS enthält die Anzahl der Monome in der $n$-ten Zeile (= Anzahl der Partitionen einer Vorlage:Nowrap Menge).

Mit der Substitution $D(X):=a_{1}B(X)$ ergibt sich

$X=C(D(X))=A(a_1^{-1}D(X))=A(B(X))$,

so dass $B(X)=a_1^{-1}D(X)$ die gesuchte Umkehrfunktion von $A(X)$ ist. Sie hat die Koeffizienten

$b_k=a_1^{-1}{d}_k$,

die allesamt ganzzahlige Polynome in $a_1^{-1}$ und den $a_n$ ($n\ge 2$) sind.

Formel von Lagrange-Bürmann

Ein Sonderfall der Lagrangeschen Inversionsformel, die in der Kombinatorik benutzt wird, gilt für $f(w)=w/\varphi (w)$ mit analytischem $\varphi (w)$ und $\varphi (0)\ne 0.$ Durch die Setzung $a=0$ wird $f(a)=f(0)=0.$ Dann ist für die Inverse $g(z):=f^{-1}(z)$

$\begin{array}{cc}g(z)& =\sum _{n=1}^{\mathrm{\infty }}\left(\underset{w\to 0}{\lim }\left(\frac{{\mathrm{d}}^{n-1}}{\mathrm{d}{w}^{n-1}}{\left(\frac{w}{w/\varphi (w)}\right)}^n\right)\frac{z^n}{n!}\right)\\ & =\sum _{n=1}^{\mathrm{\infty }}\frac{1}{n}\left(\frac{1}{(n-1)!}\underset{w\to 0}{\lim }(\frac{{\mathrm{d}}^{n-1}}{\mathrm{d}{w}^{n-1}}\varphi (w{)}^n)\right)z^n,\end{array}$

welches auch als

$[z^n]g(z)=\frac{1}{n}[w^{n-1}]\varphi (w{)}^n$

geschrieben werden kann mit dem Operator $[w^r]$, der den Koeffizienten des Terms $w^r$ in der rechts davon stehenden formalen Potenzreihe in $w$ extrahiert.

Eine nützliche Verallgemeinerung ist bekannt als Formel von Lagrange-Bürmann:

$[z^n]H(g(z))=\frac{1}{n}[w^{n-1}](H^{\prime }(w)\varphi (w{)}^n)$

mit einer beliebigen analytischen Funktion $H$.

Die Ableitung $H^{\prime }(w)$ kann eine sehr komplizierte Form annehmen, wann es durch $H(w)(1-{\varphi }^{\prime }(w)/\varphi (w))$ ersetzt werden kann, um

$[z^n]H(g(z))=[w^n]H(w)\varphi (w{)}^{n-1}(\varphi (w)-w{\varphi }^{\prime }(w))$

zu erhalten, welches auf ${\varphi }^{\prime }(w)$ anstelle von $H^{\prime }(w)$ Bezug nimmt.

Anwendungen

Die Lambertsche W-Funktion

Die Lambertsche W-Funktion ist die durch die implizite Gleichung

$W(z)e^{W(z)}=z$

definierte Funktion $W(z)$ .

Mithilfe der Lagrangesche Inversionsformel errechnet man für die Taylor-Reihe von $W(z)$ am Punkt $z=0$ wegen $f(w)=w{\mathrm{e}}^w$ und $a=b=0$ zuerst

$\frac{{\mathrm{d}}^n}{\mathrm{d}{x}^n}{\mathrm{e}}^{\alpha x}={\alpha }^n{\mathrm{e}}^{\alpha x}$ ,

woraus

$W(z)=\sum _{n=1}^{\mathrm{\infty }}\underset{w\to 0}{\lim }\left(\frac{{\mathrm{d}}^{n-1}}{\mathrm{d}{w}^{n-1}}{\mathrm{e}}^{-nw}\right)\frac{z^n}{n!}=\sum _{n=1}^{\mathrm{\infty }}(-n{)}^{n-1}\frac{z^n}{n!}=x-x^2+\frac{3}{2}{x}^3-\frac{8}{3}{x}^4+\frac{125}{24}{x}^5-\cdots .$

Der Konvergenzradius dieser Reihe ist ${\mathrm{e}}^{-1}$ .

Einen größeren Konvergenzradius erhält man auf ähnliche Weise:
Die Funktion $f(z)=W(e^z)-1$ erfüllt die Gleichung

$1+f(z)+\ln (1+f(z))=z$ .

Entwickelt man $z+\ln (1+z)$ in eine Potenzreihe und invertiert, dann erhält man für $f(z+1)=W(e^{z+1})-1$ :

$W(e^{1+z})=1+\frac{z}{2}+\frac{z^2}{16}-\frac{z^3}{192}-\frac{z^4}{3072}+\frac{13z^5}{61440}-\frac{47z^6}{1474560}-\frac{73z^7}{41287680}+\frac{2447z^8}{1321205760}+𝒪(x^9).$

Man kann daraus $W(x)$ ableiten, indem man $\ln x-1$ durch $z$ in dieser Reihe substituiert. Bspw. findet man $W(1)=0,567143\dots$ bei Vorlage:Nowrap

Binärbäume

Sei $ℬ$ die Menge der Binärbäume mit NIL-Knoten. Ein Baum aus $ℬ$ ist entweder ein NIL-Knoten oder ein Knoten mit zwei Teilbäumen.

Die Anzahl solcher Binärbäume mit $n$ (echten) Knoten sei mit $B_n$ bezeichnet.

Die Entfernung der Wurzel spaltet den Binärbaum in zwei kleinere Teilbäume. Daraus folgt für die erzeugende Funktion $B(z)=\sum _{n=0}^{\mathrm{\infty }}{B}_n{z}^n$:

$B(z)=1+zB(z{)}^2.$

Nun sei $C(z)=B(z)-1$, und damit $C(z)=z(C(z)+1{)}^2$ .

Die Anwendung der Lagrangeschen Inversionsformel mit $\varphi (w)=(w+1{)}^2$ ergibt:

$B_n=[z^n]C(z)=\frac{1}{n}[w^{n-1}](w+1{)}^{2n}=\frac{1}{n}(\genfrac{}{}{0ex}{}{2n}{n-1})=\frac{1}{n+1}(\genfrac{}{}{0ex}{}{2n}{n}).$

und das ist die $n$-te Catalan-Zahl.

Siehe auch

• Formel von Faà di Bruno

Weblinks

• Vorlage:MathWorld
• Vorlage:MathWorld
• Vorlage:MathWorld
• en:Kepler's equation#Inverse Kepler equation Umkehrung der Kepler-Gleichung mithilfe der Lagrangeschen Inversionsformel