Next Article in Journal
Evaluation of Hydrodynamic Performance of New Random Packing Structure Using CFD
Next Article in Special Issue
The Design and the Development of a Biped Robot Cooperation System
Previous Article in Journal
Green Chemistry Metrics, A Review
Previous Article in Special Issue
Solving the Formation and Containment Control Problem of Nonlinear Multi-Boiler Systems Based on Interval Type-2 Takagi–Sugeno Fuzzy Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Fuzzy System Based Iterative Learning Control for Nonlinear Discrete-Time Systems with Iteration-Varying Uncertainties

Department of Technology for Smart Living, Huafan University, New Taipei City 223, Taiwan
*
Author to whom correspondence should be addressed.
Processes 2022, 10(7), 1275; https://doi.org/10.3390/pr10071275
Submission received: 19 June 2022 / Revised: 25 June 2022 / Accepted: 26 June 2022 / Published: 29 June 2022
(This article belongs to the Special Issue Application of Fuzzy Control in Computational Intelligence)

Abstract

:
In this paper, we consider an iterative learning control problem for a class of unknown discrete-time nonlinear systems with iteration-varying initial error, iteration-varying system parameters, iteration-varying external disturbance, iteration-varying desired output, and iteration-varying control direction. These iteration-varying uncertainties are not required to take any particular structure such as the high-order internal model and only need to satisfy certain boundedness conditions. We propose an iterative learning control law with an adaptive iteration-varying fuzzy system to overcome all the uncertainties and achieve the learning control objective. Furthermore, we present a sufficient condition for designing adaptive gains and prove the convergence of the learning error to a small value as the trial number is large enough. Finally, we use two simulation examples to demonstrate all the theoretical results.

1. Introduction

For controlling a system with a given repetitive task over a finite time interval, we often considered using the iterative learning control (ILC) [1,2] in the past three decades. Being an interesting and attractive learning strategy, ILC has achieved plenty of research works in the literature. Several important theoretical results [3,4,5] and practical applications [6,7,8,9,10,11] can still be found in recent years. Basically, the ILC laws use the previous input and output data stored in the memory to construct a current control input for the system to track the desired output perfectly as the number of learning trials is large enough. A basic requirement of the control system when designing the ILC laws is that all the environmental conditions, such as the system parameters or external disturbances, are invariant from trial to trial. This requirement seems reasonable but is generally restrictive as most conditions may vary in the real physical world.
There are some possible iteration-varying uncertain sources existing in the iterative learning control systems. The most important ones include: (a) iteration-varying initial error, (b) iteration-varying external disturbance, (c) iteration-varying desired output, (d) iteration-varying system parameters, and (e) iteration-varying control direction (sign of control gain). To solve the design of ILC laws under the consideration of iteration-varying uncertainties is in general difficult. Some research works have been proposed in the area of ILC. For example, the ILC laws were studied under the uncertainty (a) in [12,13], under uncertainty (c) in [14], under uncertainties (a), (c) in [15,16], under uncertainties (a), (b), (c) in [17,18], under uncertainties (c), (d) in [19], under uncertainties (a), (c), (e) in [20], under uncertainties (a), (c), (d), (e) in [21], under uncertainties (a),(b), (c), (d) in [22,23,24], and under uncertainty (e) in [25,26]. However, to the best of the authors’ knowledge, no research works have been published to deal with the design and analysis of ILC laws under all the uncertainties of (a), (b), (c), (d), (e) at the same time.
The motivation of this paper is to find a realizable iterative learning control law for a class of discrete-time nonlinear systems under the five iteration-varying uncertainties mentioned above. Furthermore, this paper still has some more new contributions when compared with the existing ILC laws when dealing with iteration-varying uncertainties. Firstly, the uncertainties do not have to take any particular forms, such as the strict high-order internal model (HOIM) structure. Secondly, the uncertainties are not required to converge to certain iteration-invariant ones along the iteration domain in order to guarantee a satisfying learning performance. Thirdly, the upper bounds of these uncertainties can be unknown and not necessarily small.
To achieve this paper’s iterative learning control objective, we first apply a fuzzy system with iteration-varying consequent parameters as an approximator to approximate the unknown nonlinear function with iteration-varying parameters. Based on this iteration-varying fuzzy system, an adaptive strategy is adopted to design the ILC law to solve the ILC problem with all the existing iteration-varying uncertainties. The control structure is simple and easy to implement. We first reformulate the error dynamics into an iteration-invariant nominal part and an iteration-varying uncertain part, respectively. After a further derivation, the error dynamics becomes a parameterized linear combination of fuzzy basis function vector and control input as well as a lumped iteration-varying uncertainty. An iteration-varying dead-zone-like auxiliary error is proposed to deal with the lumped iteration-varying uncertainty. Then a set of adaptive laws is presented using the auxiliary error to update the fuzzy and control parameters from trial to trial in this ILC law. A rigorous technical proof is given to guarantee the boundedness of the learning system for each iteration and discrete-time instant. Furthermore, we show that the norm of output error can asymptotically converge along the iteration axis to a value that depends on the width of the dead zone. The main contributions can be summarized as follows.
(1)
This is the first work that considers the design and analysis of iterative learning control law for nonlinear unknown systems with all the five kinds of iteration-varying uncertainties.
(2)
The mentioned iteration-varying uncertainties are allowed to be unknown and without any special structure.
(3)
The upper bounds of the iteration-varying uncertainties can be unknown and not necessarily small.
(4)
A new concept of using a fuzzy system with iteration-varying consequent parameters as an approximator is proposed.
This paper is organized as follows. Section 2 gives a problem formulation and definition of the iteration-varying uncertainties for this work. The iterative learning control law and the parameter adaptive laws are proposed in Section 3. Based on the derived error dynamics, we analyzed the closed-loop stability and learning performance in Section 4. In Section 5, two simulation examples are then given to show the effectiveness of the iterative learning control law. Finally, we made a conclusion in Section 6.

2. Problem Formulation

In this paper, we consider a class of repetitive nonlinear discrete-time systems as follows,
y j ( k + 1 ) = f ( a j ( k ) , X j ( k ) ) + b j ( k ) u j ( k ) + d j ( k )
where k { 0 , 1 , 2 , , N } is the index of discrete time, j Z + is the index of learning control iteration, X j ( k ) = [ y j ( k ) , , y j ( k 1 + n ) ] R n is the state vector, f ( a j ( k ) , X j ( k ) ) R is an unknown nonlinear function of a j ( k ) and X j ( k ) . y j ( k ) R is the output, u j ( k ) R is the input and without loss of generality, we let y j ( k ) = 0 , k < 0 . The iteration-varying uncertainties appeared in this system are defined as follows: (1) the initial output y j ( 0 ) is iteration-varying, (2) the system parameters a j ( k ) R n and b j ( k ) R 1 are unknown and iteration-varying, (3) the disturbance d j ( k ) R is unknown and iteration-varying, (4) the desired output y d j ( k ) is iteration-varying, (5) the sign of control gain b j ( k ) is unknown and iteration-varying. The control objective is to design an iterative learning control law u j ( k ) such that y j ( k ) will track y d j ( k ) as close as possible, k { 1 , 2 , , N } , even all the above-mentioned uncertainties exist. Several assumptions are required for the controller design as follows.
(A1)
The iteration-varying initial output error e j ( 0 ) = y j ( 0 ) y d j ( 0 ) is bounded j 1 .
(A2)
The iteration-varying system parameters a j ( k ) , b j ( k ) are bounded k { 0 , 1 , 2 , , N } and j 1 .
(A3)
The iteration-varying disturbance d j ( k ) is bounded with | d j ( k ) | d U where d U is an unknown positive constant for all k { 0 , 1 , 2 , , N } and j 1 .
(A4)
The iteration-varying desired output y d j ( k ) is bounded with | y d j ( k ) | y d U where y d U is an unknown positive constant for all k { 0 , 1 , 2 , , N } and j 1 .
(A5)
The unknown nonlinear function f ( a j ( k ) , X j ( k ) ) is smooth and bounded if a j ( k ) and X j ( k ) are bounded.
Remark 1.
Note that all the iteration-varying uncertainties are required to be bounded. This is reasonable in a real control environment. But the upper bounds of the uncertainties are not required to be small enough so as to guarantee a satisfied performance. Furthermore, these upper bounds are allowed to be unknown. We will design suitable adaptive laws to overcome these unknown bounds. Compared with all the related works in the area of ILC dealing with similar problems, this paper can handle the more general class of iteration-varying uncertainties and needs less knowledge on them.

3. Design of the Iterative Learning Control Law

To overcome the unknown smooth nonlinear function f ( a j ( k ) , X j ( k ) ) , we proposed a fuzzy system to approximate f ( a j ( k ) , X j ( k ) ) as most of the design approaches used in the past two decades in the literature. However, we are going to present a different concept to the fuzzy system for an optimal function approximation since the system parameter vector a j ( k ) is iteration-varying. We suggest that there exists an optimal iteration-varying fuzzy system Θ j ( k ) ξ ( X j ( k ) ) for approximation of f ( a j ( k ) , X j ( k ) ) , where Θ j ( k ) R m is the iteration-varying fuzzy consequent parameter vector, and ξ j ( k ) R m is the fuzzy basis function vector. This iteration-varying fuzzy system will in general gives a better approximation to the iteration-varying nonlinear function f ( a j ( k ) , X j ( k ) ) , and the approximation error over a compact set will be smaller or the compact set for a reasonable approximation will be larger than those traditional designs using fixed consequent parameter fuzzy systems. More precisely, the approximation behavior between the optimal fuzzy system Θ j ( k ) ξ ( X j ( k ) ) and the nonlinear function f ( a j ( k ) , X j ( k ) ) will satisfy
f ( a j ( k ) , X j ( k ) ) = Θ j ( k ) ξ ( X j ( k ) + ε ( X j ( k ) )
where | ε ( X j ( k ) ) | ε U on a compact set A R n with ε U being an unknown positive constant.
Then we can rewrite (1) by using (2) as
y j ( k + 1 ) = Θ j ( k ) ξ ( X j ( k ) + b j ( k ) u j ( k ) + d j ( k ) + ε ( X j ( k ) )
We now separate the iteration-varying parameters Θ j ( k ) and b j ( k ) into the iteration-invariant nominal parts Θ * ( k ) , b * ( k ) , and the iteration-varying uncertain parts Θ ¯ j ( k ) , b ¯ j ( k ) , respectively. That is, Θ j ( k ) = Θ * ( k ) + Θ ¯ j ( k ) and b j ( k ) = b * ( k ) + b ¯ j ( k ) .
Remark 2.
It is reasonable that Θ * ( k ) , b * ( k ) , and Θ ¯ j ( k ) , b ¯ j ( k ) are all bounded according to assumption (A2). Furthermore, we assume | Θ ¯ j ( k ) | θ U and | b ¯ j ( k ) | b U for some positive unknown constants θ U and b U .
Based on (3), we derive the output error e j ( k + 1 ) as follows,
e j ( k + 1 ) = y j ( k + 1 ) y d j ( k + 1 ) = Θ j ( k ) ξ ( X j ( k ) ) + b j ( k ) u j ( k ) + d j ( k ) + ε ( X j ( k ) ) y d j ( k + 1 ) = Θ * ( k ) ξ ( X j ( k ) ) + b * ( k ) u j ( k ) + Θ ¯ j ( k ) ξ ( X j ( k ) ) + b ¯ j ( k ) u j ( k ) + d j ( k ) + ε ( X j ( k ) ) y d j ( k + 1 )
The iterative learning control law is designed as
u j ( k ) = b ^ j ( k ) λ + b ^ j ( k ) 2 Θ ^ j ( k ) ξ ( X j ( k ) ) + y d j ( k + 1 )
where λ is a small positive constant, Θ ^ j ( k ) and b ^ j ( k ) are adaptive control parameters which will be updated along the iteration domain in order to ensure the learning error convergence and closed-loop stability. Define the control parameter errors as Θ ˜ j ( k ) = Θ ^ j ( k ) Θ * ( k ) , b ˜ j ( k ) = b ^ j ( k ) b * ( k ) and substitute the iterative learning control law (5) into the output error (4). Then we have
e j ( k + 1 ) = Θ * ( k ) ξ ( X j ( k ) ) Θ ^ j ( k ) ξ ( X j ( k ) ) + Θ ^ j ( k ) ξ ( X j ( k ) ) + b * ( k ) u j ( k ) b ^ j ( k ) u j ( k ) + b ^ j ( k ) u j ( k ) + Θ ¯ j ( k ) ξ ( X j ( k ) ) + b ¯ j ( k ) u j ( k ) + d j ( k ) + ε ( X j ( k ) ) y d j ( k + 1 ) = Θ ˜ j ( k ) ξ ( X j ( k ) ) b ˜ j ( k ) u j ( k ) + λ λ + b ^ j ( k ) 2 Θ ^ j ( k ) ξ ( X j ( k ) ) y d j ( k + 1 ) + Θ ¯ j ( k ) ξ ( X j ( k ) ) + b ¯ j ( k ) u j ( k ) + d j ( k ) + ε ( X j ( k ) ) Θ ˜ j ( k ) ξ ( X j ( k ) ) b ˜ j ( k ) u j ( k ) + δ j ( k )
where
δ j ( k ) = Θ ¯ j ( k ) ξ ( X j ( k ) ) + b ¯ j ( k ) u j ( k ) + d j ( k ) + ε ( X j ( k ) ) + λ λ + b ^ j ( k ) 2 Θ ^ j ( k ) ξ ( X j ( k ) ) y d j ( k + 1 ) = Θ ¯ j ( k ) ξ ( X j ( k ) ) + d j ( k ) + ε ( X j ( k ) ) + b ¯ j ( k ) b ^ j ( k ) λ + b ^ j ( k ) 2 Θ ^ j ( k ) ξ ( X j ( k ) ) y d j ( k + 1 ) + λ λ + b ^ j ( k ) 2 Θ ^ j ( k ) ξ ( X j ( k ) ) y d j ( k + 1 )
It is easily shown that
b ¯ j ( k ) b ^ j ( k ) λ + b ^ j ( k ) 2 b U 2 λ
so that the term δ j ( k ) in (7) will satisfy the bounding condition as follows,
| δ j ( k ) | θ U | ξ ( X j ( k ) ) | + d U + ε U + b U 2 λ + 1 m | Θ ^ j ( k ) | + y d U ψ * R j ( k )
where ψ * is a suitably defined unknown positive constant and R j ( k ) = | Θ ^ j ( k ) | + 1 . To overcome all the unknown Θ * ( k ) , b * ( k ) , ψ * caused by iteration-varying uncertainties, a dead-zone like auxiliary error e ϕ j ( k + 1 ) is firstly introduced:
e ϕ j ( k + 1 ) = e j ( k + 1 ) ϕ j ( k + 1 ) s a t e j ( k + 1 ) ϕ j ( k + 1 )
for k { 0 , 1 , 2 , , N 1 } . It is noted that e ϕ j ( 0 ) is not defined since it will not be utilized in the later design. The notation sat in (10) is a saturation function defined as
s a t e j ( k + 1 ) ϕ j ( k + 1 ) = 1 i f e j ( k + 1 ) > ϕ j ( k + 1 ) e j ( k + 1 ) ϕ j ( k + 1 ) i f | e j ( k + 1 ) | ϕ j ( k + 1 ) 1 i f e j ( k + 1 ) < ϕ j ( k + 1 )
ϕ j ( k + 1 ) can be considered as a width of dead zone which is iteration-time varying and it is defined as
ϕ j ( k + 1 ) = ψ ^ j ( k ) R j ( k )
where ψ ^ j ( k ) is another control parameter. According to the definition, it is easy to prove that e ϕ j ( k + 1 ) s a t e j ( k + 1 ) ϕ j ( k + 1 ) = | e ϕ j ( k + 1 ) | .
In this iterative learning controller, Θ ^ j ( k ) , b ^ j ( k ) in (5) and ψ ^ j ( k ) in (11) are used to compensate for the unknown iteration-invariant parameters Θ * ( k ) , b * ( k ) and ψ * , respectively. As the uncertainties are unknown, a direct strategy is to construct certain adaptive laws to automatically search for the optimal parameters. The adaptive laws are designed to ensure the closed-loop stability and improve the learning performance as the number of iterations is large enough. The parameter adaptive laws for Θ ^ j ( k ) , b ^ j ( k ) and ψ ^ j ( k ) at ( j + 1 )th iteration are proposed in the following:
Θ ^ j + 1 ( k ) = Θ ^ j ( k ) + β 1 e ϕ j ( k + 1 ) ξ ( X j ( k ) ) 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2
b ^ j + 1 ( k ) = b ^ j ( k ) + β 2 e ϕ j ( k + 1 ) u j ( k ) 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2
ψ ^ j + 1 ( k ) = ψ ^ j ( k ) + β 3 | e ϕ j ( k + 1 ) | R j ( k ) 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2
for k { 0 , 1 , , N 1 } , where β 1 , β 2 , β 3 > 0 are the adaptive gains. In the first trial, we will set Θ ^ 1 ( k ) = Θ 1 , b ^ 1 ( k ) = b 1 to be any constant vector, and ψ ^ 1 ( k ) = ψ 1 > 0 to be a small number. According to (14), we have ψ ^ j ( k ) > 0 , k { 0 , 1 , , N 1 } and j 1 . If we define the parameter errors as Θ ˜ j ( k ) = Θ ^ j ( k ) Θ * ( k ) , b ˜ j ( k ) = b ^ j ( k ) b * ( k ) , ψ ˜ j ( k ) = ψ ^ j ( k ) ψ * and subtract Θ * ( k ) , b * ( k ) and ψ * on both sides of (12)–(14) respectively, then it will yield
Θ ˜ j + 1 ( k ) = Θ ˜ j ( k ) + β 1 e ϕ j ( k + 1 ) ξ ( X j ( k ) ) 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2
b ˜ j + 1 ( k ) = b ˜ j ( k ) + β 2 e ϕ j ( k + 1 ) u j ( k ) 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2
ψ ˜ j + 1 ( k ) = ψ ˜ j ( k ) + β 3 | e ϕ j ( k + 1 ) | R j ( k ) 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2

4. The Main Results for Stability and Convergence

This section will study the system convergent performance and the boundedness of all the internal signals. The main results are summarized in the following theorem.
Theorem 1.
Consider the class of discrete-time nonlinear systems (1) which satisfies the assumptions (A1)–(A5) and design iterative learning control law as in (5), (10) and (12)–(14). If β = max { β 1 , β 2 , β 3 } satisfies
2 β > 0 ,
then we can get the following results.
(t1) 
The adaptive parameters Θ ^ j ( k ) , b ^ j ( k ) , ψ ^ j ( k ) are bounded k { 0 , 1 , , N } , j 1 .
(t2) 
The dead-zone-like auxiliary error e ϕ j ( k + 1 ) , output error e j ( k ) and control input u j ( k ) are bounded k { 0 , 1 , , N 1 } , j 1 . Furthermore,
lim j e ϕ j ( k + 1 ) = 0 lim j | e j ( k + 1 ) | ψ ^ ( k ) R ( k ) .
Proof. 
(t1) To show the boundedness of adaptive parameters, we define a positive function as follows,
V j ( k ) = 1 β 1 Θ ˜ j ( k ) Θ ˜ j ( k ) + 1 β 2 b ˜ j ( k ) 2 + 1 β 3 ψ ˜ j ( k ) 2
By using (15)–(17), we can derive V j + 1 ( k ) V j ( k ) as follows,
V j + 1 ( k ) V j ( k ) = 1 β 1 Θ ˜ j + 1 ( k ) Θ ˜ j + 1 ( k ) Θ ˜ j ( k ) Θ ˜ j ( k ) + 1 β 2 b ˜ j + 1 ( k ) 2 b ˜ j ( k ) 2 + 1 β 3 ψ ˜ j + 1 ( k ) 2 ψ ˜ j ( k ) 2 2 e ϕ j ( k + 1 ) Θ ˜ j ( k ) ξ ( X j ( k ) ) 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 + β 1 e ϕ j ( k + 1 ) 2 | ξ ( X j ( k ) ) | 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 2 + 2 e ϕ j ( k + 1 ) b ˜ j ( k ) u j ( k ) 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 + β 2 e ϕ j ( k + 1 ) 2 u j ( k ) 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 2 + 2 | e ϕ j ( k + 1 ) | ψ ˜ j ( k ) R j ( k ) 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 + β 3 e ϕ j ( k + 1 ) 2 R j ( k ) 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 2
According to (6), we can find that
Θ ˜ j ( k ) ξ ( X j ( k ) ) + b ˜ j ( k ) u j ( k ) = e j ( k + 1 ) + δ j ( k )
This implies that
e ϕ j ( k + 1 ) Θ ˜ j ( k ) ξ ( X j ( k ) ) + e ϕ j ( k + 1 ) b ˜ j ( k ) u j ( k ) = e j ( k + 1 ) e ϕ j ( k + 1 ) + e ϕ j ( k + 1 ) δ j ( k )
If we substitute (21) into (19), then we have
V j + 1 ( k ) V j ( k ) 2 e j ( k + 1 ) e ϕ j ( k + 1 ) + 2 e ϕ j ( k + 1 ) δ j ( k ) 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 + β 1 e ϕ j ( k + 1 ) 2 | ξ ( X j ( k ) ) | 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 2 + β 2 e ϕ j ( k + 1 ) 2 u j ( k ) 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 2 + 2 | e ϕ j ( k + 1 ) | ψ ˜ j ( k ) R j ( k ) 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 + β 3 e ϕ j ( k + 1 ) 2 R j ( k ) 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 2
Substituting (10) into (22) and using the fact that | δ j ( k ) | ψ * R j ( k ) in (9), we can derive that
V j + 1 ( k ) V j ( k ) 2 e ϕ j ( k + 1 ) 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 2 | e ϕ j ( k + 1 ) | ψ ^ j ( k ) R j ( k ) 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 + 2 | e ϕ j ( k + 1 ) | ψ * R j ( k ) 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 + 2 | e ϕ j ( k + 1 ) | ψ ˜ j ( k ) R j ( k ) 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 + β 1 e ϕ j ( k + 1 ) 2 | ξ ( X j ( k ) ) | 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 2 + β 2 e ϕ j ( k + 1 ) 2 u j ( k ) 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 2 + β 3 e ϕ j ( k + 1 ) 2 R j ( k ) 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 2 2 e ϕ j ( k + 1 ) 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 + β e ϕ j ( k + 1 ) 2 | ξ ( X j ( k ) ) | 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 2 + β e ϕ j ( k + 1 ) 2 u j ( k ) 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 2 + β e ϕ j ( k + 1 ) 2 R j ( k ) 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 2 ( 2 β ) e ϕ j ( k + 1 ) 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2
where β = max { β 1 , β 2 , β 3 } . If we choose β such that p 2 β > 0 , then we have
V j + 1 ( k ) V j ( k ) p e ϕ j ( k + 1 ) 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 0
for j 1 . Note that for the first iteration, Θ ˜ 1 ( k ) = Θ ^ 1 ( k ) Θ * ( k ) , b ˜ 1 ( k ) = b ^ 1 ( k ) b * ( k ) and ψ ˜ 1 ( k ) = ψ ^ 1 ( k ) ψ * are bounded k { 0 , 1 , 2 , , N } since the initial settings of Θ ^ 1 ( k ) , b ^ 1 ( k ) and ψ ^ 1 ( k ) are bounded. This implies V 1 ( k ) is bounded k { 0 , 1 , 2 , , N } , and hence V j ( k ) , Θ ˜ j ( k ) , b ˜ j ( k ) and ψ ˜ j ( k ) are bounded j 1 according to the result of (24). This proves (t1) of the theorem.
(t2) Summing (24) from 1 to j, it yields
V j ( k ) V 1 ( k ) i = 1 j 1 p e ϕ i ( k + 1 ) 2 1 + | ξ ( X i ( k ) ) | 2 + u i ( k ) 2 + R i ( k ) 2
Because V 1 ( k ) is bounded and V j ( k ) is nonnegative, we can conclude that
lim j e ϕ j ( k + 1 ) 2 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 = 0
k { 0 , 1 , 2 , , N 1 } .
To prove that e ϕ j ( k + 1 ) will be bounded and converge to zero, we firstly note that the fuzzy basis function vector ξ ( X j ( k ) ) is bounded, and u j ( k ) given in (5) as well as R j ( k ) = | Θ ^ j ( k ) | + 1 are bounded according to (t1). We can guarantee 1 + | ξ ( X j ( k ) ) | 2 + u j ( k ) 2 + R j ( k ) 2 are bounded j 1 . Hence, (25) implies that
lim j e ϕ j ( k + 1 ) 2 = 0
j 1 and k { 0 , , N 1 } . When j , the output error e ( k + 1 ) satisfies
lim j | e j ( k + 1 ) | ϕ ( k + 1 ) = ψ ^ ( k ) R ( k )
k { 0 , 1 , 2 , , N 1 } . This concludes the result of (t2) in the theorem. □
Remark 3.
In this theorem, we show that the output error e j ( k + 1 ) will converge to a residual set which is bounded by ψ ^ ( k ) R ( k ) . Therefore, we require that ψ ^ ( k ) a ( k ) is as small as possible for all k { 0 , 1 , 2 , , N } . This is the reason for us to set the initial value of ψ ^ 1 ( k ) to be a small number and the adaptive gain β 3 in (14) as a small one such that ψ ^ j ( k ) and hence, ψ ^ ( k ) R ( k ) , k { 0 , 1 , 2 , , N } will remain in a reasonable small value for all j 1 . Of course it is no problem to choose β 3 as small as possible due to the requirement of convergent condition in (18).

5. Simulation Examples

Example 1.
In this example, we use a simple parameterized nonlinear system with iteration-varying uncertainties to study the learning performance of the iterative learning control system. The nonlinear system considered is given as
y j ( k + 1 ) = Θ j ( k ) f ( X j ( k ) ) + b j ( t ) u j ( k ) + d j ( t )
where
Θ j ( k ) = 2 + 0.5 sin ( k ) 0.1 sin ( 2 π j / 100 )
f ( X j ( k ) ) = 2 sin ( y j ( k ) ) 1 + 0.01 y j ( k 1 ) 2
b j ( k ) = s ( 1 + 0.3 e j k + 0.3 sin ( k π / 100 ) )
d j ( k ) = 0.1 r a n d
y j ( 0 ) = 0.5 + 0.1 r a n d
here s denotes the sign of control gain b j ( k ) , rand is a uniform distribution on the interval ( 0 , 1 ) and f ( X j ( k ) ) is known for simplicity. We use the following reference model to generate the desired output,
y d j ( k + 1 ) = 0.3 y d j ( k ) + r j ( k ) , y d j ( 0 ) = 0.1 + 0.1 r a n d r j ( k ) = 3 sin ( 2 π k / 15 ) + 0.1 ( 1 + 0.1 sin ( j / 20 ) )
The iterative learning control objective in this example is to force the system output y j ( k ) to track the desired output y d j ( k ) for all k { 1 , , 50 } as close as possible. The proposed iterative learning control law (5), (10) and (12)–(14) are chosen with ξ ( X j ( k ) ) = f ( X j ( k ) ) , λ = 0.01 , β 1 = 1.5 , β 2 = 0.5 , β 3 = 0.001 so that 2 β = 2 max { β 1 , β 2 , β 3 } = 0.5 > 0 . The control parameter values at first iteration are given as Θ ^ 1 ( k ) = 1.5 , b ^ 1 ( k ) = 0.2 , ψ ^ 1 ( k ) = 0.01 for all k { 0 , 1 , , 50 } .
Case 1: s = 1 (unknown positive control gain)
Figure 1a shows max k { 1 , , 50 } | e ϕ j ( k ) | versus iteration j = 1 , , 100 . The simulation result proves that (t2) of the theorem is correct. Figure 1b shows the output error e 100 ( k ) , width of iteration-varying dead zone ϕ 100 ( k ) and ϕ 100 ( k ) at the 100th trial. The profile of e 100 ( k ) roughly lies between ϕ 100 ( k ) and ϕ 100 ( k ) , for k { 1 , , 50 } . The system output y 100 ( k ) and desired output y d 100 ( k ) at the 100th iteration are shown in Figure 1c. Finally, the bounded control input u 100 ( k ) at the 100th iteration is given in Figure 1d.
Case 2: s = 1 (unknown negative control gain)
In order to check if the learning performance still remains when the sign of the control gain b j ( k ) changes, we let s of b j ( k ) in (31) be replaced by s = 1 and make a simulation again without changing any parameters in the proposed iterative learning control law in case 1. It is clear from Figure 2 that the learning performances are almost the same as those in case 1. But the control input u 100 ( k ) in Figure 2d is reversed when compared with that in Figure 1d. This proves that the proposed iterative learning control law can work without the knowledge of the control gain sign.
Case 3: s = 1 at the 15th, 30th, 60th iterations and s = 1 for the others (iteration-varying control gain)
In this case, we study the learning effect when the sign of control gain b j ( k ) is iteration-varying. We let s = 1 for most of the iterations except at the 15th, 30th, and 60th iterations. The simulation results are now shown in Figure 3. It is interesting to find that the learning system has to re-learn when the sign of control gain suddenly changes. However, it can overcome the variation quickly after several trials and capture the control direction again.
Example 2.
In this example, we use the proposed fuzzy iterative learning control law to control a non-BIBO nonlinear unknown system similar to but more complex than that in [27]. The difference equation of the nonlinear system with iteration-varying uncertainties is given as
y j ( k + 1 ) = 0.2 + 0.1 sin k / j y j ( k ) 2 + 0.2 y j ( k 1 ) + 0.4 sin 0.5 ( 1 e j k 100 ) y j ( k 1 ) + y j ( k ) cos 0.5 y j ( k 1 ) + y j ( k ) + 1.2 + 0.3 cos ( y j ( k ) u j ( k ) + d j ( k ) , y j ( 0 ) = 0.5 + 0.1 r a n d n
where y j ( k ) is the system output, u j ( k ) is the control input and d j ( k ) = 0.01 r a n d is a non-repeatable random disturbance. Here the reference model is chosen as
y d j ( k + 1 ) = 0.6 y d j ( k ) + r j ( k ) , y d j ( 0 ) = 0
where r j ( k ) = 0.2 sin ( 2 π k / 100 ) + 0.1 sin ( 2 π j / 50 ) is an iteration-varying bounded reference input. We choose λ = 0.1 in (5) and the initial value of control parameters at the first iteration as Θ ^ 1 ( k ) = [ 0.1 , 0.1 , 0.1 , 0.1 , 0.1 ] , b ^ 1 ( k ) = 1 and ψ ^ 1 ( k ) = 0.0015 , respectively. The adaptive gains are chosen as β 1 = 1.5 , β 2 = 1 , β 3 = 0.01 . Five fuzzy rules are chosen to construct the fuzzy basis function vector ξ ( y j ( k ) , y j ( k 1 ) ) . The detailed fuzzy membership functions are given as follows,
μ 11 ( y j ( k ) ) = 1 i f y j ( k ) < 2 e ( y j ( k ) ( 2 ) ) 2 i f y j ( k ) 2 μ 12 ( y j ( k ) ) = e ( y j ( k ) ( 1 ) ) 2 μ 13 ( y j ( k ) ) = e ( y j ( k ) ( 0 ) ) 2 μ 14 ( y j ( k ) ) = e ( y j ( k ) ( 1 ) ) 2 μ 15 ( y j ( k ) ) = e ( y j ( k ) ( 2 ) ) 2 i f y j ( k ) < 2 1 i f y j ( k ) 2 μ 21 ( y j ( k 1 ) ) = 1 i f y j ( k 1 ) < 2 e ( y j ( k 1 ) ( 2 ) ) 2 i f y j ( k 1 ) 2 μ 22 ( y j ( k 1 ) ) = e ( y j ( k 1 ) ( 1 ) ) 2 μ 23 ( y j ( k 1 ) ) = e ( y j ( k 1 ) ( 0 ) ) 2 μ 24 ( y j ( k 1 ) ) = e ( y j ( k 1 ) ( 1 ) ) 2 μ 25 ( y j ( k 1 ) ) = e ( y j ( k 1 ) ( 2 ) ) 2 i f y j ( k 1 ) < 2 1 i f y j ( k 1 ) 2
Hence, the fuzzy basis functions are given as
ξ i ( y j ( k ) , y j ( k 1 ) ) = μ 1 i ( y j ( k ) ) μ 2 i ( y j ( k 1 ) )
i = 1 , 2 , 3 , 4 , 5 . In Figure 4a, we show the evolution of max k { 1 , , 200 } | e ϕ j ( k ) | with respective to iteration j. After 100 iterations, the converged learning error e 100 ( k ) is shown in Figure 4b. Same as example 1, the trajectory of e 100 ( k ) roughly lies between ϕ 100 ( k ) and ϕ 100 ( k ) , for k { 1 , , 200 } . After the 100th trial, a satisfied output tracking performance is achieved. We then show the relation between y 100 ( k ) and y d 100 ( k ) in Figure 4c for k { 0 , 1 , 2 , , 200 } . Finally, Figure 4d shows the control input u 100 ( k ) which is clearly bounded.

6. Conclusions

A robustness problem caused by iteration-varying uncertainties in studying iterative learning control is always a practical but difficult challenge. This paper aims to design an iterative learning control law for a class of unknown discrete-time nonlinear systems with five kinds of iteration-varying uncertainties. These uncertainties include iteration-varying initial error, iteration-varying external disturbance, iteration-varying desired output, iteration-varying system parameters, and iteration-varying control direction. This paper is the first ILC work that can deal with these uncertainties at the same time. In addition to this main contribution, the five kinds of iteration-varying uncertainties can take a general structure. They don’t have to satisfy the assumption of the high-order internal model or any other special formulation. We not only show that all the internal signals are bounded for each iteration and discrete-time instant but also guarantee that the output error asymptotically converges to a small value. Two simulation results are studied with several scenarios to show that the fuzzy iterative learning control law is practical and feasible.

Author Contributions

Conceptualization, C.-J.C., Y.-C.W.; methodology, C.-J.C.; software, Y.-C.W.; validation, C.-J.C., Y.-C.W.; formal analysis, C.-J.C.; writing–original draft preparation, C.-J.C.; writing–review and editing, C.-J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by Ministry of Science and Technology, Taiwan, under Grants MOST110-2221-E-211-002.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bristow, D.A.; Tharayil, M.; Allyne, A.G. A survey of iterative learning. IEEE Control Syst. Mag. 2006, 26, 96–114. [Google Scholar]
  2. Ahn, H.S.; Chen, Y.; Moore, K.L. Iterative learning control: Brief survey and categorization. IEEE Trans. Syst. Man Cybern. Part C Rev. 2007, 37, 1099–1121. [Google Scholar] [CrossRef]
  3. Xiong, W.; Xu, L.; Ho, D.W.C.; Cao, J.; Huang, T. Synchronous and asynchronous iterative learning strategies of T-S fuzzy systems with measurable and unmeasurable state information. IEEE Trans. Fuzzy Syst. 2018, 26, 3042–3053. [Google Scholar] [CrossRef]
  4. Bu, X.; Hou, Z.; Yu, Q.; Yang, Y. Quantized data driven iterative learning control for a class of nonlinear systems with sensor saturation. IEEE Trans. Syst. Man Cybern. Syst. 2020, 50, 5119–5129. [Google Scholar] [CrossRef]
  5. Chen, Q.; Shi, H.; Sun, M. Echo state network-based backstepping adaptive iterative learning control for strict-feedback systems: An error-tracking approach. IEEE Trans. Cybern. 2020, 50, 3009–3022. [Google Scholar] [CrossRef]
  6. Li, X.; Liu, Y.H.; Yu, H.Y. Iterative learning impedance control for rehabilitation robots driven by series elastic actuators. Automatica 2018, 90, 1–7. [Google Scholar] [CrossRef]
  7. Jian, Y.; Huang, D.; Liu, J.; Min, D. High-precision tracking of piezoelectric actuator using iterative learning control and direct inverse compensation of hysteresis. IEEE Trans. Ind. Electron. 2019, 66, 368–377. [Google Scholar] [CrossRef]
  8. Cobb, M.K.; Barton, K.; Fathy, H.; Vermillion, C. Iterative learning-based path optimization for repetitive path planning, with application to 3-D crosswind flight of airborne wind energy systems. IEEE Trans. Control Syst. Technol. 2020, 28, 1447–1459. [Google Scholar] [CrossRef]
  9. Huang, D.; Min, D.; Jian, Y.; Li, Y. Current-cycle iterative learning control for high-precision position tracking of piezoelectric actuator system via active disturbance rejection control for hysteresis compensation. IEEE Trans. Ind. Electron. 2020, 67, 8680–8690. [Google Scholar] [CrossRef]
  10. Liu, G.; Hou, Z. Adaptive iterative learning control for subway trains using multiple-point-mass dynamic model under speed constraint. IEEE Trans. Intell. Transp. Syst. 2021, 22, 1388–1400. [Google Scholar] [CrossRef]
  11. Wang, L.; Yu, J.; Zhang, R.; Li, P.; Gao, F. Iterative learning control for multiphase batch processes with asynchronous switching. IEEE Trans. Syst. Man Cybern. Syst. 2021, 51, 2536–2549. [Google Scholar] [CrossRef]
  12. Sun, M.; Wang, D. Iterative learning control with initial rectifying action. Automatica 2002, 38, 1177–1182. [Google Scholar] [CrossRef]
  13. Xu, J.-X.; Yan, R. On initial conditions in iterative learning control. IEEE Trans. Autom. Control 2005, 50, 1349–1354. [Google Scholar]
  14. Xu, J.-X.; Xu, J. On iterative learning from different tracking tasks in the presence of time-varying uncertainties. IEEE Trans. Syst. Man Cybern. Part B Cybern. 2004, 34, 589–597. [Google Scholar] [CrossRef] [PubMed]
  15. Chi, R.; Hou, Z.; Xu, J.-X. Adaptive ILC for a class of discrete-time systems with iteration-varying trajectory and random initial condition. Automatica 2008, 44, 2207–2213. [Google Scholar] [CrossRef]
  16. Yu, M.; Huang, D.; He, W. Robust adaptive iterative learning control for discrete-time nonlinear systems with both parametric and nonparametric uncertainties. Int. J. Adapt. Control Signal Process. 2016, 30, 972–985. [Google Scholar] [CrossRef]
  17. Chien, C.-J. A combined adaptive law for fuzzy iterative learning control of nonlinear systems with varying control tasks. IEEE Trans. Fuzzy Syst. 2008, 16, 40–51. [Google Scholar] [CrossRef]
  18. Chi, R.; Hou, Z.; Jin, S.; Wang, D. Discrete-time adaptive ILC for non-parametric uncertain nonlinear systems with iteration-varying trajectory and random initial condition. Asian J. Control 2013, 15, 562–570. [Google Scholar] [CrossRef]
  19. Yin, C.; Xu, J.-X.; Hou, Z. A high-order internal model based iterative learning control scheme for nonlinear systems with time-iteration-varying parameters. IEEE Trans. Autom. Control 2010, 55, 2665–2670. [Google Scholar] [CrossRef]
  20. Yu, M.; Wang, J.; Qi, D. Discrete-time adaptive iterative learning control for high-order nonlinear systems with unknown control directions. Int. J. Control. 2013, 86, 299–308. [Google Scholar] [CrossRef]
  21. Zhao, D.; Yang, Y. An iterative learning control design method for nonlinear discrete-time systems with unknown iteration-varying parameters and control direction. Math. Probl. Eng. 2016, 7, 8971407. [Google Scholar] [CrossRef] [Green Version]
  22. Yin, C.; Xu, J.-X.; Hou, Z. An ILC scheme for a class of nonlinear continuous-time systems with time-iteration-varying parameters subjected to second-order internal model. Asian J. Control 2011, 13, 126–135. [Google Scholar] [CrossRef]
  23. Meng, D.; Moore, K.L. Robust iterative learning control for nonrepetitive uncertain systems. IEEE Trans. Autom. Control 2017, 62, 907–913. [Google Scholar] [CrossRef]
  24. Yu, M.; Li, C. Robust adaptive iterative learning control for discrete-time nonlinear systems with time-iteration-varying parameters. IEEE Trans. Syst. Man Cybern. Syst. 2017, 47, 1737–1745. [Google Scholar] [CrossRef]
  25. Yan, W.; Sun, M. Adaptive iterative learning control of discrete-time varying systems with unknown control direction. Int. J. Adapt. Control Signal Process. 2013, 27, 340–348. [Google Scholar] [CrossRef]
  26. Xu, Q.-Y.; Li, X.-D. Adaptive fuzzy ILC of nonlinear discrete-time systems with unknown dead zones and control directions. Int. J. Syst. Sci. 2018, 49, 1878–1894. [Google Scholar] [CrossRef]
  27. Wang, Y.C.; Chien, C.J. Design and analysis of fuzzy-neural discrete adaptive iterative learning control for nonlinear plant. Int. J. Fuzzy Syst. 2013, 15, 149–158. [Google Scholar]
Figure 1. (a) max k { 1 , , 50 } | e ϕ j ( k ) | (∗) versus iteration j. (b) e 100 ( k ) (solid line) and ϕ 100 ( k ) , ϕ 100 ( k ) (dotted lines) versus time k. (c) y 100 ( k ) (solid line) and y d 100 ( k ) (dotted line) versus discrete-time k. (d) u 100 ( k ) versus discrete-time k.
Figure 1. (a) max k { 1 , , 50 } | e ϕ j ( k ) | (∗) versus iteration j. (b) e 100 ( k ) (solid line) and ϕ 100 ( k ) , ϕ 100 ( k ) (dotted lines) versus time k. (c) y 100 ( k ) (solid line) and y d 100 ( k ) (dotted line) versus discrete-time k. (d) u 100 ( k ) versus discrete-time k.
Processes 10 01275 g001
Figure 2. (a) max k { 1 , , 50 } | e ϕ j ( k ) | (∗) versus iteration j. (b) e 100 ( k ) (solid line) and ϕ 100 ( k ) , ϕ 100 ( k ) (dotted lines) versus discrete-time k. (c) y 100 ( k ) (solid line) and y d 100 ( k ) (dotted line) versus discrete-time k. (d) u 100 ( k ) versus discrete-time k.
Figure 2. (a) max k { 1 , , 50 } | e ϕ j ( k ) | (∗) versus iteration j. (b) e 100 ( k ) (solid line) and ϕ 100 ( k ) , ϕ 100 ( k ) (dotted lines) versus discrete-time k. (c) y 100 ( k ) (solid line) and y d 100 ( k ) (dotted line) versus discrete-time k. (d) u 100 ( k ) versus discrete-time k.
Processes 10 01275 g002
Figure 3. (a) max k { 1 , , 50 } | e ϕ j ( k ) | (∗) versus iteration j. (b) e 100 ( k ) (solid line) and ϕ 100 ( k ) , ϕ 100 ( k ) (dotted lines) versus discrete-time k. (c) y 100 ( k ) (solid line) and y d 100 ( k ) (dotted line) versus discrete-time k. (d) u 100 ( k ) versus discrete-time k.
Figure 3. (a) max k { 1 , , 50 } | e ϕ j ( k ) | (∗) versus iteration j. (b) e 100 ( k ) (solid line) and ϕ 100 ( k ) , ϕ 100 ( k ) (dotted lines) versus discrete-time k. (c) y 100 ( k ) (solid line) and y d 100 ( k ) (dotted line) versus discrete-time k. (d) u 100 ( k ) versus discrete-time k.
Processes 10 01275 g003
Figure 4. (a) max k { 1 , , 200 } | e ϕ j ( k ) | (∗) versus iteration j. (b) e 100 ( k ) (solid line) and ϕ 100 ( k ) , ϕ 100 ( k ) (dotted lines) versus discrete-time k. (c) y 100 ( k ) (solid line) and y d 100 ( k ) (dotted line) versus discrete-time k. (d) u 100 ( k ) versus discrete-time k.
Figure 4. (a) max k { 1 , , 200 } | e ϕ j ( k ) | (∗) versus iteration j. (b) e 100 ( k ) (solid line) and ϕ 100 ( k ) , ϕ 100 ( k ) (dotted lines) versus discrete-time k. (c) y 100 ( k ) (solid line) and y d 100 ( k ) (dotted line) versus discrete-time k. (d) u 100 ( k ) versus discrete-time k.
Processes 10 01275 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chien, C.-J.; Wang, Y.-C. A Fuzzy System Based Iterative Learning Control for Nonlinear Discrete-Time Systems with Iteration-Varying Uncertainties. Processes 2022, 10, 1275. https://doi.org/10.3390/pr10071275

AMA Style

Chien C-J, Wang Y-C. A Fuzzy System Based Iterative Learning Control for Nonlinear Discrete-Time Systems with Iteration-Varying Uncertainties. Processes. 2022; 10(7):1275. https://doi.org/10.3390/pr10071275

Chicago/Turabian Style

Chien, Chiang-Ju, and Ying-Chung Wang. 2022. "A Fuzzy System Based Iterative Learning Control for Nonlinear Discrete-Time Systems with Iteration-Varying Uncertainties" Processes 10, no. 7: 1275. https://doi.org/10.3390/pr10071275

APA Style

Chien, C. -J., & Wang, Y. -C. (2022). A Fuzzy System Based Iterative Learning Control for Nonlinear Discrete-Time Systems with Iteration-Varying Uncertainties. Processes, 10(7), 1275. https://doi.org/10.3390/pr10071275

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