Abstract

We solve the problem of regulating the output of the Pan system (2010), which is one of the recently discovered three-dimensional chaotic attractors. Pan system has many interesting complex dynamical behaviours, and it has potential applications in secure communication. In this paper, we construct explicit state feedback control laws for regulating the output of the Pan system so as to track constant reference signals. The state feedback control laws are derived using the regulator equations of Byrnes and Isidori (1990). The simulation results are provided to illustrate the effectiveness of the regulation schemes derived for the output regulation of the Pan system.

1. Introduction

Regulating the output of nonlinear control systems is one of the central problems in nonlinear control theory. Essentially, the output regulation problem is the problem of controlling a fixed linear or nonlinear plant so that the output of the plant tracks the reference signals produced by some external input generator (the exosystem). For linear control systems, the output regulation problem was solved by Francis and Wonham [1]. For nonlinear control systems, the output regulation problem was solved by Byrnes and Isidori [2], who generalized the internal model principle obtained by Francis and Wonham and proved their results using Centre Manifold Theory for flows [3]. Byrnes and Isidori made an important assumption in their work [2] that the exosystem is a neutrally stable system (Lyapunov stable system in both forward and backward time). This class of exosystem signals includes the important special cases of constant reference signals as well as periodic reference signals.

The output regulation problem for linear and nonlinear control systems has been the focus of many important studies in the last two decades [4–14]. In [4], Mahmoud and Khalil derived results for the problem of asymptotic regulation of minimum phase nonlinear systems using output feedback. In [5], Fridman derived results on the output regulation problem for nonlinear control systems with delay. In [6, 7], Chen and Huang derived results on the robust output regulation problem for output feedback systems with nonlinear exosystems. In [8], Liu and Huang derived results on the global robust output regulation problem for lower triangular nonlinear systems with unknown control direction. In [9], Immonen derived results on the practical output regulation for bounded linear infinite-dimensional state space systems. In [10], Pavlov et al. derived results on the global nonlinear output regulation using convergence-based controller design. In [11], Xi and Ding studied the global adaptive output regulation of a class of nonlinear systems with nonlinear exosystems. In [12–14], Serrani et al. derived results on the semiglobal and global output regulation problem for minimum phase nonlinear systems.

In this paper, we solve the output regulation problem for the chaotic Pan system (2010). The Pan system [15] is an important model of three-dimensional chaotic system discovered by Pan et al. From the control engineering point of view, the Pan system provides an interesting framework for advanced control techniques since it is a very complex system than the Lorenz’s system. In this paper, we use the regulator equations for nonlinear control systems [2] to derive state feedback control laws for regulating the output of the Pan system so as to track constant reference input signals (set-point signals).

This paper is organized as follows. In Section 2, we present a review of the output regulation problem for nonlinear control systems and its solution using the regulator equations derived by Byrnes and Isidori [2]. In Section 3, we detail our solution of the output regulation problem for the Pan system. In Section 4, we discuss the simulation results illustrating the effectiveness of the state feedback control laws derived in Section 3. In Section 5, we summarize the main results obtained in this paper.

2. Review of the Output Regulation for Nonlinear Control Systems

In this section, we consider a multivariable nonlinear control system given byΜ‡π‘₯=𝑓(π‘₯)+𝑔(π‘₯)𝑒+𝑝(π‘₯)πœ”,(2.1a)Μ‡πœ”=𝑠(πœ”),(2.1b)𝑒=β„Ž(π‘₯)βˆ’π‘ž(πœ”).(2.2) Here, the differential equation (2.1a)) describes the plant dynamics with state π‘₯ defined in a neighbourhood 𝑋 of the origin of 𝑅𝑛, and the input 𝑒 takes values in π‘…π‘š subject to the effect of a disturbance represented by the vector field 𝑝(π‘₯)πœ”. The differential equation ((2.1b) describes the autonomous system, known as the exosystem, defined in a neighbourhood π‘Šof the origin of π‘…π‘˜, which models the class of disturbance and reference signals taken into consideration.

We also assume that all the constituent mappings of the system (2.1a), and (2.1b) and the error equation (2.2), namely, 𝑓, 𝑔, 𝑝, 𝑠, β„Ž, and π‘ž are𝐢1 mappings vanishing at the origin, that is,𝑓(0)=0,𝑔(0)=0,𝑝(0)=0,𝑠(0)=0,β„Ž(0)=0,π‘ž(0)=0.(2.3) Thus, for 𝑒=0, the composite system (2.1a), and (2.1b) has an equilibrium state (π‘₯,πœ”)=(0,0) with zero error (2.2).

A state feedback controller for the composite system (2.1a), and (2.1b) has the form𝑒=𝛼(π‘₯,πœ”),(2.4) where 𝛼 is a 𝐢1 mapping defined on π‘‹Γ—π‘Š such that 𝛼(0,0)=0.

Upon the substitution of the feedback law (2.4) into the composite system (2.1a), and (2.1b), we get the closed-loop system given byΜ‡π‘₯=𝑓(π‘₯)+𝑔(π‘₯)𝛼(π‘₯,πœ”)+𝑝(π‘₯)πœ”,Μ‡πœ”=𝑠(πœ”).(2.5) The purpose of designing the state feedback controller (2.4) is to achieve both internal stability and output regulation. Internal stability means that when the input is disconnected from (2.5) (i.e. when πœ”=0), the closed-loop system (2.5) has an exponentially stable equilibrium at π‘₯=0. Output regulation means that for the closed-loop system (2.5), for all initial states (π‘₯(0),πœ”(0)) sufficiently close to the origin, 𝑒(𝑑)β†’0 asymptotically as π‘‘β†’βˆž. Formally, we can summarize the requirements as follows.

2.1. State Feedback Regulator Problem [2]

Find, if possible, a state feedback control law 𝑒=𝛼(π‘₯,πœ”) such that the following two conditions are satisfied.(OR1) (Internal Stability) The equilibrium π‘₯=0 of the dynamics Μ‡π‘₯=𝑓(π‘₯)+𝑔(π‘₯)𝛼(π‘₯,0)(2.6) is locally exponentially stable.(OR2)(Output Regulation) There exists a neighbourhood π‘ˆ of (π‘₯,πœ”)=(0,0) contained in π‘‹Γ—π‘Š such that for each initial condition (π‘₯(0),πœ”(0))βˆˆπ‘ˆ, the solution (π‘₯(𝑑),πœ”(𝑑)) of the closed-loop system (2.5) satisfies limπ‘‘β†’βˆž[]β„Ž(π‘₯(𝑑))βˆ’π‘ž(πœ”(𝑑))=0.(2.7)

Byrnes and Isidori [2] solved this problem under the following assumptions.(H1) The exosystem dynamics Μ‡πœ”=𝑠(πœ”) is neutrally stable at πœ”=0, that is, the system is Lyapunov stable in both forward and backward time at πœ”=0.(H2) The pair (𝑓(π‘₯),𝑔(π‘₯)) has a stabilizable linear approximation at π‘₯=0, that is, if 𝐴=πœ•π‘“ξ‚Ήπœ•π‘₯π‘₯=0ξ‚Έ,𝐡=πœ•π‘”ξ‚Ήπœ•π‘₯π‘₯=0,(2.8) then (𝐴,𝐡) is stabilizable, which means that we can find a gain matrix 𝐾 so that 𝐴+𝐡𝐾 is Hurwitz.

Next, we state the solution of the output regulation problem derived by Byrnes and Isidori [2].

Theorem 2.1 (see [2]). Under the assumptions (H1) and (H2), the state feedback regulator problem is solvable if and only if there exist 𝐢1 mappings π‘₯=πœ‹(πœ”) with πœ‹(0)=0 and 𝑒=πœ‘(πœ”) with πœ‘(0)=0, both defined in a neighbourhood π‘Š0βŠ‚π‘Š of πœ”=0 such that following equations (called the regulator equations ) are satisfied: (1)(πœ•πœ‹/πœ•πœ”)𝑠(πœ”)=𝑓(πœ‹(πœ”))+𝑔(πœ‹(πœ”))πœ‘(πœ”)+𝑝(πœ‹(πœ”))πœ”,(2)β„Ž(πœ‹(πœ”))βˆ’π‘ž(πœ”)=0. When the regulator equations (2.1a), (2.1b), and (2.2) are satisfied, a control law solving the state feedback regulator problem is given by []𝑒=πœ‘(πœ”)+𝐾π‘₯βˆ’πœ‹(πœ”),(2.9) where 𝐾 is any gain matrix such that 𝐴+𝐡𝐾 is Hurwitz.

3. Output Regulation of the Pan System

The Pan system is a new three-dimensional chaotic system discovered by Pan et al. [15] and described byΜ‡π‘₯1ξ€·π‘₯=π‘Ž2βˆ’π‘₯1ξ€Έ,Μ‡π‘₯2=𝑐π‘₯1βˆ’π‘₯1π‘₯3+𝑒,Μ‡π‘₯3=βˆ’π‘π‘₯3+π‘₯1π‘₯2,(3.1) where π‘Ž>0,𝑏>0,𝑐>0 are the parameters of the system and 𝑒 is the control.

The unforced Pan system (i.e. with 𝑒=0) is chaotic when π‘Ž=10,𝑏=8/3, and 𝑐=16. The chaotic state portrait of the unforced Pan system is illustrated in Figure 1.

In this paper, we solve the problem of output regulation for Pan system (3.1) for the tracking of constant reference signals (set-point signals).

The constant reference signals are generated by the exosystem dynamicsΜ‡πœ”=0.(3.2) It is important to note that the exosystem given by (3.2) is neutrally stable.

This follows simply because the differential equation (3.2) admits only constant solutions, that is,πœ”(𝑑)β‰‘πœ”(0)=πœ”0,βˆ€π‘‘βˆˆπ‘….(3.3) Thus, the assumption (H1) of Theorem 2.1 holds trivially.

Linearizing the dynamics of the Pan system (3.1) at the equilibrium (π‘₯1,π‘₯2,π‘₯3)=(0,0,0), we get the following system matrices:⎑⎒⎒⎒⎣⎀βŽ₯βŽ₯βŽ₯⎦⎑⎒⎒⎒⎣010⎀βŽ₯βŽ₯βŽ₯⎦𝐴=βˆ’π‘Žπ‘Ž0𝑐0000βˆ’π‘,𝐡=.(3.4) Using Kalman’s rank test for controllability [16], it can be easily seen that the pair (𝐴,𝐡) is not controllable. However, we shall show that the pair (𝐴,𝐡) is stabilizable, that is, we can find a gain matrix 𝐾 such that 𝐴+𝐡𝐾 is Hurwitz.

Suppose thatξ€Ίπ‘˜πΎ=1π‘˜2π‘˜3ξ€»,(3.5) where π‘˜1,π‘˜2, and π‘˜3 are real numbers to be determined so that 𝐴+𝐡𝐾 is Hurwitz.

The characteristic equation of the matrix⎑⎒⎒⎒⎣𝐴+𝐡𝐾=βˆ’π‘Žπ‘Ž0𝑐+π‘˜1π‘˜2π‘˜3⎀βŽ₯βŽ₯βŽ₯⎦00βˆ’π‘(3.6) is easily obtained asξ€Ίπœ†(πœ†+𝑏)2ξ€·+πœ†π‘Žβˆ’π‘˜2ξ€Έξ€·βˆ’π‘Žπ‘+π‘˜1+π‘˜2ξ€Έξ€»=0.(3.7) Since 𝑏>0, we note that πœ†=βˆ’π‘ is a stable eigenvalue of 𝐴+𝐡𝐾, and the other two eigenvalues of 𝐴+𝐡𝐾 will be stable if and only ifπ‘Žβˆ’π‘˜2ξ€·>0,βˆ’π‘Žπ‘+π‘˜1+π‘˜2ξ€Έ>0.(3.8) Since π‘Ž>0, the above conditions are equivalent toπ‘˜2<π‘Ž,𝑐+π‘˜1+π‘˜2<0.(3.9) Since π‘˜3 does not play any role in the above eigenvalue calculations, we can take π‘˜3=0.

Thus, we shall assume that [𝐾=𝐾1𝐾20], where π‘˜1 and π‘˜2 are chosen so that the inequalities (3.9) are satisfied. Thus, 𝐴+𝐡𝐾 is Hurwitz. Hence, (𝐴,𝐡) is stabilizable. Hence, the assumption (H2) of Theorem 2.1 also holds.

Hence, we can apply Theorem 2.1 to completely solve the output regulation problem for the Pan system for the tracking of constant reference signals (set-point signals).

3.1. Constant Tracking Problem for π‘₯1

Here, we consider the problem of the state π‘₯1 of the Pan system tracking the constant signal πœ”. Thus, the nonlinear system under consideration isΜ‡π‘₯1ξ€·π‘₯=π‘Ž2βˆ’π‘₯1ξ€Έ,Μ‡π‘₯2=𝑐π‘₯1βˆ’π‘₯1π‘₯3+𝑒,Μ‡π‘₯3=βˆ’π‘π‘₯3+π‘₯1π‘₯2,𝑒=π‘₯1βˆ’πœ”.(3.10) The regulator equations for the system (3.10) are obtained by Theorem 2.1 asπœ‹2(πœ”)βˆ’πœ‹1(πœ”)=0,π‘πœ‹1(πœ”)βˆ’πœ‹1(πœ”)πœ‹3(πœ”)+πœ‘(πœ”)=0,βˆ’π‘πœ‹3(πœ”)+πœ‹1(πœ”)πœ‹2πœ‹(πœ”)=0,1(πœ”)βˆ’πœ”=0.(3.11) Solving the regulator equations (3.11), we obtain the unique solutionπœ‹1(πœ”)=πœ”,πœ‹2(πœ”)=πœ”,πœ‹3πœ”(πœ”)=2π‘πœ”,πœ‘(πœ”)=βˆ’π‘πœ”+3𝑏.(3.12) Using Theorem 2.1 and the solution (3.12) of the regulator equations (3.11) for the nonlinear system (3.10), we obtain the following result which provides a state feedback control law solving the output regulation problem for the system (3.10).

Theorem 3.1. A state feedback control law solving the output regulation for the system (3.10) is given by []𝑒=πœ‘(πœ”)+𝐾π‘₯βˆ’πœ‹(πœ”),(3.13) where πœ‘(πœ”) and πœ‹(πœ”) are as given in (3.12), and the gain matrix [𝐾=𝐾1𝐾20] is chosen so as to satisfy the inequalities (3.9).

3.2. Constant Tracking Problem for π‘₯2

Here, we consider the problem of the state π‘₯2 of the Pan system tracking the constant signal πœ”. Thus, the nonlinear system under consideration isΜ‡π‘₯1ξ€·π‘₯=π‘Ž2βˆ’π‘₯1ξ€Έ,Μ‡π‘₯2=𝑐π‘₯1βˆ’π‘₯1π‘₯3+𝑒,Μ‡π‘₯3=βˆ’π‘π‘₯3+π‘₯1π‘₯2,𝑒=π‘₯2βˆ’πœ”.(3.14) The regulator equations for the system (3.14) are obtained by Theorem 2.1 asπœ‹2(πœ”)βˆ’πœ‹1(πœ”)=0,π‘πœ‹1(πœ”)βˆ’πœ‹1(πœ”)πœ‹3(πœ”)+πœ‘(πœ”)=0,βˆ’π‘πœ‹3(πœ”)+πœ‹1(πœ”)πœ‹2πœ‹(πœ”)=0,2(πœ”)βˆ’πœ”=0.(3.15) Solving the regulator equations (3.15), we obtain the unique solutionπœ‹1(πœ”)=πœ”,πœ‹2(πœ”)=πœ”,πœ‹3πœ”(πœ”)=2π‘πœ”,πœ‘(πœ”)=βˆ’π‘πœ”+3𝑏.(3.16) Using Theorem 2.1 and the solution (3.16) of the regulator equations (3.15) for the nonlinear system (3.14), we obtain the following result which provides a state feedback control law solving the output regulation problem for the system (3.14).

Theorem 3.2. A state feedback control law solving the output regulation for the system (3.14) is given by []𝑒=πœ‘(πœ”)+𝐾π‘₯βˆ’πœ‹(πœ”),(3.17) where πœ‘(πœ”) and πœ‹(πœ”) are as given in (3.16), and the gain matrix [𝐾=𝐾1𝐾20] is chosen so as to satisfy the inequalities (3.9).

3.3. Constant Tracking Problem for π‘₯3

Here, we consider the problem of the state π‘₯3 of the Pan system tracking the constant signal πœ”. Thus, the nonlinear system under consideration isΜ‡π‘₯1ξ€·π‘₯=π‘Ž2βˆ’π‘₯1ξ€Έ,Μ‡π‘₯2=𝑐π‘₯1βˆ’π‘₯1π‘₯3+𝑒,Μ‡π‘₯3=βˆ’π‘π‘₯3+π‘₯1π‘₯2,𝑒=π‘₯3βˆ’πœ”.(3.18) The regulator equations for the system (3.18) are obtained by Theorem 2.1 asπœ‹2(πœ”)βˆ’πœ‹1(πœ”)=0,π‘πœ‹1(πœ”)βˆ’πœ‹1(πœ”)πœ‹3(πœ”)+πœ‘(πœ”)=0,βˆ’π‘πœ‹3(πœ”)+πœ‹1(πœ”)πœ‹2πœ‹(πœ”)=0,3(πœ”)βˆ’πœ”=0.(3.19) Solving the regulator equations (3.19), we obtain the unique solutionπœ‹1√(πœ”)=π‘πœ”,πœ‹2√(πœ”)=π‘πœ”,πœ‹3√(πœ”)=πœ”,πœ‘(πœ”)=π‘πœ”(πœ”βˆ’π‘).(3.20) Using Theorem 2.1 and the solution (3.20) of the regulator equations (3.19) for the nonlinear system (3.18), we obtain the following result which provides a state feedback control law solving the output regulation problem for the system (3.18).

Theorem 3.3. A state feedback control law solving the output regulation for the system (3.18) is given by []𝑒=πœ‘(πœ”)+𝐾π‘₯βˆ’πœ‹(πœ”),(3.21) where πœ‘(πœ”) and πœ‹(πœ”) are as given in (3.20), and the gain matrix [𝐾=𝐾1𝐾20] is chosen so as to satisfy the inequalities (3.9).

4. Numerical Results

For numerical simulations, the parameters are chosen as in the chaotic case of Pan’s system, namely, 8π‘Ž=10,𝑏=3,𝑐=16.(4.1) For achieving internal stability of the state feedback regulator problem, a gain matrix 𝐾 must be chosen so that 𝐴+𝐡𝐾 is Hurwitz. Since 𝑏=8/3, one eigenvalue of 𝐴+𝐡𝐾 is always given by πœ†1=βˆ’π‘=βˆ’8/3. The other two eigenvalues of 𝐴+𝐡𝐾 are given by the quadratic equation (see (3.7))πœ†2ξ€·+πœ†π‘Žβˆ’π‘˜2ξ€Έξ€·βˆ’π‘Žπ‘+π‘˜1+π‘˜2ξ€Έ=πœ†2ξ€·+πœ†10βˆ’π‘˜2ξ€Έξ€·βˆ’1016+π‘˜1+π‘˜2ξ€Έ=0.(4.2) Suppose that we wish to choose 𝐾 so that it has two stable eigenvalues πœ†2=πœ†3=βˆ’3.

Then, we must haveπœ†2ξ€·+πœ†10βˆ’π‘˜2ξ€Έξ€·βˆ’1016+π‘˜1+π‘˜2ξ€Έ=πœ†2+6πœ†+9=0.(4.3) Comparing the coefficients of the two quadratic equations, we have10βˆ’π‘˜2=6,π‘˜1+π‘˜2+16=βˆ’0.9.(4.4) A simple calculation yieldsξ€Ίπ‘˜πΎ=1π‘˜20ξ€»=ξ€Ίξ€»βˆ’20.940.(4.5)

4.1. Constant Tracking Problem for π‘₯1

Here, the initial conditions are taken asπ‘₯1(0)=18,π‘₯2(0)=10,π‘₯3(0)=12,πœ”=2.(4.6) The simulation graph is depicted in Figure 2 from which it is clear that the state trajectory π‘₯1(𝑑) tracks the constant reference signal πœ”=2 in 2 seconds.

4.2. Constant Tracking Problem for π‘₯2

Here, the initial conditions are taken asπ‘₯1(0)=20,π‘₯2(0)=30,π‘₯3(0)=6,πœ”=2.(4.7) The simulation graph is depicted in Figure 3 from which it is clear that the state trajectory π‘₯2(𝑑) tracks the constant reference signal πœ”=2 in 3 seconds.

4.3. Constant Tracking Problem for π‘₯3

Here, the initial conditions are taken asπ‘₯1(0)=12,π‘₯2(0)=8,π‘₯3(0)=24,πœ”=2.(4.8) The simulation graph is depicted in Figure 4 from which it is clear that the state trajectory π‘₯3(𝑑) tracks the constant reference signal πœ”=2 in 2 seconds.

5. Conclusions

In this paper, we have derived new results for the output regulation problem of the Pan chaotic system (2010). Basically, we have applied the regulator equations of Byrnes and Isidori (1990) to derive state feedback control laws for regulating the output of the Pan system (2010) so as to track constant reference signals. We have provided simulation results to illustrate the effectiveness of the regulation schemes for the three cases of the constant regulation problem for the chaotic Pan system.