Computeralgebra

Die Computeralgebra ist das Teilgebiet der Mathematik und Informatik, das sich mit der automatisierten symbolischen Manipulation algebraischer Ausdrücke beschäftigt.

Hauptziel der Computeralgebra ist es, durch konservative Rechnungen algebraische Ausdrücke umzuformen und eine möglichst kompakte Darstellung zu erzielen. Die Wertigkeit und Präzision der Gleichung wird dabei nicht angetastet. Rundungen oder Näherungen werden nicht zugelassen. Eine Nebenbedingung ist hierbei, die verwendeten Algorithmen und Methoden effizient zu gestalten.

Die Disziplinen Mathematik und Informatik sind in der Computeralgebra eng miteinander verwoben, einerseits über die Komplexitätstheorie für die Analyse, andererseits über die Softwaretechnik für die praktische Umsetzung computeralgebraischer Algorithmen.

Einen Schwerpunkt bildet das exakte Rechnen mit ganzen, rationalen und algebraischen Zahlen sowie mit Polynomen über diesen Zahlenräumen. Eine weitere Anwendung ist das symbolische Lösen von Gleichungen aller Arten, das symbolische Summieren von Reihen, die symbolische Berechnung von Grenzwerten und das symbolische Differenzieren und Integrieren (auch Algebraische Integration genannt).

Praktische Anwendung erfahren solche Ergebnisse durch den Einsatz effizienter Algorithmen bei der Entwicklung und Verbesserung von Computeralgebrasystemen, die die rechnergestützte Manipulation algebraischer Ausdrücke ermöglichen. Diese Systeme sind auch ein immer wichtiger werdendes Werkzeug für Mathematiker und Naturwissenschaftler verschiedener Fachrichtungen.

Anders als in der Numerik, wo mit Gleitkomma-Approximationen der auftretenden Zahlen gerechnet wird, hat die Computeralgebra den Anspruch, stets exakt zu rechnen. Entsprechend ist es grundsätzlich notwendig, Anforderungen an die zugrundeliegenden Strukturen (in der Regel Gruppen, Ringe oder Körper) zu spezifizieren.

Alle endlichen Gruppen lassen sich im Computer darstellen, für unendliche Gruppen gibt es unter bestimmten Voraussetzungen Algorithmen, etwa für polyzyklische Gruppen.

Ein Ring heißt berechenbar (oder „effektiv“), wenn folgende Bedingungen gelten:

• Die Elemente des Rings können auf einem Computer dargestellt werden, insbesondere muss die Darstellung der Elemente endlich sein,
• Gleichheit zwischen zwei Elementen kann in endlicher Zeit entschieden werden,
• es gibt Algorithmen, mit denen die Ringoperationen „+“ und „·“ durchgeführt werden können.

Berechenbar sind etwa:

• die natürlichen Zahlen N {\displaystyle \mathbb {N} }
• die ganzen Zahlen Z {\displaystyle \mathbb {Z} }
• die rationalen Zahlen Q {\displaystyle \mathbb {Q} }
• alle Formen von endlichen Körpern.

Aus einem berechenbaren Ring R {\displaystyle R} lassen sich weitere berechenbare Ringe konstruieren:

• Polynomringe über R {\displaystyle R}
• Rationale Funktionen über R {\displaystyle R}
• Matrizen über R {\displaystyle R}
• Alle endlichen algebraischen Erweiterungen der obigen Körper,
• Alle endlichen transzendenten Erweiterungen der obigen Körper.

In der Computeralgebra werden neben den Elementen der zugrundeliegenden Bereiche noch weitere „formale“ Objekte betrachtet wie etwa

• Integrale,
• Reihen,
• (formale) Potenzreihen,
• Differential- und Differenzengleichungen (bzw. -operatoren).

Hierbei geht es in der Regel nicht um die Berechnung von Zahlwerten, sondern beispielsweise um die Bestimmung von „geschlossenen Formeln“ als Lösungen.

Bei den Anwendungen mathematischer Methoden in Naturwissenschaft und Technik stehen traditionell numerische Methoden im Vordergrund. Mit den symbolischen Methoden der Computeralgebra haben sich neue Anwendungsgebiete eröffnet, bei denen es auf exakte Lösungen ankommt und bei denen strukturmathematische Überlegungen, z. B. zur Beschreibung von Symmetrien, eingehen; ferner die Behandlung von Problemen, die von unbestimmten Parametern abhängen.

Dazu gehören etwa:

• In der Physik und ihren technologischen Anwendungen: Gemischt symbolisch-numerische Berechnungen komplexer Probleme in der Himmelsmechanik, der Hochenergiephysik (Feynman-Integrale) und der Relativitätstheorie (Differentialgeometrie); Integration und Lösung von Differentialgleichungen in geschlossener Form; symbolische Berechnungen in den Algebren der Quantenmechanik; Klassifikation höherdimensionaler kristallographischer Gruppen zur Beschreibung von inkommensurabel modulierten Strukturen, Quasikristallen und magnetischen Strukturen.
• In der Chemie: Anwendungen der Darstellungstheorie auf die Klassifikation von Graphen chemischer Verbindungen, insbesondere von Isomeren; Lösung großer Gleichungssysteme zur Bestimmung chemischer Reaktionsgleichgewichte bei variablen Reaktionsbedingungen, z. B. bei Verbrennungsprozessen und der Abgasregulierung.
• In der Informationssicherheit: Algebraische Methoden der Fehlererkennung und -korrektur bei der Nachrichtenübertragung; kryptographische Kodierung von vertraulichen Nachrichten durch Public-Key-Verfahren mit Methoden der Zahlentheorie und der algebraischen Gruppen; Verifikation von Sicherheitsmechanismen (Protokollen).
• In der Robotik: Bewegungsplanung und Regulierung autonomer Roboter z. B. in der Raumfahrt; zylindrisch algebraische Zellenzerlegung des R n {\displaystyle \mathbb {R} ^{n}} ; geometrische Bildverarbeitung beim Maschinensehen.
• Im computergestützten Entwurf (CAD): Flexible Inferenzsysteme für die geometrische Modellierung parametrisierter Probleme; Konstruktion von Übergangsflächen.
• In der Kontrolltheorie: Untersuchung der Stabilität und Sicherheit von Kontrollsystemen mit Rückkopplung.
• In der Genforschung: Klassifikation von DNA-Strukturen.
• In der Ausbildung: Computeralgebrasysteme versprechen eine Verbesserung des Mathematikunterrichts, da es durch ihren Einsatz möglich wird, sich mehr auf die Unterrichtsinhalte zu konzentrieren. Ferner ist die Behandlung realistischerer anwendungsbezogener Aufgaben möglich.

Die Methoden der Computeralgebra erlauben in diesen Anwendungsbereichen eine automatische Behandlung von Problemen, die sonst nur mühsam mit Ad-hoc-Ansätzen angegangen werden konnten. Das große Potential dieser Methoden ist dabei noch lange nicht ausgeschöpft. Kontinuierliche Förderung dieser Anwendungen und der zugrundeliegenden Algorithmen kombiniert mit immer leistungsstärkerer Hardware werden dem Gebiet Computeralgebra weitere Entwicklungsmöglichkeiten geben.

Will man die Zeitkomplexität von Aufgaben und Algorithmen zur exakten Arithmetik mit ganzen Zahlen klassifizieren, so muss zunächst ein Rechnermodell zugrunde gelegt werden. Eine Diskussion diverser Rechnermodelle findet sich im Kapitel Komplexität. Ein relativ eingängiges Modell ist die Mehrband-Turingmaschine, eine Variante der klassischen Turingmaschine, die mehrere Bänder mit je einem Schreib-/Lesekopf besitzt. Für Komplexitätsabschätzungen mit der Landau-Notation wird bei Bedarf unter der Bezeichnung log {\displaystyle \operatorname {log} } ein Logarithmus zu einer nicht spezifizierten Basis B > 1 {\displaystyle B>1} verwendet. Als Maß für die Zeit wird die Zahl der benötigten Bitoperationen gewählt, die in Landau-Notation von der Bitlänge des Inputs abhängig gemacht wird.

Die präzise mathematische Angabe von (Bit-)Komplexitäten für die exakte Arithmetik mit ganzen Zahlen muss zunächst mit der genauen Festlegung der Bitlänge einer ganzen Zahl starten: Ist die Zahl a ∈ Z {\displaystyle a\in \mathbb {Z} } nicht null, so wird

L ( a ) := ⌊ log 2 ⁡ | a | ⌋ {\displaystyle L(a):=\lfloor \log _{2}|a|\rfloor }

gesetzt; zusätzlich wird L ( 0 ) := 1 {\displaystyle L(0):=1} definiert. Für die konkrete Speicherung einer ganzen Zahl wird zusätzlich mindestens noch ein Bit für das Vorzeichen benötigt.

Die Aufgaben der Vorzeichenbestimmung sgn ⁡ ( a ) , {\displaystyle \operatorname {sgn} (a),} der Berechnung des Negativen − a {\displaystyle -a} sowie der Betragsbildung | a | {\displaystyle |a|} sind alle in linearer Zeit O ( l ) {\displaystyle O(l)} mit l = L ( a ) {\displaystyle l=L(a)} durchführbar; die Addition a + b {\displaystyle a+b} sowie der Vergleich zweier Zahlen a < b {\displaystyle a<b} sind in linearer Zeit O ( l ) {\displaystyle O(l)} mit l = max ( L ( a ) , L ( b ) ) {\displaystyle l=\max(L(a),L(b))} zu bewältigen. Der n-Shift 2 n ⋅ a {\displaystyle 2^{n}\cdot a} ist in O ( n + l ) {\displaystyle O(n+l)} durchführbar.

Ein nicht-triviales Ergebnis der Computeralgebra ist die Erkenntnis, dass die Multiplikation a ⋅ b {\displaystyle a\cdot b} wesentlich schneller als in O ( l 2 ) {\displaystyle O\left(l^{2}\right)} (was dem naiven Multiplikationsalgorithmus entspricht) lösbar ist. Eine Beschleunigung erreichte zunächst Anatoli Karazuba mit dem Karazuba-Algorithmus; dieser wurde dann als ein Spezialfall einer noch allgemeineren Algorithmenfamilie erkannt, die unter den Begriff Toom-Cook-Algorithmus subsumiert werden. Bahnbrechend war dann der von Arnold Schönhage und Volker Strassen 1971 vorgestellte auf diskreten Fourier-Transformationen basierende Schönhage-Strassen-Algorithmus, für den die Autoren selbst eine Komplexität von

O ( l ⋅ log l ⋅ log log l ) {\displaystyle O(l\cdot \log \,l\cdot \log \,\log \,l)}

nachwiesen. Bedenkt man, wie aufwendig der „naive“ Multiplikationsalgorithmus ist, so erscheint diese Komplexität unglaublich „schnell“. Da der Algorithmus allerdings ziemlich komplex und schwierig programmierbar ist, gibt es bis heute keine effiziente Implementierung in einem Computeralgebrasystem.

Da die Komplexität der Integer-Multiplikation in der gesamten Computeralgebra von absolut tragender Bedeutung ist, wurde hierfür eine Kurznotation

ψ ( l ) := O ( l ⋅ log l ⋅ log log l ) {\displaystyle \psi (l):=O(l\cdot \log \,l\cdot \log \,\log \,l)}

eingeführt. Ausgestattet mit dieser „schnellen“ Integer-Multiplikation kann nun der Katalog der Grundrechenarten für die Arithmetik in Z {\displaystyle \mathbb {Z} } wie folgt vervollständigt werden: Die Aufgabe der Berechnung von a n {\displaystyle a^{n}} ist in O ( ψ ( n l ) ) {\displaystyle O\left(\psi (nl)\right)} durchführbar; für die (simultane) Berechnung der Binomialkoeffizienten ( n 0 ) ⋯ ( n n ) {\displaystyle {n \choose 0}\cdots {n \choose n}} wird O ( n ⋅ ψ ( n ) ) {\displaystyle O\left(n\cdot \psi (n)\right)} benötigt. Die ganzzahlige Division a / b {\displaystyle a/b} (mit Quotient und Rest als Ergebnis) benötigt

O ( l l m ⋅ ψ ( l m ) ) {\displaystyle O\left({\frac {l}{l_{m}}}\cdot \psi (l_{m})\right)} .

Die Berechnung des größten gemeinsamen Teilers gcd ⁡ ( a , b ) {\displaystyle \operatorname {gcd} (a,b)} benötigt

O ( ( l l m + log ⁡ l m ) ⋅ ψ ( l m ) ) {\displaystyle O\left(\left({\frac {l}{l_{m}}}+\log l_{m}\right)\cdot \psi (l_{m})\right)} .

In gleicher Komplexität ist auch gcdex ⁡ ( a , b ) {\displaystyle \operatorname {gcdex} (a,b)} berechenbar, d. h. die Kofaktoren u , v {\displaystyle u,v} mit gcd ⁡ ( a , b ) = u a + v b {\displaystyle \operatorname {gcd} (a,b)=ua+vb} werden mitberechnet.

Bevor exakte Arithmetik in Q {\displaystyle \mathbb {Q} } konkret durchgeführt werden kann, muss erst eine kanonische Darstellung (Repräsentation) rationaler Zahlen gefunden werden; dieses Problem tauchte bei der exakten Arithmetik der ganzen Zahlen noch nicht auf. Rationale Zahlen sind Äquivalenzklassen „bedeutungsgleicher“ Brüche aus ganzen Zahlen; zum Beispiel sind 1 2 {\displaystyle {\tfrac {1}{2}}} und 2 4 {\displaystyle {\tfrac {2}{4}}} unterschiedliche Repräsentanten der gleichen rationalen Zahl.

Die gängigste kanonische Darstellung rationaler Zahlen wird festgelegt, indem alle gemeinsamen Teiler aus Zähler und Nenner herausgekürzt werden: Jede rationale Zahl ist dann eindeutig durch einen gekürzten Bruch

p q {\displaystyle {\frac {p}{q}}\quad } mit p ∈ Z , {\displaystyle \quad p\in \mathbb {Z} ,} q ∈ N {\displaystyle q\in \mathbb {N} } und ggT ⁡ ( p , q ) = 1 {\displaystyle \operatorname {ggT} (p,q)=1}

repräsentierbar. Ist diese Festlegung einmal getroffen, so beinhaltet jede elementare Operation in Q {\displaystyle \mathbb {Q} } wie Addition und Multiplikation auch notwendigerweise die Aufgabe des Herauskürzens des größten gemeinsamen Teilers aus Zähler und Nenner des Ergebnisbruches.

Dank der Resultate der exakten Arithmetik in Z {\displaystyle \mathbb {Z} } sind damit die Operationen a + b ,   a − b ,   a ⋅ b ,   a / b {\displaystyle a+b,\ a-b,\ a\cdot b,\ a/b} alle in der Komplexität

O ( ψ ( l ) ⋅ log ⁡ l ) {\displaystyle O\left(\psi (l)\cdot \log l\right)}

durchführbar. Von der Hoffnung, die Addition rationaler Zahlen in linearer Komplexität bewerkstelligen zu können, muss man hierbei Abschied nehmen.

Glücklicherweise kann die Berechnung des größten gemeinsamen Teilers sehr effizient mit Hilfe des Euklidischen Algorithmus berechnet werden. Der Euklidische Algorithmus spielt an vielen Stellen in der Computeralgebra in wechselnden Varianten eine tragende Rolle.

Es genügt hierbei, die Arithmetik in Z [ x ] {\displaystyle \mathbb {Z} \left[x\right]} zu betrachten, da Operationen mit rationalen Polynomen in naheliegender Weise durch jeweiliges Ausklammern der Hauptnenner auf Operationen mit ganzzahligen Polynomen zurückgeführt werden können. Für Polynome f ∈ Z [ x ] {\displaystyle f\in \mathbb {Z} [x]} definiert man die (Koeffizienten-)Länge L ( f ) {\displaystyle L(f)} als das Maximum der Längen der einzelnen Koeffizienten.

Für die folgende Laufzeitentabelle wichtiger Berechnungsprobleme in Z [ x ] {\displaystyle \mathbb {Z} [x]} setzen wir folgendes voraus:

• f , g ∈ Z [ x ] {\displaystyle f,g\in \mathbb {Z} \left[x\right]}
• vom Grad d f = deg ⁡ f {\displaystyle d_{f}=\deg f} bzw. d g = deg ⁡ g , {\displaystyle d_{g}=\deg g,} ferner d = max ( d f , d g ) , d m = min ( d f , d g ) {\displaystyle d=\max(d_{f},d_{g}),\,d_{m}=\min(d_{f},d_{g})} ,
• von der Länge l f = L ( f ) {\displaystyle l_{f}=L(f)} bzw. l g = L ( g ) , {\displaystyle l_{g}=L(g),} ferner l = max ( l f , l g ) {\displaystyle l=\max(l_{f},l_{g})} und l m = min ( l f , l g ) {\displaystyle l_{m}=\min(l_{f},l_{g})} ,
• daneben sei a ∈ Z {\displaystyle a\in \mathbb {Z} } mit Bitlänge l a = L ( a ) {\displaystyle l_{a}=L(a)} .

Dann führen die (gemäß Bitkomplexität) schnellsten bislang bekannten Algorithmen zu folgender Laufzeittabelle:

Operation	Komplexität als O ( … ) {\displaystyle O(\dots )}	bei ungleichen Eingabegrößen
f + g {\displaystyle f+g}	d l {\displaystyle d\,l}	d m l {\displaystyle d_{m}\,l}
f ⋅ g {\displaystyle f\cdot g}	ψ ( d ( l + log ⁡ d ) ) {\displaystyle \psi (d(l+\log d))}	d d m ψ ( d m ( ( d − d m + 1 ) l + log ⁡ d ) ) {\displaystyle {\frac {d}{d_{m}}}\,\psi (d_{m}((d-d_{m}+1)l+\log d))}
f / g {\displaystyle f/g} Division mit Rest	ψ ( d ( l + log ⁡ d ) ) {\displaystyle \psi (d(l+\log d))}	d l d m l m ψ ( d m ( l m + log ⁡ d m ) ) {\displaystyle {\frac {dl}{d_{m}l_{m}}}\,\psi (d_{m}(l_{m}+\log d_{m}))}
prem ⁡ ( f , g ) {\displaystyle \operatorname {prem} (f,g)}	− {\displaystyle -}	d d m ψ ( d m ( ( d − d m + 1 ) l + log ⁡ d ) ) {\displaystyle {\frac {d}{d_{m}}}\,\psi (d_{m}((d-d_{m}+1)l+\log d))}
f k {\displaystyle f^{k}}	ψ ( k d ( l + log ⁡ d ) ) {\displaystyle \psi (kd(l+\log d))}	− {\displaystyle -}
f ( a x ) {\displaystyle f(ax)} Skalierung	d ψ ( l + d l a ) {\displaystyle d\,\psi (l+dl_{a})}	− {\displaystyle -}
a d f ( x a ) {\displaystyle a^{d}f\left({\frac {x}{a}}\right)} Skalierung	d ψ ( l + d l a ) {\displaystyle d\,\psi (l+dl_{a})}	− {\displaystyle -}
f ( 2 x ) {\displaystyle f(2x)} Skalierung	d ( l + d ) {\displaystyle d(l+d)}	− {\displaystyle -}
2 d f ( x 2 ) {\displaystyle 2^{d}f\left({\frac {x}{2}}\right)} Skalierung	d ( l + d ) {\displaystyle d(l+d)}	− {\displaystyle -}

• Symbolische Mathematik

Deutschsprachige Literatur:

• Michael Kaplan: Computeralgebra. Springer, Berlin u. a. 2005, ISBN 3-540-21379-1.
• Wolfram Koepf: Computeralgebra. Eine algorithmisch orientierte Einführung. Springer, Berlin u. a. 2006, ISBN 3-540-29894-0.
• Attila Pethö: Algebraische Algorithmen. Herausgegeben von Michael Pohst. Vieweg, Braunschweig u. a. 1999, ISBN 3-528-06598-2.

Englischsprachige Literatur:

• Peter Bürgisser, Michael Clausen, Mohammad Amin Shokrollahi: Algebraic Complexity Theory (= Die Grundlehren der mathematischen Wissenschaften in Einzeldarstellungen. 315). Springer, Berlin u. a. 1997, ISBN 3-540-60582-7 (Rezension).
• Arjeh M. Cohen, Hans Cuypers, Hans Sterk: Some Tapas of Computer Algebra (= Algorithms and Computation in Mathematics. 4). Springer, Berlin u. a. 1999, ISBN 3-540-63480-0.
• David Cox, John Little, Donal O’Shea: Ideals, varieties, and algorithms. 2nd Edition. Springer, New York NY ua. a. 1997, ISBN 0-387-94680-2.
• Joachim von zur Gathen, Jürgen Gerhard: Modern Computer Algebra. 2nd Edition. Cambridge University Press, Cambridge u. a. 2003, ISBN 0-521-82646-2.
• Keith O. Geddes, Stephen R. Czapor, George Labahn: Algorithms for Computer Algebra. Kluwer, Boston MA u. a. 1992, ISBN 0-7923-9259-0.

• Fachgruppe Computeralgebra, gemeinsame Fachgruppe von GI, DMV und GAMM
• Sieveking: Algorithmen der Computeralgebra. (PDF; 666 kB) Skript, Goethe-Universität Frankfurt a. M., Archivlink, abgerufen am 26. Februar 2022