Research Article | Open Access

# A Symbolic Computation Approach to Parameterizing Controller for Polynomial Hamiltonian Systems

**Academic Editor:**Yingwei Zhang

#### Abstract

This paper considers controller parameterization method of control for polynomial Hamiltonian systems (PHSs), which involves internal stability and external disturbance attenuation. The aims of this paper are to design a controller with parameters to insure that the systems are stable and propose an algorithm for solving parameters of the controller with symbolic computation. The proposed parameterization method avoids solving Hamilton-Jacobi-Isaacs equations, and thus the obtained controllers with parameters are relatively simple in form and easy in operation. Simulation with a numerical example shows that the controller is effective as it can optimize control by adjusting parameters. All these results are expected to be of use in the study of control for nonlinear systems with perturbations.

#### 1. Introduction

A lot of control problems can be modeled as or transformed approximately to polynomial control systems, which are of great significance in the nonlinear theory. In recent years, a series of achievements have been obtained via symbolic computation and engineering applications [1–3]. The symbolic computation is an efficient and strong tool for polynomial control systems. Meanwhile, generalized Hamiltonian systems are such important nonlinear systems that they have been widely applied in engineering, physical and life science, and so forth. Many scholars studied generalized Hamiltonian systems in recent decades [4–6]. In control design, generalized Hamiltonian systems have good structural properties, which are clear physical expression, good structure and the Hamiltonian function in such a system can be used as good candidate of Lyapunov function. So we consider the polynomial Hamiltonian systems through combining their advantages for control theory. However, the current studies of polynomial Hamiltonian systems (PHSs) merely remain in mathematical theory [7, 8], and robust control and controller parameterization are hardly mentioned. This paper is concerned with control problems of such systems.

As we all know, the internal stability and the external disturbance attenuation are the basic constraints of control systems. However, it is critical to satisfy some desired control objectives in designing a practical control system. It is a sophisticated and efficient way to find the parameterized controller to solve complex control problems, namely, the control, which can satisfy the additional design objectives. Therefore, controller parameterization is a fundamental problem in the control theory and has aroused considerable attention in recent decades [9–19]. Lu and Doyle [10] and Isidori and Astolfi [11] proposed a family of nonlinear controllers via output feedback. Astolfi [12] presented a family of nonlinear state-feedback controllers, in which the system state and the external disturbance are measurable. Yung et al. [13, 14] extended the state-space formulas and presented a family of state-feedback controllers for -dimensional nonlinear system. Fu et al. [15, 16] proposed a family of reliable nonlinear controllers via solving the Hamilton-Jacobi-Isaacs (HJI) inequality (or equations). Xu and Hou [17, 18] studied the generalized Hamiltonian system and proposed a family of parameterized controllers in control and adaptive control. Furthermore, the controllers have been applied to synchronous generators with steam valve in [19].

The controllers obtained in [10–16] are intended to solve a class of HJI inequalities (or equations), which have actually imposed a considerable difficulty. To obtain the parameterized controllers, [17–19] avoided solving HJI inequalities (or equations) applying good structure and clear physical expression of general Hamiltonian systems. However, there are so many restrictions in the process of obtaining a controller in which the form would be very complex even it could not be achieved.

In order to avoid solving HJI inequalities (or equations) and obtain a parameterized controller, this paper presents a novel, straightforward, and convenient method to parameterizing controller for PHSs and gives an algorithm for solving parameters of the controller by symbolic computation. The main results of the paper are as follows.(1)A controller with parameters is obtained via the approach to parameterizing controller for PHSs. Because the Hamiltonian function can be used to build a Lyapunov function of a dissipative Hamiltonian system, the approach obtained here avoids solving HJI inequalities (or equations). Moreover, the obtained controller with parameters has a simpler form than the controller obtained in [17].(2)Parameters ranges are obtained via an algorithm based on symbolic computation. The controller with parameters is effective to control and it can optimize control performance of systems by adjusting parameters.

The remainder of this paper is organized as follows. In Section 2, the problem of control for PHSs is formulated. The main contribution of this paper is then given in Section 3, in which a controller with parameters and an algorithm for solving parameters are provided, respectively. We present a numerical example for illustrating effectiveness and feasibility of controller in Section 4 and conclusions follow in Section 5.

#### 2. Problem Formulation and Preliminaries

Consider the following polynomial Hamiltonian system with dissipation [20]: where is the state vector; is the controller with parameters; is the disturbances; is a skew-symmetric matrix; is a positive semidefinite matrix; and are sufficiently smooth functions; is the output; is the penalty; is a weighting matrix; is the Hamiltonian function which has a local minimum at the equilibrium of system, and . The Hamiltonian function must be a positive semidefinite polynomial and the following assumption holds.

*Assumption 1. * and the Hessian matrix .

*Remark 2. *Note that has a local minimum at the equilibrium of system (1). It is straightforward that in Assumption 1, guarantees the existence of and guarantees that is strict convex on some neighborhood of equilibrium .

The problem considered in this paper is to propose an approach to parameterizing controller for systems (1), which can be described as follows: given a disturbance attenuation level , we can obtain a controller with parameters of the form such that the gain of the closed-loop system (from to ) is bounded by , and the closed-loop system is locally asymptotically stable when .

In the end, we give a definition and a lemma required in the next section.

*Definition 3 (see [20]). *System (1) is called zero-energy-gradient (ZEG) observable with respect to if and , , implies , ; system (1) is called ZEG detectable with respect to if and , , implies ; system (1) is called generalized ZEG observable (detectable) if , , , , implies , ().

Lemma 4 (see [21]). *Consider a nonlinear system
**
where is the state vector, is the disturbances, and is the penalty. If there exists the function such that HJI inequality
**
holds, it is implied that the gain of the closed-loop system (2) (from to ) is bounded by (); that is,
*

#### 3. Main Results

In this section, we propose an controller with parameters for system (1) and an algorithm for solving parameters. The parameterization methods suggest a framework to solve the control problem of polynomial Hamiltonian systems.

##### 3.1. Parameterizing Controller

Theorem 5. *Suppose Assumption 1 holds and also suppose system (1) is generalized ZEG detectable (when ) and
**
hold simultaneously (i.e., the matrix is positive semidefinite and is negative semidefinite). Then, control of system (1) can be realized by the following controller:
**
where is item containing parameters, which is part of the controller parameterization, and is an unit matrix.*

*Proof. *Consider the candidate Lyapunov function (). From Lemma 4 and controller (7), we have
which implies that the gain of the closed-loop system (1) controlled by controller (7) (from to ) is bounded by . Next, we prove that the closed-loop system is asymptotically stable at , when . From system (1), controller (7), and , it follows that

Hence, the closed-loop system converges to the largest invariant set, which is contained in

From Assumption 1 and the fact that system (1) is generalized ZEG detectable, we obtain

On the other hand, since has a local minimum at the equilibrium , it shows . From Assumption 1, there exists a neighborhood of , in which implies . Therefore, from (11), it follows that the trajectory of is strict convex on the equilibrium . According to the LaSalle invariant principle, the closed-loop system (1) controlled by controller (7) is locally asymptotically stable at . This completes the proof.

*Remark 6. *(1) As compared with the controller proposed in [13–19], the controller (7) has a much simpler form and is easier to realize.

(2) Parameters ranges of polynomial vector can be obtained by solving condition (6).

(3) Wang et al. [22, 23] presented the dissipative Hamiltonian realization method for general nonlinear systems that can be transformed to PHSs. So the controller obtained in this paper can be applied to general nonlinear systems.

##### 3.2. Solving Parameters Algorithm

From condition (5), when , we can obtain the . Let such that condition (5) holds. Then we propose an algorithm to find parameters ranges of controller (7) via solving the parameters of in condition (6). The algorithm now proceeds as follows. S1. Set and suppose a positive integer , which is the degree of polynomial vector . Write , where and is the number of state variable. S2. Let . S3. The influence of high order items can be ignored because this paper considers locally asymptotically stable for system. Choose all terms of and from and let the coefficients of these terms be zero. So obtain a set of equations . S3.1. Observe equations . When the right-hand side is only one item with parameters and the left-hand side is zero, let these parameters be zero and substitute them into . Then obtain simplified equations . S3.2. Obtain a set of parameters solution via solving by using cylindrical algebraic decompositions (CAD) algorithm [24]. S3.3. Substitute into and obtain a new polynomial , which is a quadratic form. S4. Rewrite as coefficient matrix , and all principal minors of must be positive semidefinite [25]. Choose all principal minors of and obtain inequalities . S4.1. Observe inequalities . Let some parameters be zero and substitute them into . Then obtain the simplified inequalities . S4.2. Obtain a set of parameters solution via solving by using CAD algorithm. S5. Let and substitute into controller (7), and obtain the polynomial parameterized controller. This completes the algorithm.

*Remark 7. *(1) The algorithm starts from normally.

(2) The CAD algorithm is given by Semi-Algebraic-Set-Tools of Regular-Chains in Maple 16.

(3) It is merely to simplify computation that we let some parameters be zero before using CAD algorithm. However, these parameters are not necessarily zero. So the set of parameters solution obtained by the algorithm is a subset of solutions.

#### 4. Numerical Experiments

Consider a polynomial Hamiltonian system with dissipation (1)

##### 4.1. Controller Design and Solving Parameters

From system (12), it is easy to get So Assumption 1 holds.

Then, we check that condition (5) holds for all and given . From system (12), we have Let . To ensure that condition (5) holds, the following statement should be satisfied

Next, we consider condition (6) such that system (12) satisfies robustness in control. Set the candidate Lyapunov function . It follows from controller (7) that where .

We know in system (12) and let . We have where are the parameters.

From system (12), we obtain that .

Let ; we have is a quadratic form and can be rewritten as a coefficient matrix (multiply constant 2 for simplifying computation),

All principal minors of must be positive semidefinite. We have inequalities from . From , we can easy obtain that , and . Substitute into inequalities to simplify computation; we obtain simplified inequalities . Solving inequalities by using CAD algorithm, we obtain a series of sets. Choose some sets, which satisfy inequalities , and organize them. We have Substitute into controller (7), where , , . So we have the controller with parameters for system (12). The controller (21) has a rather simple form.

##### 4.2. Simulations and Results

In order to evaluate the robustness of the controller (21), we set the parameters of system (12) as , , and the parameters of controller as , , . We obtain the following controller:

Suppose that is the preassigned operating point of system (12); we impose an external disturbance on system (12) during the time period 2~2.2 s. In order to evaluate the robust performance of controller (22), we compare it with the saturated controller proposed in [26]. The parameters of the saturated controller are chosen as follows: , , and , , .

The simulation results are shown in Figures 1, 2, 3, and 4, where is the response of the state under controller (22) and is the response of the state under saturated controller in [26]. From Figures 1–4, we can clearly see that under controller (22), it takes only 6.8 seconds for the system to return back to the equilibrium point (circle point in figures), while under saturated controller, it takes about 24 seconds. The simulation shows that controller (22) is much more efficient and feasible and has stronger robustness than saturated controller.

In order to evaluate the robustness optimization of the system by adjusting the parameters of controller (21), we choose the parameters of system (12) as , and the parameters of controller as . We have the following controller:

Suppose that is the preassigned operating point of system (12); we impose an external disturbance on system (12) during the time period 2~2.2 s. In Figure 5, is the response of the state under controller (22), and in Figure 6, is the response of the state under controller (23).

From Figure 5, we can clearly see that under controller (22), it takes 6.8 seconds for the system to return back to the equilibrium point (circle point in figures), while under controller (23), it takes about 4.8 seconds in Figure 6. The time of state vectors and back to the equilibrium point is much shorter in Figure 5 than in Figure 6. The simulation shows that controller (21) can effectively optimize the robustness of system by adjusting parameters.

#### 5. Conclusion

In this paper, an approach to parameterizing controller for polynomial Hamiltonian systems has been considered. A controller with parameters has been obtained using Hamiltonian function method and an algorithm for solving parameters of the controller has been proposed with symbolic computation. The proposed parameterization method avoids solving Hamilton-Jacobi-Isaacs equations and thus the obtained controllers with parameters are easier as compared to some existing ones. The numerical experiment and simulations show that the controller is efficient in control; meanwhile the controller can optimize the robustness of the system by adjusting parameters. The approach proposed in this paper is also applicable to general nonlinear systems, which should be transformed to the polynomial Hamiltonian systems through the available dissipative Hamiltonian realization methods.

#### Conflict of Interests

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

#### Acknowledgment

This work is supported by the National Natural Science Foundation of China (nos. 61074189 and 61374001).

#### References

- Y. Huang, Y. Shang, and W. Yuan, “Symbolic computation and the extended hyperbolic function method for constructing exact traveling solutions of nonlinear PDEs,”
*Journal of Applied Mathematics*, vol. 2012, Article ID 716719, 19 pages, 2012. View at: Publisher Site | Google Scholar | MathSciNet - J. F. Navarro, “The asymptotic expansion method via symbolic computation,”
*Journal of Applied Mathematics*, vol. 2012, Article ID 201340, 24 pages, 2012. View at: Publisher Site | Google Scholar | MathSciNet - H. Deng and J. Wu, “Approximate bisimulation and optimization of software programs based on symbolic-numeric computation,”
*Mathematical Problems in Engineering*, vol. 2013, Article ID 421926, 19 pages, 2013. View at: Publisher Site | Google Scholar | MathSciNet - R. Ortega, A. van der Schaft, B. Maschke, and G. Escobar, “Interconnection and damping assignment passivity-based control of port-controlled Hamiltonian systems,”
*Automatica*, vol. 38, no. 4, pp. 585–596, 2002. View at: Publisher Site | Google Scholar | MathSciNet - Y. Wang, S. S. Ge, and D. Cheng, “Observer and observer-based ${H}_{\infty}$ control of generalized Hamiltonian systems,”
*Science in China F*, vol. 48, no. 2, pp. 211–224, 2005. View at: Publisher Site | Google Scholar | MathSciNet - I. G. Polushin, “Stabilization of invariant tori in Hamiltonian systems under persistently acting disturbances,”
*International Journal of Robust and Nonlinear Control*, vol. 11, no. 3, pp. 253–265, 2001. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - H. Zhou, W. Xu, S. Li, and Y. Zhang, “On the number of limit cycles of a cubic polynomials Hamiltonian system under quintic perturbation,”
*Applied Mathematics and Computation*, vol. 190, no. 1, pp. 490–499, 2007. View at: Publisher Site | Google Scholar | MathSciNet - F. Iavernaro and D. Trigiante, “High-order symmetric schemes for the energy conservation of polynomial Hamiltonian problems,”
*Journal of Numerical Analysis, Industrial and Applied Mathematics*, vol. 4, no. 1-2, pp. 87–101, 2009. View at: Google Scholar | MathSciNet - J. Gadewadikar, F. L. Lewis, L. Xie, V. Kucera, and M. Abu-Khalaf, “Parameterization of all stabilizing H∞ static state-feedback gains: application to output-feedback design,”
*Automatica*, vol. 43, no. 9, pp. 1597–1604, 2007. View at: Publisher Site | Google Scholar - W.-M. Lu and J. C. Doyle, “${\mathcal{H}}_{\infty}$ control of nonlinear systems via output feedback: controller parameterization,”
*IEEE Transactions on Automatic Control*, vol. 39, no. 12, pp. 2517–2521, 1994. View at: Publisher Site | Google Scholar | MathSciNet - A. Isidori and A. Astolfi, “Disturbance attenuation and ${H}_{\infty}$-control via measurement feedback in nonlinear systems,”
*IEEE Transactions on Automatic Control*, vol. 37, no. 9, pp. 1283–1293, 1992. View at: Publisher Site | Google Scholar | MathSciNet - A. Astolfi, “Parameterization of output feedback controller that satisfy an norm bound,” in
*Proceedings of the 2nd European Control Conference*, pp. 74–78, Groningen, The Netherlands, 1993. View at: Google Scholar - C. Yung, Y. Lin, and F. Yeh, “A family of nonlinear H
_{∞}output feedback controllers,”*Institute of Electrical and Electronics Engineers. Transactions on Automatic Control*, vol. 41, no. 2, pp. 232–236, 1996. View at: Publisher Site | Google Scholar | MathSciNet - C. F. Yung, J. L. Wu, and T. T. Lee, “Parameterization of nonlinear ${H}_{\infty}$ state-feedback controllers,”
*Automatica*, vol. 33, no. 8, pp. 1587–1590, 1997. View at: Publisher Site | Google Scholar | MathSciNet - Y. Fu, Z. Tian, and S. Shi, “A family of reliable nonlinear H∞ state-feedback controllers,”
*Control Theory and Applications*, vol. 18, no. 3, pp. 447–452, 2001. View at: Google Scholar | MathSciNet - Y. S. Fu, Z. H. Tian, and S. J. Shi, “A family of reliable state-feedback controllers for nonlinear systems with strictly redundant actuators: the full information case,”
*Journal of Systems Science and Systems Engineering*, vol. 10, no. 3, pp. 350–358, 2001. View at: Google Scholar - S. Xu and X. Hou, “A family of ${H}_{\mathit{\infty}}$ controllers for dissipative Hamiltonian systems,”
*International Journal of Robust and Nonlinear Control*, vol. 22, no. 11, pp. 1258–1269, 2012. View at: Publisher Site | Google Scholar - S. Xu and X. Hou, “A family of adaptive H
_{∞}controllers with full information for dissipative Hamiltonian systems,”*International Journal of Automation and Computing*, vol. 8, no. 2, pp. 209–214, 2011. View at: Publisher Site | Google Scholar - S. Xu and X. Hou, “A family of robust adaptive excitation controllers for synchronous generators with steam valve via Hamiltonian function method,”
*Journal of Control Theory and Applications*, vol. 10, no. 1, pp. 11–18, 2012. View at: Publisher Site | Google Scholar | MathSciNet - A. van der Schaft and A. J. Schaft,
*L2-Gain and Passivity in Nonlinear Control*, Springer, Berlin, Germany, 1999. - T. Shen,
*Control Theory and its Applications*, Tsinghua University Press, Beijing, China, 1996. - Y. Wang, D. Cheng, C. Li, and Y. Ge, “Dissipative Hamiltonian realization and energy-based ${L}_{2}$-disturbance attenuation control of multimachine power systems,”
*IEEE Transactions on Automatic Control*, vol. 48, no. 8, pp. 1428–1433, 2003. View at: Publisher Site | Google Scholar | MathSciNet - Y. Wang, D. Cheng, and S. S. Ge, “Approximate dissipative Hamiltonian realization and construction of local Lyapunov functions,”
*Systems & Control Letters*, vol. 56, no. 2, pp. 141–149, 2007. View at: Publisher Site | Google Scholar | MathSciNet - F. B. Caviness and R. J. Johnson,
*Quantifier Elimination and Cylindrical Algebraic Decomposition*, Springer, New York, NY, USA, 1998. - K. N. Swamy, “Sylvester’s criterion for positive-semidefinite matrices,”
*IEEE Transactions on Automatic Control*, vol. 18, no. 3, p. 306, 1973. View at: Google Scholar - Z. Xi, G. Feng, D. Cheng, and Q. Lu, “Nonlinear decentralized saturated controller design for power systems,”
*IEEE Transactions on Control Systems Technology*, vol. 11, no. 4, pp. 539–547, 2003. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2014 Zhong Cao and Xiaorong Hou. 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.