# LMI-based model predictive control for underactuated surface vessels with input constraints.

1. Introduction

Over the past decade, control of underactuated systems has been one of active research areas in control society. One feature of underactuated systems is that the number of independent actuators of the system is less than that of the degree of freedom. The challenging control problem is how to design a stabilizing control law such that the state of the closed-loop system asymptotically converges to the origin. In this note, we consider the dynamic positioning-control problem for a ship that has no side thruster, but two independent main thrusters located in a distance from the center line in order to provide both surge force and yaw moment. The control problem considered in this note is to find a feedback law that stabilizes both the position variables and the orientation, using only the two available controls. Since we attempt to control three degrees of freedom with only two independent controls, we have an underactuated control problem.

The dynamics of this underactuated system is complex enough to yield a rich source of control problems, yet simple enough to permit a complete mathematical analysis. It has been shown in  that underactuated vessels do not satisfy Brockett's necessary condition  if the unactuated dynamics contain no gravitational field component. In this case, these vessels are not asymptotically stabilizable to a given equilibrium configuration via time-invariant continuous feedback laws.

Stabilization of underactuated surface vessels has been tackled in a number of research studies in the last few years [3-6]. The controller in , which can be viewed as a time-varying controller devised to circumvent Brockett's limitation, is proposed. It consists of a supervising logic unit and subordinate time-invariant controllers. The asymptotic stability is rigorously proven using Lyapunov analysis. A sliding mode control (SMC) law  is presented by introducing a first-order sliding surface in terms of surge tracking errors and a second-order surface in terms of lateral motion tracking errors. In , a discontinuous feedback control law has been discussed for the exponential convergence of the equilibrium point under certain assumptions on the variables initial value. Contrastively, several transformations in  are introduced to represent the system into a nonlinear cascade form. One of the desirable features of the suggested approach in  is that it provides a systematic procedure to transform a class of nonholonomic dynamical systems into a nonlinear cascade form. Such representation opens the possibility of developing more constructive control tools than those already existing in . However, it is a pity that the control design in  is not complete and lacked a further step in the backstepping procedure. In the comment letter , the control laws in  were revised and the states decayed asymptotically to zero. In , a discontinuous control approach with two stage control laws switched on at given time is proposed based on the stability analysis of the global transformed system. With the aid of terminal sliding mode (TSM) method, the two asymptotically stabilizing control laws for the transformed subsystems are designed independently to asymptotically stabilize the underactuated surface vessel. A controller designed by using Lyapunov direct method and the popular backstepping techniques  is utilized to stabilize the USV system to the origin.

In realistic implementations, the performance of control systems is often limited due to constraints on inputs or states that naturally arise. None of the references cited above has taken those constraints explicitly into account. MPC is a control technique which embeds optimization within feedback to deal with systems subject to constraints on inputs and states . Over the last few decades, MPC has proven to be successful for a wide range of applications including chemical, food processing, automotive, and aerospace systems . The fact that the rudder actuation is limited in amplitude and rate makes MPC approach a natural choice for the design of the path following controller. By taking these physical constraints into account, control actions that respect actuators limits can be generated. As the pioneer work of MPC application in control of surface vessels, authors in  consider rudder saturation in their MPC controller and adopt a 1DoF yaw dynamical model. By using model predictive control algorithm, a rudder roll stabilization (RRS) system is presented in . The control objective is to regulate the heading to a desired value and to reduce the roll angle under various sea conditions. In , a nonlinear MPC control law with terminal invariant manifolds constraints is designed through coordinate transformation and state feedback transformation based on diffeomorphism and Lyapunov stability theory. Based on MPC, the states of the USV are steered into a terminal manifolds set. After the terminal manifolds set is reached, a linear feedback control is used to stabilize the system. However, just as the authors said in their paper, the control design proposed is developed with constant parameters without considering parameter perturbation. In , the focus is on satisfying all the input (rudder) and state (roll) constraints by using MPC to achieve satisfactory path following performance. Authors in  present an experimental implementation of a model predictive control (MPC) strategy for path following on a model ship. The objective of the proposed algorithm is to control cross-tracking error and heading angle using both linearized and nonlinear models. To implement MPC with a nonlinear model, the integrated perturbation analysis and sequential quadratic programming (InPA-SQP) method are introduced to solve the constrained optimal control problem.

For notational convenience, the ship dynamics in [19, 20] are written into linearized and nonlinear form based on the assumption that surge velocity is constant and the yaw moment is proportional to the rudder angle. However, the former assumption is hardly admissible in engineering.

This paper illuminates the stabilization approach for underactuated surface vessels with only a surge force and a yaw moment. The technique in this paper leads to convex LMIs based online optimization problem when the USV model is described by SDC in terms of diffeomorphism and SDRE theory. Moreover, the constrained control input, index, and stability of the considered system are solved by convex optimization based on LMIs. The simulation results involving an actual underactuated ship Northern Clipper  verified the effectiveness of the proposed method.

Notation. The symbol * will be used in some matrix expressions to induce a symmetric structure. I denotes identity matrix. For example, when H and R are symmetric matrices, then

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (1)

2. Problem Formulation

In this paper, we consider the trajectory tracking control problem of a surface vessel shown in Figure 1. There is no side thruster, but two independent main thrusters located in a distance from the center line in order to provide both surge force and yaw moment.

The dynamics and kinematics of an underactuated surface vessel are described as 

M[??] + C (v) v + Dv = [tau] (2)

[??] = J([eta]) v. (3)

The inertia matrix M = diag{[m.sub.11], [m.sub.22], [m.sub.33]} and the damping matrix D = diag{[d.sub.11], [d.sub.22], [d.sub.33]} are constant and positive definite. The vector [tau] = [[[tau].sub.1], [[tau].sub.2], [[tau].sub.3]] denotes the control forces in surge and sway and control torque in yaw. In this paper, the surface vessel is considered as the common thruster configuration that has no side thruster, such as [[tau].sub.2] = 0. So, the second component equation of (2) behaves as a nonholonomic constraint, which is a nonintegrable relation involving not only the generalized coordinates and velocities but also the generalized accelerations 23]. C(v) is the matrix of Coriolis and centripetal terms also including added mass. J([eta]) is the rotation matrix for the transformation between body-fixed and earth-fixed coordinates. Consider

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (4)

The vector [eta] = [[x, y, [psi]].sup.T] denotes the positions and orientation of the underactuated surface vessel in the earth-fixed coordinate system. The vector v = [[u, v, r].sup.T] denotes the linear velocities in surge and sway and the angular velocity in yaw.

As a general accepted conclusion (see [24, 25]), there is no continuous time-invariant feedback control law that makes the zero origin an asymptotically stable equilibrium of the system (2)-(3), for the system does not satisfy Brockett's condition . Then time-varying and discontinuous control approaches are only taken into account in this paper.

Neglecting the motions in heave, roll, and pitch, the simplified kinematic model which describes the geometrical relationship between the earth-fixed (E-frame) and the body-fixed (B-frame) motion is given as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (5)

The following global coordinate transformation and feedback transformation are adopted before control design. Define

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (6)

It is proved that the state transformation (6) is a global diffeomorphism (see ).

The feedback transformation is

[w.sub.1] = ([d.sub.11]/[d.sub.22] - 1)u - [z.sub.3][z.sub.6] - [[tau].sub.1]/[d.sub.22]

[w.sub.2] = ([m.sub.11] - [m.sub.22])uv/[m.sub.33] - [d.sub.33]r/[m.sub.33] + [[tau].sub.3]/[m.sub.33]. (7)

With the state and feedback transformation (6)-(7), the system (2)-(3) is eventually transformed to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (8)

The system (8) has the same diffeomorphism properties as the system (2)-(3) (see 22]); that is, if [lim.sub.t[right arrow][infinity]] [z.sub.i] = 0, (1 [less than or equal to] i [less than or equal to] 6) then (x, y, [psi], u, v, r) converges to zero as t [right arrow] [infinity].

3. Control Design

MPC is a popular technique for the control of dynamical system subject to input and state constraints. At any time instant, MPC requires the online solution of an optimization problem to compute optimal control inputs over a fixed number of future time instants, known as the finite horizon or quasi-infinite horizon. Using MPC, it is possible to handle inequality constraints on the manipulated and controlled variables in a systematic manner during the design and implementation of the controller (see ).

We will consider a nonlinear plant with input constraints, where the evolution of the state after time is predicted by the following model:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (9)

The data of this model comprise a set W [member of] [R.sup.n] containing all possible terminal states with constraints, a vector X that is the set of the state constraints, and a given vector U of possible set of control constraints. Our objective is to obtain a feedback law that (asymptotically) drives the state of our plant to the origin. In this paper, we consider the following quadratic objective:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10)

The domain of this optimization problem is the set of admissible processes, namely, pairs (x, u), comprising a measurable control function u and the corresponding absolutely continuous state trajectory x which satisfy the constraints of min J(x, t, u), where function F(*) [greater than or equal to] 0 is selected as terminal costs. We show that we can guarantee stability of the resultant closed-loop system, by choosing the design terminal constraints set W and terminal costs F(*) to satisfy a certain stability condition. However, in the above technique, selecting nonlinear affine system (8) as the predictive model may give rise to a lot of problems such as complicated computation and difficulty solving the optimization; the problems inherent to nonlinear affine model are inevitable. In this work, we will make efforts to design MPC controllers involving SDC model to overcome the previous problem (see ). System (8) is equivalent to the following nonlinear affine model:

[??] = f (x) + g (x) u. (11)

It is well-known that the nonlinear dynamics (8) can be represented by the following linear structure having SDC (state-dependent coefficients) by using SDRE theory (see ):

[??] = A(x) x + B (x) u, (12)

where f(x) = xA(x) and B(x) = g(x). It is also known  that, in the multivariable case, there are an infinite number of ways to bring the nonlinear system to SDC form. It is clear that the modeling errors caused by the traditional linear methods are avoided.

Definition 1. A(x) is a controllable (stabilizable) parameterization of the nonlinear system (in a region X) if the pair (A(x), B(x)) is pointwise controllable (stabilizable) in the linear sense for all x [member of] X (see ).

Construct the nonlinear feedback controller via

u (x) = -K (x) x = -[R.sup.-1] [B.sup.T] (x) P (x) x. (13)

Solve the state-dependent Riccati equation (SDRE)

[A.sup.T] (x) P(x) + P (x) A (x) + Q (x) -P(x)B (x) [R.sup.-1] (x) [B.sup.T] (x) P(x) = 0 (14)

to obtain P(x) [greater than or equal to] 0. Note that P(x) is a solution of x.

Physical limitations inherent in ship power transmission system invariably impose hard constraints on the manipulated variable input. In this section, we show how limits on the control signal can be incorporated into our MPC algorithm as sufficient LMI constraints. Predictive control is an effective optimization control method to deal with constrains. The basic idea of the discussion that follows can be found in  in the context of discrete-time robust MPC. We present it here to clarify its application in our underactuated surface vessels setting and also for completeness of exposition.

The nonlinear dynamics (8) can be represented by the following linear structure having state-dependent coefficients:

[??] = A(z) z + B(z) w, (15)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (16)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (17)

The Euler discrete system of (15) is given by

z(k + 1) = A(z, k) z(k) + B (z, k) w (k) [[parallel]w[parallel].sub.2] [less than or equal to] [w.sub.max]. (18)

We define the nominal cost as

[J.sub.[infinity]] = [[infinity].summation over (i=0)] [z.sup.T] (k + i | k)Qz (k + i | k) + [w.sup.T] (k + i | k) Rw (k + i|k), (19)

where Q = [Q.sup.T] > 0 and R = [R.sup.T] > 0. MPC is an attractive strategy for systems subject to input constraints. Due to the systems with constrained control input, we obtain the control input by minimizing a nominal cost [J.sub.[infinity]] over an infinite predictive horizon as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (20)

It is very difficult to resolve the optimization (20) directly. We can see that, based on LMI constraints, it is easy to get the optimization in system (20).

Theorem 2. Consider the time-varying system (18)-(20), where z(k | k) is the state of the system (15) measured at sampling time k. [w.sub.max] is the maximum value of the input w. When A(x) is a controllable parameterization of the system (18), if there exist F, Q, [gamma], Q = [Q.sup.T] [member of] [R.sup.6x6] >0, [gamma] [member of] [R.sup.1x1] > 0, and Y [member of] [R.sup.2x6], the state feedback matrix F in the control law w = Fz(k) that minimizes the upper bound on the performance objective function at sampling time k is given by F = Y[Q.sup.-1] z(k), where Q and Y are obtained from the solution (if it exists) of the following linear objective minimization problem. If the optimization problem in (21)-(23) is feasible, the closed-loop underactuated system (11) with input constraints [[parallel]w[parallel].sub.2] [less than or equal to] [w.sub.max] is asymptotically stable. Consider

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (21)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (22)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (23)

Proof. Consider the following.

(1) Proof for Optimization in Performance Index. Using quadratic Lyapunov-Krasovskii function, the upper bound on the objective function [J.sub.[infinity]] is given as

V (z (k)) = [z.sup.T] (k) P (k) z (k) P (k) > 0. (24)

To guarantee the existence of the upper bound on the performance index, underlying inequalities must be given as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (25)

with the conditions V(x([infinity] | k)) = 0 and x([infinity] | k) = 0. Summing (22) from 0 to [infinity] gives

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (26)

where Q = [gamma][P.sup.-1] > 0.

Using the Schur-complement , we get (21).

(2) Proof for the Stability in System (18)-(19). Substituting (18) to (19) gives

[z.sup.T](k) [P - [(A + BF).sup.T] P(k) (A + BF) -[Q.sub.1] - [F.sup.T] RF] z (k) > 0, (27)

where w = Fz(k) = Y[Q.sup.-1]z(k); we get

P - [(A + BF).sup.T]P (k) (A + BF) - [Q.sub.1] - [F.sup.T]RF > 0. (28)

Using the Schur-complement, we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (29)

Substituting F = Y[Q.sup.-1] and pre- and postmultiplying by Q give LMI (22).

(3) Proof for Input Constraints. Consider the Euclidean norm constraint [[parallel]w[parallel].sub.2] [less than or equal to] [w.sub.max]. The constraint is imposed on the present and the entire horizon of future manipulated variables. Following , we have that

V (z(k + i | k)) = [z.sup.T] (k + i | k) [Q.sup.-1]z (k + i | k) [less than or equal to] [z.sup.T] (k|k)[Q.sup.-1]z(k|k) [less than or equal to] 1 (30)

is state-invariant ellipsoid. So, we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (31)

Using the Schur-complement, we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (32)

Remark 3. The applicable conditions of Theorem 2 is that A(x) is a controllable parameterization of the system (18). This means that the condition [z.sub.6] [not equal to] 0 must be satisfied. When[ [z.sub.6] = 0, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (33)

where the pair (A(x), B(x)) is not pointwise controllable. To solve the prevous problem, we define a constant [z.sup.*.sub.6] which is close to zero but not equal to zero. When [z.sub.6] is close to zero, let [z.sub.6] = [z.sup.*.sub.6]; the time-varying state matrix A(z) is substituted as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (34)

Substituting (34) to 22) gives

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (35)

So, we get the following Theorem 4.

Theorem 4. Consider the time-varying system (18)-(20), where z(k | k) is the state of the system (15) measured at sampling time k. [w.sub.max] is the maximum value of the input w. If there exist F, Q, [gamma], Q = [Q.sup.T] [member of] [R.sup.6x6] >0, [gamma] [member of] [R.sup.1x1] > 0, and Y [member of] [R.sup.2x6], the state feedback matrix F in the control law w = Fz(k) that minimizes the upper bound on the performance objective function at sampling time k is given by F = Y[Q.sup.-1]z(k), where Q and Y are obtained from the solution of the following linear objective minimization problem. Simultaneously, input constraints [[parallel]w[parallel].sub.2] [less than or equal to] [w.sub.max] are satisfied.

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (36)

Remark 5. After obtaining the MPC controller w, we can get the control input for the origin system (5) by inverse diffeomorphism

[[tau].sub.1] = ([d.sub.11] - [d.sub.22]) u - [d.sub.22][z.sub.3][z.sub.6] - [d.sub.22][w.sub.1] [[tau].sub.3] = [m.sub.33][w.sub.2] - ([m.sub.11] - [m.sub.22]) uv + [d.sub.33]r. (37)

4. Simulation Results

In , an underactuated actual ship named "Northern Clipper" is introduced. Consider the simulation model with parameters as Northern Clipper in : [m.sub.11] = 5.312 x [10.sup.6] kg, [m.sub.22] = 8.283 x [10.sup.6] kg, [m.sub.33] = 3.745 x [10.sup.6] kg, [d.sub.11] = 5.024 x [10.sup.4] kg/s, [d.sub.22] = 2.722 x [10.sup.5] kg/s, and [d.sub.33] = 1.189 x [10.sup.8] kg/s. Length of "Northern Clipper" L = 76.2 m, and mass m = 4.6 x [10.sup.6] kg. This suggests that the system matrices should be scaled with respect to forward speed, ship length L, and displacement (see ). The model parameters can be made nondimensional as [m".sub.11] = 1.1274, [m".sub.2] = 1.8902, [m".sub.33] = 0.1278, [d".sub.11] = 0.0358, [d".sub.22] = 0.1183, and [d".sub.33] = 0.0308. In this section, the effectiveness of the proposed MPC control law is verified by following simulation. The control law is selected as Theorem 4, and the initial values are selected as x(0) = -2, y(0) = -2, [psi](0) = 1, u(0) = v(0) = 0 m/s, r(0) = 1, [w.sub.max] = 1, [z.sup.*.sub.6] = 0.001, simulation time T = 100 s, and sampling time 0.05 s.

Simulation results are shown in the following figures, and the simulation time is set to 100 s. Figure 2 gives the time response of the positions v, y and orientation [psi]. Figure 3 gives the time response of the velocities w, v and orientation velocity r; Figure 4 gives the time response of [w.sub.1], [w.sub.2]. The responses of the control inputs [[tau].sub.1] and [[tau].sub.3] are shown in Figure 5. The trajectory is shown in Figure 6. The switching motion of the control signal has been highlighted as red "o" on the simulation Figures 2-5. Because [z.sup.*.sub.6] is close to zero, the switching procedure is relatively smooth. Figures 2-6 show that the MPC control laws asymptotically stabilize the underactuated surface vessel to the zero origin. Moreover, the constraint on the control input is satisfied.

5. Conclusion

This paper proposes the stabilization approach for underactuated surface vessels with only a surge force and a yaw moment. The stability theory of MPC controller is further developed, which is applied to the stabilization control of underactuated surface vessel with input constraints. For the stabilization control of underactuated surface vessel, a nonlinear MPC control law is designed through coordinate transformation and SDRE based on LMIs. The simulation results show that the proposed control law can effectively deal with the problem of stabilization control of underactuated surface vessel.

It should be noted that solving LMI itself is a challenge in computation. Many research efforts in LPV community are directed to address the computational issues. For LPV, those large LMIs are only solved once. But for MPC, the LMIs have to be solved repeatedly at each sampling time. In the future work, efforts will be made to analyze that whether the algorithms are fast enough to do this at the rate required by the control sampling. As is well known, for a long time, efficient setting of the large set of tunable parameters has been a hard problem for MPC. Many available methods to deal with uncertain systems have already been developed (see [30-33] and the references therein for more details). However, it should be noticed that MPC algorithm may perform very poorly when model mismatch occurs in spite of the inherent robustness provided by the feedback strategy based on the plant measurement at the next sampling time. MPC involving data-driven technique is suitable to overcome the previous problem [34-38]. In this paper, the robustness of the proposed approach was not discussed in detail. Fortunately, MPC requires the online solution of an optimization problem to compute optimal control inputs over a fixed number of future time instants, which partly guarantee the robustness of the dynamics. For the perturbed impulsive USV system with uncertainties, it is the further work that will solve the robust control problem in future.

http://dx.doi.org/10.1155/2014/673256

Conflict of Interests

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

Acknowledgments

This work is partially supported by National Natural Science Foundation of China (61201410, 51379044), Fundamental Research Funds for the Central Universities (HEUCFX41304).

References

 K. Y. Wichlund, Q. J. Sordalen, and O. Egeland, "Control of vehicles with second holonomic constraints: underactuated vehicles," in Proceedings of the European Control Conference, pp. 3086-3091, 1995.

 R. W. Brockett, "Asymptotic stability and feedback stabilization," in Proceedings of Conference Progress in Mathematics, vol. 27, pp. 181-208, 1983.

 I. Fantoni, R. Lozano, F. Mazenc, and K. Y. Pettersen, "Stabilization of a nonlinear underactuated hovercraft," International Journal of Robust and Nonlinear Control, vol. 10, no. 8, pp. 645-654, 2000.

 F. Mazenc, K. Y. Pettersen, and H. Nijmeijer, "Global uniform asymptotic stabilization of an underactuated surface vessel," in Proceedings of the 41st IEEE Conference on Decision and Control, pp. 510-515, December 2002.

 K. Y. Pettersen and H. Nijmeijer, "Global practical stabilization and tracking for an underactuated ship: a combined averaging and backstepping approach," in Proceedings of the IFAC Conference on Systems Structure and Control (SSC '98), pp. 59-64, Nantes, France, July 1998.

 K. Y. Pettersen and O. Egeland, "Exponential stabilization of an underactuated surface vessel," in Proceedings of the 35th IEEE Conference on Decision and Control, vol. 1, pp. 967-972, Kobe, Japan, December 1996.

 T.-H. Kim and T. Basar, "Asymptotic stabilization of an underactuated surface vessel via logic-based control," in Proceedings of the American Control Conference, vol. 6, pp. 4678-4683, May 2002.

 R. Yu, Q. Zhu, G. Xia, and Z. Liu, "Sliding mode tracking control of an underactuated surface vessel," IET Control Theory & Applications, vol. 6, no. 3, pp. 461-466, 2012.

 M. Reyhanoglu, "Exponential stabilization of an underactuated autonomous surface vessel," Automatica, vol. 33, no. 12, pp. 2249-2254, 1997.

 J. Ghommam, F. Mnif, A. Benali, and N. Derbel, "Asymptotic backstepping stabilization of an underactuated surface vessel," IEEE Transactions on Control Systems Technology, vol. 14, no. 6, pp. 1150-1157, 2006.

 Z. Liu, R. Yu, and Q. Zhu, "Comments on asymptotic back-stepping stabilization of an underactuated surface vessel," IEEE Transactions on Control Systems Technology, vol. 20, no. 1, pp. 286-288, 2012.

 J. Cheng, J. Yi, and D. Zhao, "Stabilization of an underactuated surface vessel via discontinuous control," in Proceedings of the American Control Conference, pp. 206-211, July 2007.

 G. Jawhar, M. Faical, and K. A. El-Metwally, "Backstepping technique for tracking control of an underactuated surface

vessel with unmeasured thrust," in Proceedings of the IFAC, vol. 17, 2008.

 D. Q. Mayne, J. B. Rawlings, C. V. Rao, and P. O. M. Scokaert, "Constrained model predictive control: stability and optimality," Automatica, vol. 36, no. 6, pp. 789-814, 2000.

 S. J. Qin and T. A. Badgwell, "A survey of industrial model predictive control technology," Control Engineering Practice, vol. 11, no. 7, pp. 733-764, 2003.

 A. Wahl and E. D. Gilles, "Track-keeping on waterways using model predictive control," in Proceedings of the Workshop on Control Applications in Marine Systems, pp. 149-154, Elsevier, Oxford, UK, 1998.

 T. Perez, Ship Motion Control: Course Keeping and Roll Stabilisation Using Rudder and Fins, Springer, Berlin, Germany, 2005.

 L. Liu, Z. Liu, and J. Zhang, "Nonlinear model predictive control with terminal invariant manifolds for stabilization of underactuated surface vessel," Abstract and Applied Analysis, vol. 2013, Article ID 846389, 8 pages, 2013.

 L. Zhen, Path following with roll constraints for marine surface vessels in wave fields [Ph.D. thesis], The University of Michigan, 2009.

 S. Oh and J. Sun, "Path following of underactuated marine surface vessels using line-of-sight based model predictive control," Ocean Engineering, vol. 37, no. 2-3, pp. 289-295, 2010.

 T. I. Fossen, Marine Control Systems, Marine Cybernetics, Trondheim, Norway, 2002.

 F Mazenc, K. Y. Pettersen, and H. Nijmeijer, "Global uniform asymptotic stabilization of an underactuated surface vessel," in Proceedings of the 41st IEEE Conference on Decision and Control, pp. 510-515, Las Vegas, Nev, USA, December 2002.

 K. Y. Wichlund, O. J. Sordalen, and O. Egeland, "Control properties of underactuated vehicles," in Proceedings of the IEEE International Conference on Robotics and Automation, vol. 2, pp. 2009-2014, Nagoya, Japan, May 1995.

 C. I. Byrnes and A. Isidori, "On the attitude stabilization of rigid spacecraft," Automatica, vol. 27, no. 1, pp. 87-95, 1991.

 K. Y. Pettersen and O. Egeland, "Exponential stabilization of an underactuated surface vessel," in Proceedings of the 35th IEEE Conference on Decision and Control, pp. 967-972, Kobe, Japan, December 1996.

 M. V. Kothare, V. Balakrishnan, and M. Morari, "Robust constrained model predictive control using linear matrix inequalities," Automatica, vol. 32, no. 10, pp. 1361-1379, 1996.

 J. R. Cloutier, C. N. D' Souza, and C. P. Mracek, "Nonlinear regulation and nonlinear HM control via the state-dependent Riccati equation technique," in Proceedings of the International Conference on Nonlinear Problems in Aviation and Aerospace, University Press, Embry-Riddle Aeronautical University, Daytona Beach, Fla, USA, May 1996.

 J. R. Cloutier, "State-dependent Riccati equation techniques: an overview," in Proceedings of the American Control Conference, pp. 932-936, June 1997.

 S. Boyd, L. El Ghaoui, E. Feron, and V. Balakrishnan, Linear Matrix Inequalities in System and Control Theory, SIAM, Philadelphia, Pa, USA, 1994.

 T. Yang and H. R. Karimi, "LMI-based model predictive control for a class of constrained uncertain fuzzy Markov jump systems," Mathematical Problems in Engineering, vol. 2013, Article ID 963089, 13 pages, 2013.

 H. R. Karimi, "Adaptive HM synchronization of master-slave systems with mixed time-varying delays and nonlinear perturbations: an LMI approach," International Journal of Automation and Computing, vol. 8, no. 4, pp. 381-390, 2011.

 H. R. Karimi, "Robust HM filter design for uncertain linear systems over network with network-induced delays and output quantization," Modeling, Identification and Control, vol. 30, no. 1, pp. 27-37, 2009.

 Y. Yin, Y. Liu, and H. R. Karimi, "A simplified predictive control of constrained Markov jump system with mixed uncertainties," Abstract and Applied Analysis, vol. 2014, Article ID 475808, 7 pages, 2014.

 S. Yin, H. Luo, and S. Ding, "Real-time implementation of fault-tolerant control systems with performance optimization," IEEE Transactions on Industrial Electronics, vol. 64, no. 5, pp. 2402-2411, 2014.

 S. Yin, G. Wang, and H. Karimi, "Data-driven design of robust fault detection system for wind turbines," Mechatronics, vol. 24, no. 4, pp. 298-306, 2014.

 S. Yin, S. X. Ding, A. H. A. Sari, and H. Hao, "Data-driven monitoring for stochastic systems and its application on batch process," International Journal of Systems Science, vol. 44, no. 7, pp. 1366-1376, 2013.

 S. Yin, S. X. Ding, A. Haghani, H. Hao, and P. Zhang, "A comparison study of basic data-driven fault diagnosis and process monitoring methods on the benchmark Tennessee Eastman process," Journal of Process Control, vol. 22, no. 9, pp. 1567-1581, 2012.

 S. Yin, X. Yang, and H. R. Karimi, "Data-driven adaptive observer for fault diagnosis," Mathematical Problems in Engineering, vol. 2012, Article ID 832836, 21 pages, 2012.

Lutao Liu, (1) Zhilin Liu, (2) and Jun Zhang (3)

(1) College of Information and Telecommunication, Harbin Engineering University, Harbin, Heilongjiang 150001, China

(2) College of Automation, Harbin Engineering University, Harbin, Heilongjiang 150001, China

(3) School of Electrical and Information Engineering, Jiangsu University, Zhenjiang, Jiangsu 212013, China

Correspondence should be addressed to Zhilin Liu; liuzhilin@hrbeu.edu.cn

Received 23 April 2014; Revised 27 June 2014; Accepted 27 June 2014; Published 14 July 2014

Title Annotation: Printer friendly Cite/link Email Feedback Research Article Liu, Lutao; Liu, Zhilin; Zhang, Jun Abstract and Applied Analysis Report Jan 1, 2014 5225 Nonlinear pulse vaccination in an SIR epidemic model with resource limitation. Dynamics of a predator-prey system with a mate-finding Allee effect. Control systems Inequalities (Mathematics) Mathematical models Ships