PDA

View Full Version : Exponential of (Markov Chain) Transition matrix


NewStudent200
Oct6-09, 01:59 PM
Hi,

I have a (markov chain) transition matrix X which I understand. In particular each row of this matrix sums to 1.

I have used this transition matrix to construct it's generator, Y. I.e. Y is the continuously compounded transition matrix,

X = exp(Y)
X*X = exp(2Y), etc

both X and Y are matrices.

I am told that the sums of Y must sum to 0, but I can not see why this should be the case. Is it obvious?

Many Thanks.

trambolin
Oct7-09, 03:35 PM
Possible hint maybe? e^0=1

NewStudent200
Oct8-09, 02:53 AM
Thanks.

But I have trouble visualizing this for a matrix. Is there aproof somewhere, or a text that you can recommend which gives examples and talks about the applications of taking exponentials or logs of matrices?

Many thanks,

trambolin
Oct8-09, 05:15 AM
It is a matrix exponential (http://en.wikipedia.org/wiki/Matrix_exponential) right ? If you write down the power series for it you will get a pattern.


X = I + Y + \frac{Y^2}{2!} + \ldots

Now if you sum up the rows of X it is 1. On the right hand side you already get 1 from the identity matrix. So all contributions from the remaining terms must be zero right? So I will let you think if your condition is sufficient or necessary.

NewStudent200
Oct8-09, 08:58 AM
Cool. Thanks a lot!

Thinking about matrices a little further. How does one raise a matrix to a decimal power. I know that in the case of an integer power:

X^n = S.M^n.S^-1

where S is the eigen vector matrix and M is the matrix with eigen values along the diagonal. Now if n is non integer, then does this still hold? Could we also do it via:

Y = X^n
ln(Y) = n.ln(X)
Y = exp(n.ln(X))?

much appreciated.

g_edgar
Oct8-09, 09:10 AM
Cool. Thanks a lot!

Thinking about matrices a little further. How does one raise a matrix to a decimal power. I know that in the case of an integer power:

X^n = S.M^n.S^-1

where S is the eigen vector matrix and M is the matrix with eigen values along the diagonal. Now if n is non integer, then does this still hold? Could we also do it via:

Y = X^n
ln(Y) = n.ln(X)
Y = exp(n.ln(X))?

much appreciated.

yes, yes

also:
if X = S M S^{-1} then e^X = S e^M S^{-1}, \log(X) = S \log(M) S^{-1}, etc.