Next Article in Journal
Analytical Study of a Circular Thin Plate Contacting with an Elastic Sphere
Next Article in Special Issue
Dynamic Analysis of Fuzzy Systems
Previous Article in Journal
Composition Component Influence on Concrete Properties with the Additive of Rubber Tree Seed Shells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discrete-Time Design of Dual Internal Model-Based Repetitive Control Systems

1
Research Center for Photonics, National Research and Innovation Agency, Tangerang Selatan 15314, Indonesia
2
Directorate of Agricultural Irrigation, Ministry of Agriculture, Jakarta Selatan 12550, Indonesia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(22), 11746; https://doi.org/10.3390/app122211746
Submission received: 13 October 2022 / Revised: 10 November 2022 / Accepted: 16 November 2022 / Published: 18 November 2022
(This article belongs to the Special Issue Advances in Intelligent Control and Engineering Applications)

Abstract

:
This paper presents a novel design of discrete-time dual internal model-based repetitive control systems. The design strategy is accomplished by combining general and high-order modified repetitive control schemes for simultaneous tracking repetitive tasks and rejection of uncertain periodic disturbances. The proposed controller is constructed from two different discrete-time internal models, rendering a dual internal model-based repetitive controller (DIMRC). The first internal model is intended to track repetitive commands with a fixed fundamental frequency. Meanwhile, the second internal model is coupled to compensate for an exogenous periodic disturbance with an uncertain frequency. The controller structure, stability conditions, and convergence analysis are discussed in this paper. The performance of the proposed controller is validated through simulation studies showing accurate tracking and excellent disturbance rejection simultaneously.

1. Introduction

Repetitive control (RC) is a learning control scheme utilizing an internal model principle by Francis and Wohnam [1] for accurate reference tracking or good disturbance rejection of the periodic signals. The inclusion of internal model inside the feedback loop establishes a generator of periodic control signal, enabling a null steady-state tracking error. The early applications of RC for reference tracking and disturbance rejection problems were listed by Hillerstrom and Walgama [2], followed by Kurniawan et al. [3]. Recently, RC was developed and applied in many different applications such as wind-turbines [4], inverter compressor refrigeration fields [5], magnetically suspended rotor systems [6], functional electrical stimulation [7], nano-positioning systems [8], centrifugal compressors [9], and many others.
In principle, RC design comprises two parts, i.e., a stabilizing controller and an internal model. In the RC-controlled system, the stabilizing controller, sometimes called the compensator/learning function, is required to stabilize the closed-loop model. In addition, the stabilizing controller also determines the convergence speed of the system error. On the other hand, the internal model represents the reference/disturbance model, which later behaves as a periodic signal generator. This internal model makes accurate reference tracking/good disturbance rejection possible.
Some important works related to the internal model design can be found in [10,11,12,13,14,15]. The first digital internal model obtained by discretizing the continuous-time delay was given in [10]. The internal model [10] has a finite-dimensional structure and provides a null steady-state tracking error for any periodic signal at a fundamental frequency and its harmonics up to Nyquist frequency components. The internal model [10] was successfully applied to the control of a single-axis electric servomechanism using a 16-bit microprocessor. In [11], a low-pass filter q ( z ) is added to discrete-time internal model [10] intended to improve the system robustness against high-frequency models. Hillerstrom and Sternby [12] presented a low-order internal model targeting specific band-limited periodic references/disturbances. The discrete-time internal model specifically used to track/compensate odd-harmonic periodic references/disturbances was proposed in [13]. In [14], the discrete-time design of the internal model was developed to track or suppress periodic signals consisting of two or more fundamental frequencies. A discrete-time higher-order internal model for tracking/rejecting periodic signals with uncertain periods was introduced in [15]. In addition to some works on the internal models, several recent works on the development and application of RCs can be found in [16,17,18,19]. A robust RC with an optimized-band-stop filter for achieving high accuracy in nanopositioning stages was proposed in [16]. In [17], the discrete-time RC with an active disturbance rejection control was implemented to realize high-precision permanent magnet synchronous motor drives. A design of higher-order RC with a phase lead stabilizing controller applied for a two-level grid-connected inverter was presented in [18]. Lu et al. [19] developed the high-order selective harmonic RC scheme to improve the steady-state tracking accuracy of the pulse-width modulation converters. The works in [10,11,12,13,14,15,16,17,18,19] aimed to develop and utilize the internal model-based RC for either periodic reference tracking or disturbance rejection. This also implies that the proposed RC designs are not intended to simultaneously track and reject periodic signals.
In this work, we focus on the internal model design used for simultaneous reference tracking and disturbance rejection. More specifically, we develop a control strategy for tracking repetitive reference and rejection of uncertain periodic disturbance of a discrete-time linear system. This is motivated by the fact that tracking and rejection are common problems in control systems. Tracking repetitive trajectories can be found in many applications, such as robotics, power sources, precision gantry systems, engine valve systems, etc. [3]. Here, we also consider a disturbance rejection problem, where the disturbance period is uncertain or varying. Such disturbance appeared in many practical situations, such as rotating machines [20,21], active suspension systems [22,23], steel casting processes [24,25], and many others. When the actual disturbance period is subject to variation, the general internal model will operate with a period mismatch. This condition makes the RC gains at the disturbance’s frequencies drastically drop to low-level magnitudes. Consequently, the perfect disturbance attenuation is no longer applicable and the tracking accuracy is degraded. Therefore, the higher-order internal model is adopted to improve the attenuation performance when the disturbance period is varied. In this paper, a systematic approach to integrating two different internal models, i.e., general and higher-order internal models, is presented. Note that the basis periods of reference and disturbance signals are different and uncorrelated. The stability analysis, the controller structure, and controller realization are also discussed. The proposed method renders a controller called a dual internal model-based RC (DIMRC). Simulation results and comparison studies validate the efficiency of the proposed DIMRC. The main contributions of this work are summarized as follows:
  • A novel dual internal model-based RC is constructed by using the denominator parts of the general and higher-order internal models.
  • Stability conditions of the plug-in DIMRC system are presented. The stability conditions are then used to determine the stabilizing controller.
  • The structure and realization of DIMRC for the non-causal stabilizing controller are developed.
The remainder of this paper is structured as follows. Section 2 describes the problem statement and several assumptions used in the design. In Section 3, the fundamental notions of the discrete-time general RC and high-order RC are presented. Section 4 discusses the proposed method covering the controller structure, the closed-loop system stability, and the controller realization. Simulation results covering minimum-phase and non-minimum phase cases, followed by discussion, are provided in Section 5. Lastly, concluding remarks are drawn in Section 6.

2. Problem Statement

In this work, we consider the following discrete-time LTI system:
y ( k ) = P ( z ) u ( k ) + v ( k ) ,
where y ( k ) , u ( k ) , v ( k ) R are, respectively, the discrete-time plant output, control input, and external disturbance, and P ( z ) is the plant model. Throughout this paper, we use notation z denoting the Z -transform variable, and also the forward shifting operator, e.g., z u ( k ) = u ( k + 1 ) . Hence, the notation z 1 denotes a backward shifting operator, e.g., z 1 u ( k ) = u ( k 1 ) . To ease the writing, the operation P ( z ) u ( k ) is used to equally represent Z 1 { P ( z ) u ( z ) } . Here, P ( z ) can be either a polynomial or transfer function in the Z domain, and u ( z ) is Z -transform of the discrete-time control signal u ( k ) .
The plant model of the LTI system (1) is expressed as
P ( z ) = N ( z ) D ( z ) ,
where N ( z ) , D ( z ) C are numerator and denominator polynomials of the plant. The desired trajectory to be tracked is referred to as a reference signal r ( k ) , and the tracking error e ( k ) is defined by
e ( k ) = r ( k ) y ( k ) .
In the design of the proposed controller, the following assumptions are used:
Assumption 1. 
The polynomials N ( z ) and D ( z ) in (2) are known. In addition, the plant model P ( z ) (2) is stabilized with a conventional controller C ( z ) , resulting in a stable transfer function. The stabilized plant model can be either a minimum or non-minimum phase system.
Assumption 2. 
The reference r ( k ) is a repetitive signal with a fixed and known fundamental frequency f r . The disturbance v ( k ) is periodic with an uncertain frequency. However, its nominal fundamental frequency f ¯ v is known. The actual frequency of disturbance v ( k ) may slightly vary from its nominal value.
Assumption 3. 
The reference frequency f r and disturbance frequency f ¯ v are uncorrelated. In other words, they are not harmonics, i.e., f r n f ¯ v , where n is an arbitrary positive integer number.
The research objective is to synthesize control law u ( k ) to simultaneously track the repetitive reference r ( k ) and suppress the effect of the uncertain periodic disturbance v ( k ) of the LTI system (1), such that the tracking error e ( k ) (3) converges to a small steady-state value, and the resulting closed-loop system is stable. In addition, the proposed controller is realizable and applicable for both minimum and non-minimum phase systems.

3. Discrete-Time RC

3.1. A General Modified RC

Suppose that the LTI system (1) is subject to the repetitive reference r ( k ) with no presence of disturbance v ( k ) . Let the reference r ( k ) have a fundamental frequency of f r . This gives a basis reference period as T r = 1 / f r . Based on the information of plant model P ( z ) and reference period T r , a general discrete-time RC can be constructed to accurately track the reference signal r ( k ) . The general modified RC (GMRC) has a transfer function as follows:
U g ( z ) E ( z ) = q r ( z ) z N r 1 q r ( z ) z N r F ( z ) ,
where U g ( z ) and E ( z ) are the control input and the tracking error in Z -domain, respectively, q r ( z ) is a zero-phase low pass filter used to cut off higher frequencies of reference r ( k ) , F ( z ) is a stabilized controller used for stabilizing the RC-closed loop system, and N r is an integer number of samples per reference period.
An integer N r in (4) is obtained from
N r = T r T s = 1 f r T s .
Here, T r is a reference period, and T s is a sampling time. Equation (5) emphasizes that Assumption 2 is compulsory in the design of GMRC. The term modified in GMRC refers to the inclusion of q-filter q r ( z ) to the time-delay z N r . The use of q r ( z ) improves the robustness of RC against high-frequency components. In addition, the q r ( z ) smooths the control signal generated by the traditional RC. The q r ( z ) is generally designed as a zero-phase low-pass filter given by
q r ( z ) = α 0 + i = 1 h r α i z i + z i .
Note that h r is the order of filter q r ( z ) , and the filter coefficients are chosen such that
α 0 + 2 i = 1 h r α i = 1 .
The condition (7) ensures the unity gain of the filter at the frequencies below the filter bandwidth ( ω q r ); that is, q r ( ω ) = 1 for ω < ω q r ( ω q r = 2 π f q r ). In addition, the q-filter (6) contributes zero-phase for all frequencies, i.e, θ q r ( ω ) = 0 o 0 < ω < π / T s . Here, q r ( ω ) and θ q r ( ω ) correspond to the magnitude and phase responses of q r ( z ) , respectively.
The F ( z ) in (4) is also a crucial part of RC. In addition to stabilizing the RC closed-loop system, the F ( z ) also determines the convergence rate of the tracking error e ( k ) . The F ( z ) is often designed as an exact inverse of the stabilized plant model. This design requires that the stabilized plant is in a class of minimum-phase systems. To design F ( z ) for non-minimum phase systems, one can refer to a zero-phase error tracking controller (ZPETC) design proposed in [26].
Next, we define an internal model of the controller (4) as follows:
I g ( z ) = N g ( z ) D g ( z ) = q r ( z ) z N r 1 q r ( z ) z N r ,
where N g ( z ) and D g ( z ) are the numerator and denominator parts of the general internal model. The denominator D g ( z ) in (8) is used later in the synthesis of the proposed controller.

3.2. A Higher-Order Modified RC

Here, the control problem for the LTI system (1) is disturbance rejection only, implying that the reference r ( k ) is set to zero. The following input–output relation expresses a discrete-time higher-order modified repetitive controller (HOMRC) [15]:
U h ( z ) E ( z ) = q v ( z ) j = 1 m w j z j N v 1 q v ( z ) j = 1 m w j z j N v F ( z ) ,
where U h ( z ) is a control input generated by HOMRC, q v ( z ) is the q-filter used for cutting off higher frequencies of the disturbance v ( k ) , F ( z ) is the stabilizing controller, and N v is the number of samples per disturbance period. In contrast with (5), N v is an integer number given by
N r = T ¯ v T s = 1 f ¯ v T s ,
where T ¯ v is a nominal disturbance period which is also known based on Assumption 2.
Similar to (4), we also determine the internal model of the HOMRC given by
I h ( z ) = N h ( z ) D h ( z ) = q v ( z ) j = 1 m w j z j N v 1 q v ( z ) j = 1 m w j z j N v ,
It can be seen from (11) that HOMRC employs a sum of the weighted delays, i.e, j = 1 m w j z j N v . The weights w j for j = 1 , 2 , , m are chosen to satisfy conditions as follows [15]:
j = 1 m w j = 1 ,
j = 1 m j m 1 w j = 0 .
The conditions (12) and (13) are added to ensure that the high gains of HOMRC at the fundamental frequency and its harmonics remain infinite. In addition, the high-gain peaks of HOMRC at the harmonics are extended to wider regions. This behavior can be seen from the magnitude responses of GMRC and HOMRC at the targeted fundamental frequency 1 Hz shown in Figure 1. Hence, the disturbance-rejection performances at the neighboring frequencies are improved by using the HOMRC.

4. Proposed Dual Internal Model-Based RC (DIMRC)

4.1. DIMRC Structure

Observing the internal models (8) and (11), we notice that both models show a similar form as follows:
I g , h ( z ) = S g , h ( z ) 1 S g , h ( z ) ,
where S g , h ( z ) corresponds to the delay term of GMRC/HOMRC, and formulated by
S g ( z ) = q r ( z ) z N r ,
S h ( z ) = q v ( z ) j = 1 m w j z j N v .
Hence, the proposed DIMRC is designed to have the following transfer function:
C d ( z ) = U d ( z ) E ( z ) = S d ( z ) 1 S d ( z ) F ( z ) ,
which is equivalent to the DIMRC law as follows:
U d ( z ) = S d ( z ) U d ( z ) + S d ( z ) F ( z ) E ( z ) .
Note that S d ( z ) in (17) is the delay term constructed based on the denominators D g ( z ) and D h ( z ) given in (8) and (11), respectively. The S d ( z ) is also in the form of polynomial and is proposed as follows:
S d ( z ) = 1 D g ( z ) D h ( z ) .
Substituting D g ( z ) in (8) and D g ( z ) in (11) to (19), we obtain
S d ( z ) = 1 1 q r ( z ) z N r 1 q v ( z ) j = 1 m w j z j N v ,
which can be further expressed as
S d ( z ) = q r ( z ) z N r + q v ( z ) j = 1 m w j z j N v q r ( z ) q v ( z ) j = 1 m w j z j N v N r .
Finally, the block diagram of the proposed DIMRC closed-loop system is illustrated in Figure 2. In addition, the controller realization is described in Figure 3. The weighted delays shown in Figure 3, is then detailed in Figure 4. Note that the controller realization shown in Figure 3 is implementable for DIMRC with causal or proper stabilizing controller F ( z ) . Some F ( z ) designs resulting in causal form can be seen in [27,28]. Here, the causal stabilizing controller F ( z ) can be implemented separately without being merged with the internal model. Other designs resulting in non-causal/improper form can be found in the references [26,29,30]. Thus, it is obvious that the stabilizing controller F ( z ) can be designed in either causal or non-causal form. Suppose that the stabilizing controller F ( z ) is non-causal by a factor of d F . This means that deg { N f ( z ) } > deg { D f ( z ) } and deg { N f ( z ) } deg { D f ( z ) } = d F , where deg { N f ( z ) } and deg { D f ( z ) } represent the degrees of F ( z ) ’s numerator N f ( z ) and denominator D f ( z ) , respectively. To reduce the controller’s complexity, the realization of DIMRC with non-causal F ( z ) can be achieved by modifying Figure 3 into Figure 5.

4.2. Stability of the Plug-In RC System

In this subsection, we develop a sufficient stability condition for the closed-loop system regulated by plug-in DIMRC. The sensitivity function of the plug-in DIMRC system shown in Figure 2 is given by
E ( z ) R ( z ) = 1 1 + { 1 + C d ( z ) } C ( z ) P ( z )
= 1 S d ( z ) 1 + P ( z ) C ( z ) 1 + P ( z ) C ( z ) S d ( z ) + P ( z ) C ( z ) F ( z ) S d ( z ) .
Equation (23) can be factorized to
E ( z ) R ( z ) = 1 S d ( z ) 1 + P ( z ) C ( z ) 1 S d ( z ) + { P ( z ) C ( z ) 1 + P ( z ) C ( z ) } F ( z ) S d ( z )
Let define the closed-loop plant model P c ( z ) as follows:
P c ( z ) = P ( z ) C ( z ) 1 + P ( z ) C ( z ) .
The sensitivity function (24) can be further derived to
E ( z ) R ( z ) = 1 S d ( z ) 1 + P ( z ) C ( z ) 1 { 1 P c ( z ) F ( z ) } S d ( z ) .
This leads to the following characteristic equation (CE):
CE : 1 + P ( z ) C ( z ) 1 { 1 P c ( z ) F ( z ) } S d ( z ) = 0 .
Based on (27), the closed-loop DIMRC system in Figure 2 is internally stable if the following two conditions are fulfilled:
  • The first factor 1 + P ( z ) C ( z ) in (27) has stable roots. This also implies that the closed-loop plant model P c ( z ) in (25) is a stable transfer function.
  • The second factor 1 { 1 P c ( z ) F ( z ) } S d ( z ) in (27) has stable roots. Following the work in [13], the stability of the second part can be assured by
    { 1 P c ( z ) F ( z ) } S d ( z ) 1 P c ( z ) F ( z ) S d ( z ) < 1 .
    Here, the notation T f ( z ) corresponds to the H -norm infinity of the transfer function T f ( z ) , which can be defined as the peak of the magnitude response of T f ( z ) for all frequency components.
Based on Assumption 1 and picking the conventional controller C ( z ) , P c ( z ) can be computed. Meanwhile, S d ( z ) can be calculated by using Assumption 2, and by choosing the q-filters q r ( z ) and q v ( z ) . Finally, the stabilizing controller F ( z ) can be designed to satisfy the condition (28). It has been mentioned before that several works addressing stabilizing controller design methods can be found in [26,27,28,29,30,31]. Here, we omit the discussion about the stabilizing controller design as we focus on the internal model structure, the controller structure and realization, and the stability of the plug-in DIMRC system.
To summarize, the proposed design is accomplished according to the following steps:
  • Obtain the open-loop plant model P ( z ) (2).
  • Pick the conventional controller C ( z ) , ensuring a stable P c ( z ) (25).
  • Obtain the reference frequency f r and disturbance frequency f ¯ v to calculate the integer numbers N r (5) and N v (10).
  • Pick the q-filters q r ( z ) and q v ( z ) according to (6) and (7).
  • Choose the order of HOMRC m, and determine the weights based on (12) and (13).
  • Construct the S d ( z ) by using (21).
  • Determine the stabilizing controller F ( z ) satisfying (28).
  • Synthesize the complete controller according to (17).
  • Realize the controller based on Figure 3 for causal F ( z ) , and Figure 5 for non-causal F ( z ) .

5. Simulation Results

5.1. Case 1: Minimum Phase Stabilized System

To validate the effectiveness of the proposed controller, the following discrete-time model of Quanser SRV02 servo plant [32] is used in the simulation:
P ( z ) = θ ( z ) V ( z ) = 10 4 7.6337 z + 7.173 z 2 1.83 z + 0.8289 ,
where P ( z ) is an open-loop plant model, θ ( z ) is an angle position (rad), and V ( z ) is an input voltage (Volt). The open-loop plant (29) is associated with the sampling time T s = 0.005 s. The conventional controller used to stabilize the open-loop plant (29) is a proportional controller with the gain 10, i.e., C ( z ) = 10 . The stabilized plant model P c ( z ) (25) is obtained as follows:
P c ( z ) = 10 3 7.634 z + 7.173 z 2 1.822 z + 0.837 .
The stabilized plant model (30) has a stable zero and two stable poles located at z 1 = 0.94 , p 1 = 0.911 + i 0.083 , and p 2 = 0.911 i 0.083 . Thus, the model (30) is classified as the minimum-phase system. The control objective is to track a periodic triangle reference r ( k ) with the fixed fundamental frequency f r = 0.4 Hz and to suppress the uncertain periodic disturbance with known nominal frequency as f ¯ v = 1 Hz. We can notice that f r and f ¯ v are uncorrelated, which is inline with Assumption 3. The reference r ( k ) is illustrated in Figure 6. The integer numbers N r and N v can be calculated as follows:
N r = 1 0.4 × 0.005 = 500 , N v = 1 1 × 0.005 = 200 .
Suppose that the q-filters q r ( z ) and q v ( z ) are equally chosen as
q r ( z ) = q v ( z ) = 0.25 z 1 + 0.5 + 0.25 z .
The chosen filter above gives the tracking/rejection performance with a bandwidth of 36.35 Hz. Note that the q-filter here determines the bandwidth, which selects the reference frequencies to be passed and the disturbance frequencies to be attenuated. The coefficients of the q-filter are calculated according to (6) and (7) to give a unity gain at the targeted bandwidth and a zero phase to all frequency components. The filter behaviors can be seen in Figure 7. Increasing the filter’s degree does not guarantee that the bandwidth is enlarged. The coefficients must be carefully tuned to give the preferred bandwidth. Here, the chosen q-filter already accommodates the targeted reference and disturbance frequencies, as specified before.
Let us pick the order of HOMRC m as 2. Then, the following weights are chosen to satisfy conditions (12) and (13):
m = 2 w 1 = 2 , w 2 = 1 .
Based on (31) and (33), the polynomial S d ( z ) in (21) can now be synthesized. Next, we need to obtain F ( z ) to stabilize the closed-loop DIMRC system. A straightforward design for F ( z ) is
F ( z ) = k r P c 1 ( z ) = k r z 2 1.822 z + 0.837 10 3 7.634 z + 7.173 ,
where k r is the learning gain and P c 1 ( z ) is the inverse model of (30). It can be seen that F ( z ) is non-causal by a factor of 1 ( d F = 1 ). Hence, the proposed DIMRC can be realized by following the controller structure shown in Figure 5. The learning gain k r satisfying condition (28) is chosen as 0.95 . As we have S d ( z ) and F ( z ) , the proposed DIMRC (17) can now be constructed.
To highlight the tracking performance of the proposed DIMRC, a comparison to multiple-periods RC (MPRC) [14] is conducted for two different disturbance models as follows:
v 1 ( t ) = 0.05 sin ( 2 π t ) + 0.03 sin ( 4 π t ) ,
v 2 ( t ) = 0.05 sin ( 2 π 0.95 t ) + 0.03 sin ( 2 π 1.95 t ) .
The MPRC law [14] is expressed as
U m ( z ) = M ( z ) U m ( z ) + M ( z ) F ( z ) E ( z ) ,
where N r and N v are provided in (31), F ( z ) is given in (34), and M ( z ) is a polynomial given by
M ( z ) = q r ( z ) z N r + z N v z N r N v .
Tracking-error performance of both DIMRC and MPRC for the minimum-phase system (30) with disturbance model v 1 ( t ) is shown in Figure 8. We can notice that the tracking errors of both DIMRC and MPRC converge to zero steady-state value. This is obvious because the disturbance model v 1 ( t ) represents the periodic disturbance signal with nominal fundamental frequency as 1 Hz. Figure 8 also indicates that the MPRC outperforms the DIMRC in terms of convergence rate. This phenomenon can be understood because the DIMRC applies a longer control delay due to the use of the higher-order internal model. The control delay (also referred as delay term) of DIMRC is calculated according to (21), and given by
S d ( z ) = q r ( z ) 2 z 200 z 400 + z 500 2 q r ( z ) z 700 + q r ( z ) z 900
Meanwhile, the control delay of MPRC is obtained from (38), and expressed as
M ( z ) = q r ( z ) z 200 + z 500 z 700
We can observe from (39) and (40) that DIMRC has more series of delays and gives longer delays to the control input and the tracking error. Consequently, the DIMRC takes longer to generate the correct control input, so the tracking error converges to a steady state. However, we mainly focus on improving the tracking accuracy, especially during the steady state of the plant under time-varying periodic disturbances. Therefore, the superior performance of DIMRC can be seen from the next simulation results.
Next, we examine the performance of both controllers when the disturbance frequencies vary from their nominal values. The model (36) represents the uncertain exogenous periodic disturbance v ( k ) where the frequency components are shifted to 0.95 Hz and 1.95 Hz. Figure 9 illustrates the tracking errors of the minimum-phase system (30) with disturbance model v 2 ( t ) . As we observe in Figure 9, the transient errors of both controllers indicate similar patterns to the ones shown in Figure 8. In contrast, the steady-state errors are significantly different, showing that zero-error steady state is no longer achieved. These results illustrate that RC systems’ tracking accuracy is greatly affected when the disturbance’s frequency is slightly changed from its nominal value. However, the steady-state error of DIMRC is significantly smaller than that of MPRC. This implies that DIMRC exploiting the higher-internal model offers better robustness against the uncertain periodic disturbance, especially during the steady-state period. This performance can also be assessed from the magnitude responses of the simulated DIMRC and MPRC shown in Figure 10. We notice that the magnitude responses at the reference frequency ( 0.4 Hz) are almost similar between DIMRC and MPRC. However, the magnitude response of the DIMRC is extended to a larger region around the nominal disturbance frequency (1 Hz). As a result, the DIMRC improves the rejection performance at the neighboring disturbance frequencies. Therefore, when the disturbance frequencies vary to 0.95 Hz and 1.95 Hz, the DIMRC gives a better attenuation than the MPRC, resulting in a smaller steady-state error.

5.2. Case 2: Non-Minimum Phase Stabilized System

In this case, the open-loop plant model (29) is stabilized using the conventional controller as follows:
C ( z ) = z + 1.25 z + 0.75 .
The resulting close-loop plant model is
P c ( z ) = 10 3 0.763 z 2 + 0.145 z + 0.896 z 3 1.079 z 2 0.541 z + 0.623 .
The closed-loop plant model (42) has three poles given by p 1 = 0.988 , p 2 = 0.841 , p 3 = 0.750 , and two zeros located at z 1 = 1.25 , and z 2 = 0.939 . It is obvious that the model (42) has stable poles and one unstable zero ( z 1 = 1.25 ) confirming that (42) is a stable non-minimum phase system. Note that the stabilizing controller F ( z ) designed according to (34) is not applicable to the non-minimum phase system. The design method (34) results in unstable stabilizing controller F ( z ) . In this case, we employ a zero-phase tracking error controller (ZPETC) design technique developed in [26]. The ZPETC-based stabilizing controller has the following transfer function:
F ( z ) = k r D c ( z ) N c u ( z 1 ) k p N c s ( z ) .
Here, k r is the learning gain similar to the ones in (34), k p is the stabilized-plant gain, D c ( z ) denotes the denominator of (42), N c s ( z ) and N c u ( z 1 ) represent stable and unstable factors of N c ( z ) , where N c ( z ) is the numerator of (42). The term N c u ( z 1 ) is obtained from N c u ( z ) with the operator z replaced by the backward shift operator z 1 . Using the Equation (43) and choosing the learning gain as k r = 0.2 , we obtain F ( z ) expressed as
F ( z ) = 0.2 1.252 z 4 0.351 z 3 1.756 z 2 + 0.239 z + 0.623 10 3 0.763 z 2 + 0.716 z
We notice that the F ( z ) above is non-causal by a factor of 2. This gives d F as 2, and the DIMRC structure can be adjusted according to Figure 5.
The tracking-error performance of the non-minimum phase system (42) for two different disturbance models v 1 ( t ) and v 2 ( t ) are, respectively, depicted in Figure 11 and Figure 12. From Figure 11 and Figure 12, we notice that the DIMRC system converges more slowly and exhibits larger peaks of the tracking error during transience compared to the MPRC system. However, when the disturbance model is uncertain, the DIMRC system offers better steady-state error. These results indicate that, in term of the transient responses, the MPRC outperforms the DIMRC for both minimum and non-minimum phase systems when there is no variation in the disturbance model. Nevertheless, when the disturbance model slightly changes from its nominal value, the DIMRC significantly outperforms the MPRC, especially during the steady-state period.
In order to highlight the controllers’ performance during the steady state, the detailed comparison in terms of the steady-state time ( t s s ) and root-mean-square of the steady-state error (rms- e s s ) are summarized in Table 1. The t s s is defined as the time needed to reach 2.5 % error ( ± 0.025 ). Meanwhile, the rms- e s s is calculated by
rms - e s s : = 1 ( n t s n s s ) k = n s s n t s e 2 ( k ) ,
where n t s = t s / T s , n s s = t s s / T s , and t s s is the simulation period. It can be seen from Table 1 that the proposed method has a smaller mean square error for both minimum and non-minimum phase systems under uncertain periodic disturbance v 2 ( t ) . However, the proposed controller offers longer steady-state time compared to the MPRC. All these results indicate that the proposed controller provides better robustness against uncertain periodic disturbances, but at the expense of a slower transient response.

6. Conclusions

In this paper, the dual internal model-based repetitive controller has been developed to simultaneously track repetitive tasks and reject uncertain periodic disturbances. The proposed design combines two internal models originating from the general modified and high-order modified repetitive control schemes. The internal model of general modified repetitive control is adopted to handle periodic reference tracking with a fixed frequency. Then, the internal model of high-order modified repetitive control is added to deal with the uncertain periodic disturbance. The controller structure, stability analysis, and controller realization are discussed in this article. Simulation and comparison studies have been conducted to highlight the tracking performance of the proposed controller. However, the proposed design is still limited for the single-input single-output system, implying that the controller is not applicable for the multivariable system. Moreover, the proposed control scheme is unsuitable for handling non-periodic (aperiodic) disturbances. Extending the proposed work for the multivariable system and for handling non-periodic disturbances are challenging tasks and have become our future research investigation.

Author Contributions

Conceptualization, J.A.P. and P.P.; methodology, E.K. and J.A.P.; software, P.P.; validation, E.B.P. and D.R.; formal analysis, E.K.; investigation, E.B.P. and P.P.; resources, S.W. and T.M.; data curation, J.A.P.; writing—original draft preparation, J.A.P.; writing—review and editing, E.K.; visualization, D.R.; supervision, P.P. and S.W.; project administration, T.M. and F.S.; funding acquisition, F.S. All authors have read and agreed to the published version of the manuscript.

Funding

The authors received financial support for publication from the UPLAND Project, with the support of IsDB and IFAD, that focus to increase smallholders’ agriculture productivity, incomes, livelihoods, and resilience in the targeted area. Meanwhile, the research was supported by Research Center for Photonics, National Research and Innovation Agency (BRIN), Indonesia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to express sincere thanks, particularly to Rahmanto, the Director of Agriculture Irrigation at the Ministry of Agriculture, Republic of Indonesia.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DIMRCDual internal model-based repetitive controller.
GMRCGeneral modified repetitive controller.
HOMRCHigher-order modified repetitive controller.
LTILinear time invariant.
MPRCMulti-Periods Repetitive Control.
RCRepetitive controller.
ZPETCZero-phase tracking error controller.

References

  1. Francis, B.A.; Wonham, W.M. The internal model principle for linear multivariable regulators. Appl. Math. Optim. 1975, 2, 170–194. [Google Scholar] [CrossRef]
  2. Hillerstrom, G.; Walgama, K. Repetitive control theory and applications—A survey. In Proceedings of the 13th IFAC World Congress, San Francisco, CA, USA, 30 June–5 July 1996; pp. 1–6. [Google Scholar]
  3. Kurniawan, E.; Cao, Z.; Mahendra, O.; Wardoyo, R. A survey on robust Repetitive Control and applications. In Proceedings of the 2014 IEEE International Conference on Control System, Computing and Engineering, Penang, Malaysia, 28–30 November 2014; pp. 524–529. [Google Scholar]
  4. Liu, Y.; Wu, P.; Ferrari, R.M.G.; Van Wingerden, J.W. Fast adaptive fault accommodation in floating offshore wind turbines via model-based fault diagnosis and subspace predictive repetitive control. IFAC-PapersOnLine 2020, 53, 12650–12655. [Google Scholar] [CrossRef]
  5. Meng, F.; Zhang, X.; Li, Z.; Wen, X.; You, L. Speed Fluctuation Suppression for the Inverter Compressor Based on the Adaptive Revised Repetitive Controller. Energies 2020, 13, 6342. [Google Scholar] [CrossRef]
  6. Cui, P.; Zhang, G. Modified Repetitive Control for Odd-Harmonic Current Suppression in Magnetically Suspended Rotor Systems. IEEE Trans. Ind. Electron. 2019, 66, 8008–8018. [Google Scholar] [CrossRef]
  7. Page, A.P.; Freeman, C.T. Point-to-point repetitive control of functional electrical stimulation for drop-foot. Control Eng. Pract. 2020, 96, 1–10. [Google Scholar] [CrossRef]
  8. Li, L.; Aphale, S.S.; Zhu, L. Enhanced odd-harmonic repetitive control of nanopositioning stages using spectrum-selection filtering scheme for high-speed raster scanning. IEEE Trans. Automat. Sci. Eng. 2021, 18, 1087–1096. [Google Scholar] [CrossRef]
  9. Cai, K.; Deng, Z.; Peng, C.; Li, K. Suppression of Harmonic Vibration in Magnetically Suspended Centrifugal Compressor Using Zero-Phase Odd-Harmonic Repetitive Controller. IEEE Trans. Ind. Electron. 2020, 67, 7789–7797. [Google Scholar] [CrossRef]
  10. Nakano, M.; Hara, S. Microprocessor-based repetitive control. Microprocess. Based Control Syst. 1986, 4, 279–296. [Google Scholar]
  11. Chew, K.-K.; Tomizuka, M. Steady-state and stochastic performance of a modified discrete-time prototype repetitive controller. Trans. ASME J. Dyn. Syst. Meas. Contr. 1990, 112, 35–41. [Google Scholar] [CrossRef]
  12. Hillerstrom, G.; Sternby, J. Repetitive control using low order models. In Proceedings of the American Control Conference, Baltimore, MD, USA, 29 June–1 July 1994; Volume 2, pp. 1873–1878. [Google Scholar]
  13. Grino, R.; Costa-Castello, R. Digital repetitive plug-in controller for odd-harmonic periodic references and disturbances. Automatica 2005, 41, 153–157. [Google Scholar] [CrossRef]
  14. Kurniawan, E.; Afandi, M.I.; Suryadi, S. Repetitive control system for tracking and rejection of multiple periodic signals. In Proceedings of the International Conference on Robotics, Automation and Sciences (ICORAS), Melaka, Malaysia, 27–29 November 2017; pp. 1–5. [Google Scholar]
  15. Steinbuch, M.; Weiland, S.; Singh, T. Design of Noise and Period-time Robust High-order Repetitive Control, with Application to Optical Storage. Automatica 2007, 43, 2086–2095. [Google Scholar] [CrossRef]
  16. Huang, W.-W.; Hu, C.; Zhu, L.-M. Robust Repetitive Control of Nanopositioning Stages Using the Spectrum-Selection Filter with Narrow Passbands. IEEE/ASME Trans. Mechatron. 2022, 27, 4211–4216. [Google Scholar] [CrossRef]
  17. Tian, M.; Wang, B.; Yu, Y.; Dong, Q.; Xu, D. Discrete-Time Repetitive Control-Based ADRC for Current Loop Disturbances Suppression of PMSM Drives. IEEE Trans. Ind. Inform. 2022, 18, 3138–3149. [Google Scholar] [CrossRef]
  18. Jamil, M.; Waris, A.; Gilani, S.O.; Khawaja, B.A.; Khan, M.N.; Raza, A. Design of Robust Higher-Order Repetitive Controller Using Phase Lead Compensator. IEEE Access 2020, 8, 30603–30614. [Google Scholar] [CrossRef]
  19. Lu, W.; Wang, W.; Zhou, K.; Fan, Q. General High-Order Selective Harmonic Repetitive Control for PWM Converters. IEEE J. Emerg. Sel. Top. Power Electron. 2022, 10, 1178–1191. [Google Scholar] [CrossRef]
  20. Duan, C.; Gu, G.; Du, C.; Chong, T.C. Robust Compensation of Periodic Disturbances by Multirate Control. IEEE Trans. Magn. 2008, 44, 413–418. [Google Scholar] [CrossRef]
  21. Yang, R.; Wang, M.; Li, L.; Wang, G.; Zhong, C. Robust Predictive Current Control of PMLSM With Extended State Modeling Based Kalman Filter: For Time-Varying Disturbance Rejection. IEEE Trans. Power Electron. 2020, 35, 2208–2221. [Google Scholar] [CrossRef]
  22. Landau, I.D.; Constantinescu, A.; Rey, D. Adaptive narrow band disturbance rejection applied to an active suspension—An internal model principle approach. Automatica 2005, 41, 563–574. [Google Scholar] [CrossRef]
  23. Zhao, R.; Xie, W.; Yu, G.; Wang, G.; Wong, P.K.; Silvestre, C. Adaptive ride height controller design for vehicle active suspension systems with uncertain sprung mass and time-varying disturbance. Int. J. Robust Nonlinear Control 2022, 32, 5950–5966. [Google Scholar] [CrossRef]
  24. Manayathara, T.J.; Tsao, T.-C.; Bentsman, J. Rejection of unknown periodic load disturbances in continuous steel casting process using learning repetitive control approach. IEEE Trans. Control Syst. Technol. 1996, 4, 259–265. [Google Scholar] [CrossRef]
  25. Ma, Z.; Fang, Y.; Zheng, H.; Liu, L. Active Disturbance Rejection Control With Self-Adjusting Parameters for Vibration Displacement System of Continuous Casting Mold. IEEE Access 2019, 7, 52498–52507. [Google Scholar] [CrossRef]
  26. Tomizuka, M. Zero phase error tracking algorithm for digital control. Trans. ASME J. Dyn. Syst. Meas. Contr. 1987, 109, 65–68. [Google Scholar] [CrossRef]
  27. Kurniawan, E.; Adinanta, H.; Harno, H.G.; Prakosa, J.A.; Suryadi, S.; Purwowibowo, P. On the synthesis of a stable and causal compensator for discrete-time high-order repetitive control systems. Int. J. Dynam. Control 2021, 9, 727–736. [Google Scholar] [CrossRef]
  28. Kurniawan, E.; Cao, Z.; Man, Z. Design of Robust Repetitive Control with Time-Varying Sampling Periods. IEEE Trans. Ind. Electron. 2014, 61, 2834–2841. [Google Scholar] [CrossRef]
  29. Xu, K.; Longman, R.W. Use of Taylor Expansions of The Inverse Model to Design FIR Repetitive Controllers. In Spaceflight Mechanics; Advances in the Astronautical Sciences Series; American Astronautical Society: Springfield, VA, USA, 2009; Volume 134, pp. 1073–1088. [Google Scholar]
  30. Zhang, B.; Wang, D.; Zhou, K.; Wang, Y. Linear Phase Lead Compensation Repetitive Control of a CVCF PWM Inverter. IEEE Trans. Ind. Electron. 2008, 55, 1595–1602. [Google Scholar] [CrossRef]
  31. Prasitmeeboon, P.; Longman, R.W. Repetitive Control Compensator Design for Frequency Response near Singularities. J. Astronaut. Sci. 2021, 68, 916–945. [Google Scholar] [CrossRef]
  32. Kurniawan, E. Robust Repetitive Control and Applications. Ph.D. Thesis, Swinburne University of Technology, Hawthorn, Australia, 2013. [Google Scholar]
Figure 1. Magnitude responses of the GMRC and the HOMRC.
Figure 1. Magnitude responses of the GMRC and the HOMRC.
Applsci 12 11746 g001
Figure 2. Block diagram of the plug−in DIMRC system.
Figure 2. Block diagram of the plug−in DIMRC system.
Applsci 12 11746 g002
Figure 3. Realization of the proposed DIMRC.
Figure 3. Realization of the proposed DIMRC.
Applsci 12 11746 g003
Figure 4. Realization of the weighted multiple delays.
Figure 4. Realization of the weighted multiple delays.
Applsci 12 11746 g004
Figure 5. Realization of the DIMRC with non−causal stabilizing controller.
Figure 5. Realization of the DIMRC with non−causal stabilizing controller.
Applsci 12 11746 g005
Figure 6. Repetitive reference r ( k ) .
Figure 6. Repetitive reference r ( k ) .
Applsci 12 11746 g006
Figure 7. Magnitude and phase responses of the q−filter (32).
Figure 7. Magnitude and phase responses of the q−filter (32).
Applsci 12 11746 g007
Figure 8. Tracking errors for the minimum phase system (30) with the disturbance model v 1 ( t ) .
Figure 8. Tracking errors for the minimum phase system (30) with the disturbance model v 1 ( t ) .
Applsci 12 11746 g008
Figure 9. Tracking errors for the minimum phase system (30) with the disturbance model v 2 ( t ) .
Figure 9. Tracking errors for the minimum phase system (30) with the disturbance model v 2 ( t ) .
Applsci 12 11746 g009
Figure 10. Magnitude responses of the simulated DIMRC and MPRC.
Figure 10. Magnitude responses of the simulated DIMRC and MPRC.
Applsci 12 11746 g010
Figure 11. Tracking errors for the non-minimum phase system (42) with the disturbance model v 1 ( t ) .
Figure 11. Tracking errors for the non-minimum phase system (42) with the disturbance model v 1 ( t ) .
Applsci 12 11746 g011
Figure 12. Tracking errors for the non-minimum phase system (42) with the disturbance model v 2 ( t ) .
Figure 12. Tracking errors for the non-minimum phase system (42) with the disturbance model v 2 ( t ) .
Applsci 12 11746 g012
Table 1. Steady-state performance of the DIMRC and MPRC.
Table 1. Steady-state performance of the DIMRC and MPRC.
Minimum Phase SystemNon-Minimum Phase System
Method v 1 ( t ) v 2 ( t ) v 1 ( t ) v 2 ( t )
t ss rms- e ss t ss rms- e ss t ss rms- e ss t ss rms- e ss
DIMRC5.380.00225.390.002613.930.003913.880.0046
MPRC3.6250.00084.140.00297.1950.00197.1350.0065
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Prakosa, J.A.; Purwowibowo, P.; Kurniawan, E.; Wijonarko, S.; Maftukhah, T.; Sari, F.; Pratiwi, E.B.; Rustandi, D. Discrete-Time Design of Dual Internal Model-Based Repetitive Control Systems. Appl. Sci. 2022, 12, 11746. https://doi.org/10.3390/app122211746

AMA Style

Prakosa JA, Purwowibowo P, Kurniawan E, Wijonarko S, Maftukhah T, Sari F, Pratiwi EB, Rustandi D. Discrete-Time Design of Dual Internal Model-Based Repetitive Control Systems. Applied Sciences. 2022; 12(22):11746. https://doi.org/10.3390/app122211746

Chicago/Turabian Style

Prakosa, Jalu A., Purwowibowo Purwowibowo, Edi Kurniawan, Sensus Wijonarko, Tatik Maftukhah, Farakka Sari, Enggar B. Pratiwi, and Dadang Rustandi. 2022. "Discrete-Time Design of Dual Internal Model-Based Repetitive Control Systems" Applied Sciences 12, no. 22: 11746. https://doi.org/10.3390/app122211746

APA Style

Prakosa, J. A., Purwowibowo, P., Kurniawan, E., Wijonarko, S., Maftukhah, T., Sari, F., Pratiwi, E. B., & Rustandi, D. (2022). Discrete-Time Design of Dual Internal Model-Based Repetitive Control Systems. Applied Sciences, 12(22), 11746. https://doi.org/10.3390/app122211746

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop