MathGroup Archive 2002

[Date Index] [Thread Index] [Author Index]

Search the Archive

RE: How to compute a MatrixPower using: A^n = P D^n Inverse[P]

  • To: mathgroup at smc.vnet.net
  • Subject: [mg34989] RE: [mg34976] How to compute a MatrixPower using: A^n = P D^n Inverse[P]
  • From: "DrBob" <majort at cox-internet.com>
  • Date: Tue, 18 Jun 2002 02:48:29 -0400 (EDT)
  • Reply-to: <drbob at bigfoot.com>
  • Sender: owner-wri-mathgroup at wolfram.com

MatrixForm[A = {{3, 1, 0}, {1, 2, -1}, {0, -1, 3}}]
{evalues, evectors} = Eigensystem[A]
MatrixForm[ch = Transpose[evectors]]
Inverse[ch].A.ch // MatrixForm
MatrixForm[d = DiagonalMatrix[evalues]]
ch.d.Inverse[ch] == A

This only works when A has a basis of eigenvectors, of course.

Bobby Treat

-----Original Message-----
From: J. Guillermo Sanchez [mailto:guillerm at usal.es] 
To: mathgroup at smc.vnet.net
Subject: [mg34989] [mg34976] How to compute a MatrixPower using: A^n = P D^n
Inverse[P] 


I have the matrix

A == {{3,1,0},{1,2,-1},{0,-1,3}}

For educational purpose I would like to evaluate

A^n (* I mean MatrixPower[A,n]*)

using the following matrix property

A^n == P D^n Inverse[P]  (*D mean Diagonal Matrix *)

How can I do with Mathematica? (Methods to obtain P and D)

Thanks

Guillermo Sanchez






  • Prev by Date: Colors
  • Next by Date: RE: How to compute a MatrixPower using: A^n = P D^n Inverse[P]
  • Previous by thread: Re: How to compute a MatrixPower using: A^n = P D^n Inverse[P]
  • Next by thread: RE: How to compute a MatrixPower using: A^n = P D^n Inverse[P]