Suppression of Wind Generator Speed Vibration Based on the Internal Model Control with Three Degrees of Freedom

: When the wind power system is working in the maximum wind energy tracking area, or when its mechanical drivetrain system vibrates and the vibration frequency is relatively high, the system cannot suppress the speed vibration of the generator in the drivetrain system by adjusting the pitch angle. In this paper, a generator speed control system based on the three degrees of freedom internal model control (3-DOF-IMC) is established to control the generator speed. Thus, a scheme of using the feedback ﬁlter in the 3-DOF-IMC to suppress generator speed vibration caused by drivetrain shaft elasticity and gear clearance in the drivetrain system is proposed. Firstly, the vibration problem and waveform of the two-mass wind power drivetrain systems are discussed, and the generator’s vector control and speed control systems are analyzed. Secondly, the principle of the 3-DOF-IMC is described, and the inﬂuence of the three controllers on the speed tracking performance and anti-interference performance of the generator is discussed. The suppression ability of the feedback ﬁlter for different forms of disturbance signals is emphasized. Finally, the feedback ﬁlter controls the generator speed and eliminates the inﬂuence of drivetrain vibration on the generator speed. To verify the superiority of the proposed method, the vibration suppression effect, tracking performance and anti-interference performance of the proposed method are compared with the engineering design method (EDM) and conventional IMC (1-DOF-IMC) method. When the parameter α / β = 0.66, the generator speed amplitude overshoot of the proposed method is the same as the EDM. When α / β = 1, it is only 4% of the amplitude overshoot of the EDM. In addition, the tracking performance and anti-interference performance of the proposed method can be adjusted independently, and it is better than the EDM and the 1-DOF-IMC method. tracking performance, anti-interference performance, and robustness of the speed loop.


Introduction
The wind turbine, gearbox, generator, and drivetrain shaft constitute the wind power drivetrain system. The essence of vibration in the wind drivetrain system is caused by the torque difference between the wind turbine's pneumatic torque and the generator's electromagnetic torque [1]. The main method for suppressing the vibration of the wind drivetrain system is reducing the torque difference by adjusting the pitch angle of the wind turbine or by controlling the electromagnetic torque of the generator.
There are two specific suppression methods. In the first method, a mechanical drivetrain system based on the wind power system adjusts the aerodynamic torque of the wind turbine through variable pitch control. In the second method, the generator is controlled by adjusting its electromagnetic torque. The vibration frequency caused by the elasticity of the drivetrain shaft is relatively low in the wind mechanical drivetrain system. The pitch-speed closed loop of the wind power system can be used to suppress the speed vibration of the generator. However, the vibration frequency caused by gear clearance is still high [2,3], and the variable pitch actuator cannot keep up with the rapid changes imposed by the variable the necessary gearbox equivalence and conversion of the gearbox's rotational inertia and elastic coefficient.
In addition, the wind turbine's PID control cannot suppress the vibration resulting from gear clearance. Therefore, the electromagnetic torque or speed of the wind generator can be controlled by the torque-speed closed loop to eliminate the wind drivetrain vibration with clearance. Figure 1 depicts a schematic diagram of the wind drivetrain system composition for a doubly fed wind-powered system [24]. The doubly fed generator is connected to the drivetrain system via a high-speed shaft. Thus, the drivetrain system's performance affects the generator's output performance. In addition, the generator is connected to the grid through the transformer, while the generator's electrical control system controls the converter that regulates the electromagnetic torque and speed of the generator. A doubly fed wind generator is one of the core components of the wind-powered system. The vibration of the drivetrain system causes uncertainty in the generator model, which increases the controlling difficulty of the electric drivetrain system. system serving as the basis. (2) The research conclusions of the two-mass syst extended to the multi-mass system. When using the two-mass model, it is ne carry out the necessary gearbox equivalence and conversion of the gearbox's inertia and elastic coefficient.

Composition and Model of Two-Mass Wind-Powered Mechanical Drivetrain System
In addition, the wind turbine's PID control cannot suppress the vibration from gear clearance. Therefore, the electromagnetic torque or speed of the wind can be controlled by the torque-speed closed loop to eliminate the wind drivet tion with clearance. Figure 1 depicts a schematic diagram of the wind drivetrain system comp a doubly fed wind-powered system [24]. The doubly fed generator is connec drivetrain system via a high-speed shaft. Thus, the drivetrain system's perfor fects the generator's output performance. In addition, the generator is connec grid through the transformer, while the generator's electrical control system co converter that regulates the electromagnetic torque and speed of the generator. fed wind generator is one of the core components of the wind-powered syste bration of the drivetrain system causes uncertainty in the generator model, creases the controlling difficulty of the electric drivetrain system. The gearbox in Figure 1 is regarded as an equivalent gear. The mass of th gear is converted to the wind turbine side, and the mass of the driven gear is co the generator side. The gear clearance influence of the three-stage gear system i box is represented by an equivalent total clearance to ensure that a two-mass model with gear clearance can be obtained ( Figure 2). In Figure 2, W is the driving torque of the wind turbine, F is the rotatio of the generator, L is the load resistance torque, W is the wind turbine rotat F is the generator rotation speed, W is the equivalent moment of inertia of turbine mass block (including the converted value of the rotational inertia of t gear), F is the equivalent moment of inertia of the generator mass block (inc converted value of the rotational inertia of the driven gear), and K is the equiva coefficient of the wind drivetrain that includes the converted elastic coefficient o The gearbox in Figure 1 is regarded as an equivalent gear. The mass of the driving gear is converted to the wind turbine side, and the mass of the driven gear is converted to the generator side. The gear clearance influence of the three-stage gear system in the gearbox is represented by an equivalent total clearance to ensure that a two-mass drivetrain model with gear clearance can be obtained ( Figure 2). carry out the necessary gearbox equivalence and conversion of the gearbox's inertia and elastic coefficient.

Composition and Model of Two-Mass Wind-Powered Mechanical Drivetrain Syst
In addition, the wind turbine's PID control cannot suppress the vibratio from gear clearance. Therefore, the electromagnetic torque or speed of the wind can be controlled by the torque-speed closed loop to eliminate the wind drivet tion with clearance. Figure 1 depicts a schematic diagram of the wind drivetrain system comp a doubly fed wind-powered system [24]. The doubly fed generator is conne drivetrain system via a high-speed shaft. Thus, the drivetrain system's perfo fects the generator's output performance. In addition, the generator is conne grid through the transformer, while the generator's electrical control system c converter that regulates the electromagnetic torque and speed of the generator fed wind generator is one of the core components of the wind-powered syste bration of the drivetrain system causes uncertainty in the generator model, creases the controlling difficulty of the electric drivetrain system. The gearbox in Figure 1 is regarded as an equivalent gear. The mass of t gear is converted to the wind turbine side, and the mass of the driven gear is co the generator side. The gear clearance influence of the three-stage gear system box is represented by an equivalent total clearance to ensure that a two-mass model with gear clearance can be obtained ( Figure 2). In Figure 2, W is the driving torque of the wind turbine, F is the rotatio of the generator, L is the load resistance torque, W is the wind turbine rota F is the generator rotation speed, W is the equivalent moment of inertia o turbine mass block (including the converted value of the rotational inertia of t gear), F is the equivalent moment of inertia of the generator mass block (inc converted value of the rotational inertia of the driven gear), and K is the equiva coefficient of the wind drivetrain that includes the converted elastic coefficient box. Lastly, the clearance module N represents the total equivalent clearance o box. In Figure 2, T W is the driving torque of the wind turbine, T F is the rotational torque of the generator, T L is the load resistance torque, ω W is the wind turbine rotation speed, ω F is the generator rotation speed, J W is the equivalent moment of inertia of the wind turbine mass block (including the converted value of the rotational inertia of the driving gear), J F is the equivalent moment of inertia of the generator mass block (including the converted value of the rotational inertia of the driven gear), and K is the equivalent elastic coefficient of the wind drivetrain that includes the converted elastic coefficient of the gearbox. Lastly, the clearance module N represents the total equivalent clearance of the gearbox.

Mathematical Model of the Two-Mass Wind Turbine Drivetrain
The gear clearance is not considered in Figure 2. First, Equation (1) can be obtained. Then, the Laplace transform of Equation (1) is used to obtain Equation (2): ( After obtaining the transfer function structure diagram of the wind drivetrain, the nonlinear clearance unit N is added according to Equation (2). The transfer function structure diagram of the two-mass wind drivetrain with gear clearance is shown in Figure 3a. The input/output characteristics of gear clearance are depicted in Figure 3b, and the mathematical model is shown in Equation (3).
where the input is the position v of the equivalent driving gear, the output is the position q of the driven gear, 2b is the total gear clearance, and K is the characteristic slope.

Mathematical Model of the Two-Mass Wind Turbine Drivetrain
The gear clearance is not considered in Figure 2. First, Equation (1) can be obtained. Then, the Laplace transform of Equation (1) is used to obtain Equation (2): After obtaining the transfer function structure diagram of the wind drivetrain, the nonlinear clearance unit N is added according to Equation (2). The transfer function structure diagram of the two-mass wind drivetrain with gear clearance is shown in Figure 3a. The input/output characteristics of gear clearance are depicted in Figure 3b, and the mathematical model is shown in Equation (3).
where the input is the position v of the equivalent driving gear, the output is the position q of the driven gear, 2b is the total gear clearance, and K is the characteristic slope.

Wind Drivetrain Vibration Analysis and Suppression
In the wind drivetrain system shown in Figure 1, the elasticity of the drive shaft and gear clearance of the gearbox causes the vibration of the drivetrain system and the electromagnetic torque or speed of the generator.

Elastic Vibration Analysis of a Two-Mass Wind Drivetrain
The generator speed for the rotational torque T W is ω W . Thus, Equation (4) can be obtained as follows: According to Equation (4), a resonant point in the drivetrain [25] indicates that vibration exists, and the resonant frequency is ω r = J W +J F J W J F K. In Equation (4), the parameter in  The system output of the same parameters without clearance is subtracted from the output of the two-mass wind drivetrain transfer function based on Figure 3a. Thus, an experimental system capable of verifying clearance vibration is created, as shown in Figure 4a. The speed deviation ∆ω F between the two output speed waveforms, which represent the vibration generated by the nonlinear gear clearance unit, is shown in Figure 4b.
According to Equation (4), a resonant point in the drivetrain [25] indicates that vibr tion exists, and the resonant frequency is r = √ W + F W F . In Equation (4), the paramet in Table 1 [5] is substituted to obtain the resonant frequency as r = 0.863 ( rad s ), whi is a relatively low vibration frequency. The system output of the same parameters without clearance is subtracted from the ou put of the two-mass wind drivetrain transfer function based on Figure 3a. Thus, an expe mental system capable of verifying clearance vibration is created, as shown in Figure 4a. Th speed deviation ∆ F between the two output speed waveforms, which represent the v bration generated by the nonlinear gear clearance unit, is shown in Figure 4b.

Vibration Suppression Strategy
The vibration frequency caused by the elasticity of the drivetrain shaft is relative low. As such, it can be suppressed by the wind turbine's PID speed closed-loop contr and the variable pitch of the wind power system. The speed wave of the generator aft suppression is denoted as Δω F1 , as shown in Figure 5a. However, the wind turbine's PI control cannot suppress the vibration resulting from gear clearance ∆ F2 , as shown Figure 5b. Furthermore, the vibration frequency caused by gear clearance is relative high [2,3,26]. If variable pitch control is employed, then the variable pitch actuator cann

Vibration Suppression Strategy
The vibration frequency caused by the elasticity of the drivetrain shaft is relatively low. As such, it can be suppressed by the wind turbine's PID speed closed-loop control and the variable pitch of the wind power system. The speed wave of the generator after suppression is denoted as ∆ω F1 , as shown in Figure 5a. However, the wind turbine's PID control cannot suppress the vibration resulting from gear clearance ∆ω F2 , as shown in Figure 5b. Furthermore, the vibration frequency caused by gear clearance is relatively high [2,3,26]. If variable pitch control is employed, then the variable pitch actuator cannot follow the controller's signal to change rapidly. Therefore, the electromagnetic torque or speed of the wind generator can be controlled by the torque-speed closed loop to eliminate the wind drivetrain vibration. follow the controller's signal to change rapidly. Therefore, the electromagnetic torque or speed of the wind generator can be controlled by the torque-speed closed loop to eliminate the wind drivetrain vibration.

Velocity Vector Control of a Doubly Fed Wind Generator
A doubly fed induction generator (DFIG) is a nonlinear multivariable system with a high order and strong coupling. The conventional control method is not only complicated but also ineffective. Vector control can simplify the control method and obtain the same control effect as a DC motor in some aspects. At present, the research on vector control in an AC motor has made a lot of achievements and has been widely used in production, but the application of vector control technology in a generator is still in the research stage [24]. The vector control of an AC motor mainly adopts the rotor flux orientation, and the vector control of a DFIG adopts the stator flux orientation.

The Mathematical Model of a Doubly Fed Wind Generator Vector Control Based on the Stator Flux Orientation
A doubly fed generator usually adopts a vector control method to achieve the torque component of the current and excitation component decoupling control of a doubly fed generator. An asynchronous rotating d-q coordinate system model based on the stator flux orientation can be used while assuming a negligible effect of the doubly fed generator stator resistance. Under this condition, the voltage and current of a doubly fed generator can be derived as follows: where Δ is the deviation of the angular frequency of a doubly fed generator, i.e., Δ = − , is the rotating magnetic field that synchronizes the angular frequency, and is the angular frequency of the rotor speed. The parameters and are the stator and rotor resistances, respectively, and , , and are the stator inductance, rotor self-inductance, and mutual inductance between stator and rotors, respectively. The parameters 1 , 1 , 2 , and 2 are the d-and q-axis components of the stator and rotor voltages; 1 , 1 , 2 , 2 are the d-and q-axis components of the stator and rotor currents; 1 , 1 , 2 , 2 are the d-and q-axis components of the stator and rotor flux; p is the differentiation operator.

Velocity Vector Control of a Doubly Fed Wind Generator
A doubly fed induction generator (DFIG) is a nonlinear multivariable system with a high order and strong coupling. The conventional control method is not only complicated but also ineffective. Vector control can simplify the control method and obtain the same control effect as a DC motor in some aspects. At present, the research on vector control in an AC motor has made a lot of achievements and has been widely used in production, but the application of vector control technology in a generator is still in the research stage [24]. The vector control of an AC motor mainly adopts the rotor flux orientation, and the vector control of a DFIG adopts the stator flux orientation.

The Mathematical Model of a Doubly Fed Wind Generator Vector Control Based on the Stator Flux Orientation
A doubly fed generator usually adopts a vector control method to achieve the torque component of the current and excitation component decoupling control of a doubly fed generator. An asynchronous rotating d-q coordinate system model based on the stator flux orientation can be used while assuming a negligible effect of the doubly fed generator stator resistance. Under this condition, the voltage and current of a doubly fed generator can be derived as follows: where ∆ω r is the deviation of the angular frequency of a doubly fed generator, i.e., ∆ω r = ω s − ω r , ω s is the rotating magnetic field that synchronizes the angular frequency, and ω r is the angular frequency of the rotor speed. The parameters R s and R r are the stator and rotor resistances, respectively, and L s , L r , and L m are the stator inductance, rotor self-inductance, and mutual inductance between stator and rotors, respectively. The parameters u d1 , u q1 , u d2 , and u q2 are the dand qaxis components of the stator and rotor voltages; i d1 , i q1 , i d2 , i q2 are the dand q-axis components of the stator and rotor currents; ψ d1 , ψ q1 , ψ d2 , ψ q2 are the dand q-axis components of the stator and rotor flux; p is the differentiation operator. After transformation, the state equation of the rotor voltage can finally be obtained as follows: where σ is the total magnetic leakage coefficient equal to σ = 1 − L 2 m /L s L r , and i m1 is the stator excitation current.
In Equation (6), the second term on the left of the equal sign of the two equations is the constant related to the stator flux that can be regarded as a disturbance of the system and is compensated by the controller. The compensations terms added to the controller are is the total magnetic leakage coefficient equal to = 1 − 2 ∕ , a the stator excitation current.
In Equation (6) Rotor Model Figure 6. Dynamic structure diagram of DFIG in synchronous rotation d-q coordinate s A cross-coupling term between the d-axis and q-axis components of the ro in the generator can be seen in Figure 6. This indicates that the two components independently adjusted and affects the electromagnetic torque and speed perfo the generator. Hence, it is necessary to decouple the cross-coupling of the wind

IMC Decoupling of Cross-Coupling Voltage of a Doubly Fed Wind Generator
The common decoupling methods of cross-coupled voltage include volta ward decoupling, diagonal matrix decoupling, and voltage feedback decouplin coupling effect of the methods mentioned above depends on the accurate ma model of the controlled object and the stability of the parameters. The IMC does depend on the accurate mathematical model of the controlled object. Moreove has few model accuracy requirements. Thus, it can be easily achieved in engi this section, the conventional IMC decoupling method is adopted. The structu that implements this decoupling control is shown in Figure 7. A cross-coupling term between the d-axis and q-axis components of the rotor current in the generator can be seen in Figure 6. This indicates that the two components cannot be independently adjusted and affects the electromagnetic torque and speed performance of the generator. Hence, it is necessary to decouple the cross-coupling of the wind generator.

IMC Decoupling of Cross-Coupling Voltage of a Doubly Fed Wind Generator
The common decoupling methods of cross-coupled voltage include voltage feedforward decoupling, diagonal matrix decoupling, and voltage feedback decoupling. The decoupling effect of the methods mentioned above depends on the accurate mathematical model of the controlled object and the stability of the parameters. The IMC does not overly depend on the accurate mathematical model of the controlled object. Moreover, the IMC has few model accuracy requirements. Thus, it can be easily achieved in engineering. In this section, the conventional IMC decoupling method is adopted. The structure diagram that implements this decoupling control is shown in Figure 7.
When the controlled object matches the forecast model, i.e., G(s) =Ĝ(s), the two components of the rotor current are decoupled [20].
The following parameters of the wind generator required for simulation are obtained from Table 7-3 in [24]. The rotor resistance R r = 0.1442 Ω, the rotor inductance L r = 35.7 mH, and the calculated parameters ∆ω r σL r = 0.35 Ω. The simulation model was built according to Figure 7. The provided input value i * d2 jumped from 0 A to 5.1 A at t = 1 s, and i * q2 from 0 A to 5.1 A at t = 1.5 s. The following can be proven: (1) When the inductance and resistance parameters of the wind generator rotor increase, decrease or remain constant in the same proportion, the cross-coupled voltage can achieve full decoupling by the 1-DOF-IMC, and the corresponding decoupling effect is good. (2) When the wind generator's rotor inductance and resistance parameters change in different proportions and reverse directions, the cross-coupled voltage cannot achieve full decoupling by the 1-DOF-IMC due to few adjustable parameters. However, the cross-coupled voltage can achieve full decoupling by the 3-DOF-IMC. Limited by the length of the paper, the corresponding discussion will not be carried out.
The following parameters of the wind generator required for simulation ar from Table 7-3 in [24]. The rotor resistance = 0.1442 Ω, the rotor inductance mH, and the calculated parameters ∆ = 0.35 Ω. The simulation model w cording to Figure 7. The provided input value 2 * jumped from 0 A to 5.1 A at 2 * from 0 A to 5.1 A at t = 1.5 s. The following can be proven: (1) When the inductance and resistance parameters of the wind generator roto decrease or remain constant in the same proportion, the cross-coupled v achieve full decoupling by the 1-DOF-IMC, and the corresponding decoup is good.
(2) When the wind generator's rotor inductance and resistance parameters different proportions and reverse directions, the cross-coupled volta achieve full decoupling by the 1-DOF-IMC due to few adjustable parame ever, the cross-coupled voltage can achieve full decoupling by the 3-DOFited by the length of the paper, the corresponding discussion will not be c

Design of Control System for a Doubly Fed Generator
The electric control system of a doubly fed wind generator usually adop and current double closed-loop structure. The inner current loop is designe lowed by the outer speed loop. Generally, the inner current loop is designed as order optimal system" of the EDM to improve its dynamic performance. The 3

Design of Control System for a Doubly Fed Generator
The electric control system of a doubly fed wind generator usually adopts a speed and current double closed-loop structure. The inner current loop is designed first, followed by the outer speed loop. Generally, the inner current loop is designed as a "second-order optimal system" of the EDM to improve its dynamic performance. The 3-DOF-IMC method is used to design various controllers for the speed loop that represent the performance of the speed-regulating system. Thus, the speed loop's tracking, anti-interference, and robustness performances can be flexibly adjusted.
In addition, the EDM is also used to design the speed loop, and it is compared with the system designed by the 3-DOF-IMC.

Design of the Rotor's Current Control System for a Doubly Fed Generator
According to the rotor's current control block diagram shown in Figure 7, the dand q-axis channels of the inner current loop are completely independent after decoupling and have the same structure. Therefore, the two current controllers are the same, and only a single controller needs to be designed. The simplified dynamic structure of the rotor current loop of a doubly fed generator is shown in Figure 8. Equation (7) represents the current feedback filtering unit, β is the current feedback coefficient, and T oi is the filtering time constant. Usually, a filter unit is also added to a given part of the current.  The current loop is designed according to the "second-order optimal system" method, and the obtained current controller ACR is a PI controller: The transfer function of the current loop is approximated to an inertial unit with a small time constant via the closed-loop control of the rotor current: Ws Ts An important function of inner current loop control is transforming the controlled object and accelerating the current tracking.

Design of Rotor Speed Control System for a Doubly Fed Generator
The reactive power control channel usually adopts a single closed-loop control to ensure that the wind generator captures the maximum wind energy. This is achieved by controlling the current of the d-axis of the rotor within the d-q coordinate system of the synchronous rotation of two phases. The speed control channel adopts a double closedloop control. The inner loop is the current loop of the doubly fed generator rotor, which is achieved by controlling the current of the q-axis of the rotor. The outer loop is the generator speed loop. The block diagram of a doubly fed generator's complete rotor electrical control system is shown in Figure 9. (1) The controlled object of the speed loop The simplified rotor inner current loop, i.e., Equation (9), is substituted into Figure 9. After further simplification, the dynamic structure diagram of the speed loop is shown in Figure 10. The current loop is designed according to the "second-order optimal system" method, and the obtained current controller ACR is a PI controller: The transfer function of the current loop is approximated to an inertial unit with a small time constant via the closed-loop control of the rotor current: An important function of inner current loop control is transforming the controlled object and accelerating the current tracking.

Design of Rotor Speed Control System for a Doubly Fed Generator
The reactive power control channel usually adopts a single closed-loop control to ensure that the wind generator captures the maximum wind energy. This is achieved by controlling the current of the d-axis of the rotor within the d-q coordinate system of the synchronous rotation of two phases. The speed control channel adopts a double closedloop control. The inner loop is the current loop of the doubly fed generator rotor, which is achieved by controlling the current of the q-axis of the rotor. The outer loop is the generator speed loop. The block diagram of a doubly fed generator's complete rotor electrical control system is shown in Figure 9. The current loop is designed according to the "second-order optima method, and the obtained current controller ACR is a PI controller: The transfer function of the current loop is approximated to an inertial u small time constant via the closed-loop control of the rotor current: An important function of inner current loop control is transforming the object and accelerating the current tracking.

Design of Rotor Speed Control System for a Doubly Fed Generator
The reactive power control channel usually adopts a single closed-loop ensure that the wind generator captures the maximum wind energy. This is ac controlling the current of the d-axis of the rotor within the d-q coordinate sys synchronous rotation of two phases. The speed control channel adopts a doub loop control. The inner loop is the current loop of the doubly fed generator ro is achieved by controlling the current of the q-axis of the rotor. The outer loop erator speed loop. The block diagram of a doubly fed generator's complete roto control system is shown in Figure 9. (1) The controlled object of the speed loop The simplified rotor inner current loop, i.e., Equation (9), is substituted int After further simplification, the dynamic structure diagram of the speed loop is Figure 10. (1) The controlled object of the speed loop The simplified rotor inner current loop, i.e., Equation (9), is substituted into Figure 9. After further simplification, the dynamic structure diagram of the speed loop is shown in Figure 10. The output speed of the wind generator requires good tracking and anti-in performances; it is also robust to the changes in system parameters. Most of controllers designed by the EDM are PI controllers, making it difficult to simu obtain good tracking and anti-interference performance. In this paper, a 3 method is used to design the speed loop, and three controllers are used to adjus ing performance, anti-interference performance, and robustness of the speed lo

The Principle of a 3-DOF-IMC
(1) Decoupling effect: In the early stage, the 3-DOF-IMC was mainly used for pling of a multivariable process control system, and in the late 1990s, it w to the decoupling of a motor model [27]. With the development of power and power grid technology, the method of IMC decoupling has also pene this field. For example, the cross-coupling problem of active and reactiv solved in the PWM rectifier. Active and reactive currents are decoupled DC transmission systems, stationary synchronous series compensators a source converters. The IMC decoupling method also has many applicatio wind power generation system. The principle of the conventional IMC (1decoupling is discussed in the reference, and the methods of the 2-DOF-I DOF-IMC decoupling and their applications in DFIG vector control are pr (2) The multiple degrees of freedom IMC is used to achieve the independent of multiple performance indexes of the control system.
The performance index of the control system includes tracking perform interference performance and robustness performance. The controller designe DOF-IMC is difficult to achieve excellent tracking performance and anti-interfe formance at the same time, which cannot meet the control requirements of a h mance system. The 2-DOF-IMC can adjust the tracking performance and anti-in performance of the system, respectively, and can mainly control the anti-interfe formance of the system. However, when 2-DOF-IMC controllers are used to re tracking performance and anti-interference performance, respectively, then th ward IMC controller needs to take both the anti-interference and robustness in In the 3-DOF-IMC, the given filter regulates the tracking performance of the s The declaring constants can be expressed as follows: where J F is the generator's moment of inertia and p m is the pole number of a doubly fed generator. The output speed of the wind generator requires good tracking and anti-interference performances; it is also robust to the changes in system parameters. Most of the speed controllers designed by the EDM are PI controllers, making it difficult to simultaneously obtain good tracking and anti-interference performance. In this paper, a 3-DOF-IMC method is used to design the speed loop, and three controllers are used to adjust the tracking performance, anti-interference performance, and robustness of the speed loop.

The Principle of a 3-DOF-IMC
(1) Decoupling effect: In the early stage, the 3-DOF-IMC was mainly used for the decoupling of a multivariable process control system, and in the late 1990s, it was applied to the decoupling of a motor model [27]. With the development of power electronics and power grid technology, the method of IMC decoupling has also penetrated into this field. For example, the cross-coupling problem of active and reactive power is solved in the PWM rectifier. Active and reactive currents are decoupled in flexible DC transmission systems, stationary synchronous series compensators and voltage source converters. The IMC decoupling method also has many application uses in a wind power generation system. The principle of the conventional IMC (1-DOF-IMC) decoupling is discussed in the reference, and the methods of the 2-DOF-IMC and 3-DOF-IMC decoupling and their applications in DFIG vector control are presented. (2) The multiple degrees of freedom IMC is used to achieve the independent regulation of multiple performance indexes of the control system.
The performance index of the control system includes tracking performance, antiinterference performance and robustness performance. The controller designed by the 1-DOF-IMC is difficult to achieve excellent tracking performance and anti-interference performance at the same time, which cannot meet the control requirements of a high-performance system. The 2-DOF-IMC can adjust the tracking performance and anti-interference performance of the system, respectively, and can mainly control the anti-interference performance of the system. However, when 2-DOF-IMC controllers are used to regulate the tracking performance and anti-interference performance, respectively, then the feedforward IMC controller needs to take both the anti-interference and robustness into account. In the 3-DOF-IMC, the given filter regulates the tracking performance of the system, the feedforward IMC controller mainly regulates the anti-interference performance, and the feedback filter regulates the robustness performance. From the research status, the 3-DOF-IMC is mainly used in the control system performance index regulation and dynamic decoupling. In addition to analyzing the tracking performance and anti-interference performance of the 3-DOF-IMC, this section also tries to use it to suppress the vibration of the drive train.

The Structure of a 3-DOF-IMC
The structure diagram of a 3-DOF-IMC is shown in Figure 11 [23].

The Structure of a 3-DOF-IMC
The structure diagram of a 3-DOF-IMC is shown in Figure 11 [23]. Figure 11. Block diagram of a 3-DOC-IMC.
In Figure 11, IMC I (s), IMC II (s), and f ( ) constitute a 3-DOF-IMC f ( ) that regulates the system's robustness, f ( ), which is also known as a feedback filter. When f ( ) = 1, the system degenerates to a 2-DOF-IMC with IMC II (s) = 1. Then, the system is reduced to a 1-DOF-IMC. The following expression can be obtained according to Figure  11: When the controlled object matches the forecast model, i.e., ( ) = ̂( s): The 1-DOF-IMC method is used to design IMC I (s), where

The Analysis of the Tracking Performance and Anti-Interference Performance of a 3-DOF-IMC
According to Equation (12), the following expression can be obtained: The tracking and anti-interference performances of the system can be independently adjusted by changing the parameters of 1 ( ), 2 ( ), and f ( ).

The Robustness Analysis of a 3-DOF-IMC
The mechanical vibration of the wind drivetrain system will cause uncertainty in the generator model, which conforms to the multiplicative description of Equation (15) [22].
Usually, the controlled object containing uncertainty can be described in one of the following two ways: (1) Additive description: (2) Multiplicative description: Figure 11. Block diagram of a 3-DOC-IMC.
In Figure 11, C I IMC (s), C II IMC (s), and F f (s) constitute a 3-DOF-IMC F f (s) that regulates the system's robustness, F f (s), which is also known as a feedback filter. When F f (s) = 1, the system degenerates to a 2-DOF-IMC with C II IMC (s) = 1. Then, the system is reduced to a 1-DOF-IMC. The following expression can be obtained according to Figure 11: When the controlled object matches the forecast model, i.e., G(s) =Ĝ(s): The 1-DOF-IMC method is used to design C I IMC (s), where C I IMC (s) =Ĝ −1 − (s)L 1 (s), and controller C II IMC (s) = L 2 (s) L 1 (s) is selected. Since the decoupled system is approximately a first-order system, it is desirable that F f (s) = αs+1 βs+1 . In the above correlation formula, α and β represent filter parameters of the filter F f (s). The internal model of the controlled object isĜ(s) =Ĝ + (s)Ĝ − (s).

The Analysis of the Tracking Performance and Anti-Interference Performance of a 3-DOF-IMC
According to Equation (12), the following expression can be obtained: The tracking and anti-interference performances of the system can be independently adjusted by changing the parameters of L 1 (s), L 2 (s), and F f (s).

The Robustness Analysis of a 3-DOF-IMC
The mechanical vibration of the wind drivetrain system will cause uncertainty in the generator model, which conforms to the multiplicative description of Equation (15) [22].
Usually, the controlled object containing uncertainty can be described in one of the following two ways: (1) Additive description: (2) Multiplicative description: Equations (14) and (15) are equal to each other, thus obtaining the following equation: where ∆(jω) ∞ ≤ 1 and W(s) is the weighting function that expresses the degree to which the uncertainty of the controlled object depends on the frequency ω. The controlled object G(s) in Figure 11 is expressed in the form of Equation (15), and its structure is transformed to obtain the equivalent structure in the dashed box, as shown in Figure 12. In this section, the robustness of an IMC system is mainly discussed. Therefore, only the influence of model uncertainty on the system characteristics is considered, while the influence of disturbance D(s) on the system performance is temporarily neglected.
where ( ) 1 j    and ( ) is the weighting function that expresses the deg which the uncertainty of the controlled object depends on the frequency ω.
The controlled object ( ) in Figure 11 is expressed in the form of Equatio and its structure is transformed to obtain the equivalent structure in the dashed b shown in Figure 12. In this section, the robustness of an IMC system is mainly disc Therefore, only the influence of model uncertainty on the system characteristics is c ered, while the influence of disturbance ( ) on the system performance is temp neglected. The following relationships can be obtained according to Figures 11 and 12: I  I  I  IMC  1  f  IMC  1  IMC  f   I  I  I  IMC 1 Figure 11 can be equalized to Figure 12. Thus, the 3-DOF-IMC is an equivalent loop system.

The Robustness Stability Analysis of the 3-DOF-IMC
According to the IMC requirements, the parameters IMC II (s) and ( ) are sta shown in Figure 12. Therefore, the necessary and sufficient condition for the stab the open-loop control system shown in Figure 12 is that the local feedback unit the figure is stable. Furthermore, the necessary and sufficient conditions for ( ) ity can be expressed as follows: Equation (16) can be substituted into Equation (22) to obtain: The following relationships can be obtained according to Figures 11 and 12: Figure 11 can be equalized to Figure 12. Thus, the 3-DOF-IMC is an equivalent open-loop system.

The Robustness Stability Analysis of the 3-DOF-IMC
According to the IMC requirements, the parameters C II IMC (s) and G(s) are stable, as shown in Figure 12. Therefore, the necessary and sufficient condition for the stability of the open-loop control system shown in Figure 12 is that the local feedback unit P(s) in the figure is stable. Furthermore, the necessary and sufficient conditions for P(s) stability can be expressed as follows: Equation (16) can be substituted into Equation (22) to obtain: According to Figure 12 and Equation (22), the IMC system is an open-loop control. The zero poles of C II IMC (s) and G(s) can exactly counteract to achieve the system output tracking input given value. When uncertainty exists in the object model, the internal model of the object is introduced asĜ(s). The model uncertainty is separated from the object model. Then, the local feedback unit P(s) is used. The feedback unit compensates to achieve robust control of the system.
The actual object may have a variety of structures. Regardless of the structure, the robust stability of the 3-DOF-IMC system is related to the feedback filter related to F f (s). Due to an import filter F f (s), C I IMC (s)F f (s)∆G(s) ∞ is lower than C I IMC (s)∆G(s) ∞ . Therefore, the degree of uncertainty ∆G(s) ∞ that the system will allow is higher than the disturbance from the outer environment.

The Robustness Performance Analysis of a 3-DOF-IMC
The system should be robust, stable, and characterized by good dynamic performance in engineering. The system should completely overcome the influence of the model uncertainty; hence, the following equation should be established according to Equation (19): Substituting Equation (20) into Equation (23) yields: The sensitivity function S = 1−F f (s).C I IMC (s).Ĝ(s) 1+C I IMC (s).F f (s).∆G(s) is then defined. The sensitivity function S's value reflects the system's robust tracking performance. The smaller S is, the better the robust tracking performance.
Similarly, regardless of the structure of the controlled object, the robustness of the 3-DOF-IMC system is also related to the feedback filter F f (s), whose introduction decreases the sensitivity S and improves the robustness of the system.

The 3-DOF-IMC of Wind Generator Speed Loop
In this section, the MATLAB (2019b) simulation method is used to verify the tracking performance, anti-interference performance and vibration suppression ability of the IMC to control the DFIG. Compared with the system designed based on the EDM, the experiment proves the superiority of the 3-DOF-IMC method.
The dynamic structure diagram of the generator speed loop based on the 3-DOF-IMC is shown in Figure 13. The following parameters comprise Figure 13: The system should be robust, stable, and characterized by good dynam mance in engineering. The system should completely overcome the influence o uncertainty; hence, the following equation should be established according t (19): Substituting Equation (20) into Equation (23) yields: The sensitivity function is then defined. The function S's value reflects the system's robust tracking performance. The smal better the robust tracking performance. Similarly, regardless of the structure of the controlled object, the robustne DOF-IMC system is also related to the feedback filter f ( ), whose introductio the sensitivity S and improves the robustness of the system.

The 3-DOF-IMC of Wind Generator Speed Loop
In this section, the MATLAB (2019b) simulation method is used to verify t performance, anti-interference performance and vibration suppression ability to control the DFIG. Compared with the system designed based on the EDM, ment proves the superiority of the 3-DOF-IMC method.
The dynamic structure diagram of the generator speed loop based on the 3 is shown in Figure 13. The following parameters comprise  (1) The controlled object: (2) Internal model controller C I IMC (s) designed via the 1-DOF-IMC method: L 1 (s) is introduced for the physical realization of the controller C I IMC (s).
The system output can be obtained according to Equation (13): Since the controlled object contains no right-half plane zero,Ĝ + (s) = 1. (2) Simulation parameters The parameters of the doubly fed wind generator used in the simulation are as follows: the moment of inertia J F is 0.776 kg·m 2 , the polar logarithm P m is 3, the constant K Ψ /β is 0.388, and the filter time constant T oi = 0.001 s.
(3) Model of the controlled object and the controller The mathematical model of the controlled object can be expressed as:

(4) Simulation waveform and conclusion
The given input is set to 40 × 1 (t), the load disturbance input is 5 × 1 (t−15), and the 3-DOF-IMC is adopted. The ratio of α/β = 1 in the filter is kept constant, and three parameters (a), (b), and (c) in Table 2 are selected to obtain the tracking performance and anti-interference performance waveforms of the generator output speed, as shown in Figure 14a-c. Table 2. Selection schemes and conclusions for the three parameters in Figure 14.

Conclusion
When the parameter λ n1 is constant, the anti-interference performance is also constant The tracking performance is improved with a decrease in parameter λ n2 Anti-interference performance The tracking performance The close-up of the anti-interference performance curve is shown in Figure 14. When the parameter n1 is constant, the anti-interference performance is also constant, while the tracking performance is improved with a decrease in the parameter n2 .  The close-up of the anti-interference performance curve is shown in Figure 14. When the parameter λ n1 is constant, the anti-interference performance is also constant, while the tracking performance is improved with a decrease in the parameter λ n2 .
6.2. The Anti-Interference Performance of the Wind Generator Based on C I I MC (s) Three parameter-selection schemes (a), (b), and (c) in Table 3 are taken for α/β = 1 within the filter. It can be concluded that the waveform of scheme (a) is the same as that shown in Figure 14b. The tracking performance and anti-interference waveforms of the generator speed corresponding to schemes (b) and (c) are shown in Figure 15a,b. Table 3. Selection schemes and conclusions for the three parameters in Figure 15. The close-up of the anti-interference performance curve is shown in Figure 14. When the parameter n1 is constant, the anti-interference performance is also constant, while the tracking performance is improved with a decrease in the parameter n2 .

The Anti-Interference Performance of the Wind Generator Based on ( )
Three parameter-selection schemes (a), (b), and (c) in Table 3 are taken for α/β = 1 within the filter. It can be concluded that the waveform of scheme (a) is the same as tha shown in Figure 14b. The tracking performance and anti-interference waveforms of the generator speed corresponding to schemes (b) and (c) are shown in Figure 15a,b. Table 3. Selection schemes and conclusions for the three parameters in Figure 15.  The close-up of the anti-interference performance curve is shown in Figure 15. It can be concluded that the anti-interference performance of the speed loop can be adjusted by modifying the parameter λ n1 . The parameters are small, and the anti-interference performance is good. When the parameter λ n2 is constant, the tracking performance is unchanged.

The Anti-Interference Performance of the Generator Based on the Feedback Filter F f (s)
According to Y(s) = L 2 (s) · R(s) + [1 − F f (s) · L 1 (s)] · D(s), the parameter F f (s) can also adjust the anti-interference performance of the generator speed with respect to the disturbance.
(1) Anti-interference performance analysis of generator speed against step disturbance when parameters of F f (s) are changed. The provided input R(t) = 40 × 1 (t) and the disturbance input D(t) = −5 × 1(t − 20) are set, and a 3-DOF-IMC is adopted. (1) Two schemes (a) and (b) in Table 4 are employed to obtain the tracking performance and anti-interference performance waveforms of generator speed, as shown in Figure 16a,b. The close-up of the disturbance point waveform is indicated within the figure.
(2) According to the anti-interference performance analysis, it can be concluded that the tracking performance of the speed loop is the same when the parameter λ n2 is unchanged. For the same parameter λ n1 , the maximum dynamic descent of the anti-interference performance remains constant. The recovery time varies with the value of α/β. The closer the value of α/β is to 1, the better the anti-interference performance. (2) Analysis of the anti-interference performance of generator speed to periodic disturbance signal when F f (s) is changed.
The following inputs are employed: R(t) = 40*1 (t), and a load disturbance input in a periodic sine wave signal equal to D(t) = 8sinωt. Lastly, a 3-DOF-IMC is also employed. (1) Two schemes (a) and (b) in Table 5 are considered to obtain the anti-interference performance waveforms of the generator speed to periodic disturbance signals, as shown in Figure 17a,b.
Anti-interference performance analysis: It can be concluded that different α/β values of the filter can affect the disturbance suppression ability of periodic disturbance signal. The closer the value of α/β is to 1, the stronger the ability of the system to suppress periodic disturbances. Table 4. Two-step disturbance schemes and corresponding conclusions.

No Tracking Performance Anti-Interference Performance
Conclusion λ n2 is unchanged, tracking performance is unchanged λ n1 is unchanged, the maximum dynamic descent of anti-interference performance is unchanged The recovery time of anti-interference performance varies with the value of α/β Table 5. Two schemes and conclusions of periodic perturbation with sine wave signal.

2022, 15, x FOR PEER REVIEW 17 of 26
The close-up of the anti-interference performance curve is shown in Figure 15. It can be concluded that the anti-interference performance of the speed loop can be adjusted by modifying the parameter n1 . The parameters are small, and the anti-interference performance is good. When the parameter n2 is constant, the tracking performance is unchanged.

The Anti-Interference Performance of the Generator Based on the Feedback Filter ( )
According to , the parameter f ( ) can also adjust the anti-interference performance of the generator speed with respect to the disturbance.
(1) Anti-interference performance analysis of generator speed against step disturbance when parameters of f ( ) are changed.
(1) Two schemes (a) and (b) in Table 4 are employed to obtain the tracking performance and anti-interference performance waveforms of generator speed, as shown in Figure 16a,b. The close-up of the disturbance point waveform is indicated within the figure.  (2) According to the anti-interference performance analysis, it can be concluded that the tracking performance of the speed loop is the same when the parameter n2 is unchanged. For the same parameter n1 , the maximum dynamic descent of the anti-interference performance remains constant. The recovery time varies with the value of α/β. The closer the value of α/β is to 1, the better the anti-inter-
(a) (b) (2) Anti-interference performance analysis: It can be concluded that different α/β values of the filter can affect the disturbance suppression ability of periodic disturbance signal. The closer the value of α/β is to 1, the stronger the ability of the system to suppress periodic disturbances.

Suppression of Drivetrain Vibration by Generator Speed Control Based on 3-DOF-IMC
Wind generator speed control based on the 3-DOF-IMC can suppress periodic disturbances, as demonstrated in Section 6.3. Moreover, generator speed vibration can be suppressed by adjusting the parameters of the filter f ( ). According to the analysis in Sections 2.3 and 2.4, the vibration waveform of the generator in the wind mechanical drivetrain system is also periodic. Therefore, the influence of the drivetrain vibration on the generator speed fluctuation can be eliminated by a 3-DOF-IMC of the wind generator.

Experimental Analysis of Drivetrain Vibration Suppression Based on the 3-DOF-IMC of Generator
The experimental simulation model of suppressing drivetrain vibration based on the 3-DOF-IMC is shown in Figure 18.

Suppression of Drivetrain Vibration by Generator Speed Control Based on 3-DOF-IMC
Wind generator speed control based on the 3-DOF-IMC can suppress periodic disturbances, as demonstrated in Section 6.3. Moreover, generator speed vibration can be suppressed by adjusting the parameters of the filter F f (s). According to the analysis in Sections 2.3 and 2.4, the vibration waveform of the generator in the wind mechanical drivetrain system is also periodic. Therefore, the influence of the drivetrain vibration on the generator speed fluctuation can be eliminated by a 3-DOF-IMC of the wind generator.

Experimental Analysis of Drivetrain Vibration Suppression Based on the 3-DOF-IMC of Generator
The experimental simulation model of suppressing drivetrain vibration based on the 3-DOF-IMC is shown in Figure 18. Parameters of the two-mass drivetrain are taken from Table 3 in [5]. The 3-DOF-IMC generator parameters are described in this section. Furthermore, the effect of the f ( ) parameters on the system is discussed.
(1) Analysis of the suppression of generator speed vibration for a variable f ( ) and constant α = β ratio.
The input is set to ( ) = 40×1 (t) and the 3-DOF-IMC is adopted. Two schemes (a) and (b) in Table 6 are taken to obtain the anti-interference performance waveforms of the generator speed to drivetrain vibration signals, as shown in Figure 19a,b. In the table, n1 = 0.08, n2 = 0.5, and α/β = 1. However, the values of α or β are different. The vibration waveform of the two-mass drivetrain is shown in Figure 19a. The generator speed waveform with drivetrain vibration input signal is shown in Figure 19b. Table 6. Two schemes and conclusions of periodic perturbation.

No
Tracking Performance Anti-Interference Performance α/β Parameters of the two-mass drivetrain are taken from Table 3 in [5]. The 3-DOF-IMC generator parameters are described in this section. Furthermore, the effect of the F f (s) parameters on the system is discussed.
(1) Analysis of the suppression of generator speed vibration for a variable F f (s) and constant α = β ratio. The input is set to R(t) = 40×1 (t) and the 3-DOF-IMC is adopted. Two schemes (a) and (b) in Table 6 are taken to obtain the anti-interference performance waveforms of the generator speed to drivetrain vibration signals, as shown in Figure 19a,b. In the table, λ n1 = 0.08, λ n2 = 0.5, and α/β = 1. However, the values of α or β are different. The vibration waveform of the two-mass drivetrain is shown in Figure 19a. The generator speed waveform with drivetrain vibration input signal is shown in Figure 19b.
(1) According to the anti-interference performance analysis, it can be concluded that the optimal anti-interference performance of the drivetrain is obtained for the value of α/β equal to 1. Here, the output speed of the generator is not characterized by vibration.
(2) Analysis of the suppression of generator speed vibration for a varying F f (s) parameter with respect to α = β.
(1) Two schemes (a) and (b) in Table 7 are taken to obtain the anti-interference performance waveforms of the generator speed to drivetrain vibration signals, as shown in Figure 20a,b. In the table, λ n1 = 0.08, λ n2 = 0.5, and the α/β values are mutually different. The vibration waveform of the two-mass drivetrain remains constant, as shown in Figure 19a. The generator speed waveform with the drivetrain vibration signal is shown in Figure 20a,b.
According to the anti-interference performance analysis, it can be concluded that the disturbance suppression ability of F f (s) of the periodic disturbance signal differs when α = β. The closer the value of α/β is to 1, the stronger the ability of the system to suppress periodic disturbances.
(3) Analysis of generator speed vibration for the constant α/β ratio and modified values of α and β.
(1) Three schemes (a), (b), and (c) in Table 8 are considered to obtain the antiinterference performance waveforms of the generator speed to the drivetrain vibration signals, as shown in Figure 21a-c. The value of α/β is different for λ n1 = 0.08 and λ n2 = 0.5 in the table. However, the specific values of α and β are also different. The vibration waveform of the two-mass drivetrain remains constant. The waveform of the generator speed interference of the periodic vibration of the drivetrain is shown in Figure 21a-c. Table 6. Two schemes and conclusions of periodic perturbation.

Conclusion
The value of α/β is equal to 1, and good suppression of periodic disturbance is observed. The outputs of schemes (a) and (b) are the same. Table 7. Two schemes and conclusions of periodic vibration of drivetrain.

No. Tracking Performance Anti-Interference Performance
The closer the value of α/β is to 1, the better the suppression performance of the periodic disturbance. Table 8. Two schemes of periodic perturbation of drivetrain vibration with corresponding conclusions. (1) Analysis of the suppression of generator speed vibration for a variable f ( ) and constant α = β ratio.

No. Tracking Performance Anti-Interference Performance
The input is set to ( ) = 40×1 (t) and the 3-DOF-IMC is adopted. Two schemes (a) and (b) in Table 6 are taken to obtain the anti-interference performance waveforms of the generator speed to drivetrain vibration signals, as shown in Figure 19a,b. In the table, n1 = 0.08, n2 = 0.5, and α/β = 1. However, the values of α or β are different. The vibration waveform of the two-mass drivetrain is shown in Figure 19a. The generator speed waveform with drivetrain vibration input signal is shown in Figure 19b. Table 6. Two schemes and conclusions of periodic perturbation.

Conclusion
The value of α/β is equal to 1, and good suppression of periodic disturbance is observed. The outputs of schemes (a) and (b) are the same. (1) According to the anti-interference performance analysis, it can be concluded that the optimal anti-interference performance of the drivetrain is obtained for the value of α/β equal to 1. Here, the output speed of the generator is not characterized by vibration.
(2) Analysis of the suppression of generator speed vibration for a varying f ( ) parameter with respect to α ≠ β.  Table 7 are taken to obtain the anti-interference performance waveforms of the generator speed to drivetrain vibration signals, as shown in Figure 20a,b. In the table, n1 = 0.08, n2 = 0.5, and the α/β values are mutually different. The vibration waveform of the two-mass drivetrain remains constant, as shown in Figure 19a. The generator speed waveform with the drivetrain vibration signal is shown in Figure 20a,b.

Conclusion
The closer the value of α/β is to 1, the better the suppression performance of the periodic disturbance.
According to the anti-interference performance analysis, it can be concluded that the disturbance suppression ability of f ( ) of the periodic disturbance signal differs when α ≠ β. The closer the value of α/β is to 1, the stronger the ability of the system to suppress periodic disturbances.  Table 8 are considered to obtain the anti-interference performance waveforms of the generator speed to the drivetrain vibration signals, as shown in Figure 21a-c. The value of α/β is different for n1 = 0.08 and n2 = 0.5 in the table. However, the specific values of α and β are also different. The vibration waveform of the two-mass drivetrain remains constant. The waveform of the generator speed interference of the periodic vibration of the drivetrain is shown in Figure 21a-c.   According to the vibration suppression analysis, it can be concluded that different values of α and β have different effects on the suppression of generator speed vibration when α/β = 0.3. The smaller the values of α and β, the smaller the speed vibration is, and the stronger the periodic disturbance suppression ability. This is because the smaller the values of α and β are, the smaller the influence of S in the transfer function F f (s) is, i.e., it is closer to 1.

The Vibration Suppression Ability Comparison of the 3-DOF-IMC and the EDM
(1) Experimental scheme: The experimental scheme of the generator speed control is shown in Figure 22, and the drivetrain parameters are from Table 1. The speed detection and speed given signal have filtering parts.
The dashed box in Figure 22 is the speed loop, which is generally designed by the EDM, and the transfer function of the speed regulator ASR is chosen as where: K n is the proportional gain of the speed controller; τ n is the integration time constant of the speed controller. By substituting Equation (10) into Figure 22, the small time constants T on and 2T oi are combined and approximated into an inertial part with a time constant T Σn = T on + 2T oi , so as to obtain the integral time constant of the speed controller Gain coefficient of speed controller where: α is the feedback coefficient of speed, using unit feedback; h is the intermediate frequency bandwidth, which is generally set as h = 5 in engineering [28]. The generator moment of inertia J F and parameters K ψ /β have been defined previously.   Experimental results and analysis: The generator's speed waveform is shown in Figure 23a, with parameters in scheme (a). Line 1 is the generator's speed waveform based on the 3-DOF-IMC; Line 2 is the speed waveform based on the EDM; Line 3 is the reference speed waveform in the diagram. The generator's speed waveforms of different values of α/β are shown in Figure 23b,c. In the figure, the parameters of the PI controller designed according to the EDM are shown in Equations (30) and (31), and its vibration amplitude overshoot is 5%; however, different values of α/β are characterized by different performances in suppressing periodic disturbances. With the increase in α/β, the overshoot of the vibration amplitude decreases. When α/β = 0.66, the amplitude overshoot of the generator speed using the 3-DOF-IMC method is the same as the amplitude overshoot using the EDM. When α/β = 0.82, the amplitude overshoot of the generator using the 3-DOF-IMC method is reduced to 50%. When α/β = 1, the amplitude overshoot of the generator speed using the 3-DOF-IMC method is only 4% of the amplitude overshoot of the EDM, which is also a unit feedback, and the vibration amplitude of the 3-DOF-IMC method is greatly reduced. However, when the value of α/β is less than 0.66, the vibration suppression effect of the EDM is better than that of the 3-DOF-IMC method. In Table 9, the vibration suppression ability comparison of the 3-DOF-IMC and the EDM with corresponding conclusions is shown.

Generator's Speed Tracking and Anti-Interference Performance Based on 3-DOF-IMC, Compared with EDM
(1) Experimental scheme: The experimental scheme of the generator's speed tracking and anti-interference performance based on the EDM is obtained according to reference [28]. Based on the EDM, the controlled object's mathematical model can be expressed as Equation (32):   [28]. Based on the EDM, the controlled object's mathematical model can be expressed as Equation (32): Thus, the integral time of the speed controller is τ n = hT Σn = 1; the gain coefficient of the speed controller is k n = J F β αp m k ψ · h+1 2hT Σn = 2.
(2) Experimental results and analysis: The generator's speed waveform is shown in Figure 24. The solid Line-A in the figure is the generator's output speed of the tracking performance and the anti-interference performance based on the 3-DOF-IMC. The dashed Line-B is the generator's speed of the tracking performance and the anti-interference performance based on the EDM. The tracking and anti-interference performance controlled by the 3-DOF-IMC is better than that of the PI controller designed by the EDM.

Generator's Speed Tracking and Anti-Interference Performance Based on 3-DOF-IMC, Compared with 1-DOF-IMC
(1) Experimental scheme: The experimental scheme of the generator's speed tracking and anti-interference performance based on the 1-DOF-IMC method is obtained ac cording to reference [28], the mathematical model of the controlled object can be ex pressed as Equation (28).
The 1-DOF-IMC of the generator is as follows: (1) Experimental scheme: The experimental scheme of the generator's speed tracking and anti-interference performance based on the 1-DOF-IMC method is obtained according to reference [28], the mathematical model of the controlled object can be expressed as Equation (28). The 1-DOF-IMC of the generator is as follows: C IMC−F (s) = G −1 − (s)L(s) = s(0.002s + 1) 1.5 · 1 λs + 1 (2) Experimental results and analysis: The generator's speed waveform is shown in Figure 25. The solid Line 1 in the figure is the tracking and anti-interference performance when the parameter λ is 0.4. The dashed Line 2 is the tracking and antiinterference performance of the generator's output speed when the parameter λ is 1.6. It can be seen that when the parameter λ is changed, the tracking and the antiinterference performance of the generator's speed change at the same time, which cannot be adjusted independently. However, the tracking performance and antiinterference performance controlled by the 3-DOF-IMC can be adjusted flexibly.

Compared with 1-DOF-IMC
(1) Experimental scheme: The experimental scheme of the generator's speed tracking and anti-interference performance based on the 1-DOF-IMC method is obtained ac cording to reference [28], the mathematical model of the controlled object can be ex pressed as Equation (28).
The 1-DOF-IMC of the generator is as follows: (2) Experimental results and analysis: The generator's speed waveform is shown in Figure 25 The solid Line 1 in the figure is the tracking and anti-interference performance when the parameter λ is 0.4. The dashed Line 2 is the tracking and anti-interference perfor mance of the generator's output speed when the parameter λ is 1.6. It can be seen that when the parameter λ is changed, the tracking and the anti-interference perfor mance of the generator's speed change at the same time, which cannot be adjusted independently. However, the tracking performance and anti-interference perfor mance controlled by the 3-DOF-IMC can be adjusted flexibly.

Conclusions
In this paper, the method of the suppression of vibration was proposed to adjust the aerodynamic torque of the wind turbine via variable pitch control and to control the elec tromagnetic torque and speed of the generator in order to suppress its speed vibration Moreover, vector control and IMC were used for the doubly fed wind generator. A double closed-loop control scheme of wind generator speed based on 3-DOF-IMC was proposed The following conclusions can be drawn: Figure 25. When the parameter λ is changed, the tracking and anti-interference performance of the generator's speed based on 1-DOF-IMC.

Conclusions
In this paper, the method of the suppression of vibration was proposed to adjust the aerodynamic torque of the wind turbine via variable pitch control and to control the electromagnetic torque and speed of the generator in order to suppress its speed vibration. Moreover, vector control and IMC were used for the doubly fed wind generator. A double closed-loop control scheme of wind generator speed based on 3-DOF-IMC was proposed. The following conclusions can be drawn: The 3-DOF-IMC has three controllers with adjustable parameters and flexible adjustments that can independently control the tracking performance, anti-interference performance, and robustness performance of the wind generator speed loop.
The feedback filter has a suppression effect on the external periodic disturbance signal and can be used to suppress the periodic disturbance from the drivetrain. Thus, the influence of drivetrain vibration on the wind generator speed can be eliminated.
A 3-DOF-IMC plays an important role in decoupling the cross-coupling voltage of the generator. When the controlled object and forecast model are matched, the wind generator rotor resistance and inductance parameters are proportionally increased, reduced, or remain constant, and the decoupling effect of the 1-DOF-IMC is good. However, when the wind generator's rotor inductance and resistance parameters increase, decrease, remain constant, or change in reverse according to different proportions, the 1-DOF-IMC can hardly achieve full decoupling due to a lower number of adjustable parameters. Full decoupling can be achieved through multi-parameter coordination and by employing the 3-DOF-IMC. Further work is needed in the future.