基于混沌理论的蛋白质序列特性的研究

蛋白质序列特性的研究对于蛋白质的结构及功能具有重要意义。该文为了研究蛋白质序列是否具有混沌行为,先将蛋白质序列通过氨基酸电子离子相互作用势(electron interaction potential,EIIP)转化为时间序列,再根据混沌理论对其进行相空间重构,利用去偏自相关系数,经典G-P算法确定系统的时间延迟t和嵌入维数m,系统的最大Lyapunov指数则用改进的最大Lyapunov指数计算方法计算,其结果绝大多数为正,从而确认了蛋白质时间序列的混沌行为,并对特例进行了说明。     

具有共生作用两种群离散耦合Logistic模型混沌控制

针对具有共生作用的离散耦合Logistic模型,首先采用Lyapunov指数方法验证了混沌现象的存在.然后详细地分析了系统随参数变化的分岔图,发现了系统中存在更复杂的现象.最后应用混沌跟踪控制方法控制系统的混沌现象,使得种群稳定到正不动点轨道上,消除了种群中存在的混沌现象.仿真结果验证了控制方法的有效性.     

Research on cascading high-dimensional isomorphic chaotic maps

In order to overcome the security weakness of the discrete chaotic sequence caused by small Lyapunov exponent and keyspace, a general chaotic construction method by cascading multiple high-dimensional isomorphic maps is presented in this paper. Compared with the original map, the parameter space of the resulting chaotic map is enlarged many times. Moreover, the cascaded system has larger chaotic domain and bigger Lyapunov exponents with proper parameters. In order to evaluate the effectiveness of the presented method, the generalized 3-D Hénon map is utilized as an example to analyze the dynamical behaviors under various cascade modes. Diverse maps are obtained by cascading 3-D Hénon maps with different parameters or different permutations. It is worth noting that some new dynamical behaviors, such as coexisting attractors and hyperchaotic attractors are also discovered in cascaded systems. Finally, an application of image encryption is delivered to demonstrate the excellent performance of the obtained chaotic sequences.     

Can noise induce chaos?

An important component of the mathematical definition of chaos is sensitivity to initial conditions. Sensitivity to initial conditions is usually measured in a deterministic model by the dominant Lyapunov exponent (LE), with chaos indicated by a positive LE. The sensitivity measure has been extended to stochastic models; however, it is possible for the stochastic Lyapunov exponent (SLE) to be positive when the LE of the underlying deterministic model is negative, and vice versa. This occurs because the LE is a long-term average over the deterministic attractor while the SLE is the long-term average over the stationary probability distribution. The property of sensitivity to initial conditions, uniquely associated with chaotic dynamics in deterministic systems, is widespread in stochastic systems because of time spent near repelling invariant sets (such as unstable equilibria and unstable cycles). Such sensitivity is due to a mechanism fundamentally different from deterministic chaos. Positive SLE's should therefore not be viewed as a hallmark of chaos. We develop examples of ecological population models in which contradictory LE and SLE values lead to confusion about whether or not the population fluctuations are primarily the result of chaotic dynamics. We suggest that "chaos" should retain its deterministic definition in light of the origins and spirit of the topic in ecology. While a stochastic system cannot then strictly be chaotic, chaotic dynamics can be revealed in stochastic systems through the strong influence of underlying deterministic chaotic invariant sets.     

Birhythmicity, trirhythmicity and chaos in bursting calcium oscillations

We have analyzed various types of complex calcium oscillations. The oscillations are explained with a model based on calcium-induced calcium release (CICR). In addition to the endoplasmic reticulum as the main intracellular Ca2+ store, mitochondrial and cytosolic Ca2+ binding proteins are also taken into account. This model was previously proposed for the study of the physiological role of mitochondria and the cytosolic proteins in gene rating complex Ca2+ oscillations [1]. Here, we investigated the occurrence of different types of Ca2+ oscillations obtained by the model, i.e. simple oscillations, bursting, and chaos. In a bifurcation diagram, we have shown that all these various modes of oscillatory behavior are obtained by a change of only one model parameter, which corresponds to the physiological variability of an agonist. Bursting oscillations were studied in more detail because they express birhythmicity, trirhythmicity and chaotic behavior. Two different routes to chaos are observed in the model: in addition to the usual period doubling cascade, we also show intermittency. For the characterization of the chaotic behavior, we made use of return maps and Lyapunov exponents. The potential biological role of chaos in intracellular signaling is discussed.     

Chaotic attractor in two-prey one-predator system originates from interplay of limit cycles

We investigate the appearance of chaos in a microbial 3-species model motivated by a potentially chaotic real world system (as characterized by positive Lyapunov exponents (Becks et al., Nature 435, 2005). This is the first quantitative model that simulates characteristic population dynamics in the system. A striking feature of the experiment was three consecutive regimes of limit cycles, chaotic dynamics and a fixed point. Our model reproduces this pattern. Numerical simulations of the system reveal the presence of a chaotic attractor in the intermediate parameter window between two regimes of periodic coexistence (stable limit cycles). In particular, this intermediate structure can be explained by competition between the two distinct periodic dynamics. It provides the basis for stable coexistence of all three species: environmental perturbations may result in huge fluctuations in species abundances, however, the system at large tolerates those perturbations in the sense that the population abundances quickly fall back onto the chaotic attractor manifold and the system remains. This mechanism explains how chaos helps the system to persist and stabilize against migration. In discrete populations, fluctuations can push the system towards extinction of one or more species. The chaotic attractor protects the system and extinction times scale exponentially with system size in the same way as with limit cycles or in a stable situation.     

Long food chains are in general chaotic

The question whether chaos exists in nature is much debated. In this paper we prove that chaotic parameter regions exist generically in food chains of length greater than three. While nonchaotic dynamics is also possible, the presence of chaotic parameter regions indicates that chaotic dynamics is likely. We show that the chaotic regions survive even at high exponents of closure. Our results have been obtained using a general food chain model that describes a large class of different food chains. The existence of chaos in models of such generality can be deduced from the presence of certain bifurcations of higher codimension.     

Parameter-dependent transitions and the optimal control of dynamical diseases

Many theoretical studies in biological and physical sciences consider the dynamical behavior of ann-dimensional ordinary differential equation that contains a large number of independent parameters. A frequently asked question is, are there permissible parameter sets that result in periodic or chaotic behavior? The large number of distinct parameters often limits the feasibility of trial and error calculations. The large dimension and nonlinearity of the system make application of analytic methods at best difficult and at worst effectively impossible. It is shown here that a computational search for parameter-dependent transitions of attractor topology can be effected by constrained optimization of quantitative measures of dynamical behavior (Hurwitz polynomials, Floquet coefficients, Lyapunov exponents and correlation dimension). As an example, we examine a three-dimensional nonlinear ordinary differential equation containing seven parameters that was constructed by Goldbeter and Segel to model periodic synthesis of cyclic AMP inDictyostelium. A search for bifurcations to periodic solutions is made by minimizing Hurwitz coefficients subject to parameter constraints. By comparing four optimization algorithms, the defects and advantages of the procedure are identified. It is also argued that it may be possible to use this characterization of dynamics to construct optimal responses to dynamical diseases (those disorders that result from parameter-dependent bifurcations in physiological control systems).     

Variability of bursting patterns in a neuron model in the presence of noise

Spiking and bursting patterns of neurons are characterized by a high degree of variability. A single neuron can demonstrate endogenously various bursting patterns, changing in response to external disturbances due to synapses, or to intrinsic factors such as channel noise. We argue that in a model of the leech heart interneuron existing variations of bursting patterns are significantly enhanced by a small noise. In the absence of noise this model shows periodic bursting with fixed numbers of interspikes for most parameter values. As the parameter of activation kinetics of a slow potassium current is shifted to more hyperpolarized values of the membrane potential, the model undergoes a sequence of incremental spike adding transitions accumulating towards a periodic tonic spiking activity. Within a narrow parameter window around every spike adding transition, spike alteration of bursting is deterministically chaotic due to homoclinic bifurcations of a saddle periodic orbit. We have found that near these transitions the interneuron model becomes extremely sensitive to small random perturbations that cause a wide expansion and overlapping of the chaotic windows. The chaotic behavior is characterized by positive values of the largest Lyapunov exponent, and of the Shannon entropy of probability distribution of spike numbers per burst. The windows of chaotic dynamics resemble the Arnold tongues being plotted in the parameter plane, where the noise intensity serves as a second control parameter. We determine the critical noise intensities above which the interneuron model generates only irregular bursting within the overlapped windows.     

Routes to chaos in a model of a bursting neuron.

Chaotic regimens have been observed experimentally in neurons as well as in deterministic neuronal models. The R15 bursting cell in the abdominal ganglion of Aplysia has been the subject of extensive mathematical modeling. Previously, the model of Plant and Kim has been shown to exhibit both bursting and beating modes of electrical activity. In this report, we demonstrate (a) that a chaotic regime exists between the bursting and beating modes of the model, and (b) that the model approaches chaos from both modes by a period doubling cascade. The bifurcation parameter employed is the external stimulus current. In addition to the period doubling observed in the model-generated trajectories, a period three "window" was observed, power spectra that demonstrate the approaches to chaos were generated, and the Lyaponov exponents and the fractal dimension of the chaotic attractors were calculated. Chaotic regimes have been observed in several similar models, which suggests that they are a general characteristic of cells that exhibit both bursting and beating modes.     

When can noise induce chaos and why does it matter: a critique

Noise‐induced chaos illustrates how small amounts of exogenous noise can have disproportionate qualitative impacts on the long term dynamics of a nonlinear system. This property is particularly clear in chaotic systems but is also important for the majority of ecological systems which are nonchaotic, and has direct implications for analyzing ecological time series and testing models against field data. Dennis et al. point out that a definition of chaos which we advocated allows a noise‐dominated system to be classified as chaotic when its Lyapunov exponent λ is positive, which misses what is really going on. As a solution, they propose to eliminate the concept of noise‐induced chaos: chaos “should retain its strictly deterministic definition”, hence “ecological populations cannot be strictly chaotic”. Instead, they suggest that ecologists ask whether ecological systems are strongly influenced by “underlying skeletons with chaotic dynamics or whatever other dynamics”– the skeleton being the hypothetical system that would result if all external and internal noise sources were eliminated. We agree with Dennis et al. about the problem – noise‐dominated systems should not be called chaotic – but not the solution. Even when an estimated skeleton predicts a system's short term dynamics with extremely high accuracy, the skeleton's long term dynamics and attractor may be very different from those of the actual noisy system. Using theoretical models and empirical data on microtine rodent cycles and laboratory populations of Tribolium, we illustrate how data analyses focusing on attributes of the skeleton and its attractor – such as the “deterministic Lyapunov exponent”λ₀ that Dennis et al. have used as their primary indicator of chaos – will frequently give misleading results. In contrast, quantitative measures of the actual noisy system, such as λ, provide useful information for characterizing observed dynamics and for testing proposed mechanistic explanations.     

Exciting chaos with noise: unexpected dynamics in epidemic outbreaks

 In this paper, we identify a mechanism for chaos in the presence of noise. In a study of the SEIR model, which predicts epidemic outbreaks in childhood diseases, we show how chaotic dynamics can be attained by adding stochastic perturbations at parameters where chaos does not exist apriori. Data recordings of epidemics in childhood diseases are still argued as deterministic chaos. There also exists noise due to uncertainties in the contact parameters between those who are susceptible and those who are infected, as well as random fluctuations in the population. Although chaos has been found in deterministic models, it only occurs in parameter regions that require a very large population base or other large seasonal forcing. Our work identifies the mechanism whereby chaos can be induced by noise for realistic parameter regions of the deterministic model where it does not naturally occur. Received: 13 October 2000 / Revised version: 15 May 2001 / Published online: 7 December 2001     

Nonlinear dynamics of two-dimensional cardiac action potential duration mapping model with memory

The aim of this work is the analysis of the nonlinear dynamics of two-dimensional mapping model of cardiac action potential duration (2D-map APD) with memory derived from one dimensional map (1D-map). Action potential duration (APD) restitution, which relates APD to the preceding diastolic interval (DI), is a useful tool for predicting cardiac arrhythmias. For a constant rate of stimulation the short action potential during alternans is followed by a longer DI and inversely. It has been suggested that these differences in DI are responsible for the occurrence and maintenance of APD alternans. We focus our attention on the observed bifurcations produced by a change in the stimulation period and a fixed value of a particular parameter in the model. This parameter provides new information about the dynamics of the APD with memory, such as the occurrence of bistabilities not previously described in the literature, as well as the fact that synchronization rhythms occur in different ways and in a new fashion as the stimulation frequency increases. Moreover, we show that this model is flexible enough as to accurately reflect the chaotic dynamics properties of the APD: we have highlighted the fractal structure of the strange attractor of the 2D-map APD, and we have characterized chaos by tools such as the calculation of the Lyapunov exponents, the fractal dimension and the Kolmogorov entropy, with the next objective of refining the study of the nonlinear dynamics of the duration of the action potential and to apply methods of controlling chaos.

     

Explosive route to chaos through a fractal torus in a generalized Lotka-Volterra model

The behavior of a model that generalizes the Lotka-Volterra problem into three dimensions is presented. The results show the analytic derivation of stability diagrams that describe the system's qualitative features. In particular, we show that for a certain value of the bifurcation parameter the system instantly jumps out of a steady state solution into a chaotic solution that portrays a fractal torus in the three-dimensional phase space. This scenario, is referred to as the explosive route to chaos and is attributed to the non-transversal saddle connection type bifurcation. The stability diagrams also present a region in which the Hopf type bifurcation leads to periodic and chaotic solutions. In addition, the bifurcation diagrams reveal a qualitative similarity to the data obtained in the Texas and Bordeaux experiments on the Belousov-Zhabotinskii chemical reaction. The paper is concluded by showing that the model can be useful for representing dynamics associated with biological and chemical phenomena.     

Analysis of cedar pollen time series: no evidence of low-dimensional chaotic behavior

Much of the current interest in pollen time series analysis is motivated by the possibility that pollen series arise from low-dimensional chaotic systems. If this is the case, short-range prediction using nonlinear modeling is justified and would produce high-quality forecasts that could be useful in providing pollen alerts to allergy sufferers. To date, contradictory reports about the characterization of the dynamics of pollen series can be found in the literature. Pollen series have been alternatively described as featuring and not featuring deterministic chaotic behavior. We showed that the choice of test for detection of deterministic chaos in pollen series is difficult because pollen series exhibit power spectra. This is a characteristic that is also produced by colored noise series, which mimic deterministic chaos in most tests. We proposed to apply the Ikeguchi–Aihara test to properly detect the presence of deterministic chaos in pollen series. We examined the dynamics of cedar (Cryptomeria japonica) hourly pollen series by means of the Ikeguchi–Aihara test and concluded that these pollen series cannot be described as low-dimensional deterministic chaos. Therefore, the application of low-dimensional chaotic deterministic models to the prediction of short-range pollen concentration will not result in high-accuracy pollen forecasts even though these models may provide useful forecasts for certain applications. We believe that our conclusion can be generalized to pollen series from other wind-pollinated plant species, as wind speed, the forcing parameter of the pollen emission and transport, is best described as a nondeterministic series that originates in the high dimensionality of the atmosphere.     

Nonlinear population dynamics: complication in the age structure influences transition-to-chaos scenarios

The work continues a series of studies on the evolution of a natural population of explicitly seasonal organisms. Model analyses have revealed relationships between the duration of ontogenesis and the pattern of temporal dynamics in size of an isolated population (i.e., the structure and dimensionality of the chaotic attractors). For nonlinear models of age-structured population dynamics (under long-lasting ontogenesis), increase in the reproductive potential is shown to result in the chaotic attractors whose structure and dimensionality changes in response to variations in the model parameters. When the ontogenesis becomes longer and more complicated, it does not, "on the average", augment the level of chaos in the attractors observed. There are wide enough regions in the space of the birth and death parameter values that provide for windows in the chaotic dynamics where the total or partial regularization occurs.     

Subharmonic resonance and chaos in forced excitable systems

 Forced excitable systems arise in a number of biological and physiological applications and have been studied analytically and computationally by numerous authors. Existence and stability of harmonic and subharmonic solutions of a forced piecewise-linear Fitzhugh-Nagumo-like system were studied in Othmer ad Watanabe (1994) and in Xie et al. (1996). The results of those papers were for small and moderate amplitude forcing. In this paper we study the existence of subharmonic solutions of this system under large-amplitude forcing. As in the case of intermediate-amplitude forcing, bistability between 1 : 1 and 2 : 1 solutions is possible for some parameters. In the case of large-amplitude forcing, bistability between 2 : 2 and 2 : 1 solutions, which does not occur in the case of intermediate-amplitude forcing, is also possible for some parameters. We identify several new canonical return maps for a singular system, and we show that chaotic dynamics can occur in some regions of parameter space. We also prove that there is a direct transition from 2 : 2 phase-locking to chaos after the first period-doubling bifurcation, rather than via the infinite sequence of period doublings seen in a smooth quadratic interval map. Coexistence of chaotic dynamics and stable phase-locking can also occur. Received: 6 July 1998 / Revised version: 2 October 1998     

Complex Generalized Synchronization and Parameter Identification of Nonidentical Nonlinear Complex Systems

In this paper, generalized synchronization (GS) is extended from real space to complex space, resulting in a new synchronization scheme, complex generalized synchronization (CGS). Based on Lyapunov stability theory, an adaptive controller and parameter update laws are designed to realize CGS and parameter identification of two nonidentical chaotic (hyperchaotic) complex systems with respect to a given complex map vector. This scheme is applied to synchronize a memristor-based hyperchaotic complex Lü system and a memristor-based chaotic complex Lorenz system, a chaotic complex Chen system and a memristor-based chaotic complex Lorenz system, as well as a memristor-based hyperchaotic complex Lü system and a chaotic complex Lü system with fully unknown parameters. The corresponding numerical simulations illustrate the feasibility and effectiveness of the proposed scheme.