Polynomdivision

Im Folgenden seien ${\displaystyle n}$  und ${\displaystyle m}$  natürliche Zahlen einschließlich Null (${\displaystyle n,m\in \mathbb {N} _{0}}$ ) und der Einfachheit halber die Größen ${\displaystyle a_{i}\ (0\leq i\leq n)}$  und ${\displaystyle b_{j}\ (0\leq j\leq m)}$  stets ganze Zahlen, also Elemente von ${\displaystyle \mathbb {Z} }$ . Hat man nun zwei Polynome, etwa

${\displaystyle p(x)=a_{n}x^{n}+a_{n-1}x^{n-1}+\cdots +a_{2}x^{2}+a_{1}x+a_{0}\quad {\text{mit}}\ a_{n}\neq 0}$ 

und

${\displaystyle q(x)=b_{m}x^{m}+b_{m-1}x^{m-1}+\cdots +b_{2}x^{2}+b_{1}x+b_{0}\quad {\text{mit}}\ b_{m}\neq 0,}$ 

so kann man sie unter gewissen formalen Voraussetzungen ähnlich wie ganze Zahlen durcheinander dividieren, also die Rechenaufgabe

${\displaystyle p(x):q(x)=\;?}$ 

lösen. Im Ergebnis finden sich dann zwei Polynome: Ein Polynom ${\displaystyle s(x)}$ , das dem Ganzzahlquotienten in der Zahlendivision mit Rest entspricht, und ein Polynom ${\displaystyle r(x)}$ , das sich nicht mehr weiter durch ${\displaystyle q(x)}$  teilen lässt und das dem Rest in der Zahlendivision entspricht:

${\displaystyle p(x)=q(x)s(x)+r(x)}$ 

oder in Analogie zur Schulschreibweise

${\displaystyle p(x):q(x)=s(x){\text{ Rest }}r(x).}$ 

Das Verfahren zum Auffinden dieser Lösung, bestehend aus ${\displaystyle s(x)}$  und ${\displaystyle r(x)}$ , ist die Polynomdivision.

Dass sich hiernach ${\displaystyle r(x)}$  nicht weiter durch ${\displaystyle q(x)}$  teilen lässt, ist gleichbedeutend damit, dass der Polynomgrad von ${\displaystyle r(x)}$  kleiner ist als der von ${\displaystyle q(x)}$ , weshalb dies in der formalen Definition der Rechenvorschrift (Algorithmus) auch als Abbruchbedingung gefordert wird. In der Zahlendivision mit Rest wird stattdessen gefordert, dass der Rest kleiner als der Divisor ist. Beide Nebenbedingungen sorgen im jeweiligen Verfahren dafür, dass der Rest eindeutig bestimmt ist.

Bei der formalen Definition des Verfahrens werden einige zusätzliche Bedingungen beachtet. Das kommt daher, dass man Polynome im Allgemeinen viel weitläufiger definieren kann, als es hier zur einfacheren Erklärung geschehen ist oder man es zum Beispiel aus der Schule kennt. Die Koeffizienten eines Polynoms etwa können dann aus beliebigen Ringen stammen. Dann dürfen aber wiederum die Koeffizienten der beiden Polynome nicht aus verschiedenen Ringen stammen. Daher definiert man, dass die Polynome in einem gemeinsamen Polynomring liegen müssen. Auch reicht es nicht mehr zu fordern, dass der „höchste“ Koeffizient (Leitkoeffizient) ${\displaystyle b_{m}}$  von ${\displaystyle q(x)}$  nur ungleich Null sein müsse. Vielmehr muss man fordern, dass er zudem eine Einheit des Ringes sein muss. Oder es wird das unten beschriebene Verfahren der Pseudo-Division angewendet.

Bei der Polynomdivision sind zwei Polynome ${\displaystyle p(x)}$  und ${\displaystyle q(x)}$  eines Polynomringes ${\displaystyle R[x]}$  gegeben, wobei ${\displaystyle R}$  ein kommutativer Ring mit ${\displaystyle 1\neq 0}$  und der Leitkoeffizient von ${\displaystyle q(x)}$  eine Einheit in ${\displaystyle R}$  ist, und es wird die Gleichung

${\displaystyle p(x)=s(x)q(x)+r(x),}$ 

nach den gesuchten Polynomen ${\displaystyle s(x)}$  und ${\displaystyle r(x)}$  gelöst, und zwar so, dass der Grad von ${\displaystyle r(x)}$  kleiner als der von ${\displaystyle q(x)}$  ist.

• Wegen ${\displaystyle \operatorname {grad} \,r(x)<\operatorname {grad} \,q(x)}$  sind die Polynome ${\displaystyle s(x)}$  und ${\displaystyle r(x)}$  in ${\displaystyle R[x]}$  eindeutig bestimmt.
• Die Polynomdivision ist im Allgemeinen keine innere Verknüpfung auf ${\displaystyle R[x]}$ , da sich als Ergebnis der Division zweier Polynome im Allgemeinen nicht ein einzelnes, sondern zwei Polynome in ${\displaystyle R[x]}$  ergeben und sich somit keine Zuordnung der Form ${\displaystyle R[x]\times R[x]\rightarrow R[x]}$  machen lässt. Ist das jedoch im Einzelfall möglich, so wird ${\displaystyle R[x]}$  zu einem Körper mit der Polynomdivision als Umkehrung der Polynommultiplikation.
• Ist die Polynomdivision für jedes beliebige Paar von zwei Polynomen aus ${\displaystyle R[x]}$  möglich, so wird ${\displaystyle R[x]}$  zu einem euklidischen Ring bzgl. der Polynomgrad-Funktion. Das ist genau dann der Fall, wenn ${\displaystyle R}$  ein Körper ist.

• Eine Anwendung ist das Lösen von Gleichungen höheren Grades. Wenn eine Lösung (Nullstelle) ${\displaystyle x_{n}}$  bekannt ist, findet die Polynomdivision Anwendung, um den Grad der Gleichung um Eins zu senken. Diese Vorgehensweise wird „Abspalten einer Nullstelle“ genannt.
• Eine weitere Anwendung findet die Polynomdivision bei der Kurvendiskussion mit der Bestimmung der Näherungskurven einer rationalen Funktion.
• Bei der Partialbruchzerlegung rationaler Funktionen wird die Polynomdivision ebenfalls benötigt.
• Bei der Berechnung von Prüfsummen findet die Polynomdivision über dem Ring der ganzen Zahlen modulo 2 Anwendung, siehe CRC-Polynom.
• Nach erfolgter Polynomdivision kann man dasselbe Verfahren auf Divisor und Rest erneut anwenden und so einen weiteren Rest berechnen, und so weiter. Man erhält dann eine sogenannte Polynomrestfolge.

Das Verfahren funktioniert für Polynome mit ganzzahligen Koeffizienten genau so wie die schriftliche Division ganzer Zahlen mit Rest und kann mit dem gleichen Schema (Verfahren, Vorgehensweise) gelöst werden. Hier werden die einzelnen Schritte am Beispiel

${\displaystyle {\frac {p(x)}{q(x)}}={\frac {4\cdot x^{5}-x^{4}+2\cdot x^{3}+x^{2}-1}{x^{2}+1}}}$ 

erläutert:

• Wie bei der Division ganzer Zahlen wird zuerst der Summand höchsten Grades des Polynoms ${\displaystyle p}$  eliminiert. Dazu wird zunächst der Summand höchsten Grades von ${\displaystyle p}$  durch den Summanden höchsten Grades von ${\displaystyle q}$  dividiert. Das Ergebnis ist ${\displaystyle {\frac {4x^{5}}{x^{2}}}=4x^{3}}$ . Danach wird ${\displaystyle q}$  mit ${\displaystyle 4x^{3}}$  multipliziert und von ${\displaystyle p}$  subtrahiert.

${\displaystyle {\begin{matrix}(4x^{5}&-x^{4}&+2x^{3}&+x^{2}&-1)&:&(x^{2}&+1)&=&4x^{3}+{\frac {-x^{4}-2x^{3}+x^{2}-1}{x^{2}+1}}\\{\underline {-4x^{5}}}&&{\underline {-4x^{3}}}\\&-x^{4}&-2x^{3}\end{matrix}}}$ 

Es bleibt der Rest ${\displaystyle -x^{4}-2x^{3}+x^{2}-1}$ . Sein Grad ist nicht kleiner als der des Divisors.

• Im nächsten Schritt wird von diesem Rest der Summand höchsten Grades eliminiert, bis in mehreren solchen Schritten ein Rest entsteht (hier: ${\displaystyle 2x-3}$ ), der nicht mehr eliminiert werden kann, weil sein Grad kleiner als der Grad des Divisors ${\displaystyle q}$  ist.

${\displaystyle {\begin{array}{rrrrl}(-x^{4}&-2x^{3}&+x^{2}&&-1)&:&(x^{2}&&+1)=-x^{2}&-2x&+2&+{\frac {2x-3}{x^{2}+1}}\\{\underline {+x^{4}}}&&{\underline {+x^{2}}}\\&-2x^{3}&+2x^{2}\\&{\underline {+2x^{3}}}&&{\underline {+2x}}\\&&+2x^{2}&+2x\\&&{\underline {-2x^{2}}}&&{\underline {-2}}\\&&&+2x&-3\end{array}}}$ 

Weitere Beispiele

Das folgende Code-Fragment in BASIC zeigt den Kern der Berechnung:

For i = GradZ - GradN To 0 Step -1
    Quotient(i) = Zähler(i + GradN) / Nenner(GradN)
    For j = GradN To 0 Step -1
        Zähler(i + j) = Zähler(i + j) - Nenner(j) * Quotient(i)
    Next j
Next i
For j = GradN - 1 To 0 Step -1
    Rest(j) = Zähler(j)
Next j

Die Variable Zähler() ist ein Feld (Array), welches die Koeffizienten des Zählerpolynoms enthält, so dass Zähler(i) den Koeffizienten der Potenz ${\displaystyle x^{i}}$  enthält. Entsprechend ist Nenner() ein weiteres Feld, welches in gleicher Art die Koeffizienten des Nennerpolynoms enthält. Das Ergebnis sind zwei Polynome, welche in Quotient() und Rest() ausgegeben werden. Die Variablen GradN und GradZ enthalten den jeweiligen Polynomgrad von Zähler und Nenner.

In einem optimierten Programm könnte man die innere Schleife von 0 bis GradN-1 laufen lassen und die Ergebnisse in Zähler() zurückschreiben, so dass die Variablen Quotient() und Rest() entfallen würden. Der Einfachheit halber wurde hier darauf verzichtet.

Die oben beschriebene Methode zur Polynomdivision ist nur dann anwendbar, wenn der Leitkoeffizient des Divisors ${\displaystyle q(x)}$  eine Einheit im Grundring ist. Das ist immer der Fall, wenn der Grundring ein Körper ist. Über allgemeinen Grundringen muss das jedoch nicht immer der Fall sein. Deswegen wird eine sogenannte Pseudo-Division definiert, die über allen Integritätsringen funktioniert. Gelöst wird dabei nicht die obige Gleichung, sondern die leicht variierte Gleichung

${\displaystyle \alpha p(x)=s(x)q(x)+r(x),}$ 

wobei die Polynome ${\displaystyle p(x)}$  und ${\displaystyle q(x)}$  vorgegeben sind und eine Konstante ${\displaystyle \alpha }$  sowie Polynome ${\displaystyle s(x)}$  und ${\displaystyle r(x)}$  gesucht werden. Auch hier soll wieder der Grad von ${\displaystyle r(x)}$  kleiner als derjenige von ${\displaystyle q(x)}$  sein.

Das Vorgehen ist ähnlich der normalen Polynomdivision. Allerdings werden im Divisionsschritt nicht nur das Polynom ${\displaystyle q(x)}$ , sondern auch ${\displaystyle p(x)}$  mit geeigneten Faktoren multipliziert, um zu erreichen, dass sich die Leitkoeffizienten gegenseitig herauslöschen.

Als Beispiel soll eine Pseudo-Division im Polynomring ${\displaystyle \mathbb {Z} [x]}$  über den ganzen Zahlen durchgeführt werden. Seien

${\displaystyle p(x)=2x^{2}+1\quad {\text{und}}\quad q(x)=5x+5.}$ 

Eine normale Polynomdivision ist hier nicht möglich, da ${\displaystyle 5}$ , der Leitkoeffizienten von ${\displaystyle q}$ , in ${\displaystyle \mathbb {Z} }$  nicht invertierbar ist. Wir können aber ${\displaystyle p}$  mit ${\displaystyle 5}$  multiplizieren. Nun kann man ${\displaystyle q}$  mit ${\displaystyle 2x}$  multipliziert abziehen und erhält

${\displaystyle 5p(x)-2xq(x)=10x^{2}+5-10x^{2}-10x=-10x+5.}$ 

Der Grad von ${\displaystyle -10x+5}$  ist dabei kleiner als derjenige von ${\displaystyle p}$  aber noch nicht kleiner als der von ${\displaystyle q}$ . Ziehen wir nun von diesem Zwischenergebnis ${\displaystyle -2}$ -mal ${\displaystyle q}$  ab, erhalten wir

${\displaystyle -10x+5-(-2)q(x)=-10x+5+10x+10=15.}$ 

Da ${\displaystyle 15}$  als konstantes Polynom einen kleineren Grad als ${\displaystyle q}$  besitzt, sind wir hier fertig. Rückwärts einsetzen ergibt

${\displaystyle 15=(5p(x)-2xq(x))-(-2)q=5p(x)-(2x-2)q(x)}$ 

oder umgeformt

${\displaystyle 5p(x)=(2x-2)q(x)+15.}$ 

Eine Probe bestätigt dies.

Das Vorgehen soll nun noch durch den Algorithmus illustriert werden. Dieser rekursive Algorithmus hat als Argumente zwei Polynome ${\displaystyle p}$  und ${\displaystyle q}$ , wobei ${\displaystyle q}$  nicht das Nullpolynom sein darf, sowie die Variable ${\displaystyle x}$ , bezüglich der die Pseudodivision zu erfolgen hat. Das Ergebnis ist ein Tripel ${\displaystyle (c,s,r)}$  bestehend aus Polynomen ${\displaystyle s}$  und ${\displaystyle r}$  sowie einer Konstanten ${\displaystyle c}$ , so dass ${\displaystyle cp=sq+r}$  und ${\displaystyle \operatorname {grad} (r)<\operatorname {grad} (q)}$  gilt.

   pseudoDivision(p, q, x) =
     if d < 0
     then (1, 0, p)
     else (c * a, c * t + s, r) where
       d       = grad(p, x) - grad(q, x)
       a       = lcoeff(q, x)
       b       = lcoeff(p, x)
       t       = b*xd
       (c,q,r) = pseudoDivision(a*p - t*q, q, x)

Hierbei liefert ${\displaystyle \operatorname {grad} (f,x)}$  den Grad sowie ${\displaystyle \operatorname {lcoeff} (f,x)}$  den Leitkoeffizienten eines Polynomes. Man kann noch weitere Verbesserungen am Algorithmus vornehmen, indem man etwa wie im Beispiel die Multiplikation mit x unterlässt, wenn sie nicht notwendig ist.

Will man bei einer Gleichung

${\displaystyle a_{n}x^{n}+a_{n-1}x^{n-1}+\cdots +a_{2}x^{2}+a_{1}x+a_{0}=0\quad {\text{mit}}\ a_{n}\neq 0}$ 

den Linearfaktor ${\displaystyle (x-x_{1})}$  einer Lösung ${\displaystyle x_{1}}$  abspalten, so ergibt sich das um ein Grad reduzierte Polynom

${\displaystyle b_{n-1}x^{n-1}+b_{n-2}x^{n-2}+\cdots +b_{2}x^{2}+b_{1}x+b_{0}=0}$ 

mit den Koeffizienten

${\displaystyle b_{n-1}=a_{n}}$ 

${\displaystyle b_{n-2}=a_{n-1}+a_{n}\cdot x_{1}}$ 

${\displaystyle b_{n-3}=a_{n-2}+a_{n-1}\cdot x_{1}+a_{n}\cdot {x_{1}}^{2}}$ 

${\displaystyle \cdots }$ 

${\displaystyle b_{0}=a_{1}+a_{2}\cdot x_{1}+a_{3}\cdot {x_{1}}^{2}+\cdots +a_{n}\cdot {x_{1}}^{n-1}}$ 

Beispiel:

Die Polynomgleichung

${\displaystyle 2x^{5}-4x^{4}+4x^{3}+3x^{2}+1{,}5x+0,75=0}$ 

hat die Lösung

${\displaystyle x_{1}=-0{,}4841657}$ 

Das Restpolynom hat also die Koeffizienten

${\displaystyle b_{4}=a_{5}=2}$ 

${\displaystyle b_{3}=a_{4}+a_{5}\cdot x_{1}=-4{,}968331}$ 

${\displaystyle b_{2}=a_{3}+a_{4}\cdot x_{1}+a_{5}\cdot {x_{1}}^{2}=6{,}405496}$ 

${\displaystyle b_{1}=a_{2}+a_{3}\cdot x_{1}+a_{4}\cdot {x_{1}}^{2}+a_{5}\cdot {x_{1}}^{3}=-0{,}101321}$ 

${\displaystyle b_{0}=a_{1}+a_{2}\cdot x_{1}+a_{3}\cdot {x_{1}}^{2}+a_{4}\cdot {x_{1}}^{3}+a_{5}\cdot {x_{1}}^{4}=1{,}549056}$ 

und lautet:

${\displaystyle 2x^{4}-4{,}968331x^{3}+6{,}405496x^{2}-0{,}101321x+1{,}549056=0}$ 

Mit Leitkoeffizient 1 kann schneller mit dem Horner-Schema (zur Funktionswert-Berechnung eines Polynoms) gearbeitet werden. Interessant ist die Umkehrung: man kann mit der Polynomdivision auch Funktionswerte bestimmen. Beispiel: ${\displaystyle p(x)=x^{3}-2x+1}$  mit ${\displaystyle p(3)=22}$ 

Polynomdivision liefert: ${\displaystyle \left(x^{3}-2x+1\right)\colon \left(x-3\right)=x^{2}+3x+7+{\frac {22}{x-3}}}$ 

Nach Multiplikation mit ${\displaystyle (x-3)}$  sieht man, dass der Rest 22 der Funktionswert ${\displaystyle p(3)}$  ist.

Es existiert eine verallgemeinerte Polynomdivision in multivariablen Polynomringen ${\displaystyle K[x_{1},x_{2},\ldots ,x_{n}]}$ , wenn ${\displaystyle K}$  ein Körper ist. Dabei werden einige Abstriche in Kauf genommen, wie beispielsweise die Eindeutigkeit.

• Peter Hartmann: Mathematik für Informatiker. Vieweg+Teubner, 2006, ISBN 3-8348-0096-1, S. 88–90 (Auszug in der Google-Buchsuche)
• Schülerduden Mathematik II. Dudenverlag, 2004, ISBN 3-411-04275-3, S. 327–328
• Charles D. Miller, Margaret L. Lial, David I. Schneider: Fundamentals of College Algebra. 3. überarbeitete Auflage. Scott & Foresman / Little & Brown Higher Education, 1990, ISBN 0-673-38638-4, S. 24–26

• JavaScript berechnet Polynomdivision und erzeugt Übungsaufgaben
• Polynomdivision Rechner mit Rechenweg