Gröbnerbasis

Eine Gröbnerbasis (nach Bruno Buchberger, 1965) bzw. Standardbasis (nach Heisuke Hironaka, 1964) ist ein endliches Erzeugendensystem zu einem Ideal I {\displaystyle I} im Polynomring K [ X 1 , , X n ] {\displaystyle K[X_{1},\ldots ,X_{n}]} über dem Körper K {\displaystyle K} , das besonders gut dafür geeignet ist, zu entscheiden, ob ein gegebenes Polynom zum Ideal gehört oder nicht.

Buchberger entwickelte diese 1965 in seiner Dissertation bei Wolfgang Gröbner und benannte sie nach ihm.

Motivation: Idealzugehörigkeitsproblem und verallgemeinerte Polynomdivision

Das Idealzugehörigkeitsproblem

Es sei I = ( g 1 , , g n ) K [ X 1 , , X k ] {\displaystyle I=(g_{1},\ldots ,g_{n})\subseteq K[X_{1},\dots ,X_{k}]} , also das Ideal I {\displaystyle I} von den Polynomen g 1 , , g n {\displaystyle g_{1},\ldots ,g_{n}} erzeugt, dann gehört ein Polynom f {\displaystyle f} genau dann zu I {\displaystyle I} , wenn sich f {\displaystyle f} als Linearkombination

f = a 1 g 1 + + a n g n {\displaystyle f=a_{1}g_{1}+\dots +a_{n}g_{n}}

mit Polynomen a 1 , , a n K [ X 1 , , X k ] {\displaystyle a_{1},\dots ,a_{n}\in K[X_{1},\dots ,X_{k}]} schreiben lässt.

Man kann nun versuchen, so eine Darstellung mit Hilfe der Polynomdivision mit Rest zu finden, indem man so lange dividiert, bis der erhaltene Rest verschwindet und damit die Division aufgeht:

f = a 1 g 1 + r 1 {\displaystyle f=a_{1}g_{1}+r_{1}}
f = a 1 g 1 + a 2 g 2 + r 2 {\displaystyle f=a_{1}g_{1}+a_{2}g_{2}+r_{2}}

Dabei treten aber zwei Probleme auf:

  1. Bei multivariaten Polynomen p K [ X 1 , , X k ] {\displaystyle p\in K[X_{1},\ldots ,X_{k}]} lassen sich die Terme nicht mehr ohne Weiteres „der Größe nach“ ordnen, was für die Polynomdivision aber notwendig ist.
  2. Können wir denn überhaupt sicher sein, dass durch wiederholte Polynomdivision immer der Rest 0 erscheint?

Das erste Problem lässt sich recht einfach lösen, indem eine Monomordnung gewählt wird. Dann können die Terme jedes Polynomes nach dieser Ordnung angeordnet werden; vor allem können wir nun vom Leitterm LT {\displaystyle \operatorname {LT} } eines Polynoms sprechen, also des bezüglich der Monomordnung größten Monoms mit seinem Koeffizienten. Als Nachteil bleibt, dass die Anordnung der Monome und auch das Ergebnis der Polynomdivisionen von der Wahl der Monomordnung abhängen.

Das zweite Problem ist schwieriger, denn es ist tatsächlich nicht lösbar, wenn die erzeugenden Polynome fest vorgegeben sind. Es kann nur gelöst werden, indem das Erzeugendensystem passend geändert wird. Eine Gröbnerbasis stellt sich als ein geeignetes Erzeugendensystem heraus.

Verallgemeinerte Polynomdivision

Die Aufgabe der verallgemeinerten Polynomdivision ist nun also: Für ein Polynom f {\displaystyle f} und mehrere Polynome g 1 , , g n {\displaystyle g_{1},\dots ,g_{n}} sollen Polynome a 1 , , a n , r {\displaystyle a_{1},\dots ,a_{n},r} gefunden werden, die die Gleichung

f = a 1 g 1 + + a n g n + r {\displaystyle f=a_{1}g_{1}+\dots +a_{n}g_{n}+r}

erfüllen.

Dazu bietet sich folgendes Vorgehen an:

  1. Beginne mit a 1 = = a n = r = 0 {\displaystyle a_{1}=\dots =a_{n}=r=0} .
  2. Falls f = 0 {\displaystyle f=0} , brich ab, sonst vergleiche der Reihe nach alle LT ( g i ) {\displaystyle \operatorname {LT} (g_{i})} , ob sie LT ( f ) {\displaystyle \operatorname {LT} (f)} teilen.
  3. Falls etwa a LT ( g i ) = LT ( f ) {\displaystyle a\operatorname {LT} (g_{i})=\operatorname {LT} (f)} , so ersetze a i {\displaystyle a_{i}} durch a i + a {\displaystyle a_{i}+a} und f {\displaystyle f} durch f a g i {\displaystyle f-ag_{i}} und gehe zu Schritt 2.
  4. Falls LT ( f ) {\displaystyle \operatorname {LT} (f)} von keinem LT ( g i ) {\displaystyle \operatorname {LT} (g_{i})} geteilt wird, so ersetze r {\displaystyle r} durch r + LT ( f ) {\displaystyle r+\operatorname {LT} (f)} und f {\displaystyle f} durch f LT ( f ) {\displaystyle f-\operatorname {LT} (f)} und gehe zu Schritt 2.

Dann ist in jedem Schritt die Gleichung

f urspr = a 1 g 1 + + a n g n + r + f aktuell {\displaystyle f_{\text{urspr}}=a_{1}g_{1}+\dots +a_{n}g_{n}+r+f_{\text{aktuell}}}

erfüllt, und schließlich, wenn f aktuell = 0 {\displaystyle f_{\text{aktuell}}=0} , ist die gewünschte Darstellung gefunden.

Mit dieser Division haben wir das Problem wie gewünscht auf ein kleineres Polynom reduziert, denn es ist f I {\displaystyle f\in I} genau dann, wenn r I {\displaystyle r\in I} . Falls r = 0 {\displaystyle r=0} , ist das klar, und f I {\displaystyle f\in I} . Ist aber r 0 {\displaystyle r\neq 0} , können wir auf diesem Weg nicht entscheiden, ob r I {\displaystyle r\in I} oder r I {\displaystyle r\not \in I} :

Beispiel: Seien g 1 = x {\displaystyle g_{1}=x} , g 2 = x 2 + 1 {\displaystyle g_{2}=x^{2}+1} und f = g 1 + g 2 = x 2 + x + 1 {\displaystyle f=g_{1}+g_{2}=x^{2}+x+1} . Testet man die erzeugenden Polynome der Reihe nach (also in aufsteigender Reihenfolge der Indizes) und wendet die beschriebene Division an, erhält man:

f = x g 1 + ( x + 1 ) = ( x + 1 ) g 1 + 1 {\displaystyle f=xg_{1}+(x+1)=(x+1)g_{1}+1}

Offenbar gilt aber f = g 1 + g 2 I {\displaystyle f=g_{1}+g_{2}\in I} , also auch r = 1 I {\displaystyle r=1\in I} (nämlich 1 = x g 1 + g 2 {\displaystyle 1=-xg_{1}+g_{2}} ). Man kann also im Allgemeinen aus r 0 {\displaystyle r\neq 0} nicht folgern, dass r I {\displaystyle r\not \in I} und damit f I {\displaystyle f\not \in I} gilt.

Definition

Ein Erzeugendensystem G = ( g 1 , , g n ) {\displaystyle G=(g_{1},\ldots ,g_{n})} von I {\displaystyle I} ist eine Gröbnerbasis (bezüglich einer Monomordnung < {\displaystyle <} ) von I {\displaystyle I} , falls für jedes Polynom f I { 0 } {\displaystyle f\in I\setminus \{0\}} ein g i {\displaystyle g_{i}} existiert, dessen Leitmonom das Leitmonom von f {\displaystyle f} teilt.

Eine Gröbnerbasis G {\displaystyle G} heißt reduziert, falls alle g G {\displaystyle g\in G} normiert sind, und kein Monom von g {\displaystyle g} durch die Leitterme der anderen Gröbnerbasispolynome dargestellt werden kann, also kein Monom von g {\displaystyle g} in { LT ( g ) : g G { g } } {\displaystyle \langle \{\operatorname {LT} (g'):g'\in G\setminus \{g\}\}\rangle } liegt. Man kann zeigen, dass jedes Ideal (für eine gegebene Monomordnung) eine eindeutig bestimmte reduzierte Gröbnerbasis besitzt.

Anwendungen

Das Konzept von Gröbnerbasen gibt zunächst eine Lösung des Idealzugehörigkeitsproblems. Damit verbunden lassen sich aber auch andere Probleme lösen (oder zumindest vereinfachen).

Lösung des Idealzugehörigkeitsproblems

Wird mit dem oben beschriebenen Verfahren eine nicht weiter reduzierbare Darstellung

f = k 1 g 1 + + k n g n + r {\displaystyle f=k_{1}g_{1}+\ldots +k_{n}g_{n}+r}

bezüglich einer Gröbnerbasis G = ( g 1 , , g n ) {\displaystyle G=(g_{1},\ldots ,g_{n})} des Ideals I {\displaystyle I} gefunden, so gilt f I {\displaystyle f\in I} genau dann, wenn r I {\displaystyle r\in I} . Da aber G {\displaystyle G} eine Gröbnerbasis ist, gilt das wiederum genau dann, wenn r = 0 {\displaystyle r=0} , da nach Annahme kein Leitmonom eines g i {\displaystyle g_{i}} das Leitmonom von r {\displaystyle r} teilt.

Mit dem Buchberger-Algorithmus können Gröbnerbasen berechnet werden. Damit ist das Problem, ob ein Polynom zu einem Ideal gehört oder nicht, von Computeralgebrasystemen lösbar.

Beispiel: Wenn wie im Beispiel oben die erzeugenden Polynome g 1 = x , g 2 = x 2 + 1 {\displaystyle g_{1}=x,g_{2}=x^{2}+1} gegeben sind, sowie f = g 1 + g 2 = x 2 + x + 1 {\displaystyle f=g_{1}+g_{2}=x^{2}+x+1} , dann hatte die Polynomdivision Rest 1 {\displaystyle 1} geliefert.

Wenden wir zunächst den Buchberger-Algorithmus auf dieses Beispiel an, so erhalten wir die (nicht reduzierte) Gröbnerbasis ( x , x 2 + 1 , 1 ) {\displaystyle (x,x^{2}+1,-1)} von I {\displaystyle I} . Bezüglich dieses Divisors ist die Division hier noch nicht abgeschlossen, denn es ist mit g 3 = 1 {\displaystyle g_{3}=-1} :

f = x g 1 + ( x + 1 ) = ( x + 1 ) g 1 + 1 = ( x + 1 ) g 1 + ( 1 ) g 3 {\displaystyle f=xg_{1}+(x+1)=(x+1)g_{1}+1=(x+1)g_{1}+(-1)g_{3}}

Wir sehen hier auch, dass die Darstellung bezüglich der Gröbnerbasis nicht eindeutig ist (da f = g 1 + g 2 = ( x + 1 ) g 1 g 3 {\displaystyle f=g_{1}+g_{2}=(x+1)g_{1}-g_{3}} ), sondern von der Reihenfolge der erzeugenden Polynome und der gewählten Monomordnung abhängt.

Polynomiale Gleichungssysteme

Ein polynomiales Gleichungssystem besteht aus endlich vielen Gleichungen f 1 ( x ) = 0 , , f k ( x ) = 0 {\displaystyle f_{1}(x)=0,\ldots ,f_{k}(x)=0} , wobei f 1 , , f k K [ x 1 , , x n ] {\displaystyle f_{1},\ldots ,f_{k}\in K[x_{1},\ldots ,x_{n}]} gegebene Polynome in n {\displaystyle n} Unbestimmten über einem Körper K {\displaystyle K} und deren gemeinsame Nullstellen x K n {\displaystyle x\in K^{n}} gesucht sind. Die Lösungsmenge { x K n : f 1 ( x ) = = f k ( x ) = 0 } {\displaystyle \{x\in K^{n}:f_{1}(x)=\cdots =f_{k}(x)=0\}} eines solchen Gleichungssystems beschreibt eine algebraische Varietät. Diese ist offenbar gleich { x K n : f ( x ) = 0 f I } {\displaystyle \{x\in K^{n}:f(x)=0\quad \forall f\in I\}} , wobei I = ( f 1 , , f k ) {\displaystyle I=(f_{1},\ldots ,f_{k})} das von den Polynomen erzeugte Ideal ist.

Soll ein bestimmtes polynomiales Gleichungssystem gelöst werden, genügt es also das Ideal zu betrachten, das von den Polynomen erzeugt wird. Dann kann es hilfreich sein, mit Hilfe des Buchberger-Algorithmus und einer geeigneten lexikographischen Monomordnung eine reduzierte Gröbnerbasis zu finden. Dann bleibt zwar das Problem, die Nullstellen dieser Polynome zu finden (z. B. näherungsweise durch numerische Verfahren), aber die zu untersuchenden Polynome haben immerhin eine kleinere Anzahl an Variablen und kleineren Grad.

Beispiel: Welche Lösungen ( x , y , z ) R 3 {\displaystyle (x,y,z)\in \mathbb {R} ^{3}} hat das folgende Gleichungssystem?

x 2 + y 2 + z 2 1 = 0 {\displaystyle x^{2}+y^{2}+z^{2}-1=0}
x 2 y + z 2 = 0 {\displaystyle x^{2}-y+z^{2}=0}
x z = 0 {\displaystyle x-z=0}

Mit Hilfe des Computers erhält man für die lexikographische Monomordnung x > y > z {\displaystyle x>y>z} die (reduzierte) Gröbnerbasis g 1 = x z , g 2 = y + 2 z 2 {\displaystyle g_{1}=x-z,g_{2}=-y+2z^{2}} und g 3 = z 4 + 1 2 z 2 1 4 {\displaystyle g_{3}=z^{4}+{\frac {1}{2}}z^{2}-{\frac {1}{4}}} . Die gesuchten Lösungen sind also genau die Lösungen des einfacheren Gleichungssystems

x = z {\displaystyle x=z}
y = 2 z 2 {\displaystyle y=2z^{2}}
z 4 + 1 2 z 2 1 4 = 0 {\displaystyle z^{4}+{\frac {1}{2}}z^{2}-{\frac {1}{4}}=0}

Und wir sehen, dass die Lösungsmenge aus nur zwei Punkten besteht: { ( z , 2 z 2 , z ) : z = ± 1 2 5 1 } {\displaystyle \left\{(z,2z^{2},z):z=\pm {\frac {1}{2}}{\sqrt {{\sqrt {5}}-1}}\right\}} (die zwei weiteren Lösungen ± 1 2 5 1 {\displaystyle \pm {\frac {1}{2}}{\sqrt {-{\sqrt {5}}-1}}} sind in R {\displaystyle \mathbb {R} } nicht enthalten).

Gleichheit von Idealen und algebraischen Varietäten

Da (zu einer gegebenen Monomordnung) die reduzierte Gröbnerbasis eines Ideals eindeutig bestimmt ist, kann die Gleichheit von Idealen untersucht werden, indem (zu irgendeiner Monomordnung) die reduzierten Gröbnerbasen gebildet werden. Die Ideale sind genau dann gleich, wenn diese reduzierten Gröbnerbasen gleich sind.

Auf diese Weise kann man auch rein rechnerisch die Gleichheit von algebraischen Varietäten untersuchen, da diese durch ihre erzeugenden Ideale eindeutig bestimmt sind. Sind die reduzierten Gröbnerbasen gleich, dann sind die erzeugenden Ideale und damit auch die erzeugten Varietäten gleich.

Literatur

  • Heisuke Hironaka: Resolution of singularities of an algebraic variety over a field of characteristic zero, In: Annals of Mathematics 79 (1964), Nr. 1, S. 109–326
  • Bruno Buchberger: Ein Algorithmus zum Auffinden der Basiselemente des Restklassenringes nach einem nulldimensionalen Polynomideal (PDF; 1,8 MB). Österreich, Universität Innsbruck, Diss., 1965
  • Bruno Buchberger: Ein algorithmisches Kriterium für die Lösbarkeit eines algebraischen Gleichungssystems (PDF; 475 kB). Aequationes Mathematicae 4 (1970): 374–383.
  • T. Becker, B. V. Weispfenning: Gröbner bases: a computational approach to commutative algebra. Graduate Texts in Mathematics: readings in mathematics. Bd. 141, New York: Springer Verlag, 1993
  • David Cox; Little, John; O’Shea, Donald: Ideals, varieties, and algorithms. New York: Springer-Verlag, 1997
  • Joachim von z. Gathen, Jürgen Gerhard: Modern Computer Algebra. Cambridge University Press, 1999