imported>B wik: /* Faulhaber-Polynome */ Übersichtlichere Formeln

2025-09-23T20:51:06Z

Faulhaber-Polynome: Übersichtlichere Formeln

Neue Seite

Die '''Faulhabersche Formel''' beschreibt, wie sich die Summe der ersten <math>n</math> <math>k</math>-ten [[Potenz (Mathematik)|Potenzen]] mit einem [[Polynom]] <math>P_{k+1}(n) </math> in <math>n</math> vom Grad <math>k+1</math> berechnen lässt.

:<math>\sum_{h=1}^n h^k = 1^k + 2^k + 3^k + \cdots + n^k=P_{k+1}(n) \quad \text{mit}\; k\in\N_0 \; \text{und} \; n\in\N</math>

Die Koeffizienten des Polynoms <math>P_{k+1}(n)</math> können dabei mit Hilfe der [[Bernoulli-Zahl]]en berechnet werden.

Der Name „Faulhabersche Formel“ geht auf [[Donald Knuth]] zurück, der sie nach [[Johannes Faulhaber]] benannte.

== Darstellung des Polynoms mittels der Bernoulli-Zahlen ==

Zur Berechnung der Koeffizienten dieses Polynoms werden die [[Bernoulli-Zahlen]] benötigt. Im Folgenden bezeichne <math>B_j</math> die <math>j</math>-te Bernoulli-Zahl erster Art und die <math>\overline{B}_j</math> die Bernoulli-Zahlen zweiter Art, dann sieht die Faulhabersche Formel wie folgt aus:

:<math>\begin{align}
\sum_{h=1}^n h^k &= \frac{1}{k+1} \sum_{j=0}^k (-1)^j {k+1 \choose j} B_j n^{k+1-j}
= \frac{1}{k+1} \sum_{j=0}^k {k+1 \choose j} \overline{B}_j n^{k+1-j} \\
&= \frac{n^k}{2} + \frac{1}{k+1} \sum_{i=0}^{\lfloor k/2 \rfloor} {k+1 \choose 2i} B_{2i} n^{k+1-2i}
\end{align}</math>

Wenn man statt der ersten <math>n</math> nur die ersten <math>n-1</math> Potenzen betrachtet, so kann man die Faulhabersche Formel auch „ohne die Ausnahme“ für <math>B_1</math> beschreiben und erhält
:<math>\sum_{h=0}^{n-1} h^k = \frac{1}{k+1} \sum_{j=0}^k {k+1 \choose j} B_j n^{k+1-j}</math>

== Rekursive Darstellung ==

Betrachtet man das Polynom <math>P_{k+1}(n)</math> als Funktion der reellen Veränderlichen <math>n</math> anstelle <math>n\in\N_0</math> ergibt sich aus der Darstellung mittels Bernoulli-Zahlen auch die folgende elementar beweisbare rekursive Darstellung in aufsteigender Reihenfolge mit den Bernoulli-Zahlen zweiter Art

:<math>\begin{align}
P_1(n) &= B_0 n = B_0^* n = n \\
P_{k+1}(n) &= k\int_0^n P_k(n')dn' + n B_k^* \quad\text{falls}\; k \geq 1
\end{align}</math>

und in absteigender Folge gemäß:
:<math>k P_{k}(n) = \frac{d}{dn}P_{k+1}(n) - B_k^* \quad\text{falls}\; k \geq 1</math>

== Explizite Darstellungen ==
:<math>
\begin{array}{lllr}
1^0 + 2^0 + 3^0 + \cdots + n^0 &= n \\
1^1 + 2^1 + 3^1 + \cdots + n^1 &= \tfrac12 n^2 + \tfrac12 n \qquad \text{(Gaußsche Summenformel)} \\
1^2 + 2^2 + 3^2 + \cdots + n^2 &= \tfrac13 n^3 + \tfrac12 n^2 + \tfrac16 n \qquad \text{(Quadratische Pyramidalzahl)} \\
1^3 + 2^3 + 3^3 + \cdots + n^3 &= \tfrac14 n^4 + \tfrac12 n^3 + \tfrac14 n^2 \qquad \text{(Quadrat der Summe der ersten n Zahlen)}\\
1^4 + 2^4 + 3^4 + \cdots + n^4 &= \tfrac15 n^5 + \tfrac12 n^4 + \tfrac13 n^3 - \tfrac1{30} n \\
1^5 + 2^5 + 3^5 + \cdots + n^5 &= \tfrac16 n^6 + \tfrac12 n^5 + \tfrac5{12} n^4 - \tfrac1{12} n^2 \\
1^6 + 2^6 + 3^6 + \cdots + n^6 &= \tfrac17 n^7 + \tfrac12 n^6 + \tfrac12 n^5 - \tfrac16 n^3 + \tfrac1{42} n \\
1^7 + 2^7 + 3^7 + \cdots + n^7 &= \tfrac18 n^8 + \tfrac12 n^7 + \tfrac7{12} n^6 - \tfrac7{24} n^4 + \tfrac1{12} n^2 \\
1^8 + 2^8 + 3^8 + \cdots + n^8 &= \tfrac19 n^9 + \tfrac12 n^8 + \tfrac23 n^7 - \tfrac7{15} n^5 + \tfrac29 n^3 - \tfrac1{30} n \\
1^9 + 2^9 + 3^9 + \cdots + n^9 &= \tfrac1{10} n^{10} + \tfrac12 n^9 + \tfrac34 n^8 - \tfrac7{10} n^6 + \tfrac12 n^4 - \tfrac3{20} n^2 \\
1^{10} + 2^{10} + 3^{10} +\cdots + n^{10} &= \tfrac1{11} n^{11} + \tfrac12 n^{10} + \tfrac56 n^9 - n^7 + n^5 - \tfrac12 n^3 + \tfrac5{66} n \\
1^{11} + 2^{11} + 3^{11}+ \cdots + n^{11} &= \tfrac1{12} n^{12} + \tfrac12 n^{11} + \tfrac{11}{12} n^{10} - \tfrac{11}8 n^8 + \tfrac{11}6 n^6 - \tfrac{11}8 n^4 + \tfrac5{12} n^2 \\
1^{12} + 2^{12} + 3^{12} + \cdots + n^{12} &= \tfrac1{13} n^{13} + \tfrac12 n^{12} + n^{11} - \tfrac{11}6 n^9 + \tfrac{22}7 n^7 - \tfrac{33}{10} n^5 + \tfrac53 n^3 - \tfrac{691}{2730} n \\
\end{array}</math>

Die niedrigen Koeffizienten in Stammbrüchen, wie man sie bei kleinem <math>k</math> aus der Schulmathematik kennt, sind aber für den weiteren Verlauf überhaupt nicht typisch. Bereits bei <math>k=11</math> tritt zum ersten Mal ein Koeffizient > 1 auf; bei noch höheren Potenzen wird das zur Regel. Grund dafür sind die [[Bernoulli-Zahlen]], die nach einer Reihe von niedrigen Werten stark ansteigen, sogar stärker als jede Exponentialfunktion, und gegen Unendlich gehen. Sie selbst bilden die Koeffizienten der linearen Glieder, und da sie bei ungeraden Exponenten ungleich 1 Null werden, fehlen diese Glieder auch dementsprechend in den Summenformeln.

Allgemein gilt:
:<math>1^k + 2^k + 3^k + \dotsb + n^k = \tfrac1{k+1} n^{k+1} + \tfrac12 n^k + \mathcal{O}(n^{k-1}) </math>

(Das <math>\mathcal O</math> bezeichnet die [[Landau-Symbole|O-Notation]].) Hier sieht man auch den Zusammenhang mit [[Bonaventura Cavalieri|Cavalieris]] Integralformel; eine Summe von Potenzen ist eine Potenz mit einem um 1 höheren Grad. Das gilt auch für den trivialen Sonderfall von <math>k=0</math>, denn das Integral einer konstanten Funktion ist eine lineare.

Bei der Erweiterung von <math>k</math> auf <math>\Z</math> erhält man zunächst bei <math>k=-1</math> die divergente [[harmonische Reihe]], aber bei allen <math>k \leq -2</math> konvergente Potenzsummen. Ihre Grenzwerte sind definitionsgemäß die Funktionswerte der [[Riemannsche Zeta-Funktion|Riemannschen Zeta-Funktion]].

All das sind Spezialfälle der allgemeinen [[Euler-Maclaurin-Formel]] angewandt auf die Funktion <math>x^k</math> mit beliebigem reellem Exponenten <math>k</math>.

== Zusammenhang mit Bernoulli-Polynomen ==
Die Summe der ersten <math>n</math> <math>k</math>-ten Potenzen lässt sich auch mit Hilfe von [[Bernoulli-Polynome]]n ausdrücken:
:<math>\sum_{h=0}^{n} h^k = \frac{\text{B}_{k+1}(n+1)-\text{B}_{k+1}(0)}{k+1},</math>

Hierbei bezeichnet <math>\text{B}_j</math> das <math>j</math>-te Bernoulli-Polynom.

== Faulhaber-Polynome ==
Die Summen ungerader Potenzen
:<math>\sum_{h=1}^n h^{2k+1} = 1^{2k+1} + 2^{2k+1} + 3^{2k+1} + \cdots + n^{2k+1}\qquad (k,n\in\N_0)</math>

lassen sich auch als Polynom in <math> m=1+2+\dotsb+n =\tfrac 12 n(n+1)</math> darstellen. Solche Polynome in <math> m</math> statt in <math>n</math> werden auch als ''Faulhaber-Polynome'' bezeichnet.
Johannes Faulhaber selbst kannte nur die Formel in der folgenden beschriebenen Form und berechnete lediglich die ungeraden Fälle <math>k=1,3,5,\dotsc,17 </math> als Polynom in <math>m </math> und vermutete, dass für alle ungeraden Zahlen <math>k</math> ein entsprechendes Polynom existiere, ohne jedoch einen Beweis dafür zu geben. Das Konzept der Bernoulli-Zahlen war ihm nicht bekannt.

Einige Beispiele für kleinen Exponenten:

:<math>1^3 + 2^3 + 3^3 + \cdots + n^3 = m^2</math>     ({{OEIS|A000537}}, [[Satz von Nikomachos]])

:<math>1^5 + 2^5 + 3^5 + \cdots + n^5 = \frac{1}{3}(4m^3 - m^2)</math>     ({{OEIS|A000539}})

:<math>1^7 + 2^7 + 3^7 + \cdots + n^7 = \frac{1}{3}(6m^4 -4m^3 + m^2)</math>     ({{OEIS|A000541}})

:<math>1^9 + 2^9 + 3^9 + \cdots + n^9 = \frac{1}{5}(16m^5 - 20m^4 +12m^3 - 3m^2)</math>     ({{OEIS|A007487}})

:<math>1^{11} + 2^{11} + 3^{11} + \cdots + n^{11} = \frac{1}{3}(16m^6 - 32m^5 + 34m^4 - 20m^3 + 5m^2)</math>     ({{OEIS|A123095}})

Allgemein gilt für alle <math>k \in \N</math>:
:<math>
1^{2k-1} + 2^{2k-1} + 3^{2k-1} + \cdots + n^{2k-1} = \frac{1}{2^{2k}(2k)} \sum_{h=0}^{k-1} B_{2h} \binom{2k}{2h}
(2-2^{2h})\left((8m+1)^{k-h}-1\right)
</math>
was ein Polynom vom Grad <math>k</math> in <math>8m+1</math> darstellt oder explizit als Polynom in <math>m</math>
:<math> \sum_{h=1}^n h^{2k-1} = \frac{1}{2k} \sum_{i=1}^{k} a_{i,k}\, m^{i}
\qquad \text{mit} \quad
a_{i,k} = (-2)^{i} \sum_{j=0}^{i} {2k \choose i-j} {i+j \choose j} \frac{i-j}{i+j} B_{2k-i+j}
</math>

== Historisches ==
Faulhaber selbst kannte die Formel in ihrer heutigen allgemeinen Form nicht, auch waren die Bernoullizahlen zu seiner Zeit noch nicht bekannt. Er kannte jedoch zumindest die ersten 17 Fälle und die Konstruktionen der nach ihm benannten Polynome.<ref name="Knuth">[[Donald Knuth|Donald E. Knuth]]: ''Johann Faulhaber and Sums of Powers''. Math. Comp. 61 (1993), no. 203, S. 277–294. {{arXiv|math/9207222}}.</ref> Im Jahre 1834 veröffentlichte [[Carl Gustav Jacob Jacobi]] den ersten bekannten Beweis der Faulhaberschen Formel und verwendete dazu die [[Euler-Maclaurin-Formel]].<ref>Carl Gustav Jacob Jacobi: [http://visualiseur.bnf.fr/StatutConsulter?N=sorel2.bnf.fr-1296074706623&B=1&E=PDF&O=NUMM-77778 ''De usu legitimo formulae summatoriae Maclaurinianae.''] In: Crelles Journal für die reine und angewandte Mathematik 12 (1834), S. 263–272.</ref> Weitere Beweise wurden unter anderem 1923 von L. Tits und 1986 von A. W. F. Edwards publiziert.<ref>L. Tits: ''Sur la sommation des puissances numériques.'' In: Mathesis 37 (1923), S. 353–355.</ref><ref>Anthony William Fairbank Edwards: ''A quick route to sums of powers.'' In: American Mathematical Monthly 93 (1986), S. 451–455. {{JSTOR|2323466}}</ref> Donald Ervin Knuth untersuchte Verallgemeinerungen, Darstellungen der Summen als Polynome in <math>n(n+1)\cdots(n+r)</math> mit festem <math>r</math>,<ref name="Knuth" /> und trug zur Popularisierung der Faulhaberschen Polynome bei.

== Literatur ==
* [[Donald Knuth|Donald E. Knuth]]: ''Johann Faulhaber and Sums of Powers.'' Math. Comp. 61 (1993), no. 203, S. 277–294. {{arXiv|math/9207222}}.
* [[John H. Conway]], [[Richard Guy]]: ''The Book of Numbers.'' Copernicus (Springer), New York 1996, ISBN 0-387-97993-X, S. 106, ({{Google Buch|BuchID=0--3rcO7dMYC|Seite=106|Linktext=Auszug (Google)|KeinText=ja}}).
* [[Johann Faulhaber]]: ''Academia Algebrae. Darinnen die miraculosische Inventiones, zu den höchsten Cossen weiters continuirt und profitiert werden. Dergleichen zwar vor 15. Jahren den Gelehrten auff allen Vniversiteten in gantzem Europa proponiert, darauff continuiert, auch allen Mathematicis inn der gantzen weiten Welt dediciert, aber bißhero, noch nie so hoch, biß auff die regulierte Zensicubiccubic Coß, durch offnen Truck publiciert worden. Welcher vorgesetzet ein kurtz Bedencken, Was einer für Authores nach ordnung gebrauchen solle, welcher die Coß fruchtbarlich, bald, auch fundamentaliter lehrnen vnd ergreiffen will.'' Augsburg: Johann Ulrich Schönig, 1631. ({{Google Buch|BuchID=0pw_AAAAcAAJ|Seite=|Linktext=Online-Kopie (Google)|KeinText=ja}})
* {{Literatur |Autor=Greg Orosi |Titel=A Simple Derivation Of Faulhaber's Formula |Sammelwerk=Applied Mathematics E-Notes |Band=18 |Datum=2018 |Seiten=124–126 |Online=[http://www.math.nthu.edu.tw/~amen/2018/AMEN-170803.pdf Online]}}

== Weblinks ==
* {{MathWorld|id=FaulhabersFormula|title=Faulhaber's formula}}
* {{MacTutor|id=Faulhaber|title=Johann Faulhaber|page=bio}}
* Helmut Richter, Bernhard Schiekel: ''Potenzsummen, Bernoulli-Zahlen und Eulersche Summenformel.'' [[doi:10.18725/OPARU-1819]]
* Tobias Krähling: [http://www.semibyte.de/wp/download/maths/ermittlung_der_potenzsummenformel.pdf ''Aussagenlogik.''] (PDF; 148 kB) – Herleitung der Potenzsummenformel/des Polynoms

== Einzelnachweise ==
<references />

[[Kategorie:Diskrete Mathematik]]
[[Kategorie:Theorie der Polynome]]
[[Kategorie:Folgen und Reihen]]

Faulhabersche Formel - Versionsgeschichte

imported>B wik: /* Faulhaber-Polynome */ Übersichtlichere Formeln