A Piecewise Constant Method for Frobenius-Perron Operators via Delta Function Approximations
Let S : [0, 1] -> [0, 1] be a nonsingular transformation such that the corresponding Frobenius-Perron operator P-S : L-1(0, 1) -> L-1(0, 1) has a stationary density f*. We develop a piecewise constant method for the numerical computation of f*, based on the approximation of Dirac's delta function via pulse functions. We show that the numerical scheme out of this new approach is exactly the classic Ulam's method. Numerical results are given for several one dimensional test mappings. (C) 2012 Elsevier Inc. All rights reserved.
Applied Mathematics and Computation
(2012). A Piecewise Constant Method for Frobenius-Perron Operators via Delta Function Approximations. Applied Mathematics and Computation, 219(3), 1047-1052.
Available at: http://aquila.usm.edu/fac_pubs/7594