Kurs:Mathematik für Anwender (Osnabrück 2011-2012)/Teil I/Vorlesung 26

Stammfunktionen zu rationalen Funktionen

Wir möchten zeigen, wie man zu einer rationalen Funktion ${}P/Q$ (gegeben durch zwei Polynome $P,Q$, $Q\neq 0$) eine Stammfunktion finden kann. Im Allgemeinen gehen bei der Bestimmung einer solchen Stammfunktion verschiedene Techniken ein, die wir im Laufe der Vorlesung kennengelernt haben: die Division mit Rest für Polynome, das Lösen inhomogener linearer Gleichungssysteme und Integrationsregeln.

Wenn das Nennerpolynom ${}Q=1$ ist, so handelt es sich einfach um ein Polynom ${}P$, das problemlos zu integrieren ist. Für die Funktion ${}1/x$ ist der natürliche Logarithmus eine Stammfunktion. Damit ist auch eine Funktion vom Typ

${\frac {1}{ax+b}}$

(mit $a\neq 0$) integrierbar, eine Stammfunktion ist ${}{\frac {1}{a}}\ln(ax+b)$. Damit kann man überhaupt beliebige rationale Funktionen der Form

${\frac {P}{ax+b}}$

integrieren. Die Division mit Rest führt zu einer Darstellung

$P=H(ax+b)+c,$

mit einem weiteren Polynom ${}H$, und wobei das Restpolynom ${}c$ konstant ist, da sein Grad kleiner als der Grad des linearen Polynoms ist, durch das die Division durchgeführt wird. Aus dieser Gleichung erhält man die Darstellung

${\frac {P}{ax+b}}=H+{\frac {c}{ax+b}},$

wobei wir für die beiden Summanden Stammfunktionen angeben können. Die Division mit Rest wird auch im allgemeinen Fall entscheidend sein. Davor betrachten wir aber noch den Fall eines quadratischen Nennerpolynoms mit Zähler ${}1$, also

${\frac {1}{ax^{2}+bx+c}}$

mit ${}a,b,c\in \mathbb {R} ,\,a\neq 0$. Durch Multiplikation mit ${}a$ kann man den Koeffizienten vor ${}x^{2}$ zu ${}1$ normieren. Durch quadratisches Ergänzen kann man

$x^{2}+bx+c=(x+d)^{2}+e$

schreiben. Mit der neuen Variablen ${}u=x+d$ (bzw. mit der Substitution $u=x+d$) schreibt sich dies als ${}u^{2}+e$. Mit einer weiteren Substitution unter Verwendung der Quadratwurzel von ${}e$ bzw. von ${}-e$ gelangt man zu

${\frac {1}{1+v^{2}}}{\text{ oder }}{\frac {1}{1-v^{2}}}.$

Im ersten Fall gilt

$\int {\frac {1}{1+v^{2}}}dv=\arctan v$

und im zweiten Fall gilt

$\int {\frac {1}{1-v^{2}}}dv={\frac {1}{2}}\ln {\frac {1+v}{1-v}},$

wie früher gezeigt wurde. Für die inversen Funktionen zu Potenzen von quadratischen nullstellenfreien Polynomen werden die Stammfunktionen durch folgende Rekursionsformel bestimmt.

Lemma

Es sei ${}x^{2}+bx+c$ (mit ${}b,c\in \mathbb {R}$) ein quadratisches Polynom ohne reelle Nullstelle (d.h., dass ${}\triangle :={\frac {b^{2}-4c}{4}}<0$ ist).

Dann ist

${\frac {1}{\sqrt {-\triangle }}}\arctan {\frac {1}{\sqrt {-\triangle }}}{\left(x+{\frac {b}{2}}\right)}$

eine Stammfunktion von

${\frac {1}{x^{2}+bx+c}}$

und für ${}n\geq 1$ gilt die Rekursionsformel

$\int _{r}^{s}{\frac {1}{{\left(x^{2}+bx+c\right)}^{n+1}}}dx={\frac {1}{n{\left(4c-b^{2}\right)}}}{\left({\frac {2z+b}{{\left(z^{2}+bz+c\right)}^{n}}}{|}_{r}^{s}+(4n-2)\int _{r}^{s}{\frac {1}{{\left(x^{2}+bx+c\right)}^{n}}}dx\right)}\,.$

Beweis

Ableiten ergibt

${}{\begin{aligned}{\left({\frac {1}{\sqrt {-\triangle }}}\arctan {\left({\frac {1}{\sqrt {-\triangle }}}{\left(x+{\frac {b}{2}}\right)}\right)}\right)}'&={\frac {1}{\sqrt {-\triangle }}}\cdot {\frac {1}{\sqrt {-\triangle }}}\cdot {\frac {1}{1+{\frac {1}{-\triangle }}{\left(x+{\frac {b}{2}}\right)}^{2}}}\\&={\frac {1}{-\triangle +{\left(x+{\frac {b}{2}}\right)}^{2}}}\\&={\frac {1}{c-{\frac {b^{2}}{4}}+x^{2}+bx+{\frac {b^{2}}{4}}}}\\&={\frac {1}{x^{2}+bx+c}}.\end{aligned}}$

Zum Beweis der Rekursionsformel setzen wir ${}q(x):=x^{2}+bx+c$ und leiten ab.

${}{\begin{aligned}{\left((2x+b)q^{-n}\right)}'&=2q^{-n}-n(2x+b)q^{-n-1}(2x+b)\\&=2q^{-n}-nq^{-n-1}(2x+b)^{2}\\&=2q^{-n}-nq^{-n-1}{\left(4x^{2}+4xb+b^{2}\right)}\\&=2q^{-n}-nq^{-n-1}{\left(4q-4c+b^{2}\right)}\\&=2q^{-n}-4nq^{-n}+n{\left(4c-b^{2}\right)}q^{-n-1}\\&=(2-4n)q^{-n}+n{\left(4c-b^{2}\right)}q^{-n-1}.\end{aligned}}$

Division durch ${}n{\left(4c-b^{2}\right)}$ und Umstellen ergibt

$q^{-n-1}={\frac {1}{n{\left(4c-b^{2}\right)}}}{\left((2x+b)q^{-n}\right)}'+{\frac {4n-2}{n{\left(4c-b^{2}\right)}}}q^{-n}\,.$

Dies ist die Behauptung.

$\Box$

Bemerkung  

Mit Lemma 26.1 kann man auch rationale Funktionen der Form

${\frac {rx+s}{(x^{2}+bx+c)^{n}}}$

(mit ${}r,s\in \mathbb {R} ,\,r\neq 0$,) integrieren, wo also das Zählerpolynom linear und das Nennerpolynom eine Potenz eines quadratischen Polynoms ist. Bei ${}n=1$ ist

${\left({\frac {r}{2}}\ln \left(x^{2}+bx+c\right)\right)}'={\frac {r}{2}}\cdot {\frac {2x+b}{x^{2}+bx+c}}={\frac {rx+{\frac {rb}{2}}}{x^{2}+bx+c}}\,.$

D.h., dass die Differenz zwischen dieser Ableitung und der zu integrierenden Funktion vom Typ

${\frac {u}{x^{2}+bx+c}}$

ist, was wir aufgrund von Lemma 26.1 integrieren können. Bei ${}n\geq 2$ ist

${}{\begin{aligned}{\left({\frac {-r}{2(n-1)}}\cdot {\frac {1}{{\left(x^{2}+bx+c\right)}^{n-1}}}\right)}'&={\frac {-r}{2(n-1)}}\cdot (-n+1)\cdot (2x+b)\cdot {\frac {1}{{\left(x^{2}+bx+c\right)}^{n}}}\\&={\frac {rx+{\frac {rb}{2}}}{{\left(x^{2}+bx+c\right)}^{n}}}\,\end{aligned}}$

und wieder ist das Integral auf eine schon behandelte Situation zurückgeführt.

Wir möchten für beliebige rationale Funktionen ${}f={\frac {P}{Q}}$ mit ${}P,Q\in \mathbb {R} [X]$ Stammfunktionen bestimmen. Dies geht grundsätzlich immer, vorausgesetzt, dass man eine Faktorzerlegung des Nennerpolynoms besitzt. Aufgrund der reellen Version des Fundamentalsatzes der Algebra gibt es eine Faktorzerlegung

$Q=(X-b_{1})\cdots (X-b_{r})\cdot q_{1}\cdots q_{s},$

wobei die ${}q_{j}$ quadratische Polynome ohne reelle Nullstellen sind. Das Bestimmen der Stammfunktionen zu rationalen Funktionen beruht auf der Partialbruchzerlegung von rationalen Funktionen, die wir zuerst besprechen.

Partialbruchzerlegung

Die Partialbruchzerlegung liefert eine wichtige Darstellungsform für eine rationale Funktion ${}P/Q$, bei der die Nenner besonders einfach werden. Wir beginnen mit dem Fall ${}K={\mathbb {C} }$, wo wir den Fundamentalsatz der Algebra zur Verfügung haben.

Satz

Es seien ${}P,Q\in {\mathbb {C} }[X]$, ${}Q\neq 0$, Polynome und es sei

${}Q=(X-a_{1})^{r_{1}}\cdots (X-a_{s})^{r_{s}}\,$

mit verschiedenen ${}a_{i}\in {\mathbb {C} }$.

Dann gibt es ein eindeutig bestimmtes Polynom ${}H\in {\mathbb {C} }[X]$ und eindeutig bestimmte Koeffizienten ${}c_{ij}\in {\mathbb {C} }$, ${}1\leq i\leq s$, ${}1\leq j\leq r_{i}$, mit

${}{\frac {P}{Q}}=H+{\frac {c_{11}}{X-a_{1}}}+{\frac {c_{12}}{(X-a_{1})^{2}}}+\cdots +{\frac {c_{1r_{1}}}{(X-a_{1})^{r_{1}}}}+\cdots +{\frac {c_{s1}}{X-a_{s}}}+{\frac {c_{s2}}{(X-a_{s})^{2}}}+\cdots +{\frac {c_{sr_{s}}}{(X-a_{s})^{r_{s}}}}\,.$

Beweis

Wir verzichten auf den Beweis, obwohl uns die dazu benötigten Methoden zur Verfügung stehen.

$\Box$

Wir wenden uns nun der reellen Situation zu.

Korollar

Es seien ${}P,Q\in \mathbb {R} [X]$, ${}Q\neq 0$, Polynome und es sei

$Q=(X-a_{1})^{r_{1}}\cdots (X-a_{s})^{r_{s}}Q_{1}^{t_{1}}\cdots Q_{u}^{t_{u}}$

mit verschiedenen ${}a_{i}\in \mathbb {R}$ und verschiedenen quadratischen Polynomen ${}Q_{k}$ ohne reelle Nullstellen.

Dann gibt es ein eindeutig bestimmtes Polynom ${}H\in \mathbb {R} [X]$ und eindeutig bestimmte Koeffizienten ${}c_{ij}\in \mathbb {R}$, ${}1\leq i\leq s$, ${}1\leq j\leq r_{i}$, und eindeutig bestimmte lineare Polynome ${}L_{k\ell }=d_{k\ell }X+e_{k\ell }$, ${}1\leq k\leq u$, ${}1\leq \ell \leq t_{k}$, mit

${}{\begin{aligned}{\frac {P}{Q}}&=H+{\frac {c_{11}}{X-a_{1}}}+{\frac {c_{12}}{(X-a_{1})^{2}}}+\cdots +{\frac {c_{1r_{1}}}{(X-a_{1})^{r_{1}}}}+\cdots +{\frac {c_{s1}}{X-a_{s}}}+{\frac {c_{s2}}{(X-a_{s})^{2}}}+\cdots +{\frac {c_{sr_{s}}}{(X-a_{s})^{r_{s}}}}\\&\,\,\,+{\frac {L_{11}}{Q_{1}}}+{\frac {L_{12}}{Q_{1}^{2}}}+\cdots +{\frac {L_{1t_{1}}}{Q_{1}^{t_{1}}}}+\cdots +{\frac {L_{u1}}{Q_{u}}}+{\frac {L_{u2}}{Q_{u}^{2}}}+\cdots +{\frac {L_{ut_{u}}}{Q_{u}^{t_{u}}}}.\end{aligned}}$

Beweis

Auch hier verzichten wir auf den Beweis, der auf der komplexen Partialbruchzerlegung beruht.

$\Box$

Neben dem Umweg über die komplexe Partialbruchzerlegung gibt es weitere Methoden, in Beispielen die reelle Partialbruchzerlegung zu bestimmen. Grundsätzlich bedeutet das Bestimmen der (reellen oder komplexen) Koeffizienten in der Partialbruchzerlegung, ein (inhomogenes) lineares Gleichungssystem zu lösen, wobei man sowohl durch Koeffizientenvergleich als auch durch das Einsetzen von bestimmten Zahlen zu hinreichend vielen linearen Gleichungen kommt.

Beispiel

Wir betrachten die rationale Funktion

${}{\frac {1}{X^{3}-1}}={\frac {1}{(X-1){\left(X^{2}+X+1\right)}}}\,.$

wobei der Faktor rechts reell nicht weiter zerlegbar ist. Daher muss es eine eindeutige Darstellung

${}{\frac {1}{X^{3}-1}}={\frac {a}{X-1}}+{\frac {bX+c}{X^{2}+X+1}}\,$

geben. Multiplikation mit dem Nennerpolynom führt auf

${}{\begin{aligned}1&=a{\left(X^{2}+X+1\right)}+(bX+c)(X-1)\\&=(a+b)X^{2}+(a+c-b)X+a-c.\end{aligned}}$

Koeffizientenvergleich führt auf das inhomogene lineare Gleichungssystem

$a+b=0{\text{ und }}a+c-b=0{\text{ und }}a-c=1$

mit den eindeutigen Lösungen

$a={\frac {1}{3}},\,b=-{\frac {1}{3}},\,c=-{\frac {2}{3}}.$

Die Partialbruchzerlegung ist also

${}{\frac {1}{X^{3}-1}}={\frac {\frac {1}{3}}{X-1}}+{\frac {-{\frac {1}{3}}X-{\frac {2}{3}}}{X^{2}+X+1}}={\frac {1}{3}}\cdot {\frac {1}{X-1}}-{\frac {1}{3}}\cdot {\frac {X+2}{X^{2}+X+1}}\,.$

Beispiel

Wir betrachten die rationale Funktion

${}{\frac {X^{3}-X+5}{X^{4}+X^{2}}}={\frac {X^{3}-X+5}{X^{2}{\left(X^{2}+1\right)}}}\,,$

wo die Faktorzerlegung des Nennerpolynoms sofort ersichtlich ist. Der Ansatz

${}{\frac {X^{3}-X+5}{X^{2}{\left(X^{2}+1\right)}}}={\frac {a}{X}}+{\frac {b}{X^{2}}}+{\frac {cX+d}{X^{2}+1}}\,$

führt durch Multiplikation mit dem Nennerpolynom auf

${}{\begin{aligned}X^{3}-X+5&=aX{\left(X^{2}+1\right)}+b{\left(X^{2}+1\right)}+(cX+d)X^{2}\\&=aX^{3}+aX+bX^{2}+b+cX^{3}+dX^{2}\\&=(a+c)X^{3}+(b+d)X^{2}+aX+b.\end{aligned}}$

Koeffizientenvergleich führt auf das inhomogene lineare Gleichungssystem

$a+c=1{\text{ und }}b+d=0{\text{ und }}a=-1{\text{ und }}b=5$

mit der Lösung

$b=5,\,a=-1,\,d=-5,\,c=2.$

Insgesamt ist die Partialbruchzerlegung also gleich

${}{\frac {X^{3}-X+5}{X^{2}(X^{2}+1)}}=-{\frac {1}{X}}+{\frac {5}{X^{2}}}+{\frac {2X-5}{X^{2}+1}}\,.$

Integration rationaler Funktionen

Verfahren  

Es sei eine rationale Funktion

${}f={\frac {P}{Q}}\,$

gegeben, für die eine Stammfunktion gefunden werden soll. Dabei seien ${}P$ und ${}Q$ reelle Polynome. Man geht folgendermaßen vor.

1. Bestimme die reelle Faktorzerlegung des Nennerpolynoms ${}Q$.
2. Finde die Partialbruchzerlegung

   ${}{\frac {P}{Q}}=H+\sum _{i=1}^{r}{\left(\sum _{j=1}^{s_{i}}{\frac {c_{ij}}{(X-a_{i})^{j}}}\right)}+\sum _{k=1}^{u}{\left(\sum _{\ell =1}^{t_{k}}{\frac {d_{k{\ell }}X+e_{k{\ell }}}{Q_{k}^{\ell }}}\right)}\,.$

3. Bestimme für ${}H$, für jedes

   ${\frac {c_{ij}}{(X-a_{i})^{j}}}$

   und für jedes

   ${\frac {d_{k{\ell }}X+e_{k{\ell }}}{Q_{k}^{\ell }}}$

   eine Stammfunktion.

Beispiel

Wir möchten eine Stammfunktion zu

${}f(x)={\frac {1}{x^{3}-1}}\,$

bestimmen. Nach Beispiel 26.5 ist die reelle Partialbruchzerlegung gleich

${}{\begin{aligned}{\frac {1}{x^{3}-1}}&={\frac {1}{3}}\cdot {\frac {1}{x-1}}-{\frac {1}{3}}\cdot {\frac {x+2}{x^{2}+x+1}}\\&={\frac {1}{3}}\cdot {\frac {1}{x-1}}-{\frac {1}{6}}\cdot {\frac {2x+4}{x^{2}+x+1}}\\&={\frac {1}{3}}\cdot {\frac {1}{x-1}}-{\frac {1}{6}}\cdot {\frac {2x+1}{x^{2}+x+1}}-{\frac {1}{6}}\cdot {\frac {3}{x^{2}+x+1}}\\&={\frac {1}{3}}\cdot {\frac {1}{x-1}}-{\frac {1}{6}}\cdot {\frac {2x+1}{x^{2}+x+1}}-{\frac {1}{2}}\cdot {\frac {1}{x^{2}+x+1}}.\end{aligned}}$

Als Stammfunktion ergibt sich daher

${\frac {1}{3}}\ln \vert {x-1}\vert -{\frac {1}{6}}\cdot \ln \left(x^{2}+x+1\right)-{\frac {1}{\sqrt {3}}}\cdot \arctan {\frac {2}{\sqrt {3}}}{\left(x+{\frac {1}{2}}\right)},$

wobei wir für den rechten Summanden Lemma 26.1 verwendet haben.

Beispiel

Wir möchten eine Stammfunktion zu

${}f(x)={\frac {x^{3}-x+5}{x^{4}+x^{2}}}\,$

bestimmen. Nach Beispiel 26.6 ist die reelle Partialbruchzerlegung gleich

${}{\frac {x^{3}-x+5}{x^{2}(x^{2}+1)}}=-{\frac {1}{x}}+{\frac {5}{x^{2}}}+{\frac {2x-5}{x^{2}+1}}\,.$

Als Stammfunktion ergibt sich daher

$-\ln \vert {x}\vert -5x^{-1}+\ln \left(x^{2}+1\right)-5\arctan x.$

Stammfunktionen zu rationalen Funktionen in der Exponentialfunktion

Nachdem wir nun rationale Funktionen integrieren können, können wir auch für eine ganze Reihe von Funktionen eine Stammfunktion finden, die wir durch gewisse Standardsubstitution auf eine rationale Funktion zurückführen können. Wir führen dies exemplarisch für Funktionen durch, die sich als rationale Funktionen in der Exponentialfunktion schreiben lassen.

Lemma

Es sei ${}f$ eine rationale Funktion in der Exponentialfunktion, d.h. es gebe Polynome ${}P,Q\in \mathbb {R} [X]$, ${}Q\neq 0$, derart, dass

${}f(t)={\frac {P{\left(e^{t}\right)}}{Q{\left(e^{t}\right)}}}\,$

gilt.

Dann kann man durch die Substitution ${}t=\ln s$ das Integral ${}\int _{}^{}f(t)\,dt$ auf das Integral einer rationalen Funktion zurückführen.

Beweis

Bei der Substitution ${}t=\ln s$ ist

${}dt={\frac {1}{s}}ds\,,$

und für die Polynome ${}P{\left(e^{t}\right)}$ und ${}Q{\left(e^{t}\right)}$ ergeben sich

$P{\left(e^{t}\right)}=P{\left(e^{\ln s}\right)}=P(s)\,\,{\text{ und }}\,\,Q{\left(e^{t}\right)}=Q{\left(e^{\ln s}\right)}=Q(s).$

Insgesamt ergibt sich also die rationale Funktion ${}{\frac {P(s)}{sQ(s)}}$. In deren Stammfunktion muss man dann ${}s=e^{t}$ einsetzen.

$\Box$

Beispiel

Wir wollen eine Stammfunktion für die Funktion

${}f(t)={\frac {1}{e^{t}+e^{3t}}}\,$

finden. Das in Lemma 26.10 beschriebene Verfahren führt auf die rationale Funktion

${}{\frac {1}{{\left(s+s^{3}\right)}s}}={\frac {1}{s^{2}{\left(s^{2}+1\right)}}}={\frac {1}{s^{2}}}-{\frac {1}{s^{2}+1}}\,,$

sodass die Partialbruchzerlegung direkt vorliegt. Die Stammfunktion von dieser rationalen Funktion ist

$-{\frac {1}{s}}-\arctan s.$

Die Stammfunktion von ${}f$ ist daher

$-{\frac {1}{e^{t}}}-\arctan e^{t}.$

Fußnoten

<< | Kurs:Mathematik für Anwender (Osnabrück 2011-2012)/Teil I | >> PDF-Version dieser Vorlesung Arbeitsblatt zur Vorlesung (PDF) (PDF englisch)