Paper
1 September 2009 Ulam's method to estimate invariant measures and Lyapunov exponents for one-dimensional discretely randomized photonic structures
Author Affiliations +
Abstract
In the one-dimensional optical analog to Anderson localization, a periodically layered medium has one or more parameters randomly disordered. Such a medium can be modeled by an infinite product of 2x2 random transfer matrices with the upper Lyapunov exponent of the matrix product identified as the localization factor (inverse localization length). Furstenberg's integral formula for the Lyapunov exponent requires integration with respect to both the probability measure of the random matrices and the invariant probability measure of the direction of the vector propagated by the random matrix product. This invariant measure is difficult to find analytically, so one of several numerical techniques must be used in its calculation. Here, we focus on one of those techniques, Ulam's method, which sets up a sparse matrix of the probabilities that an entire interval of possible directions will be transferred to some other interval of directions. The left eigenvector of this sparse matrix forms the estimated invariant measure. While Ulam's method is shown to produce results as accurate as others, it suffers from long computation times. The Ulam method, along with other approaches, is demonstrated on a random Fibonacci sequence having a known answer, and on a quarter-wave stack model with discrete disorder in layer thickness.
© (2009) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.
Glen J. Kissel "Ulam's method to estimate invariant measures and Lyapunov exponents for one-dimensional discretely randomized photonic structures", Proc. SPIE 7392, Metamaterials: Fundamentals and Applications II, 73920S (1 September 2009); https://doi.org/10.1117/12.826439
Lens.org Logo
CITATIONS
Cited by 2 scholarly publications.
Advertisement
Advertisement
RIGHTS & PERMISSIONS
Get copyright permission  Get copyright permission on Copyright Marketplace
KEYWORDS
Matrices

Monte Carlo methods

Binary data

Analog electronics

Detection and tracking algorithms

MATLAB

Numerical integration

Back to Top