• 0 Vote(s) - 0 Average
• 1
• 2
• 3
• 4
• 5
 Matrix Operator Method Gottfried Ultimate Fellow Posts: 758 Threads: 117 Joined: Aug 2007 08/13/2007, 09:22 AM (This post was last modified: 08/13/2007, 02:53 PM by Gottfried.) *** I just tried my first latex-code in the hope to improve readability *** Matrix method for tetration Well - there is not much special here. 1) ----------------------------------------------------- Assume, you denote the summation of a powerseries $\hspace{24} s(x) = \sum_{k=0}^\infty x^k * a_k$ as a vector-product: $\hspace{24} s(x) = rowvector(1,x,x^2,x^3,...) * colvector(a_0,a_1,a_2,a_3,...)$ then it may be useful for the further analysis, to give names for such vectors. I say $\hspace{24} V(x) =colvector(1,x,x^2,x^3,...)$ Let the other vector be denoted as $\hspace{24} A = colvector(a_0,a_1,a_2,a_3,...)$ then $\hspace{24} s(x) = V(x)\sim * A$ Example: if $\hspace{24} A = colvector(a_0,a_1,a_2,a_3,...) = (1,1,\frac1{2!},\frac1{3!},\frac1{4!},...)$ then V(x)~ * A represents simply the exponential-series in x, so $\hspace{24} e^x = V(x)\sim * A$ 2) ------------------------------------------------------- Now, if we want iteration, to get e^(e^x) that means, that we simply need to set e^x=y and use y as paraemter for the "powerseries"-vector (or better: "vandermonde"-vector, thus the letter "V" at V(x)) $\hspace{24} e^y = V(y)\sim * A = e^{e^x}$ But what is V(y) now? It contains the powers of y, or the powers of e^x. But by the previous formula $\hspace{24} y = e^x = V(x)\sim * A$ we only have the "first power" of e^x. How to obtain the other powers too? To make in in one formula, we should -instead of a single vectorproduct- have a full matrix-product, which provides all required powers of e^x as well, and in the same one shot. Something like Code:.   V(y)~ = V(x)~ [A0,  A1,    A2,     A3    ,  ... ]         =       [1 , e^x , (e^x)^2, (e^x)^3,  ... ] where A1 is our already known vector A. A2, for instance, is then simply $\hspace{24} A2 = colvector(1, 2, \frac{2^2}{2!}, \frac{2^3}{3!}, \frac{2^4}{4!},...)$ since Code:.    (e^x)^2 = e^(2x) = sum (k=0..inf) (2x)^k/k!                     = sum (k=0..inf) 2^k * x^k / k!                     = sum (k=0..inf) 2^k/k! * x^kand written as a vector-product it is $\hspace{24} \left(e^x\right)^2 = V(x)\sim * colvector(\frac{2^0}{0!},\frac{2^1}{1!},\frac{2^2}{2!}, ....)$ This is completely analoguous for all powers of (e^x). So Code:.   A0 = colvector(0^0/0!, 0^1/1!, 0^2/2!,...)   A1 = colvector(1^0/0!, 1^1/1!, 1^2/2!,...)   A2 = colvector(2^0/0!, 2^1/1!, 2^2/2!,...)   A3 = colvector(3^0/0!, 3^1/1!, 3^2/2!,...) ...     A_k= colvector(k^0/0!, k^1/1!, k^2/2!,...) Call this collection of vectors B $\hspace{24} B:=[A0,A1,A2,...]$ 3) ---------------------------------------------------------- If we consider the collection of all A_k into the matrix B, then we may also extract the common factorial denominators into a diagonal "factorial"-vector F Code:.     F     = diagonal( 0! , 1! , 2! , 3! ,... )   F^-1  = diagonal(1/0!,1/1!,1/2!,1/3!,... ) then Code:.    A0 = F^-1  * V(0)    A1 = F^-1  * V(1)    A2 = F^-1  * V(2)    A3 = F^-1  * V(3)     ...   A_k = F^-1  * V(k)     ... so we have Code:.   V(y)~ = V(x) * F^-1 * [V(0), V(1),   V(2) ,   V(3) , ....]         =               [1   , e^x , (e^x)^2, (e^x)^3,  ... ] and for obvious reasons I denote the above collection of V-vectors of subsequent parameters VZ so we have: Code:.      V(y)~ = V(x) * F^-1 * VZ            = V(e^x)~ and for more brevity I call the product F-1 * VZ $\hspace{24} B := F^{-1} * VZ$ and we have the iterable matrix-operator: $\hspace{24} V(x)\sim * B = V(e^x)\sim$ 4) ------------------------------------------------------------- Introducing another diagonal Vandermonde-vector, with powers of logs of a parameter s parametrizes this for s and I call the resulting matrix Bs: $\hspace{24} diag(V(\log{(s)})) * B = B_s$ and $\hspace{24} V(x)\sim * B_s = V(s^x)\sim$ the iterable matrix-expression for obtaining $\hspace{24} x,s^x,s^s^{\small x},s^s^{^{\small{s^x}}},\dots$ I call the matrix B a "constant", since it is independent of the parameter s. ==================================================================== The rationale in short: . We need a powerseries-vector in x and a exponential-series-vector . to get one scalar result for s^x. . . But to make it iterable the result should again be not a scalar . alone, but again a full vector of powers of s^x, we need not only . one exponential-series-vector, but one for each power, thus a . full matrix of such vectors. . . That matrix is B (resp. Bs) in the tetration case. ------------------- Once having introduced a matrix as an operator for tetration (the analoguous is valid for other iterated operations and functions as well, just take the appropriate matrices), we are in a situation to discuss iterations in terms of powers of Bs, and if Bs has an accessible eigensystem, we are also in the situation to define the fractional iteration by fractional powers and even complex powers of Bs, thus the continuous version of tetration. In my first heuristic approach I used the matrix-logarithm of Bs instead, and defined arbitrary powers by $\hspace{24} B_s ^y = \exp \left( y * \log(B_s) \right)$ I found that this provides also numerical stable approximates (with the same results, of course) but didn't investigate this path deeper, since I thought, that the eigensystem-decomposition is a more general or fundamental approach. ------------------- Hope I made it readable/understandable... 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 336 09/09/2019, 10:55 PM Last Post: tommy1729 A fundamental flaw of an operator who's super operator is addition JmsNxn 4 6,926 06/23/2019, 08:19 PM Last Post: Chenjesu Tetration and Sign operator tetration101 0 425 05/15/2019, 07:55 PM Last Post: tetration101 2 fixpoints , 1 period --> method of iteration series tommy1729 0 1,481 12/21/2016, 01:27 PM Last Post: tommy1729 Tommy's matrix method for superlogarithm. tommy1729 0 1,676 05/07/2016, 12:28 PM Last Post: tommy1729 [split] Understanding Kneser Riemann method andydude 7 7,996 01/13/2016, 10:58 PM Last Post: sheldonison [2015] New zeration and matrix log ? tommy1729 1 3,180 03/24/2015, 07:07 AM Last Post: marraco Kouznetsov-Tommy-Cauchy method tommy1729 0 2,024 02/18/2015, 07:05 PM Last Post: tommy1729 Problem with cauchy method ? tommy1729 0 1,866 02/16/2015, 01:51 AM Last Post: tommy1729 Regular iteration using matrix-Jordan-form Gottfried 7 8,643 09/29/2014, 11:39 PM Last Post: Gottfried

Users browsing this thread: 1 Guest(s)