Biquadratische Gleichung
aus Wikipedia, der freien Enzyklopädie
Eine biquadratische Gleichung oder polynomiale Gleichung 4. Grades (auch quartische Gleichung genannt) hat die Form
mit komplexen Koeffizienten und
.
(Im Fall läge nur eine Gleichung niedrigeren Grades vor.)
Nach dem Fundamentalsatz der Algebra lässt sich die Gleichung bis auf die Reihenfolge eindeutig in die Form
bringen, wobei x1,x2,x3 und x4 die – nicht notwendigerweise verschiedenen – vier Lösungen der Gleichung sind.
Ist B = 0 und D = 0, dann lässt sich die Gleichung durch Substitution auf eine quadratische Gleichung zurückführen. Diese Spezialform wird manchmal in Lehrbüchern als biquadratische Gleichung bezeichnet, was fälschlicherweise den Eindruck erweckt, dass nur dieser Spezialfall biquadratisch heißt. Es ist jedoch jede Gleichung 4. Grades eine biquadratische (also auch solche mit kubischem und linearem Glied), denn der Name "biquadratisch" bedeutet lediglich, dass die höchste vorkommende Potenz die 4. Potenz ist (nämlich das Biquadrat x4).
Inhaltsverzeichnis |
[Bearbeiten] Geschichte
Die erste geschlossene Lösung der biquadratischen Gleichung fand der italienische Mathematiker Lodovico Ferrari (1522-1565). Diese Lösung veröffentlichte sein Lehrer Gerolamo Cardano 1545 in dem Werk Ars magna de Regulis Algebraicis. Eine weitere Lösungsmethode mit unterschiedlichem Ansatz wurde von Leonhard Euler 1738 in Sankt Petersburg publiziert, in dem Bestreben, eine allgemeine Lösungsformel auch für Gleichungen höherer Grade zu finden. Daß dies unmöglich ist, wurde von Niels Henrik Abel 1824 bewiesen.
[Bearbeiten] Satz
Voraussetzung: Gegeben sei eine quartische Gleichung Ax4 + Bx3 + Cx2 + Dx + E = 0 mit .
Aussage: Dann kann man ihre Lösungen wie folgt angeben:
- Algorithmus 1 (frei nach Ferrari; entnommen aus der englischen Wikipedia Quartic equation):
- Zunächst wird die Gleichung mit der Substitution
- dahingehend vereinfacht, dass der kubische Koeffizient verschwindet.
- Mit den Festlegungen
- reduziert sich die Gleichung zu
- Falls β = 0, dann löse u4 + αu2 + γ = 0 und berechne
.
- Sonst:
wählen, um alle Lösungen zu erhalten
- q.e.f.
[Bearbeiten] Beweis
(konstruktiv)
bis zur Erstellung der deutschen Übersetzung möge die englische Version der Herleitung hinreichen: en:Quartic equation.
[Bearbeiten] Hilfssatz A
Sei P und Q wie oben, dann gilt:
Beweis: Limes-Rechenregeln (man beachte die höheren Potenzen im Zähler; man beachte, dass P und Q unabhängig von einander gegen 0 streben können):
- 0
Es folgt die Richtigkeit der Fallunterscheidung bezüglich der Bedingung U=0.
[Bearbeiten] Spezialformen
[Bearbeiten] B=0 und D=0
In dem Fall lässt sich die Gleichung einfach auf eine quadratische Gleichung zurückführen:
-
- Substituieren: x2 durch z:
,
- Finden: Die Lösungen zu Gleichung (99) z1 und z2 finden: siehe Quadratische Gleichung,
- Rück-subsitituieren:
:
,
,
,
.
- Substituieren: x2 durch z:
[Bearbeiten] E=0
In diesem Fall ist x1 = 0 eine Lösung der Gleichung. Dann kann man den Faktor (x − x1) also (x − 0) ausklammern und erhält die Gleichung
Die Lösungen der quartischen Gleichung sind dann 0 und die drei Lösungen der kubischen Gleichung
-
.
[Bearbeiten] reelle Koeffizienten
Sind alle Koeffizienten reell, lassen sich Fallunterscheidungen für die möglichen Lösungen angeben. Dies beruht auf folgender Tatsache: Ist die nicht-reelle Zahl Nullstelle eines beliebigen Polynoms mit reellen Koeffizienten, so ist es auch die konjugiert komplexe Zahl
(Beweis). Bei der Zerlegung des zugehörigen Polynoms ergibt das Produkt der beiden Faktoren
ein quadratisches Polynom mit reellen Koeffizienten, nämlich
.
Also lässt sich jedes Polynom mit reellen Koeffizienten unabhängig von seinem Grad in lineare und quadratische Faktoren mit reellen Koeffizienten zerlegen. Es gibt für die quartische Gleichung also drei Möglichkeiten:
- Die Gleichung hat vier reelle Lösungen. Sie zerfällt in vier Linearfaktoren mit reellen Koeffizienten.
- Die Gleichung hat zwei reelle und zwei konjugiert komplexe Lösungen. Sie zerfällt in zwei Linearfaktoren und einen quadratischen Faktor mit reellen Koeffizienten.
- Die Gleichung hat zwei Paare konjugiert komplexer Lösungen. Sie zerfällt in zwei quadratische Faktoren mit reellen Koeffizienten.
[Bearbeiten] vier reelle Lösungen
Unter den Lösungen können einzelne Lösungen oder solche mit einer Vielfachheit 2, 3 oder 4 sein. (Erläuterung).
Im einzelnen gibt es diese Möglichkeiten:
- eine Lösung mit Vielfachheit 4
- Beispiel:
, zerlegt
- hat die vierfache Lösung
- eine Lösung mit Vielfachheit 3 und eine einfache Lösung
- Beispiel:
, zerlegt
- hat die dreifache Lösung
und die einfache Lösung
- zwei Lösungen, jeweils mit Vielfachheit 2
- Beispiel:
, zerlegt
- hat die zweifache Lösung
und die zweifache Lösung
- eine Lösung mit Vielfachheit 2 und zwei einfache Lösungen
- Beispiel:
, zerlegt
- hat die zweifache Lösung
und die einfachen Lösungen
- vier einfache Lösungen
- Beispiel:
, zerlegt
- hat die einfachen Lösungen
[Bearbeiten] zwei reelle und zwei konjugiert komplexe Lösungen
Auch hier kann die reelle Lösung mit Vielfachheit 2 auftreten. Es gibt also diese beiden Möglichkeiten:
- eine reelle Lösung mit Vielfachheit 2 und zwei konjugiert komplexe Lösungen
- Beispiel:
, zerlegt
- oder mit reellem quadratischen Faktor
- hat die zweifache Lösung
und die konjugiert komplexen Lösungen
- zwei einfache reelle Lösungen und zwei konjugiert komplexe Lösungen
- Beispiel:
, zerlegt
- oder mit reellem quadratischen Faktor
- hat die einfachen Lösungen
und die konjugiert komplexen Lösungen
[Bearbeiten] zwei Paare konjugiert komplexer Lösungen
Hier gibt es diese beiden Möglichkeiten:
- zwei konjugiert komplexe Lösungen mit Vielfachheit 2
- Beispiel:
, zerlegt
- oder mit zwei reellen quadratischen Faktoren
- hat die zweifachen konjugiert komplexen Lösungen
- zwei Paare einfacher konjugiert komplexer Lösungen
- Beispiel:
, zerlegt
- oder mit zwei reellen quadratischen Faktoren
- hat die konjugiert komplexen Lösungen
und
[Bearbeiten] Pseudo-Code
A, B, C, D, E, a, b, c, P, Q, R, U, y, x1, x2, x3 und x4 können komplexe Zahlen sein.
Demzufolge müssen die Funktionen „sqrt“, „exp“ und „ln“ auf ebensolchen Zahlen definiert sein.
Für B=0 und D=0
- dann
-
- Algorithmus aus dem Artikel Quadratische Gleichung auf Az2 + Cz + E = 0 anwenden
- und x1,2,3,4 berechnen
-
-
- x1 = sqrt(z1);
- x2 = -sqrt(z1);
- x3 = sqrt(z2);
- x4 = -sqrt(z2);
-
- sonst
-
- a = -3*B*B/(8*A*A) + C/A;
- b = B*B*B/(8*A*A*A) - B*C/(2*A*A) + D/A;
- c = -3*B*B*B*B/(256*A*A*A*A) + C*B*B/(16*A*A*A) - B*D/(4*A*A) + E/A;
- wenn b = 0
- dann
-
- löse u4 + au2 + c = 0
- und berechne
-
- x1=u1-B/(4*A);
- x2=u2-B/(4*A);
- x3=u3-B/(4*A);
- x4=u4-B/(4*A);
- sonst
-
-
- P = -a*a/12 - c;
- Q = -a*a*a/108 + a*c/3 - b*b/8;
- R = Q/2 + sqrt(Q*Q/4+P*P*P/27);
- wenn R = 0
- dann U = 0;
- sonst U = exp(ln(R)/3);
- wenn U = 0
- dann (prüfe, ob P=0 gilt) y = -5*a/6;
- sonst y = -5*a/6 + P/(3*U) - U;
- W = sqrt(a+2*y)
- x1 = -B/(4*A) + (W + sqrt(-(a+2*y) - 2*(a + b/W)))/2;
- x2 = -B/(4*A) + (-W + sqrt(-(a+2*y) - 2*(a - b/W)))/2;
- x3 = -B/(4*A) + (W - sqrt(-(a+2*y) - 2*(a + b/W)))/2;
- x4 = -B/(4*A) + (-W - sqrt(-(a+2*y) - 2*(a - b/W)))/2;
-
Es empfiehlt sich die Probe zu machen...
[Bearbeiten] Beispiele
Beispiel für
> calc C-style arbitrary precision calculator (version 2.11.5t4.5) Calc is open software. For license details type: help copyright [Type "exit" to exit, or "help" for help.] > A=1 > B=10 > C=-6 > D=60 > E=36 > a=-3*B^2/A^2/8 + C/A > a -43.5 > b=B^3/8/A^3 - B*C/2/A^2 + D/A > b 215 > c=-3*B^4/256/A^4 + C*B^2/16/A^3 - B*D/4/A^2 + E/A > c -268.6875 > P = -a*a/12 - c > P 111 > Q = -a*a*a/108 + a*c/3 - b*b/8 > Q -1120 > R = Q/2 + sqrt(Q*Q/4+P*P*P/27) > R 43.53376044758258411234 > U = exp(ln(R)/3) > U 3.51783442380909981526 > y = -5*a/6 + P/(3*U) - U > y ~43.24999999999999999999 > x1 = -B/(4*A) + (sqrt(a+2*y) + sqrt(-(a+2*y) - 2*(a + b/sqrt(a+2*y) )))/2 > x1 .77871926215100032617+2.32241174444907628892i > A*x1^4+B*x1^3+C*x1^2+D*x1+E ~.00000000000000000007-~.00000000000000000073i > x2 = -B/(4*A) + (-sqrt(a+2*y) + sqrt(-(a+2*y) - 2*(a - b/sqrt(a+2*y) )))/2 > x2 ~-.54483004754633870880 > A*x2^4+B*x2^3+C*x2^2+D*x2+E ~.00000000000000000014 > x3 = -B/(4*A) + (sqrt(a+2*y) - sqrt(-(a+2*y) - 2*(a + b/sqrt(a+2*y) )))/2 > x3 .77871926215100032617-2.32241174444907628892i > A*x3^4+B*x3^3+C*x3^2+D*x3+E ~.00000000000000000007+~.00000000000000000073i > x4 = -B/(4*A) + (-sqrt(a+2*y) - sqrt(-(a+2*y) - 2*(a - b/sqrt(a+2*y) )))/2 > x4 ~-11.01260847675566194353 > A*x4^4+B*x4^3+C*x4^2+D*x4+E ~-.00000000000000000335
Beispiel für
> calc C-style arbitrary precision calculator (version 2.11.5t4.5) Calc is open software. For license details type: help copyright [Type "exit" to exit, or "help" for help.] > A=1 > B=0 > C=1 > D=sqrt(-8/27) > E=-1/12 > a=-3*B^2/A^2/8 + C/A > a 1 > b=B^3/8/A^3 - B*C/2/A^2 + D/A > b .54433105395181735515i > c=-3*B^4/256/A^4 + C*B^2/16/A^3 - B*D/4/A^2 + E/A > c ~-.08333333333333333333 > P = -a*a/12 - c > P 0 > Q = -a*a*a/108 + a*c/3 - b*b/8 > Q ~-.00000000000000000000 > R = Q/2 + sqrt(Q*Q/4+P*P*P/27) > R ~-.00000000000000000000 > # U = exp(ln(R)/3) > U=0 > U 0 > # y = -5*a/6 + P/(3*U) - U > y = -5*a/6 > y ~-.83333333333333333333 > x1 = -B/(4*A) + (sqrt(a+2*y) + sqrt(-(a+2*y) - 2*(a + b/sqrt(a+2*y) )))/2 > x1 ~1.22474487139158904909i > A*x1^4+B*x1^3+C*x1^2+D*x1+E ~-.00000000000000000000 > x2 = -B/(4*A) + (-sqrt(a+2*y) + sqrt(-(a+2*y) - 2*(a - b/sqrt(a+2*y) )))/2 > x2 -~.40824829041868291618i > A*x2^4+B*x2^3+C*x2^2+D*x2+E ~-.00000000000000000000 > x3 = -B/(4*A) + (sqrt(a+2*y) - sqrt(-(a+2*y) - 2*(a + b/sqrt(a+2*y) )))/2 > x3 -~.40824829046386301636i > A*x3^4+B*x3^3+C*x3^2+D*x3+E ~-.00000000000000000000 > x4 = -B/(4*A) + (-sqrt(a+2*y) - sqrt(-(a+2*y) - 2*(a - b/sqrt(a+2*y) )))/2 > x4 -~.40824829050904311654i > A*x4^4+B*x4^3+C*x4^2+D*x4+E ~-.00000000000000000000
[Bearbeiten] Siehe auch
[Bearbeiten] Literatur
siehe Cardanische Formeln
[Bearbeiten] Links
TU-Freiberg – Verschiedene Lösungswege für Gleichungen 4. Grades – (pdf) gut verständlich