Singular
https://www.singular.uni-kl.de/forum/

newbie: how to find minimal polynomial of a matrix?
https://www.singular.uni-kl.de/forum/viewtopic.php?f=10&t=1898
Page 1 of 1

Author:  alhutch [ Mon Jan 03, 2011 11:20 pm ]
Post subject:  newbie: how to find minimal polynomial of a matrix?

Thank you for Singular.
I want to do what seems a simple job: to find the minimal polynomial of a square matrix.
The matrix entries are polynomials in the ring Q[x], so the minimal polynomial will also be in this ring.
The algorithmic complications seem to be
- although all the entities involved are in commutative rings, the calculations of powers of the original matrix are made in the ring of all such matrices over Q[x], which is not commutative;
- the tools of Singular seem designed (for good reason) to hide the information I want. They tell the user whether an element is in an ideal, but not (obviously) how the element can be expressed in terms of the ideal's generators.
Any guidance welcome.

Author:  malex [ Tue Jan 04, 2011 8:52 pm ]
Post subject:  Re: newbie: how to find minimal polynomial of a matrix?

1. don't you need to adjoin another commutative variable to your ring (e.g. Q[x] -> Q[x][T])?
2. matrix multiplication is generally non-commutative, no?
3. one uses 'lift' for "how the element can be expressed in terms of the ideal's generators"

4. i'd suggest you to start with characteristic polynomial (of M) as follows:
a. compute D = det(M) in our original ring,
b. switch to enriched ring (e.g. new variable T), and imap D and M there
c. compute det(D * E * T - M) (where E is the identity matrix)
it can be used for testing and checking your results...


ps: Happy New Year

Author:  alhutch [ Wed Jan 05, 2011 11:16 am ]
Post subject:  Re: newbie: how to find minimal polynomial of a matrix?

Thanks, malek.
This makes things much clearer.
Alan.

Page 1 of 1 All times are UTC + 1 hour [ DST ]
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group
http://www.phpbb.com/