Using Special-Purpose Computing to Examine Chaotic Behavior in Nonlinear Mappings

Studying chaotic behavior in nonlinear systems requires numerous computations in order to simulate the behavior of such systems. The Standard Map Machine was designed and implemented as a special computer for performing these intensive computations with high-speed and high-precision. Its impressive performance is due to its simple architecture specialized to the numerical computations required of nonlinear systems. This report discusses the design and implementation of the Standard Map Machine and its use in the study of nonlinear mappings; in particular, the study of the standard map.

One-dimensional Fermi accelerator model with moving wall described by a nonlinear van der Pol oscillator

A modification of the one-dimensional Fermi accelerator model is considered in this work. The dynamics of a classical particle of mass m, confined to bounce elastically between two rigid walls where one is described by a nonlinear van der Pol type oscillator while the other one is fixed, working as a reinjection mechanism of the particle for a next collision, is carefully made by the use of a two-dimensional nonlinear mapping. Two cases are considered: (i) the situation where the particle has mass negligible as compared to the mass of the moving wall and does not affect the motion of it; and (ii) the case where collisions of the particle do affect the movement of the moving wall. For case (i) the phase space is of mixed type leading us to observe a scaling of the average velocity as a function of the parameter (χ) controlling the nonlinearity of the moving wall. For large χ, a diffusion on the velocity is observed leading to the conclusion that Fermi acceleration is taking place. On the other hand, for case (ii), the motion of the moving wall is affected by collisions with the particle. However, due to the properties of the van der Pol oscillator, the moving wall relaxes again to a limit cycle. Such kind of motion absorbs part of the energy of the particle leading to a suppression of the unlimited energy gain as observed in case (i). The phase space shows a set of attractors of different periods whose basin of attraction has a complicated organization. © 2013 American Physical Society.

Finding invariant tori in the problem of a periodically corrugated waveguide

Some dynamic properties for a light ray suffering specular reflections inside a periodically corrugated waveguide are studied. The dynamics of the model is described in terms of a two dimensional nonlinear area preserving map. We show that the phase space is mixed in the sense that there are KAM islands surrounded by a large chaotic sea that is confined by two invariant spanning curves. We have used a connection with the Standard Mapping near a transition from local to global chaos and found the position of these two invariant spanning curves limiting the size of the chaotic sea as function of the control parameter.

On the statistical and transport properties of a non-dissipative Fermi-Ulam model

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Detection and Classification of Whale Acoustic Signals

This dissertation focuses on two vital challenges in relation to whale acoustic signals: detection and classification.

In detection, we evaluated the influence of the uncertain ocean environment on the spectrogram-based detector, and derived the likelihood ratio of the proposed Short Time Fourier Transform detector. Experimental results showed that the proposed detector outperforms detectors based on the spectrogram. The proposed detector is more sensitive to environmental changes because it includes phase information.

In classification, our focus is on finding a robust and sparse representation of whale vocalizations. Because whale vocalizations can be modeled as polynomial phase signals, we can represent the whale calls by their polynomial phase coefficients. In this dissertation, we used the Weyl transform to capture chirp rate information, and used a two dimensional feature set to represent whale vocalizations globally. Experimental results showed that our Weyl feature set outperforms chirplet coefficients and MFCC (Mel Frequency Cepstral Coefficients) when applied to our collected data.

Since whale vocalizations can be represented by polynomial phase coefficients, it is plausible that the signals lie on a manifold parameterized by these coefficients. We also studied the intrinsic structure of high dimensional whale data by exploiting its geometry. Experimental results showed that nonlinear mappings such as Laplacian Eigenmap and ISOMAP outperform linear mappings such as PCA and MDS, suggesting that the whale acoustic data is nonlinear.

We also explored deep learning algorithms on whale acoustic data. We built each layer as convolutions with either a PCA filter bank (PCANet) or a DCT filter bank (DCTNet). With the DCT filter bank, each layer has different a time-frequency scale representation, and from this, one can extract different physical information. Experimental results showed that our PCANet and DCTNet achieve high classification rate on the whale vocalization data set. The word error rate of the DCTNet feature is similar to the MFSC in speech recognition tasks, suggesting that the convolutional network is able to reveal acoustic content of speech signals.

Block and parallel modelling of broad domain nonlinear continuous mapping based on NN

The necessity of the use of the block and parallel modeling of the nonlinear continuous mappings with NN is firstly expounded quantitatively. Then, a practical approach for the block and parallel modeling of the nonlinear continuous mappings with NN is proposed. Finally, an example indicating that the method raised in this paper can be realized by suitable existed software is given. The results of the experiment of the model discussed on the 3-D Mexican straw hat indicate that the block and parallel modeling based on NN is more precise and faster in computation than the direct ones and it is obviously a concrete example and the development of the large-scale general model established by Tu Xuyan.

The Specialized Mappings Architecture

A probabilistic, nonlinear supervised learning model is proposed: the Specialized Mappings Architecture (SMA). The SMA employs a set of several forward mapping functions that are estimated automatically from training data. Each specialized function maps certain domains of the input space (e.g., image features) onto the output space (e.g., articulated body parameters). The SMA can model ambiguous, one-to-many mappings that may yield multiple valid output hypotheses. Once learned, the mapping functions generate a set of output hypotheses for a given input via a statistical inference procedure. The SMA inference procedure incorporates an inverse mapping or feedback function in evaluating the likelihood of each of the hypothesis. Possible feedback functions include computer graphics rendering routines that can generate images for given hypotheses. The SMA employs a variant of the Expectation-Maximization algorithm for simultaneous learning of the specialized domains along with the mapping functions, and approximate strategies for inference. The framework is demonstrated in a computer vision system that can estimate the articulated pose parameters of a human’s body or hands, given silhouettes from a single image. The accuracy and stability of the SMA are also tested using synthetic images of human bodies and hands, where ground truth is known.

Specialized Mappings Architecture with Applications to Vision-based Estimation of Articulated Body Pose

A fundamental task of vision systems is to infer the state of the world given some form of visual observations. From a computational perspective, this often involves facing an ill-posed problem; e.g., information is lost via projection of the 3D world into a 2D image. Solution of an ill-posed problem requires additional information, usually provided as a model of the underlying process. It is important that the model be both computationally feasible as well as theoretically well-founded. In this thesis, a probabilistic, nonlinear supervised computational learning model is proposed: the Specialized Mappings Architecture (SMA). The SMA framework is demonstrated in a computer vision system that can estimate the articulated pose parameters of a human body or human hands, given images obtained via one or more uncalibrated cameras. The SMA consists of several specialized forward mapping functions that are estimated automatically from training data, and a possibly known feedback function. Each specialized function maps certain domains of the input space (e.g., image features) onto the output space (e.g., articulated body parameters). A probabilistic model for the architecture is first formalized. Solutions to key algorithmic problems are then derived: simultaneous learning of the specialized domains along with the mapping functions, as well as performing inference given inputs and a feedback function. The SMA employs a variant of the Expectation-Maximization algorithm and approximate inference. The approach allows the use of alternative conditional independence assumptions for learning and inference, which are derived from a forward model and a feedback model. Experimental validation of the proposed approach is conducted in the task of estimating articulated body pose from image silhouettes. Accuracy and stability of the SMA framework is tested using artificial data sets, as well as synthetic and real video sequences of human bodies and hands.

Semi-hyperbolic mappings in Banach spaces

The definition of semi-hyperbolic dynamical systems generated by Lipschitz continuous and not necessarily invertible mappings in Banach spaces is presented in this thesis. Like hyperbolic mappings, they involve a splitting into stable and unstable spaces, but a slight leakage from the strict invariance of the spaces is possible and the unstable subspaces are assumed to be finite dimensional. Bi-shadowing is a combination of the concepts of shadowing and inverse shadowing and is usually used to compare pseudo-trajectories calculated by a computer with the true trajectories. In this thesis, the concept of bi-shadowing in a Banach space is defined and proved for semi-hyperbolic dynamical systems generated by Lipschitz mappings. As an application to the concept of bishadowing, linear delay differential equations are shown to be bi-shadowing with respect to pseudo-trajectories generated by nonlinear small perturbations of the linear delay equation. This shows robustness of solutions of the linear delay equation with respect to small nonlinear perturbations. Complicated dynamical behaviour is often a consequence of the expansivity of a dynamical system. Semi-hyperbolic dynamical systems generated by Lipschitz mappings on a Banach space are shown to be exponentially expansive, and explicit rates of expansion are determined. The result is applied to a nonsmooth noninvertible system generated by delay differential equation. It is shown that semi-hyperbolic mappings are locally φ-contracting, where -0 is the Hausdorff measure of noncompactness, and that a linear operator is semi-hyperbolic if and only if it is φ-contracting and has no spectral values on the unit circle. The definition of φ-bi-shadowing is given and it is shown that semi-hyperbolic mappings in Banach spaces are φ-bi-shadowing with respect to locally condensing continuous comparison mappings. The result is applied to linear delay differential equations of neutral type with nonsmooth perturbations. Finally, it is shown that a small delay perturbation of an ordinary differential equation with a homoclinic trajectory is ‘chaotic’.

Saddle points and rare collisions under scaling approach in a Fermi accelerator with two nonlinear terms

Rare collisions of a classical particle bouncing between two walls are studied. The dynamics is described by a two-dimensional, nonlinear and area-preserving mapping in the variables velocity and time at the instant that the particle collides with the moving wall. The phase space is of mixed type preventing diffusion of the particle to high energy. Successive and therefore rare collisions are shown to have a histogram of frequency which is scaling invariant with respect to the control parameters. The saddle fixed points are studied and shown to be scaling invariant with respect to the control parameters too. © 2012 Elsevier B.V. All rights reserved.

Scaling invariance of the diffusion coefficient in a family of two-dimensional Hamiltonian mappings

We consider a family of two-dimensional nonlinear area-preserving mappings that generalize the Chirikov standard map and model a variety of periodically forced systems. The action variable diffuses in increments whose phase is controlled by a negative power of the action and hence effectively uncorrelated for small actions, leading to a chaotic sea in phase space. For larger values of the action the phase space is mixed and contains a family of elliptic islands centered on periodic orbits and invariant Kolmogorov-Arnold-Moser (KAM) curves. The transport of particles along the phase space is considered by starting an ensemble of particles with a very low action and letting them evolve in the phase until they reach a certain height h. For chaotic orbits below the periodic islands, the survival probability for the particles to reach h is characterized by an exponential function, well modeled by the solution of the diffusion equation. On the other hand, when h reaches the position of periodic islands, the diffusion slows markedly. We show that the diffusion coefficient is scaling invariant with respect to the control parameter of the mapping when h reaches the position of the lowest KAM island. © 2013 American Physical Society.

A rescaling of the phase space for Hamiltonian map: Applications on the Kepler map and mappings with diverging angles in the limit of vanishing action

A rescale of the phase space for a family of two-dimensional, nonlinear Hamiltonian mappings was made by using the location of the first invariant Kolmogorov-Arnold-Moser (KAM) curve. Average properties of the phase space are shown to be scaling invariant and with different scaling times. Specific values of the control parameters are used to recover the Kepler map and the mapping that describes a particle in a wave packet for the relativistic motion. The phase space observed shows a large chaotic sea surrounding periodic islands and limited by a set of invariant KAM curves whose position of the first of them depends on the control parameters. The transition from local to global chaos is used to estimate the position of the first invariant KAM curve, leading us to confirm that the chaotic sea is scaling invariant. The different scaling times are shown to be dependent on the initial conditions. The universality classes for the Kepler map and mappings with diverging angles in the limit of vanishing action are defined. © 2013 Published by Elsevier Inc. All rights reserved.

Nonlinear Normalization in Limit Theorems for Extremes

2000 Mathematics Subject Classification: 60G70, 60F05.