Maximum Power Tracking Control of Wind Turbines Based on a New Prescribed Performance Function

: The primary control goals of a wind turbine (WT) are structural load shedding, maximum wind energy capture in the underpowered situation, and consistent power production in the full power condition. A crucial component of the control problem for wind turbines with varying speeds is maximum power tracking control. Conventional maximum power tracking control tracks the ideal blade tip speed ratio to provide the most wind power at the speciﬁed wind speeds. However, because of the wind turbine’s great nonlinearity and the signiﬁcant external disturbances it encounters, it is difﬁcult to react quickly to variations in wind speed, and the tracking speed is sluggish, which lowers the amount of electricity produced annually. In light of this, this work develops a novel preset performance controller for a wind power system maximum power tracking control. With this technique, the convergence rate and tracking precision may be set. In particular, based on the concept of time-varying feedback, a time-varying function, known as the preset performance function, is ﬁrst created to allow the convergence speed and accuracy to be predetermined; then this time-varying function is used to transform the actual speciﬁed time problem of the original system into a bounded time problem of the new system; ﬁnally, a direct robust controller design strategy with pre-deﬁned performance is suggested based on the design concept of the backstepping technique. The plan may maximize the rotor power coefﬁcient by altering the wind turbine speed, track the ideal blade tip speed ratio for a given tracking accuracy and speed, and get the most wind power to produce the most power with the strongest robustness. The simulation results show that the recommended control technique works.


Introduction
Wind energy is a green and renewable energy source, and the process of generating electricity from wind does not produce any carbon emissions.It is one of the power generation methods with the greatest potential [1].Wind power generation will unavoidably play a significant role in the worldwide development of clean energy in order to fulfill the objectives of peak carbon and carbon neutrality, accelerate the proportion of renewable energy consumption, decrease the use of fossil fuels, and build a green and hospitable environment.The control of wind turbines has also received a great deal of attention in Europe and in some other countries around the world.
Generally speaking, there are two control loops to control a wind turbine, namely torque control loop and pitch control loop [2]. Figure 1 shows that, depending on the wind speed, the control of a wind turbine may be roughly split into three phases, i.e., the standby segment, the optimal Cp segment, and the rated power segment.The green line depicts the turbine's actual operating curve, while the blue line shows that the turbine is running at its ideal blade tip speed ratio.Region 1 denotes that the turbine is in standby mode, since the wind speed is too low to generate electricity, with the rotational speed Energies 2023, 16, 4022 2 of 21 rising and the torque at 0. Region 2 is the optimal Cp section, where the blade tip speed ratio is ideal and the turbine has not yet reached its rated output.The turbine is controlled to maintain constant power by altering the pitch angle and generator torque in Region 3, which shows that the turbine has reached rated power and that the rotational speed is inversely proportional to the torque.Area 1 1/2 indicates the start-up state and area 2 1/2 the transition state where the rated power is about to be reached, these two areas account for very little.The region 2 is under power operation.To utilise wind energy as much as possible in this operating region is the control aim.At this stage, the blade pitch angle is locked at the optimal value in order to optimize power production, and the torque control circuit is used to adjust the rotor speed to run relatively near to the ideal speed.Numerous control techniques, including tip speed ratio (TSR) control [3], optimal torque control [4], power signal feedback control [5], and extremum seeking control [6], etc., have been researched and used in this sector.However, their characteristics include a slow tracking speed and that tracking accuracy cannot be set in advance.
Energies 2023, 16, x FOR PEER REVIEW 2 of 22 running at its ideal blade tip speed ratio.Region 1 denotes that the turbine is in standby mode, since the wind speed is too low to generate electricity, with the rotational speed rising and the torque at 0. Region 2 is the optimal Cp section, where the blade tip speed ratio is ideal and the turbine has not yet reached its rated output.The turbine is controlled to maintain constant power by altering the pitch angle and generator torque in Region 3, which shows that the turbine has reached rated power and that the rotational speed is inversely proportional to the torque.Area 11/2 indicates the start-up state and area 21/2 the transition state where the rated power is about to be reached, these two areas account for very little.The region 2 is under power operation.To utilise wind energy as much as possible in this operating region is the control aim.At this stage, the blade pitch angle is locked at the optimal value in order to optimize power production, and the torque control circuit is used to adjust the rotor speed to run relatively near to the ideal speed.Numerous control techniques, including tip speed ratio (TSR) control [3], optimal torque control [4], power signal feedback control [5], and extremum seeking control [6], etc., have been researched and used in this sector.However, their characteristics include a slow tracking speed and that tracking accuracy cannot be set in advance.In control engineering, a common control goal is trajectory tracking [7].The convergence rate and the final tracking error range, two crucial metrics for assessing the controller, must frequently be taken into account when developing a controller.The tracking error needs to converge to a specified precision within a given period for several application scenarios that have higher criteria for safety and reliability, such as target interception, spacecraft docking, and auto parts assembly [1,8].The latest results of prescribed time control [9][10][11][12][13] solve the problem of system stability in a limited time.The designer can predefine the convergence time, which is a parameter that is independent of the initial situation.Using this method, the control gain is engineered to grow infinitesimally as time moves closer to the terminal time T, resulting in the desired temporal stability [14].It is challenging to build controllers for systems with the majority of uncertainties and external disturbances in order to achieve zero error tracking.In fact, in practical engineering applications, the tracking error is allowed to have an appropriate accuracy range rather than zero error.Due to this, a few scholars have proposed a useful theory of finite-time control in which the tracking error converges in limited time to a local neighborhood (rather than zero) of the origin [15][16][17][18][19][20][21][22][23].The tracking error after the stabilization time is uncertain since the size of this neighborhood is unknown, and its value actually depends on some unknowable factors.In the field of engineering, this is quite undesirable.
In order to obtain the maximum power output of a wind turbine at low wind speeds, we attempt to create a new prescribed performance controller (PPC) based on the discussion above.No matter what the initial conditions are, users can fully specify the convergence speed and tracking accuracy in advance.In order to allow the designer to set the convergence rate and accuracy in advance, we first create a time-varying function, known as a prescribed performance function (PPF), with embedded stability time and tracking In control engineering, a common control goal is trajectory tracking [7].The convergence rate and the final tracking error range, two crucial metrics for assessing the controller, must frequently be taken into account when developing a controller.The tracking error needs to converge to a specified precision within a given period for several application scenarios that have higher criteria for safety and reliability, such as target interception, spacecraft docking, and auto parts assembly [1,8].The latest results of prescribed time control [9][10][11][12][13] solve the problem of system stability in a limited time.The designer can predefine the convergence time, which is a parameter that is independent of the initial situation.Using this method, the control gain is engineered to grow infinitesimally as time moves closer to the terminal time T, resulting in the desired temporal stability [14].It is challenging to build controllers for systems with the majority of uncertainties and external disturbances in order to achieve zero error tracking.In fact, in practical engineering applications, the tracking error is allowed to have an appropriate accuracy range rather than zero error.Due to this, a few scholars have proposed a useful theory of finite-time control in which the tracking error converges in limited time to a local neighborhood (rather than zero) of the origin [15][16][17][18][19][20][21][22][23].The tracking error after the stabilization time is uncertain since the size of this neighborhood is unknown, and its value actually depends on some unknowable factors.In the field of engineering, this is quite undesirable.
In order to obtain the maximum power output of a wind turbine at low wind speeds, we attempt to create a new prescribed performance controller (PPC) based on the discussion above.No matter what the initial conditions are, users can fully specify the convergence speed and tracking accuracy in advance.In order to allow the designer to set the convergence rate and accuracy in advance, we first create a time-varying function, known as a prescribed performance function (PPF), with embedded stability time and tracking accuracy.Then, we use this time function to transform the actual specified time problem of the original system into the bounded time problem of the new system.Finally, derived from the design idea of the backstepping method [24], we propose a direct robust controller (DRC) design scheme.
The contributions of this paper include: (1) The convergence rate and tracking precision can both be predetermined by the designer utilizing the new error transformation based on PPF; (2) Not only the tracking error variable, but also all other intermediate error variables will reach the predetermined stable region within a given time; (3) In spite of external disturbance, the WT can respond to the change of wind speed quickly and achieve maximum power control.

Problem Formulation
Take into account the following uncertain nonlinear strict-feedback systems: where R i is the system state; u ∈ R and y are the system input and the system output, respectively; f i (x i ) and g i (x i ) represent the uncertain nonlinear continuous function and unknown gain function, respectively; d i (.) denotes the unknown time-varying disturbances.
Designing a correct u to achieve the following objectives is the main task of this paper: Certain assumptions are made in order to accomplish the aforementioned control aims.
Assumption 1.The gain functions g i (x i ) are positive and there are unidentified constants 0 < g i < g i that case 0 < g i ≤ g i (x i ) ≤ g i .
Assumption 2. For lumped uncertainties there exist unknown constants a i and known functions Assumption 3. The unknown time-varying disturbances d i (.) are bounded, and there exist unknown positive constant D i satisfy |d i (.)| ≤ D i .
Assumption 4. The desired trajectory y d and its i(i = 1, . . ., n)-order derivatives are known, piecewise continuous and bounded.

Some Lemmas and Definitions
Lemma 1.For any |r| < 1, the following inequality holds:

Lemma 2 ([11]
).For constants l > 0 and time-varying functions v(t) > 0, there are: x = f (x) with x being system state, is said to be reached if for any initial condition x(0) = x 0 , there exist constant ε > 0 and finite time T < ∞ (in which T is independent of systems initial conditions and other design parameter vector), such that |x| < ε for all t > T.
It can be seen from Definition 1 that the accuracy ε in the semi-global practical prescribed time (PPT) stability is not a specific value.In other words, there is no explicit need for the size of this domain; just that the system's state converge to a finite error domain within a certain amount of time.The following practical predefined time and precision stability is defined as the tracking error being restricted to a predetermined compact set around zero within a stipulated period T.

Definition 2 ([26]
).If a time-varying function ρ(t) holds the following generalized properties, then this function is called a new prescribed performance function (NPPF): ρ(t) is a continuous and non-increasing C n function from an initial ρ(0) = ρ 0 to a terminal value ρ(T) = ε, where ρ 0 , ε, T are given as constants greater than 0. .
Immediate examples of such ρ(t) include the subsequent formulae with t ≥ 0: where the constant m meets the m ≥ 2 requirement.It should be noted that there are many functions that satisfy Definition 2, which are not limited to the above examples.
We build the following lemma to continue the control design because it will be vital to our later technological development.

Definition of a New Variable
Lemma 3. We firstly define a new variable h = r ln 1 − r 2 , where r = e/ρ, e = y − y d .If we can design an appropriate controller u(.) so that h bounded with the condition of the initial state satisfy |e(0)| < ρ 0 , then we can get −ρ < e < ρ.

Modelling of Wind Turbines
A wind turbine is made up of four basic parts, including aerodynamics, mechanics, a generator, and a pitch actuator, as depicted in Figure 2. The two discs on the left and right in Figure 3 are the determinant of inertias, J r and J g .In other words, J r and J g represent the rotor's and the generator's respective rotational inertia.The rotation of the rotor shaft at high torques is described as a mass-spring system with damping D s and a spring constant K s .In addition, T r and T g stand for the torque on each side of the rotor and generator, respectively.The transmitted torques T rs and T gs are interconnected via the gear ratio N g .
A wind turbine is made up of four basic parts, including aerodynamics, mechanics, a generator, and a pitch actuator, as depicted in Figure 2. The two discs on the left and right in Figure 3    The aerodynamic model of a wind turbine is where ρ represents air density; R represents rotor radius; and v represents wind speed; the rotor power coefficient, denoted by the symbol , is as follows [27]: where ω d is the optimal generator speed.The aerodynamic model of a wind turbine is where ρ represents air density; R represents rotor radius; and v represents wind speed; the rotor power coefficient, denoted by the symbol C p (λ, β), is as follows [27]: ) where ω d is the optimal generator speed.
Energies 2023, 16, 4022 The single mass block model of the wind turbine drive train is as follows where J is the total moment of inertia of the wind turbine; K a = 1/2ρπR 5 C p /λ 3 is wind turbine's operational state coefficient; T g is generator torque; B is the total damping coefficient of the transmission system; r g is the gear ratio.You may rewrite System (6) as where

Error Transformation
Based on Assumption 4 in Section 2.1, the tracking error is defined as: Obviously, the control objective (b) holds if the controller can be designed so that the tracking error satisfies the following inequality: Define the intermediate error variable as: where ∂ i−1 stand for output signal of virtual controller ∂ i−1 to be designed through the first-order filter with variable time constant as follows: Filtering errors are defined as follows: To guarantee that the error variables w i are constrained to the specified precision ε at the same rate, and that the system's state x i does not go beyond a predetermined bound, we define the transformation error variable: where r i = vw i , v = 1/ρ.Take the derivative of h i , the following transformation is performed on system (1): where and controller u(.) are designed to make the conversion error variable h = [h 1 • • • h n ] bounded, then we will get −ρ < w i < ρ, i = 1, . . ., n.The PPT control issue of ( 1) is simplified to building a controller such that the state of the system shown in ( 15) is limited, hence simplifying the previously complex control problem.

Direct Robust Control (DRC) Scheme
The backstepping method divides a complex nonlinear system into subsystems no larger in size than the system itself, designs partial Lyapunov functions for each subsystem, "backstepping" all the way to the whole system, and then integrates them to finish designing the entire control law.
Then, the time derivative of 1 2 h 1 2 along ( 17) is It is simple to obtain that with the aid of Assumptions 1-3 in Section 2.1 and Young's inequality [28]: Adding the aforementioned disparities' two sides gives us where Hence the virtual control is constructed as where k 1 > 0. Using Lemma 1 in Section 2.2, we can get: Energies 2023, 16, 4022 8 of 21 Combining ( 16) and ( 21), we get A candidate for a Lyapunov function is now defined as Differentiating (24) and combining ( 20)-( 23), we have where φ 1 is the filtering error, and its derivative is: .
Each variables of the function κ 1 are in the compact sets, at the same time, κ 1 is a smooth function, it follows that the largest |κ 1 | is κ 1 , i.e., |κ 1 | < κ 1 with κ 1 being an unknown constant.Thus, we have: Substituting ( 27) into (25), we have Step i (i = 2, . . ., n − 1).Similar to step 1, combined with the system (15), 1  2 h i 2 derivative of time is: The virtual controller is designed as Define a positive definite Lyapunov function Differentiating (31) and combining (29), we have .
With the help of Young's inequality and also considering the approach described in (19), we have where Using a similar approach as in Step 1, we have .
and the largest |κ 1 | is κ 1 , i.e., |κ i | < κ i with κ i being an unknown constant.Therefore, it is possible to obtain that: Substituting ( 33) and ( 35) into (32), we get: . where Step n.Select the following choice for the Lyapunov function Then, we obtain .
V n as .
Design the actual control u as Consider (37) and Assumptions 1-3 in Section 2.1, with the help of Young's inequality, there is Substituting (39) and (40) into (38), we have . where .
Figure 4 shows the structure of the DRC controller designed in this paper.Where, ω is the actual generator speed, ω d is the desired generator speed, T g stand for the electromagnetic torque, ε and T are the pre-settable tracking accuracy and convergence speed respectively.

System Feasibility Analysis
Theorem 1. Think about the closed-loop system (1) that complies with Assumptions 1-3 in Section 21), (30) with the first-order variable time constant filter (12) and the actual controller u in (39) are constructed, then the subsequent goals are accomplished: (a) The closed-loop system is stable.
(b) The tracking error converges to the designated area By combining ( 28), (36) and (41), it follows that where To solve the above differential equation, we have: Using Lemma 2 in Section 2.2, we can get:

System Feasibility Analysis
Theorem 1. Think about the closed-loop system (1) that complies with Assumptions 1-3 in Section 2.1.If the virtual controllers ∂ i , i = 1, . . ., n − 1 in ( 21), ( 30) with the first-order variable time constant filter ( 12) and the actual controller u in (39) are constructed, then the subsequent goals are accomplished: (a) The closed-loop system is stable.
(b) The tracking error converges to the designated area Ω = {w i ∈ R, i = 1, . . ., n : |w i (t)| < ε} within the designated periodT, whereεandTare predetermined values provided by the user.(c) Each and every intermediate signal has a limit, i.e., bounded.
Proof of Theorem 1.Consider the following Lyapunov function: By combining ( 28), (36) and (41), it follows that . where To solve the above differential equation, we have: Using Lemma 2 in Section 2.2, we can get: We are now prepared to establish the following theorem's conclusions.Firstly, we show that objective (a) is achieved.Since we may deduce from (45) that V ∈ L ∞ holds for every constrained beginning condition, it follows that h i ∈ L ∞ , φ i ∈ L ∞ .Thus, the closed-loop system is stable, the signals h i and φ i are ultimately uniformly bounded.
Next, we show that objective (b) is achieved.Since h i is bounded, we can derive −ρ < w i < ρ according to Lemma 3 in Section 2.3.The tracking error obeys the performance constraint |w i (t)| < ρ and converges to the specified region Ω = {w i ∈ R, i = 1, . . ., n : |w i (t)| < ε} within the predetermined time T, as shown by the fact that ρ(t) as defined in Definition 2 of Section 2.2 will fall to ε within finite time T.

Simulation Study
In this paper, Matlab (Version 9.12, R2022a) is used for simulation studies and Simulink (Version 10.5, R2022a) is used for controller design; the simulation model built is shown in Figure 5. Table 1 displays the precise wind turbine specifications [29].The simulation model includes wind speed disturbance uncertainty, and the control effect of the designed direct robust controller (DRC) is compared with that of the sliding mode controller (SMC) commonly used in maximum power tracking control.The turbulent wind model is produced using the Turbsim simulator, the use of which can be found in the OpenFAST|TurbSim User Manual.The resulting turbulent wind has an effective wind speed of 10 m/s and is used to simulate natural winds., it follows that i w are constrained, which also suggests that λ i and π i are constrained.Note that i r is bounded, then it fellows from ( 9), ( 11), ( 13), ( 21), ( 30), ( 34) and (39) that the virtual con- and actual controller u are bounded.The proof is com- pleted.□

Simulation Study
In this paper, Matlab (Version 9.12, R2022a) is used for simulation studies and Simulink (Version 10.5, R2022a) is used for controller design; the simulation model built is shown in Figure 5. Table 1 displays the precise wind turbine specifications [29].The simulation model includes wind speed disturbance uncertainty, and the control effect of the designed direct robust controller (DRC) is compared with that of the sliding mode controller (SMC) commonly used in maximum power tracking control.The turbulent wind model is produced using the Turbsim simulator, the use of which can be found in the OpenFAST | TurbSim User Manual.The resulting turbulent wind has an effective wind speed of 10 m/s and is used to simulate natural winds.We now confirm that the suggested direct robust controller can provide the needed tracking precision |e(t)| = |x(t) − y d (t)| < 0.1 radians within the required time T = 2 s.The WT is first set up with x(0) = 0.5 radians, since x(0) − y d (0) = 0.5, the prescribed performance functions is chosen as: where ρ 0 = 2, ε = 0.1, T = 2.In simulation, the direct robust control algorithm (32) is utilized to deduce the control input u.
Figure 6 shows the turbulent winds with effective speed being 10 m/s. Figure 7 displays the controller's tracking performance for tracking accuracy ε = 0.1 and convergence time T = 2.In Figure 7, (a) shows that the control method achieves sensitive tracking of the generator speed ω against the optimal generator speed ω d under random fluctuations in wind speed; (b) demonstrates that the tracking error achieves the control goal of error tracking by converging to the required accuracy range in the provided convergence period t = 2 s; (c) demonstrates how the tip speed ratio is kept to vary about the ideal value even after the tracking error has converged to the defined range; (d) demonstrates how the rotor power coefficient is maximized and the greatest amount of wind energy is captured by altering the tracking wind turbine's ideal speed; (e) and (f) show that the control method enables fast tracking of its optimal trajectory by aerodynamic torque and wind power; (g) shows that the speed is varied by adjusting the generator torque to cope with wind speed variations.Fast tracking is accomplished via the control approach.
gence time = 2 T .In Figure 7, (a) shows that the control me ing of the generator speed ω against the optimal genera fluctuations in wind speed; (b) demonstrates that the track goal of error tracking by converging to the required accura vergence period t = 2 s; (c) demonstrates how the tip speed ideal value even after the tracking error has converged to t strates how the rotor power coefficient is maximized and energy is captured by altering the tracking wind turbine's that the control method enables fast tracking of its optim torque and wind power; (g) shows that the speed is varie torque to cope with wind speed variations.Fast tracking is approach.The impact of various controller settings on the system is then confirmed.The effect of convergence time T on the system is first verified.Figure 8 shows the control performance of the controller for a fixed tracking accuracy ε = 0.1 and varying the controller convergence time T to 1, 2, and 3 respectively.(a) shows that the magnitude of the convergence time T affects how fast the wind turbine speed tracks to the desired speed trajectory; the smaller T, the faster the tracking speed, where ω 1 , ω 2 , and ω 3 are the wind turbine speeds at convergence times T = 1, 2, and 3, respectively, and ω d stands for desired speed; (b) demonstrates how the size of the convergence time T affects how quickly the tracking error converges to within a certain accuracy; the smaller the T, the quicker the convergence, where ε 1 , ε 2 , and ε 3 correspond to convergence times T of 1, 2, and 3 respectively; (c) and (d) show that the smaller the convergence time T, the shorter the time taken to reach the ideal blade tip speed ratio and the highest rotor power coefficient.Panels (e), (f) and (g) show that the smaller the convergence time T, the shorter the time to track the optimal aerodynamic torque and wind power trajectory, and the faster the generator torque is involved in the regulation process.The impact of various controller settings on the system is then confirmed.The effect of convergence time T on the system is first verified.Figure 8 shows the control perfor- mance of the controller for a fixed tracking accuracy ε = 0.1 and varying the controller convergence time T to 1, 2, and 3 respectively.(a) shows that the magnitude of the con- vergence time T affects how fast the wind turbine speed tracks to the desired speed tra- jectory; the smaller T , the faster the tracking speed, where ω 1 , ω 2 , and ω 3 are the wind turbine speeds at convergence times T = 1, 2, and 3, respectively, and ω d for de- sired speed; (b) demonstrates how the size of the convergence time T affects how quickly the tracking error converges to within a certain accuracy; the smaller the T , the quicker the convergence, where ε , ε , and ε correspond to convergence times T magnitude of the tracking accuracy ε affects the final convergence range of the tracking error, tracking accuracy is increased by decreasing the ε, where the values of ε 1 , ε 2 , and ε 3 are 0.05, 0.1, and 0.15, respectively, and ρ 1 , ρ 2 , and ρ 3 are the corresponding performance bounds at different tracking accuracies; (c) and (d) show the different fluctuations of the blade tip speed ratio and rotor power coefficient around the optimal value at different tracking accuracies, the larger the ε, the greater the fluctuation amplitude, and vice versa, the flatter the fluctuation amplitude; (e) and (f) show that the smaller the tracking accuracy ε, the higher the coincidence of optimal aerodynamic torque and wind power trajectory tracking, i.e., the higher the tracking accuracy; (g) shows how the smaller the tracking accuracy ε, the more sensitive the regulation of the generator torque and the greater the regulation when the wind speed changes rapidly.
Energies 2023, 16, x FOR PEER REVIEW 16 of 22 ε = 0.05, 0.1, and 0.15 respectively, and ω d is the desired turbine speed; from (b) it is clear that the magnitude of the tracking accuracy ε affects the final convergence range of the tracking error, tracking accuracy is increased by decreasing the ε , where the values of ε 1 , ε 2 , and ε 3 are 0.05, 0.1, and 0.15, respectively, and ρ 1 , ρ 2 , and ρ 3 are the corresponding performance bounds at different tracking accuracies; (c) and (d) show the different fluctuations of the blade tip speed ratio and rotor power coefficient around the optimal value at different tracking accuracies, the larger the ε , the greater the fluctuation amplitude, and vice versa, the flatter the fluctuation amplitude; (e) and (f) show that the smaller the tracking accuracy ε , the higher the coincidence of optimal aerodynamic torque and wind power trajectory tracking, i.e., the higher the tracking accuracy; (g) shows how the smaller the tracking accuracy ε , the more sensitive the regulation of the generator torque and the greater the regulation when the wind speed changes rapidly.Next, we verify that preset performance tracking (PPT) can be achieved.We demonstrate that, given various beginning circumstances, the suggested controller can successfully perform PPT tracking.Figure 10 shows the tracking performance when ε = 0.01, T = 1.In Figure 10, compared to Figure 7a, (a) shows that the generator speed ω tracks the ideal speed ω d trajectory more sensitively after changing the initial conditions; the tracking error converges to a smaller error range more quickly, as seen in (b), improving the tracking accuracy; (c) and (d) show that after changing the initial conditions, the rotor power coefficient and tip speed ratio fluctuate less and stay close to their ideal levels; (e) and (f) show that the tracking sensitivity of the aerodynamic torque and wind power to their desired values becomes higher after changing the combination of controller parameters; (g) shows how the improved control accuracy results in a more rapid and sensitive adjustment of the generator torque to cope with rapid changes in wind speed.The control objective is achieved.Next, we verify that preset performance tracking (PPT) can be achieved.We demonstrate that, given various beginning circumstances, the suggested controller can successfully perform PPT tracking.Figure 10 shows the tracking performance when ε = 0.01 , = 1 T .In Figure 10, compared to Figure 7a, (a) shows that the generator speed ω tracks the ideal speed ω d trajectory more sensitively after changing the initial conditions; the tracking error converges to a smaller error range more quickly, as seen in (b), improving the tracking accuracy; (c) and (d) show that after changing the initial conditions, the rotor power coefficient and tip speed ratio fluctuate less and stay close to their ideal levels; (e) and (f) show that the tracking sensitivity of the aerodynamic torque and wind power to their desired values becomes higher after changing the combination of controller parameters; (g) shows how the improved control accuracy results in a more rapid and sensitive adjustment of the generator torque to cope with rapid changes in wind speed.The control objective is achieved.Subsequently, we varied the parameters of the controller so that ε = 0.15, T = 3 to verify the effect of different parameter combinations on the performance of the controller.Figure 11 shows the tracking performance when ε = 0.01, T = 1. Figure 11 compares with Figure 7 as well as Figure 10.The corresponding panel (a) show a decrease in the sensitivity of the turbine speed ω to the desired trajectory ω d tracking and a decrease in trajectory coincidence after increasing the tracking accuracy and convergence time T. Panel (b) shows a slowdown in the convergence of the tracking error to within a given accuracy range, with a larger accuracy range leading to a decrease in control accuracy; panels (c) and (d) show controller parameter combinations show a significant increase in the fluctuation of the leaf tip speed ratio and rotor power coefficient around the optimum value.It can be seen that the selection of suitable controller parameters is crucial to improve the control effect; (e) and (f) show a reduction in sensitivity to optimal trajectory tracking of aerodynamic torque and wind power after increasing the convergence time T and tracking accuracy ε; (g) shows the slower and less regulated involvement of the generator torque in the regulation.Choosing the right controller parameters is crucial.
In this paper the controllers are selected as direct robust controllers and sliding mode controllers, and the simulation curves are compared.The SMC controllers used are as follows [30]: First define a sliding surface: s = .
where y represents the system status, the control laws of the system are as follows produces .s = 0.For a zero starting condition, the system converges exponentially s(t 0 ).By introducing a discontinuous term as KSgn(s) in the control u 0 for K > 0, the system's closed-loop dynamics are brought closer to a sliding surface.The system will arrive to the sliding surface in a limited amount of time at s = 0.If a parametric estimate of f , f is provided in the presence of uncertainties, the control rule changes to û = KSgn(s) + u.Thus there are: This results in .s = f − f − KSgn(s).It is clear from Figure 12 that the DRC controller is more effective.(a) and (b) show that under the strong disturbance of external perturbation, the DRC control makes the tracking error converge quickly to the specified accuracy range, achieving sensitive tracking of the required speed trajectory and ensuring the best speed of the generator, while the tracking error under the SMC control gradually maintains stability after 35 s and the fluctuation amplitude exceeds the given accuracy range; (c) and (d) show that the DRC controller is able to track the TSR more accurately and quickly, and achieves the tracking of the optimal tip speed ratio within the preset time of 2 s, so that the rotor power coefficient reaches the maximum value to obtain the best wind power and achieve the control targets, while the SMC controller has a slower tracking speed and the control effect is unstable, and the TSR and C p fluctuate more around the optimal value; (e) and (f) show that the DRC controller is more responsive to rapid changes in wind speed than the SMC controller, enabling maximum capture of wind energy by maximising aerodynamic torque and wind power within a pre-given time of T = 2 s.
the fluctuation of the leaf tip speed ratio and rotor power coefficient around the optimum value.It can be seen that the selection of suitable controller parameters is crucial to improve the control effect; (e) and (f) show a reduction in sensitivity to optimal trajectory tracking of aerodynamic torque and wind power after increasing the convergence time In this paper the controllers are selected as direct robust controllers and sliding mode controllers, and the simulation curves are compared.The SMC controllers used are as follows [30]: First define a sliding surface: where y represents the system status, the control laws of the system are as follows of the optimal tip speed ratio within the preset time of 2 s, so that the rotor power coeffi-cient reaches the maximum value to obtain the best wind power and achieve the control targets, while the SMC controller has a slower tracking speed and the control effect is unstable, and the TSR and Cp fluctuate more around the optimal value; (e) and (f) show that the DRC controller is more responsive to rapid changes in wind speed than the SMC controller, enabling maximum capture of wind energy by maximising aerodynamic torque and wind power within a pre-given time of T = 2 s.

Conclusions
In this study, a DRC control strategy for WT systems' maximum power tracking is suggested.The control scheme can be pre-given arbitrary parameters and can make all tracking errors converge to within a specified accuracy range ε within a time T given by the designer.Experiments prove that the size of the convergence time T determines how fast the tracking error converges to within the given accuracy range; the smaller the T , the faster the speed.While the tracking accuracy ε determines the final convergence range of the tracking error, i.e., it affects the accuracy of the control, the smaller the ε , the higher the control accuracy.Choosing the right tracking accuracy ε and conver- gence time T can achieve sensitive tracking of the error.For wind turbine systems, the control scheme can track the optimal blade tip speed ratio for a given tracking accuracy

Conclusions
In this study, a DRC control strategy for WT systems' maximum power tracking is suggested.The control scheme can be pre-given arbitrary parameters and can make all tracking errors converge to within a specified accuracy range ε within a time T given by the designer.Experiments prove that the size of the convergence time T determines how fast the tracking error converges to within the given accuracy range; the smaller the T, the faster the speed.While the tracking accuracy ε determines the final convergence range of the tracking error, i.e., it affects the accuracy of the control, the smaller the ε, the higher the control accuracy.Choosing the right tracking accuracy ε and convergence time T can achieve sensitive tracking of the error.For wind turbine systems, the control scheme can track the optimal blade tip speed ratio for a given tracking accuracy and speed, and achieve fast and accurate tracking of the optimal turbine speed trajectory, e.g., within a convergence time T = 1 s, so that the tracked speed error converges to an accuracy range of 0.01 for the upper and lower limits, and the rotor power coefficient reaches its maximum.Further, the tracking accuracy can be changed to reduce the fluctuation of wind energy utilization coefficient and TSR near the ideal value, with the goal to achieve the greatest power tracking control and the most wind energy collection, with strong robustness in the face of complex wind speed changes, and the control input is bounded and smooth.However, the wind turbine model is simplified for research, and a complicated non-linear wind turbine model will be incorporated in further work to confirm the controller's effectiveness in terms of control.

( a )
In closed-loop systems, the intermediate signals are constrained; (b) The preset performance tracking (PPT) is achieved, i.e., there are matching parameter values for every ε > 0 and T > 0 ensuring |w 1 (t)| = |y(t) − y d (t)| < ε when t ≥ T, where y d (t) is the target reference signal.
are the determinant of inertias, r J and g J .In other words, r J and g J represent the rotor's and the generator's respective rotational inertia.The rotation of the rotor shaft at high torques is described as a mass-spring system with damping s D and a spring constant s K .In addition, r T and g T stand for the torque on each side of the rotor and generator, respectively.The transmitted torques rs T and gs T are interconnected via the gear ratio g N .

Figure 2 .
Figure 2. Schematic diagram of the wind turbine structure.Figure 2. Schematic diagram of the wind turbine structure.

Figure 2 . 22 Figure 3 .
Figure 2. Schematic diagram of the wind turbine structure.Figure 2. Schematic diagram of the wind turbine structure.Energies 2023, 16, x FOR PEER REVIEW 6 of 22

Figure 4
Figure4shows the structure of the DRC controller designed in this paper.Where, ω is the actual generator speed, ω d is the desired generator speed, g T stand for the electromagnetic torque, ε and T are the pre-settable tracking accuracy and conver- gence speed respectively.

Figure 4 .
Figure 4. Diagram of the controller.
period T , where ε and T are predetermined values provided by the user.(c) Each and every intermediate signal has a limit, i.e., bounded.Proof of Theorem 1.Consider the following Lyapunov function:

Figure 4 .
Figure 4. Diagram of the controller.

Figure 10 .
Figure 10.Tracking performance under ε = 0.01, T = 1.(a) Generator speed; (b) Tracking error; (c) Tip speed ratio; (d) Rotor power coefficient; (e) Aerodynamic torque; (f) Wind power; (g) Generator torque.Subsequently, we varied the parameters of the controller so that ε = 0.15 , = 3 T to verify the effect of different parameter combinations on the performance of the controller.

Figure 11
Figure11shows the tracking performance when ε = 0.01 , = 1 T. Figure11compares with Figure7as well as Figure10.The corresponding panel (a) show a decrease in the
1, , : t as defined in Definition2of Section 2.2 will fall to ε within finite time ρ i w

Table 1 .
Specific parameters of wind turbines.