Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
Matrix Operator Method
#28
One encouraging result for fractional iterates with non-real base

Well , perhaps I was too much scared because of the difficult approximations of the matrix-operator-method for the cases when not e^-e < b < eta.

With a careful computation of the halfiterate of the complex base "I" I got
Code:
. y1= {I,1}^^0.5                ~ 1.16729812784 + 0.735996102206*I

  y2= {I,y1}^^0.5 = {I,1}^^1 ~  0.000150635188062 + 1.00000615687*I
which is near the expected result, without need to change my hypothesis.

Remember my hypothetical formula for continuous tetration

y = {b,x}^^h

implemented by

V(y)~ = V(x)~ * (dV(log(b))*B)^h = V(x)~ * Bb^h
and
y = V(y)[1]

To arrive at the desired result, Bb^h must be constructed by the analytical description:
Let
W^-1 * D * W = Bb
be the eigendecomposition of Bb and

W^-1 * D^h * W = Bb^h

the h'th power of Bb
This can simply approximated with any eigensystem-solver, when fed with the matrix Bb. Exponentiate the eigenvalues and recompute...
But the result will then be only heuristical and we don't know the degree of approximation.

-------------

Following my hypothese about the further structure of the eigen-matrices, we can compute this structure analytically.

I assumed:
W^-1 = dV(1/t)*P^-1 ~ * X^-1
and
W = X * P~ * dV(t)

P is the known lower triangular pascal-matrix, X is a lower triangular matrix depending on t and u, and dV(t) is a diagonalmatrix containing the consecutive powers of t. The eigenvalues in D are the consecutive powers of u.

Here t and u are dependend on the base-parameter b, such that
t^(1/t) = b
u = log(t)

With my fixpoint-tracer I can first find a solution for t and b (at least) for some values b outside the range e^-e < b <e^(1/e), and it seems even for complex values of b (I have to make this sure for the general case of complex b)

The entries of X can be composed by finite polynomials in t and u and since X is triangular its inverse as well as X^-1 * D^h * X can be computed exactly/symbolically without need of approximation.

As I have already indicated in a previous post, the direct computation of W^-1 is *not* possible without approximation, since P^-1 ~ is not rowfinite, but the order of computation of the whole matrix expression can be optimized to have only one situation, where we need approximative summation, giving inexact values. This can be the last step.

The whole formula in the analytical decomposition is

V(y)~ = V(x)~ * Bb^h
V(y)~ = V(x)~ * W^-1 * D^h * W
= V(x)~ * (dV(1/t)*P^-1 ~ * X^-1 ) * D^h *(X * P~ * dV(t)

Exploiting associativity we implement this by

= (V(x)~ * dV(1/t)*P^-1 ~ )* X^-1 * D^h * X * (P~ * dV(t)

V(y/t-1)~ = V(x/t-1)~ * X^-1 * D^h * X

Denote M(t,h) = X^-1 * D^h * X
then

V(y/t-1)~ = V(x/t-1)~ * M(t,h)

M(t,h) is computed with exact terms, and only its second column is needed for the final computation. Call this terms m_r using r for the row-index.

Let (x/t-1) = z then

y/t-1 = sum(r=0..inf) m_r * z^r

and finally

y = ((sum(r=0..inf) m_r * z^r ) + 1)*t

The terms of the series m_r * z^r still do not converge over the short run of 32 terms, which I have available. But they are far better summable by Euler-summation than the terms, which occur, if I applied the naive form of summation using the second column of the precomputed Bb^h.

Possibly the problems, which remained for the matrix-method, can be reduced to mere acceleration of convergence (or regular summation of alternating divergent series), and are then just object of technically improvements of the numerical methods. That would be a very nice outcome, because we had then a firm base, from where only optimizations are required...

(This is worth a good glass of wine tonight :-) )

Gottfried
Gottfried Helms, Kassel
Reply


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

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



Users browsing this thread: 1 Guest(s)