Experiment Research on Complex Optimization Algorithm-Based Adaptive Iterative Learning Control for Electro-Hydraulic Shaking Tables

: The adaptive iterative learning control method for electro-hydraulic shaking tables based on the complex optimization algorithm was proposed to overcome the potential stability problem of the traditional iteration control method. The system identiﬁcation precision’s inﬂuence on convergence was analyzed. Based on the real optimization theory and the mapping relationship between real vector space and complex vector space, the complex Broyden optimization iterative algorithm was proposed, and its stability and convergence was analyzed. To improve the stability and accelerate the convergence of the proposed algorithm, the complex steepest descent algorithm was proposed to cooperate with the complex Broyden optimization algorithm, which can adaptively optimize the complex steepest gradient iterative gain and update the system impedance in real time during the control process. The shaking tables experiment system was designed, applying xPC target rapid prototype control technology, and a series of experimental tests were performed. The results indicated that the proposed control method can quickly and stably converge to the optimal solution no matter whether the system identiﬁcation error is small or large, and, thus, veriﬁed that validity and feasibility of the proposed adaptive iterative learning method.


Introduction
The electro-hydraulic shaking table is the most used and reliable earthquake simulation test equipment in the field of structure anti-seismic research, especially in recent decades [1][2][3].However, for earthquake experiments with structures with a large span such as trains, bridges and pipelines, etc., a single shaking table is not at all sufficient; thus, a shaking table array system composed of multiple identical or different shaking tables is proposed to solve this problem [4][5][6][7].The shaking table array is more complex compared with the single shaking table system due to the tracking precision needed and also the synchronization precision needed to protect the test specimen [8].Usually, the large structure used to carry out earthquake test experiments is very complex and represents obvious nonlinearity, and the shaking table array system also encounters various disturbances, etc., all of which increase the difficulty of the control system and, thus, influence the experiment effect [9,10].
For the shaking table control problem, many experts have devoted much effort to improving the response performance from different perspectives.Guo [11] adopts the latest advances in deep learning to develop a physics-guided long short-term memory network for system identification of the shaking table.Huang [12] analyzed the coupling phenomena and reasons for the multiaxis hydraulic servo table based on a dynamic model of a multiaxis hydraulic servo shaking table.Li [13] analyzed the effects of the characteristics of the dual shaking tables-specimen interaction on the system performance.Soleymani [14] designed a fuzzy sliding mode supervisory controller for an electric seismic shaking table with variable payload.Zhao [15] proposed modal space three-state feedback and feedforward control to solve the coupling problem.Cheng [16] proposed a control algorithm based on the frequency-domain filtered-x least mean square adaptive algorithm to achieve power spectral density replication.
How to improve the control performance of shaking tables with system nonlinearity and parametric uncertainty has been of great interest.Yang [11] proposed a combined controller with ultra-local model principle-based adaptive SMC to improve the sweptsine control precision of an electro-hydraulic shaking table considering the nonlinearities, uncertainties and unexpected disturbance.In [12], an adaptive dynamic surface control method with nonlinear disturbance observers was proposed to improve the tracking precision of an electro-hydraulic shaking table with unknown time-varying disturbances.To obtain a good acceleration replication effect for a hydraulic shaking table, a feed-forward inverse model controller based on the RELS algorithm was proposed by Shen [13].In [14], a novel stochastic gradient algorithm was presented to estimate the harmonic information of a hydraulic shaking table.Giacomo [15] compared the results of two tests of automotive parts, respectively using a single-axis shaker and a three-axis shaker, which indicated that the MIMO control scheme had a better control effect.In [16], an acceleration waveform replication scheme combining an offline feed-forward compensator with an online adaptive controller was proposed to improve the control performance of a hydraulic shaking table.
The offline iterative control method was proposed for the electro-hydraulic shaking tables system, and has been become the most widely and successfully applied compensation control method.Offline iterative control is an incausal feedforward control technique based on the inverse model of the system linear frequency model.Offline iterative control can realize the non-minimum phase system to obtain very high precision, which is unrealizable using the real-time feedforward control method.In 1976, Cryer first proposed the offline iterative control method and applied this technique to the road's waveform replication of a four-channel road simulator [17].In 1977, MTS introduced the remote parameter control technique to an electro-hydraulic shaking table control system [18].Plummer reviewed and summarized the development history of the offline iterative control technique, and indicated its two main applications: vehicle experiments and seismic experiments [19].Underwood researched the application of the iterative control technique in random, sine and shock vibration experiment in depth [20].Bladh [21] carried out a virtual vehicle persistence vibration test combining offline iterative control and a vehicle software model.Eksteen [22] applied offline iterative control based on the nonlinear inverse model to time domain waveform replication.Yang [23] realized acceleration time domain waveform replication using offline iterative control based on singular value threshold control.
The traditional iterative control mentioned above usually adopts the same small iterative gain to ensure the system is stable, which reduces the convergence efficiency and requires more iterative times to obtain high control precision.Liu [24] proposed adaptive variable gain iterative control based on the frequency Fx-LMS algorithm and applied it to a shock vibration test.Cornelis [25] improved traditional offline iterative control, increasing the identification precision by using the optimal exciting signals and introducing a weighting correction factor to update system frequency response function.Moten [26] combined adaptive inverse control and offline iterative control to improve the control precision and convergence efficiency of shaking tables' time domain waveform replication.Cuyper [27] proposed hybrid control, combining offline iterative control and a real-time feedback controller based on H∞ robust control.Tang [28] proposed inner model control combined with offline iterative control to improve the iterative convergence and control precision.
Aiming at the phenomenon that a nonlinear mechanical system may be a non-minimum phase system and the system can be non-convergent using traditional iterative control, Smolders [29] proposed nonlinear iterative control based on the Newton optimization algorithm to improve the control precision and convergence speed of a four-channel road simulator.Daley [30] applied iterative control based on the gradient optimization algorithm to a road simulator to reduce the iterative times and improve the system robustness.Underwood [31] adopted the quasi-Newton optimization algorithm to update system impedance and iterative gain to realize sine sweep adaptive iterative control.Yang [32] proposed a single sine sweep control method based on the extended quasi-Newton optimization algorithm to correct system FRF to solve the potential stability problem.
Roberts [33] proposed iterative control based on the Broyden optimization algorithm and analyzed the global convergence to solve the slow convergence problem of simulation test equipment using traditional iterative control.Gunnarsson [34] proposed the frequency domain optimization iterative learning control method to solve the problem of minimizing the control error of a nonlinear system.Adali [35] researched the equivalence relation of the optimization algorithm in a real domain and a complex domain and the advantage in the complex domain, and then applied the complex domain optimization algorithm to the nonlinear function of signal processing.Brandwood [36] proposed an optimization algorithm based on a complex vector gradient operator to solve the real valve scalar optimization of the complex vector function.Hjorungnes [37] studied the complex domain optimization iterative algorithm and its application in signal processing.Kreutzdelgado [38] deduced the gradient operator and related complex differential algorithm in a complex domain.Sorber [39] studied the unconstrained optimization algorithm of complex vector real valve function and its related theory.
This paper presents a novel adaptive iterative control method based on the complex domain optimization algorithm to improve the system convergence rate and system stability.The complex domain quasi-Newton optimization algorithm is deduced based on the real domain quasi-Newton optimization algorithm.The complex steepest gradient iterative gain is proposed to obtain a better convergence rate, which can better compensate the phase error of the control system.The proposed adaptive iterative control can make the system stable even if the identified system model error is quite large, and make the system converge to the desired signal more quickly than the traditional iterative control.The proposed method can be introduced to the shaking tables system to efficiently realize the precise acceleration time domain waveform replication.
The remainder of this paper is organized as follows.Section 2 presents the investigated electro-hydraulic shaking tables system and the basic principle of the shaking tables' traditional iterative control method.In Section 3, the proposed complex domain adaptive iterative control scheme is introduced.Section 4 shows the application of the proposed complex domain optimization algorithm to shaking tables' adaptive iterative control.Experiments are carried out to validate the performance of the designed control method in Section 5, followed by Section 6, which briefly concludes the paper.

The Researched Shaking Tables System
The shaking tables system researched in this paper has two electro-hydraulic shaking tables, which are shown in Figure 1.The large-span structure test specimen is fixed on the two platforms.Each shaking table has eight servo valves, eight hydraulic cylinders, eight displacement sensors, eight accelerator sensors, one upper platform and one lower platform fixed in the base.Both the shaking tables can realize six DOF motions.Each hydraulic cylinder has one upper hinge connected with the upper platform and one lower hinge connected with the lower platform.The control system sends the control command to the servo valves, and then causes the hydraulic cylinder to move, then the upper platform drive test specimen to move according to the expected trajectory.Meanwhile, the feedback displacement and acceleration information is translated to the control system in real time to compensate the control signal.
one lower hinge connected with the lower platform.The control system sends the con trol command to the servo valves, and then causes the hydraulic cylinder to move, then the upper platform drive test specimen to move according to the expected trajectory Meanwhile, the feedback displacement and acceleration information is translated to the control system in real time to compensate the control signal.

Offline Iterative Learning Control for Shaking Tables
The traditional offline iterative learning control method for shaking tables is shown in Figure 2. It can be seen that the total process is mainly composed of three phases: ex periment preparation, which includes experiment data preprocessing, software launch and the preparation of the hydraulic power supply and hardware equipment; system model identification, which builds the linear system FRF model and then computes the inverse model; and target simulation, which iteratively computes the driven signal ac cording to identified FRF model and desired target signal.Generally, the input signal of the identified system is accurately known, so the sys tem identification method usually applies the H1 identification method to evaluate the system FRF model: where [ ] ( ) xx f S is the self-spectrum matrix of the input signal and

S
is the cross-spectrum matrix of the input signal and output signal.The process of the target simulation stage is as follows:

Offline Iterative Learning Control for Shaking Tables
The traditional offline iterative learning control method for shaking tables is shown in Figure 2. It can be seen that the total process is mainly composed of three phases: experiment preparation, which includes experiment data preprocessing, software launch and the preparation of the hydraulic power supply and hardware equipment; system model identification, which builds the linear system FRF model and then computes the inverse model; and target simulation, which iteratively computes the driven signal according to identified FRF model and desired target signal.
one lower hinge connected with the lower platform.The control system sends the con trol command to the servo valves, and then causes the hydraulic cylinder to move, then the upper platform drive test specimen to move according to the expected trajectory Meanwhile, the feedback displacement and acceleration information is translated to the control system in real time to compensate the control signal.

Offline Iterative Learning Control for Shaking Tables
The traditional offline iterative learning control method for shaking tables is shown in Figure 2. It can be seen that the total process is mainly composed of three phases: ex periment preparation, which includes experiment data preprocessing, software launch and the preparation of the hydraulic power supply and hardware equipment; system model identification, which builds the linear system FRF model and then computes the inverse model; and target simulation, which iteratively computes the driven signal ac cording to identified FRF model and desired target signal.Generally, the input signal of the identified system is accurately known, so the sys tem identification method usually applies the H1 identification method to evaluate the system FRF model: where [ ]  Generally, the input signal of the identified system is accurately known, so the system identification method usually applies the H1 identification method to evaluate the system FRF model: where [S xx ( f )] is the self-spectrum matrix of the input signal and S xy ( f ) is the crossspectrum matrix of the input signal and output signal.The process of the target simulation stage is as follows: (1).Apply FFT to the time domain reference signal {r(t)} to obtain the frequency domain reference signal {R( f )}, compute the driven signal in the frequency domain according to the identified system FRF model [H 0 ( f )], then carry out IFFT to obtain the original time domain driven signal {d 0 (t)}: (2).Excite the system with {d 0 (t)}, acquire the response signal {c 0 (t)} and compute the control error signal {e 0 (t)}: (3).Determine whether {e 0 (t)} meets the precision requirement.If the condition is met, the iterative process will terminate; if the condition is not met, carry out FFT to {e 0 (t)} and compute the frequency domain compensated driven signal {∆D( f )}: (4).Carry out IFFT to {∆D( f )} and compute the 1st iteratively driven signal {d 1 (t)} with iterative gain α: (5).Excite the system with the new driven signal {d 1 (t)}, acquire the response signal {c 1 (t)} and compute the control error signal {e 1 (t)}: . Determine whether {e 1 (t)} meets the precision requirement.If the condition is met, the iterative process will terminate; if the condition is not met, repeat steps (3)~( 5).The iterative process will terminate until the control error signal {e k (t)} meets the desired precision requirement, where k = 1, 2 . . . .

Influence of Identification Precision on Control System
It can be seen from the above that the identified FRF model directly participates in the control process and plays an important role in the control system, so this section will study the influence of system identification precision on offline iterative control.The identified FRF model is just the approximation of the true system due to nonlinearity, uncertainty and identification method error.From this, the simplified iterative process shown in Figure 2 considering the identified FRF model error [∆H( f )] is shown in Figure 3.

Influence of Identification Precision on Control System
It can be seen from the above that the identified FRF model directly participates the control process and plays an important role in the control system, so this section w study the influence of system identification precision on offline iterative control.Th identified FRF model is just the approximation of the true system due to nonlinearit uncertainty and identification method error.From this, the simplified iterative proce shown in Figure 2 considering the identified FRF model error In Figure 3, the solid line indicates the real-time control part and the dotted line indicates the offline control part.{R( f )} is the frequency domain reference signal vector, {C k ( f )} is the kth iterative frequency domain response signal vector, {D k ( f )} is the kth iterative frequency domain driven signal vector, {∆D( f )} is the frequency domain compensated driven signal vector, [H( f )] is the true system model, [H 0 ( f )] is the identified system model and its inverse model is the impedance matrix where [I n ] is the n dimensional unit matrix and n is the dimension of the control system.From Figure 3, the driven signal update function can be obtained: Multiplying both sides of Equation ( 8) by [H( f )], it can be derived that From Equation (10), the following can be derived: From Equation (11), it can be seen that if the norm of control error {E k ( f )} during the iteration process decreases with the increase in iterative time, then the iterative control is convergent, and the following condition of {E k ( f )} tending to zero can be expressed as The equivalent condition of Equation ( 12) is From Equation (13), it can be seen that the convergence and convergent rate are determined by the iterative gain α and identified FRF error [∆H( f )].The iteration process will be unstable if [∆H( f )] is quite large in some frequency points.As a result, to improve the iterative control performance, the optimal iterative gain needs to be chosen and the identified FRF needs to be updated, which makes the iterative control have a faster convergent rate while keeping the system stable.

The Complex Optimization Algorithm-Based Iterative Learning Control
As is mentioned above, to solve problem of offline iterative learning control, it is necessary to adapt different iterative gains in different frequency points.Thus, to obtain a better control effect, the optimal iterative gain and system impedance update are developed during the iteration process.The purpose of iterative control is to minimize the control error {R( f )} − {C k ( f )} .To solve the potential stability and convergence problem of iterative control, this paper proposed an adaptive iterative control based on the complex optimization algorithm, which makes the response signal {C k ( f )} close to the desired signal {R( f )} to reduce {R( f )} − {C k ( f )} through optimize iterative gain {α k } and update system impedance [J 0 ( f )].Due to the fact that the computation is realized in the complex domain, which contains magnitude and phase, this method belongs to complex optimization theory.
The essence of the MIMO iterative control problem for shaking tables is the minimum of the object function described using the L2 norm of control error: The purpose of iterative control is to find the optimal complex driving signal D( f ) * , which makes Equation ( 14) the minimum.However, the electro-hydraulic shaking tables system is a typical nonlinear system, which means [H( f )] is nonlinear and is the function of the driving signal {D( f )}.Thus, the problem of this paper can be considered to solve the nonlinear norm function shown in Equation ( 14), whose optimization is the function of complex vector {D( f )}, so this paper will research the optimization theory in the complex domain.

The Real Optimization Algorithm Theory
Assuming the dimensions of the input and output signals are both n, the identified system FRF matrix [H 0 ( f )] and the corresponding impedance matrix [J 0 ( f )] are both n × n complex matrixes.The real domain optimization theory is relatively complete and mature.However, the complex domain optimization algorithm is still incomplete; thus, this paper will research the complex optimization algorithm based on the real optimization algorithm and the mapping relationship between real vector space and complex vector space.
The local unconstrained minimum of nonlinear function F(x) in real vector space is given by where {x} is the n dimension real vector in R n space and g(•) is the nonlinear mapping of R n → R n .Assume that g(•) is twice continuously differentiable to {x}, then the goal of the optimization problem is to find optimal {x * } to minimize the objective function F(x).This problem belongs to optimization theory in real vector space and can be solved by applying the following theorem.
Theorem 1. Letting G : Ω → R 1 be a real valued function defined on a closed and bounded region Ω in R n , (1) if {x 0 } is an inner point of Ω, at which G(•) achieves an external value, and if G(•) has a continuous first partial derivative, then ∂G/∂x 1 = ∂G/∂x 2 = . . .= ∂G/∂x n = 0, that is, the partial derivative of G(•) is zero at the optimal vector {x 0 }.
(2) if G(•) has a continuous second-order partial derivative in a neighborhood N 0 of any inner point {x 0 } of Ω, then for an arbitrary {x} in N 0 , there is the quadratic Taylor representation of G(•): Theorem 1 is the basic principle to solve the unconstrained optimization problem.The common method is using the quadratic function given by Theorem 1(2) and the optimal condition given by Theorem 1(1) to define optimal vector, in which there is ∂G/∂x t {x 0 } = 0, for Equation (15), F(•) = G(•).At an arbitrary vector near the optimal vector {x 0 }, exe- cuting the Taylor series expansion to ∂G/∂x t , then the approximate linear equation is given by Equation ( 16) is commonly used to solve the optimization problem by generating a sequence of vectors {x k }, which converge to the optimal vector to minimize F(•).The methods used to solve Equation ( 16) are Newton-Gauss, Newton-Raphson and quasi-Newton, etc.The Newton-Raphson is an efficient algorithm that uses the gradient and Hessian matrix of an objective function to find the minimization sequence.The minimizing sequence {x k } is given by the following recursion: where α k is the iteration step to minimize F({x k } + α k {P k }) and {P k } is the Newton direction for the nonlinear optimization problem and is described as where [H k ] is the Hessian matrix of F(•) to {x} and [J k ] is the Jacobian matrix of g(•) to {x}.
For the Newton iterative optimization method, there is the following local convergence theorem.
Then, there is an ε > 0 such that for all x 0 ∈ N(x * , ε), the sequence generated by The G(•) in Theorem 2 is the gradient of F(•) in Equation ( 2) and the Jacobian of G(•) is the Hessian matrix of F(•).Submit F(•) into the Newton direction of the nonlinear norm function shown in Equation (17), and it can be obtained that It can be seen from Equation ( 18) that {P k } is the product of the inverse matrix of the Hessian matrix and the gradient matrix of F(•).If the first item of the Hessian matrix can be ignored, then submit the approximate Hessian matrix into Equation (19), and there is Submitting Equation ( 20) into (17), there is the Newton-Gauss optimization algorithm, which has global convergence.
From the above analysis, it can be seen that it is necessary to compute the gradient and Hessian of the function to solve the optimization problem.To simplify this problem, some approximate Secant and Hessian techniques are studied.The most successful approximate method is the Broyden method, whose basic idea is to use an estimate of the Jacobian matrix in Theorem 2 instead of the actual Jacobian matrix, as is shown in following.
Let G : R n → R n and x 0 ∈ R n , A 0 ∈ R n×n , then there are the following vector and matrix update equations: It is similar to the Newton-Gauss method, but it uses [A k ] instead of the Jacobian matrix [J(x k )] at {x k }.There is a similar global convergence theorem for this method.Theorem 3. Let the assumptions in Theorem 2 exist, and there is ε, δ > 0, which makes x 0 − x * 2 ≤ ε and A 0 − J(x * ) ≤ δ, then sequence {x k } generated by Equation ( 22) con- verges q-superlinearly to {x * }.
If G : R n → R n is affine, then the quasi-Newton method given by Equation ( 22) converges finitely and {x 2k+1 } = {x * }, that is, this method converges to optimal vector {x * } in no more than 2k + 1 iterations.

The Complex Optimization Algorithm
This paper will research the optimization algorithm in complex space and the local unconstrained minimum to nonlinear function F : C n → R 1 given by where Let z ∈ C 1 be the complex plane z = x r + jx i and g : C 1 → C 1 be continuously differentiable as a function of R 2 , then g : C 1 → C 1 is called the C 1 function and its derivative with respect to complex z and conjugate complex z is, respectively, Assuming that G(z 1 , . . ., z i , . . ., According to Equation ( 25), the complex gradient and Hessian matrix with respect to {z} and {z} of function F(•) : C n → R , which is twice continuously differentiable, it can be deduced that For function g(•) : C n → R , its Taylor series of complex derivatives at complex vector {z a } is Equations ( 29) and ( 30) are, respectively, the holomorphic and non-holomorphic Hessian.For g(•) : C n → R , if the non-holomorphic Hessian is positive definite at {z a }, then g(•) is strictly plurisubharmonic.
Let U ⊂ C n be open and g(•) : C n → R be strictly plurisubharmonic in U, and for w ⊆ U, there is c > 0, ε > 0 and the following inequation exists: Equation (32) indicates that when L(•, z) = 0 and g(•) is strictly plurisubharmonic, there is g(z a + λ) ≥ g(z a ) + c λ 2 , λ < ε, which means g(•) has a strict local minimum at {z a }.
According to Theorem 2, it can be seen that when the start point is close enough to optimum vector, it can be used to find the optimal complex vector to minimize the second-order approximation iteratively shown in Equation (28).Equation (32) shows that the holomorphic Hessian can be ignored in a neighborhood of optimal vectors, which makes the Newton-Raphson update method for complex vector space.
According to Theorem 1, the complex and conjugate complex gradients vanish if the real gradient vanishes, which is the classic necessary condition that optimal vector must satisfy.Describe F(•) using the approximate second order polynomial in Equation ( 28) ignoring the last item o(|λ| 3 ), compute the complex and conjugate complex gradients of Equation ( 30) and set the results to zero, then the complex Newton-Raphson update algorithm is obtained, which is similar to the real algorithm, in which the update sequence {z} satisfies Equations ( 33) and ( 34) and converges to the optimal point.
Assume the holomorphic Hessian of F(•) is zero, then there is According to the mapping relationship of real vector space and complex vector space shown in Equations ( 24)-( 27), and Theorems 1 and 2, the local convergence of the complex vector Newton-Raphson update algorithm shown in Equation ( 35) can be proven.
The modification to the complex Newton-Raphson algorithm making the algorithm have global convergence is shown in the following: where F(•) is the steepest descent gain to minimize F(z k+1 (α k )).
Combining Equation (36) with the minimum norm shown in Equation ( 23), in which g : C n → C n is the holomorphic function, the gradient of F(z) = g(z) 2 can be obtained: The holomorphic and conjugate non-holomorphic Hessian is given by From Equation (38), it can be seen that the holomorphic Hessian is zero when the holomorphic Hessian of g(•) is zero or g(•) is zero, or the dot product in Equation ( 38) is zero.For the iterative control problem of multi-dimension shaking tables shown in Equation ( 14), g(•) = 0 can be satisfied at the optimal point of F(•), because g(•) is the control error vector, and it is zero at the global optimal vector of F(•).
The analysis above indicates the Broyden method can solve the optimization problem of the nonlinear system well; thus, we will extend the real vector space Broyden method to the complex vector space Broyden.
Let g : C n → C n , z 0 ∈ C n and A 0 ∈ C n×n , and there is the complex Broyden update algorithm: The proof of the global convergence of the complex Broyden update algorithm can follow from converting the complex formulation to the real vector formulation and real Broyden method convergence theorem shown in Theorem 3. When g(•) is conjugate affine, the complex Broyden method converges to the optimal vector {z m } within 2k + 1 iterative steps, that is, the algorithm is finitely convergent.

The Complex Steepest Descent Algorithm
The traditional method for the optimization algorithm is using the positive steepest descent gain to ensure convergence.However, for the complex quasi-Newton method, the complex steepest descent gain would be more satisfactory.Compared with the classical necessary condition [∂F/∂z]| {z m } = 0 at optimal vector {z m }, Equation (25) has the simpler necessary condition g(z)| {z m } = 0. To utilize the nature of the optimization problem, the descent condition on the root finding algorithm in the condition is used to obtain the sequence {z k }: The way to ensure that condition (40) satisfied is to use the steepest descent approach.For Equation (39), the necessary condition is Consider F(•) as function of α, and describe it with Taylor series, then there is where The complex quantities of the steepest descent condition shown in Equation (41) are Solve Equation (43), and the steepest descent gain can be obtained: If the holomorphic Hessian of f (•) is approximately zero, Equation ( 44) is simplified as According to the optimal iterative step s k = −α k A −1 k g(z k ), applying the complex steepest descent, there is Equation (47) indicates that the complex steepest descent and change in g(•) are orthogonal, that is, the complex descent technique produces a monotonically decreasing sequence of values for the objective function.

The Adaptive Iteration Learning Control for Shaking Tables
The shaking tables' iterative driven signal update function shown in Equation ( 8) is similar to the steepest descent method of the complex Newton-Raphson algorithm, whose object is to minimize the object function The difference between Equation ( 8) and traditional iterative control is that the complex matrix [J 0 ( f )] and complex error {R( f )} − {C( f )} replace the real positive Hessian and real gradient.Due to this, F(•) in equation ( 48) is a plurisubharmonic function and is the function of {D( f )} and D( f ) .Thus, this paper can apply the complex optimization algorithm proposed previously to solve the problem of shaking tables' iterative control.

Application of Complex Optimization Algorithm
This paper researches the optimization problem of the complex vector realvalued function, and the problem can be described using the complex vector error 48) is a plurisubharmonic function.Therefore, we can use the optimization theory in complex vector space proposed before to obtain the driven signal vector sequence {D( f )}, which converges to optimal vector {D m ( f )} to minimize the objective function F(•).
The nonlinear function g(•) of Combining Equations ( 38) and ( 49), the conjugate-transposed gradient of F(•) can be obtained, so that Combining Equations ( 39) and (50), the complex conjugate of the non-holomorphic Hessian of F(•) can be obtained, so that Combining Equations ( 40) and ( 51), the holomorphic Hessian of F(•) can be obtained, so that According to Equation (52), the following algorithm shown in Equations ( 53) and ( 54), which converges to {D m ( f )}, can be obtained: According to Equations ( 50)-(54), the global Newton-Raphson update algorithm of shaking tables' iteration control can be obtained, so that Assume [H( f )] is a holomorphic function of {D( f )}, and the complex conjugate gradient can be obtained according to Equation ( 38): Additionally, the non-holomorphic Hessian can be obtained, so that Let the holomorphic Hessian be approximately zero, and the simplified complex steepest descent Newton update algorithm can be obtained according to Equations ( 53)-(57): Apply the proposed complex Broyden update algorithm to the iterative control to obtain the optimal vector that makes the objective function vanish.In this case, the approximate Jacobian becomes an approximate Hessian, and the approximate inverse Hessian is described as Instead the real Hessian with approximate Hessian in Equation ( 59), there is Comparing Equations ( 8) and (60), we note that the traditional iterative control algorithm is similar to the proposed complex Broyden update algorithm.However, it is only similar in form, and the difference is the approximate impedance sequence [J k ( f )] instead of initial [J 0 ( f )], and complex gain sequence α k instead of the real non-optimized iterative gain α.The complex gain can compensate phase error in impedance estimation.
According to Equation ( 39), the approximate impedance [J k ( f )] and control system {y k ( f )} can be obtained, so that where The above algorithm produces an adaptive iterative control, but this method cannot ensure the approximate Hessian [A n ( f )] positive definite, which is needed to ensure the descent direction for this method.The aforementioned complex steepest descent can solve this problem, which can ensure the descent direction, that is, the decrease direction of the control error.
To effectively use the steepest descent method, we need to find a complex α k minimizing the next control error F({D k+1 ( f )}), and there is F({D k+1 ( f )}) ≤ F({D k ( f )}).The shaking tables system is holomorphic, and combining Equation (60): Combining Equations ( 42), ( 43) and (47), it can be derived that If the objective function given by Equation ( 47) is approximately Hermitian quadratic and the system is holomorphic, then calculating the complex derivative of F({D k+1 ( f )}) in Equation (64) with respect to α k and setting it to zero, the optimal α k can be obtained, so that The control process described in Equation (62) using complex steepest descent gain reduces the control error at each iterative control loop, and the process converges in, at most, 2k + 1 iteration steps for a linear system.The algorithm given by Equations ( 60) and ( 61) is a locally q-superlinear convergent for a nonlinear system, and in this case, the algorithm will track system impedance as it varies with drive and control response amplitudes if the control and impedance errors are small.

Realization of Adaptive Iteration Learning Control for Shaking Tables
From Equation (64), it can be seen that {ξ k ( f )} and [H( f )] are unknown, and the optimal gain α k depends on unknown [H( f )].The two-step control loop iteration method is applied to solve this problem, where a learning loop is used to measure {ξ k ( f )} to calculate α k .The obtained α k is used to correct the control error in a control loop, and then to update the impedance matrix estimation, followed by the next learning loop iteration, and so forth until the control precision is satisfactory.
The traditional method given is used to obtain the driven vector of the learning loop: Due to learning driven vector change being relatively small, Equation (66) holds.According to Equation (63), the approximate response in the learning loop is From Equation (67), it can be seen that unknown {ξ k ( f )} can be described using a measurable system response.Combining Equation (46), the optimal gain α k only depending on measurable control system quantities can be obtained, so that The problem of complex optimization iteration control for shaking tables has been solved and all of the quantities used to update [J k ( f )] and {D k ( f )} are available in the control process.The method given by Equations ( 60), ( 61), ( 66) and (68) forms an adaptive iteration learning control method based on the complex optimization algorithm and complex steepest descent for shaking tables, and the total control process diagram is shown in Figure 4.
Electronics 2023, 12, x FOR PEER REVIEW 17 of 2  From Figure 4, it can be seen that the control process contains three main parts: an initial iterative step, iterative loop and learning loop.In this process, the iterative loop and learning loop are performed alternately until the control precision at each frequency is satisfactory.After the control precision is satisfied at a particular frequency, then th driven spectrum at this frequency does not change at all and the iteration repeat is only still performed at other frequency points.The repeat control can realize the fast and sta From Figure 4, it can be seen that the control process contains three main parts: an initial iterative step, iterative loop and learning loop.In this process, the iterative loop and learning loop are performed alternately until the control precision at each frequency is satisfactory.After the control precision is satisfied at a particular frequency, then the driven spectrum at this frequency does not change at all and the iteration repeat is only still performed at other frequency points.The repeat control can realize the fast and stable adaptive optimization iteration control of the shaking tables system and accurately reproduce the acceleration time domain waveform.

Experiment Results and Analysis
The experiment system was built to test the researched content above.The total control system is shown in Figure 5, and mainly contains a hydraulic power manager unit, manager PC, motion controller, emergency button, target PC, signal condition unit, PLC, hydraulic power, shaking tables and test specimen.The main parameters of the experiment system are listed in Table 1.The real-time control experiment system is based on the xPC rapid prototype control technology, in which the system sample time is 0.001 s.The hydraulic power is controlled by PLC, which is communicated with through the RS485.The host manager PC is communicated with through the network.The real-time control system is built applying Matlab/Simulink, which is compiled to the real-time control

Tests Using Exactly Identified FRF Model
To make comparisons, two series of experimental tests were performed through adopting two kinds of identified system FRF models: small and large estimation errors.Each series of tests compares the control effects of the traditional offline iterative control method usually used in shaking table control in practice and the proposed adaptive iterative learning control method.The iteration control algorithm of the shaking tables system is realized in the frequency domain, whose control essence is the reduction in control error in the frequency domain at each frequency.The signal described in the frequency domain is a complex that contains two factors: a real part and imaginary part.Table 1.Main parameters of the shaking tables system.

Parameters
Values/Units mass of platform 4000/kg mass of load 13,000/kg effective area of cylinder 0.0075/m 2 density of hydraulic oil 845/kg/m 3 nature frequency of servo valve 120/Hz rated flow of servo valve 400 L/min pressure of hydraulic source 25 × 106/Pa

Tests Using Exactly Identified FRF Model
To make comparisons, two series of experimental tests were performed through adopting two kinds of identified system FRF models: small and large estimation errors.Each series of tests compares the control effects of the traditional offline iterative control method usually used in shaking table control in practice and the proposed adaptive iterative learning control method.The iteration control algorithm of the shaking tables system is realized in the frequency domain, whose control essence is the reduction in control error in the frequency domain at each frequency.The signal described in the frequency domain is a complex that contains two factors: a real part and imaginary part.Thus, some signals and control parameters during the control process will be described in complex vector form in this section.The experiment results and detailed comparison analysis are shown in the following part of the paper.
It should be noted that shaking tables are mainly used to simulate earthquake vibration environments, whose motion waveforms are usually three translations along the X, Y and Z axes.Therefore, the tests and analysis are also implemented in three translation directions in this paper.Due to the shaking tables test being mainly an earthquake waveform test, which performs the characteristic of randomness, this section presents tests for various cases using an earthquake waveform random signal, whose amplitude is about ±0.4 g, frequency range is 0.1~20 Hz, duration time is 20 s and control channel is six, corresponding to directions X1, Y1, Z1, X2, Y2, Z2 of the shaking tables, as shown in Figure 6.This section presents the tests using the traditional iteration learning control and proposed adaptive iteration learning control method with relatively accurate FRF esti mation, respectively.The amplitude and phase of the real and estimated FRFs are, re spectively, shown in Figures 7 and 8, in which blue curves are the real FRF model and red curves are the estimated FRF model.It can be seen that the errors between the rea FRF model and estimated FRF model are relatively small, especially in the frequency range 0.1~20 Hz, which is the main frequency of the reference signal.This section presents the tests using the traditional iteration learning control and proposed adaptive iteration learning control method with relatively accurate FRF estimation, respectively.The amplitude and phase of the real and estimated FRFs are, respectively, shown in Figures 7 and 8, in which blue curves are the real FRF model and red curves are the estimated FRF model.It can be seen that the errors between the real FRF model and estimated FRF model are relatively small, especially in the frequency range 0.1~20 Hz, which is the main frequency of the reference signal.This section presents the tests using the traditional iteration learning control and proposed adaptive iteration learning control method with relatively accurate FRF estimation, respectively.The amplitude and phase of the real and estimated FRFs are, respectively, shown in Figures 7 and 8, in which blue curves are the real FRF model and red curves are the estimated FRF model.It can be seen that the errors between the rea FRF model and estimated FRF model are relatively small, especially in the frequency range 0.1~20 Hz, which is the main frequency of the reference signal.

Experiments Using Traditional Iteration Learning Control Method
This section will present experiments using traditional iteration control with dif ferent iteration gains.The system control errors at the 300th frequency point (which is

Experiments Using Traditional Iteration Learning Control Method
This section will present experiments using traditional iteration control with different iteration gains.The system control errors at the 300th frequency point (which is 14.947 Hz) with iteration gain 0.2 are shown in Figure 9.The time domain control errors after 20 iterations are shown in Figure 10.The system RMSEs (root mean square errors) with the increase in iteration time are shown in Figure 11.The absolute RMSE and the relative RMSE are, respectively, defined as follows:  The absolute RMSE and the relative RMSE are, respectively, defined as follows: where r(i) is the reference signal sequence, c(i) is the response signal sequence and RMS{•} and mean{•}, respectively, describe the root mean square and mean value of {•}.
From Figures 9-11, it can be seen that the convergence rates are slow when the iteration gain is 0.2; the system RMSEs are about 10% after 10 iterations and are about 0.01 g after 20 iterations, which is about 5% relative to the reference signal RMSs about 0.2 g.Furthermore, it can be seen from Figure 11 that the difference in control errors between different channels in each iteration step are large.Note that the control errors after the 15th iteration do not decrease at all, which is because the system inevitably includes various factors such as electrical noise, pressure fluctuation of the hydraulic power supply and so on.
The control errors at the 300th frequency point with iteration gain 0.7 are shown in Figure 12.The time domain control errors after 20 iterations are shown in Figure 13.The RMSEs with the increase in iteration time are shown in Figure 14.From Figures 12-14, it can be seen that the convergence rates become fast when the iteration gain is 0.7; the system RMSEs are about 0.01 g after five iterations, which are about 5% relative to the reference signals.The system control errors converge to the minimized control errors after five iterations and the errors will not decrease at all with the increase in iteration time.From Figures 9-11, it can be seen that the convergence rates are slow when the iteration gain is 0.2; the system RMSEs are about 10% after 10 iterations and are about 0.01 g after 20 iterations, which is about 5% relative to the reference signal RMSs about 0.2 g.Furthermore, it can be seen from Figure 11 that the difference in control errors between different channels in each iteration step are large.Note that the control errors after the 15th iteration do not decrease at all, which is because the system inevitably includes various factors such as electrical noise, pressure fluctuation of the hydraulic power supply and so on.
The control errors at the 300th frequency point with iteration gain 0.7 are shown in Figure 12.The time domain control errors after 20 iterations are shown in Figure 13.The RMSEs with the increase in iteration time are shown in Figure 14.From Figures 12-14, it can be seen that the convergence rates become fast when the iteration gain is 0.7; the system RMSEs are about 0.01 g after five iterations, which are about 5% relative to the reference signals.The system control errors converge to the minimized control errors after five iterations and the errors will not decrease at all with the increase in iteration time.

Experiments Using Adaptive Iteration Learning Control Method
This section will present tests applying the proposed adaptive optimization iteration learning control.The learning loop iteration gains are set to the same value, 0.2, in six directions, and the threshold value of control error during the iteration process is set to 10 −5 , that is, when the ratio of the control error and reference signal is less than 10 −5 , the system impedance and the complex steepest descent gain stop updating.The control errors, complex steepest descent gain and the main diagonal elements of the system impedance matrix at the 300th frequency point with the increase in iteration time are, respectively, shown in Figures 15-17.The time domain control errors after 20 iterations are shown in Figure 18 and the system RMSEs with the increase in iteration time are shown in Figure 19.

. Experiments Using Adaptive Iteration Learning Control Method
This section will present tests applying the proposed adaptive optimization iteration learning control.The learning loop iteration gains are set to the same value, 0.2, in six directions, and the threshold value of control error during the iteration process is set to 10 −5 , that is, when the ratio of the control error and reference signal is less than 10 −5 , the system impedance and the complex steepest descent gain stop updating.The control errors, complex steepest descent gain and the main diagonal elements of the system impedance matrix at the 300th frequency point with the increase in iteration time are, respectively, shown in Figures 15-17         From Figures 15-19, it can be seen that the control errors converge to about 0.01 g after only three iterations applying the proposed adaptive optimization iteration learning control method.Compared with traditional iteration control, with an iteration gain of 0.7, the proposed control method has a faster convergence rate, which is because it From Figures 15-19, it can be seen that the control errors converge to about 0.01 g after only three iterations applying the proposed adaptive optimization iteration learning control method.Compared with traditional iteration control, with an iteration gain of 0.7, the proposed control method has a faster convergence rate, which is because it adopts different optimal iteration gains at each frequency point.Furthermore, from Figures 11,14 and 19, it can be seen that the differences between each channel's RMSEs are relatively small compared with traditional iteration control.This is because the former treats all channels of shaking tables as an entirety, and the iteration gain is calculated based on the entire system optimization algorithm.However, the traditional iteration control method applies the same not-optimized iteration gain and the gain makes one channel converge easily and the others converge poorly.Furthermore, with the convergence of control errors, the real part of the optimized iteration gain tends to one and the imaginary part tends to zero.The system impedance matrix does not converge to the real value, but the actual value will tend to the real value with the increase in iteration time, which is because the purpose of the proposed adaptive iterative learning control method is to minimize the system control error and not to obtain a precise system impedance value.The experiment results indicate that the proposed adaptive optimization iteration learning control method based on the proposed complex optimization algorithm and complex steepest descent algorithm has a faster convergence rate compared with the traditional iteration control method.

Tests Using Large Error-Identified FRF Model
This section presents tests using the traditional iteration learning control method and proposed adaptive optimization iteration learning control method, respectively, when the system FRF model estimation error is relatively large.The amplitude and phase of the real and estimated FRFs are, respectively, shown in Figures 20 and 21, in which blue curves are the real FRF model and red curves are the estimated FRF model.It can be seen that the errors between the real FRF model and estimated FRF model are large in the frequency range 0.1~20 Hz, which is the main frequency of the reference signal.In particular, the maximum of the phase error is about 180 • , so it is very easy to make the control system unstable, as mentioned above.
Electronics 2023, 12, x FOR PEER REVIEW large in the frequency range 0.1~20 Hz, which is the main frequency of th signal.In particular, the maximum of the phase error is about 180°, so it is v make the control system unstable, as mentioned above.

Experiments Using Traditional Iteration Learning Control Method
Corresponding to Section 5.1.1,this section will also present tests using the traditional iteration learning control method.The system control errors at the 300th frequency point with iteration gain 0.1 are shown in Figure 22.The system time domain control errors after 20 iterations are shown in Figure 23.The system RMSEs with the increase in iteration time are shown in Figure 24.Note that the iteration gain 0.1 (not 0.2 or 0.7 in this part) is chosen to ensure that the test system is relatively stable and controllable and is not destroyed.

Experiments Using Traditional Iteration Learning Control Method
Corresponding to Section 5.1.1,this section will also present tests using tional iteration learning control method.The system control errors at the 300t cy point with iteration gain 0.1 are shown in Figure 22.The system time dom errors after 20 iterations are shown in Figure 23.The system RMSEs with the iteration time are shown in Figure 24.Note that the iteration gain 0.1 (not 0 this part) is chosen to ensure that the test system is relatively stable and contro is not destroyed.

Experiments Using Traditional Iteration Learning Control Method
Corresponding to Section 5.1.1,this section will also present tests using the traditional iteration learning control method.The system control errors at the 300th frequency point with iteration gain 0.1 are shown in Figure 22.The system time domain control errors after 20 iterations are shown in Figure 23.The system RMSEs with the increase in iteration time are shown in Figure 24.Note that the iteration gain 0.1 (not 0.2 or 0.7 in this part) is chosen to ensure that the test system is relatively stable and controllable and is not destroyed.From Figures 22-24, it can be seen that the control system is not convergent even if using the relatively small iteration gain of 0.1, which is because the error between the actual and estimated system FRF is relatively large, and the phase error is larger than 180° at some frequency points.In this case, the applied iteration gain and the system FRF estimation error cannot satisfy the stable condition of the iteration control shown in Equation (13).The compensation control in this case is equivalent to a reversed compensation, that is, this compensation control makes the control errors increase along with the increase in the iteration time.The experiment results indicate the validity of the traditional iteration control stability analysis.

Experiments Using Adaptive Iteration Learning Control Method
Corresponding to Section 5.1.2,this section presents tests applying the proposed adaptive optimization iteration learning control method.The learning loop iteration gains are set to the same value, 0.2, in six directions and the threshold value of the control error during iteration process is set to 10 −5 .The control errors at the 300th frequency point with the increase in iteration time are shown in Figure 25, the complex steepest descent gain at the 300th frequency point is shown in Figure 26 and the main diagonal elements of the system impedance matrix at the 300th frequency point are shown in Figure 27.The time domain control errors after 20 iterations are shown in Figure 28 and the system RMSEs with the increase in iteration time are shown in Figure 29.From Figures 22-24, it can be seen that the control system is not convergent even if using the relatively small iteration gain of 0.1, which is because the error between the actual and estimated system FRF is relatively large, and the phase error is larger than 180 • at some frequency points.In this case, the applied iteration gain and the system FRF estimation error cannot satisfy the stable condition of the iteration control shown in Equation ( 13).The compensation control in this case is equivalent to a reversed compensation, that is, this compensation control makes the control errors increase along with the increase in the iteration time.The experiment results indicate the validity of the traditional iteration control stability analysis.

Experiments Using Adaptive Iteration Learning Control Method
Corresponding to Section 5.1.2,this section presents tests applying the proposed adaptive optimization iteration learning control method.The learning loop iteration gains are set to the same value, 0.2, in six directions and the threshold value of the control error during iteration process is set to 10 −5 .The control errors at the 300th frequency point with the increase in iteration time are shown in Figure 25, the complex steepest descent gain at the 300th frequency point is shown in Figure 26 and the main diagonal elements of the system impedance matrix at the 300th frequency point are shown in Figure 27.The time domain control errors after 20 iterations are shown in Figure 28 and the system RMSEs with the increase in iteration time are shown in Figure 29.From Figures 25-29, it can be seen that the system RMSEs converge to about 0.023 g, which is about 15% relative to the reference signal, after five iterations using the proposed adaptive optimization iteration learning control method.The system RMSEs converge to about 0.02 g after 15 iterations, which is about 10%, and the control errors converge with the increase in iteration time.Compared with the results applying the traditional iteration learning control method shown in Figures 22-24, the proposed control method has successfully solved the potentially unstable problem of the traditional iteration learning control method when the system FRF estimation error is relatively large in some frequency points.Furthermore, with the convergence of control errors, the real part of the optimized iteration gain tends to one and the imaginary part tends to zero during the iteration process.The system impedance matrix does not converge to the real value, but the actual value will tend to the real value with the increase in iteration time, which is similar to the results in Section 5.1.2.From Figures 25-29, it can be seen that the system RMSEs converge to about 0.023 g, which is about 15% relative to the reference signal, after five iterations using the proposed adaptive optimization iteration learning control method.The system RMSEs converge to about 0.02 g after 15 iterations, which is about 10%, and the control errors converge with the increase in iteration time.Compared with the results applying the traditional iteration learning control method shown in Figures 22-24, the proposed control method has successfully solved the potentially unstable problem of the traditional iteration learning control method when the system FRF estimation error is relatively large in some frequency points.Furthermore, with the convergence of control errors, the real part of the optimized iteration gain tends to one and the imaginary part tends to zero during the iteration process.The system impedance matrix does not converge to the real value, but the actual value will tend to the real value with the increase in iteration time, which is similar to the results in Section 5.1.2.
Combining the experiment results in this and the above sections, it can be seen that the proposed adaptive optimization iteration learning control method based on the complex optimization algorithm and complex steepest descent algorithm makes the control system convergent no matter if the system FRF model estimation error is small or large.The proposed control method has a stronger robustness to system FRF model estimation error, and it has a faster convergence rate when applying the same system FRF estimation model compared with the traditional iteration learning control method, which indicates that the proposed adaptive optimization iteration learning control method is effective and feasible, and can accurately realize the acceleration time domain waveform replication of shaking tables systems.

Conclusions
In this paper, the adaptive iterative learning control method based on the complex optimization algorithm was proposed to overcome the problem of the traditional iteration learning control method.The complex Broyden optimization algorithm combined with complex steepest descent gain is deduced according to the real Broyden optimization algorithm and the mapping relationship between real vector space and complex vector space.The shaking tables experiment system is built based on xPC target rapid prototyping technology.The tests performed use the proposed adaptive iterative learning control and traditional iterative learning control, respectively, under two conditions: with relatively accurate system FRF estimation and with a large system FRF estimation error.The percent RMSE of the reference signal (RMS is about 0.2 g) can decrease to 5% (about 0.01 g) within three iterations using the proposed control method and at least 20 iterations using the traditional method with iterative gain 0.2, and can only decrease to a minimum of 10% within 10 iterations with iterative gain 0.7 when the estimated FRF error is relatively small.When the estimated FRF error is large, the system is not stable Combining the experiment results in this and the above sections, it can be seen that the proposed adaptive optimization iteration learning control method based on the complex optimization algorithm and complex steepest descent algorithm makes the control system convergent no matter if the system FRF model estimation error is small or large.The proposed control method has a stronger robustness to system FRF model estimation error, and it has a faster convergence rate when applying the same system FRF estimation model compared with the traditional iteration learning control method, which indicates that the proposed adaptive optimization iteration learning control method is effective and feasible, and can accurately realize the acceleration time domain waveform replication of shaking tables systems.

Conclusions
In this paper, the adaptive iterative learning control method based on the complex optimization algorithm was proposed to overcome the problem of the traditional iteration learning control method.The complex Broyden optimization algorithm combined with complex steepest descent gain is deduced according to the real Broyden optimization algorithm and the mapping relationship between real vector space and complex vector space.The shaking tables experiment system is built based on xPC target rapid prototyping technology.The tests performed use the proposed adaptive iterative learning control and traditional iterative learning control, respectively, under two conditions: with relatively accurate system FRF estimation and with a large system FRF estimation error.The percent RMSE of the reference signal (RMS is about 0.2 g) can decrease to 5% (about 0.01 g) within three iterations using the proposed control method and at least 20 iterations using the traditional method with iterative gain 0.2, and can only decrease to a minimum of 10% within 10 iterations with iterative gain 0.7 when the estimated FRF error is relatively small.When the estimated FRF error is large, the system is not stable and convergent using the traditional method, even with an iteration gain of 0.1, and the RMSE can still decrease to 15% (0.023 g) within five iterations and the control error will keep decreasing with the increase in iteration time.The experiment results indicate that the proposed control method makes the control system convergent no matter if the system FRF estimation error is small or large, and it has a faster convergence rate and stronger robust stability.

Figure 2 .
Figure 2. Offline iterative learning control schematic for shaking tables.

Figure 2 .
Figure 2. Offline iterative learning control schematic for shaking tables.
is the self-spectrum matrix of the input signal and ( ) matrix of the input signal and output signal.The process of the target simulation stage is as follows:

Figure 2 .
Figure 2. Offline iterative learning control schematic for shaking tables.

Figure 3 .
Figure 3. Offline iterative control schematic considering model error.

Figure 4 .
Figure 4. Adaptive iteration learning control process diagram for shaking tables.

Figure 4 .
Figure 4. Adaptive iteration learning control process diagram for shaking tables.

Figure 5 .Table 1 .
Figure 5. Shaking tables experiment system.(a) Real-time control system, (b) Built experiment system.sTable 1.Main parameters of the shaking tables system.Parameters Values/Units mass of platform 4000/kg mass of load 13,000/kg effective area of cylinder 0.0075/m 2 density of hydraulic oil 845/kg/m 3 nature frequency of servo valve 120/Hz rated flow of servo valve 400 L/min pressure of hydraulic source 25 × 106/Pa

Figure 7 .
Figure 7. Amplitude frequency characteristics of shaking tables.

Figure 6 .
Figure 6.Time-waveform curves of six reference signals.

Figure 7 .
Figure 7. Amplitude frequency characteristics of shaking tables.

Figure 7 .
Figure 7. Amplitude frequency characteristics of shaking tables.

Figure 7 .
Figure 7. Amplitude frequency characteristics of shaking tables.

Figure 8 .
Figure 8. Phase frequency characteristics of shaking tables.

Figure 8 .
Figure 8. Phase frequency characteristics of shaking tables.

Figure 9 .
Figure 9. Errors in 300th frequency under traditional iterative control with gain 0.2.(a) Real part, (b) Imaginary part.

Figure 10 .
Figure 10.Errors under traditional iterative control with gain 0.2 after 20 iterations.

Figure 10 .
Figure 10.Errors under traditional iterative control with gain 0.2 after 20 iterations.

Figure 10 .
Figure 10.Errors under traditional iterative control with gain 0.2 after 20 iterations.
 , respectively, describe the root mean square and mean value of { }  .

Figure 13 .
Figure 13.Errors under traditional iterative control with gain 0.7 after 20 iterations.

Figure 17 .
Figure 17.Impedance matrix update in 300th frequency under adaptive iterative control.(a) Amplitude of impedance matrix, (b) Phase of impedance matrix.

Figure 16 .
Figure 16.Optimal adaptive gain in 300th frequency under adaptive iterative control.(a) Real part, (b) Imaginary part.

Figure 17 .
Figure 17.Impedance matrix update in 300th frequency under adaptive iterative control.(a) Amplitude of impedance matrix, (b) Phase of impedance matrix.

Figure 20 .
Figure 20.Amplitude frequency characteristics of shaking tables.

Figure 20 .
Figure 20.Amplitude frequency characteristics of shaking tables.

Figure 20 .
Figure 20.Amplitude frequency characteristics of shaking tables.

Figure 21 .
Figure 21.Phase frequency characteristics of shaking tables.

Figure 22 .
Figure 22.Errors in 300th frequency under traditional iterative control with gain 0.1.(a (b) Imaginary part.

Figure 21 .
Figure 21.Phase frequency characteristics of shaking tables.

Figure 20 .
Figure 20.Amplitude frequency characteristics of shaking tables.

Figure 21 .
Figure 21.Phase frequency characteristics of shaking tables.

Figure 22 .
Figure 22.Errors in 300th frequency under traditional iterative control with gain 0.1.(a) Real part, (b) Imaginary part.

Figure 23 .
Figure 23.Errors under traditional iterative control with gain 0.1 after 20 iterations.

Figure 26 .
Figure 26.Optimal adaptive gain in 300th frequency under adaptive iterative control.(a) Real part, (b) Imaginary part.

Figure 27 .
Figure 27.Impedance matrix update in 300th frequency under adaptive iterative control.(a) Amplitude of impedance matrix, (b) Phase of impedance matrix.