Dynamic High-Type Interval Type-2 Fuzzy Logic Control for Photoelectric Tracking System

: This paper proposes a dynamic high-type control (DHTC) method based on an interval type-2 fuzzy logic controller (IT2FLC), which is used in the photoelectric tracking system to improve the steady-state accuracy and response speed. Adding integrators to the traditional multi-loop feedback control loop can increase the system type, thereby speeding up the response speed and improving the steady-state accuracy, but there is a risk of integral saturation. Switching the type dynamically according to the system state can avoid integral saturation while retaining the advantages of the high-type. Fuzzy logic control (FLC) can dynamically change the output value according to the input change and has the advantages of fast response speed and strong ability to handle uncertainties. Therefore, in this paper, the FLC is introduced into the high-type control system, and the output of the FLC is used as the gain of the integrator to control the on-off to achieve the goal of dynamic switching type, which is successfully veriﬁed in the experiment. IT2FLC introduces a three-dimensional membership function, which further improves the FLC’s ability to handle uncertainties. From the experimental results, compared with T1FLC, IT2FLC’s ability to handle uncertainties is signiﬁcantly improved. In addition, in order to speed up the calculation speed of IT2FLC, this paper proposes an improved type-reduction algorithm, which is called weighted-trapezoidal Nie-Tan (WTNT). Compared with the traditional type-reduction algorithm, WTNT has faster calculation speed and better steady-state accuracy, and has been successfully applied to real-time control systems, which has good engineering application value. Finally, in order to reduce the interference of human factors and improve the automation level of the system, a multi-population genetic algorithm (MPGA) is used to iteratively optimize the parameters of the FLC, which improves the output accuracy. On the experimental platform of the ﬂexible fast steering mirror (FFSM), the control effects of the traditional controller, T1FLC and IT2FLC are compared, which proves that the IT2FLC-DHTC system has a faster response performance, higher steady-state accuracy, and stronger ability to handle uncertainties.


Introduction
Photoelectric servo tracking equipment is widely used in beam control systems, such as adaptive optics, free-space communication, line of sight stabilization, etc., and has important application value in military and civilian fields [1][2][3]. Its tracking performance is the most important index to evaluate the control system. Commonly used tracking control methods include PID [4,5], speed/acceleration lag compensation [6,7], and disturbance observer based on a multi-loop feedback control [8]. Among them, the PID method is simple in design, stable in performance, and the most widely used. However, due to its 1.
A DHT-IT2FLC system is constructed to realize dynamic high-type control. As a judging mechanism, the FLC changes the system type in real time according to the system status, which not only achieves the purpose of improving the steady-state accuracy and speeding up the response speed, but also avoids the hidden dangers of system oscillation and integral saturation caused by the increase of type.

2.
Compared with the traditional controller and the T1FLC, the IT2FLC can significantly improve the ability of the system to handle uncertainty, make the system run more smoothly and stably, and enhance the robustness of the system. 3.
Aiming at the TR problem of T2FLC in real-time system applications, the improved WTNT TR algorithm is proposed to reduce computational complexity. By comparison, WTNT has faster response speed and calculation accuracy than EKM and EIASC in this project.
The structure of this paper is as follows. Section 2 systematically introduces the background knowledge of DHTC, FLC, and MPGA; Section 3 introduces the design of the control system; Section 4 presents the experimental results, and make conclusions in Section 5.

Problem Formulation and Preliminaries
This chapter systematically introduces DHTC, FLC, and MPGA to prepare the basics for controller design.

DHTC
The commonly used photoelectric tracking control method is the classic type-I double closed-loop feedback control, as shown in Figure 1. In the figure, Gv is the mathematical model identified from the FSM experimental platform, and Cv is the speed controller. These two parts form an inner loop called the speed loop; the speed loop and the position controller Cp in the series form an outer loop called the position loop. The open-loop transfer function of the system is G(s) = a·k·(T 3 s + 1) s·(T 4 s + 1)·(a + (T 1 s + 1)·(T 2 s + 1)) = K ∏ m i=1 (τ i s + 1) s ν ∏ n−ν j=1 T j s + 1 (1) where a and k are the gains of the inner loop and the outer loop, respectively, and T 1 , T 2 , T 3 , and T 4 are the time constants of the transfer function of each controller. According to the definition in the principle of automatic control, since there is an integrator in the system, namely υ = 1, it is a type-I system. The systematic error e is defined as e(s) = r(s) − y(s) = (1 − Φ(s))r(s) = E(s)r(s), where r is reference input, y is systematic output, E(s) is the error transfer function, Φ(s) is the closed-loop transfer function of the system Φ(s) = y(s) r(s) = b m s m + b m−1 s m−1 + . . . + b 1 s + b 0 a n s n + a n−1 s n−1 + . . . + a 1 s + a 0 The step signal, velocity signal, and acceleration signal existing in the physical system are expressed as ( ) = , ( ) = , ( ) = . For the superimposed input of these three signals, e(∞) will remain in a certain range only when ≥ 2; that is, the system can reach a stable state only if the system type is greater than 2.
To improve the steady-state accuracy of the system, multiple integrators are added to the system to change the system into a high-type state. However, the increase of the Therefore, e(s) can be expressed as e(s) = a n s n + . . . + c m s m + c m−1 s m−1 + . . . + c 1 s + c 0 a n s n + a n−1 s n−1 + . . . + a 1 s + a 0 (4) where c i = (a i − b i ), i = 0, . . . , m. According to the final value theorem, when the system is in a stable state, the steady-state error e(∞) can be expressed as e(∞) = lim s→0 s a n s n + . . . + c m s m . . . + c 1 s + c 0 a n s n + a n−1 s n−1 + . . . + a 1 s + a 0 r(s) = lim s→0 s ν+1 r(s) The step signal, velocity signal, and acceleration signal existing in the physical system are expressed as r(s) = R 1 s , r(s) = R 2 s 2 , r(s) = R 3 s 3 . For the superimposed input of these three signals, e(∞) will remain in a certain range only when υ ≥ 2; that is, the system can reach a stable state only if the system type is greater than 2.
To improve the steady-state accuracy of the system, multiple integrators are added to the system to change the system into a high-type state. However, the increase of the system type will increase the step response overshoot and increase the system oscillation; as the gain of the integrator becomes larger, the oscillation will be greater, and in some extreme cases it will even reach the integral saturation state, causing the system to lose stability completely.
To solve the above-mentioned problems, the 'dynamic high-type' method was proposed. As shown in Figure 2, multiple integrators are connected in parallel in the forward path, and the on and off of the integrators are judged according to the system error state, and the type of the system is dynamically changed. However, there are two problems to be solved: The step signal, velocity signal, and acceleration signal existing in the physical system are expressed as ( ) = , ( ) = , ( ) = . For the superimposed input of these three signals, e(∞) will remain in a certain range only when ≥ 2; that is, the system can reach a stable state only if the system type is greater than 2.
To improve the steady-state accuracy of the system, multiple integrators are added to the system to change the system into a high-type state. However, the increase of the system type will increase the step response overshoot and increase the system oscillation; as the gain of the integrator becomes larger, the oscillation will be greater, and in some extreme cases it will even reach the integral saturation state, causing the system to lose stability completely.
To solve the above-mentioned problems, the 'dynamic high-type' method was proposed. As shown in Figure 2, multiple integrators are connected in parallel in the forward path, and the on and off of the integrators are judged according to the system error state, and the type of the system is dynamically changed. However, there are two problems to be solved: Question 1. The standard for type switching is not defined. When to increase the type and when to decrease it, no theoretical guidance is given in the current literature. Besides, the manual switching method used in most of the literature is inaccurate and not practical. Question 2. The jitter caused by type switching is difficult to eliminate. In the existing research, every moment of type switching will bring greater jitter to the system, increase instability, and reduce its robustness.

IT2FLC
The common point of the T1FLS and T2FLS fuzzy systems is that they are based on the fuzzy rule library constructed by expert experience for fuzzy inference, and the input is processed to obtain the output. Generally, the structure of two inputs and one output is commonly used. The two inputs are the system error e and the error rate of change ∆e, and the output is the control signal. The difference between the T1FLS and T2FLS lies in the membership function, defuzzification method, and the ability to deal with uncertainties. Question 1. The standard for type switching is not defined. When to increase the type and when to decrease it, no theoretical guidance is given in the current literature. Besides, the manual switching method used in most of the literature is inaccurate and not practical. Question 2. The jitter caused by type switching is difficult to eliminate. In the existing research, every moment of type switching will bring greater jitter to the system, increase instability, and reduce its robustness.

IT2FLC
The common point of the T1FLS and T2FLS fuzzy systems is that they are based on the fuzzy rule library constructed by expert experience for fuzzy inference, and the input is processed to obtain the output. Generally, the structure of two inputs and one output is commonly used. The two inputs are the system error e and the error rate of change ∆e, and the output is the control signal. The difference between the T1FLS and T2FLS lies in the membership function, defuzzification method, and the ability to deal with uncertainties.

Construct of IT2FLC
Let ∼ A be a type-2 fuzzy set (T2FS), u ∈ [0, 1] is the primary membership degree of ∼ A, and u ∈ [0, 1] is the secondary membership degree of ∼ A. When u is always 1, the Processes 2022, 10, 562 5 of 20 calculation amount of fuzzy inference can be greatly simplified. At this time, T2FLC is called interval T2FLC (IT2FLC). When x ∈ X is input, ∼ A can be expressed by Equation (6): The membership function of IT2FLC is different from that of T1FLC. As shown in Figure 3, the value of the main membership function (MF) on the two-dimensional plane is an area (shaded area), which is contained by LMF (Lower Membership Function) and UMF (Upper Membership Function) and is called FOU (Footprint of Uncertainty). LMF and UMF are represented by µ∼ A (x) and µ∼ A (x), respectively, so ∼ A can also be expressed as The membership function of IT2FLC is different from that of T1FLC. As shown in Figure 3, the value of the main membership function (MF) on the two-dimensional plane is an area (shaded area), which is contained by LMF (Lower Membership Function) and UMF (Upper Membership Function) and is called FOU (Footprint of Uncertainty). LMF and UMF are represented by ( ) and ( ), respectively, so can also be expressed as  The introduction of FOU can include the uncertainties in the input of the system, so IT2FLC has a stronger ability to handle uncertainties. Figure 4 shows the structure of IT2FLC. Compared with T1FLC, there is one more step: type reduction. Its function is to convert the T2FS output by fuzzy inference into T1FS, and then perform defuzzification calculation. Other processes are similar to T1FLC.  The introduction of FOU can include the uncertainties in the input of the system, so IT2FLC has a stronger ability to handle uncertainties. Figure 4 shows the structure of IT2FLC. Compared with T1FLC, there is one more step: type reduction. Its function is to convert the T2FS output by fuzzy inference into T1FS, and then perform defuzzification calculation. Other processes are similar to T1FLC.
The membership function of IT2FLC is different from that of T1FLC. As shown in Figure 3, the value of the main membership function (MF) on the two-dimensional plane is an area (shaded area), which is contained by LMF (Lower Membership Function) and UMF (Upper Membership Function) and is called FOU (Footprint of Uncertainty). LMF and UMF are represented by ( ) and ( ), respectively, so can also be expressed as  The introduction of FOU can include the uncertainties in the input of the system, so IT2FLC has a stronger ability to handle uncertainties. Figure 4 shows the structure of IT2FLC. Compared with T1FLC, there is one more step: type reduction. Its function is to convert the T2FS output by fuzzy inference into T1FS, and then perform defuzzification calculation. Other processes are similar to T1FLC.
The activation intensity of the l rule is calculated by Equation (9): The most commonly used center-of-gravity (COG) method reducer is shown in Equation (10): where y l , y r can be obtained by the TR algorithm, and the final defuzzification obtains the output of IT2FLC as y = y l + y r 2 (11)

The Principle of FLC Realizes DHTC Technology
According to the introduction in Part 2 of this article, FLC can solve the two problems of DHTC.
For Question 1: As shown in Figure 3, the input of FLC is a range. By setting a reasonable membership function, a weight can be set for each precise value in the input range. We take the e and ∆e of the photoelectric tracking system as the input of the FLC, set up reasonable fuzzy rules, and let the FLC determine the timing of switching.
For Question 2: Use the output value of FLC as the gain of the integrator. In the same way, by setting a reasonable membership function, the output value of the fuzzy controller is set to a gradual change, and the traditional "0" and "1" switch structure is changed to an integrator gain gradual change structure, which can avoid the shock caused by hard switching. When the fuzzy output is 0, the integrator is disconnected to realize the dynamic switching type.

WTNT TR Algorithm
Since the input membership function of IT2FLC is type-2, we need to transform T2-FS into T1-FS after the fuzzy inference to get the accurate output. Commonly used TR methods include KM, EKM, IASC, EIASC, etc., based on iterative convergence, which has the advantage of high calculation accuracy, but the disadvantage is slow calculation speed. There is also a kind of direct TR method without iteration, such as NT, BMM, UB algorithm, etc., which has the advantage of fast calculation speed, but in some cases, it is difficult to guarantee the accuracy of the degradation. The WEKM algorithm proposed by Mendel et al. is based on numerical integration to perform the weighted calculation of EKM. After simulation comparison, the calculation accuracy can be improved without affecting the speed [30]; based on Mendel research, CHEN et al. proposed GT2 FLS' WNT TR method; it is proved through simulation that the combination of numerical integration method with NT algorithm can also effectively improve the accuracy of the fuzzy defuzzification [31][32][33]. This provided inspiration to the author and constructed the WTNT reduction algorithm based on the compound trapezoidal rule in IT2FLC and applied it to the physical experimental system for the first time.

NT and CNT
The most recent studies on the type-reduction of interval type-2 fuzzy set have shown that the continuous version of NT (CNT) algorithms is an accurate method for defuzzifing [34]. Intuitively, we can feel that if the sampling size approaches infinity, random sampling will become an accurate TR method. The proof is as follows. Theorem 1. When the number sampling is infinite, the random sampling algorithm will calculate the accurate centroid of IT2FS.

Proof Theorem 1. Suppose there is a discrete IT2FS
∼ A, whose domain is X. Let M represent the number of vertical slices along the y-axis, and L represents the number of horizontal slices along the µ-axis. Therefore, the total number of T1FS embedded in ∼ A is LM, because T1FS contains M vertical slices [35].
Let µ∼ A (x) and µ∼ A (x) denote the LMF and UMF of FOU, respectively. Select N embedded T1FS of ∼ A randomly. Let µ i (x) be the MF of the i-th embedding set. Then, we summarized N randomly embedded T1-FS, Hence, Hence, For Equation (14), the left side is a random sampling of N → ∞ , and the right side is exactly the set of all ∼ A embedded MFs. Note that, since L N is only a constant, it does not affect calculating the centroid of ∼ A.
When ∼ A is a continuous IT2FS form, it is easy to get Equation (15) by replacing integral with summation End of proof. Therefore, when the step size of random sampling satisfies certain conditions, the NT algorithm in the discrete domain approximates the CNT algorithm, which is sufficient to meet the accuracy of the reduction.

Newton-Cotes Quadrature Formulas
Applying the Newton-Cotes formula, the definite integral of f (x) in the interval [a, b] can be estimated by using a limited number of sampling points. This provides a mathematical basis for the WTNT method through systematic sampling.
Then the above equation is called numerical integration or quadrature formula, where x k is the integration node, ω k is the weight coefficient, and E( f ) is the truncation error, that is, the error term of the integration formula.
If the function f is second-order continuous and differentiable on [a, b], then the remainder where a < ζ < b. O h 2 means that when the calculation step size is reduced by 1/2, the error should be reduced to (1/2) 2 = 0.25 times. That means that when the sampling step is smaller, the estimation error is smaller, which is sufficient to meet the accuracy requirements.

WTNT TR Algorithm
The classic NT algorithm was proposed by Nie and Tan [26]. It does not require iterative calculations, such as EM, which greatly simplifies the computational complexity and is more suitable for real-time systems. The NT algorithm is as Equation (21), and the calculated y l , y r are used in Equation (11) to calculate the exact value of fuzzy output.
Through the discussion of Theorems 1 and 2, it can be concluded that the NT algorithm in the discrete domain can meet the accuracy requirements of the defuzzification. Set a = y 1 < y 2 < . . . < y N = b, construct the WTNT TR algorithm where ω j is the weighting coefficient, which is positioned as Based on the principle of numerical integration, the WTNT algorithm can theoretically approximate the CNT algorithm, which is more accurate than the NT algorithm. In fact, the NT algorithm is only a special case of the WTNT algorithm when ω j = 1 (j = 1, 2, . . . , N). Since the output of the IT2FLC set in this paper is discretized, it is not necessary to discuss the reduction of the continuous domain as in the literature [28,29], which simplifies the calculation. The sampling step is determined by fuzzy rules, rather than artificial settings.
Since the number of fuzzy rules is generally more than 9, the estimation error of the WTNT algorithm is very small, and the reduction accuracy can be guaranteed theoretically.

DHT-IT2FLC Design
Based on the above discussion, this paper constructs an IT2FLC-DHTC system as shown in Figure 5. Based on the original double closed-loop control system, a branch of FLC with integrator is connected in parallel. The FLC judges the system error and the rate of error change. The output of FLC is the gain of the integrator. The input and output scaling factors of FLC are iteratively optimized by MPGA.
Processes 2022, 10, x FOR PEER REVIEW 9 of 20 of the WTNT algorithm is very small, and the reduction accuracy can be guaranteed theoretically.

DHT-IT2FLC Design
Based on the above discussion, this paper constructs an IT2FLC-DHTC system as shown in Figure 5. Based on the original double closed-loop control system, a branch of FLC with integrator is connected in parallel. The FLC judges the system error and the rate of error change. The output of FLC is the gain of the integrator. The input and output scaling factors of FLC are iteratively optimized by MPGA.

TI-FLC Design
1. The MF of T1FLC is set as shown in Figure 6. For two inputs e, ∆e, and one output k, the five language variables NB, N, Z, P, PB are set. The NB language value of k is set to a constant of 0 instead of an interval, representing the disconnection of the integrator.
Before setting up the fuzzy rule base, the judgment mechanism of type switching should be analyzed. Through a large number of previous experiments, it is found that when e is large, the access integrator can easily lead to system divergence. When ∆e is large, the system error will change rapidly; at this time, the interruption of integrator should be handled carefully,

1.
The MF of T1FLC is set as shown in Figure 6. For two inputs e, ∆e, and one output k, the five language variables NB, N, Z, P, PB are set. The NB language value of k is set to a constant of 0 instead of an interval, representing the disconnection of the integrator.

2.
The rule base setting.
Processes 2022, 10, x FOR PEER REVIEW 9 of 20 of the WTNT algorithm is very small, and the reduction accuracy can be guaranteed theoretically.

DHT-IT2FLC Design
Based on the above discussion, this paper constructs an IT2FLC-DHTC system as shown in Figure 5. Based on the original double closed-loop control system, a branch of FLC with integrator is connected in parallel. The FLC judges the system error and the rate of error change. The output of FLC is the gain of the integrator. The input and output scaling factors of FLC are iteratively optimized by MPGA.

TI-FLC Design
1. The MF of T1FLC is set as shown in Figure 6. For two inputs e, ∆e, and one output k, the five language variables NB, N, Z, P, PB are set. The NB language value of k is set to a constant of 0 instead of an interval, representing the disconnection of the integrator.
Before setting up the fuzzy rule base, the judgment mechanism of type switching should be analyzed. Through a large number of previous experiments, it is found that when e is large, the access integrator can easily lead to system divergence. When ∆e is large, the system error will change rapidly; at this time, the interruption of integrator should be handled carefully, Before setting up the fuzzy rule base, the judgment mechanism of type switching should be analyzed. Through a large number of previous experiments, it is found that when e is large, the access integrator can easily lead to system divergence. When ∆e is large, the system error will change rapidly; at this time, the interruption of integrator should be handled carefully, otherwise the system divergence will easily occur. When e is near the steady-state range, the access of the integrator should be avoided as much as possible; because there is noise in the actual system, its ∆e is irregular. In this case, it is easy to connect the integrator and cause the system to deviate from the balanced position again. When e is in P and N state, it is the best time to connect the integrator, which can achieve a significant correction effect under the premise of system stability. Finally, the mechanism for determining type switching is shown in Figure 7, when e and ∆e are in six states in the figure: • State 1: e is P, ∆e is PB, and e will continue to increase. At this time, the integrator should be connected to intervene to make the output k value of FLC be PB so that the system can quickly return to the reference value. • State 2: e is PB, ∆e is P. To prevent integral saturation, disconnect the integrator. • State 3: e is PB, ∆e is Z. The system is about to approach the equilibrium position, and the integrator should be connected to speed up the response. • State 4: e is PB, and ∆e is N. The system approaches the equilibrium position at a faster speed. At this time, the integrator can be properly connected, but the gain cannot be too large. • State 5: e is P and ∆e is NB. It is safer to connect to the integrator in this situation. • State 6: e is Z. As mentioned earlier, the integrator should be disconnected to prevent the system from deviating from the equilibrium position. • Similarly, the same judgment mechanism is used when e is N and NB. Therefore, set the fuzzy rules of T1FLC as shown in Table 1.
2022, 10, x FOR PEER REVIEW 10 of 20 otherwise the system divergence will easily occur. When e is near the steady-state range, the access of the integrator should be avoided as much as possible; because there is noise in the actual system, its ∆e is irregular. In this case, it is easy to connect the integrator and cause the system to deviate from the balanced position again. When e is in P and N state, it is the best time to connect the integrator, which can achieve a significant correction effect under the premise of system stability. Finally, the mechanism for determining type switching is shown in Figure 7, when e and ∆e are in six states in the figure: • State 1: e is P, Δe is PB, and e will continue to increase. At this time, the integrator should be connected to intervene to make the output k value of FLC be PB so that the system can quickly return to the reference value. • State 2: e is PB, Δe is P. To prevent integral saturation, disconnect the integrator. • State 3: e is PB, Δe is Z. The system is about to approach the equilibrium position, and the integrator should be connected to speed up the response. • State 4: e is PB, and Δe is N. The system approaches the equilibrium position at a faster speed. At this time, the integrator can be properly connected, but the gain cannot be too large.

•
State 5: e is P and Δe is NB. It is safer to connect to the integrator in this situation. • State 6: e is Z. As mentioned earlier, the integrator should be disconnected to prevent the system from deviating from the equilibrium position.

•
Similarly, the same judgment mechanism is used when e is N and NB. Therefore, set the fuzzy rules of T1FLC as shown in Table 1.

3.
The defuzzification uses the COG method:

IT2FLC Settings
The MF setting is shown in Figure 8. Since the membership function of IT2FLC is two-dimensional, the calculation complexity is relatively high. To speed up the calculation, we reduce the number of language variables. From the experimental comparison, the calculation accuracy can still meet actual needs. Set three language variables N, Z, P for e and ∆e, and each variable has upper and lower bounds. The output variable k follows the form of T1FLC.
Processes 2022, 10, x FOR PEER REVIEW 11 of 20 1 The dark blue font in the table is FLC output k.
3. The defuzzification uses the COG method:

IT2FLC Settings
The MF setting is shown in Figure 8. Since the membership function of IT2FLC is two-dimensional, the calculation complexity is relatively high. To speed up the calculation, we reduce the number of language variables. From the experimental comparison, the calculation accuracy can still meet actual needs. Set three language variables N, Z, P for e and Δe, and each variable has upper and lower bounds. The output variable k follows the form of T1FLC.
The design of the fuzzy rule library follows the principle of T1FLC, but the number of rules is changed from 5 × 5 to 3 × 3, as shown in the dotted box in Table 1.
The reduction algorithm selects WTNT, EKM, and EIASC for comparison.

MPGA Advantages
As Figure 9 shows, the MPGA is constructed based on the coding, selection, crossover, and mutation of the standard genetic algorithm (SGA) which adds an immigration operator and multiple populations. MPGA breaks through the framework of SGA that only relies on a single population for genetic evolution, introducing multiple populations for simultaneous optimization searches. Different populations are assigned different crossover and mutation probabilities, taking into account global and local search performance. The various populations are connected through the migration operator, and its function is to replace the best value of the former group with the worst quality of the latter group to realize the co-evolution of multiple groups. The best individuals in each evolutionary generation of various groups are preserved by manual selection operators to form the elite populations. The elite populations do not perform selection, crossover, mutation, and other genetic operations to ensure that the best individuals generated by various groups are not destroyed or lost during the evolution process. At the same time, the elite populations are also the basis for judging the termination of the algorithm. In MPGA, the minimum algebra of the optimal individual is used as the termination criterion. This criterion makes full use of the knowledge accumulation of the genetic algorithm in the evolution process and is more reasonable than the maximum genetic algebra criterion [36,37]. The design of the fuzzy rule library follows the principle of T1FLC, but the number of rules is changed from 5 × 5 to 3 × 3, as shown in the dotted box in Table 1.
The reduction algorithm selects WTNT, EKM, and EIASC for comparison.

MPGA Advantages
As Figure 9 shows, the MPGA is constructed based on the coding, selection, crossover, and mutation of the standard genetic algorithm (SGA) which adds an immigration operator and multiple populations. MPGA breaks through the framework of SGA that only relies on a single population for genetic evolution, introducing multiple populations for simultaneous optimization searches. Different populations are assigned different crossover and mutation probabilities, taking into account global and local search performance. The various populations are connected through the migration operator, and its function is to replace the best value of the former group with the worst quality of the latter group to realize the co-evolution of multiple groups. The best individuals in each evolutionary generation of various groups are preserved by manual selection operators to form the elite populations. The elite populations do not perform selection, crossover, mutation, and other genetic operations to ensure that the best individuals generated by various groups are not destroyed or lost during the evolution process. At the same time, the elite populations are also the basis for judging the termination of the algorithm. In MPGA, the minimum algebra of the optimal individual is used as the termination criterion. This criterion makes full use of the knowledge accumulation of the genetic algorithm in the evolution process and is more reasonable than the maximum genetic algebra criterion [36,37].

Determine the Optimized Parameters
There are many ways to optimize FLC with GA, which can optimize parameters, such as rule base, membership function, scale factor, and so on. Through experimental comparison, it is found that the optimization rule base or membership function has a large amount of calculation, complex parameter settings, and low efficiency; the optimization scale factor is simple, the calculation speed is fast, and it is convenient to modify. Therefore, this project keeps the rule base and membership function unchanged and uses the three scaling factors Ke, Kec, and Ku ( Figure 5) as the parameters to optimize the chromosome coding.
The larger the Ke value, the faster the step response, but at the same time, the overshoot and adjustment time will increase, and even the system will not work steadily. Contrary to the Ke effect, the smaller the Kec, the larger the overshoot and adjustment time. Ku corresponds to the total magnification of the system. Generally, Ku increases and the rise speed is faster, but when Ku is too large, it will produce a larger overshoot which affects the steady-state operation.
Therefore, using MPGA to optimize three parameters at the same time can give full play to the advantages of global search and find the optimal solution of the three parameters.

Determine the Fitness Function
The GA does not use external information in the optimization search. It only uses the fitness function as the basis and uses the fitness value of each individual in the population to search. In the control system, our commonly used objective function is the integral error criterion (ITAE): where t is the running time of the system, and |e(t)| is the absolute value of the system error. The smaller the ITAE value, the better the performance. It can comprehensively evaluate the control system's response time, overshoot, and other dynamic and static performance. Since the genetic operation is carried out according to the size of the fitness value, and the optimization direction of the objective function should correspond to the direction in which the fitness value increases, the performance index function is appropri- Figure 9. The MPGA structure.

MPGA Design Process Determine the Optimized Parameters
There are many ways to optimize FLC with GA, which can optimize parameters, such as rule base, membership function, scale factor, and so on. Through experimental comparison, it is found that the optimization rule base or membership function has a large amount of calculation, complex parameter settings, and low efficiency; the optimization scale factor is simple, the calculation speed is fast, and it is convenient to modify. Therefore, this project keeps the rule base and membership function unchanged and uses the three scaling factors Ke, Kec, and Ku ( Figure 5) as the parameters to optimize the chromosome coding.
The larger the Ke value, the faster the step response, but at the same time, the overshoot and adjustment time will increase, and even the system will not work steadily. Contrary to the Ke effect, the smaller the Kec, the larger the overshoot and adjustment time. Ku corresponds to the total magnification of the system. Generally, Ku increases and the rise speed is faster, but when Ku is too large, it will produce a larger overshoot which affects the steady-state operation.
Therefore, using MPGA to optimize three parameters at the same time can give full play to the advantages of global search and find the optimal solution of the three parameters.

Determine the Fitness Function
The GA does not use external information in the optimization search. It only uses the fitness function as the basis and uses the fitness value of each individual in the population to search. In the control system, our commonly used objective function is the integral error criterion (ITAE): where t is the running time of the system, and |e(t)| is the absolute value of the system error.
The smaller the ITAE value, the better the performance. It can comprehensively evaluate the control system's response time, overshoot, and other dynamic and static performance.
Since the genetic operation is carried out according to the size of the fitness value, and the optimization direction of the objective function should correspond to the direction in which the fitness value increases, the performance index function is appropriately transformed to obtain the fitness function: Determine the Operating Parameters of the Genetic Algorithm Operating parameters include population size M, population number N, optimal individual minimum retention algebra MAXGEN, elite population ratio ELITE, mutation probability Pm, crossover probability Pc, etc. Among them, the population size, number, and the minimum number of optimal individuals have the greatest impact on the results and calculation time. If it is too small, the optimization result is not ideal, and if it is too large, the iteration will take more time. The mutation probability and the crossover probability determine the system's global and local search capabilities. When they are larger, the recombined individuals will have a high probability of appearing and converge quickly; but at the same time, the replacement of the old and the new is too fast, making some superior individuals may be eliminated prematurely [38][39][40].
In this project, the MPGA parameters are set as M = 20, N = 10, MAXGEN = 10, ELITE = 20%, Pm is a random number between 0.7-0.9, Pc is a random number between 0.001-0.05, and different populations take different probabilities to prevent precocity.

Experimental Verification
On the flexible fast steering mirror (FFSM) experimental platform, the control effects of Cp, T1FLC, and IT2FLC were verified, respectively, and the advantages and disadvantages of the control effects were compared from the four aspects of step response adjustment time, steady-state error, mean square error, and ITAE index. To verify the ability of IT2FLC to deal with uncertainties, three repeated experiments were performed on three control controllers for comparison. In addition, NT, EKM, and EIASC are selected to compare the effects of IT2FLC, and the best TR algorithm in the real-time control system is obtained.

Experimental Platform Characteristics and Cp Design
As shown in Figure 10, the FFSM experiment platform consists of an adjustable base, a flexible mirror, a voice coil motor, a position sensor, a control cabinet, and a communication bus. By adjusting the adjustable base, different characteristics of the controlled object can be obtained. The flexible mirror is driven by a voice coil motor to track the reference input signal. The mirror is made of metal with softer material properties, and it will be slightly deformed during the movement, which will bring uncertainties to the system. The position sensor outputs the mirror position in real time. The communication bus is responsible for the signal transmission, and the control cabinet implements the control operation.
Through frequency sweep processing, the model of the controlled object is obtained. The transfer function G 0 (s) (Equation (26)) consists of two second-order oscillation elements (SOOE), one second-order differentiation element (SODE), one first-order inertial element (FOIE), and a delay element (DE). The data in Table 2 are obtained after identification. The identification effect comparison is shown in Figure 11a.  Through frequency sweep processing, the model of the controlled object is obtained. The transfer function G0(s) (Equation (26)) consists of two second-order oscillation elements (SOOE), one second-order differentiation element (SODE), one first-order inertial element (FOIE), and a delay element (DE). The data in Table 2 are obtained after identification. The identification effect comparison is shown in Figure 11a.
The position controller Cp is designed with a simple and practical zero-pole cancellation method. The open-loop transfer function of the forward path Gcp (s) is obtained as Equation (27).
The measured closed-loop bandwidth under Cp control is 37.83 Hz (Figure 11b), which meets the bandwidth requirements under actual application conditions. ( ) = * + 2 + 1

Cp and T1 and T2 Step Response Comparison
In the experimental platform, Cp, T1FLC, and T2FLC are used for step response experiments. Figure 12 shows the comparison of the response curves of the three control modes. Table 3 shows the settling time (TS), steady-state error (ESS), variance (VAR), and ITAE value of each mode. It can be seen that in the FLC control mode, the four aspects of TS, ESS, VAR, and ITAE index are better than the system with only Cp. T2FLC is better than T1FLC, which reaches a stable state faster and has the least jitter. The measured closed-loop bandwidth under Cp control is 37.83 Hz (Figure 11b), which meets the bandwidth requirements under actual application conditions.

Cp and T1 and T2 Step Response Comparison
In the experimental platform, Cp, T1FLC, and T2FLC are used for step response experiments. Figure 12 shows the comparison of the response curves of the three control modes. Table 3 shows the settling time (TS), steady-state error (ESS), variance (VAR), and ITAE value of each mode. It can be seen that in the FLC control mode, the four aspects of TS, ESS, VAR, and ITAE index are better than the system with only Cp. T2FLC is better than T1FLC, which reaches a stable state faster and has the least jitter.

Cp and T1 and T2 Step Response Comparison
In the experimental platform, Cp, T1FLC, and T2FLC are used for step response experiments. Figure 12 shows the comparison of the response curves of the three control modes. Table 3 shows the settling time (TS), steady-state error (ESS), variance (VAR), and ITAE value of each mode. It can be seen that in the FLC control mode, the four aspects of TS, ESS, VAR, and ITAE index are better than the system with only Cp. T2FLC is better than T1FLC, which reaches a stable state faster and has the least jitter.  It can be seen from Figure 13 that under the guidance of the fuzzy rules when the system starts a step response, the output of the FLC changes accordingly, so the output of  It can be seen from Figure 13 that under the guidance of the fuzzy rules when the system starts a step response, the output of the FLC changes accordingly, so the output of the integrator changes accordingly, and the system type changes dynamically. The step output curve quickly reached a stable state, indicating that the FLC access is effective and the goal of dynamically changing the type has been achieved.

Comparison of the Ability to Handle Uncertainties
Next, compare the uncertainties performance of each control mode. The uncertainties of the experimental platform mainly come from the following aspects:   Table 4 show the step response state in the mode of Cp only. In the case of the same control parameters, the shape of the step response curve is quite different, and the indicators are also very different, which shows that the Cp has a poor ability to handle uncertainties. the integrator changes accordingly, and the system type changes dynamically. The step output curve quickly reached a stable state, indicating that the FLC access is effective and the goal of dynamically changing the type has been achieved. Figure 13. FLC dynamic response.

Comparison of the Ability to Handle Uncertainties
Next, compare the uncertainties performance of each control mode. The uncertainties of the experimental platform mainly come from the following aspects:   Table 4 show the step response state in the mode of Cp only. In the case of the same control parameters, the shape of the step response curve is quite different, and the indicators are also very different, which shows that the Cp has a poor ability to handle uncertainties.

Comparison of the Ability to Handle Uncertainties
Next, compare the uncertainties performance of each control mode. The uncertainties of the experimental platform mainly come from the following aspects:   Table 4 show the step response state in the mode of Cp only. In the case of the same control parameters, the shape of the step response curve is quite different, and the indicators are also very different, which shows that the Cp has a poor ability to handle uncertainties.   Figure 14b and Table 5 show the step response under the T1FLC mode. Its various indicators are more stable and better than the Cp mode, but the dynamic curve shape is   Figure 14b and Table 5 show the step response under the T1FLC mode. Its various indicators are more stable and better than the Cp mode, but the dynamic curve shape is still quite different.  Figure 14c and Table 6 show the step response in T2FLC control mode. The dynamic curve shape difference is small, the steady-state error is reduced to 4 decimal places, and the fluctuation range is not large. In addition, the steady-state accuracy is improved, the adjustment time can be maintained at 2 decimal places, and VAR and ITAE values keep within a small fluctuation range, which is more excellent than the other two control modes. In summary, IT2FLC can better handle system uncertainties than Cp and T1FLC. Figure 15 and Table 7 are comparisons of the three TR algorithms of NT, WTNT, EKM, and EIASC. EKM and EIASC apply the principle of iterative optimization to find the optimal solution, which is more complicated than the WTNT algorithm. From the response results, under the same parameters, both NT, WTNT, and EKM can ensure the stability of the system, but EKM will bring greater jitter, and the values of VAR and ITAE are large and unstable, while EIASC diverges the step response and does not apply to this experimental platform. Comparing NT and WTNT, we can see that the WTNT algorithm has a faster response speed and higher steady-state accuracy. Meanwhile, the VAR and ITAE values of WTNT are smaller than NT, indicating that the WTNT algorithm has stronger robustness. stability of the system, but EKM will bring greater jitter, and the values of VAR and ITAE are large and unstable, while EIASC diverges the step response and does not apply to this experimental platform. Comparing NT and WTNT, we can see that the WTNT algorithm has a faster response speed and higher steady-state accuracy. Meanwhile, the VAR and ITAE values of WTNT are smaller than NT, indicating that the WTNT algorithm has stronger robustness.  On the whole, for this experimental platform, the WTNT algorithm results are the best, and the algorithm itself does not require iteration, the calculation speed is the fastest, and it is more suitable for real-time control systems.

Conclusions
A IT2FLC-DHTC is constructed and verified in an FFSM platform. On the premise of ensuring the stability of the system, the dynamic and steady-state performance is improved, and the goal of dynamic high-type control is achieved.
Adding an integrator to a traditional closed-loop control system can effectively improve the steady-state accuracy of the system. However, the increase of integrators will aggravate the oscillation of the system. In order to avoid system divergence, a dynamic adjustment mechanism is introduced to change the number of integrators with the state of the system, that is, the type of the system. From the experimental results, the dynamic high-type controller can effectively improve the steady-state accuracy of the system while avoiding the system oscillation.
Compared with traditional position controllers, fuzzy logic controller shows great advantages in handling uncertainties, which benefits from the inherent characteristics of FLC. In this paper, a FLC is introduced as the on-off switch of the integrator, and also control the gain of the integrator. By reasonably designing the fuzzy rule base and membership function, the FLC can adjust the on-off and gain of the integrator in real time according to the system state, which effectively improves the system's response speed, steady-state accuracy, and ability to handle uncertainties. On the whole, for this experimental platform, the WTNT algorithm results are the best, and the algorithm itself does not require iteration, the calculation speed is the fastest, and it is more suitable for real-time control systems.

Conclusions
A IT2FLC-DHTC is constructed and verified in an FFSM platform. On the premise of ensuring the stability of the system, the dynamic and steady-state performance is improved, and the goal of dynamic high-type control is achieved.
Adding an integrator to a traditional closed-loop control system can effectively improve the steady-state accuracy of the system. However, the increase of integrators will aggravate the oscillation of the system. In order to avoid system divergence, a dynamic adjustment mechanism is introduced to change the number of integrators with the state of the system, that is, the type of the system. From the experimental results, the dynamic high-type controller can effectively improve the steady-state accuracy of the system while avoiding the system oscillation.
Compared with traditional position controllers, fuzzy logic controller shows great advantages in handling uncertainties, which benefits from the inherent characteristics of FLC. In this paper, a FLC is introduced as the on-off switch of the integrator, and also control the gain of the integrator. By reasonably designing the fuzzy rule base and membership function, the FLC can adjust the on-off and gain of the integrator in real time according to the system state, which effectively improves the system's response speed, steady-state accuracy, and ability to handle uncertainties.
IT2FLC introduces a three-dimensional membership function, which further improves the fuzzy controller's ability to handle uncertainties. From the experimental results, compared with T1FLC, IT2FLC's ability to handle uncertainties is significantly improved. In addition, in order to speed up the calculation speed of IT2FLC, the WTNT type-reduction algorithm constructed in this paper has faster calculation speed and better steady-state accuracy than the traditional TR algorithm, and has been successfully applied to real-time control systems, with good engineering application value.
Future work will continue on the combination of FLC and DHTC. DHTC is a relatively new field. The current mainstream control structure type does not exceed type-3, and the combination of FLC and DHTC does not exceed type-2. If the type can be broken to a higher level, the steady-state accuracy of the system can undoubtedly be further improved, which has important research significance in the engineering field that requires high precision. In addition, the research of the type-2 FLC is also in the preliminary stage, the application field is almost blank, and there are many worthy research areas waiting for us to explore. I believe that in the near future, fuzzy control and high-type control can shine in the field of control, let us wait and see.