20.05.2014 Views

link to my thesis

link to my thesis

link to my thesis

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

2.5. THE BUCHBERGER ALGORITHM 21<br />

The algorithm above gives a set of genera<strong>to</strong>rs G = {g 1 ,...,g t } which is a Gröbner<br />

basis for the ideal I. If there is any member of G whose leading term can be divided<br />

by the leading term of any other member of G, one can eliminate this member from<br />

the basis <strong>to</strong> get a minimal size of the basis.<br />

Example 2.12 (continued). Suppose we start with the set of polynomials {f 1 ,f 2 } =<br />

{x 2 − 2y 2 ,xy − 3} of Example 2.12 where the monomials are ordered with respect<br />

<strong>to</strong> a lexicographic ordering. In Example 2.12 the Gröbner basis for I = 〈f 1 ,f 2 〉 is<br />

given by {g 1 ,g 2 } = {2y 4 − 9, 3x − 2y 3 }. We now show step-by-step how <strong>to</strong> obtain this<br />

Gröbner basis for the ideal I using the Buchberger algorithm as shown in Algorithm<br />

1 (each line starts with a number denoting the line of pseudocode in Algorithm 1).<br />

1 The set G is initialized as G = F = {f 1 ,f 2 }.<br />

2 The set M is initialized as M = {{f i ,f j } : f i ,f j ∈ G, f i f j } = {{f 1 ,f 2 }}.<br />

Note that {f 2 ,f 1 } need not be included.<br />

3 The set M is not empty: new iteration starts.<br />

4 The pair {p, q} can only be chosen as {f 1 ,f 2 }.<br />

5 The set M becomes: M = M −{f 1 ,f 2 } = ∅.<br />

6 The S-polynomial of {f 1 ,f 2 } is computed using LT (f 1 )=x 2 , LT (f 2 )=xy and<br />

J = LCM(LM(f 1 ),LM(f 2 )) = LCM(x 2 ,xy)=x 2 y as: S = S-polynomial(f 1 ,<br />

f 2 )=<br />

J<br />

LT (f f 1) 1 −<br />

J<br />

LT (f f 2) 2 = x2 y<br />

x<br />

f 2 1 − x2 y<br />

xy f 2 = yf 1 − xf 2 =3x − 2y 3 = g 2 .<br />

7 The NormalForm of g 2 is computed as h = NormalForm(g 2 ,G)=g 2 =3x−2y 3<br />

because:<br />

3x is not divisible by LM(f 1 )=x 2 nor divisible by LM(f 2 )=xy.<br />

−2y 3 is not divisible by LM(f 1 )=x 2 nor divisible by LM(f 2 )=xy.<br />

Therefore the polynomial g 2 is already in normal form.<br />

8 Because h = g 2 0, the sets M and G are extended.<br />

9 The set M is extended <strong>to</strong> M = {{f 1 ,g 2 }, {f 2 ,g 2 }}.<br />

10 The set G becomes G = {f 1 ,f 2 ,g 2 } = {x 2 − 2y 2 ,xy− 3, 3x − 2y 3 }.<br />

3 The set M is not empty: new iteration starts.<br />

4 The pair {p, q} is chosen as {f 1 ,g 2 }.<br />

5 The set M becomes: M = M −{f 1 ,g 2 } = {{f 2 ,g 2 }}.<br />

6 The S-polynomial of {f 1 ,g 2 } is computed using LT (f 1 )=x 2 , LT (g 2 )=3x and<br />

J = LCM(LM(f 1 ),LM(g 2 )) = LCM(x 2 ,x)=x 2 as: S = S-polynomial(f 1 ,<br />

g 2 )=<br />

J<br />

LT (f 1) f 1 −<br />

2<br />

3 xy3 − 2y 2 .<br />

2<br />

7 The NormalForm of<br />

J<br />

LT (g 2) g 2 = x2<br />

x 2 f 1 − x2<br />

3x (3x − 2y3 )=x 2 − 2y 2 − x 2 + 2 3 xy3 =<br />

3 xy3 − 2y 2 is computed as h = NormalForm( 2 3 xy3 −<br />

2y 2 ,G)= 2 3 xy3 − 2y 2 − 2 3 (xy − 3)y2 =0.<br />

8 h =0.<br />

9 Because h = 0 the set M is not extended.<br />

10 Because h = 0 the set G is not extended.<br />

3 The set M is not empty: new iteration starts.<br />

4 The pair {p, q} is chosen as {f 2 ,g 2 }.<br />

5 The set M becomes: M = M −{f 2 ,g 2 } = ∅.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!