Data-Driven Nonlinear Iterative Inversion Suspension Control

: The commercial operation of the maglev train has strict requirements for the reliability and safety of the suspension control system. However, due to a large number of unmodeled dynamics of the suspension system, it is difﬁcult to obtain the precise mathematical model of the suspension system. After the suspension system has been operated for a long time with high load, the system model will change due to the wear, aging and failure of components, as well as the settlement of the line and track. The control performance is degraded. Therefore, this paper proposes a data-driven nonlinear iterative inversion suspension control algorithm, which can achieve high-precision tracking performance recovery control after control performance degradation without depending on the suspension system model. The control performance of the suspension system is improved by learning the measured data of the historical suspension system, and the fast convergence of the tracking error and high-precision stable suspension control are realized in the presence of unmodeled dynamics and external noise interference. Based on the historical suspension data of the maglev train suspension control system, the inverse dynamics model of the suspension system is identiﬁed by iterative inversion learning based on data drive, and the suspension control framework based on iterative inversion is designed. Then, the nonlinear input update strategy is used to realize the rapid convergence of the learning process. Finally, the simulation experiment of the maglev train suspension system and the physical experiment of the maglev system experimental platform are combined. It is veriﬁed that the proposed levitation control algorithm can achieve high-precision fast tracking performance recovery control after the system control performance degrades under noise environment.


Introduction
As a new rail transit tool with mature technology, quiet operation, comfortable ride and high degree of localization, maglev train has become a popular means of transportation for citizens.Commercial operation lines, including Beijing Maglev Line S1, Changsha Maglev Express Line, Shanghai High speed Maglev Demonstration Line and Phoenix Culture Tourism Maglev Line, have been well integrated into people's daily travel.In order to meet the strict requirements of the commercial operation reliability and operation performance of the maglev train, it is necessary to ensure that the suspension control system can always maintain a high performance working state.However, when the maglev train on the commercial operation line runs under a high load for a long time, the model parameters of the suspension system will change due to the wear, aging and failure of system components, as well as the track settlement of the line, which will lead to the deterioration of the control performance of the suspension control system.As the allowable variation range of the suspension gap is small, the maglev train will be prone to hit the track.This will affect the passenger's riding experience.In serious cases, it may even affect the normal operation of the maglev train and endanger the safety of passengers.Therefore, in order to ensure the stable operation of the maglev train at the set suspension working gap of 10 mm, it is necessary to design a suspension control algorithm with high tracking performance.
The iterative learning control (ILC) method, which is used to implement repetitive control by learning from historical data to achieve high tracking performance, is widely popular for such problems as high-performance tracking control of suspension system clearance setting [1].In this case, in order to improve the tracking performance of the suspension control, it is necessary to conduct accurate inverse dynamic modeling of the suspension system explicitly or implicitly.Accurate inverse dynamic model is very important to the performance of tracking control.The methods of constructing inverse models include model-based method and modeless method.Model-based iterative learning control uses the inverse model of the system as a learning filter for learning.If an accurate system model can be identified and the learning filter approximates the inverse dynamic model of the controlled object, the learning process can achieve rapid convergence [2,3].Model-based methods usually use stable approximation [4,5] or numerical optimization [6,7] to obtain the inverse model.
However, model-based inversion of ILC requires sophisticated modeling, which is not practical in actual industrial sites.Zhongsheng Hou et al. also pointed out that due to the development of information science and technology, the complexity of industrial systems has increased, and it is difficult to model based on the first principle or identification method [8].Therefore, it is very important to directly design controllers based on a large amount of data generated in industrial processes, which also promotes the development of data-driven control methods.Ronghu Chi et al. proposed a data-driven optimal terminal ILC method for linear and nonlinear discrete-time systems for fixed-input constant-value control systems [9,10].The controller design only depends on the input and output measurement data of the system.Bolder et al. designed a data-driven steepest descent ILC method for the practical application of printing system control [11].For these data-driven learning methods, the system model is constructed by model identification in the iterative learning process.Therefore, a modeless ILC method based on input-output data to identify the inverse dynamic model of the system is increasingly welcomed by researchers [12].
In addition, there are also ILCs based on compression mapping methods, such as D-type ILC, P-type ILC, PD-type ILC and PID-type ILC, and modeless methods based on time reversal ILC.The D-type ILC uses the derivative of the previous error signal to update the input, introducing high-frequency noise.Type P ILC avoids this problem.PD-type ILC can effectively ensure the convergence of tracking error.However, the tuning of PD parameters does not always meet the monotonic convergence conditions of the system, and the tuning process of PD parameters is particularly time consuming and may damage the system.Moreover, such methods based on contractive mapping require that the nonlinear terms of the system meet the global Lipschitz condition and have the same initial conditions.More significantly, such methods can only ensure that the system tracking error converges in the sense of the norm, and the weight of the norm is an exponential that decreases with time, so the tracking error will certainly decrease rapidly with time, which cannot accurately describe the dynamic characteristics of the system on the time axis.The transient problem caused by overshoot in the sense of a 2-norm cannot be avoided, and the fixed learning law form cannot adapt to the degradation of the system well.The ILC based on time reversal uses the adjoint operator of the system as the learning filter, and the iterative learning speed is slow.
Therefore, in order to achieve rapid error convergence, Kim et al. proposed a model free of iterative inversion-based control (MFIIC) method to estimate the inverse frequency response function (FRF) model based on the historical input and output data of the system [13].It is unnecessary to determine the FRF model in advance, reducing the requirements for the model.However, when the system reference value is equivalent to the output interference, the performance will be significantly reduced.Pieter Janssens et al. proposed a data-driven constrained-norm optimal iterative learning control framework for linear time invariant systems for tracking and point to point motion problems [14].They estimated the impulse response of the system using the input and output measurements of the previous iteration, eliminating time-consuming identification experiments.
The inverse dynamic model identification of the system can be obtained by solving the least squares problem in the time domain [15,16] or by point-by-point division in the frequency domain [17][18][19][20].However, because there is no sustained excitation, these methods are usually sensitive to output disturbances, and the inverse model identification may not converge due to output disturbances.Therefore, ILC needs to balance the effective compensation of repetitive interference with the amplification effect of non-repetitive interference [21,22].The key reason for ILC amplification of non-repetitive interference is that the effectiveness of repetitive interference is largely attributed to the non-causality in the time domain, while the essence of feedback control is causality, and ILC has strict causality in the iterative domain.
Since fast convergence can be obtained when the learning filter approaches the inverse system dynamics [23,24], a typical learning filter design method is based on the inversion of the parametric object model [3], which requires taking measures to ensure that the inverse of the non-minimum phase model is stable.The zero point of the non-minimum phase model means that the causal and stable inverse of the system does not exist, and the stable inverse is anti-causal [25].For inverse causal inversion, the input required by the system can be determined according to the stable inversion process [26].Because these methods are essentially model-based, requirements are put forward for parameter identification and approximate inversion of learning filters.
At present, although the system inverse dynamics estimation based on iterative learning reduces the requirement of prior modeling, the unmeasured interference of the system will still affect the convergence of iterative learning when the system inverse dynamics identification is conducted using historical test data.
Therefore, a data-driven nonlinear iterative inversion control method is designed in this paper to achieve high-precision fast-tracking suspension control of the maglev train so as to ensure that the suspension system of the maglev train can also ensure highperformance suspension tracking control when disturbed by track irregularities.

Maglev Control Problem Description
The main interference of the suspension control system is the track irregularity caused by the settlement of the maglev operation line track and the deformation of the track beam, which directly affects the measured value of the suspension gap.The vertical deflection of the track beam increases with the increase of the weight of the train and the span of the track beam.The track irregularity caused by the settlement of the line foundation or supporting pier will also affect the safe operation of the train.As shown in Figure 1.The lifting of the track outside the curve section of the maglev commercial operation line will also bring the influence of track irregularity to the maglev train suspension system.As shown in Figure 2.  In addition, for the control of maglev train suspension system, the following issues need to be considered: (1) The control input of the suspension system is limited (input saturation problem), mainly including the limitation of the output duty cycle of the suspension chopper and the limitation of the input current of the suspension controller.Working for a long time with an excessive current will cause the suspension electromagnet to overheat, which may affect the safe operation of the maglev train in serious cases.(2) The state of the maglev train suspension system is also limited.According to the structure of the suspension frame of the medium and low speed maglev train, the suspension gap in normal operation needs to be maintained at 8-12 mm.(3) Considering the actual operation of the maglev train, such as lifting, lowering, ramps and curves, the control of the suspension system needs to ensure passenger comfort on the maglev train during lifting and lowering, as well as the safe operation of the maglev train on ramps and curves under the influence of track irregularity.(4) Because the mileage of the maglev operation line is limited, the running time of the maglev train is limited, and the actual control performance of the maglev train cannot be guaranteed by the levitation control algorithm that requires the running time to reach infinity to obtain the gradual convergence effect.(5) The dynamic model of the actual suspension system of the maglev train has highorder nonlinear characteristics, and there are a lot of unmodeled dynamics, which makes it impossible to obtain an accurate model of the suspension system.At the same time, because the components of the suspension system will age or even fail, and the mechanical components will wear or even break after the maglev train runs for a long time, the model of the suspension system will also change.This means that the suspension control of the maglev train cannot simply rely on the accurate system model.
Therefore, aiming at the suspension control problem of maglev train suspension system in the presence of track irregularity interference, this paper makes full use of the characteristics of maglev train's cyclic and repetitive operation and designs a data-driven nonlinear iterative inversion suspension control framework, as shown in Figure 3. Where, R is the reference value of the suspension gap, E i is the tracking error, U i is the control input of the suspension system, D i is the output interference of the suspension system, Y i is the output suspension gap, F is the learning filter, F i is the output signal of the i-th iteration of the learning filter, G(z) is the z-domain form of the transfer function of the SISO discrete LTI suspension system that is gradually stable under zero initial conditions, and i is the number of iterations.
Assume that the time domain signal length of the system is infinite to meet the requirements of frequency domain analysis and avoid leakage in frequency domain calculation, that is N → ∞ .The entire impulse response is included in the object model [27].
where k is the time index of the discrete time signal, g(k) ∈ R is the impulse response of object G, and u(k) ∈ R and y(k) ∈ R are the input and output signals of the object respectively.Output signal y(k) is obtained by convolving g(k) and u(k) in the time domain y(k) = G(z)u(k).Therefore, the time domain update law of the general ILC algorithm is where r(k) is the reference signal of the system, e i (k) is the error signal of the i-th iteration of the system, and u i (k) and y i (k) are the input and output signals of the i-th iteration of the system respectively.The general ILC algorithm will reset the initial condition u 0 (k) = F(z)r(k) before each iterative learning.Therefore, it can be obtained that the system tracking error is Iterative tracking control needs to ensure that the tracking error converges to zero after multiple iterations.The goal is to minimize tracking error by updating feedforward using measurement data generated by previous control tasks.That is, for system G(z), the ILC stability condition in frequency domain shown in Equation (4) needs to be satisfied.
where • ∞ denotes H∞-norm.Therefore, the iterative tracking control problem needs to meet the following assumptions [27]: is the DFT form of the output suspension gap, and the DFT form of the tracking error is where the output of the suspension system can be obtained by measurement, and its expression is where D i (ω k ) is the DFT form of bounded track irregularity interference, and δ(ω k ) is the boundary, i.e., ∀i, so that This assumption is to achieve robustness to the worst-case estimation error.For the general ILC algorithm, the error signal is minimized by updating the feedforward with the measurement data from the previous task, thus improving the performance of repetitive control.For example, in the norm optimal ILC algorithm, suppose r i+1 (k) = r i (k) and the tracking error of the i-th iteration is [28]: where S(z) = (I + G(z)C(z)) −1 is the sensitivity function, and f i (k) is the output signal of the learning filter of the i-th iteration.The system minimizes the tracking error based on historical measurement data.Then the tracking error of the i + 1 iteration is where the feedforward signal of the i + 1 iteration can be determined as V( f i+1 (k)) is a data-based performance standard, W e > 0 is a positive definite weight matrix, and W f ≥ 0 and W ∆ f ≥ 0 are semi positive definite weight matrices, which are used to balance the robustness and convergence of iterative learning algorithms.However, due to the influence of track irregularity, the reference signal r(k) of the suspension control system will change, making the key assumptions of the ILC algorithm no longer meet, which will lead to the decline of the tracking performance of the ILC algorithm.

Model Free Nonlinear Iterative Inverse Learning Control for Suspension System
Considering the existence of track irregularity in levitation control, a data-driven iterative inversion learning control update law is determined as [29] where ρ ∈ R is the gain learning function, ρ ≥ 0 and i ≥ 1.The basic principle is to use ρ(|Y i |) to avoid unbounded amplification of input.In addition, if the system is disturbed by changes, the tracking error will not converge to 0. Define the tracking error change rate as

Tracking Performance Conditions Based on Nonlinear IIL
Consider the update law given in Equation (12).If for α > 1, the following constraints are satisfied [29]: where E 0 is the initial tracking error.Then we can obtain When ∆δ = 0 and sequence {E i } converge monotonically to E ∞ = 0. Therefore, when the following conditions are met, the error asymptotic convergence can be guaranteed: (1) The learning gain ρ(|Y i |) has a non-zero strict lower bound, which does not exceed 1; (2) |R| is much larger than δ and ∆δ; (3) The initial tracking error is far less than |R| and δ; and (4) The tracking error converges to the neighborhood of 0, which is proportional to ∆δ.This conclusion can be proved by the following lemma [29]: i=0 satisfies the following boundary conditions, where k ∈ R, k ≥ 0 and ∀i satisfies k i ≤ k. and where k i is given by Equation (20).
In order to obtain this conclusion, suppose that Equation ( 21) is true for i, according to Lemma 2, k i ≤ 1/α, combined with Equation ( 15 Therefore, if Equation ( 21) is true for i, then it is also true for i + 1.Since Equation ( 21) is also true when i = 0, it is proven that Equation ( 21) is true for ∀i [29].
According to Lemma 2, ∀i makes This also proves Equation (17).In addition, if ∆δ = 0, we can directly draw the conclusion that {E i } ∞ i=0 converges monotonically to 0 from Equation ( 18) in Lemma 1.In addition, for the iterative inversion learning control update law of Equation ( 18), the input can converge to the neighborhood of the limit point.
According to the following Theorem 1, if the error converges to the neighborhood of 0, the input will converge to the neighborhood of U ∞ [29]: Theorem 1.Consider Equation ( 12) under the condition of Assumption 1.If Equations ( 14)-( 16) hold for α > 1 then In addition, when the reference is relatively small relative to the interference, Theorem 2 can be obtained [29].
According to Theorem 2, if |R| is smaller than |D i |, then U i will converge to 0. This is obviously different from the iterative learning method based on the update law of the inverse model.In this method, if the convergence condition is not satisfied, the unstable iterative dynamics may lead to unbounded input.However, in NLIIC, the iterative update will not result in unbounded input, and the input will converge to 0. Theorem 2 can be proven by Lemmas 3 and 4 [29]: Lemma 3.Under the condition of Assumption 1, consider Equation ( 12), and each sequence that is, each sequence {U i } ∞ i=0 meets the following restrictions where ∀i, k ∈ R and k ≥ 0, meet k i ≤ k.
The proof of Theorem 2 is as follows: Given 0 < ρ(|Y i |) ≤ 1 for Y i = 0 and for β > 1, ∀i let β|R| ≤ |D i | and |U 0 | < (β − 1) G −1 R hold according to Lemma 4. For A, it is equivalent to B. Because of C, according to D, E can be obtained.Then according to Lemma 4, we can know k 1 < 1. Repeat the above process to obtain ∀i, k i < 1.Therefore, Theorem 2 proves thatTherefore, Theorem 2 proves that {U i } ∞ i=0 converges to U ∞ = 0.This section mainly analyzes the system input convergence of the control algorithm.In addition, according to the reference [20], the conditions that need to be satisfied to obtain the bounded input amplification ρ(|Y i |) can be derived, so as to ensure the global input convergence of the system.According to Equation (27), where τ sup Therefore, ∀x ∈ (0, ∞), if and only if f (x) ρ(x)/x is bounded, the third term on the right side of Equation ( 28) is bounded.Under the assumption of 1, if ρ(x) is continuous in the neighborhood of 0, and lim x→0 ρ(x) = 0 and (dρ/d|Y i |) |Y i |=0 ∈ R exist, then ∀x ∈ (0, ∞), the real function f (x) ρ(x)/x is bounded, and |R| is bounded.
Therefore, under the assumption of 1, if ρ(|Y i |) is continuous in the neighborhood of 0, lim This conclusion shows that if ρ(|Y i |) is continuous near 0, when |Y| → 0 , ρ(|Y i |) → 0 , and ρ(|Y i |) is differentiable near 0, then k satisfies an upper bound.Therefore, for the iterative inverse learning control update law designed in this paper, the amplification of the system input is bounded.This conclusion also provides a global guarantee for the convergence of the system input [29].
In this section, the input amplification is adjusted by using the learning gain function ρ(|Y i |).If ρ(|Y i |) and |R| meet specific attributes, NLIIC can solve the problem of data- driven iterative tracking.
In the next section, we will study the limitations of traditional MFIIC to help design specific gain functions, so as to design and complete the data-driven iterative inversion levitation control algorithm.

Data-Driven Nonlinear Iterative Inversion Suspension Control Algorithm
In this section, the IIC algorithm and the limitations of MFIIC algorithm are reviewed to help design a data-driven nonlinear iterative inversion levitation control algorithm.
According to the literature [13], when the modeling error is not too large and the iterative coefficient is selected, IIC algorithm can achieve accurate tracking, but the algorithm is sensitive to the phase uncertainty of system dynamics, especially when the phase uncertainty is close to π/2.
MFIIC updates the inverse model by using input and output measurement data in each iteration, and makes the iteration coefficient 1, which not only eliminates IIC's modeling process and its related model quality and modeling time consuming problems but also improves the quality of the inverse model through iteration [13].Considering the case of NLIIC ρ(|Y i |) = 1, NLIIC will be simplified to the traditional MFIIC.Obviously, when Y i → 0 , lim the boundedness of input amplification cannot be guaranteed.
According to MFIIC [30], , ∀i ≥ 1 (30) For random noise and non-random effects (such as nonlinear dynamics not modeled), as long as D i /R is small enough, the iterative process will make the noise correlation terms P i and S i tend to zero.For ∀i, when ∃ε makes |D i /R| ≤ ε < 0.5, the relative tracking error is bounded, This ensures MFIIC's ability to improve output tracking.Therefore, MFIIC can improve the tracking performance and input convergence performance of the system.Therefore, the gain function designed in this section is [29] when γ > 0, |Y i | > γ, ρ|Y i | = 1, and γ = |R| are often used.Therefore, when |R| > δ and MFIIC is used for reference, the tracking performance of the system increases monotonously, and the system input is strictly reduced during iteration.
When |Y i | ≤ γ is selected, NLIIC is used for reference to ensure that the amplification of system input is bounded.
where µ ∈ (0, π], when µ ≈ 2.33, Therefore, the input is bounded. In addition, in order to improve the robustness of the control system, especially considering the influence of high-frequency model uncertainty and to achieve a low tracking error in the steady state of the system, the adjoint-based learning method is used for reference.Design the initial switching learning filter F 0 (z) = αG * (z), and make the gain α > 0 small enough to ensure the stability condition, Among them, αG * (z)e i (k) can be obtained based on modeless time inversion filtering technology [28], as shown in Algorithm 1.

S1
Reverse the error signal e1 i (k) = e i (N − k).

S2
Apply the reverse error signal e1 i (k) to the system to obtain the output signal:e2 i (k) = G(z)e1 i (k).

S3
Then reverse the output signal and multiply by a factor small enough:

S4
Finally, the learning law based on adjoint is: However, a small gain α will result in slow convergence.Therefore, this paper combines the advantages of MFIIC's fast convergence, low steady-state tracking error based on time inversion filtering, and the advantages of NLIIC's high tracking performance in the presence of interference and designs a data-driven nonlinear iterative inversion levitation control algorithm.The computational cost of the algorithm inherits the advantages of MFIIC.Its time complexity is O(n) and space complexity is O(1).Therefore, the iterative learning speed of the algorithm is fast, as shown in Figure 4.The algorithm flow is shown in Algorithm 2.

S2
Apply input u i (k) to the system to obtain output y i (k), and then calculate the tracking error e i (k) = r(k) − y i (k) of the system.If the tracking error is small enough and within the acceptable range, stop the iterative learning process or go to the next step.

S3
If i = 0, go to step 4, otherwise go to step 5.

S4
Execute the ILC based on time inversion, u i+1 (k) = u i (k) + αG * (z)e i (k), and then let i = i + 1. Go to step 2 for execution.

S5
Discrete Fourier transform is used to obtain frequency domain signals

S6
Then combine learning gain function ρ(|Y i |) and obtain U i+1 (ω k ) according to Equation (11).Finally, perform IDFT to obtain time domain u i+1 (k), and then turn to step 2 for i = i + 1 to execute.
In the iterative learning process, the inverse model is applied to the learning filter of the nonlinear iterative inversion levitation control algorithm to accelerate the convergence of the tracking error [27,31].The learning filter update structure is shown in Figure 5. Once the tracking error converges, the inverse model-based optimal learning filter can be used to track any trajectory.The learning filter update algorithm flow based on the inverse model is shown in Algorithm 3. Finally, the updated learning filter is replaced by the inverse model.

Suspension Control Experiment and Analysis
This section takes Changsha medium-and low-speed maglev trains as the research object, combines a large number of maglev train suspension gap data to realize the identification of the inverse model of the suspension system, and selects historical data with a high signal-to-noise ratio to improve the estimation accuracy.
The suspension clearance data of the maglev train is obtained through the data acquisition system.The data acquisition system is mainly composed of suspension gap sensor, acceleration sensor, current sensor, temperature sensor and corresponding data receiving and processing module.It can collect suspension gap, vertical acceleration, lateral acceleration, electromagnet working current and electromagnet coil surface temperature data during the operation of the maglev train, with a sampling rate of 1 kHz, as shown in Figure 6.The prior information of the suspension system is used to improve the convergence speed of the data-driven nonlinear iterative inversion suspension control algorithm.The suspension system of maglev train after mechanical decoupling of suspension bogie is simplified as a single point suspension system, as shown in Figure 7. Considering the single point suspension system of the maglev train, after linear processing at the suspension balance point, the state variable is ∆z(t) ∆ .z(t) , and the state space description of the suspension system model can be obtained as follows [32]: where ∆i(t) is the change of coil current; ∆z(t) is the change of air gap; ∆u(t) is the change of the voltage at both ends of the solenoid coil.According to the data of Changsha maglev trains, vacuum permeability µ 0 = 12.6 × 10 −7 H/m, coil current i 0 = 22.0A in steady state, time gap value in steady state z 0 = 8.0 mm, coil turns N = 360 turns, electromagnet sectional area S = 0.038 m 2 , coil resistance R = 0.92 Ω, and electromagnet equivalent mass m = 535 kg.Therefore, the state space of the suspension system model describes the system matrix A = 0 1 5467.520 , control matrix B = 0 −1.99 , and output matrix C = 1 0 .
Then, a controller based on PID algorithm is designed to make the closed-loop stable.On the basis of this closed-loop single point suspension system model, considering the influence of track irregularity, a simulation experiment of nonlinear iterative inversion suspension control based on data drive is carried out [33,34].After the suspension system simulation experiment, the iterative learning error curve is obtained, as shown in Figure 8.The system converges after finite iterative learning.It can be observed from the simulation experiment result curve that the error convergence speed of ILC based on time inversion is the slowest under the output noise interference environment with standard deviation of 1 mm and mean value of 0 mm.Compared with the ILC method based on time inversion, Accelerated Convergence Interleaving further improves the learning rate.The error convergence rate of MFIIC method is the fastest, but the method may have unpredictable learning transient for the case with output interference, as shown in Figure 9. NLIIC solves the robustness problem of MFIIC in the presence of output interference and has a fast error convergence rate.However, the steady-state tracking performance of NLIIC is worse than that of MFIIC and the methods proposed in this paper.As shown in Figure 10, at 2.5 s of the steady-state tracking error curve, NLIIC has a larger tracking error than other algorithms.The NLIIC method based on time inversion proposed in this paper has good steady-state tracking performance, and the iterative learning error convergence rate is fast, close to the error convergence rate of NLIIC.The time domain tracking error after iterative learning is shown in Figure 10.According to Figure 11, when the number of updates of the learning filter increases, the tracking error of the system pulse response decreases.When the seventh learning filter update is completed, it can be considered that the FIR inverse filter at this time is superior to the previous results.The FIR inverse model is used to replace the learning filter, which achieves a faster error convergence rate [31].Then the tracking control problem of the suspension system in the presence of track irregularity interference is studied.A step signal with Gaussian white noise interference is applied at the output position of the suspension system.After the suspension system simulation experiment, the step response simulation curve of the suspension system is obtained, as shown in Figure 12.It can be seen from Figure 12 that the suspension system can realize fast tracking control, and the tracking error of the system is within ±0.05 mm, as shown in Figure 13.The experimental results show that the data-driven nonlinear iterative inversion suspension control based on time inversion can achieve high-precision tracking control.Finally, the tracking control experiment of the levitation system is carried out on the experimental platform of the magnetic levitation system.Verify the tracking control performance of the control algorithm designed in this paper in the presence of track irregularity interference.The control circuit module of the magnetic suspension vibration experimental platform is shown in Figure 14.The magnetic suspension vibration experiment platform includes the baseline series real-time target machine of Speedgoat Company (Bern, Switzerland), which realizes the real-time operation control and monitoring of the suspension system.The levitation electromagnet module and its corresponding control circuit board, chopper circuit board, upper computer and power supply module can simulate the levitation system of the maglev train.It also includes a vibrating table base module capable of applying external interference excitation.The suspension electromagnet module of the magnetic suspension system experimental platform is shown in Figure 15.For the tracking control experiment of the suspension system, first of all, build the most commonly used closed-loop suspension control structure based on PID in engineering, adjust the PID parameters to make the system suspend stably, and obtain the experimental results of the step response output suspension gap curve of the closed-loop suspension system experiment based on PID algorithm without feedforward filter link, as shown in Figure 16.The control error curve of the step response experiment of the corresponding closed-loop suspension system is shown in Figure 17.At this time, the suspension control system can better follow the step input reference signal of the system but has a large overshoot.Then, on the basis of the closed-loop suspension control structure based on PID, the feedforward filter link of the data-driven nonlinear iterative inversion suspension control based on time inversion is added, and the suspension gap curve of the step response output of the data-driven nonlinear iterative inversion suspension control experiment with the feedforward filter link is obtained as shown in Figure 18.It can be seen from Figure 18 that after the results are updated iteratively, the system can track the step input reference signal well and has a small overshoot.The experimental results verify the effectiveness of the proposed method.

Figure 1 .
Figure 1.Vertical section of maglev commercial operation line project.

Figure 2 .
Figure 2. Cross section diagram of track beam structure in curved section of maglev commercial operation line.

Figure 3 .
Figure 3. Suspension control block diagram of nonlinear IIC based on data-driven.

Figure 4 .
Figure 4. Schematic diagram of gain switching structure of control algorithm learning law proposed in this paper.

Figure 5 .
Figure 5. Schematic diagram of control algorithm learning filter updating structure proposed in this paper.Algorithm 3: Update learning filter S1 Initialize the reference pulse signal to r(k) = αδ(k).Execute the algorithm initialization in Algorithm 2. S2 Perform the algorithm iteration update steps in Algorithm 2. S3 Construct the inverse model as G(z) = ∑ N/2 k=−N/2 u(k)z −k .S4Finally, the updated learning filter is replaced by the inverse model.

Figure 6 .
Figure 6.Schematic diagram of maglev data system.

Figure 7 .
Figure 7. Schematic diagram of single point suspension system of maglev train.

Figure 8 .
Figure 8. Error curve of iterative learning process.

Figure 9 .
Figure 9. Error curve of iterative learning process.

Figure 11 .
Figure 11.Pulse response tracking error curve after learning filter updating.

Figure 12 .
Figure 12.Simulation curve of suspension system step response.

Figure 13 .
Figure 13.Suspension system tracking control error curve.

Figure 14 .
Figure 14.Schematic diagram of control circuit module of magnetic suspension system experimental platform.

Figure 15 .
Figure 15.Schematic diagram of suspension electromagnet module of magnetic suspension system experimental platform.

Figure 16 .
Figure 16.Suspension clearance curve of step response output of closed-loop suspension system experiment based on PID algorithm without feedforward filter link.

Figure 17 .
Figure 17.Control error curve of step response of closed-loop suspension system experiment based on PID algorithm without feedforward filter link.

Figure 18 .
Figure 18.Suspension gap curve of step response output of suspension control experiment based on time driven nonlinear iterative inversion with feedforward filter link.