Wallissches Produkt

Das wallissche Produkt, auch Wallis-Produkt, ist eine Produktdarstellung der Kreiszahl $\pi$, das heißt, es handelt sich um ein Produkt mit unendlich vielen Faktoren, dessen Grenzwert $\pi$ ist. Es wurde 1655 von dem englischen Mathematiker John Wallis entdeckt. Dazu nutzte er eine schachbrettartige 'Interpolation' zwischen den (in ganzen Dimensionen) figurierten Zahlenfolgen des Pascalschen Dreiecks zur Bestimmung von 4/Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \pi} als mittleren Binomialkoeffizienten zwischen nullter und erster Dimension. Im Jahr 2015 wurde erstmals ein Zusammenhang mit quantenmechanischen Berechnungen bezüglich des Wasserstoffatoms festgestellt.

Formel

Üblich ist die Darstellung des Produktes in der Form:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \frac{\pi}{2} = \frac{2}{1} \cdot \frac{2}{3} \cdot \frac{4}{3} \cdot \frac{4}{5} \cdot \frac{6}{5} \cdot \frac{6}{7} \cdot \dots}

Über eine Umformung ergibt sich die Kurzschreibweise des Wallisproduktes wie folgt:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \frac{\pi}{2}= \left( \frac{2}{1} \cdot \frac{2}{3} \right) \cdot \left( \frac{4}{3} \cdot \frac{4}{5} \right) \cdot \left( \frac{6}{5} \cdot \frac{6}{7} \right) \cdot\ \dots = \prod_{i=1}^{\infty} \frac{(2i)(2i)}{(2i-1)(2i+1)} = \prod_{i=1}^\infty \frac{4 \, i^2}{4 \, i^2 - 1} = \prod_{i=1}^\infty\left(1+\frac{1}{4i^2-1}\right)}

Für den Kehrwert folgt:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \frac 2{\pi} =\prod_{i=1}^\infty \bigg(1- \frac1{4 \, i^2}\bigg)}

Die Konvergenz dieses Produktes folgt aus der Konvergenz der unendlichen Reihe

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \sum_{i=1}^\infty \frac{-1}{4 \, i^2}} bzw. $\sum _{i=1}^{\infty }{\frac {1}{i^{2}}}$

Konvergenzgeschwindigkeit

n	2·Produkt	2·Produkt / Pi	relativer Fehler
1	2,7	0,85	15 %
2	2,8	0,91	09 %
3	2,9	0,93	07 %
10	3,07	0,976	02,4 %
100	3,134	0,9975	00,25 %
1000	3,1408	0,99975	00,025 %
10000	3,14151	0,999975	00,0025 %
100000	3,141585	0,9999975	00,00025 %
1000000	3,1415918	0,99999975	00,000025 %
Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \lim_{n\to\infty}}	3,14159265…	1,00000000…	00 %

Zur effizienten Berechnung einer Näherung von Pi ist die Formel nicht geeignet. Berechnet man etwa die ersten 5 Terme des Wallischen Produkts und verdoppelt das Ergebnis, so erhält man als Näherung für Pi:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle 2\cdot \prod_{i=1}^{5} \frac{4 \cdot i^2}{4 \cdot i^2 - 1} \approx 3{,}002}

Mit dieser Näherung konnte nicht einmal die erste Nachkommastelle korrekt bestimmt werden.

Nach Ausmultiplizieren der ersten 50 Terme ergibt sich ein Quotient aus zwei 160-stelligen Zahlen, der aber für Pi nur die Näherung 3,126 liefert, also nicht einmal 2 Nachkommastellen korrekt angibt. Da 3,126/3,14159 = 0,9950 ist, ist der relative Fehler etwa 0,5 %. Die Konvergenzgeschwindigkeit ist langsamer als linear.

Die nebenstehende Tabelle gibt für einige ausgewählte Werte von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle N} an, wie gut die Approximation von Pi ist, die man nach Ausmultiplizieren von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle N} Termen im wallisschen Produkt erhält. Die Tabelle legt die Vermutung nahe, dass der Fehler nach Ausmultiplizieren von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle N} Termen in etwa Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \tfrac{25}{N}\%} beträgt (z. B. nach 100 Termen: 0,25 % = ${\tfrac {1}{400}}$).

Dies kann man auch durch folgende mathematische Überlegung beweisen: Der Quotient zwischen der Approximation und dem gewünschten Wert ist gleich dem unendlichen Produkt

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \prod_{i=N+1}^\infty \bigg(1- \frac1{4 \, i^2}\bigg)}

Mit Hilfe der Rechenregeln für Logarithmen, der Abschätzung Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \log(1+x) \approx x } (für kleine Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle x} ) sowie durch Approximation einer unendlichen Summe durch ein Integral sieht man, dass dieses Produkt ungefähr den folgenden Wert hat:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \exp\left(\int\limits_N^\infty \log\left(1- \frac1{4 \, x^2}\right) \,dx\right) \approx \exp\left(\int\limits_N^\infty -\frac1{4 \, x^2}\,dx\right) \approx \exp\left(-\frac1{4 \, N}\right) \approx 1-\frac1{4 \, N}} .

Damit die ersten beiden Nachkommastellen richtig sind, braucht man demzufolge eine Genauigkeit von ca. 0,3 % (3,13/3,14 = 0,997), also etwa Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle N=60.} Für 3 Nachkommastellen braucht man Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle N=600,} für 4 Nachkommastellen $N=6000$ etc.

Beweisskizze

Man definiert Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle C_n(x) := \int\limits_0^x \sin^n(t) dt} , für welche die Rekursionsformel Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle (n+1)C_{n+1}(x) = -\cos(x)\sin^n(x) + nC_{n-1}(x)} gilt. Insbesondere erhält man für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle c_n := C_n(\pi/2)} die Formel Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle c_{n+1}=\frac{n}{n+1} c_{n-1}} .

Man berechnet Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle c_{2m} = \frac{\pi}{2} \frac{1}{4^m} \binom{2m}{m}} und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle c_{2m+1} = \frac{4^m}{2m+1} / \binom{2m}{m} } . Nun gilt ${\frac {n}{n+1}}={\frac {c_{n+1}}{c_{n}}}{\frac {c_{n}}{c_{n-1}}}<{\frac {c_{n+1}}{c_{n}}}<1$, und daher Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle 1 < \frac{c_{2m}}{c_{2m+1}} = \pi (m+\frac{1}{2}) \binom{2m}{m}^2 4^{-2m} < 1+\frac{1}{2m} } .

Insbesondere also Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): {\displaystyle \sqrt{\pi} = \lim_{m \rightarrow \infty} \frac{1}{\sqrt{m}} \frac{2 \cdot 4 \cdot \cdot \cdot (2m)}{1 \cdot 3 \cdot \cdot \cdot (2m-1)} = \lim_{m \to \infty} \frac{1}{\sqrt{m}} \frac{(2m)!!}{(2m-1)!!}} , aus der man durch quadrieren die übliche Formel erhält.

Physik

Von Tamar Friedmann, C. R. Hagen und Studenten der Uni. Rochester (USA) wurde 2015 eine Anwendung dieses Produkts bei der Berechnung des Fehlers der quantenmechanischen Variationsrechnung der Energieeigenzustände im angeregten Wasserstoffatom relativ zur Lösung im Bohr’schen Atommodell entdeckt.

Literatur

• John Wallis: The arithmetic of infinitesimals (Übersetzung vom Latein ins Englische mit einem Vorwort von Jacqueline A. Stedall). 1. Auflage. Springer Verlag, Heidelberg / Berlin / New York 2004, ISBN 0-387-20709-0.
• Kurt Endl, Wolfgang Luh: Analysis. Eine integrierte Darstellung; Studienbuch für Studierende der Mathematik, Physik und anderer Naturwissenschaften ab 1. Semester. Band 2. 7. überarbeitete Auflage. Aula-Verlag, Wiesbaden 1989, ISBN 3-89104-455-0, S. 343.
• Fridtjof Toenniessen: Das Geheimnis der transzendenten Zahlen: Eine etwas andere Einführung in die Mathematik. Springer 2009, ISBN 978-3-8274-2274-3, S. 321–322 (Auszug (Google)).

Weblinks

Wikibooks

Wikibooks: Herleitung des Wallis-Produktes – Lern- und Lehrmaterialien

• Johan Wästlund: An Elementary Proof of the Wallis Product Formula for Pi. (PDF; 59 kB) MAA Monthly, Dezember 2007
• Eric W. Weisstein: Wallis Formula. In: MathWorld (englisch).