home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #18 / NN_1992_18.iso / spool / sci / math / 10259 < prev    next >
Encoding:
Internet Message Format  |  1992-08-13  |  1.2 KB
  1. Path: sparky!uunet!psinntp!kepler1!andrew
  2. From: andrew@rentec.com (Andrew Mullhaupt)
  3. Newsgroups: sci.math
  4. Subject: Re: Help - non-integral power of a matrix?
  5. Keywords: matrix
  6. Message-ID: <1173@kepler1.rentec.com>
  7. Date: 13 Aug 92 14:08:06 GMT
  8. References: <Aug.10.15.45.34.1992.26563@clam.rutgers.edu>
  9. Organization: Renaissance Technologies Corp., Setauket, NY.
  10. Lines: 20
  11.  
  12. In article <Aug.10.15.45.34.1992.26563@clam.rutgers.edu> gonzalez@clam.rutgers.edu (Ralph Gonzalez) writes:
  13. >
  14. >Hi.  Does anyone know of an algorithm to find a non-integral
  15. >power of a matrix, e.g. A^.5 or A^1.3? Thus, A^2 is the same
  16. >as AxA and A^0 is the identity.
  17.  
  18. Yes. _Run_ do not walk, to the nearest copy of _Matrix Computations_ by
  19. Golub and van Loan. Read carefully the part about the matrix exponential.
  20. Do not fail to notice the ominous title of one of the more important
  21. references on this subject: "Nineteen Dubious Ways to Compute the Matrix
  22. Exponential", and thus chastened, implement one of the recommended algorithms
  23. (Shur decomposition or Squaring, scaling and Pade).
  24.  
  25. I've had good success with the squaring, scaling & Pade approach.
  26.  
  27. Later,
  28. Andrew Mullhaupt
  29.  
  30.  
  31. P.S. this should be in a FAQ, right?
  32.