Partialbruchzerlegung

Die Partialbruchzerlegung oder Partialbruchentwicklung ist eine standardisierte Darstellung rationaler Funktionen. Sie wird in der Mathematik verwendet, um das Rechnen mit solchen Funktionen zu erleichtern. Insbesondere kommt sie bei der Integration der rationalen Funktionen zur Anwendung.

Hier liegt die Tatsache zugrunde, dass jede rationale Funktion als Summe einer Polynomfunktion und Brüchen der Form

a i ( x x i ) j {\displaystyle {\frac {a_{i}}{(x-x_{i})^{j}}}}

dargestellt werden kann. Die x i {\displaystyle x_{i}} sind dabei die Polstellen der Funktion.

Sind die Polstellen bereits bekannt, so ist die Bestimmung der Zähler a i {\displaystyle a_{i}} die eigentliche Aufgabe der Partialbruchzerlegung.

Bei reellwertigen Funktionen müssen die Polstellen x i {\displaystyle x_{i}} und infolgedessen auch die Zahlen a i {\displaystyle a_{i}} nicht unbedingt reell sein, denn die reellen Zahlen sind nicht algebraisch abgeschlossen. Man kann das Rechnen mit komplexen Zahlen aber vermeiden, weil mit jeder komplexen Nullstelle z i {\displaystyle z_{i}} auch die konjugiert komplexe Zahl z i ¯ {\displaystyle {\overline {z_{i}}}} Nullstelle ist.

Statt a 1 ( x z i ) j {\displaystyle {\tfrac {a_{1}}{(x-z_{i})^{j}}}} und a 2 ( x z i ¯ ) j {\displaystyle {\tfrac {a_{2}}{(x-{\overline {z_{i}}})^{j}}}} verwendet man dann einen Term b i + c i x ( x 2 + p i x + q i ) j {\displaystyle {\tfrac {b_{i}+c_{i}x}{(x^{2}+p_{i}x+q_{i})^{j}}}} , wobei x 2 + p i x + q i = ( x z i ) ( x z i ¯ ) {\displaystyle {x^{2}+p_{i}x+q_{i}}=(x-z_{i})\cdot (x-{\overline {z_{i}}})} ein reelles quadratisches Polynom ist und auch b i {\displaystyle b_{i}} und c i {\displaystyle c_{i}} reell sind.

Geschichte

Die Partialbruchzerlegung wurde ab 1702 in Arbeiten zur Infinitesimalrechnung von Gottfried Wilhelm Leibniz und Johann I Bernoulli entwickelt. Beide Gelehrten nutzten diese Methode zur Integration von gebrochenrationalen Funktionen. Da zu dieser Zeit der Fundamentalsatz der Algebra noch nicht bewiesen war – er wurde damals aber schon vermutet –, behauptete Leibniz, dass es für das Nennerpolynom x 4 + a 4 = ( x x a a 1 ) ( x x + a a 1 ) {\displaystyle x^{4}+a^{4}=\left(xx-aa{\sqrt {-1}}\right)\left(xx+aa{\sqrt {-1}}\right)} keine Partialbruchzerlegung gebe. Johann Bernoulli schloss sich dieser Meinung nicht an. Dieses Beispiel wurde in den Folgejahren von verschiedenen Mathematikern diskutiert und um 1720 erschienen mehrere Arbeiten, die das Beispiel als fehlerhaft nachwiesen und das (unbestimmte) Integral

d x x 4 + a 4 {\displaystyle \int {\frac {\mathrm {d} x}{x^{4}+a^{4}}}}

korrekt berechneten.[1]

Verfahren

Die Partialbruchzerlegung einer reellen rationalen Funktion R {\displaystyle R} wird in mehreren Schritten bestimmt:

  1. Man vergleicht den Grad des Zählers mit dem des Nenners von R {\displaystyle R} :
    • Ist der Zählergrad größer oder gleich dem Nennergrad, so dividiert man den Zähler durch den Nenner. Man erhält daraus das Polynom P {\displaystyle P} und möglicherweise eine rationale Restfunktion R = Z N {\displaystyle R^{*}={\tfrac {Z^{*}}{N^{*}}}} , sodass gilt: R ( x ) = P ( x ) + R ( x ) {\displaystyle R(x)=P(x)+R^{*}(x)} .
      • Ist R 0 {\displaystyle R^{*}\equiv 0} , ist das Verfahren abgeschlossen.
      • Andernfalls hat der Zähler Z {\displaystyle Z^{*}} von R {\displaystyle R^{*}} einen kleineren Grad als der Nenner N {\displaystyle N^{*}} . Man arbeitet dann nur mehr mit der Restfunktion R {\displaystyle R^{*}} weiter.
    • Ist der Zählergrad kleiner als der Nennergrad, so kann man die Funktion R {\displaystyle R} direkt betrachten. Um im Folgenden eine einheitliche Bezeichnungsweise zu ermöglichen, setzen wir in diesem Fall R := R {\displaystyle R^{*}:=R} .
  2. Anschließend betrachtet man die Nullstellen von N {\displaystyle N^{*}} . Abhängig von der Art der Nullstellen wird ein geeigneter Ansatz verwendet.
  3. Die Konstanten a i j {\displaystyle a_{ij}} , b i j {\displaystyle b_{ij}} und c i j {\displaystyle c_{ij}} erhält man dann zum Beispiel durch Koeffizientenvergleich nach Multiplikation der Zerlegung mit dem Nennerpolynom.

Die beiden letzten Schritte sollen nun im Detail erläutert werden.

Ansatz

Vorausgesetzt wird hier, dass R {\displaystyle R^{*}} in der Form R ( x ) = Z ( x ) N ( x ) {\displaystyle R^{*}(x)={\tfrac {Z^{*}(x)}{N^{*}(x)}}} gegeben ist, wobei der Grad von Z {\displaystyle Z^{*}} kleiner als der Grad des Nennerpolynoms N {\displaystyle N^{*}} ist und sämtliche Nullstellen von N {\displaystyle N^{*}} bekannt sind. Sind, wie oben angenommen, die n {\displaystyle n} verschiedenen Nullstellen x i {\displaystyle x_{i}} und ihr jeweiliger Grad r i {\displaystyle r_{i}} bekannt, so kann das Nennerpolynom auf folgende Form gebracht werden:

N ( x ) = ( x x 1 ) r 1 ( x x 2 ) r 2 ( x x n ) r n {\displaystyle N^{*}(x)=(x-x_{1})^{r_{1}}\cdot (x-x_{2})^{r_{2}}\dotsm (x-x_{n})^{r_{n}}}

Zu beachten ist, dass einige der x i {\displaystyle x_{i}} nicht-reell sein können.

Der Ansatz ist nun folgendermaßen aufgebaut:

  • Für jede einfache reelle Nullstelle x i {\displaystyle x_{i}} enthält der Ansatz einen Summanden a i 1 x x i {\displaystyle {\tfrac {a_{i1}}{x-x_{i}}}} .
  • Für jede r i {\displaystyle r_{i}} -fache reelle Nullstelle x i {\displaystyle x_{i}} enthält der Ansatz r i {\displaystyle r_{i}} Summanden a i 1 x x i + a i 2 ( x x i ) 2 + + a i r i ( x x i ) r i {\displaystyle {\tfrac {a_{i1}}{x-x_{i}}}+{\tfrac {a_{i2}}{(x-x_{i})^{2}}}+\dotsb +{\tfrac {a_{ir_{i}}}{(x-x_{i})^{r_{i}}}}} .

Da R {\displaystyle R} reell ist, gehört zu jeder nicht-reellen Nullstelle z i {\displaystyle z_{i}} notwendigerweise auch die konjugiert komplexe Nullstelle z i ¯ {\displaystyle {\overline {z_{i}}}} . Sei x 2 + p i x + q i {\displaystyle x^{2}+p_{i}x+q_{i}} das quadratische Polynom mit den Nullstellen z i {\displaystyle z_{i}} und z i ¯ {\displaystyle {\overline {z_{i}}}} , also x 2 + p i x + q i := ( x z i ) ( x z i ¯ ) {\displaystyle x^{2}+p_{i}x+q_{i}:=(x-z_{i})(x-{\overline {z_{i}}})} .

  • Für jede einfache nicht-reelle Nullstelle z i {\displaystyle z_{i}} enthält der Ansatz nun einen Summanden b i x + c i x 2 + p i x + q i {\displaystyle {\tfrac {b_{i}x+c_{i}}{x^{2}+p_{i}x+q_{i}}}} .
  • Entsprechend enthält der Ansatz für jede s i {\displaystyle s_{i}} -fache nicht-reelle Nullstelle z i {\displaystyle z_{i}} (und die zugehörige, ebenfalls s i {\displaystyle s_{i}} -fache, konjugiert komplexe Nullstelle z i ¯ {\displaystyle {\overline {z_{i}}}} ) die s i {\displaystyle s_{i}} Terme b i 1 x + c i 1 x 2 + p i x + q i + b i 2 x + c i 2 ( x 2 + p i x + q i ) 2 + + b i s i x + c i s i ( x 2 + p i x + q i ) s i {\displaystyle {\tfrac {b_{i1}x+c_{i1}}{x^{2}+p_{i}x+q_{i}}}+{\tfrac {b_{i2}x+c_{i2}}{(x^{2}+p_{i}x+q_{i})^{2}}}+\dotsb +{\tfrac {b_{is_{i}}x+c_{is_{i}}}{(x^{2}+p_{i}x+q_{i})^{s_{i}}}}} .

Jeder Ansatz enthält somit genau g {\displaystyle g} unbekannte Koeffizienten a i 1 , , a i r i , b i 1 , , b i s i , c i 1 , , c i s i {\displaystyle a_{i1},\dotsc ,a_{ir_{i}},b_{i1},\dots ,b_{is_{i}},c_{i1},\dotsc ,c_{is_{i}}} .

Bestimmung der Konstanten

Um die Konstanten a i j {\displaystyle a_{ij}} , b i j {\displaystyle b_{ij}} und c i j {\displaystyle c_{ij}} zu ermitteln, wird R {\displaystyle R^{*}} mit dem Ansatz gleichgesetzt und diese Gleichung mit dem Nennerpolynom N {\displaystyle N^{*}} multipliziert.

Auf der einen Seite der Gleichung steht dann nur noch das Zählerpolynom Z {\displaystyle Z^{*}} , auf der anderen ein Ausdruck mit allen Unbekannten, der ebenfalls ein Polynom in x {\displaystyle x} ist und entsprechend nach den Potenzen von x {\displaystyle x} geordnet werden kann. Ein Koeffizientenvergleich der linken mit der rechten Seite ergibt dann ein lineares Gleichungssystem, aus dem sich die unbekannten Konstanten berechnen lassen. Alternativ kann man bis zu g {\displaystyle g} beliebige verschiedene Werte für x {\displaystyle x} in diese Gleichung einsetzen, was wie der Koeffizientenvergleich zu einem aus g {\displaystyle g} Gleichungen bestehenden linearen Gleichungssystem führt. Sinnvoll ist das Einsetzen der zuvor berechneten (reellen) Nullstellen, was sofort jeweils einen Koeffizientenwert liefert.

Diese beiden Möglichkeiten können auch kombiniert werden.

Beispiele

Einfache Polstellen

Gegeben sei die rationale Funktion

R ( x ) = x x 2 1 {\displaystyle R(x)={\frac {x}{x^{2}-1}}} .

Es gibt zwei einfache Polstellen x 1 = 1 {\displaystyle x_{1}=1} und x 2 = 1 {\displaystyle x_{2}=-1} . Der Ansatz lautet also

x x 2 1 = a 1 x 1 + a 2 x + 1 {\displaystyle {\frac {x}{x^{2}-1}}={\frac {a_{1}}{x-1}}+{\frac {a_{2}}{x+1}}} ,

wobei a 1 {\displaystyle a_{1}} und a 2 {\displaystyle a_{2}} unbekannte, noch zu ermittelnde Konstanten sind. Multipliziert man beide Seiten der Gleichung mit ( x 2 1 ) = ( x + 1 ) ( x 1 ) {\displaystyle (x^{2}-1)=(x+1)(x-1)} , erhält man

x = a 1 ( x + 1 ) + a 2 ( x 1 ) {\displaystyle x=a_{1}(x+1)+a_{2}(x-1)} .

Sortiert man die rechte Seite nach Gliedern mit x {\displaystyle x} und Gliedern ohne x {\displaystyle x} , so ergibt sich

x = ( a 1 + a 2 ) x + ( a 1 a 2 ) {\displaystyle x=(a_{1}+a_{2})x+(a_{1}-a_{2})} .

Koeffizientenvergleich: Der Koeffizient von x {\displaystyle x} ist Eins: a 1 + a 2 = 1 {\displaystyle a_{1}+a_{2}=1} und das absolute Glied Null: a 1 a 2 = 0 {\displaystyle a_{1}-a_{2}=0} . Hieraus lässt sich berechnen: a 1 = a 2 = 1 2 . {\displaystyle a_{1}=a_{2}={\tfrac {1}{2}}.} Die gesuchte Partialbruchzerlegung ist also

x x 2 1 = 1 2 1 x 1 + 1 2 1 x + 1 . {\displaystyle {\frac {x}{x^{2}-1}}={\frac {1}{2}}\cdot {\frac {1}{x-1}}+{\frac {1}{2}}\cdot {\frac {1}{x+1}}.}

Doppelte Polstellen

Gegeben sei die rationale Funktion

R ( x ) = x 2 x 2 2 x + 1 {\displaystyle R(x)={\frac {x^{2}}{x^{2}-2x+1}}} .

Mittels Polynomdivision und Faktorenzerlegung des Nenners folgt

R ( x ) = x 2 x 2 2 x + 1 = 1 + 2 x 1 ( x 1 ) 2 {\displaystyle R(x)={\frac {x^{2}}{x^{2}-2x+1}}=1+{\frac {2\;x-1}{(x-1)^{2}}}} .

Die einzige, allerdings doppelte Nullstelle des Nenners ist x 0 = 1 {\displaystyle x_{0}=1} . Ansatz:

2 x 1 ( x 1 ) 2 = a 1 x 1 + a 2 ( x 1 ) 2 | ( x 1 ) 2 {\displaystyle {\frac {2\;x-1}{(x-1)^{2}}}={\frac {a_{1}}{x-1}}+{\frac {a_{2}}{(x-1)^{2}}}\quad |\cdot (x-1)^{2}}
2 x 1 = a 1 ( x 1 ) + a 2 {\displaystyle 2\;x-1=a_{1}(x-1)+a_{2}}
2 x 1 = a 1 x a 1 + a 2 {\displaystyle 2\;x-1=a_{1}x-a_{1}+a_{2}}

Koeffizientenvergleich:

a 1 = 2 {\displaystyle a_{1}=2}
a 1 + a 2 = 1 {\displaystyle -a_{1}+a_{2}=-1}

Lösung:

a 1 = 2 , a 2 = 1 {\displaystyle a_{1}=2,\quad a_{2}=1} ,

also erhalten wir die Partialbruchzerlegung

x 2 x 2 2 x + 1 = 1 + 2 x 1 + 1 ( x 1 ) 2 {\displaystyle {\frac {x^{2}}{x^{2}-2x+1}}=1+{\frac {2}{x-1}}+{\frac {1}{(x-1)^{2}}}} .

Komplexe Polstellen

Gegeben sei die rationale Funktion

R ( x ) = 5 x 2 + 2 x + 1 x 3 + x {\displaystyle R(x)={\frac {5x^{2}+2x+1}{x^{3}+x}}} .

Der Nenner hat hier die reelle Nullstelle x 1 = 0 {\displaystyle x_{1}=0} , die komplexe Nullstelle x 2 = z 1 = i {\displaystyle x_{2}=z_{1}=\mathrm {i} } und deren konjugiert komplexe x 3 = z 1 ¯ = i {\displaystyle x_{3}={\overline {z_{1}}}=-\mathrm {i} } . Das quadratische Polynom mit den Nullstellen z 1 {\displaystyle z_{1}} und z 1 ¯ {\displaystyle {\overline {z_{1}}}} ist ( x z 1 ) ( x z 1 ¯ ) = ( x i ) ( x + i ) = x 2 + 1 {\displaystyle (x-z_{1})(x-{\overline {z_{1}}})=(x-\mathrm {i} )(x+\mathrm {i} )=x^{2}+1}

Ansatz:

5 x 2 + 2 x + 1 x 3 + x = a 1 x + b 1 x + c 1 x 2 + 1 {\displaystyle {\frac {5x^{2}+2x+1}{x^{3}+x}}={\frac {a_{1}}{x}}+{\frac {b_{1}x+c_{1}}{x^{2}+1}}}
5 x 2 + 2 x + 1 = a 1 x 2 + a 1 + b 1 x 2 + c 1 x {\displaystyle 5x^{2}+2x+1=a_{1}x^{2}+a_{1}+b_{1}x^{2}+c_{1}x}
5 x 2 + 2 x + 1 = ( a 1 + b 1 ) x 2 + c 1 x + a 1 {\displaystyle 5x^{2}+2x+1=(a_{1}+b_{1})x^{2}+c_{1}x+a_{1}}

Koeffizientenvergleich:

a 1 + b 1 = 5 c 1 = 2 a 1 = 1 {\displaystyle {\begin{aligned}a_{1}+b_{1}&=5\\c_{1}&=2\\a_{1}&=1\end{aligned}}}

Lösung:

a 1 = 1 , b 1 = 4 , c 1 = 2 {\displaystyle a_{1}=1,\quad b_{1}=4,\quad c_{1}=2} ,

Partialbruchzerlegung:

5 x 2 + 2 x + 1 x 3 + x = 1 x + 4 x + 2 x 2 + 1 {\displaystyle {\frac {5x^{2}+2x+1}{x^{3}+x}}={\frac {1}{x}}+{\frac {4x+2}{x^{2}+1}}}

Kubische und quartische Nenner

Kubische Nenner:

Für Brüche mit kubischem Nenner gilt unter der Bedingung a²e + b²c − abd ≠ 0 folgende Partialbruchzerlegung:

t x 2 + u x + v ( a x + b ) ( c x 2 + d x + e ) = a 2 v + b 2 t a b u ( a 2 e + b 2 c a b d ) ( a x + b ) + [ a ( e t c v ) b ( d t c u ) ] x + a ( e u d v ) b ( e t c v ) ( a 2 e + b 2 c a b d ) ( c x 2 + d x + e ) {\displaystyle {\frac {tx^{2}+ux+v}{(ax+b)(cx^{2}+dx+e)}}={\frac {a^{2}v+b^{2}t-abu}{(a^{2}e+b^{2}c-abd)(ax+b)}}+{\frac {[a(et-cv)-b(dt-cu)]x+a(eu-dv)-b(et-cv)}{(a^{2}e+b^{2}c-abd)(cx^{2}+dx+e)}}}

Beispielsweise kann dieser Bruch mit der genannten Formel zerlegt werden:

1 x 3 + 1 = 1 ( x + 1 ) ( x 2 x + 1 ) = 1 3 ( x + 1 ) + x + 2 3 ( x 2 x + 1 ) = 1 3 ( x + 1 ) 2 x 1 6 ( x 2 x + 1 ) + 1 2 ( x 2 x + 1 ) {\displaystyle {\frac {1}{x^{3}+1}}={\frac {1}{(x+1)(x^{2}-x+1)}}={\frac {1}{3(x+1)}}+{\frac {-x+2}{3(x^{2}-x+1)}}={\frac {1}{3(x+1)}}-{\frac {2x-1}{6(x^{2}-x+1)}}+{\frac {1}{2(x^{2}-x+1)}}}

Hiermit kann ein kubisches Analogon zur Leibniz-Reihe ermittelt werden:

k = 0 ( 1 ) k 3 k + 1 = 0 1 1 x 3 + 1 d x = 0 1 ( 1 3 ( x + 1 ) 2 x 1 6 ( x 2 x + 1 ) + 1 2 ( x 2 x + 1 ) ) d x = {\displaystyle \sum _{k=0}^{\infty }{\frac {(-1)^{k}}{3k+1}}=\int _{0}^{1}{\frac {1}{x^{3}+1}}\mathrm {d} x=\int _{0}^{1}\left({\frac {1}{3(x+1)}}-{\frac {2x-1}{6(x^{2}-x+1)}}+{\frac {1}{2(x^{2}-x+1)}}\right)\mathrm {d} x=}
= [ 1 3 ln | x + 1 | 1 6 ln ( x 2 x + 1 ) + 1 3 3 arctan [ 1 3 3 ( 2 x 1 ) ] ] x = 0 x = 1 = 1 9 3 π + 1 3 ln ( 2 ) 0,835 648848264721 {\displaystyle ={\biggl [}{\frac {1}{3}}\ln |x+1|-{\frac {1}{6}}\ln(x^{2}-x+1)+{\frac {1}{3}}{\sqrt {3}}\arctan {\bigl [}{\frac {1}{3}}{\sqrt {3}}(2x-1){\bigr ]}{\biggr ]}_{x=0}^{x=1}={\frac {1}{9}}{\sqrt {3}}\pi +{\frac {1}{3}}\ln(2)\approx 0{,}835648848264721}

Quartische Nenner:

Die Partialbruchzerlegung von Brüchen mit quartischem Nenner kann mit einer Matrix ermittelt werden:

t x 3 + u x 2 + v x + w ( a x 2 + b x + c ) ( d x 2 + e x + f ) = y 1 x + y 2 a x 2 + b x + c + y 3 x + y 4 d x 2 + e x + f {\displaystyle {\frac {tx^{3}+ux^{2}+vx+w}{(ax^{2}+bx+c)(dx^{2}+ex+f)}}={\frac {y_{1}x+y_{2}}{ax^{2}+bx+c}}+{\frac {y_{3}x+y_{4}}{dx^{2}+ex+f}}}

Für diese Form muss folgendes Produkt von reziproker Matrix und Vektor ermittelt werden:

[ y 1 y 2 y 3 y 4 ] = [ d 0 a 0 e d b a f e c b 0 f 0 c ] 1 [ t u v w ] {\displaystyle {\begin{bmatrix}y_{1}\\y_{2}\\y_{3}\\y_{4}\end{bmatrix}}={\begin{bmatrix}d&0&a&0\\e&d&b&a\\f&e&c&b\\0&f&0&c\end{bmatrix}}^{-1}{\begin{bmatrix}t\\u\\v\\w\end{bmatrix}}}

Beispielsweise soll folgender Bruch zerlegt werden:

x 3 + 2 x 2 + 4 x + 6 ( 3 x 2 + 5 x + 7 ) ( 11 x 2 + 13 x + 17 ) {\displaystyle {\frac {x^{3}+2x^{2}+4x+6}{(3x^{2}+5x+7)(11x^{2}+13x+17)}}}

Hierfür muss nach diesem Verfahren folgende Rechnung durchgeführt werden:

[ y 1 y 2 y 3 y 4 ] = [ 11 0 3 0 13 11 5 3 17 13 7 5 0 17 0 7 ] 1 [ 1 2 4 6 ] = 1 290 [ 82 57 204 387 ] {\displaystyle {\begin{bmatrix}y_{1}\\y_{2}\\y_{3}\\y_{4}\end{bmatrix}}={\begin{bmatrix}11&0&3&0\\13&11&5&3\\17&13&7&5\\0&17&0&7\end{bmatrix}}^{-1}{\begin{bmatrix}1\\2\\4\\6\end{bmatrix}}={\frac {1}{290}}{\begin{bmatrix}82\\-57\\-204\\387\end{bmatrix}}}

Daraus folgt:

x 3 + 2 x 2 + 4 x + 6 ( 3 x 2 + 5 x + 7 ) ( 11 x 2 + 13 x + 17 ) = 82 x 57 290 ( 3 x 2 + 5 x + 7 ) + 204 x + 387 290 ( 11 x 2 + 13 x + 17 ) {\displaystyle {\frac {x^{3}+2x^{2}+4x+6}{(3x^{2}+5x+7)(11x^{2}+13x+17)}}={\frac {82x-57}{290(3x^{2}+5x+7)}}+{\frac {-204x+387}{290(11x^{2}+13x+17)}}}

Der Hauptsatz der Partialbruchzerlegung

Reellwertige Funktionen

Jede rationale Funktion R : D R R {\displaystyle R\colon D\subset \mathbb {R} \to \mathbb {R} } mit den m {\displaystyle m} verschiedenen reellen Polstellen x i {\displaystyle x_{i}} der Ordnung r i {\displaystyle r_{i}} und den n {\displaystyle n} bis auf Konjugation verschiedenen komplexen Polstellen z i {\displaystyle z_{i}} der Ordnung s i {\displaystyle s_{i}} hat eine eindeutig bestimmte Darstellung

R ( x ) = P ( x ) + i = 1 m j = 1 r i a i j ( x x i ) j + i = 1 n j = 1 s i b i j x + c i j ( x z i ) j ( x z i ¯ ) j {\displaystyle R(x)=P(x)+\sum _{i=1}^{m}\sum _{j=1}^{r_{i}}{\frac {a_{ij}}{(x-x_{i})^{j}}}+\sum _{i=1}^{n}\sum _{j=1}^{s_{i}}{\frac {b_{ij}x+c_{ij}}{(x-z_{i})^{j}(x-{\overline {z_{i}}})^{j}}}}

mit einer Polynomfunktion P {\displaystyle P} und reellen Konstanten a i j {\displaystyle a_{ij}} , b i j {\displaystyle b_{ij}} und c i j {\displaystyle c_{ij}} . Diese wird die Partialbruchzerlegung (abgekürzt PBZ) von R {\displaystyle R} genannt.

Die Brüche a i j ( x x i ) j {\displaystyle {\tfrac {a_{ij}}{(x-x_{i})^{j}}}} heißen Partial- oder Teilbrüche 1. Art, die Brüche b i j x + c i j ( x z i ) j ( x z ¯ i ) j {\displaystyle {\tfrac {b_{ij}x+c_{ij}}{(x-z_{i})^{j}(x-{\overline {z}}_{i})^{j}}}} Partial- oder Teilbrüche 2. Art.

Komplexwertige Funktionen

Jede rationale Funktion R : D C C {\displaystyle R\colon D\subset \mathbb {C} \to \mathbb {C} } mit den n {\displaystyle n} verschiedenen Polstellen z i {\displaystyle z_{i}} der Ordnung r i {\displaystyle r_{i}} hat eine eindeutig bestimmte Darstellung

R ( z ) = P ( z ) + i = 1 n j = 1 r i a i j ( z z i ) j {\displaystyle R(z)=P(z)+\sum _{i=1}^{n}\sum _{j=1}^{r_{i}}{\frac {a_{ij}}{(z-z_{i})^{j}}}}

mit einer Polynomfunktion P {\displaystyle P} und komplexen Konstanten a i j {\displaystyle a_{ij}} .

Dieser Satz lässt sich für Polynome über jedem anderen algebraisch abgeschlossenen Schiefkörper verallgemeinern.

Anwendungen

Die Partialbruchzerlegung wird unter anderem zum Integrieren rationaler Funktionen benutzt. Da die Integrale sämtlicher Partialbrüche bekannt sind, ist die Integration immer möglich, wenn sich die Polstellen der betrachteten Funktion angeben lassen.[2]

Des Weiteren wird die Partialbruchzerlegung bei der Laplace- und der z-Transformation verwendet. Die Transformierten der einzelnen Partialbrüche können in Tabellen nachgeschlagen werden. Somit erspart man sich eine analytische Berechnung, wenn der zu transformierende Term in entsprechende Summanden zerlegt werden kann.

Integration der Partialbrüche

Beim Auffinden der Stammfunktionen von Partialbrüchen lassen sich sechs Fälle unterscheiden, je nachdem, ob der Zählergrad 0 oder 1 ist, ob die Polstellen, also die Nullstellen des Nenners, reell oder nicht reell sind und ob sie einfach oder mehrfach sind.

Partialbrüche mit reellen Polstellen

Bei Partialbrüchen mit reellen Polstellen gibt es zwei Fälle, da der Zähler nur den Grad 0 haben kann.

Damit ergibt sich bei reellen und einfachen Polstellen

( 1 ) A x a d x = A ln | x a | + c {\displaystyle (1)\qquad \int {\frac {A}{x-a}}dx=A\cdot \ln |x-a|+c}

und bei reellen und mehrfachen Polstellen ( n 1 {\displaystyle n\neq 1} )

( 2 ) A ( x a ) n d x = A ( x a ) ( n 1 ) n 1 + c {\displaystyle (2)\qquad \int {\frac {A}{{(x-a)}^{n}}}dx={\frac {-A\cdot (x-a)^{-(n-1)}}{n-1}}+c} .

Partialbrüche mit komplexen Polstellen

Bei Partialbrüchen mit komplexen Polstellen gibt es vier Fälle, da der Zählergrad sowohl 0 als auch 1 sein kann.

Damit ergibt sich bei komplexen und einfachen Polstellen und Zählergrad 0

( 3 ) B x 2 + p x + q d x = 2 B 4 q p 2 arctan ( 2 x + p 4 q p 2 ) + c {\displaystyle (3)\qquad \int {\frac {B}{x^{2}+px+q}}dx={\frac {2B}{\sqrt {4q-p^{2}}}}\cdot \arctan \left({\frac {2x+p}{\sqrt {4q-p^{2}}}}\right)+c} .

Der Fall mit komplexen und einfachen Polstellen und Zählergrad 1 lässt sich auf (3) zurückführen:

( 4 ) B x + C x 2 + p x + q d x = B 2 ln ( x 2 + p x + q ) + ( C p B 2 ) 1 x 2 + p x + q d x + c {\displaystyle (4)\qquad \int {\frac {Bx+C}{x^{2}+px+q}}dx={\frac {B}{2}}\cdot \ln(x^{2}+px+q)+\left(C-{\frac {pB}{2}}\right)\cdot \int {\frac {1}{x^{2}+px+q}}dx+c}

Für die beiden Fälle mit mehrfachen Polstellen lassen sich nicht direkt Stammfunktionen bestimmen, es lassen sich jedoch Rekursionsvorschriften finden. Damit ergibt sich für den Fall mit komplexen und mehrfachen Polstellen ( n 0 {\displaystyle n\neq 0} ) und Zählergrad 0

( 5 ) B ( x 2 + p x + q ) n + 1 d x = B ( 4 q p 2 ) n 2 x + p ( x 2 + p x + q ) n + 2 4 q p 2 2 n 1 n B ( x 2 + p x + q ) n d x + c {\displaystyle (5)\qquad \int {\frac {B}{(x^{2}+px+q)^{n+1}}}dx={\frac {B}{(4q-p^{2})\cdot n}}\cdot {\frac {2x+p}{(x^{2}+px+q)^{n}}}+{\frac {2}{4q-p^{2}}}\cdot {\frac {2n-1}{n}}\cdot \int {\frac {B}{(x^{2}+px+q)^{n}}}dx+c} .

Der Fall mit komplexen und mehrfachen Polstellen und Zählergrad 1 lässt sich auf (5) zurückführen ( n 1 {\displaystyle n\neq 1} )

( 6 ) B x + C ( x 2 + p x + q ) n d x = B 2 ( n 1 ) 1 ( x 2 + p x + q ) n 1 + ( C p B 2 ) 1 ( x 2 + p x + q ) n d x + c {\displaystyle (6)\qquad \int {\frac {Bx+C}{(x^{2}+px+q)^{n}}}dx=-{\frac {B}{2(n-1)}}\cdot {\frac {1}{{(x^{2}+px+q)}^{n-1}}}+\left(C-{\frac {pB}{2}}\right)\cdot \int {\frac {1}{(x^{2}+px+q)^{n}}}dx+c} .

Laurent-Reihen-Entwicklung

Ist für jede Polstelle eine Laurent-Reihen-Entwicklung der Funktion bekannt, so erhält man die Partialbruchzerlegung sehr einfach als Summe der Hauptteile dieser Laurent-Reihen. Dieser Weg steht im Zusammenhang mit dem Residuenkalkül.

Verallgemeinerung auf rationale Funktionenkörper

Die Partialbruchzerlegung lässt sich für einen Körper K {\displaystyle K} auf den rationalen Funktionenkörper K ( X ) {\displaystyle K(X)} verallgemeinern. Bezeichnet man die normierten irreduziblen Polynome im Polynomring K [ X ] {\displaystyle K[X]} mit P {\displaystyle P} , so sind die rationalen Funktionen der Form X i p j {\displaystyle {\tfrac {X^{i}}{p^{j}}}} mit p P , j N + , 0 i < deg ( p ) {\displaystyle p\in P,j\in \mathbb {N} _{+},0\leq i<\deg(p)} linear unabhängig und bilden mit den Monomen X i , i N {\displaystyle X^{i},i\in \mathbb {N} } eine K {\displaystyle K} -Basis des K {\displaystyle K} -Vektorraums K ( X ) {\displaystyle K(X)} .[3]

Literatur

  • Schülerduden Mathematik II. Bibliographisches Institut & F.A. Brockhaus, 2004, ISBN 3-411-04275-3, S. 316–317.
  • Charles D. Miller, Margaret L. Lial, David I. Schneider: Fundamentals of College Algebra. 3. Auflage. Scott & Foresman / Little & Brown Higher Education, 1990, ISBN 0-673-38638-4, S. 364–370.
  • L.D. Kudryavtsev: Undetermined coefficients, method of. In: Michiel Hazewinkel (Hrsg.): Encyclopedia of Mathematics. Springer-Verlag und EMS Press, Berlin 2002, ISBN 1-55608-010-7 (englisch, encyclopediaofmath.org). 
  • Eric W. Weisstein: Partial Fraction Decomposition. In: MathWorld (englisch).
  • Günter Scheja, Uwe Storch: Lehrbuch der Algebra. Teubner, Stuttgart 1988, ISBN 3-519-02212-5.

Einzelnachweise

  1. Heinz-Wilhelm Alten: 4000 Jahre Algebra. Geschichte, Kulturen, Menschen. Springer, Berlin u. a. 2003, ISBN 3-540-43554-9, S. 285–286. 
  2. Christoph Bock: Elemente der Analysis. (PDF; 2,2 MB) Abschnitt 8.35.
  3. Günter Scheja, Uwe Storch: Lehrbuch der Algebra. Teubner, Stuttgart 1988, ISBN 3-519-02212-5, S. 148.