<< Chapter < Page Chapter >> Page >

Here is some code computing the same Gröbner basis in Sage:

---------------------------------------------------------------------- | Sage Version 4.5.2, Release Date: 2010-08-05                       || Type notebook() for the GUI, and license() for information.        | ----------------------------------------------------------------------sage: R.<x,y,z> = PolynomialRing(QQ,3,order='deglex')  # or degrevlex, lex, etc. sage: f = x^2*y+y^2*z;sage: g = x*y^2+z^2; sage: I = (f,g)*Rsage: I Ideal (x^2*y + y^2*z, x*y^2 + z^2) of Multivariate Polynomial Ring in x, y, z over Rational Fieldsage: gI = I.groebner_basis(); gI [x^2*z^2 + y*z^3, y^3*z - x*z^2, x^2*y + y^2*z, x*y^2 + z^2]sage: h = y^4*z^2+x^3*z^2; sage: h.mod(I)0 sage: h.mod(gI)  # I'm not sure how to get Sage to just do the division algorithm0

Exercises

You may wish to use a computer to do most of the work in the following calculations. One advantage of Macaulay 2 is the getChangeMatrix command to express the generators of the Gröbner basis as linear combinations of the original generators, something I'm not sure how to do in the other programs.

    1. Determine whether of not f = x y 3 - z 2 + y 5 - z 3 is in the ideal
      I = - x 3 + y , x 2 y - z .
    2. Determine whether or not f = x 3 z - 2 y 2 is in the ideal
      I = x z - y , x y + 2 z 2 , y - z .
    1. Find the points on the variety
      V x 2 + y 2 + z 2 - 1 , x 2 + y 2 + z 2 - 2 x , 2 x - 3 y - z .
    2. Find the points on the variety
      V x 2 y - z 3 , 2 x y - 4 z - 1 , z - y 2 , x 3 - 4 z y .
    1. Find an implicit equation for the surface parametrized by:
      x = u t y = 1 - u z = u + t - u t
    2. Find an implicit equation for the surface parametrized by:
      x = t + u y = t 2 + 2 t u z = t 3 + 3 t 2 u

Reduced gröbner bases

Last time we found that

{ x 2 y + y 2 z , x y 2 + z 2 , y 3 z - x z 2 , x 3 z 2 + x y z 3 , x 2 z 2 + y z 3 }

was a Gröbner basis for I = x 2 y + y 2 z , x y 2 + z 2 R [ x , y , z ] in grlex order with x > y > z .

Of course, we'd like to be able to say our Gröbner bases are unique. As a first step, we noticed last time that one element of the Gröbner basis was redundant: x 3 z 2 + x y z 3 = x ( x 2 z 2 + y z ) , so we could remove it and still have a Gröbner basis

{ x 2 y + y 2 z , x y 2 + z 2 , y 3 z - x z 2 , x 2 z 2 + y z 3 } .

More generally, if G is a Gröbner basis with f , g G and L T ( f ) is a multiple of L T ( g ) , then f will be redundant and can be removed, i.e. the set G - { f } is still a Gröbner basis for the same ideal. (Why?)

Of course, multiplying any element of a Gröbner basis by a scalar will give a different Gröbner basis, so if we want uniqueness, we should require that each leading coefficient be 1. Since a monomial ideal certainly has a unique minimal monomial generating set, we might hope that forcing constant leading coefficients and removing redundant elements would be enough to get a Gröbner basis which is unique, but that is not quite the case. The problem is that we could still add a multiple of one generator to another. For example, for any a , b R

{ x 2 y + a x y 2 + y 2 z + a z 2 , x y 2 + z 2 , y 3 z - x z 2 , x 2 z 2 + b y 3 z + y z 3 - b x z 2 }

is another Gröbner basis for the same ideal as above. To avoid non-uniqueness arising in this way, we say that a reduced Gröbner basis G is a Gröbner basis where the leading coefficient of every f G is 1 and no term of any f G is divisible by the leading term of any g G with g f .

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, The art of the pfug. OpenStax CNX. Jun 05, 2013 Download for free at http://cnx.org/content/col10523/1.34
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'The art of the pfug' conversation and receive update notifications?

Ask