Hidden Positivity and a New Approach to Numerical Computation of Hausdorff Dimension: Higher Order Methods
share
In [14], the authors developed a new approach to the computation of the Hausdorff dimension of the invariant set of an iterated function system or IFS. In this paper, we extend this approach to incorporate high order approximation methods. We again rely on the fact that we can associate to the IFS a parametrized family of positive, linear, Perron-Frobenius operators L_s, an idea known in varying degrees of generality for many years. Although L_s is not compact in the setting we consider, it possesses a strictly positive C^m eigenfunction v_s with eigenvalue R(L_s) for arbitrary m and all other points z in the spectrum of L_s satisfy |z| ≤ b for some constant b < R(L_s). Under appropriate assumptions on the IFS, the Hausdorff dimension of the invariant set of the IFS is the value s=s_* for which R(L_s) =1. This eigenvalue problem is then approximated by a collocation method at the extended Chebyshev points of each subinterval using continuous piecewise polynomials of arbitrary degree r. Using an extension of the Perron theory of positive matrices to matrices that map a cone K to its interior and explicit a priori bounds on the derivatives of the strictly positive eigenfunction v_s, we give rigorous upper and lower bounds for the Hausdorff dimension s_*, and these bounds converge rapidly to s_* as the mesh size decreases and/or the polynomial degree increases.
READ FULL TEXT