1. Introduction
Conventional proportional–integral–derivative (PID) controllers have a simple structure, convenient adjustment, good stability, and reliable operation. Hence, they have been widely used in industrial control applications such as synchronous generator control, aircraft automatic navigation system control, computer numerical control, machine tool control, robot control, temperature control, and liquid level control [
1,
2,
3,
4,
5]. Many intelligent controllers are based on conventional PID controllers. However, when the controlled object is a high-order complex system with a time delay, the performance of PID controllers is greatly reduced, and the control ability may even be lost [
6,
7].
To overcome this problem, many researchers have combined fuzzy logic systems (FLS) with conventional PID controllers [
8,
9]. Type-1 fuzzy PID (T1F-PID) and interval type-2 fuzzy PID (IT2F-PID) controllers have recently been developed [
10,
11]. The structure of the T1-FLS is similar to the IT2-FLS, with the IT2-FLS formed by adding uncertainties to the membership function of the T1-FLS [
12,
13]. Therefore, the IT2-FLS has a stronger control ability than that of the T1-FLS and can effectively control uncertain systems with inherent modeling nonlinearities and time delays [
14,
15]. The combination of the IT2-FLS and PID controllers is a current research hotspot [
16]. For example, IT2F-PID controllers have been successfully applied to high-precision electro-optical tracking [
17], a low-cost microcontroller [
18], a highly complex nonlinear double-link robot operating system [
19], and a mobile robot [
20].
At present, the common structure of IT2F-PID controllers includes cascade and parallel types. The cascade-type IT2F-PID controllers connect the IT2-FLS with a proportional–integral controller in series and perform the proportional–integral operation on the output of the IT2-FLS [
21,
22,
23]. This method does not control the three parameters of the PID controller separately when establishing the fuzzy control rules. So it does not take advantage of the independent control of the three parameters of the PID controller. The parallel IT2F-PID controller structure combines the IT2-FLS and the PID controller in parallel [
6]. Compared with cascade-type controllers, the parallel-type IT2F-PID controller not only retains the benefits of PID controllers but also utilizes the powerful control ability of the IT2-FLS to adjust the control parameters of the PID controller in real time, thus enabling the PID controller to control high-order complex systems. Therefore, the parallel-type structure is more conducive to engineering applications. However, at present, the design of parallel IT2F-PID controllers employs an empirical method, and there is no other optimization approach [
24].
By analyzing the optimal design of cascade-type IT2F-PID controllers, it appears that the construction of IT2F-PID controllers requires appropriate input and output factors and membership function parameters [
23]. Therefore, the appropriate parameters must be selected for parallel-type IT2F-PID controllers. Several researchers have found that the input and output factors have significant effects on the performance of IT2F-PID controllers [
25]. Therefore, on the basis of a fixed membership function, only the input and output factors need to be optimized, which is an efficient method of constructing the controller [
26]. Parallel-type IT2F-PID controllers retain the structure of PID controllers, so there are eight parameters to be optimized, namely, the three coefficients of the PID controller and five scaling factors of the IT2-FLS. Cascade-type IT2F-PID controllers usually have only four parameters to be optimized; even so, it is very difficult to optimize four parameters at the same time using regular intelligent optimization algorithms [
27]. Optimizing the eight parameters, therefore, presents a real difficulty in constructing a parallel-type IT2F-PID controller.
In this paper, a practical disassembled gradational optimization (D-GO) method is proposed for the design of a parallel-type IT2F-PID controller. The main feature of this method is the optimization of eight parameters in two steps. The first step is to optimize the proportional, integral, and differential coefficients of the PID controller. The second step then optimizes the scaling factor of error input, error change rate, proportional coefficient output, integral coefficient output, and differential coefficient output.
The merit of this method is that the optimal PID controller is determined in the first step, so the IT2-FLS only controls the increment of the proportional, integral, and differential coefficients. Compared with the general optimization method, which uses an intelligent optimization algorithm to optimize the parameters at the same time, the optimization range of the parameters is greatly reduced, which significantly reduces the optimization time. The D-GO method provides support for the IT2F-PID controller to realize online real-time parameter adjustment and control time-varying nonlinear complex systems. This paper compares the D-GO method with the concurrent and general optimization methods. Under the same conditions, the controller built using the D-GO method achieves better control performance and has a shorter optimization time. Through simulations, the validity of the parallel-type IT2F-PID controller designed using the D-GO method is verified using a forced closed-loop control system.
The balance of this paper is organized as follows. In
Section 2, the IT2F-PID controller is introduced and the process of constructing a parallel IT2F-PID controller with the D-GO method is described in detail.
Section 3 presents the forced closed-loop control system and the corresponding IT2F-PID controller constructed using the D-GO method. The simulation results are analyzed in
Section 4, and the conclusions of this study are summarized in
Section 5.
2. Interval Type-2 Fuzzy PID Controller
The parallel-type IT2F-PID controller is composed of an upper-level IT2-FLS controller and a lower-level classical PID controller. The upper intelligent controller provides a parameter change selection mechanism for the lower classical PID controller. The parameters of the lower classical PID controller are adjusted using the upper controller. The structure of the parallel IT2F-PID controller is shown in
Figure 1.
The structure of the traditional PID controller is as follows [
6]:
where
kp,
ki, and
kd are the proportional, integral, and derivative gains, respectively, and
u(k),
Fref(k),
Fout(k), and
e(k) are the controller output, expected target value, system output value, and the difference between the target value and the system output value, respectively. In the classic PID controller, the values of
kp,
ki, and
kd are adjusted by technicians according to changes in the system. The control ability of the controller depends heavily on the experience of these technicians, which is not conducive to stable and accurate control [
7]. Therefore, an intelligent IT2-FLS based on PID control rules is designed to adjust the parameters of the classical PID controller online, enabling better adaptation to changes in the system without excessive operator intervention and further enhancing the application range of the traditional controller. The composition and working process of the IT2-FLS is introduced below, and the proposed D-GO method is described in detail.
2.1. Interval Type-2 Fuzzy Logic System
The structure of the upper-level IT2-FLS consists of four parts: a fuzzifier, an inference engine, a rule base, and an output processor. The specific process is shown in
Figure 2 [
28].
The working process of the IT2-FLS is as follows [
14]: the crisp input
e(k) and
de(k)/
dt from the input sensor are first multiplied by the input factors (
Ke and
Kde) and fuzzified into the input interval type-2 fuzzy sets (IT2-FSs). The input IT2-FSs then activate the inference engine and rule base to generate the output IT2-FSs. The IT2-FSs rules are the same as in the T1-FLS, but the antecedents and/or consequents are represented by the IT2-FSs. The inference engine combines the fired rules and gives the mapping from the input IT2-FSs to the output IT2-FSs. The IT2 fuzzy outputs of the inference engine are processed using a type reducer, which combines the output sets and performs a centroid calculation to obtain the T1-FSs, called the type-reduction sets. After the type-reduction process, the type-reduction sets are defuzzified to obtain a clear output. The product of the crisp output and the scale factors (
Kup,
Kui, and
Kud) gives the adjustment values (Δ
kp, Δ
ki, and Δ
kd) in the classic PID controller.
2.1.1. Interval Type-2 Fuzzy Set
An IT2-FS
can be expressed as [
29]:
where
x and
u are the primary and secondary variables, respectively. Jx is the primary membership of
x,
Jx ⊆ [0, 1]. ∬ means taking all possible values of
x and
u. The secondary grades of
are 1, fx(
u) = 1(∀
u ∈
Jx ⊆ [0, 1]).
Figure 3 describes the footprint of uncertainty (FOU) of the fuzzy set
:
where
and
are the lower membership function (LMF) and the upper membership function (UMF), respectively, and
. The LMF and UMF are type-1 membership functions. Like type-1 fuzzy sets, the interval type-2 membership function (IT2MF) can take different forms, such as triangular, trapezoidal, Gaussian, or bell. The membership function used in this paper is triangular, as shown in
Figure 3.
2.1.2. Rule Base
The rules of the IT2-FLS are generally IF-THEN rules that formulate the knowledge of technicians. The rules determine how the IT2-FLS adjusts the parameters of the PID controller. The parameters
kp,
ki, and
kd mainly influence the system through the rise time, overshoot, and settling time. The influence of each parameter is described in
Table 1.
The function of the kp is to speed up the response of the system. As kp increases, the response speed of the system becomes faster and the adjustment accuracy is enhanced. However, it is easy to overshoot the target value, potentially leading to system instability. If the value of kp is too small, the adjustment accuracy of the system will be reduced, the response speed will be slow, the adjustment time will be prolonged, and the static and dynamic characteristics of the system will deteriorate.
The function of ki is to eliminate the steady-state error of the system. With increasing ki, the static error of the system can be eliminated faster, but if the value is too large, the integral saturation phenomenon will occur at the initial stage of the response process, resulting in a large overshoot. If ki is too small, it will be difficult to eliminate the static error, thus affecting the adjustment accuracy of the system.
The function of kd is to improve the dynamic characteristics of the system. Its main function is to suppress changes in the deviation in any direction during the response process and predict the changes in deviation in advance. However, if kd is too large, the adjustment process will be overdamped, which will prolong the adjustment time and reduce the anti-interference performance of the system.
Based on the above PID parameter characteristics, the fuzzy control rules are formulated. The rule base takes the error and the change in error as the premise and the proportional, integral, and differential gains as the result. The structure of the fuzzy rules is as follows:
If E(k) is and (k) is
then Δkp is and Δki is and Δkd is
where
,
,
,
, and
are IT2-FSs, as described in
Figure 3.
2.1.3. Fuzzification and Inference
When , the vertical line at crosses the region to create an interval [, ]. When , the vertical line at crosses the region, creating an interval [, ]. After obtaining two sets of interval values, two firing levels are computed. The lower and upper firing levels are calculated as min[, ] and min[, ], respectively. After obtaining the upper and lower firing levels, the firing interval is obtained as .
2.1.4. Type Reduction and Defuzzification
The type-reduction process is part of the IT2-FLS that differs from the T1-FLS. The process of type reduction involves transforming IT2-FSs into T1-FSs [
30]. The type reduction of IT2-FSs outputs T1-FSs, which are collections of all embedded T1-FSs. In the proposed controller, the most common center-of-set method based on the Karnik–Mendel (KM) algorithm is selected to map IT2-FSs to T1-FSs [
31]. The calculation process is as follows:
where
Ycos(x) is the T1-FS,
yl and
yr are the left and right endpoints of the IT2-FS, respectively.
yl and
yr can be calculated using the KM algorithm [
32]. After the type-reduction process is applied, the clear output value of the IT2-FLS is calculated as:
2.2. Proposed D-GO Method
Based on the structure of the IT2F-PID controller (
Figure 1), at least eight parameters (
kp,
ki,
kd,
Ke,
Kde,
Kup,
Kui, and
Kud) should be optimized to establish a suitable controller. To reduce the parameter optimization time and the hardware requirements, this paper disassembles the parameter adjustment into two categories and performs the optimization gradationally as shown in
Figure 4.
The PID controller of the IT2F-PID controller plays the control role. The IT2-FLS adjusts the PID control parameters. Therefore, we optimize the PID controller parameters (kp, ki, and kd) first, and then optimize the IT2-FLS parameters. This effectively reduces the optimization time. The specific optimization process, optimization algorithm, and objective function are introduced below.
2.2.1. Specific Process of the D-GO Method
The D-GO method proceeds as follows. First, the IT2F-PID controller is simplified to a PID controller, and the PID controller is connected to the controlled system as shown in
Figure 4 (step 1).
At this time, only the PID controller parameters need to be optimized. After the PID parameters (kp, ki, and kd) have been optimized, the IT2-FLS can be established and optimized. Once the optimal PID controller parameters have been determined, the IT2-FLS can adjust Δkp, Δki, and Δkd. Compared with the cascade-type IT2F-PID controller parameters, the optimization range of these parameters is greatly reduced, which significantly reduces the optimization time.
The membership function of the IT2-FLS is uncertain. Hence, optimizing the input and output factors directly will take a lot of time, which affects work efficiency. The structure and mathematical foundations of the IT2-FLS are similar to those of the T1-FLS. Typically, the IT2-FLS is constructed by adding uncertainty to the membership function based on the T1-FLS [
33,
34,
35]. The method of expanding a T1-FLS to an IT2-FLS is theoretically feasible and takes less time.
Therefore, we first establish a T1-FLS, select an appropriate membership function for the fuzzy sets E,
, Up, Ui, and Ud, and then optimize the scale factors (
Ke,
Kde,
Kup,
Kui, and
Kud) of the T1-FLS as shown in
Figure 4 (step 2). After that, the uncertainty of the membership functions in the T1-FLS is added, and then the IT2-FLS is established (step 3). To verify the effectiveness of the proposed D-GO method,
kp,
ki,
kd,
Ke,
Kde,
Kup,
Kui, and
Kud are optimized under the condition that the T1-FLS and IT2-FLS membership functions have been determined. Various intelligent algorithms can be used to optimize the above parameters, including genetic algorithms [
36], particle swarm optimization (PSO) algorithm [
37], grey wolf optimization algorithm [
38], and imperialist competition algorithm [
39]. This paper uses the PSO algorithm as an example to demonstrate the implementation of the D-GO method and optimize the PID controller and IT2-FLS. The optimization objective function is the integral-of-time absolute error (ITAE).
2.2.2. PSO and Objective Function
PSO is a global optimization algorithm inspired by the study of birds’ predation behavior [
37]. The algorithm starts with randomly generated particles in the search space. In each iteration, particles move around the search space at the specified speed, hoping to find the optimal solution. The position of each particle is updated based on the particle’s current position, the particle’s individual optimal solution position, and the global optimal solution position. The particle’s position is then iteratively updated based on the new speed until the stop criterion is reached [
37]. The mathematical expression of the above process is:
where
and
are the current and future positions of the ith particle, respectively. The parameter
is the current speed of the particle, and
w is a weight function that controls the impact of the particle’s current speed on its future speed. The variable
pi,d represents the best solution of particle
i at iteration
t, and
pg,d is the best global solution attained so far. The weight factors c
1 and c
2 determine the importance of the current local optimal solution and the global optimal solution. Finally,
r1 and
r2 are random numbers in the range [0, 1], which further increase the randomness of the particle search in the whole search space.
The selection of a fitness function in the optimization algorithm directly affects the optimization effect; therefore, the commonly used ITAE [
17] is selected as the fitness function. The ITAE is calculated as follows:
4. Simulation Results
This section describes the application of the IT2F-PID controller established using the D-GO method to simulate and control a forced closed-loop system. At present, IT2F-PID controller construction mostly uses the general optimization method, which optimizes the parameters at the same time [
27]. Beyond that, the concurrent method is added for comparison. Comparative analysis was conducted on the optimization time and ITAE of the D-GO and the concurrent and general optimization methods. Four different tasks were performed and the performance was compared with that of a PID controller and IT2F-PID controllers established using the general optimization method and the concurrent optimization method, respectively. To ensure the effectiveness of the comparison, the optimization algorithm and the objective function are the same for each controller. The parameters of the forced closed-loop system are listed in
Table 5.
4.1. Task 1: Normal Case
First, we compare the performance of the PID controller and the IT2F-PID controllers established using the general optimization method, the D-GO method, and the concurrent optimization method, respectively, for a normal case with a 3000 N expected output force. The results are presented in
Figure 9. The parameters (
kp,
ki, and
kd) were optimized in the range [0–1, 0–500, 0–1] and the optimal PID parameters [0.04, 222, 0] were obtained after rounding. The scale factors (
Ke,
Kde,
Kup,
Kui, and
Kud) of the T1-FLS were then optimized. The initial optimization range of these five parameters was defined as [0–0.001, 0–0.000001, 0–0.1, 0–200, 0–0.01], and the optimal parameters were rounded. The IT2F-PID controller was then established as described in
Section 3.2.
The general method is to set up the IT2F-PID controller directly and then optimize its eight parameters (kp, ki, kd, Ke, Kde, Kup, Kui, and Kud) together. The general optimization method is different from the D-GO method. It does not optimize the PID controller first, so it has no preliminary understanding of the system error, error rate of change, or changes in the proportional, integral, and derivative gains. Therefore, the initial parameters (kp, ki, kd, Ke, Kde, Kup, Kui, and Kud) are estimated values. The initial optimization range is [0–1, 0–500, 0–1, 0–1, 0–1, 0–1, 0–500, 0–1] and the optimal parameters are rounded.
The concurrent optimization is to connect the PID controller and IT2-FLS into the controlled system, respectively, for parameter optimization, and then combine them. The PID parameter optimization process is the same as that of the D-GO method, so it is not described here. The IT2-FLS parameter optimization process is similar to that of the general method, which is also not described here. The parameters to be optimized are
Ke,
Kde,
Kup,
Kui, and
Kud and the initial optimization range is [0–1, 0–1, 0–1, 0–500, 0–1]. The IT2F-PID controller built using the D-GO method has obvious advantages over the controllers built using the concurrent and general optimization methods, as shown in
Figure 9. Compared with the other controllers, the controller reduces the overshoot by 80% and the fluctuation by 67%. It has only a small overshoot, no oscillation, and no steady-state error after stabilization. Compared with the IT2F-PID controller, the PID controller has a longer rise time, larger overshoot, and longer settling time. The parameters of the PID controller used are consistent with those of the IT2F-PID controller built using the D-GO method; therefore, the reason for the poor performance is that the PID controller is not connected with the IT2-FLS. Similarly, the IT2F-PID controller established using the general method also has problems like the PID controller and it exhibits sudden fluctuations even after the system has stabilized. Although the fluctuations quickly disappear, they have a certain impact on the stability of the system. The reason for the above problems is that compared with the D-GO method, the initial parameter selection is not accurate enough, resulting in the optimized parameters being not optimal and the control performance being poor. For the D-GO method, a better parameter range was preliminarily determined in the PID parameter optimization step. The controller constructed using the concurrent method does not converge. The control result of the IT2F-PID controller is poor because the PID controller and IT2-FLS are connected to the controlled system, respectively, for parameter optimization during the construction process. The combination of the two causes the IT2F-PID controller to lose control of the controlled system. This illustrates that the proposed controller achieves superior performance.
To illustrate the advantages of the D-GO method, it is compared with the concurrent and general optimization methods in
Table 6. Taking the normal case as an example, the D-GO method has a significantly shorter optimization time than the concurrent and general optimization methods in constructing the IT2F-PID controller and outperforms the concurrent and general optimization methods in terms of the ITAE. The D-GO method reduces the optimization time by over 90% and the ITAE value by about 30% compared with the general method. Compared with the concurrent optimization method, the D-GO method reduces time by over 25% and the ITAE value by about 95%.
By comparing the optimization time of the three methods, it can be found that when the IT2F-PID controller is disassembled into the PID controller and IT2-FLS, fewer parameters are optimized at the same time, and a lot of optimization time can be saved. By comparing ITAE values, it can be found that if the relationship between the PID controller and IT2-FLS is not considered and the controller is only optimized and recombined, the controller will lose control ability. The optimization method that saves time while ensuring control ability is the D-GO method, which first disassembles the controller and analyzes the relationship between the PID controller and the IT2-FLS, and then carries out gradational optimization. It ensures that fewer parameters are optimized at the same time, and the relationship between the two controllers is considered to ensure the control performance.
4.2. Task 2: Model Uncertainty
In practical control applications, the existence of uncertainties such as temperature and friction [
44] makes it difficult to obtain accurate values from the physical model of the forced closed-loop system. Therefore, there is a mismatch between the simulation results and the actual model output. To verify the effectiveness of the D-GO method, the motor inductance parameter and the winding resistance were set to 1.1 times their nominal values. All control parameters remained unchanged.
The results in
Figure 10 show that, in the case of model mismatch, the force tracking performance of the controller of the D-GO method is consistent with that of model matching. The force response curve has only a small overshoot, and there is no steady-state error after a slight oscillation, so the required value is quickly attained and accurately tracked. The PID controller and the IT2F-PID controller established using the general optimization method exhibit worse force tracking ability than the IT2F-PID controller established using the D-GO method. In addition, the IT2F-PID controller established using the general method still exhibits sudden fluctuations and it is more obvious after the system has stabilized. Although the fluctuations quickly disappear, they have a certain impact on the stability of the system. The reasons for the above phenomena are the same as that in normal cases. The PID controller is not connected with the IT2-FLS, so it does not have real-time tuning ability, and its control performance is inevitably poor. The parameters of the IT2F-PID controller established using the general method are not optimal, which leads to poor control performance when the controlled system parameters are uncertain, and the fluctuation after stability is more obvious. In normal cases, the concurrent method has not converged, and it still does not converge in the model uncertainty; therefore, it is not discussed here. This again illustrates that the proposed controller achieves superior performance.
4.3. Task 3: Target Uncertainty
In the actual force control process, there is great uncertainty about the expected output force. Thus, we performed experiments to verify the adaptability of the D-GO method to different expected forces. The expected input was a pulse force with an amplitude of 3000 N, a period of 0.1 s, and a pulse width of 50%. This input is significantly different from that used to establish the controller and is intended to reflect the effectiveness of the proposed method.
The results in
Figure 11 show that the IT2F-PID controller designed using the D-GO method achieves excellent force tracking performance under the peak of pulse force. The force tracking fluctuates at the bottom, so it can be seen that the optimization parameters obtained in the normal case cannot be applied to any situation. Compared with the other controllers, the proposed IT2F-PID controller reduces the overshoot by 80% and the fluctuation by 67%. The force response curve has a small overshoot, and quickly tracks the required value thereafter. The performance is obviously better than that of the PID controller and the IT2F-PID controller designed using the general optimization method. Because the PID controller does not have the ability to adjust its own parameters, the parameters of the IT2F-PID controller designed using the general method are far from the optimal value, resulting in an obvious overshoot and a significant increase in the stability time of both.
4.4. Task 4: External Disturbance
The forced closed-loop control system is inevitably subjected to external force interference during operation. Under external disturbance, the validity of the D-GO method is tested. The increase in disturbance force was Fd = 300 N in the normal case. The parameters of each controller remain unchanged.
The force tracking performance is shown in
Figure 12. Under external disturbance, the controller constructed using the D-GO method has a short rise time, small overshoot, and short settling time in the process of force tracking, and the overall tracking performance is obviously better. Compared with the PID controller and the controller constructed using the general method, the overshoot was reduced by 80%, the fluctuation was reduced by 67%, and especially the recovery speed after interference was obviously better. The reason for the poor performance of the controller established using the general method is similar to the normal case. The PID has no ability to adjust its own parameters, and the controller parameters established using the general method have a wide range of optimization, which leads to a long distance from the optimal value and poor control performance. The parameters of the controller constructed using the D-GO method are closer to the optimal value and the control performance is better.
5. Conclusions
This paper has described the D-GO method for designing a parallel-type IT2F-PID controller. The proposed method disassembles the parallel IT2F-PID controller into a PID controller and a T1-FLS for parameter optimization and then blurs the membership function of the T1-FLS to construct a parallel IT2F-PID controller. Compared with the general optimization method, the D-GO method reduces the optimization time by over 90% and the ITAE value by about 30%. Compared with the concurrent optimization method, the D-GO method reduces time by over 25% and the ITAE value by about 95%. The D-GO method significantly reduces the optimization time and supports the IT2F-PID controller to realize online real-time parameter adjustment and control a time-varying nonlinear complex system.
Simulations were conducted to verify the control performance of the D-GO method. Through comparisons with a PID controller and an IT2F-PID controller established using the concurrent and general optimization methods, we found that in the normal case, model uncertainty, target uncertainty, and external disturbance, the D-GO method’s controller achieves superior control performance for the forced closed-loop system. The overshoot is reduced by 80% and the fluctuation is reduced by 67% compared with the PID controller and the parallel IT2F-PID controller established using the concurrent and general optimization methods. In addition, compared with the other controllers, the D-GO method’s controller does not produce sudden fluctuations after the system has been stable and the recovery speed after interference is obviously better. This illustrates that the D-GO method’s controller achieves superior performance. The results presented in this paper are sufficient to demonstrate that the D-GO method is suitable for the establishment of IT2F-PID controllers. By extension, it also provides inspiration for the establishment of other intelligent controllers.