Research Article | Open Access

# Suboptimal Control Using Model Order Reduction

**Academic Editor:**X. Ma

#### Abstract

A Padé approximation based technique for designing a suboptimal controller is presented. The technique uses matching of both time moments and Markov parameters for model order reduction. In this method, the suboptimal controller is first derived for reduced order model and then implemented for higher order plant by partial feedback of measurable states.

#### 1. Introduction

The design of an optimal control law for large scale systems is inconvenient due to one of the main drawbacks of optimal control theory that requires feedback from all state variables that are defined to describe the dynamics of the plant. It is rare that all the states are available for measurement [1] and if one tries to design optimal controller, Kalman filter or state observer is required to reconstruct the missing state of the plant. This introduces the high order dynamics in the control function and leads to a complicated and costly controller. In order to minimize this problem, an attempt has been made for the design of incomplete state feedback suboptimal control law using measurable states [2–4]. In this paper a method for suboptimal controller design is proposed. The suboptimal control can be achieved by using model order reduction technique for higher order system.

The plant (higher order system) can be reduced by Padé approximation technique (matching of both time moments and Markov parameters) into reduced order model, and then, with the help of reduced order controller, one can design the suboptimal controller for plant (higher order systems). The usefulness of techniques for deriving reduced order approximations of high order system has been already accepted due to the advantages of reduced computational effort and increased understanding of original system [5, 6]. Consequently, a large number of time-domain (Hankel norm, balance truncation, and Krylov subspace projection) and frequency-domain (Moment matching, Padé approximation, Routh-Padé approximation, and pole clustering methods) system simplification techniques have been developed to suit the different requirements and a variety of numerical approximation are applied to obtain the reduced order model [7]. In this paper the plant (higher order system) is reduced by Padé approximation technique and suboptimal controller is derived with the reduced order model. Then the suboptimal controller is implemented for original plant (higher order system). This reduces the computational complexity of controller design. Two numerical examples are included to illustrate the procedure.

#### 2. The Design Procedure

Consider a single-input single-output linear (th-order) dynamic system described by with the quadratic cost function where is the positive semidefinite matrix and is a positive weight. , , and are matrices of appropriate dimensions. It is well know that the optimal feedback control law is linear combination of the state variables. A block diagram of the plant using control law has been shown in Figure 1.

The control action of the plant is given as where is the symmetric positive definite matrix whose elements may be obtained by solving the following matrix Riccati equation: The closed loop transfer function with controller is obtained as

##### 2.1. Model Order Reduction Procedure by Padé Approximation

Consider a th order system [5] as It is desired to obtain a stable th order approximant as It is easy to verify that the following equations hold true: and for ; ; for ; We seek a reduced order model so as to satisfy where .

Using (8)–(16) the reduced order model is obtained.

##### 2.2. Suboptimal Controller Design Procedure

The suboptimal controller [4] is designed with the help of reduced order model by choosing the optimal gain matrix in such a way that where is optimal gain matrix of reduced order model, obtained using (3) and (4). Finally this is applied to higher order system in order to control the plant (higher order system).

##### 2.3. Flow Diagram of Suboptimal Controller Design

See Figure 2.

##### 2.4. Comparison of Responses of Optimal and Suboptimal Controller

Using (5) one can try to determine the closed loop transfer functions with optimal controller and suboptimal controller with the gain matrices and , respectively being derived for comparison.

#### 3. Numerical Examples

The step-by-step procedure to obtain reduced the order model is explained with the help of the example presented below and procedure is derived for suboptimal controller.

*Example 1. *Consider the following 4th order system [3]:
Suppose that reduced order model of the form
is required.

*Reduction Steps*. The approximant can systematically be derived in the following steps for second order approximant.

*Step 1.* Expand (12) in power series around and :

*Step 2.* Expand (13) in power series around and :
The time moment and Markov parameter expansion can obtain about and .

*Step 3.* For perfect matching we have to match the () term:

*Step 4.* By considering the above assumptions we can get the second order reduced model by solving the following and get the variables of reduced order model (13),
But for better overall time response approximation, equal emphasis is given to both time moments and Markov parameters. For second order reduced model, we have to match terms and satisfy ; here and consider .

By solving (17) we get and require the reduced order model obtained as The step response of higher order system (12) and reduced order model (19) is plotted as shown in Figure 3.

Model (19) is seen to be a closed approximant to (12).

The optimal controller and suboptimal controller can be designed for a particular performance index by solving Riccati equation (4).

The optimal controller for actual higher order system (12) is obtained by using (3)-(4) as and the transfer function with optimal controller is obtained by putting the value of gain matrix (20) in (5) as Similarly, the optimal gain matrix for reduced order model (19) by using (3) and (4) is as follows: The closed loop transfer function of reduced order model (19) is obtained by putting (22) in (5) as The step response of higher order system with optimal controller (21) and reduced order model with optimal controller (23) is shown in Figure 4.

It is found that (23) is a closed approximation to (21).

The suboptimal controller (11) is obtained with the help of reduced order model using .

Thus takes the following form: The transfer function with suboptimal controller is obtained by putting the value of (24) in (5) as The step response of higher order system with optimal control (21) and with suboptimal controller (25) is shown in Figure 5.

The step responses of higher order system with optimal controller (21) and with suboptimal controller (25) are shown in Figure 5. The response of suboptimal controller (25) is seen identical to (21). This shows that we can control the plant with the help of only measurable states which reduces the complexity and cost.

*Example 2. *Consider the higher order system [6]
The optimal controller for (26) is obtained using (3) and (4) as
The closed loop transfer function using (5) is expressed as
It is assumed that one state is not available for measurement. Thus a reduced order model is desired which takes the following form:
And the optimal gain matrix for (29) is obtained with the help of (3) and (4) as
The closed loop transfer function of (29) with optimal controller is obtained by putting (30) into (5) as
The comparison of the result of step responses of higher order system with optimal controller (28) and reduced order model with optimal controller (31) is shown in Figure 6.

It is clear that the step response of higher order system with optimal controller (28) and reduced order model with optimal controller (31) is very close to identical.

The suboptimal gain matrix is obtained for the system (26) with the help of (11) and (30) and is expressed as

The closed loop transfer function of higher order system (26) with suboptimal controller is obtained using (32) and (5) as
The step responses of (28) and (33) are shown in Figure 7 and are found identical in transient as well as steady state region.

The results of the two examples signify that the suboptimal control provides a good platform for controlling the higher order system and reduces the complexity and cost of control action.

For purpose of comparison, the suboptimal controller is also derived using the technique in [4]. The result obtained in [4] shows remarkable deviation in the transient part which confirms the observation of [4] that, for better overall time response approximation, matching of both time moments and Markov parameters is required.

#### 4. Conclusion

A design technique of suboptimal controller using partial feedback is proposed. The method is based on the Padé approximation technique (matching of both time moments and Markov parameters) for model order reduction which results in better overall time response approximation. In this method, the suboptimal controller is first derived for reduced order model and then implemented for higher order plant. This reduces the design complexity and cost of control action. The method also avoids use of observer or Kalman filter for reconstruction of missing states.

In this paper equal emphasis is given for matching of time moments and Markov parameters; however, the optimal choice of Markov parameters and time moments to be matched remains unresolved and it is open to further investigation.

#### Conflict of Interests

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

#### References

- G. Bengtsson and S. Lindahl, “A design scheme for incomplete state or output feedback with applications to boiler and power system control,”
*Automatica*, vol. 10, no. 1, pp. 15–30, 1974. View at: Google Scholar - T. C. Hsia, “An approach for incomplete state feedback control systems design,”
*IEEE Transactions on Automatic Control*, vol. AC-17, pp. 383–386, 1972. View at: Google Scholar - V. P. Singh, P. Chaubey, and D. Chandra, “Model order reduction of continuous time systems using pole clustering and Chebyshev polynomials,” in
*Proceedings of the IEEE Students Conference on Engineering and Systems (SCES '12)*, 2012. View at: Google Scholar - J. Pal, “Suboptimal control using Padé approximation technique,”
*IEEE Transactions on Automatic Control*, vol. AC-25, no. 5, 1980. View at: Google Scholar - V. Singh, D. Chandra, and H. Kar, “Improved Routh-Padé approximants: a computer-aided approach,”
*IEEE Transactions on Automatic Control*, vol. 49, no. 2, pp. 292–296, 2004. View at: Publisher Site | Google Scholar - D. Chandra,
*An investigation into Routh’s Padé approximants [Ph.D. thesis]*, Department of Electrical Engineering, Motilal Nehru National Institute of Technology, Allahabad, India. - D. S. Janardhanan, “Model Order Reduction and Controller Design,” http://web.iitd.ac.in/~janas/courses/material/eel879/sp_topics_01.pdf. View at: Google Scholar

#### Copyright

Copyright © 2014 Avadh Pati et al. 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.