Re: How to compute a MatrixPower using: A^n = P D^n Inverse[P]
- To: mathgroup at smc.vnet.net
- Subject: [mg34987] Re: [mg34976] How to compute a MatrixPower using: A^n = P D^n Inverse[P]
- From: Andrzej Kozlowski <andrzej at platon.c.u-tokyo.ac.jp>
- Date: Tue, 18 Jun 2002 02:48:27 -0400 (EDT)
- Sender: owner-wri-mathgroup at wolfram.com
Since your matrix is symmetric and has distinct eigenvalues this is rather easy. Let's just do one case, say n=3. First let's compute: In[1]:= A={{3,1,0},{1,2,-1},{0,-1,3}}; In[2]:= B=MatrixPower[A,3]; Next we find the eigenvalues and eigenvectors of B. In[3]:= {eigenvals, eigenvecs} = Eigensystem[B]; We normalize the eigenvectors and compute the change of basis matrix: In[4]:= q = Map[1/Sqrt[#.#]&, eigenvecs]*eigenvecs; the diagonal matrix has the eigenvalues along the diagonal: In[5]:= d=DiagonalMatrix[eigenvals]; Its transpose = its inverse In[6]:= Transpose[q]==Inverse[q] Out[6]= True Finally we check the formula: In[7]:= B==Inverse[q].d.q Out[7]= True Andrzej Kozlowski Toyama International University JAPAN http://platon.c.u-tokyo.ac.jp/andrzej/ On Monday, June 17, 2002, at 04:26 PM, J. Guillermo Sanchez wrote: > > 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 > > > >