Abstract

This paper addresses the exponential synchronization problem of a class of master-slave distributed parameter systems (DPSs) with spatially variable coefficients and spatiotemporally variable nonlinear perturbation, modeled by a couple of semilinear parabolic partial differential equations (PDEs). With a locally Lipschitz constraint, the perturbation is a continuous function of time, space, and system state. Firstly, a sufficient condition for the robust exponential synchronization of the unforced semilinear master-slave PDE systems is investigated for all admissible nonlinear perturbations. Secondly, a robust distributed proportional-spatial derivative (P-sD) state feedback controller is desired such that the closed-loop master-slave PDE systems achieve exponential synchronization. Using Lyapunov’s direct method and the technique of integration by parts, the main results of this paper are presented in terms of spatial differential linear matrix inequalities (SDLMIs). Finally, two numerical examples are provided to show the effectiveness of the proposed methods applied to the robust exponential synchronization problem of master-slave PDE systems with nonlinear perturbation.

1. Introduction

The study on master-slave systems has become more important for theoretical and practical points in many fields, including communication, mechanical systems, robotics, chemical reactions, and biological systems [19]. Ever since the discovery of Christian Huygens in 1665 on the synchronization of two pendulum clocks [10], synchronization has received considerable attention for a long time as a typical collective behavior and a basic motion in nature with potential applications in many different areas including secure communication, chaos generators design, chemical reactions, biological systems, and information science [1126]. The theory of synchronization for master-slave systems, which aims to control the slave system so that the output of the slave system follows the output of the master system [27], is a recent research area extensively investigated nowadays in many industrial and technical processes, such as unmanned aerial vehicle (UAV) team, vehicular platoons, rendezvous of space shuttles, and many other practical control systems (see, e.g., [28, 29] and the references therein).

The existing works most considered dynamical behavior of master-slave systems described by ordinary differential equations (ODEs) or delay differential equations (DDEs), and a variety of synchronization criteria have been presented [1416]. In practice, however, the outputs, inputs, and process states with relevant parameters usually vary both temporally and spatially in nature, and thereby their behavior depending on time and spatial position could be described by distributed parameter systems (DPSs) modeled by partial differential equations (PDEs). Unfortunately, few works have investigated the synchronization of master-slave PDE systems.

As a result of the infinite-dimensional nature of master-slave PDE systems, the existing finite-dimensional control theory and techniques for the master-slave ODE systems are difficult to be directly employed for the control design of master-slave PDE systems. In this situation, it is important to study the synchronization problem of master-slave PDE systems. Therefore, some researchers have paid attention to the study of synchronization of master-slave PDE systems [1729], where “design-then-reduce” approach was employed to take the full advantage of the original PDE model for the controller design [3034]. References [3537] researched synchronization of neural networks with reaction-diffusion terms. Yuan et al. proposed synchronization of the coupled distributed parameter system with time delay via P-sD control [38]. Yang et al. proposed exponential synchronization for complex spatiotemporal networks with space-varying coefficients via P-sD control [39]. Wang et al. studied exponential synchronization for networked linear PDE systems via boundary control [40]. However, to the best of our knowledge, few results are available on the exponential synchronization of master-slave PDE systems with nonlinear perturbation and spatially variable coefficients, which motives the present investigation.

In this paper, we will deal with the problem of robust exponential synchronization for master-slave parabolic PDE systems with spatially variable coefficients and spatiotemporally variable nonlinear perturbation. Initially, a master-slave parabolic PDE model is discussed, and then a synchronization error dynamic of PDE systems is developed in spatial domain. Then, on the basis of Lyapunov’s direct method and the technique of integration by parts, a spatial differential linear matrix inequality (SDLMI) based condition for the robust exponential synchronization of the unforced semilinear master-slave PDE systems is studied. Once the unforced semilinear master-slave PDE systems cannot achieve robust exponential synchronization, distributed proportional-spatial derivative (P-sD) control design is developed to achieve that the closed-loop slave system which is exponentially synchronized with the master system with a given decay rate for all admissible nonlinear perturbations in terms of SDLMI. Furthermore, the SDLMI optimization problem can be approximately solved by the finite difference method and LMI optimization techniques [41, 42]. Finally, the simulation study on the exponentially synchronous control of a master-slave PDE systems with nonlinear perturbation is given to show the effectiveness of the proposed design method.

The remainder of this paper is organized as follows. The problem formulation and preliminaries are given in Section 2. Section 3 provides the exponential synchronization analysis of the unforced master-slave PDE systems. Robust P-sD controller is designed in Section 4. Section 5 presents two examples on master-slave PDE systems to illustrate the effectiveness of the proposed method. Finally, Section 6 offers some concluding remarks.

Notations. The following notations will be employed throughout this paper. , , and denote the set of all real numbers, n-dimensional Euclidean space, and the set of all matrices, respectively. and denote the Euclidean norm and inner product for vectors, respectively. Identity matrix, of appropriate dimensions, will be denoted by . For a symmetric matrix , (, , resp.) means that is positive definite (negative definite, negative semidefinite, resp.). The space-varying symmetric matrix function , , is positive definite (negative definite, negative semidefinite, resp.), if (, , resp.) for each . represents the eigenvalue of a matrix. The superscript “” is used for the transpose of a vector or a matrix. is a Sobolev space of absolutely continuous n-dimensional vector functions with square integrable derivatives of the order and with the norm . The symbol “” is used as an ellipsis in matrix expressions that are induced by symmetry; for example,

2. Problem Formulation

Consider the following synchronization scheme of semilinear master-slave PDE systems with spatially variable coefficients and spatiotemporally variables as follows.

The master system is as follows:and the slave system   is as follows:where , are the state vectors, the subscripts and stand for the partial derivatives with respect to , , respectively, and are the spatial position and time, respectively, and is the distributed control input. , , and , , are real known matrix functions. , , is a real known matrix function. is n-dimensional real vector-valued spatiotemporally variable nonlinear perturbation with for all and . Moreover, assume that satisfies the following Lipschitz condition as follows.

Assumption 1. is continuous with respect to its arguments and is locally Lipschitz continuous in on ; that is, there exists a scalar such that the following inequality holds for all , and : where . Therefore, we have thatDefine the synchronization error , and we have the following synchronization error system :where .
Choose as the state space and the trajectory segment as the state. Define the spatial differential operator in aswith its domainLet be control input space and the trajectory segment as the control input. Then the state-space description of the synchronization error system (6) can be rewritten as the following nonlinear abstract differential equation on the Hilbert space [43]:where , , , and .
Using (8), we get from (6) thatwhere .

We introduce the following definition of the exponential synchronization for the master-slave systems (2) and (3) in the sense of norm .

Definition 2. For a given constant , the master-slave systems (2)-(3) achieve ρ-exponential synchronization or exponential synchronization with a given decay rate  ρ, if there exists a constant such that the following inequality holds for any initial condition , :

It is easily seen from Definition 2 that the master-slave systems (2)-(3) achieve exponential synchronization with a given decay rate ρ if and only if the error system (6) is exponentially stable with a given decay rate ρ.

For simplicity, when , the master-slave systems (2)-(3) are referred to as unforced master-slave systems.

Lemma 3. For any two square integrable vector functions , , , the following inequality holds for any positive scalar function defined on interval :

Proof. It is easily found that the inequality holds for any . Hence, which impliesIntegrating inequality (14) from to , we can obtain thatwhich means that inequality (12) holds. The proof is complete.

3. Exponential Synchronization Analysis

This section aims to analyze exponential synchronization for the unforced semilinear master-slave PDE systems (2)-(3).

We consider the following Lyapunov functional for the unforced synchronization error system (10):where is a real matrix function to be determined. The time derivative of along the solution of system (10) is given byBy integrating by parts and taking into account the boundary condition in (10), we haveAccording to the definition of in (7) and considering (18), we haveBy Assumption 1 and Lemma 3, for any positive scalar function , , we have Substituting (19), (20) into (17), we obtainwhere and

Theorem 4. Consider the unforced semilinear master-slave PDE systems (2) and (3) under Assumption 1. For a given scalar , the master-slave PDE systems (2) and (3) achieve exponential synchronization with a given decay rate ρ, if there exist a matrix function and a scalar function satisfying the following SDLMI in space:where

Proof. Assume that SDLMI (23) is satisfied for matrix functions and . Using Schur complement, the following inequality is achieved if SDLMI (23) is satisfied for each :Substituting inequality (25) into (21) yields Integration of (26) from to derivesSince is a spatially continuous matrix function of defined on , it is easily observed that given by (18) satisfies the following inequality:where and are two positive constants. Using (28), we can get the following relation:Therefore, we haveThus, from (30) and Definition 2, the unforced semilinear master-slave PDE systems (2) and (3) achieve exponential synchronization with a given decay rate ρ. The proof is complete.

Theorem 4 presents an SDLMI-based condition for the robust exponential synchronization of the unforced semilinear master-slave PDE systems (2) and (3). Once the unforced semilinear master-slave PDE systems (2) and (3) cannot exponentially synchronize by themselves, it is desired to design a distributed P-sD controller in the semilinear master-slave PDE systems (2) and (3). Moreover, the SDLMI feasibility problem is approximately solved via the finite difference method and the existing LMI optimization techniques [41, 42].

4. P-sD Control Design

Once the semilinear master-slave PDE systems (2)-(3) are not exponentially synchronized by themselves, a state feedback control design is therein desired. The aim of this section is to propose a distributed P-sD state feedback control design method to achieve exponential synchronization for semilinear master-slave PDE systems (2)-(3).

This study considers the following distributed P-sD state feedback controller of the slave system (3) as shown in Figure 1:where , , , are continuous matrices to be determined.

Remark 5. The controller (31) of the slave system needs distributed sensing and actuation, whereas it has become feasible nowadays to produce integrated control circuitry with a large number of microsensors and actuators. It has been pointed out in [33, 34, 38, 39] that the controller (31) with spatial derivative term provides more spatial performance and that is why the controller (31) is chosen in this paper.
Substituting (31) into (10) gives the following closed-loop synchronization error system:where .
Hence, the objective of this section is to develop a design method of the distributed state feedback controller (31) to ensure exponential synchronization of the master-slave systems (2)-(3), that is, the exponential stability of the closed-loop system (32). To achieve the goal, the following assumption is necessary.

Assumption 6. The matrix in (2) and (3) is diagonal; that is, , .

Assumption 6 is representative to some extent and is widely used in the literature, including the FHN equation and the tubular reactor with axial dispersion and heat exchangers [33, 34, 39, 4446].

Consider the Lyapunov functional (16), and the time derivative of along the solution of system (32) is given byBy Assumption 1 and Lemma 3, for any positive scalar function , , similar to (20), we have Substituting (19), (34) into (33) and considering defined in (32), we obtainwhere is defined in (22) and

Theorem 7. Consider the closed-loop semilinear master-slave PDE systems (2) and (3) under Assumptions 1 and 6. For a given scalar , the distributed P-sD controller (31) can guarantee that the master-slave PDE systems (2) and (3) achieve exponential synchronization with a given decay rate ρ, if there exist matrix function , a scalar function , and matrix functions and satisfying the following SDLMI in space:whereIn this case, the gain matrix functions , , , of the suitable controller (31) can be given as

Proof. Assume that SDLMI (37) is satisfied for , , , and . Using the Schur complement, SDLMI (37) is equivalent to the inequality where .
Lettingpre- and postmultiplying (40) by , respectively, with application of Assumption 6, and considering the property , , and the commutative law of matrix multiplication [47], we can obtainSince , , we can get the following inequality from (40) and (42):Substituting inequality (43) into (35) yields The rest of proof of Theorem 7 is similar to that of Theorem 4, and thus we can also get the conclusion that the semilinear master-slave PDE systems (2) and (3) achieve exponential synchronization with a given decay rate ρ with the suitable controller (31). Moreover, from (41), we have (39). The proof is complete.

Theorem 7 presents an SDLMI-based condition for the existence of a distributed controller (31) for the exponential synchronization of the semilinear master-slave PDE systems (2) and (3). Explicit expressions of a desired controller are proposed when the SDLMI (37) is feasible.

Remark 8. Notice that an SDLMI control design has been more recently proposed in [33, 34] for a class of PDE systems. Different from the stable control design for PDE system in [33, 34], this paper provides an SDLMI-based sufficient condition on the exponential synchronization for master-slave PDE systems.

Remark 9. Notice that [3840] dealt with exponential synchronization of linear models, while this paper considers exponential synchronization of the model with nonlinear perturbation.

Remark 10. For a special case, when the master system : , we have , and then the exponential synchronization of the master-slave systems is equal to exponential stability of the slave system . In other words, Theorem 4 provides the sufficient condition of exponential stability of the slave system with a given decay rate ρ when the master system , while Theorem 7 also contains the way of the P-sD controller (31) design for exponential stability.

5. Numerical Simulation

In this section, in order to show the effectiveness of Theorems 4 and 7, we consider the following two examples.

Example 1. Consider semilinear master-slave PDE systems (2) and (3) with the following spatial variable coefficients: and the following initial conditions:where .
Let . The parameter satisfying Assumption 1 is chosen asDivide the spatial domain into space instances of the same distance, where . With application of Theorem 4, employing Algorithm  1 in [34] to solve SDLMI (23) and setting and , we get feasible solutions of and . Therefore, the unforced semilinear master-slave PDE systems (2) and (3) with coefficients (45) and initial conditions (46) achieve robust exponential synchronization with a given decay rate ρ, which is apparent from Figure 2.

Example 2. Consider semilinear master-slave PDE systems (2) and (3) with the following spatially variable coefficients: and the following initial conditions:where .
Figure 3 shows the unforced profiles of evolution of the semilinear master-slave PDE systems (2) and (3) with the parameters given in (40) under the initial condition given in (41). It is apparent from Figure 3 that the operating steady state of the master PDE system (2) and the operating steady state of the slave PDE system (3) are not synchronized.
Let . The parameter satisfying Assumption 1 is chosen asIn this example, the distributed P-sD state feedback controller (31) is employed, where , , , are matrix functions to be determined. Divide the spatial domain into space instances of the same distance, where . With application to Theorem 7, employing Algorithm  1 in [34] to solve SDLMI (37) and setting and , we get feasible solutions of , , , and . Thereby, we can get the approximated control gain matrix functions and by (39) as shown in Figures 4 and 5, respectively.
Applying the controller (31) with the approximated control gain matrices and shown in Figures 4 and 5 to the semilinear master-slave PDE systems (2) and (3) with the parameters given in (48) under the initial condition given in (49), Figure 6 indicates the closed-loop profiles of evolution of the semilinear master-slave PDE systems of and , respectively.
It is clear from Figure 6 that the proposed distributed state feedback controller (31) can ensure that the semilinear master-slave PDE systems (2) and (3) with the parameters given in (48) under the initial condition given in (49) achieve exponential synchronization with a given decay rate ρ. Moreover, the evolution profile of the distributed control is shown in Figure 7.

6. Conclusions

In this paper, we have addressed the exponential synchronization problem of semilinear master-slave PDE systems with spatially variable coefficients. An SDLMI-based sufficient condition for the robust exponential synchronization of the unforced semilinear master-slave PDE systems has been presented. Moreover, when the master-slave PDE systems are not exponentially synchronized, it has been proved in detail that P-sD controllers with appropriate gain parameters can guide the closed-loop master-slave PDE systems to exponentially synchronize with a given decay rate ρ if a given SDLMI condition is fulfilled. Finally, numerical examples showed the effectiveness of the proposed methods.

Conflict of Interests

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

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China under Grant nos. 61273012, 61403179, and 61304023, in part by the Natural Science Foundation of Shandong Province of China under Grant nos. ZR2014AL009 and ZR2014CP008, in part by the Science and Technology Development Plan Project of Shandong Province under Grant no. 2012YD01052, and in part by the Applied Mathematics Enhancement Program of Linyi University.