Mathematical Problems in Engineering

Volume 2016, Article ID 6924324, 10 pages

http://dx.doi.org/10.1155/2016/6924324

## Optimal Preview Control for Discrete-Time Systems in Multirate Output Sampling

School of Mathematics and Physics, University of Science and Technology Beijing, Beijing 100083, China

Received 9 August 2016; Accepted 22 September 2016

Academic Editor: Michele Brun

Copyright © 2016 Fucheng Liao and Yujian Guo. 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.

#### Abstract

This paper studies the disturbance preview optimal control problem for discrete-time systems with multirate output sampling. By constructing the error system and using the discrete lifting technique, we reduce the multirate preview control problem to a single-rate one for a formal augmented system. Then, applying preview control theory, the optimal preview control law of the augmented error system is obtained. Meanwhile, we introduce a discrete integrator to eliminate the static error. Then we study a method to design a controller with preview action for the original system. And the existence conditions of the controller are also discussed in detail. Finally, numerical simulation is included to illustrate the effectiveness of the proposed method.

#### 1. Introduction

In multirate output control systems, each component of the input vectors only changes once, and each component of the output vectors is detected several times during one sampling period. The increase in sampling frequency of the output vectors allows the controllers to obtain more information on the controlled objects and to acquire stronger control ability [1]. Reference [2] studied a “multirate output controller” and designed a multirate output feedback controller for a class of time-invariant continuous-time systems. References [3, 4] studied the design of the multirate output sampling functional observer and periodic output feedback controller based on fast output sampling, respectively. The above two papers take full advantage of the multirate output control systems’ ability to increase the number of detecting output vectors.

Optimal control theory has numerous applications in both science and engineering. Recent reference [5] develops a comprehensive optimality theory of problems described by ordinary and partial differential inclusions and presents a number of new and important results. The theoretical basis of preview control is optimal control. Preview control theory takes full advantage of the known future reference signals or disturbance signals to improve the performance of control systems. Preview control theory has generated a whole set of theories and methods during its development [6–8]. There are many studies in the literature on this theory, and the perfect one is the study based on linear time-invariant systems [9]. Lately, progress has been made in combining preview control theory with stochastic systems [10], nonlinear systems [11], and descriptor systems [12]. Moreover, preview control theory has been applied successfully to vehicle active suspensions [13], wind turbines [14], and other aspects.

The sampling rate for control systems may be limited by physical hardware. For example, the sampling rate in the control of hard-disk systems is determined by the number of sectors on circular tracks and the spindle speed. A low sampling rate can reduce the cost of hardware resources. However, an overly low sampling rate has a bad effect on system performance and system stability. Generally, the performance of control systems can be improved by increasing sampling frequency properly [15]. In [16], it has been shown that there is an advantage in updating feedforward control input at a higher rate in track-seeking control. These cases motivate the research work for preview control theory of multirate systems.

Research into preview control theory in multirate settings began with [17]. In that paper, Liao et al. presented an LQ optimal preview servomechanism design method for linear discrete-time systems with time-delay in a multirate setting. The paper eliminated the time-delay by state augmentation, similar to the technique offered in [6]. Reference [18] discussed the optimal preview control problem for a class of multirate sampled systems with state time-delay. Reference [19] applied preview control theory to solve the optimal preview control problem for discrete-time descriptor causal systems in a multirate setting. But the above papers all study the preview control problem for multirate input control systems. Compared with the frequent changes of input vectors for multirate input control systems, multirate output control systems make the structures of systems relatively simpler and suitable for industrial process control. So this paper studies optimal preview control for discrete-time systems in multirate output sampling. We will use the following lemmas repeatedly in [20].

Lemma 1 (PBH rank test). * is stabilizable if and only if, for any complex satisfying , the matrix has full row rank. is detectable if and only if, for any complex satisfying , the matrix has full column rank.*

#### 2. Problem Description and Related Assumptions

Consider the following linear discrete-time system:where , , , and represent the state vector, the control input vector, the output vector, and the disturbance vector, respectively. , , , and are known constant matrices with appropriate dimensions.

For system (1), we introduce the following assumptions.

*Assumption 2. * is stabilizable.

*Assumption 3. * is detectable.

*Assumption 4. *The input vector can be input only at , where is a positive integer.

*Remark 5. *Assumptions 2 and 3 are basic assumptions of control systems. Assumption 4 makes the system multirate. That is, the input vector can be input once during every sampling interval and the output vector can be measured at every moment.

*Assumption 6. *, ; .

*Remark 7. *Suppose the system has zero-order-hold; that is, .

*Assumption 8. *The matrixhas full row rank.

*Remark 9. *Assumption 8 is one of the basic assumptions of preview control theory in [6].

We assume that the disturbance signal is previewable.

*Assumption 10. *The future values of the disturbance , as well as the present and past values, are available at time , where is the preview length of the disturbance and , is a nonnegative integer.

Furthermore, we assume the following conditions, for simplicity of discussion.

*Assumption 11. *The disturbance converges to zero as time goes to infinity. That is,

We assume reference signal is given by the step functionAccording to Assumptions 8 and 11, there always exist constant vectors and such thatThat is,For the equation of , the rank of the coefficient matrix is the same as that of the augmented matrix, by Assumption 8. So, this equation has a solution. In general, we take one solution from the above equation and express this solution aswhere denotes the right inverse.

Define the new vectors,and the error signal,From (1) and (5), we get the error system

In this paper, we wish to design a type-1 servo controller such that the output can track reference signal without steady-state error. That is,To achieve good transient response, we introduce the quadratic performance indexwhere the weight matrices satisfy ,.

#### 3. The Derivation of the Augmented Error System

In this section, we eliminate the multirate feature of Assumption 4 using the discrete lifting technique. Using the lifting scheme of this paper, we can construct the necessary augmented error system. Hence the preview control problem for a multirate output sampling system is converted into an optimal regulator problem for the normal augmented error system. And the controller with preview action for the augmented error system can be obtained by applying optimal preview control theory.

##### 3.1. Discrete Lifting for the Multirate Output Sampling System

According to Assumptions 4 and 6, the input vector of system (10) can be input at and ().

Noticing the first equation in (10) and the characteristic of , we getThat is,Similarly, we obtain the equation of the error signal

We introduce the vectors as follows:Combining (14) and (15), we getNow we have successfully transformed the multirate output sampling system into the single-rate system by the lifting technique.

##### 3.2. Introduction of Integral Compensation

To achieve the robust servo property mentioned in the previous section and eliminate the steady-state error, we need to introduce the discrete integrator, defined bywhereEquation (18) is the discrete integrator and we can derive the following equation: That is,Noticing the second equation in (17), we obtainIf the closed-loop system is stabilized, there exists a stationary value of . Let and defineFrom (22), we haveThis enables us to add a degree of freedom to improve the transient response of the closed-loop system since can be arbitrarily assigned [17].

Combining (17) and (24) and matching the corresponding block, we obtain the augmented systemwhereAugmented system (25) is the formal system for designing the controller.

##### 3.3. Modification of the Performance Index Function

From (12) and the structure of and , we can seewhere , , and , , ,.

Noticing the second equation of (17), the performance index (27) can be rewritten as

In order to regularize the LQ optimal control problem and guarantee the existence of a stabilizing solution to the associated Riccati equation, we modify the performance index as follows:where , , , and .

Now, the preview control problem of multirate output control system (1) becomes a standard LQ optimal preview control problem of augmented system (25) with performance index (29).

Applying the results in [6–8], we can derive the following theorem.

Theorem 12. *If is stabilizable and is detectable, the optimal preview controller of system (25) minimizing the performance index (29) is given bywhereand is the unique symmetric semi-positive definite solution of the algebraic Riccati equation,*

##### 3.4. Realization of the System

In the processes of constructing the augmented system, the order of formal system (25) increases highly. As a result, the order of corresponding Riccati equation (32) is very high. We study a method to solve high-order Riccati equation (32) by solving a low-order Riccati equation. Noticing the structure of the coefficient matrices and for system (25) and the weight matrices and for performance index (29), the solution of the Riccati equation (32) can be decomposed intoSubstituting the partitioned matrix of into (32), we deriveIn fact, we can obtain (32) by repeating (34) times. Thus, the solution of high-order Riccati equation (32) is transformed into solving low-order Riccati equation (34).

And this is equivalent to solving (34) for the following system:That is,We know that (34) has a unique symmetric semi-positive definite solution if and only if is detectable. Notice that (25) can be written asand the first equation of (37) is (36). We can obtain by the solution of the Riccati equation (34) and the state vector . Substituting into the first equation of (37), is obtained. Then, substituting and into the second equation of (37), is obtained. By that analogy, is obtained. Then, the simulation can be realized.

#### 4. Existence Conditions of the Controller

It is a well-known fact that the stabilizability of guarantees that the state feedback gain in Theorem 12 exists. And the detectability of ensures that Riccati equation (32) has a unique symmetric semi-positive solution [20].

First, we examine the stabilizability of .

Lemma 13. * is stabilizable if and only if is stabilizable.*

*Proof. *By the PBH rank test,for any complex , if , we haveThus, matrix has full row rank if and only if matrix has full row rank. So, Lemma 13 is established.

Lemma 14. * is stabilizable if and only if ** is stabilizable and the matrix ** has full row rank.**Noticing the structure of ** and **, one gets this lemma in [6]. **Therefore, one obtains the following theorem on the stabilizability of *.

Theorem 15. * is stabilizable if and only if ** is stabilizable and the matrix ** has full row rank.**Theorem 15 guarantees that ** is stabilizable under Assumptions 2 and 8. *

Next, we discuss the detectability of .

Lemma 16. * is detectable if and only if ** is detectable*.*Noticing the structure of ** and *,* one can prove Lemma 16 by using the similar method in Lemma *13*. Here, one omits the proof. *

Lemma 17. * is detectable if and only if is detectable and .*

*Proof. *According to [17], is detectable if and only if is detectable and . And , so the detectability of is equivalent to that of . Then Lemma 17 holds.

From Lemmas and 17, we obtain.

Theorem 18. * is detectable if and only if is detectable and .**Consequently, we see from Theorems 15 and 18 that the solvability of the Riccati equation (32) is guaranteed under the standard servomechanism assumptions for the original system.*

#### 5. Design of an Optimal Preview Controller for the Original System

Now we return to the controller design with preview action for the original system. We can see that Riccati equation (32) can be obtained by repeating (34) times, as pointed out above. Then we consider the optimal input (30) of system (25) and associated (31) and (32) in Theorem 12.

Noticing the structure of , , and and calculating , we haveThis is the block diagonal matrix. Similarly, we notice that , , and are all block diagonal matrices and the submatrix blocks of the diagonal are identical. Calculating these block diagonal matrices and submitting them into (30) we getwhereand is the unique symmetric semi-positive definite solution of algebraic Riccati equation (34).

We decompose into

Equation (41) can be written asThat is,Further, the above equation can be written asFrom (21), we have ; that is,

By substituting (7) into the above equation, we have

Then we obtain the main theorem of this paper.

Theorem 19. *If the following conditions hold:*(1)* is stabilizable (Assumption 2 holds);*(2)* is detectable (Assumption 3 holds);*(3)*the matrix has full row rank (Assumption 8 holds);*(4)*,, .** Riccati equation (32) has a unique symmetric semi-positive definite solution, and the optimal control input of system (1) with multirate output sampling iswhere , , and are determined by (42) and (43) and where*

*, .*

*The closed-loop system of (1) is*

*6. Numerical Example*

*Consider a linear discrete-time system:where the coefficient matrices are*

*Let ; that is, the input vector can be input only at . Let the initial conditions be and . Take the weight matrices , , and .*

*Through verifying, system (52) satisfies Assumptions 2, 3, and 8. In this case, , , and are calculated to beRiccati equation (34) is To solve this, we obtainNow, let the disturbance signal be*

*For the following step reference signal, we perform MATLAB simulation results for the multirate controllers with preview action (let the preview length be , i.e., ) and the one without preview action. *

*The output responses of the linear discrete-time system in multirate output sampling are shown in Figure 1. The error signals are shown in Figure 2. We see from the figures that the preview action restrains the system disturbance effectively. And the simulation result shows better tracking in the case of disturbance preview action. Note that the error signal is asymptotically zero. Furthermore, the vibrations of output response curves are the multirate feature of this system.*