images
Computer Modeling in Engineering & Sciences
images

DOI: 10.32604/cmes.2022.020569

ARTICLE

Improved High Order Model-Free Adaptive Iterative Learning Control with Disturbance Compensation and Enhanced Convergence

Zhiguo Wang*, Fangqing Gao and Fei Liu

Key Laboratory of Advanced Process Control for Light Industry (Ministry of Education), Institute of Automation, Jiangnan University, Wuxi, 214122, China
*Corresponding Author: Zhiguo Wang. Email: zhiguowang@jiangnan.edu.cn
Received: 01 December 2021; Accepted: 22 February 2022

Abstract: In this paper, an improved high-order model-free adaptive iterative control (IHOMFAILC) method for a class of nonlinear discrete-time systems is proposed based on the compact format dynamic linearization method. This method adds the differential of tracking error in the criteria function to compensate for the effect of the random disturbance. Meanwhile, a high-order estimation algorithm is used to estimate the value of pseudo partial derivative (PPD), that is, the current value of PPD is updated by that of previous iterations. Thus the rapid convergence of the maximum tracking error is not limited by the initial value of PPD. The convergence of the maximum tracking error is deduced in detail. This method can track the desired output with enhanced convergence and improved tracking performance. Two examples are used to verify the convergence and effectiveness of the proposed method.

Keywords: Pseudo partial derivative; enhanced convergence; tracking error; disturbance compensation

1  Introduction

Generally, there are two strategies for controller design. One is based on the exact process information, and the other is based on the input-output data which is a model-free controller design strategy. Many practical industrial processes contain nonlinearity, uncertainty, and time-varying characteristics. In this case, model-based control becomes inapplicable because it is difficult to obtain an accurate mathematical model of the controlled object. In recent years, the model-free methods have attracted the attention of many researchers because of their remarkable advantages and have achieved some valuable results [1,2].

Model-free adaptive control (MFAC) is one of the typical data-driven control methods proposed by Hou [3]. The principle of the MFAC algorithm is to continuously generate a dynamic linearization model containing pseudo partial derivative (PPD), which transforms the nonlinear system into a time-varying linearized model that utilizes the input and output data of every step [4]. Compared with the model-based algorithm, model-free adaptive control has the characteristics of fewer identification parameters and fewer calculations. This method has been applied in various fields, for example, process control [5], motor system [6], and so on. Arimoto et al. [7] proposed iterative learning control (ILC), which is suitable for a class of controlled objects with repetitive running characteristics, such as robots with reciprocating operation characteristics [8] and a chemical reactor in a batch process [9]. Iterative learning control was applied to improve the tracking performance of nonlinear processes whose model and parameters are unknown [10]. Zhao et al. [11] used iterative learning to design the unbiased finite impulse response. Combining the advantages of model-free adaptive control with iterative learning control, model-free adaptive iterative learning control (MFAILC) was proposed in [12]. Unlike the iterative learning control with fixed learning gain, the learning gain of model-free adaptive iterative learning control is time-varying. It transforms the original nonlinear time-domain into a nonlinear iterative domain and can be applied to other nonlinear controlled objects whose structures and parameters are time-varying [13,14]. Thus MFAILC algorithm has better stability and adaptability for complex systems [15].

Based on the previous research results, many researchers have further investigated the model-free adaptive iterative learning control algorithm. In [16], an improved model-free adaptive iterative learning control method which based on an encoding and decoding quantization mechanism was proposed for a class of unknown nonlinear systems with data quantization. Li [17] adopted a neural network to optimize the control law and the penalty of PPD. For a class of nonlinear discrete systems with load disturbance and data loss, Hua et al. [18] proposed a control method with data compensation. Model-free adaptive iterative learning control was applied to a multi-agent system by Bu et al. [19]. Chi et al. [20] proposed a high order model-free adaptive iterative learning control, which uses the input and output information of the previous iterations to construct a new control input. It is obvious from the above observations that MFAILC has made great progress in theoretical research. In the control process, convergence is a critical index for ILC so as to MFAILC [21]. However, the convergence speed is related to the initial pseudo partial derivative, and it often takes several experiments to get the appropriate value of the initial PPD [22]. The slower convergence will deteriorate the performance of the controlled process [23].

To obtain a faster convergence, the high-order algorithm was used to estimate the value of pseudo partial derivatives, which introduced more parameters to be adjusted unexpectedly [24]. This method brings more difficulties in tuning parameters. On the basis of the former, this paper proposed an improved high-order model-free adaptive iterative control (IHOMFAILC). The main contributions of this paper are as follows:

• An improved PPD estimation algorithm, which takes the average value of the previous PPDs can achieve faster convergence without introducing more parameters.

• A differential module is introduced into the criteria function of the controller to compensate for the dynamic performance of the system under random disturbance.

• The proposed method contains fewer parameters to design and with the number of iterations increases, the maximum error gradually converges to zero.

The rest of this paper is organized as follows: Section 2 introduces the basic principle of model-free adaptive iterative learning control. The proposed method with disturbance compensation and enhanced convergence is shown in detail in Section 3. In Section 4, the convergence and stability of the proposed algorithm are proved in detail. Two examples are used to illustrate the effectiveness of the method in Section 5 and the conclusions are drawn in Section 6.

2  Basic Principle of MFAILC

Consider the following nonlinear discrete-time SISO system:

yk(t+1)=f(yk(t),yk(t−1),⋯,yk(t−ny),uk(t),uk(t−1),⋯,uk(t−nu))(1)

where uk(t) and yk(t) are the kth iteration input and output of system, respectively, ny and nu are the unknown orders of system, and f(⋅) is a nonlinear function.

The following two assumptions are given for system (1) to make the discussion rigorous:

Assumption 1: The partial derivative of the nonlinear function f(⋅) with respect to its control input uk(t) is continuous.

Assumption 2: System (1) satisfies the generalized Lipschitz condition, that is, for any time t and iteration k, if |Δuk(t)|≠0 , then the system (1) satisfies the following conditions:

|Δyk(t+1)|≤b|Δuk(t)|(2)

where Δyk(t+1)=yk(t+1)−yk−1(t+1),Δuk(t)=uk(t)−uk−1(t), and b is a positive constant. Based on the two assumptions, Lemma 1 can be achieved:

Lemma 1: For any nonlinear system that satisfies assumptions 1 and 2, if |Δuk(t)|≠0, there must exist time-varying pseudo-partial derivative ϕk(t), such that system (1) can be transformed into the following generic model:

Δyk(t+1)=ϕk(t)Δuk(t)(3)

where ϕk(t) is bounded and satisfies |ϕk(t)|≤b.

In order to design the control law, the following criteria function is considered:

J(uk(t))=|yr(t+1)−yk−1(t+1)|2+λ|uk(t)−uk−1(t)|2(4)

where λ is a weighting factor used to restrict the control input variation between adjacent iterations, and yr(t+1) is the desired output.

Substituting (3) into (4) and minimizing (4) with respect to uk(t), we get the following control law:

uk(t)=uk−1(t)+ρkϕk(t)ek−1(t+1)|ϕk(t)|2+λ(5)

where ρk is a step factor, and ek−1(t+1)=yr(t+1)−yk−1(t+1) is the output tracking error. Since ϕk(t) is a time-varying parameter which is difficult to acquire its accurate values, it needs to be estimated based on input and output data.

As with solving uk(t), the following criteria function is denoted:

J(ϕk(t))=|Δyk−1(t+1)−ϕ^k(t)Δuk−1(t)|2+μ|ϕ^k(t)−ϕ^k−1(t)|2(6)

where μ>0 is a weighting factor, ϕ^k(t) is the estimation of ϕk(t).

Taking the derivative of Eq. (4) with respect to ϕ^k(t) results in the following parameter estimation algorithm:

ϕ^k(t)=ϕ^k−1(t)+ηkΔuk−1(t)μ+|Δuk−1(t)|2(Δyk−1(t+1)−ϕ^k−1(t)Δuk−1(t))(7)

where ηk>0 is a step factor.

In order to make the algorithm estimate value of Eq. (7) more accurate, the following reset algorithm is denoted: if |ϕ^k(t)|<ε or |Δuk(t)|≤ε or sign(ϕ^k(t))≠sign(ϕ^0(t))

ϕ^k(t)=ϕ^0(t)(8)

where ε is a small positive constant, and ϕ^0(t) is the initial value of iteration. The reset algorithm ensures the condition of Lemma 1 and the stronger tracking ability of Eq. (7) for time-varying parameters.

However, the rapidity of convergence of the maximum error of MFAILC is related to the value of the initial PPD. The larger the value is, the slower the convergence rate is and the larger the error is. Therefore, an improved method is proposed to solve the above problem.

3  Improved High-Order Model-Free Adaptive Iterative Learning Control Method

In the proposed method, a differential link is added to compensate for the impact of random disturbance on the system and only use input and output data to design the control law.

Consider the following criteria function:

J(uk(t))=|yr(t+1)−yk(t+1)|2+λ|uk(t)−uk−1(t)|2+γ|ek−1(t+1)−ek−2(t+1)|2(9)

where γ≥0 is weighting factor limits the change of tracking error, and when γ=0 the criteria function is the same as Eq. (4). Substituting (3) into (9), we can get:

J(uk(t))=|ek−1(t+1)−ϕk(t)Δuk(t)|2+λ|uk(t)−uk−1(t)|2+γ|ek−1(t+1)−ϕk(t)Δuk(t)−[ek−2(t+1)−ϕk−1(t)Δuk−1(t)]|2(10)

To get the optional solution, we minimizing (10) with respect to uk(t) leads to following control law:

uk(t)=uk−1(t)+ρkϕk(t)λ+|ϕk(t)|2+γ|ϕk(t)|2(yr(t+1)−yk−1(t+1))+γϕk(t)(ek−1(t+1)−ek−2(t+1)+ϕk−1(t)Δuk(t))λ+|ϕk(t)|2+γ|ϕk(t)|2(11)

according to (11), ϕk(t) is an unknown time-varying parameter that needs to be estimated utilizing input and output data.

Unlike other high-order algorithms designed based on control input [25] and tracking error terms [26], the improved high-order algorithm uses previous iterations to estimate the PPD. Compared to the algorithm in [24], the improved algorithm does not need to take different weight coefficients for the previous PPDs.

If the number of iterations is less than the iteration learning n, the estimation algorithm of Eq. (7) is adopted:

ϕ^k(t)=ϕ^k−1(t)+ηkΔuk−1(t)μ+|Δuk−1(t)|2(Δyk−1(t+1)−ϕ^k−1(t)Δuk−1(t))(12)

when the number of iterations is larger than the iteration learning n, the following algorithm is adopted to reduce the dependence of the error convergence speed on the initial PPD value. Denote the following criteria function:

J(ϕ^k(t))=|Δyk−1(t+1)−ϕ^k(t)Δuk−1(t)|2+μ|ϕ^k(t)−1n∑i=1nϕ^k−i(t)|2(13)

Taking the partial derivative of ϕ^k(t) in Eq. (13):

∂J(ϕ^k(t))∂ϕ^k(t)=2(−ϕ^k(t))(Δyk−1(t+1)−ϕ^k(t)Δuk−1(t))+2μ(ϕ^k(t)−1n∑i=1nϕ^k−i(t))(14)

Minimizing Eq. (14) by solving ∂(ϕ^k(t))/∂ϕ^k(t)=0, the estimate algorithm is expressed:

ϕ^k(t)=Δyk−1(t+1)Δuk−1(t)μ+|Δuk−1(t)|2+μηkμ+|Δuk−1(t)|21n∑i=1nϕ^k−i(t)(15)

where n is number of iterative learning. The reset algorithm of ϕ^k(t) is the same as Eq. (8).The improved algorithm uses the PPDs of previous n iterations to estimate the current PPD. Thus the PPD can learn more from past data, making the tracking error has fast convergence. Therefore, the overall control scheme of IHOMFAILC consists of (16)–(18):

{ϕ^k(t)=ϕ^k−1(t)+ηkΔuk−1(t)μ+|Δuk−1(t)|2(Δyk−1(t+1)−ϕ^k−1(t)Δuk−1(t))2≤k<nϕ^k(t)=Δyk−1(t+1)Δuk−1(t)μ+|Δuk−1(t)|2+μηkμ+|Δuk−1(t)|21n∑i=1nϕ^k−i(t)k≥n(16)

uk(t)=uk−1(t)+ρkϕk(t)λ+|ϕk(t)|2+γ|ϕk(t)|2(yr(t+1)−yk−1(t+1))+γϕk(t)(ek−1(t+1)−ek−2(t+1)+ϕk−1(t)Δuk(t))λ+|ϕk(t)|2+γ|ϕk(t)|2(17)

ϕ^k(t)=ϕ^0(t),if|ϕ^k(t)|<εor|Δuk(t)|≤ε(18)

4  Convergence Analysis

In this section, the convergence of the proposed method is mainly proved. In order to make the following discussion more rigorous, the following assumption is given:

Assumption 3: For all k∈{0,1,2,…,N−1} and t∈{0,1,2,…,T−1}, if |Δuk(t)|≠0,ϕk(t)>0.

Assumption 3 is similar to a restriction on the direction of control input gain.

Lemma 2: If system (1) satisfies assumption 1–3, there must exist λmin>0 so that when λ>λmin, the IHOMFAILC method has the following properties:

1.    For all k∈{0,1,2,…,N−1} and t∈{0,1,2,…,T−1}, ϕ^k(t) is bounded.

2.    When k→∞, the system output tracking error converges to 0.

3.    The close-loop system is BIBO stability, for all k∈{0,1,2,…,N−1} andt∈{0,1,2,…,T−1}, uk(t) is bounded.

Proof: The first step is to prove that ϕ^k(t) is bounded.

when |ϕ^k(t)|≤ε or |uk(t)|≤ε,ϕ^k(t)=ϕ^0(t), ϕ^k(t) is bounded.

In other cases, when 2≤k<n, denote the following estimation error of the PPD parameters:

ϕ⌣k(t)=ϕ^k(t)−ϕk(t)(19)

Subtracting ϕk(t) from both sides of Eq. (12):

ϕ⌣k(t)=ϕ⌣k−1(t)−(ϕk(t)−ϕk−1(t))+ηΔuk−1(t)μ+|Δuk−1(t)|×(Δyk−1(t)−ϕ^k−1(t)Δuk−1(t))(20)

Denote Δϕk(t)=ϕk(t)−ϕk−1(t), substituting (3) into (20), we have:

ϕ⌣k(t)={1−η|Δuk−1(t)|2μ+|Δuk−1(t)|2}ϕ⌣k−1(t)−Δϕk(t)(21)

When 0<η≤1, μ>0, from Eq. (21) we can get:

0<|1−η|Δuk−1(t)|2μ+|Δuk−1(t)|2|<1(22)

From Lemma 1 we get |ϕk(t)|≤b, so |ϕk(t)−ϕk−1(t)|≤2b. Taking absolute worthwhile at both sides of Eq. (21):

ϕ⌣k(t)≤d|ϕ⌣k−1(t)|+2b≤dn|ϕ⌣0(t)|+(1−dn−1)1−d2b(23)

where d={1−η|Δuk−1(t)|2μ+|Δuk−1(t)|2}, therefore ϕ⌣k(t) is bounded.

when n≤k and |uk(t)|>ε, we have:

ϕ^k(t)=Δyk−1(t+1)Δuk−1(t)μ+|Δuk−1(t)|2+μηkμ+|Δuk−1(t)|21n∑i=1nϕk−i(t)(24)

Substituting Eq. (3) into (24) and taking absolute worthwhile at both sides, we have:

|ϕ^k(t)|=|Δϕk−1(t)Δu2k−1(t)μ+|Δuk−1(t)|2+μηkμ+|Δuk−1(t)|21n∑i=1nϕ^k−i(t)|<|ϕk−1(t)|+|ηkn|∑i=1n|ϕ^k−i(t)|(25)

Eq. (25) can finally be simplified as follows:

|ϕ^k(t)|<|ϕk−n+1(t)|+|ηkn|∑j=1k−n+1∑i=1n|ϕ^k−i(t)|(26)

From Lemma 1, we know ϕk(t) is bounded and therefore |ϕk−n(t)| is bounded. Finally from Eq. (26), we can prove that ϕ^k(t) is bounded.

Then the convergence of tracking error and the bounded control input are proved as follows:

Denote if k≤1, then |Δuk(t)|=0. Taking absolute worthwhile at both sides of Eq. (11):

|Δuk(t)|=|ρkϕ^k(t)ek−1(t+1)+γϕ^k(t)(ϕ^k−1(t)−ϕk−1(t))Δuk−1(t)λ+|ϕ^k(t)|2+γ|ϕ^k(t)|2|(27)

Since ϕk(t) and ϕ^k(t) are bounded, when λ>λmin>0, then there must exist M1,M2,M3 satisfy the following inequality:

|ϕ^k(t)λ+|ϕ^k(t)|2+γ|ϕ^k(t)|2|≤|ϕ^k(t)2λ(1+γ)|ϕ^k(t)||<|12λmin(1+γ)||=M1≤0.5b(28)

0<M2<|γϕ^k(t)ϕk(t)λ+|ϕ^k(t)|2+γ|ϕ^k(t)|2|≤b|γϕ^k(t)2λ(1+γ)|ϕ^k(t)||<|γb2λmin(1+γ)||<0.5(29)

From Eqs. (28) and (29) we get:

M1ϕk(t)≤M3<0.5(30)

M2+M3<1(31)

Since ϕk(t) and ϕ^k(t) are bounded, there exists d1 that:

|γϕ^k(t)(ϕ^k−1(t)−ϕk−1(t))λ+|ϕ^k(t)|2+γ|ϕ^k(t)|2|≤d1<|γϕ^k(t)bλ+|ϕ^k(t)|2+γ|ϕ^k(t)|2|≤M3(32)

According to Eqs. (28)–(31), we have:

Δuk(t)<ρM1|ek−1(t+1)|+d1|Δuk−1(t)|<ρM1∑i=1k−1d1k−1−i|ei(t+1)|(33)

From Eq. (33) we can get, if ek(t) converges to 0, then uk(t) is bounded.

|ek(t+1)|=|yr(t+1)−yk(t+1)|=|yr(t+1)−yk−1(t+1)−ϕk(t)Δuk(t)|<|1−ρkϕk(t)ϕ^k(t)λ+|ϕ^k(t)|2+γ|ϕ^k(t)|2|×|ek−1(t+1)|+ϕk(t)|γϕ^k(t)(ϕ^k−1(t)−ϕk−1(t))λ+|ϕ^k(t)|2+γ|ϕ^k(t)|2||Δuk−1(t)|(34)

Let d2=1−ρM2, substituting Eq. (33) into (34) we have:

|ek(t+1)|<d2|ek−1(t+1)|+d1ϕk(t)|Δuk−1(t)|(35)

Then Eq. (35) can be described as follows:

|ek(t+1)|<d2k−1|e1(t+1)|+d1∑i=1k−2d2k−2−iϕi+1(t)ρM1∑J=1i−1d1i−1−j|ej(t+1)|(36)

Let d3=ρM3, according to Eq. (29) we get M1ϕk(t)≤M3<0.5, so M1ϕi+1(t)≤M3<0.5. Eq. (36) can be transformed into:

|ek(t+1)|<d2k−1|e1(t+1)|+d1d3∑i=1k−2d2k−2−i∑j=1i−1d1i−1−j|ej(t+1)|(37)

Let gk(t+1)=d2k−1|e1(t+1)|+d1d3∑i=1k−2d2k−2−i∑j=1i−1d1i−1−j|ej(t+1)|, from Eq. (37) we have:

|ek(t+1)|<gk(t+1)(38)

where g2(t+1)=d2|e1(t+1)|. According to (38), when k→∞, if gk(t+1)converges to 0 then ek(t+1) will converge to 0. Therefore, it is possible to prove the convergence of gk(t+1) rather than ek(t+1) directly.

gk+1(t+1)=d2k|e1(t+1)|+d1d3∑i=1k−1d2k−1−i∑j=1i−1d1i−1−j|ej(t+1)|=d2gk(t+1)+d3d1k−1|e1(t+1)|+⋯+d3d12|ek−1(t+1)|+d3d1|ek−1(t+1)|<d2gk(t+1)+d3d1k−1|e1(t+1)|+⋯+d3d12|ek−1(t+1)|+d3d1|gk−1(t+1)|=d2gk(t+1)+hk−1(t+1)(39)

where hk−1(t+1)=d3d1k−1|e1(t+1)|+⋯+d3d12|ek−1(t+1)|+d3d1|gk−1(t+1)|.

According to Eq. (31): M2+M3<1, so we have:

d2=1−ρM2>ρ(M2+M3)−ρM2=ρM3=d3(40)

Substituting Eq. (40) into hk−1(t+1), we get the following inequality:

hk−1(t+1)<d3d1k−1|e1(t+1)|+⋯+d3d12|ek−1(t+1)|+d3d1|gk−1(t+1)|<d3d1k−1|e1(t+1)|+⋯d3d12|ek−1(t+1)|+d2d1(d2k−2|e1(t+1)|+d1d3∑i=1k−3d2k−3−i∑j=1i−1d1i−1−j|ej(t+1)|)(41)

Eq. (41) hk−1(t+1) can finally be expressed as:

hk−1(t+1)=d1(d2k−1|e1(t+1)|+d1d3∑i=1k−2d2k−2−i∑j=1i−1d1i−1−j|ej(t+1)|)=d1gk(t+1)(42)

Substituting (42) into (39) yields:

gk+1(t+1)<(d2+d1)gk(t+1)(43)

When d2+d1<1, gk(t+1) converges gradually. Selecting ρ∈(0,1),γ>0 to make:

1−|ρkϕk(t)ϕ^k(t)λ+|ϕ^k(t)|2+γ|ϕ^k(t)|2|+|γϕ^k(t)(ϕ^k−1(t)−ϕk−1(t))λ+|ϕ^k(t)|2+γ|ϕ^k(t)|2|<1(44)

Thus the following inequality is given:

γ<|ρϕk(t)ϕ^k−1(t)−ϕk−1(t)|(45)

Substituting d2+d1<1 into Eq. (43),

limk→∞gk+1(t+1)<limk→∞(d2+d1)gk(t+1)<limk→∞(d2+d1)k−1g2(t+1)=0(46)

According to inequality (38) and (46), with the increase of the number of iterations k, gk(t+1) converge to 0 gradually, and therefore ek(t+1) will converge to 0 finally.

From inequality (33), we have:

|uk(t)|≤∑i=1k|Δuk(t)|<ρM1∑i=1k∑j=1i−1d1i−1−j|ej(t+1)|<ρM11−d1(e1(t+1)+g2(t+1)+⋯+gk−1(t+1))(47)

Eq. (47) can be reorganized as follows:

|uk(t)|<ρM11−d1(e1(t+1)+g2(t+1)1−d1−d2)(48)

This implies that uk(t) is bounded. Lemma 2 (1)–(3) are established.

5  Simulation

Example 1: Consider the following time-varying system [14]:

y(t+1)={y(t)1+y(t)2+u(t)30≤t≤50(y(t)y(t−1)y(t−2)u(t−1))×(y(t−2)−1)+α(t)u(t)/(1+y(t−1)2+y(t−2)2)50<t≤100(49)

The desired output is:

yr(t+1)={0.5×(−1)round(t/10)0≤t≤300.5×sin⁡(tπ/10)+0.3×cos⁡(tπ/10)30<t≤700.5×(−1)round(t/10)70<t≤100(50)

where α(t)=1+round(t/50) is a time-varying parameter of the system (49). For IHOMAFILC, the controller parameters are set as: λ=0.4, ρ=1, μ=1, η=0.5, γ=0.01 and ε=10−5. The number of previous iterative learning is set as n=3.The initial value of PPD is set to ϕ0(t)=10.

Table 1 shows the ITAE performance index of the two methods. When the number of iterations is set to 10, 20, 30, and 40 while other parameters are unchanged, the IATE values of IHOMFAILC is less than that of MFAILC. As the number of iterations increase, the IATE values of both methods shows a downward trend. Therefore the proposed method performance index is better.

Fig. 1 shows the tracking results of IHOMFAILC and MFAILC with different iterations. From the 15th iteration result, the proposed method can gradually approach desired output ultimately tracks it at 40 iterations. The output of the MFAILC method at the 40th iteration still has a large error with the expected output.

Figure 1: Tracking results of IHOMFAILC and MFAILC with different iterations

The maximum tracking error of IHOMFAILC and MFAILC with different initial PPDs is shown in Fig. 2. As the picture shows, the change of the initial value of PPD has a great impact on the error convergence speed of MFAILC. The proposed method ensures fast convergence and the final error convergence to 0. Results demonstrate that IHOMFAILC can achieve fast convergence with a few attempts to select the best initial value of PPD.

Figure 2: Comparison of max tracking error of MFAILC and IHOMFAILC with different initial PPDs

Example 2: In order to verify the effectiveness of the proposed algorithm, simulation is performed on a permanent magnet DC linear motor [1]. The dynamic characteristics of the motor are described as follows:

x˙(t)=υ˙(t)υ˙(t)=u(t)−ffriction(t)−fripple(t)m(51)

where x(t) is the position of linear motor (m), υ(t) is the speed (m/s), ffriction(t) is the friction force (N), fripple(t) is the ripple force (N), u(t) is the developed force (N) , t is the continuous-time (s).

The mathematical model between friction force and ripple force is described as follows:

ffriction(t)=(fc+(fs−fc)e−(x˙/x˙δ)δ+fυx˙)sgnx˙fripple(t)=b1sin⁡(ω0x(t))(52)

where fc is the minimum value of coulomb friction and fs is the static friction, xδ is the lubrication parameters and fυ is the load parameters,δ is an additional empirical parameter. m=0.59, xδ=0.1, δ=1, fc=10, fs=20, fυ=10, b1=8.5, ω0=314s−1.

Denote x1(t)=x(t), x2(t)=υ(t). Eq. (45) can be described as follows:

[x˙1(t)x˙2(t)]=[x2(t)−ffriction(t)+fripple(t)m]+[01m]×u(t)y(t)=x2(t)+d(t)(53)

where d(t) is a zero-mean white noise with a variance 0.005. The desired output is yr(t)=10sin⁡(πt). Discretize (46) by Euler Formula and take sample time as ts=0.001. The controller parameters are set asλ=0.06, ρ=1, μ=0.3, η=1, and γ=0.005. Fig. 3 shows the maximum tracking error convergence of IHOMFAILC and MFAILC. It can be seen from the figure that the errors of the two methods are the same at the 11th and 14th iterations and the convergence errors of the IHOMFAILC in other iterations are all less than MFAILC. This simulation shows that the proposed method has fast error convergence even though the system exists interference.

Figure 3: Max tracking error convergence of MFAILC and IHOMFAILC

6  Conclusion

Based on the model-free adaptive iterative learning control, a novel IHOMFAILC method has been proposed in this paper based on compact form dynamic linearization. When the system is disturbed, the proposed method can still keep track of the desired output. Furthermore, this method can achieve enhanced convergence without many experiments to find a suitable PPD value. Two examples show that the proposed method has good trajectory tracking and disturbance resistance.

Acknowledgement: This work was supported by the Natural Science Foundation of China, under Grant No. 61773183.

Funding Statement: The authors received no specific funding for this study.

Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.

References

 1.  Madadi, E., Söffker, D. (2018). Model-free control of unknown nonlinear systems using an iterative learning concept: Theoretical development and experimental validation. Nonlinear Dynamics, 94(2), 1151–1163. DOI 10.1007/s11071-018-4415-7.

 2.  Liu, D., Yang, G. H. (2018). Model-free adaptive control design for nonlinear discrete-time processes with reinforcement learning techniques. International Journal of Systems Science, 49(11), 2298–2308. DOI 10.1080/00207721.2018.1498557.

 3.  Hou, Z. S. (1994). The parameter identification, adaptive control and model-free learning adaptive control for nonlinear systems (Ph.D. Thesis). North-eastern University

 4.  Bu, X. H., Hou, Z. S., Zhang, H. W. (2018). Data-driven multi-agent systems consensus tracking using model-free adaptive control. IEEE Transactions on Neural Network and Learning Systems, 29(5), 1514–1524. DOI 10.1109/TNNLS.2017.2673020.

 5.  Hou, Z. S., Liu, S., Yin, C. (2016). Local learning-based model-free adaptive predictive control for adjustment of oxygen concentration in the syngas manufacturing industry. IET Control Theory & Applications, 10(12), 1384–1394. DOI 10.1049/iet-cta.2015.0835.

 6.  Di, S., Li, H. (2018). Model-free adaptive speed control on the travelling wave ultrasonic motor. Journal of Electronic Engineering, 69(1), 14–23. DOI 10.1515/jee-2018-0002.

 7.  Arimoto, S., Kawamura, S., Miyazaki, F. (1984). Bettering operation of robots by learning. Journal of Robotic Systems, 1(2), 123–140. DOI 10.1002/(ISSN)1097-4563.

 8.  Treesatayapun, C. (2015). Data input-output adaptive controller based on IF-THEN rules for a class of non-affine discrete-time systems: The robotic plant. Journal of Intelligent & Fuzzy Systems, 28(2), 661–668. DOI 10.3233/IFS-141347.

 9.  Xiong, Z., Zhang, J. (2003). Product quality trajectory tracking in batch processes using iterative learning control based on time-varying perturbation models. Industrial & Engineering Chemistry Research, 42(26), 6802–6814. DOI 10.1021/ie034006j.

10. Fu, Q., Li, X. D., Du, L. L., Xu, G. Z., Wu, J. R. (2018). Consensus control for multi-agent systems with quasi-one-sided Lipschitz nonlinear dynamics via iterative learning algorithm. Nonlinear Dynamics, 91(4), 2621–2630. DOI 10.1007/s11071-017-4035-7.

11. Zhao, S., Shmaliy, Y. S., Ahn, C. K., Liu, F. (2018). Adaptive-horizon iterative UFIR filtering algorithm with applications. IEEE Transactions on Industrial Electronics, 65(8), 6393–6402. DOI 10.1109/TIE.2017.2784405.

12. Jin, S. T., Hou, Z. S., Chi, R. H., Liu, X. B. (2012). Data-driven model-free adaptive iterative learning control for a class of discrete-time nonlinear systems. Control Theory & Applications, 29(8), 1001–1009. DOI 10.7641/j.issn.1000-8152.2012.8.lcta120438.

13. Bu, X. H., Wang, S., Hou, Z. S. (2019). Model-free adaptive iterative learning control for a class of nonlinear systems with randomly varying iteration lengths. Journal of the Franklin Institute, 356(5), 2491–2504. DOI 10.1016/j.jfranklin.2019.01.003.

14. Hou, Z. S., Jin, S. T. (2013). Model-free adaptive control: Theory and Application. Beijing: Science Press.

15. Chi, R. (2006). Adaptive iterative learning control for nonlinear discrete-time systems and its applications. Beijing, China: Beijing Jiao Tong University.

16. Zhu, P. P., Bu, X. H., Yu, Q. X., Hou, Z. S., Liang, J. Q. (2020). An improved model-free adaptive iterative learning control algorithm with data quantization. Control Theory & Applications, 37(5), 1178–1184. DOI 10.1002/rnc.5107.

17. Li, J. (2010). Model-free adaptive iterative learning control based on BP algorithm. Microcomputer & Its Applications, 29(22), 83–85. DOI 10.19358/j.issn.1674-7720.2010.22.028.

18. Hua, C. C., Qiu, Y. F., Guan, X. P. (2020). Enhanced model-free adaptive iterative learning control with load disturbance and data dropout. International Journal of Systems Science, 51(11), 2057–2067. DOI 10.1080/00207721.2020.1784492.

19. Bu, X. H., Yu, Q. X., Hou, Z. S. (2019). Model free adaptive iterative learning consensus tracking control for a class of nonlinear multi-agent systems. IEEE Transactions on Systems, 49(4), 677–686. DOI 10.1109/TSMC.2017.2734799.

20. Chi, R. H., Hou, Z. S., Yu, L., Sui, S. L. (2008). Higher order model-free adaptive iterative learning control. Control and Decision, 23(7), 795–798. DOI 10.3321/j.issn:1001-0920.2008.07.015.

21. Chi, R., Hou, Z. S., Jin, S. H., Huang, B. (2018). Computationally efficient data driven higher order optimal iterative learning control. IEEE Transactions on Neural Networks and Learning Systems, 29(12), 5971–5980. DOI 10.1109/TNNLS.2018.2814628.

22. Luo, B., Wu, H. N., Huang, T. (2018). Optimal output regulation for model-free Quanser helicopter with multistep Q-learning. IEEE Transactions on Industrial Electronics, 65(6), 4953–4961. DOI 10.1109/TIE.2017.2772162.

23. Zhao, S., Huang, B., Liu, F. (2017). Linear optimal unbiased filter for time-variant systems without apriori information on initial conditions. IEEE Transactions on Automatic Control, 62(2), 882–887. DOI 10.1109/TAC.2016.2557999.

24. Ai, Q. S., Ke, D., Zuo, J., Meng, W., Liu, Q. et al. (2020). High-order model-free adaptive iterative learning control of pneumatic artificial muscle with enhanced convergence. IEEE Transactions on Industrial Electronics, 67(11), 9548–9559. DOI 10.1109/TIE.2019.2952810.

25. Chi, R. H., Huang, B., Hou, Z., Jin, S. (2017). Data-driven high-order terminal iterative learning control with a faster convergence speed. Robust Nonlinear Control, 28(1), 103–119. DOI 10.1002/rnc.3861.

26. Chi, R. H., Liu, Y., Hou, Z. S., Jin, S. T. (2015). Data-driven Terminal Iterative learning control with high-order learning law for a class of non-linear discrete time multiple-input-multiple output systems. IET Control Theory & Applications, 9(7), 1075–1082. DOI 10.1049/iet-cta.2014.0754.

images
This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.