The Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Math Forum » Discussions » sci.math.* » sci.math.num-analysis

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Computation of the matrix exponential
Replies: 6   Last Post: Feb 10, 2013 5:00 AM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]

Posts: 11
Registered: 5/19/10
Computation of the matrix exponential
Posted: May 24, 2010 8:37 PM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply


i have the following problem: given a 3x3 real matrix, compute exp(A).

I need a really fast way to do this. I have searched a bit with
google, but it seems to me that
computing the matrix exponential is not so simple, at least if your
matrix does not have a special
structure (for example A=diagonal matrix).

I have found a simple method that use the diagonalization of A. If A
has 3 distinct eigenvalues, than compute
A=PDP^-1, where P is the matrix of the eigenvectors, D is a diagonal
matrix (whose diagonal elements are
the eigenvalues of A). Than, exp(A) = P exp(D) P^-1. Since P^-1 is
fast enough and exp(D) is simple
to compute, this should be a fast method.

But, the problem is: i am not sure that the matrix A will always have
3 distinct eigenvalues...what happens
if this does not happen? Can i use that formula even if 2 (or all
three) eigenvalues are equal?

Are there any other ways to compute exp(A) in a fast way?

Thank you,

Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum at NCTM 1994-2018. All Rights Reserved.