• 0 Vote(s) - 0 Average
• 1
• 2
• 3
• 4
• 5
 Matrix Operator Method Gottfried Ultimate Fellow Posts: 757 Threads: 116 Joined: Aug 2007 08/27/2007, 12:43 PM bo198214 Wrote:Gottfried Wrote:bo198214 Wrote:Did you anyway already realize that instead of $\exp(t\cdot\log(A))$ you can directly use the binomial series for computation? $A^t = \sum_{n=0}^\infty \left(t\\n\right) (A-I)^n$No, but it looks very good. I'll give it a deeper look, thanks for the hint!But it seems as if A must have a diagonal of 1. So probably you have to convert A into a Jordan normal form first as you must do with the logarithm (or how do you compute the matrix logarithm?) Second try. You say $\hspace{24} A^t = \sum_{n=0}^\infty \left(t\\n\right) (A-I)^n$ Let $\hspace{24} A = Q * D * Q^{-1}$ if A is diagonalizable. Also, since $\hspace{24} I = Q * I * Q^{-1}$ it is also $\hspace{24} (A - I) = Q *(D - I) * Q^{-1}$ Then the above means, for example t = 3 $\hspace{24} A^3 = Q * ( 1*(D -I)^0 + 3*(D-I)^1 + 3*(D-I)^2 + 1*(D-I)^3 ) * Q^{-1}$ Let's look at the parenthese only as the diagonal matrix E. The entries e_k can directly be computed: by binomial-theorem $\hspace{24} e_0 = ((d_0 - 1) + 1)^3$ $\hspace{24} e_1 = ((d_1 - 1) + 1)^3$ $\hspace{24}$ $\hspace{24} ...$ $\hspace{24} e_0 = d_0^3$ $\hspace{24} e_1 = d_1^3$ $\hspace{24}$ $\hspace{24} ...$ and $\hspace{24} E = D^3$ So the rhs is $\hspace{24} Q * E * Q^{-1} = Q * D^3 * Q^{-1}$ and the lhs is $\hspace{24} A^3 = Q * D^3 * Q^{-1}$ and the both sides are equal. Here was no assumption needed concerning the actual eigenvalues of A, except the assumption, that it is diagonalizable at all. It is perhaps important to note, that symbolically we may write formulae with equal eigenvalues even if a numerical solver would not produce a solution because of the mentioned indeterminacy of the eigenvectors. But this does not affect the analytical description. Well, here I assumed A is diagonalizable (with no assumptions about the actual value of its eigenvalues). If it is not, then at least it is possible to derive the Jordan-canonical form, and in A = Q*J*Q^-1 we get J non-diagonal, but in its Jordan-canonical form. All what was said about the diagonal form is still valid in principle, only that powers and functions on J have not the simple form as in the true diagonal case. But still it is a simpler computation than with the original matrix, especially if infinite series are assumed. Gottfried Gottfried Helms, Kassel « Next Oldest | Next Newest »

 Messages In This Thread Matrix Operator Method - by Gottfried - 08/12/2007, 08:08 PM RE: Matrix Operator Method - by bo198214 - 08/13/2007, 04:15 AM RE: Matrix Operator Method - by jaydfox - 08/13/2007, 05:40 AM RE: Matrix Operator Method - by Gottfried - 08/13/2007, 09:22 AM RE: Matrix Operator Method - by bo198214 - 08/14/2007, 03:43 PM RE: Matrix Operator Method - by Gottfried - 08/14/2007, 04:15 PM RE: Matrix Operator Method - by bo198214 - 08/26/2007, 12:18 AM RE: Matrix Operator Method - by Gottfried - 08/26/2007, 11:24 AM RE: Matrix Operator Method - by bo198214 - 08/26/2007, 11:39 AM RE: Matrix Operator Method - by Gottfried - 08/26/2007, 04:22 PM RE: Matrix Operator Method - by Gottfried - 08/26/2007, 10:54 PM RE: Matrix Operator Method - by bo198214 - 08/27/2007, 08:29 AM RE: Matrix Operator Method - by Gottfried - 08/27/2007, 11:04 AM RE: Matrix Operator Method - by bo198214 - 08/27/2007, 11:35 AM RE: Matrix Operator Method - by Gottfried - 08/27/2007, 11:58 AM RE: Matrix Operator Method - by bo198214 - 08/27/2007, 12:13 PM RE: Matrix Operator Method - by Gottfried - 08/27/2007, 01:19 PM RE: Matrix Operator Method - by Gottfried - 08/27/2007, 02:29 PM RE: Matrix Operator Method - by bo198214 - 08/27/2007, 02:36 PM RE: Matrix Operator Method - by Gottfried - 08/27/2007, 03:09 PM RE: Matrix Operator Method - by bo198214 - 08/27/2007, 07:15 PM RE: Matrix Operator Method - by Gottfried - 08/27/2007, 08:15 PM RE: Matrix Operator Method - by bo198214 - 08/29/2007, 05:28 PM RE: Matrix Operator Method - by Gottfried - 08/27/2007, 12:43 PM RE: Matrix Operator Method - by Gottfried - 10/08/2007, 12:11 PM RE: Matrix Operator Method - by Gottfried - 10/14/2007, 09:32 PM RE: Matrix Operator Method - by Gottfried - 04/04/2008, 09:41 AM RE: Matrix Operator Method - by Gottfried - 04/17/2008, 09:21 PM RE: Matrix Operator Method - by bo198214 - 04/25/2008, 03:39 PM RE: Matrix Operator Method - by Gottfried - 04/26/2008, 06:09 PM RE: Matrix Operator Method - by bo198214 - 04/26/2008, 06:47 PM RE: Matrix Operator Method - by Gottfried - 04/18/2008, 01:55 PM RE: Matrix Operator Method - by Gottfried - 07/08/2008, 06:46 AM Diagonalization for dxp/basic facts/Pari-routine - by Gottfried - 08/08/2008, 01:12 PM Exact entries for T-tetration Bell-matrix - by Gottfried - 09/24/2008, 08:22 PM RE: Exact entries for T-tetration Bell-matrix - by bo198214 - 09/26/2008, 07:30 AM RE: Exact entries for T-tetration Bell-matrix - by Gottfried - 09/26/2008, 09:56 AM

 Possibly Related Threads... Thread Author Replies Views Last Post Half-iterates and periodic stuff , my mod method [2019] tommy1729 0 175 09/09/2019, 10:55 PM Last Post: tommy1729 A fundamental flaw of an operator who's super operator is addition JmsNxn 4 6,464 06/23/2019, 08:19 PM Last Post: Chenjesu Tetration and Sign operator tetration101 0 307 05/15/2019, 07:55 PM Last Post: tetration101 2 fixpoints , 1 period --> method of iteration series tommy1729 0 1,345 12/21/2016, 01:27 PM Last Post: tommy1729 Tommy's matrix method for superlogarithm. tommy1729 0 1,541 05/07/2016, 12:28 PM Last Post: tommy1729 [split] Understanding Kneser Riemann method andydude 7 7,336 01/13/2016, 10:58 PM Last Post: sheldonison [2015] New zeration and matrix log ? tommy1729 1 2,975 03/24/2015, 07:07 AM Last Post: marraco Kouznetsov-Tommy-Cauchy method tommy1729 0 1,917 02/18/2015, 07:05 PM Last Post: tommy1729 Problem with cauchy method ? tommy1729 0 1,758 02/16/2015, 01:51 AM Last Post: tommy1729 Regular iteration using matrix-Jordan-form Gottfried 7 8,090 09/29/2014, 11:39 PM Last Post: Gottfried

Users browsing this thread: 1 Guest(s)