Discrete Dynamics in Nature and Society
Volume 2015 (2015), Article ID 658048, 5 pages
http://dx.doi.org/10.1155/2015/658048

Research Article

Controlling the Stochastic Sensitivity in Nonlinear Discrete-Time Systems with Incomplete Information

Lev Ryashko and Irina Bashkirtseva

Ural Federal University, Lenina, 51, Ekaterinburg 620000, Russia

Received 19 April 2015; Revised 13 September 2015; Accepted 16 September 2015

Academic Editor: Zhan Zhou

Copyright © 2015 Lev Ryashko and Irina Bashkirtseva. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

For stochastic nonlinear discrete-time system with incomplete information, a problem of the stabilization of equilibrium is considered. Our approach uses a regulator which synthesizes the required stochastic sensitivity. Mathematically, this problem is reduced to the solution of some quadratic matrix equations. A description of attainability sets and algorithms for regulators design is given. The general results are applied to the suppression of unwanted large-amplitude oscillations around the equilibria of the stochastically forced Verhulst model with noisy observations.

1. Introduction

Controlling of the complex systems in nature and society is a challenging and fundamental problem of the modern mathematical theory of nonlinear dynamics and engineering. Discrete dynamic models because of widely used computer-oriented technologies, attract attention of many researchers [1, 2]. Even in simple discrete models, due to nonlinearity, a variety of dynamic regimes, both regular and chaotic, is observed [3–5]. An interplay of nonlinearity and stochasticity can generate new unexpected phenomena [6–10].

A lot of nonlinear systems operate in zones of stochastic transitions from order to chaos. After the pioneering work [11], a problem of controlling chaos is extensively studied [12–14]. Most of the reported results are based on the direct numerical simulation. A detailed theoretical description of stochastic attractor is given by the stationary probabilistic density function. For discrete systems, this function is governed by Frobenius-Perron equation [15]. Unfortunately, an analytical solution of this equation is possible only in very special cases, so, a development of the asymptotic approximations is a highly relevant area of research.

For the constructive analysis of the stochastic attractors of nonlinear discrete-time dynamical systems, a stochastic sensitivity functions technique was elaborated [16]. This technique was applied to the analysis on noise-induced intermittency [17] and neuron excitability [18]. On the base of this technique, a new approach for the solution of control problems in stochastic discrete-time systems was suggested in [19]. In these studies, it was supposed that the complete information about current state of the controlled system is known. However, in many practical situations, the system data are far from complete. For example, only some coordinates of the system state are observable, and, moreover, these observations contain stochastic errors. So, the control of stochastic systems with incomplete information is an urgent research domain [20–22].

In present paper, we further develop a theory for the synthesis of the stochastic sensitivity for the equilibria in a randomly forced control discrete system with incomplete information. Mathematically, presence of noise in the observations leads to a new algebraic analysis of quadratic matrix equations. In Section 2, we introduce the stochastic sensitivity matrix as a basic probabilistic characteristics for the randomly forced equilibria. A problem of the synthesis of this matrix is considered. A important notion of the attainability is discussed here. A problem of the stochastic sensitivity matrix synthesis is reduced to the analysis of the corresponding quadratic matrix equation. Results of this theoretical analysis in the general multidimensional case are presented in a Theorem. This Theorem gives a description of attainability sets and algorithms for regulators design.

One-dimensional case is discussed in details in Section 3. In Section 4, we apply the results to the suppression of unwanted large-amplitude oscillations around the equilibria of the stochastically forced Verhulst model with noisy observations. We show that our regulator can be used for the suppression of chaos.

2. Synthesis of Stochastic Sensitivity

Consider a nonlinear controlled discrete-time stochastic system

where , , , and is a control input. Here, is an uncorrelated random sequence with parameters and , and is the identity -matrix and is a scalar parameter of noise intensity.

It is supposed that the corresponding deterministic uncontrolled system (1) (with and therein) has an equilibrium . Stability of is not assumed.

In present paper, we consider a case of incomplete information when the measurement vector is known only:

where , . Here, is an uncorrelated random sequence with parameters and , and is the identity -matrix.

In this circumstance, we consider the following regulator:

The dynamics of the closed-loop stochastic system (1) with the regulator (3) using noisy observations (2) is governed by the following system:

For the asymptotics of the deviations of solutions of system (4) from the equilibrium , the following stochastic system can be written:

where

Due to the uncorrelatedness of random terms and , the second moments matrix is governed by the equation

where , . A set of matrices that provide an exponential stability to the equilibrium of the closed deterministic system (4) (with therein) has the following form:

where is a spectral radius of the matrix . We suppose that the set is not empty.

For any , (7) has a unique stable stationary solution satisfying the equation

This matrix is called the stochastic sensitivity matrix of the equilibrium for system (4). The stochastic sensitivity matrix approximates a limit behavior of the second moments for deviations of solutions from :

So, the matrix characterizes a dispersion of the stationary distributed random states of system (4) around the equilibrium .

For any , the regulator (3) forms a corresponding stochastic equilibrium of system (4) with the stochastic sensitivity matrix which is a solution of (9).

Consider further the following inverse problem.

Problem of Stochastic Sensitivity Synthesis. Let be a set of symmetric and positive-definite -matrices. Let be some assigned matrix. The problem is to find a feedback matrix of regulator (3) such that the equality holds. Here, is a solution of (9).

In some cases, this problem can be unsolvable. Therefore, we consider an important notion of the attainability.

Definition 1. An element is said to be attainable for system (4) if the equality holds for some .

Definition 2. The set of all attainable elements,

is called the attainability set for system (4).

As it follows from (9), the attainability analysis is reduced to the study of solvability of the quadratic matrix equation:

Rewrite (12) with respect to a new unknown matrix

in the following form:

Denote . Suppose that the matrix is positive-definite . A substitution transforms (14) into the following equation:

where

A necessary condition of (15) solvability is in the nonnegative definiteness of the matrix :

Let condition (17) be fulfilled. Then, quadratic equation (15) is equivalent to the linear equation

where is an arbitrary orthogonal -matrix. It follows from (18) that the feedback matrix of the regulator (3) which synthesizes the stochastic sensitivity matrix , satisfies to the linear matrix equation

In the following theorem, we summarize our theoretical results.

Theorem 3. Let noises in system (1) and observations (2) be nonsingular (, ).
(a) If the matrix is quadratic and nonsingular () then

and, for any matrix , (19) has a solution

(b) If then

and, for any matrix , (19) has a solution

Here, is an arbitrary orthogonal -matrix, is a projective matrix, and a “+” sign means a pseudoinversion [23].

3. Controlling of One-Dimensional Stochastic System

Consider one-dimensional discrete stochastic controlled system

with noisy observations

Here, , , are scalar variables of state, output, and control input; , are uncorrelated random scalar sequences with parameters , , , and , and , , are scalar parameters of noise intensities.

It is supposed that the corresponding deterministic uncontrolled system (24) with and therein has an equilibrium . In what follows, we use the regulator

For the synthesis of the assigned scalar stochastic sensitivity of the equilibrium , we apply theoretical results presented above.

At first describe the attainability set for the considered example. The function from (17) has here the following representation:

The attainability condition (see Theorem 3 in Section 2) is equivalent to quadratic inequality

Thus, all attainable values of have to satisfy the inequality

Note that the value is a minimal value of the stochastic sensitivity that we can provide by this regulator.

Our regulator will synthesize any assign stochastic sensitivity if we will take (see Theorem 3 in Section 2) a feedback coefficient as follows:

Note that the optimal regulator synthesizing the minimal value of the stochastic sensitivity has the feedback coefficient

4. Example: Controlling Stochastic Verhulst System

Consider stochastically forced well-known Verhulst system with control and noisy observations:

The corresponding deterministic uncontrolled system (32) has a nontrivial equilibrium . This equilibrium is stable for and unstable for .

At first consider the influence of random disturbances for system (32) without control (). Under the stochastic disturbances, for random states of this system some probabilistic distribution is formed [24]. In Figure 1(a), random states of system (32), with for , calculated by direct numerical simulation, are plotted by red color.

(a)

(b)

(c)

Figure 1: Stochastic Verhulst model for , , and without control (red) and with optimal control (blue): (a) random states; (b) dispersion; (c) largest Lyapunov exponent.

As the parameter crosses the bifurcation value from the left to right, a dispersion of these random states sharply increases. Let be a mean square deviation of random states from the equilibrium . This function is plotted in Figure 1(b) by red color.

Consider now abilities of control . The aim of the control is to stabilize the equilibrium and provide a small dispersion of random states in a wide range of the parameter . For the solution of this problem, we will use the optimal regulator that minimizes a stochastic sensitivity of the equilibrium .

For the Verhulst system, , and the minimal value of the stochastic sensitivity (see (29)) that we can provide by the regulator has the following explicit parametric representation:

For the optimal regulator that synthesizes this minimal value , the feedback coefficient

Results of the control by this regulator are shown in Figure 1 by blue color. As one can see in Figure 1(a), random states of the controlled system are well localized near the equilibrium regardless of whether this equilibrium is stable or not. This optimal control provides the smaller and uniform dispersion (see Figure 1(b)).

Along with the spatial description of the stochastic Verhulst system, consider dynamics of the mutual arrangement of random trajectories in stochastic flows. For quantitative description of this dynamics, the largest Lyapunov exponent is traditionally used [25]. In Figure 1(c), plots of are shown for uncontrolled (red) and controlled (blue) Verhulst system. A change of sign from minus to plus justifies a transition of the uncontrolled Verhulst system from order to chaos.

As one can see, the optimal regulator essentially decreases values of , and provides their negativeness for the whole zone of the parameter . So, our regulator suppresses chaos and provides a structural stabilization of the system.

5. Conclusion

In this paper, a problem of the stabilization of random states for the nonlinear discrete-time stochastic system with incomplete information was studied. To solve this problem, the stochastic sensitivity technique was developed. Mathematically, the considered control problem was reduced to the analysis of the corresponding quadratic matrix equation. For this equation, we have analyzed its solvability and suggested a constructive method of the solution. A description of attainability sets and algorithms for regulators design were accumulated in Theorem 3. A constructiveness of these general theoretical results was demonstrated for the stochastically forced Verhulst model. It is shown that our method suppresses unwanted large-amplitude oscillations around the equilibria and transforms the system dynamics from chaotic to regular. It is worth noting that the elaborated technique is readily applicable to the control of higher dimensional discrete-time stochastic models.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The work was supported by Act 211 Government of the Russian Federation, Contract 02.A03.21.0006.

References

J. T. Sandefur, Discrete Dynamical Modeling, Oxford University Press, Oxford, UK, 1993.

S. Elaydi, “Is the world evolving discretely?” Advances in Applied Mathematics, vol. 31, no. 1, pp. 1–9, 2003.

R. L. Devaney, An Introduction to Chaotic Dynamical Systems, Westview Press, 2003.

H. G. Schuster, Deterministic Chaos, Physik, Weinheim, Germany, 1984.

B. Ibarz, J. M. Casado, and M. A. F. Sanjuán, “Map-based models in neuronal dynamics,” Physics Reports, vol. 501, no. 1-2, pp. 1–74, 2011.

W. Horsthemke and R. Lefever, Noise-Induced Transitions, Springer Series in Synergetics, Springer, Berlin, Germany, 1984.

K. Matsumoto and I. Tsuda, “Noise-induced order,” Journal of Statistical Physics, vol. 31, no. 1, pp. 87–106, 1983.

F. Gassmann, “Noise-induced chaos-order transitions,” Physical Review E, vol. 55, article 2215, 1997.

J. B. Gao, S. K. Hwang, and J. M. Liu, “When can noise induce chaos?” Physical Review Letters, vol. 82, no. 6, pp. 1132–1135, 1999.

N. Hritonenko, A. Rodkina, and Y. Yatsenko, “Stability analysis of stochastic Ricker population model,” Discrete Dynamics in Nature and Society, vol. 2006, Article ID 64590, 13 pages, 2006.

E. Ott, C. Grebogi, and J. A. Yorke, “Controlling chaos,” Physical Review Letters, vol. 64, no. 11, pp. 1196–1199, 1990.

A. L. Fradkov and A. Y. Pogromsky, Introduction to Control of Oscillations and Chaos, vol. 35, World Scientific, River Edge, NJ, USA, 1998.

G. Chen and X. Yu, Chaos Control: Theory and Applications, Springer, New York, NY, USA, 2003.

M. Sanjuan and C. Grebogi, Recent Progress in Controlling Chaos, World Scientific, Singapore, 2010.

A. Lasota and M. C. Mackey, Chaos, Fractals, and Noise: Stochastic Aspects of Dynamics, Springer-Verlag, New York, NY, USA, 1994.

I. Bashkirtseva, L. Ryashko, and I. Tsvetkov, “Sensitivity analysis of stochastic equilibria and cycles for the discrete dynamic systems,” Dynamics of Continuous, Discrete and Impulsive Systems, Series A: Mathematical Analysis, vol. 17, no. 4, pp. 501–515, 2010.

I. Bashkirtseva and L. Ryashko, “Stochastic sensitivity analysis of noise-induced intermittency and transition to chaos in one—dimensional discrete-time systems,” Physica A, vol. 392, no. 2, pp. 295–306, 2013.

I. Bashkirtseva, “Stochastic phenomena in one-dimensional Rulkov model of neuronal dynamics,” Discrete Dynamics in Nature and Society, vol. 2015, Article ID 495417, 7 pages, 2015.

I. Bashkirtseva and L. Ryashko, “Control of equilibria for nonlinear stochastic discrete-time systems,” IEEE Transactions on Automatic Control, vol. 56, no. 9, pp. 2162–2166, 2011.

K. J. Astrom, Introduction to the Stochastic Control Theory, Academic Press, New York, NY, USA, 1970.

B. Shen, Z. Wang, and H. Shu, Nonlinear Stochastic Systems with Incomplete Information, Filtering and Control, Springer, 2013.

H. Dong, Z. Wang, and H. Gao, Filtering, Control and Fault Detection with Randomly Occurring Incomplete Information, John Wiley & Sons, Ltd, New York, NY, USA, 2013.

A. Albert, Regression and the Moore-Penrose Pseudoinverse, Academic Press, New York, NY, USA, 1972.

J. P. Crutchfield, J. D. Farmer, and B. A. Huberman, “Fluctuations and simple chaotic dynamics,” Physics Reports A: Review Section of Physics Letters, vol. 92, no. 2, pp. 45–82, 1982.

O. Martin, “Lyapunov exponents of stochastic dynamical systems,” Journal of Statistical Physics, vol. 41, no. 1-2, pp. 249–261, 1985.