A Tuning Method via Borges Derivative of a Neural Network-Based Discrete-Time Fractional-Order PID Controller with Hausdorff Difference and Hausdorff Sum

: In this paper, the fractal derivative is introduced into a neural network-based discrete-time fractional-order PID controller in two areas, namely, in the controller’s structure and in the parameter optimization algorithm. The ﬁrst use of the fractal derivative is to reconstruct the fractional-order PID controller by using the Hausdorff difference and Hausdorff sum derived from the Hausdorff derivative and Hausdorff integral. It can avoid the derivation of the Gamma function for the order updating to realize the parameter and order tuning based on neural networks. The other use is the optimization of order and parameters by using Borges derivative. Borges derivative is a kind of fractal derivative as a local fractional-order derivative. The chain rule of composite function is consistent with the integral-order derivative. It is suitable for updating the parameters and the order of the fractional-order PID controller based on neural networks. This paper improves the neural network-based PID controller in two aspects, which accelerates the response speed and improves the control accuracy. Two illustrative examples are given to verify the effectiveness of the proposed neural network-based discrete-time fractional-order PID control scheme with fractal derivatives


Introduction
Due to the introductions of integral-order and differential-order, the design of fractionalorder proportion-integral-derivative (PID) controllers is more flexible compared with that of integral-order PID controllers, allowing better control, which has attracted attention from many scholars [1].Due to the simple structure of fractional-order PID controllers, as effective fractional-order controllers, they have been widely used in many control systems and achieved satisfactory control performances [2].Take the examples of power systems [3], TCP networks [4], DC motors [5], supercapacitor energy storage systems [6], and robotic manipulator systems [7].In addition, a variety of mixed fractional-order PID controllers have been derived from other algorithms, such as the fuzzy PID controller [8] and fractional-order PID sliding mode controllers [9].
The tuning methods of parameters in fractional-order PID controllers mainly include the analytic method and intelligent optimization algorithms.The analytic method can give the optimal parameters of fractional-order PID controllers with the given performance indexes.Generally, the method of analytic method is to design the parameters of the fractional-order controller in the frequency domain [10].The dominant pole placement was used to tune the parameters of fractional-order PID controller to obtain a satisfactory load disturbance response in [11].The phase margin and gain crossover frequency were concerned to present a tuning method of fractional-order PID controllers for minimum phase plants [12].The robustness on the loop gain variations was guaranteed to tune the parameters of fractional-order PID controller in [13].Although the analytic method can provide the mathematical expression of the parameter tuning, this method depends on the mathematical model of the controlled plant.It is difficult to use for the parameter tuning of fractional-order PID controllers of complex nonlinear controlled plants.Therefore, many scholars have studied the parameter tuning methods of fractional-order PID controllers based on various intelligent optimization methods, such as the genetic algorithm [14], the chaotic atom search optimization algorithm [15], the quantum bacterial foraging algorithm [16], the extremum seeking algorithm [17], and the particle swarm optimization algorithm [18].Although the above intelligent optimization algorithms can obtain the effective parameter tuning results for the nonlinear controlled plant, the parameter tuning methods are offline calculation methods.
As an intelligent optimization and decision algorithm, the neural network algorithm is also applied in the design of fractional-order PID controllers.In [19], the RBF neural network algorithm was used to tune the parameters of fractional-order PID controllers and applied for a wind turbine.However, it is also an offline parameter tuning method and requires a lot of offline data.The difficulty of online parameter tuning methods of fractionalorder PID controllers is mainly the adaptive tuning of fractional-orders, due to the global operators of fractional-order derivatives in the Caputo sense or Riemann-Liouville sense, and the fractional-order difference in the Gr ünwald-Letnikov (G-L) sense.The fractionalorder operators also play a role in neural networks.Many scholars have proposed various fractional-order neural networks defined by fractional-order derivatives [20][21][22][23][24].There are many research results in the cross research of neural networks and fractional-order calculus.Hence, it is worth studying the application of neural network optimization algorithms to the parameter tuning of fractional-order PID controllers.
In order to emphasize the adaptive adjustments of fractional-orders, the parameters of fractional-order PID controller in this paper include the parameters the proportional coefficient, the integral coefficient, and the differential coefficient; and orders integral order and differential order.The neural network algorithm is also used to adjust the parameters and orders of fractional-order PID controllers.It is necessary to study the discretized method of fractional-order operators to design a discrete-time fractional-order PID controller.There are many discretization methods of fractional-order derivatives [25,26]; the simplest discretization form of a fractional-order derivative is the G-L difference [27].The discretization method via the G-L difference is simple, but the mathematical expression of the discretization result contains the Gamma function related to the order.If the orders in fractional-order PID controllers are adjusted adaptively by a neural network, the derivatives of the Gamma function with respect to the orders need to be obtained.Unfortunately, the derivative of the Gamma function is not a closed-form mathematical expression [28], and it not easy to calculate directly to hinder the adaptive updating of orders.
As a special fractional-order derivative, the fractal derivative deals with the local characteristics.The fractal derivatives are widely used in the modeling of real physical systems [29], such as anomalous diffusion [30] and viscoelastic materials [31].Meanwhile, the validity of the fractal derivative for achieving chaotic attractors was discussed in [32] to show that the nonlinear system described by fractal derivative can also produce complex, dynamic behavior.The Hausdorff derivative as one of the fractal derivative definitions does not contain the Gamma function [33].Therefore, a novel discrete-time fractional-order PID controller can be constructed by using the discretized forms of the Hausdorff derivative and the Hausdorff integral, and the adaptive tuning method of parameters and orders can be realized by using neural networks.Besides, the power-law stability of nonlinear systems via the Hausdorff derivative was investigated in [34].The fractal derivative mainly generalizes the first-order derivative with many definitions.The Borges derivative is another form of fractal derivative which is given in [35], and it can also describes the dynamic behavior of complex systems.
In this paper, fractal calculus theory is applied to the design of a neural network PID controller in two ways.According to the Hausdorff derivative and the Hausdorff integral, the definitions of the Hausdorff difference and Hausdorff sum are derived.Compared with the fractional-order PID controller via the G-L definition, the derivation of the Gamma function is avoided in parameter updating based on neural networks.In addition, the fractional-order PID controller based on fractal theory has no memory of difference, and the updating of parameters and orders is simpler.Another improvement is to introduce the Borges derivative into the optimization of parameters and orders.The derivation of a composite function based on the Caputo derivative is very complex; hence, the chain rule of the derivative of am integer-order composite function cannot be used directly.However, the chain rule of the Borges derivative of a composite function is consistent with the integral-order derivative, and the effect can also be optimized by adjusting the order.Therefore, this paper applies the Borges derivative to the training of controller parameters and orders, and discusses the problem of order selection in the Borges derivative.
The main contributions of this paper are summarized as follows: (1) The definition and calculation method of Hausdorff difference and Hausdorff sum are proposed, and their introduction into PID controllers is described.(2) A novel type of fractional-order PID controller is proposed to realize the online adjustment of parameters and orders, and the online tuning method of parameters and orders based on the first order gradient information is given.(3) The Borges derivative is introduced into the tuning of parameters and orders, and the influence of the order in the Borges derivative on the optimization performance in regard to parameters and orders is offered.(4) The tuning methods of the parameters and orders of the proposed fractional-order PID controller based on the fractal derivative is proposed, and it improves the optimization speed of the parameters and orders.
The rest of this paper is organized as follows.In Section 2, the definitions of Hausdorff difference and sum derived from Hausdorff derivative and integral are given, and the fractional-order PID controller in the sense of a Riemann-Liouville (R-L) derivative and Caputo derivative with zero-initial-condition is also given.Based on the Hausdorff difference and sum, two kinds of neural network-based tuning methods are proposed for the parameters and fractional-orders of a fractional-order PID controller via the first order and the Borges derivative in Section 3. In Section 4, two illustrative examples are offered to verify the effectiveness of this paper.The conclusion is given in Section 5.

Problem Statement
In this paper, a novel type of fractional-order PID controller is constructed based on the difference and sum operators extended by the Hausdorff derivative and Hausdorff integral.The Hausdorff derivative with α ∈ R + for t > 0 is defined in [33] as where R + is the set of positive real numbers.It is not difficult to calculate the Hausdorff derivative if the function is a differentiable function as follows: The Hausdorff difference for k > 0 derived from the Hausdorff derivative is given as follows.
Similarly, the Hausdorff integral is defined in [32] based on (2) as follows.
Consequently, the Hausdorff sum for k > 0 derived from the Hausdorff integral is determined by For α = 1 in (3) and ( 5), the Hausdorff difference and Hausdorff sum become Thus means that the Hausdorff difference and Hausdorff sum include the first order difference and sum operations.Therefore, the main investigated issue of this paper is the tuning design of parameters and orders in the fractional-order PID control scheme with neural networks based on the Hausdorff difference and Hausdorff sum.
As the fractional-order PID controller has two more important parameters that are the integral-order and the differential-order, the controller design is more flexible compared with the integer-order PID controller.Generally, the basic definitions of fractional-order derivatives are the R-L definition and Caputo definition.In a continuous-time, fractionalorder PID controller, the mathematical descriptions under these two definitions are usually used.The definitions are given as follows.
The Caputo derivative [25] is defined by where C 0 D α t is Caputo derivative from 0 to t with the order α ∈ R + , and it satisfies P − 1 < α ≤ P with P ∈ Z + ; Z + is the set of positive integer numbers, and the function Γ(•) is the gamma function.
Accordingly, the R-L derivative [25] is defined by where RL 0 D α t is R-L derivative from 0 to t with the order α ∈ R + and satisfies P − 1 ≤ α < P with P ∈ Z + .
If the initial conditions are zero for both the Caputo derivative and the R-L derivative; the Laplace transformation of fractional-order derivative with respect to the function f (t) with the fractional-order α is s α F(s) with the Laplace operator s, where F(s) is the Laplace transformation of f (t).Then, the transfer function of fractional-order PID controller [36] in the frequency domain is offered as follows. where , and µ ∈ (0, 2) are the proportional coefficient, the integral coefficient, the differential coefficient, the integral order, and the differential order; and R + is the set of positive real numbers.The fractional-order operation in the discrete-time fractional-order PID controller is described by G-L definition.The G-L difference is provided as follows.
where the factor c α j is computed by Generally, the fractional-order operators in the fractional-order controller (8) are discretized by the G-L difference, but the coefficients of the discretized fractional-order operators are all related to the orders, and the relationship between the corresponding coefficient and the order is more complex due to the Gamma function in terms of the order.Therefore, the adaptive adjustments of corresponding orders λ and µ are quite difficult via the G-L difference.As the Hausdorff derivative, and its discretization are both local operators, the orders in the discrete-time fractional-order PID controller based on the Hausdorff difference and Hausdorff sum can be adjusted adaptively.

A Discrete-Time Fractional-Order PID Controller Based on the Hausdorff Difference and Hausdorff Sum
Since the fractional-order operators in fractional-order PID controller defined by Equation ( 11) are global operators, the adaptive adjustments of orders are not easy to calculate.Therefore, a fractional order PID controller based on the Hausdorff derivative and Hausdorff integral is proposed in this paper.The corresponding controller in the time domain is offered as follows: where e(t) is the output error.By using the Hausdorff difference and Hausdorff sum described by Equations ( 3) and ( 5), the discrete-time, fractional-order PID controller is given as follows.
where e(k) is the sampling value of e(t) at t = kT with the sampling period T, and The influences of different orders on the control effect based on Hausdorff difference and the Hausdorff sum in the fractional-order controller (12) are discussed.Then, we define K λ (k) = λk λ−1 and K µ (k) = k 1−µ /µ in u(k) respectively, and the curves of K λ (k) and K µ (k) are drawn in Figures 1 and 2 respectively.It can be seen from the curves described in Figures 1 and 2 that K λ (k) increases and K µ (k) decreases with an increase of fractional-order λ or µ.If the order λ is set as a larger value, the role of the current time is larger, and the role of the historical information of the error is smaller.That is to say, more attention is paid to the error information of the current time in H Σ λ k e(k).However, a larger order λ makes the sum term too large; hence, it is necessary to choose a reasonable order λ.Correspondingly, the order µ is set as a smaller value; then the effect of the difference becomes larger.Therefore, the order µ should be made smaller with the increase of the iteration k.
Next, we give the calculation method of the parameters and orders of the discrete-time fractional-order PID controller (12) trained by the neural network method.

Tuning Method of Parameters and Orders of Discrete-Time Fractional-Order PID Controllers with Neural Networks
The neural network-based discrete-time fractional-order (NNDTFO) PID controller with the neural network identifier (NNI) is the focus, and its structure is shown in Figure 3.

Learning algorithm
Learning algorithm The nonlinear controlled plant is presented by where g(•) is the nonlinear function in terms of the control input u(k) and output y(k), n y ∈ Z + , n u ∈ Z + , and Z + is the set of positive integer numbers.
The weight matrices corresponding the connection from the input layer to the hidden layer and the connection from the hidden layer to the output layer are denoted as W (1) (k) ∈ R S 1 ×N and W (2) (k) ∈ R 1×S 1 , where N = n y + n u + 1 is the number of the neurons in the input layer.
The input of NNI is where is the output of the hidden layer, f (•) is activation function of the hidden layer, and b (1) (k) and b (2) (k) are the biases of the hidden layer and output layer.
Let the cost function ; then the updating formulas of the weights and biases in NNI are determined by where ζ is the learning rate of NNI, and s (1) (k) and s (2) (k) are determined by 1 (k))/∂r (1) 2 (k))/∂r (1) S 1 (k))/∂r (1) , and r (1) i,j (k) is the entry in W (1) (k) corresponding to the ith row and the jth column; and x j (k) is the jth entry in the input x(k).
A BP neural network with one hidden layer is adopted as the NNI for the controlled plant (13).Let the controller for k > 0 be where e(k) = r(k) − y(k), r(k) is the reference signal, and Define the cost function as J 2 (k) = 0.5e 2 2 (k) with e 2 (k) = r(k) − ŷ(k); then the tuning formulas of parameters are determined by where η k , η i , and η d are the learning rates for the parameters k p , k i , and k d .Consider that then the tuning formulas of orders are also determined by where η λ and η µ are the learning rates for the orders λ and µ.
There is the function term ŷ(k)/∂u(k) in the above updating formulas of the parameters and orders; hence, the calculation method of ŷ(k)/∂u(k) is as follows: S 1 ,n y +1 (k)] T .
Therefore, the updating formulas of parameters and orders of NNDTFO PID controller (19) are determined by (2) (k)F (1) (k)W (1) (k)e(k), Considering the effective selection intervals of parameters and orders in the controller (19), we restrict the optimization results of parameters and orders as

Tuning Method of Parameters and Orders Based on the Borges Derivative
In the updating Formulas (45)-( 49), the updating methods of parameters k p (k), k i (k), and k d (k) and orders λ(k) and µ(k) depend on the information of the first order negative gradient.We investigate the updating method based on the improved Borges gradient to improve the optimizations of parameters and fractional-orders in NNTDFO PID controller (19).
The Borges derivative is provided in [35] as follows.
If the function f (x) is differentiable, we have If we set q = 1 in (31), the Borges derivative is equivalent to the first-order derivative.Therefore, we use the Borges gradient consisting of the Borges derivatives to generalize the updating Equations ( 45)-(49).In addition, the Borges derivative is a local operator, which satisfies the chain rule; hence, the training of parameters and orders of the controller (19) based on the neural network algorithm can be applied.Let z = m(v) and v = n(x); then the Borges derivative of z = m(n(x)) is obtained as Based on (32), we have where q p , q i , q d , q λ , and q µ are orders which are positive real numbers.
In order to give the adjustment methods of these orders, we take the order of the updating formula of the proportional coefficient as an example, and the relative speed is defined as The relationship between the index L(k p (k)) and k p (k) is shown in Figure 4.It is obvious from Figure 4 that the Borges derivative of J 2 (k) with respect to k p (k) is lager if the order q k p is set to be smaller.If q k p is larger than a certain value, the sign of derivative changes, and this is an undesirable result.With the increase of parameter k p (k), the index L(k p (k)) corresponding to a smaller order becomes larger, while the index L(k p (k)) corresponding to a larger order becomes smaller, and it may be less than zero.
Let q p (k) = 0; then we have L(k p (k)) = 1 + k p (k), which is the maximum of L(k p (k)).Thereby, the range of L(k p (k)) is (0, 1 + k p (k)].If the indexes L(k p (k)) corresponding to different parameters k p (k) are the same as L m ∈ (0, 1 + k p (k)], we have We hope that the Borges derivative is less than the first-order derivative to accelerate the optimization at the beginning of parameter optimization.Meanwhile, the Borges derivative is greater than the first-order derivative to improve the accuracy of parameter optimization at the end of the optimization.Then, the order q p (k) is determined by where ; rp is a random number uniformly distributed between 0 and 1; and M is the number of iterations.
Similarly, the tuning methods of other fractional-orders are offered as follows: where ri , rd , rλ , or rµ is a random number uniformly distributed between 0 and 1; and the conditions 0 Based on the Equations ( 33)-(37), the updating formulas of parameters and orders in the controller (19) based on the Borges derivative are offered as follows. (2)(k)F (1) (k)W (1) (k)e(k), (45) (49)

Illustrative Examples 4.1. Example 1
Consider the following discrete-time nonlinear controlled plant as The reference signal r(k) is chosen as r(k) = 0.2 + 0.1 sin(0.02k),and the initial value y(0) = 0 is set.The NNI is used to identified the model of the nonlinear plant and the activation function in the hidden layer of NNI is selected as In the NNI, the number of the neurons in the hidden layer is S 1 = 6, and N = 2 is obtained with x(k) = [y(k), u(k)] T and L = 500.To verify the effectiveness of the proposed NNDTFO PID controller, the results were compared with those of a neural network discretetime (NNDT) PID controller.The learning rates in NNI were both set as ζ = 0.5 for the two types of controller.The learning rates of the parameters in two types of controllers were both set as η p = 5000, η i = 10, and η d = 10, and the learning rates of the orders were set as η λ = 10,000 and η µ = 10,000.The initial value for weights and biases in NNI was chosen as 0.01, and the weight in the NNDTFO PID controller and NNDT PID controller was set as 0.1; u(0) = 0 was used.Meanwhile, the initial value of the orders λ and µ was set as 1 in the NNDTFO PID controller.
The responses of the outputs y(k) via two types of controllers, and the responses of two types of controllers u(k), are shown in Figures 5 and 6, respectively.
It can be seen from the output responses via the different controllers described in Figure 5 that the NNDTFO PID controller with Hausdorff difference and Hausdorff sum can obtain a more accurate response compared with the output response via the NNTD PID controller.From the curves of the two types of PID controllers described in Figure 6, it can be seen that the outputs of the two types of PID controller are bounded with no significant control effects.In order to show that Hausdorff sum and difference operations can reduce the response error, the index function is defined as While letting r(k) = 0.2 + A sin(ωk), the parameters were set as ω = 0.01, 0.02, • • • , 0.5 with A = 0.2, and A = 0.1, 0.15, • • • , 0.3 with ω = 0.02 respectively; then the corresponding indexes E are offered in Tables 1 and 2. It can be seen from the data in Tables 1 and 2 that the response error of system output can be reduced by introducing the Hausdorff difference and Hausdorff sum for different reference signals.These results show that the proposed NNDTFO PID controller can further improve the accuracy of the output response.
To illustrate the improvement of the proposed controller in regard to the dynamic response of the controlled plant, the reference input signal was changed to r(k) = 0.3 with L = 100.
From the experimental results in Figure 7, it is obvious that the controller with the Hausdorff sum and difference operations can make the response of the output y(k) faster, which is another advantage of the control strategy proposed in this paper.Besides, the parameters and orders in the NNDTFO PID controller are shown in Figure 8 respectively.
From the curves of control parameters and orders in Figure 8, it can be seen that parameters k p (k), k i (k), and k d (k), and orders λ(k) and µ(k) converge to the fixed values quickly.The orders λ(k) and µ(k) in the proposed controller are greater than 1 or less than 1 for the reference signal r(k) = 0.3, and this phenomenon corresponds to the analysis of orders in Section 3.

Example 2
In this example, the Borges derivative is focused on to improve the control performance based on the NNTD PID controller.The nonlinear controlled plant used was y(k) = −0.3y(k− 1) + 0.1 cos(y(k − 1)) + u(k − 1). (53) The parameters of the NNI and NNDTFO PID controller were consistent with the corresponding ones in example 1.Other parameters were set as , and L µ 2 = 0.9 + µ(k).In consideration of the existence of random functions with the use of the Borges derivative to optimize the parameters and orders of the controller, the average result of 100 experiments was used for the evaluation in this case.The reference signal was set as r(k) = 0.2 + 0.1 sin(0.02k)with L = 600; then the response of the average value of output y(k) via the NNDTFO PID controller with Borges derivative, and the responses of outputs y(k) via the NNDTFO PID controller and NNDT PID controller were plotted in Figure 9.The corresponding errors based on (52) via the NNDTFO PID controller via Borges derivative, the NNDTFO PID controller, and the NNDT PID controller were 0.0729 (the average result for 100 times), 0.0892, 0.1296 respectively.
It can be seen from the experimental results in Figure 9 that the NNDTFO PID controller with Borges derivative achieved the best control effect in three types of PID controller.By using Borges derivative, the output of the system can track the reference signal well; the tracking errors via the three controllers were the smallest.This phenomenon illustrates that the introduction of the Borges derivative can improve the optimization effect of parameters and orders in NNTDFO controller.For this example, we can draw the conclusion that the use of the Hausdorff difference and Hausdorff sum can improve the control accuracy.On this basis, the control accuracy can be further improved by introducing the Borges derivative.It was verified that the control effect of the NNTD controller can be improved by introducing fractal derivatives into two areas, i.e., the controller structure and the optimization strategy.
Figure 10 shows the average values of the relative speeds corresponding to the parameters and orders of the NNTD PID controller with the Borges derivative from 100 experiments.It can be seen from the average values of relative speeds corresponding to different parameters and orders described in Figure 10 that all of the average relative speeds can be gradually reduced from a number greater than 1 to a number less than 1, by adopting the order adjustment strategy proposed in this paper.According to this strategy, the optimization speed can be accelerated at the beginning of optimization stage for parameters and orders in the NNTDFO PID controller, and the optimization precision and effect can be improved at the end of the optimization stage.
To compare the robustness of the three types of controllers, the coefficient 0.1 before the nonlinear function term cos(y(k − 1)) in the model was changed to 0.1 + B sin(0.5k), and the error for different values of parameter B are shown in Table 3. (The index of the NNDTFO PID controller is the average result of 100 rounds).The experimental results in Table 3 show that the three controllers based on neural networks have robustness and can adapt to disturbances in parameters.In addition, the fractal derivative can reduce the control error E to improve the control accuracy of the closed-loop system even, if parameter 0.1 is disturbed.
In order to investigate the response of the Borges derivative to the convergence speed of the control effect, the reference signal was set to r(k) = 0.25, and L = 100 was set for this case; other parameters in three types of controllers were not changed.Then, from the average value of output y(k) via the NNDTFO PID controller with Borges derivative, the responses of outputs y(k) via NNDTFO PID controller and NNDT PID controller are drawn in Figure 11.The control effect of the three controllers described in Figure 11 shows that the NNTDFO PID controller with the Borges derivative achieved the fastest convergence speed, followed by the NNTDFO PID controller without the Borges derivative, and the NNTD PID controller without the fractal derivative had the slowest convergence speed.It is shown that the Borges derivative can further accelerate the response of the controlled plant.
Next, we study the error index for a differing constant reference signal, and the results are given in Table 4. From the experimental results described in Table 4, it can be seen that the rule of the control effect is consistent for constant input signals of different values.In other words, the use of a fractal derivative to optimize controller parameters and the use of fractal difference and fractal sum can improve control effect.
To verify the fact that the designed NNI can identify the system model well based on three kinds of neural network PID controllers, the errors between the output of NNI ŷ(k) and the output y(k) as E x (k) = | ŷ(k) − y(k)| for L = 100 and r(k) = 2.5 are shown in Figure 12.It can be seen from the error curves in Figure 12 that NNIs are not affected by the three types of neural network PID controllers, and NNIs can identify the dynamic characteristics of the controlled plant effectively, and the identification error converges to zero quickly.

Conclusions
This paper mainly studied the applications of two kinds of fractal derivative in the design of a neural network-based PID controller.Based on the definition of the Hausdorff derivative and Hausdorff integral, the calculation formulas of Hausdorff difference and Hausdorff sum were given to construct NNTDFO PID controllers.By using the local operation characteristics of Hausdorff difference and Hausdorff sum, the tuning method of parameters and orders in the proposed NNTDFO PID controllers were given.The orders corresponding to the Hausdorff difference and Hausdorff sum were introduced to make the controller design more flexible, and the adaptive adjustment algorithm of orders is also provided in this paper.
In addition, the learning rules of parameters and orders were further discussed.The first derivative was replaced by the Borges derivative; then the parameters and orders with the adaptive relative velocity were adjusted to improve the optimization performance.The adaptive tuning rules of the orders in the Borges derivatives were given to gain a more satisfactory control effect via NNTDFO PID controllers.The use of a fractal derivative, like the integer derivative, satisfies the chain rule, which provides the basis for the training parameters of neural networks.The applications of fractal difference and fractal sum derived from the fractal derivative and integral can avoid the derivation of the Gamma function, which is convenient for calculations and can also increase the flexibility of controller design.In this paper, two kinds of fractional-order PID controller based on fractal derivative and derived fractal difference and fractal sum were proposed, which can improve the response speed and accuracy of the closed-loop control system.Compared with the integer-order PID algorithm, the calculation is more complex and the calculation time is longer.However, in view of the fact that many controllers are operated by computers, the calculation burden of the proposed algorithm is not large.In future research, the basic neural network learning and construction method based on fractal derivatives will be extended to a more complex neural network model and neural network control algorithm.
Funding: This research was funded by Liaoning Revitalization Talents Program under grant number XLYC1807229, the Natural Science Foundation of Liaoning Province, China under grant number 20180520009, the China Postdoctoral Science Foundation Funded Project under grant number 2019M651206, and the Scientific Research Fund of Liaoning Provincial Education Department, China under grant LJC202010.

Figure 3 .
Figure 3. Structure of the control system with a neural network-based discrete-time fractional-order (NNDTFO) PID controller and a neural network identifier (NNI).

Figure 4 .
Figure 4. Curves of L(k p (k)) with different orders q k p (k).

Figure 5 .
Figure 5. Responses of outputs y(k) for the closed loop system of example 1 via two types of controller.

Figure 6 .
Figure 6.Responses of controllers u(k) for the closed loop system of example 1 with a sinusoidal signal.

Figure 7 .
Figure 7. Responses of outputs y(k) closed loop system for example 1 with a constant signal.

Figure 9 .
Figure9.Responses of outputs y(k) via three types of PID controller of a closed loop system for example 2 with a sinusoidal signal.

Figure 10 .
Figure 10.Relative average speeds of parameters and orders.

Figure 11 .
Figure 11.Responses of outputs y(k) via three types of PID controller using the closed loop system of example 2 with a constant signal.

Figure 12 .
Figure 12.Errors E x (k) for three types of neural network PID controller.

Table 1 .
Index E for different parameters ω with A = 0.1 via two types of controller.

Table 2 .
Index E for different parameters A with ω = 0.02 via two types of controller.

Table 3 .
Index E for different parameters B via three types of controller.

Table 4 .
Index E for constant reference signal via three types of controller.