Preview Tracking Control for Discrete-Time Multirate Systems: An Internal Model-Based Approach

: This paper investigates the preview tracking problem for discrete-time multirate systems under the scenario of the external signals with amplitude non-attenuation. Since the input and the state cannot be measured at the same time, the multirate system is ﬁrst lifted into the single-rate system. Noting that the external signals are previewable and their amplitudes are unattenuated, an auxiliary system with preview information and a dynamic compensator with internal model are introduced, respectively. On this basis, the preview tracking problem can be converted to the output regulation problem of an augmented system with the aid of the state augmentation technique. Through verifying the solvability of the output regulation problem, the su ﬃ cient conditions and the corresponding controller with preview action and feedforward compensation are presented for the output regulation problem as well as the original problem. Finally, a numeric simulation result is provided to verify the validity of the proposed method.


Introduction
Multirate control system has attracted considerable attention due to a mount of sample requirements in many fields of industry. The analysis and synthesis for multirate systems originates in [1]. To date, various results have been developed for this issue [2][3][4][5][6][7][8][9][10]. In [2], a discrete-time state space equation was established for multirate systems, where four stability criteria were derived based on two forms of a transfer matrix related to the state space equation. In [3], proper algebraic conditions were provided to ensure the desired structural properties (such as, reachability, controllability and stabilizability, etc.) of multirate systems. If the saturated communication sequences are equidistant, multirate systems and networked control system with limited bandwidth were proved to be equivalent in [4], as a result, the related conclusions of the latter on controllability and stabilizability can be applied to analyze the former. In order to simulate the dynamic characteristics in sampling interval, two different discrete lifting techniques were applied in [5] to obtain the functional space model of the multirate control systems, then based on the definition of zeros and zero directions, the conditions for realizing asymptotic ripple-free tracking of multirate systems are provided.
Note that dual-rate control is suitable for the scenario that output data are sampled by a rate slower than that of control data sampling, [6] presented two adaptive control schemes for a dual-rate discretized continuous-time system. In addition, for dual-rate systems, the parameters of lifted single-rate systems were usually assumed to be known in most literature, [7] investigated the case of unknown parameters and developed an auxiliary model to identify the parameters of the lifted Mathematics 2020, 8,1443 2 of 20 fast single-rate systems. In [8], the controller design for multirate systems were split into two steps, where the slow part was used to compute the control error and the fast part was applied to get the similar results with those achieved by the fast single-rate control. [9] considered the H ∞ control synthesis problem for multirate systems and the solution of which was obtained by formulating that into a convex optimization problem. [10] studied the robust sample-data H ∞ control problem for uncertain active vehicle suspension systems, in light of an input delay method, the original system was converted to a continuous-time system with a state delay. For other results of multirate systems, interested readers can also refer to References [11][12][13].
For many practical control systems, such as wind turbine control [14]; walking pattern generation for biped robots [15], the part (even full) of the external signals are usually known in advance for the system. Preview control theory comes along with the occurrence of such phenomenon, the core issue of which is to design a controller containing the future information to improve the response performance of the closed-loop system. To date, there have been many important advances in the theoretical research in this field. Tomizuka has done many groundbreaking works in this field, such as [16,17]. In [16], a method for designing tracking controller with preview information to minimize the optimality criterion is presented, where the reference signal is subject to Gaussian white noise. Based on the linear quadratic optimal control theory, a digital controller with error integral and preview compensation for discrete-time linear system is designed in [17]. In order to enrich the results of [17], the case of the unmeasurable state is further considered in [18]. It should be pointed out that if the reference signal incorporates the white noise or other random disturbances, the linear quadratic Gaussian (LQG) control technique will not be suitable for designing the optimal controller [19]. In addition, if the reference signal is treated as an external perturbation and the tracking problem is expressed as a H ∞ optimal control problem, the effect of the preview information will not be reflected. To this end, game theorem is adopted in [20] to study the finite time H ∞ tracking control problem. Furthermore, the discrete-time counterpart of [20] was investigated in [21]. Other relevant results for the H ∞ preview control can be also found in [22,23]. Moreover, by building an auxiliary system with preview information and augmenting it with the original system, the state feedback and the preview feedforward compensation in the controller can be designed in a unified way. Moreover, relevant results along this line can be referred to [24][25][26][27].
Currently, the preview control for multirate systems has also been extensively studied. In [28,29], the multirate system subject to input delay and state delay were transformed into a single-rate system without time-delay by discrete lifting technique, which simplified the problem significantly. Recently, based on the above results, optimal preview control problem for descriptor multirate systems was solved in [30,31]. Moreover, the cases of dual-rate systems and uncertain multirate systems have also been considered in [32,33], respectively. Note that, in the theoretical research of the preview control for multirate system, it is assumed that the reference signal tends to be a constant vector (or scalar) as time approaches infinity. However, many reference signals with unattenuated amplitudes, for instant sinusoidal signal, do not meet such assumption. Recently, using the distributed internal model principle, the cooperative preview tracking problem of multi-agent system with undamped amplitude was investigated in [34].
Inspired by the research method of [34], this study intends to consider the preview control problem of multirate systems with non-decreased amplitude reference signal by introducing the internal model compensator. Consequently, the preview tracking problem is processed via transforming it into an output regulation problem of an augmented system. The innovations are summarized as follows: (1) The preview control problem of multirate systems is further improved by considering the scenario of amplitude non-decaying reference signal; (2) The stabilizability and detectability of the augmented systems are proved strictly, which guarantees the existence and uniqueness of solutions of Riccati equation and Sylvester equation; (3) The asymptotical tracking conditions for existing an multirate controller with error integral, internal model compensation and preview feedforward compensation are presented.
Notations: Let R n and R n×m be the n-dimensional Euclidean space and the n × m matrix space, respectively. I n denotes the n × n identity matrix. Λ(S) represents the spectrum of the matrix S, i.e., λ det(λI − S) = 0 . 'iff' means 'if and only if'. diag{· · · } stands for a diagonal matrix. Q > 0 denotes that Q is a symmetric positive definite matrix. '∪' means the union of different sets; 'A ⊗ B' represents the Kronecker product of matrices A and B.

Problem Formulation
Consider a discrete-time system where x(k) ∈ R n , u(k) ∈ R m , w(k) ∈ R p and y(k) ∈ R q are the state, the control input, the disturbance and the measured output, respectively. A ∈ R n×n , B ∈ R n×m , C ∈ R q×n and D ∈ R n×p are constant matrices. Additionally, assume that r(k) ∈ R q is a reference signal, which is generated by the following system: Moreover, suppose that the disturbance signal w(k) is produced by where v(k) ∈ R d denotes the state of system (2), S 1 ∈ R d×d , S 2 ∈ R p×p and F 1 ∈ R q×d are constant matrices.
To facilitate the description of the problem, the following assumptions associated with the Equations (1)-(3) are required.

Hypothesis 2 (H2). (A, B) is stabilizable and (C, A) is detectable.
Hypothesis 3 (H3). S 1 and S 2 have no eigenvalue inside the unit circle in the complex plane.
Hypothesis 4 (H4). both of r(k) and w(k) are previewable, and the preview lengths are h r = s r N and h d = s d N, respectively, which means that values of r(k + 1), r(k + 2), · · · , r(k + h r ) and w(k + 1), w(k + 2),.., w(k + h d ) are known beforehand. Furthermore, s r and s d are nonnegative integers.

Remark 1.
For assumption H4, there are many studies in the literature talking about the application of previewable reference signal, such as [14,15]. In what follows, we take the active suspension of the vehicle as an example to explain the application and significance of disturbance previewable information in real systems. As shown in [35], moving objects such as vehicles and mobile robots will vibrate when driving on uneven roads. For the vibration damping system of the mobile objects, the unevenness of the road is actually a kind of external disturbance and a natural idea is to use the active suspension device to control vibration. That is to say, the road conditions, such as concavities and convexities, in front of the mobile objects are detected by some sensors, and the collected data can be viewed as preview information and applied in the control of the suspension device. Since the front or future road information are known in advance, there is enough time for the suspension device to choose suitable control policy to against vibration.
Define the following error signal: This study aims to design a controller with preview feedforward compensation based on internal model principle, which achieves lim t→∞ e(k) = 0 when t → +∞ .

Discrete Lifting and Problem Translation
In general, state feedback takes the form of u(k) = Kx(k) + v(k), where v(k) is a new assistant input. According to H1, state vector x(k) and output vector y(k) only obtain sampling value at k = iN (i = 0, 1, · · · and N ∈ Z + ). However, the input vector u(k) has sampled data N times rapidly. Therefore, x(k) cannot be used in state feedback when k iN, which brings difficulties to design a controller for multirate systems (1).
Discrete lifting technique can usually be applied to transform the multirate system into the single-rate system. The main method is to obtain the sampling values of each variable at the same time by lifting the dimension of the signal. For the convenience of subsequent derivation, we give the following lifting process according to (1): Using (1), Formula (5) can be further written as The rest can be done in the same procedure, that is where j = 1, 2, · · · , N. Assume that . . .
then according to Formulas (5)-(7), the single-rate system is obtained: In addition, for the error signal e(iN) = y(iN) − r(iN), set . . .
then, it follows from the system (7) that For system (8), (A 1 , B 1 ) was verified to be stabilizable in [28]. Using the duality principle, it can be easily verified that: if (C, A) is detectable, then (C 1 , A 1 ) is detectable. Moreover, to address the original problem, the following assumption is necessary for system (8).

Hypothesis 5 (H5). For any
is assumed to be of full row rank.
Next, we employ the state augmentation technique to establish an augmented system with respect to the tracking error, state variable difference and preview information. A basic consideration is that if a suitable controller is designed to make the closed-loop augmented system asymptotically stable, then as the state component, the tracking error can be asymptotically stabilized to zero. Simultaneously, the preview information will be involved in the controller by means of the state feedback.
Note that A 1 X(i), B 1 U(i) and D 1 W(i) are essentially vectors in (8) and so is A 1 X(i) + B 1 U(i) + D 1 W(i). Then, utilizing the first-order forward difference operator to both sides of the first equation of system (8), we have Moreover, taking the same operation for (9) yields Combining the above two formulas gives , then the above equation can be written as Note that the reference signal and disturbance signal are previewable and the preview lengths are h r = s r N and h d = s d N, respectively. Based on the construction of ∆R(i) and ∆W(i), we know that the future values ∆R(i), ∆R(i + 1), · · · , ∆R(i + s r − 1) and ∆W(i), ∆W(i + 1), · · · , ∆W(i + s d − 1) are available for system (12). In order to design a controller with preview information, the following auxiliary systems are constructed, where . . .

Remark 2.
The benefits of constructing the augmented system (15) are twofold: (1) Since the future values of the reference and disturbance signals are incorporated in the state vector of the augmented system (15), the preview feedforward compensation can be realized through state feedback; (2) The augmented system can deal with the circumstance of existing preview information or not simultaneously, which makes the design method proposed in this study more applicable. In what follows, by exploiting systems (2) and (3), ∆R(i + s r ) and ∆W(i + s d ) can be rewritten as and W(i) was defined in Formula (8). Furthermore, by virtue of the relationships among ∆R(i), ∆W(i), ∆X r (i) and ∆X w (i), as well as the Formula (16), the augmented system (15) can be described as where Consider the non-decreased amplitude of r(k) and w(k), thus ∆ w(i) is an unstable signal. To achieve the control goal of zero-error tracking, the internal model principle will be applied to solve the problem in the sequence. Denote S = diag(S 1 , S 2 ) p+q , the definition of the least p + q-fold internal model of the matrix S is introduced first.
Definition 1 [36,37]. If (G 1 , G 2 ) admits the following form where β is a l × l constant matrix and σ is a constant column vector with the dimension of l and (1) the pair (β, σ) is completely controllable; (2) the characteristic polynomial and the minimal polynomial of β, as well as the minimal polynomial of S are equivalent, then (G 1 , G 2 ) is said to incorporate a least p + q-fold internal model of the matrix S.
Considering the characteristic of the multirate system, a dynamic compensator is proposed based on the least p + q-fold internal model and tracking error e(k), k = iN.
By employing the discrete lifting technique to the compensator mentioned above, the following dynamic state feedback for system (17) is constructed, where , moreover, X c (i), G 1 and G 2 , respectively, have the form of . . . Substituting the dynamic feedback (18) into system (17) and setting X(i) = X T (i) X T c (i) T yield the following closed-loop system where Through the above procedure, the preview tracking problem of the multirate system (1) is eventually formulated as the output regulation problem with respect to system (19). According to [36], the latter problem is solvable provided that the dynamic feedback (18) makes the closed loop system (19) satisfy the following two properties: (1) homogeneous system is exponentially stable; (2) for any initial condition X c (0) and w(i), lim i→∞ E(i) = 0 holds.
Note that the form of E(i) and the sampling time k = iN, if the controller ∆U(i) satisfies the property 2), then it follows that lim k→∞ e(k) = 0. That is to say, ∆U(i) solves the preview tracking problem of the multirate system (1). Correspondingly, by the definition of the difference operator ∆, the control law U(i) can be derived through iterative calculation, which also contains the specific expression of u(k + j) j = 0, 1, · · · , N − 1.

Main Results
In this section, sufficient conditions are given to ensure the solvability of the output regulation problem. Meanwhile, the controller needed for the output regulation problem and the original problem will be designed through the optimal control method. Before giving the main results, we first introduce the following lemma, Lemma 1 [37]. Under H3, assume that (G 1 , G 2 ) incorporates a minimal p + q-fold internal model of the matrix S. If ÂB G 2Ĉ G 1 is exponentially stable for the matricesÂ,B andĈ with appropriate dimension. Then, for any matricesÊ andF, the following regulation equation exists a unique solution pair (X, Z). Moreover, X also satisfieŝ The following theorem gives sufficient conditions for implementing the output regulation of the system (19).

Theorem 1.
Suppose that (G 1 , G 2 ) incorporates a minimal (p + q)-fold internal model of the matrix S. If A c is stable, then the dynamic feedback controller (18) makes the system (19) achieve the output regulation.
Proof. According to [36], the asymptotical stability of the system (20) implies its exponential stability. Therefore, if the matrix A c is stable, then the property 1) holds. Next, we prove lim t→∞ E(i) = 0.
In fact, it follows from the basic knowledge of linear algebra that the least polynomial of β is equal to the least common multiple with respect to the least polynomials of S 1 and S 2 . Hence the least polynomial of G N 1 is the least common multiple corresponding to the least polynomial of S N 1 and S N 2 , which implies that the least polynomial of G N 1 is equal to that of S. In addition, it is apparent that (G 1 , G 2 ) is completely controllable. Recalling the Definition 1, it follows that (G 1 , G 2 ) incorporates a minimal (p + q)N-fold internal model of the matrix S. Under the condition that A c is stable, letÊ = B,F = 0, according to Lemma 1, there exists a unique solution pair (X, Z) for the following regulation Equation with X satisfyingĈX = 0. Let Γ = [X T , Z T ] T and take the coordinate transformation X c (i) = X(i) − Γ∆ w(i), using system (19) and regulation Equation (23) gives Due to the fact that A c is stable, then X c (i) asymptotically stabilizes to the origin. Performing the coordinate transformation to E(i) = C c X(i), we have SinceĈX = 0, the result of the property 2) is obvious. Q.E.D. From Theorem 1, it can be observed that the key to achieve the output regulation is the stability of A c . Note that (20) can be redescribed as follows, Corresponding to the tracking problem, we introduce the following quadratic performance index function for system (26). where with Q = diag(Q E , 0, Q C ), then utilizing the conclusions of optimal control [38] gives the following result.

Lemma 2.
If (A 0 , B 0 ) is stabilizable and ( Q 1/2 , A 0 ) is detectable, the optimal control input of the system (27) with the minimal performance index function (28) is where K = −(H + B T 0 PB 0 ) −1 B T 0 PA 0 . P is a symmetric positive semidefinite matrix satisfying the following discrete algebraic Riccati equation and A 0 + B 0 K is stable.
Proof. From the PBH rank criterion [39], the stabilizability of (A 0 , B 0 ) will be proved if the matrix sI − A 0 B 0 is shown to be of full row rank for any |s|≥ 1 . In accordance with the structures of A 0 and B 0 , sI − A 0 B 0 can be expanded as  s  I  C  D  I  B  sI A  D  B  sI  I  sI  I   sI  sI  I  sI  I  sI  G s Therefore, it suffices to prove that ( ) M s has full row rank for any | | 1 s ≥ under the given conditions.
For the case of It can be easily seen that In view of the controllability of Because sI is nonsingular (|s|≥ 1 ), Using sI to perform the row elementary transformation on the above formula yields rank sI − A 0 B 0 = qs r N + ps d N + rankM(s) where Therefore, it suffices to prove that M(s) has full row rank for any |s|≥ 1 under the given conditions.
For the case of s = 1, It can be easily seen that In view of the controllability of (G 1 , G 2 ), the matrix G 2 I − G 1 has full row rank. In addition, by H5, it follows that has full row rank. Therefore, M(s) is of full row rank for s = 1.
For the case of s 1 and s Λ(S 1 ) ∪ Λ(S 2 ), it can be derived from H3 and the Definition 1 that the matrices (s − 1)I and sI − G 1 are nonsingular. performing column elementary transformation on M(s) leads to In light of the property that the elementary transformation does not change the rank of a matrix, which gives rankM(s) = qN + n + rank sI − A 1 B 1 In fact, we have known that (A 1 , B 1 ) is stabilizable, therefore M(s) has full row rank for s 1 and s Λ(S 1 ) ∪ Λ(S 2 ).
Finally, when s 1 and s ∈ Λ(S 1 ) ∪ Λ(S 2 ), M(s) is changed into Therefore, matrices M(s) and have the same full row rank conditions. Similar to the proof of Theorem 2 in [34], this matrix can be expressed as the product form of two matrices, namely According to the basic knowledge in linear algebraic, for any matrices A ∈ R m×n and B ∈ R n×p , the following conclusion holds: It is noted that s 1 and (G 1 , G 2 ) is controllable, hence has full row rank for any s 1. Furthermore, under assumption H5, it follows that the matrix also has full row rank for s ∈ Λ(S 1 ) ∪ Λ(S 2 ). It follows that which shows that M(s) has full row rank for s 1 and s ∈ Λ(S 1 ) ∪ Λ(S 2 ). Summarizing what we have discussed gives the conclusion of Theorem 2. Q.E.D.
The following result is the detectability of ( Q 1/2 , A 0 ). Proof. Based on the PBH rank criterion [39], Note that the detectability of ( , ) C A implies the detectability of 1 1 ( , ) C A , namely, Obviously, it follows from Theorems 2 and 3 that Lemma 2 is valid under H1-H5. For the dynamic controller (18), let [ ] , then the following conclusion associated with Theorem 1 is given without proof. , where P satisfies the algebraic Riccati equation (30).
Note that under H1-H5, the controller (18) has the ability to make the output regulation problem solvable, as a result, we have preview tracking problem of the multirate system (1) can also be addressed with the controller (18).
Due to Q e > 0, Q c > 0 and sI is an invertible matrix (|s|≥ 1 ), a routine row elementary transformation gives rise to Note that the detectability of (C, A) implies the detectability of (C 1 , A 1 ), namely, is of full column rank for |s|≥ 1 . Therefore, if (C, A) is detectable, then ( Q 1/2 , A 0 ) is detectable and vice versa. Q.E.D.
Obviously, it follows from Theorems 2 and 3 that Lemma 2 is valid under H1-H5. For the dynamic controller (18), let K 1 K 2 = K, then the following conclusion associated with Theorem 1 is given without proof.

Theorem 4.
Under H1-H5, suppose that (G 1 , G 2 ) incorporates a minimal p + q-fold internal model of the matrix S = diag(S 1 , S 2 ), then the dynamic feedback controller (18) can realize the output regulation of the system (19) and the gain matrices K 1 and K 2 are determined by K 1 K 2 = K = −(H + B T 0 PB 0 ) −1 B T 0 PA 0 , where P satisfies the algebraic Riccati equation (30).
Note that under H1-H5, the controller (18) has the ability to make the output regulation problem solvable, as a result, we have lim l→∞ E(i) = 0 and therefore lim k→∞ e(k) = 0 (k = iN). That is, preview tracking problem of the multirate system (1) can also be addressed with the controller (18). In the following, we will derive the specific form of the preview controller required for the original problem.
According to the components of X(i), we now split K 1 into the following form: Correspondingly, ∆U(i) in (18) can be expressed as Using the definition of difference operator and replacing i with i − 1, the above equation can be written as It is noted that U(i) = [u(iN) T u(iN + 1) T · · · u(iN + N − 1) T ] T , to obtain the control input for solving the original problem, we further split K E , K x , K r (l), K d (m) and K 2 into Accordingly, it can be obtained from the Formula (33) that where i = 1, 2, · · · ; j = 0, 1, · · · , N − 1. Furthermore, Denote Mathematics 2020, 8,1443 15 of 20 then it follows from the expression of u(iN + j) that Based on the Equation (34) and Theorem 4, we immediately get the conclusion about the preview tracking problem.

Remark 4.
This paper considers the situation of no delay in control input. In fact, if the disturbance can be measured and a control signal can be applied before it reaches the system, then there will exist time-delay in control input, which, according to [40], is termed as lead-lag. However, as pointed out in [40], a delay is usually needed to guarantee that the feedforward compensation control from measurable disturbances is not made too early. Therefore, it will be very important if the internal model-based preview tracking control for multirate systems with input delay and disturbances can be done.

Remark 5.
In [41], the preview tracking problem was considered for discrete-time output multirate sampling systems, where the reference signal was assumed to converge to a constant vector. This research further considers the same issue with the scenarios that the reference signal being tracked, and disturbance signal being rejected are unattenuated and previewable simultaneously, which has benefits to expand the application of preview control theory. Additionally, when S 1 = I and disturbance signal w(k) is just assumed to tend to zero as time approaches infinity, assumption H5 will degenerate to Assumption 8 in [41] and the internal model compensator (the first equation of (18)) will be reduced to the discrete integrator (18) in [41], then the method and results presented in this research can be applicable to the problem in [41].

Numeric Simulation
In this section, a numeric example is presented to demonstrate the validity of the results. Consider the following discrete time system Compared with system (1), the coefficient matrices are as follows By utilizing the PBH rank criterion, we can prove that (A, B) is stabilizable and (C, A) is detectable. Here, the state and output of the above system are assumed to be measured at k = 3i (i = 1, 2, · · · ).
The reference signal and the disturbance signal are generated by the following two external systems, respectively and A simple calculation gives that the eigenvalues corresponding to matrices (1) cos (1) are, respectively λ 1,2 = √ 2 2 ± √ 2 2 i and λ 3,4 = cos(1) ± i · sin(1). It is obvious that these eigenvalues are located on the unit circle. According to Remark 3, (1)), (1)), α 0 = 1. Since the dimension of the reference signal is 1, it follows from Definition 1 and Remark 3 that the coefficient matrices of the internal model compensator are given by A simple calculation gives that the eigenvalues corresponding to matrices . It is obvious that these eigenvalues are located on the unit circle. According to Remark 3, let In addition, through lifting the multirate system, it can be verified that the matrix Assume that the both the reference signal and the disturbance signal are previewable. In order to reflect the influence of the preview information on tracking effectiveness, the numeric simulation will be conducted under the following four scenarios, i.e.,  (27). Then, according to the above discussions, the conditions of the Theorem 5 are all satisfied. As a result, the controller designed in this paper can drive the multirate system to complete the preview tracking. Using MATLAB to perform the simulation.
As what can be seen from Figures 1-6, the larger the preview length, the better the tracking performance. That is, the controller (34) has advantage of overcoming the external disturbance in a shorter time and track the preview signal more accurately. Therefore, the controller (34) designed for the preview tracking problem is effective.
In addition, through lifting the multirate system, it can be verified that the matrix as well as I − A 1 B 1 C 1 B 2 , has full row rank, where i = 1, 2, 3, 4.
Assume that the both the reference signal and the disturbance signal are previewable. In order to reflect the influence of the preview information on tracking effectiveness, the numeric simulation will be conducted under the following four scenarios, i.e., h r = 3, h d = 6; h r = 30, h d = 6; h r = 6, h d = 3; h r = 6, h d = 15. Specifically, select the initial values x 0 = 0.01 0.01 T and x c = 0 0 0 0 T for the system and the internal model compensator, respectively. Moreover, choose the weight matrices Q E = 100 · I 3 , Q c = 10 · I 4 and H = 5000 · I 3 for the performance index function (27). Then, according to the above discussions, the conditions of the Theorem 5 are all satisfied. As a result, the controller designed in this paper can drive the multirate system to complete the preview tracking. Using MATLAB to perform the simulation.
As what can be seen from Figures 1-6, the larger the preview length, the better the tracking performance. That is, the controller (34) has advantage of overcoming the external disturbance in a shorter time and track the preview signal more accurately. Therefore, the controller (34) designed for the preview tracking problem is effective.

Conclusions
In this research, the preview tracking problem for discrete-time multirate systems was studied with external disturbance. Discrete lifting method and state augmentation technique, which together with internal model principle, were adopted to transform the original problem into the output regulation problem for a single-rate augmented system. This study is the development of preview control for discrete-time linear multirate systems [34] since the scenario of the unattenuated reference signal and disturbance signal was considered. Note that the state and output signals in this research were obtained by equal period sampling, in practical, since the bandwidth of network transmission is limited, quantized sampling is unavoidable. Future work will focus on considering the preview tracking control for discrete-time linear systems with input and output quantization.