Integrated Control for Small Power Wind Generator

: The control strategies of the small power wind generator are usually divided into the maximum power point tracking (MPPT) case, which requires the wind generator produce power as much as possible, and the power limited control (PLC) case that demands the wind generator produce a power level following the load requirement. Integration of these two operating cases responding to ﬂexible and sophisticated power demands is the main topic of this article. A small power wind generator including the sluggish mechanical dynamic phenomenon, which uses the permanent magnet synchronous generator, is introduced to validate different control methods integrating MPPT and PLC cases and based on hysteresis control. It is a matter of an indirect power control method derived from three direct methods following perturb and observe principle as well as from a look-up table. To analyze and compare the proposed power control methods, which are implemented into an emulator of a small power wind generator, a power demand proﬁle is used. This proﬁle is randomly generated based on measured rapid wind velocity data. Analyzing experimental results, from the power viewpoint, all proposed methods reveal steady-state error with big amount of peak resulting from the nature of perturb and observe.


Introduction
After decades of development, the concern level about wind power research and application is still escalating. Being a clean and renewable energy resource, the wind power generator system extracts kinetic energy transforming into electrical form. While the wind energy source is apposite to provide the utility grid, small power wind systems can be used mostly as local distributed energy sources or one part of a microgrid since it is easy to be implemented and maintained. The scale of small power is defined by the power production level. In the United States of America, the wind energy system products lower than 100 kW are named small power wind systems; but, the criterion in Europe is 50 kW [1,2]. Due to its small size and light weight as well as its advances in aspects of reliability, energy density, and efficiency, the most used small power wind generator is the permanent magnet synchronous machine (PMSM) [3,4].
The small power wind energy system usually operates at maximum power point tracking (MPPT) case when wind velocity is lower than rated wind velocity; and constant power output is demanded when wind velocity is over rated wind velocity. However, being the local distributed energy source or a part of microgrid, the small scale wind energy system can also be required for the power limited control (PLC) case, at the range lower than rated wind velocity. Consequently, a general control strategy which could deal with both MPPT and PLC cases is necessary to be implemented.
As discussed and presented in [5], the wind kinetic energy is converted into DC electrical energy by means of three possible solutions: the passive electrical structure using a three-phase diode bridge; torque still needs to be avoided. At the same time, the determination of parameters of sliding model control requires experienced designers. In [16], the authors selected Tip Speed Ratio (TSR) as the feedback parameter of the P&O method, and then compare the actual value of TSR with the ideal value calculated by mathematical model of studied system. Nevertheless, even it used the principle of P&O, and the method still performs with weak robustness. Different from the achievement discussed above, the research objective of the present work focused on the more general and flexible limited power demands, which can be applied in power balancing of a microgrid [5,10,11], when the microgrid may include the small scale wind generator.
This paper presents the modelling of a small power wind PMSM in order to validate the power control strategy methods with hysteresis control loop. Proposed power control methods are designed to cover MPPT and PLC operating cases. In addition, this study faces precisely the problem of the uncertainty of wind velocity and the demanded power amount from load. Therefore, regarding the robustness of the studied system, it should be considered beside the main focus which is the wind speed. One notes that, in this paper, the imposed wind reference is modeled on a real wind measurement with many variations in frequencies and amplitudes. The small power wind turbine is emulated by a test bench and several experiments are implemented to validate characteristics of all proposed methods. Finally, experimental results are given, compared, and analyzed, and strengths and weaknesses of each method are revealed together.
The article is structured as follows. Section 2 presents the small power conversion system including the analysis of characteristics. The formulation of the research problem and proposed power control methods integrating MPPT and PLC operating cases are presented in Section 3. In Section 4, the capability and effectiveness of all methods are evaluated based on a test bench. Conclusions and further studies are given in Section 5 and the final section presents a nomenclature.

Small Scale Wind Energy Conversion System
The study in this article concentrates on a system transforming electrical energy into DC form, allowing the output system to supply networks as follows: a DC grid, an AC grid by using an inverter, and microgrids with DC or AC bus [5,10,11]. The power conversion structure implemented for this study, shown in Figure 1 [5], comprises a three-phase diode bridge connected to a controllable DC-DC converter. Based on this combination, the optimization of energy is able to be realized by searching the maximum power operating point or limited power operating point.
Due to their simplicity and robustness in application, P&O methods are chosen to be realized and studied in this article to deal with the issue of PLC which means the wind generation system generates the power value required by loads. The test bench of the small power wind energy conversion system, based on the emulator of a small power wind turbine, is presented in this section and then the distribution of steady-state working points is displayed and used to determine the control strategy.

Overview of Test Bench
The test bench electric scheme is presented in Figure 1a. As presented in [5], a system that consists of a three-phase driver, a three-phase PMSM, and one dSPACE DS1104 control board emulates the dynamic behavior of wind and blades. A second PMSM is used as generator of whole wind energy conversion system. Generator's output connects a three-phase diode bridge transforming electric energy from AC form into DC form. A capacitor C BUS is implemented to stabilize the voltage of DC bus. An inductance L and an Insulated Gate Bipolar Transistor (IGBT) module driving by dSPACE DS1104 are combined into a DC-DC boost convertor to realize the adjustment of operating point. A programmable electronic load (PEL), shown in Figure 1a, for which the C PEL is associated, emulates the power demand of user. The PEL keeps voltage at 400 V for all operating condition [5,10]. Detailed information of equipment of the test bench is given in Table 1 and all physic equipment is shown in Figure 1b.

Electrical Power Distribution
The emulating model of the wind energy conversion system is considered as a 1 kW system by Bergey [5]. The extracted amount of electrical power is based on the aerodynamic power AERO p expressed by Equation (1): where the density of air  is 1.23 kg/m 3 , the blade radius R is 1.25 m, the v indicates the wind velocity, and the p c is the power coefficient that is approached by a 7th-order polynomial function as in Equation (2). This power coefficient is given in Figure 1c; it shows the behavior of the emulated wind turbine. The dynamic mechanical function [5] is given by Equation (3).

Electrical Power Distribution
The emulating model of the wind energy conversion system is considered as a 1 kW system by Bergey [5]. The extracted amount of electrical power is based on the aerodynamic power p AERO expressed by Equation (1): where the density of air ρ is 1.23 kg/m 3 , the blade radius R is 1.25 m, the v indicates the wind velocity, and the c p is the power coefficient that is approached by a 7th-order polynomial function as in Equation (2). This power coefficient is given in Figure 1c; it shows the behavior of the emulated wind turbine. The dynamic mechanical function [5] is given by Equation (3). 1 In Equations (2) and (3), the λ is the TSR calculated following λ = RΩ/v, in which Ω is the rotational speed used also in Equation (3), p EM is the electromagnetic power, J and F are separately the mechanical inertia equaling 1.5 kg·m 2 and the viscous damping equaling 0.06 Nm/rad.
The distribution of static state operating points is validated based on several pre-research experiments. Figure 2a,b presents the aerodynamic power p AERO and the electrical power p BUS following the DC bus voltage u BUS and the bus current i BUS as variables. The correlation of these two powers, p AERO and p BUS is displayed in Figure 2c.
In Equations (2) and (3), the λ is the TSR calculated following The analysis of the experimental results given in Figure 2 indicates that: (i) operating maximum power points (  The analysis of the experimental results given in Figure 2 indicates that: (i) operating maximum power points (MPPs) are not identical for p BUS curves and p AERO curves owing to losses of power from the input side of PMSM to the output side of three-phase diode bridge; and (ii) for limited power demand, p AERO curves supply operating points different from p BUS curves. Thus, for integrated power control strategy, p AERO curves are useless even as indirect references. Therefore, the control strategy and the regulating algorithm application have to take into account this remark. Based on the characteristics of p BUS curves, the variable chosen to realize the change of system's operating point is the DC bus voltage. For each wind speed it is axiomatic that for one value of i BUS there are two different working points within the operating range. Consequently, comparing with i BUS , the variable u BUS is more suitable to be regarded as the reference variable, since each value of u BUS presents a unique operating state.

Problem Formulation
Integration of MPPT and PLC cases is the main objective of this paper. The proposed principle of this power control method is to avoid the requirement of distinguishing the MPPT condition from the PLC condition. Based on different theories of direct and indirect MPPT methods, ways to realize the integration of MPPT and PLC methods are dissimilar and they are presented below.

Control Loop
Whatever the principle, power control methods can only describe movement of system's steady operating points. Based on the characteristic of studied wind power conversion systems, the bus voltage u BUS is selected as the control variable. Hence, before explanations of proposed power control methods, the dynamic characteristic of whole wind power conversion system with respect to the variation of u BUS needs to be analyzed first. To correctly operate the regulating algorithm, the control loop is regulated by a hysteresis controller due to its good dynamic performance and easy implementation [17]. Figure 3a,b show the response of the actual u BUS with respect to the variation of bus voltage reference u * BUS for the hysteresis controller, and respectively its enlarged version. The power response as well as the three phase electrical variables response are analyzed and presented in Figure 3c,d, with enlarge presentation of three-phase electric variables given in Figure 3e,f.
Integration of MPPT and PLC cases is the main objective of this paper. The proposed principle of this power control method is to avoid the requirement of distinguishing the MPPT condition from the PLC condition. Based on different theories of direct and indirect MPPT methods, ways to realize the integration of MPPT and PLC methods are dissimilar and they are presented below.

Control Loop
Whatever the principle, power control methods can only describe movement of system's steady operating points. Based on the characteristic of studied wind power conversion systems, the bus voltage BUS u is selected as the control variable. Hence, before explanations of proposed power control methods, the dynamic characteristic of whole wind power conversion system with respect to the variation of BUS u needs to be analyzed first. To correctly operate the regulating algorithm, the control loop is regulated by a hysteresis controller due to its good dynamic performance and easy implementation [17]. Figure 3a,b show the response of the actual BUS u with respect to the variation of bus voltage reference BUS u  for the hysteresis controller, and respectively its enlarged version. The power response as well as the three phase electrical variables response are analyzed and presented in Figure 3c,d, with enlarge presentation of three-phase electric variables given in Figure 3e,f. Then, different amplitudes of step input of u BUS were implemented into the test bench, based on 8 m/s wind velocity and the 150 V initial value of the bus voltage. In order to highlight the regulation time of p BUS response to the variation of u BUS , Figure 4 presents the relationship between the regulation time, variation value of p BUS , and the step input values of u BUS . The used values of step input of u BUS are: 10 V, 7.5 V, 5 V, and 2.5 V.
Given the results given in Figure 3 and the discussion expressed in [18], one notes that the hysteresis controller is adequate for a rapid response. Indeed, the wave forms of currents and voltages may be considered as acceptable concerning the distortion. Therefore, mutual restraint resulting from  Given the results given in Figure 3 and the discussion expressed in [18], one notes that the hysteresis controller is adequate for a rapid response. Indeed, the wave forms of currents and voltages may be considered as acceptable concerning the distortion. Therefore, mutual restraint resulting from the joint action of the power control methods and the hysteresis control algorithm has to be taken into account.
Due to the non-linearity of the system, as presented in Figure 4, the value of regulation time of output power BUS p responding the variation of BUS u is associated to the perturb step-size u  of as accurate as possible, a low-pass filter, which uses a cut-off frequency of 6000 Hz, is added before the entrance of direct methods. This is a compromise between the precision and the rapidity of tracking.
If the BUS p is used by the power control method, certain time step (which should not be smaller than 2 s) has to be added to avoid the dynamic process of BUS p response. However, for the method based on the indirect principle, this operating time step is unnecessary, because used parameters do not include the BUS p or any other parameters containing the high order dynamic characteristics.

Principle of Power Control Strategy Methods
Based on our previous studies [5,10,11], the PLC operating case means that the power output of the small scale wind generator should follow the demanded value supplied by the users' side. Thus, the power control method should focus on the absolute value of the DIFF where LIM p is the asked power value and the BUS p is the measured actual power value, with the acceptable assumption that LIM p has been obtained. Due to the non-linearity of the system, as presented in Figure 4, the value of regulation time of output power p BUS responding the variation of u BUS is associated to the perturb step-size ∆u of u BUS and the operating point. According to Figure 4, when demanded step input of u BUS is big, e.g., +10 V, +7.5 V, and even +5 V, the output power needs approximately 2.6 s to stable in the region about 5% of stable value of p BUS . For a decreased step input, e.g., +2.5 V, the value of the regulation time still rests near to 2 s.
To implement the control strategy for all methods demanding information of p BUS and u BUS as accurate as possible, a low-pass filter, which uses a cut-off frequency of 6000 Hz, is added before the entrance of direct methods. This is a compromise between the precision and the rapidity of tracking.
If the p BUS is used by the power control method, certain time step (which should not be smaller than 2 s) has to be added to avoid the dynamic process of p BUS response. However, for the method based on the indirect principle, this operating time step is unnecessary, because used parameters do not include the p BUS or any other parameters containing the high order dynamic characteristics.

Principle of Power Control Strategy Methods
Based on our previous studies [5,10,11], the PLC operating case means that the power output of the small scale wind generator should follow the demanded value supplied by the users' side. Thus, the power control method should focus on the absolute value of the p DIFF = |p BUS − p LI M |, where p LI M is the asked power value and the p BUS is the measured actual power value, with the acceptable assumption that p LI M has been obtained. Figure 5 presents the proposed power controlling strategy integrating both of MPPT and PLC cases. Dotted lines in this figure indicate the varying demanded value of power: the upper one shows the demanded values for MPPT case; the lower dotted line indicates the PLC case. Additionally, the solid line displays the power distribution for a given wind velocity. Thus, following the P&O logic principle, the proposed integrating power control methods will be committed to minimize the value of the difference p DIFF .
According to the characteristics of the wind turbine, when the demanded value of p LI M is lower than the potential maximum power value, there are theoretically two operating points supplying the required power value. Considering the transfer of operating points, if the steady component of operating currents is close to equipment's physical limitations, the dynamic components of current resulting from the transfer between operating points has high possibility for damaging devices. Therefore, in this work, the system will be forced to operate at the "low current-high voltage" side.
Aiming to realize the objective mentioned above, the operating region has been divided into four parts: I, II, III and IV which are defined by changes of variables ∆u BUS , ∆p BUS and ∆p di f f presented in Figure 6 and Table 2. for a given wind velocity. Thus, following the P&O logic principle, the proposed integrating power control methods will be committed to minimize the value of the difference DIFF p . According to the characteristics of the wind turbine, when the demanded value of LIM p is lower than the potential maximum power value, there are theoretically two operating points supplying the required power value. Considering the transfer of operating points, if the steady component of operating currents is close to equipment's physical limitations, the dynamic components of current resulting from the transfer between operating points has high possibility for damaging devices. Therefore, in this work, the system will be forced to operate at the "low currenthigh voltage" side. Aiming to realize the objective mentioned above, the operating region has been divided into four parts: I, II, III and IV which are defined by changes of variables  Figure 6 and Table 2.
Based on the Table 2  II. So, since it is desired to drive the system operating at the "low current-high voltage" side, the next perturbation step is selected as positive again.  According to the characteristics of the wind turbine, when the demanded value of LIM p is lower than the potential maximum power value, there are theoretically two operating points supplying the required power value. Considering the transfer of operating points, if the steady component of operating currents is close to equipment's physical limitations, the dynamic components of current resulting from the transfer between operating points has high possibility for damaging devices. Therefore, in this work, the system will be forced to operate at the "low currenthigh voltage" side. Aiming to realize the objective mentioned above, the operating region has been divided into four parts: I, II, III and IV which are defined by changes of variables BUS u  , BUS p  and diff p  presented in Figure 6 and Table 2.
Based on the Table 2  is also positive, the operating point can be identified at this moment in the region II. So, since it is desired to drive the system operating at the "low current-high voltage" side, the next perturbation step is selected as positive again.
Based on the Table 2 and the Figure 6a, for example, after one step of perturbation, if the change of DC bus voltage, ∆u BUS , is positive, the change of DC bus power, ∆p BUS , is still positive and the change of p DIFF is also positive, the operating point can be identified at this moment in the region II. So, since it is desired to drive the system operating at the "low current-high voltage" side, the next perturbation step is selected as positive again.
The Figure 6a and Table 2 indicate the power limited demand conditions, for which the demanded value is lower than the potential maximum power value. The Figure 6b and Table 3 present the region definition of the MPPT demand conditions. It is obvious that, to these two operating conditions, operating regions can be defined by same rules. Therefore, P&O algorithms based on these operating regions' definition can uniformly cover two conditions of MPPT and PLC. Based on trends of changes of ∆u BUS , ∆p BUS , ∆p DIFF , power control methods based on P&O theory determine the direction of change of system operating state.

Fixed Step-Size
In order to validate the control strategy mentioned above, one simple experiment was designed. Figure 7 present the selected wind velocity and the demanded power value profiles, as well as experimental results of actual electrical power, u BUS and the perturb step-size. The profile of p LI M demands 250 W, which is lower than the maximum power level for given wind velocity, and 378 W, which is the maximum power level. The first step is to confirm this strategy works with fixed perturbation step-size.

Fixed Step-Size
In order to validate the control strategy mentioned above, one simple experiment was designed. Figure 7 present the selected wind velocity and the demanded power value profiles, as well as experimental results of actual electrical power, BUS u and the perturb step-size. The profile of LIM p demands 250 W, which is lower than the maximum power level for given wind velocity, and 378 W, which is the maximum power level. The first step is to confirm this strategy works with fixed perturbation step-size. According to results in Figure 7, it is clear that the principle presented in Section 3.2 works properly for identifying the actual operating region and detecting the moving direction of perturbation. However, obviously, fixed perturbation step-size cannot drive the system producing stable power output; so, using variable perturbation step-size is the next object to be implemented. According to results in Figure 7, it is clear that the principle presented in Section 3.2 works properly for identifying the actual operating region and detecting the moving direction of perturbation. However, obviously, fixed perturbation step-size cannot drive the system producing stable power output; so, using variable perturbation step-size is the next object to be implemented.

Improved Variable Step-Size Designed with Newton-Raphson Technique
Following the logic of MPPT variable step-size method [3,18], Newton-Raphson method is one worthwhile method of calculating variable step-size for power control method integrating MPPT and PLC cases. Based on the theory of Newton-Raphson method, p DIFF is regarded as one function of bus voltage, p DIFF = f (u BUS ). Thus, the calculation of variable step-size matches the following iterative methodology for each k step: The value of Sign is determined based on Table 2. When the system works in I, II or III region, its value equals one; otherwise, the system works in IV region and the value Sign equals minus one. To avoid the convergence of variable step-size at the low voltage-high current operating conditions Energies 2018, 11, 1217 10 of 16 whose region is located at the left side of maximum power operating points, the variable step-size calculating method is used when system works in III and IV operating regions. Figure 8 presents the experimental results of Newton-Raphson method under the same experimental condition as in Figure 7.
The value of Sign is determined based on Table 2. When the system works in I, II or III region, its value equals one; otherwise, the system works in IV region and the value Sign equals minus one. To avoid the convergence of variable step-size at the low voltage-high current operating conditions whose region is located at the left side of maximum power operating points, the variable step-size calculating method is used when system works in III and IV operating regions. Figure 8 presents the experimental results of Newton-Raphson method under the same experimental condition as in Figure 7.  Figure 8). If, the maximum value of perturbation step-size is reduced, the effect of convergence of perturbation step-size is not improved significantly but the seeking speed to a new steady-state will clearly also be decreased. Based on Figure 8, when p LI M equals 250 W, which means the PLC case, the variation of p BUS is significant; even when p LI M equals 380 W, MPPT case, the variation of the stationary component of p BUS around the MPP is not obvious so much. It is clear that when the maximum value of perturbation step-size is selected big, such as 10 V, the perturbation step-size cannot significantly converge to zero (shown in Figure 8). If, the maximum value of perturbation step-size is reduced, the effect of convergence of perturbation step-size is not improved significantly but the seeking speed to a new steady-state will clearly also be decreased.
According to experimental results such as in Figure 8, direct application of Newton-Raphson method cannot converge properly. This is due to the fact that the rate of change of the power respect to the voltage, around operating points of limited power, is usually high, comparing to the response of power nearby MPP with same size of perturbation step. Therefore, one improved Newton-Raphson variable step-size iterative methodology is introduced as follows: Based on the results presented in Figure 8, the calculation iteration (Equations (4)-(7)) of variable step-size around the MPPs naturally converge smooth, since, at that place, the distribution of operating points has a large range where the change of power is not significant for at least one perturb. So, dividing classic Newton-Raphson calculation iteration by the gradient of the actual power respecting DC bus voltage could reduce the power transition if the operating goal is PLC. Taking into account the modification in Equation (10), it can be stated that with this variable step-size algorithm, this method can deal with the convergence of variable step-size for both conditions of MPPT and PLC.

Variable Step-Size with Fuzzy Logic
Another well-known technique of variable perturbation step-size is the fuzzy logic method [3,18], which can handle robust and nonlinear control requirements. However, one notes that this method requires in-depth knowledge concerning the system. Fuzzy logic is based on the input space and the output space mapping through logical actions. As usual, in this work, the process is divided into three phases: the fuzzification, the fuzzy reasoning, and the defuzzification.
Same as the method presented in Section 3.2.2, the fuzzy logic method is introduced into the calculation of variable step-size. So, gradi(k) and p DIFF (k) are chosen as two inputs of fuzzy logic method as presented in Equation (11): During the fuzzification phase, the inputs and outputs numerical values are expressed in fuzzy sets. All points in input space are translated into the "code" in form of the degree of membership which can equal any values between 0 and 1. Considering the universe of discourse of the inputs and the output normalized into [-1, 1], Figure 9b-d present the fuzzy subsets corresponding to the inputs e 1 ∈ {bn; n; z; p; bp} and e 2 ∈ {n; z; p}, and to the output s ∈ {−; 0; +}. Aiming to accomplish fuzzy logic as simple as possible, following the principle presented in Figure 9a, these input variables, the gradi(k) is normalized into e 1 by gain K 1 , and p DIFF (k) is normalized into e 2 by gain K 2 . However, there is possibility to adjust values of K 1 and K 2 , to obtain a normalized interval rather than [-1, 1], which can lead to a better regulating performance (presented in Section 4). Forms of membership function were selected as triangular and trapezoidal functions to simplify the study. For e 1 , a group of five subsets is demonstrated to separate its value range for MPPT and PLC cases. Parameters of each membership function are determined based on experimental tests.
The Figure 9e shows the fuzzy reasoning surface of the proposed fuzzy logic and the Table 4 explains the "if-then" reasoning role. The step of "if" describes a fuzzy area in the input space and the step "then" step states the output in the fuzzy area. With a logical operator, membership functions values, evaluating the degree of each activated rule, summarize the degree of satisfaction for the "if" step of each rule. Subsequently, the fuzzy set in the "then" step of each rule is determined. Finally, outputs of each rule are aggregate into one fuzzy set. The Mamdani Fuzzy Inference System using the fuzzy toolbox is applied in this work.
The defuzzification wishes to "convert" the fuzzy set resulted from the fuzzy reasoning into a numerical value between 0 and 1. Then this value is anti-normalized to be ∆u by K s . Several defuzzification operators are mentioned in [19]; however, the most used operator is the centroid, or center of gravity, which is chosen in this study, such as the implemented methods mentioned above. Figure 10 presents the experimental results of fuzzy logic method under the same operating condition as presented in Figures 7 and 8.  The Figure 9e shows the fuzzy reasoning surface of the proposed fuzzy logic and the Table 4 explains the "if-then" reasoning role. The step of "if" describes a fuzzy area in the input space and the step "then" step states the output in the fuzzy area. With a logical operator, membership functions values, evaluating the degree of each activated rule, summarize the degree of satisfaction for the "if" step of each rule. Subsequently, the fuzzy set in the "then" step of each rule is determined. Finally, outputs of each rule are aggregate into one fuzzy set. The Mamdani Fuzzy Inference System using the fuzzy toolbox is applied in this work.
The defuzzification wishes to "convert" the fuzzy set resulted from the fuzzy reasoning into a numerical value between 0 and 1. Then this value is anti-normalized to be u  by Ks. Several defuzzification operators are mentioned in [19]; however, the most used operator is the centroid, or center of gravity, which is chosen in this study, such as the implemented methods mentioned above. Figure 10 presents the experimental results of fuzzy logic method under the same operating condition as presented in Figures 7 and 8.

Analysis of Comparative Results
In order to compare above-mentioned power control methods, a real wind velocity is considered as the wind velocity input, presented in Figure 11a, from data measured by Météo France, in Compiegne, France, on 15 January 2015, during 15 min; this real wind velocity profile is the same as in [18]. Furthermore, the LIM p profile was chosen based on the potential maximum power and the physical minimum power of the test bench, as presented in Figure 11b. The selected wind profile and power limited profile highlight a strained condition with rapid changes.

Analysis of Comparative Results
In order to compare above-mentioned power control methods, a real wind velocity is considered as the wind velocity input, presented in Figure 11a, from data measured by Météo France, in Compiegne, France, on 15 January 2015, during 15 min; this real wind velocity profile is the same as in [18]. Furthermore, the p LI M profile was chosen based on the potential maximum power and the physical minimum power of the test bench, as presented in Figure 11b. The selected wind profile and power limited profile highlight a strained condition with rapid changes. (c)

Analysis of Comparative Results
In order to compare above-mentioned power control methods, a real wind velocity is considered as the wind velocity input, presented in Figure 11a, from data measured by Météo France, in Compiegne, France, on 15 January 2015, during 15 min; this real wind velocity profile is the same as in [18]. Furthermore, the LIM p profile was chosen based on the potential maximum power and the physical minimum power of the test bench, as presented in Figure 11b. The selected wind profile and power limited profile highlight a strained condition with rapid changes. Under these experimental conditions, several experimental tests are introduced to compare all above-mentioned power control methods whose characteristics are given as follow:
As key variable, p DIFF is used to analyze characteristics of each method. Referencing statistical concepts, the 'mean', defined as average value, and the 'variance', indicating the expected value of the squared deviation from the mean, are calculated to compare dynamic and steady-state characteristics of different power control methods. At first, to all methods based on the P&O theory, a sampling method is a sampling time equals two seconds, is applied to filter the peak of each perturbation step of P&O (marked as 'Mean sampling' and 'Variance sampling'). Then, the overall comparison (marked as 'Mean overall' and 'Variance overall') is calculated and analyzed also. Based on those parameters, results for the optimal one test of each kind of method are collected and presented in Figure 12. The 'Mean' and 'Variance' of those selected tests are listed in Table 5.
Summarizing all presented results, under complex conditions of wind velocity and power demands variation, all three methods based on the P&O principle matched the design target of this integrating power control method. Those methods perform weakly in the view of dynamic, resulting from the mechanical inertia which cannot be ignored; but, overall, their steady-state errors are tiny. From fixed step-size to improved Newton-Raphson method calculating variable step-size till to fuzzy logic method, the dynamic noise is reduced significantly, meanwhile the change of whole steady-state error is not much, in the 'overall' view.
Considering the complexity of implementing each power control method, fixed step-size and improved Newton-Raphson method has just one parameter needing to be determined: the value or the maximum limit of step-size. In theory, variable step-size calculated by fuzzy logic has a greater possibility to perform better, since it contains several parameters to be modified. However, this greater freedom increases the complexity of determination of those parameters, which demands a rich experience of application. a sampling method is a sampling time equals two seconds, is applied to filter the peak of each perturbation step of P&O (marked as 'Mean sampling' and 'Variance sampling'). Then, the overall comparison (marked as 'Mean overall' and 'Variance overall') is calculated and analyzed also. Based on those parameters, results for the optimal one test of each kind of method are collected and presented in Figure 12. The 'Mean' and 'Variance' of those selected tests are listed in Table 5.  Summarizing all presented results, under complex conditions of wind velocity and power demands variation, all three methods based on the P&O principle matched the design target of this integrating power control method. Those methods perform weakly in the view of dynamic, resulting from the mechanical inertia which cannot be ignored; but, overall, their steady-state errors are tiny. From fixed step-size to improved Newton-Raphson method calculating variable step-size till to fuzzy logic method, the dynamic noise is reduced significantly, meanwhile the change of whole steady-state error is not much, in the 'overall' view.
Considering the complexity of implementing each power control method, fixed step-size and improved Newton-Raphson method has just one parameter needing to be determined: the value or the maximum limit of step-size. In theory, variable step-size calculated by fuzzy logic has a greater possibility to perform better, since it contains several parameters to be modified. However, this greater freedom increases the complexity of determination of those parameters, which demands a rich experience of application.

Conclusions
Three power control methods based on the P&O principle (fixed step-size; improved Newton-Raphson variable step-size and variable step-size based on fuzzy logic) have been studied, designed, and then implemented into a test bench highlighting the sluggish mechanical inertia and dealing with the power control problem integrating MPPT and PLC cases. This study concerns the small power wind energy conversion system.
One experiment, based on designed wind and power demand profile, was used to validate the basic function of each method. In addition, another experiment, using a measured wind velocity profile and one calculated power demand profile, was implemented to compare characteristics of all proposed methods. Evaluations of all methods and several statistical indices of key variable p DIFF were calculated and analyzed. All power control methods present good steady-state characteristics, but their dynamic characteristics are limited by the sluggish mechanical inertia and the control loop's capability. Furthermore, some of the complex control algorithms, such as fuzzy logic used in this work, have the potential to achieve better performance. Moreover, the studied system uses a PMSM assembly, three-phase diode bridge, and converter with a hysteresis control, which makes it one of the most robust "drivers".
The future research direction is to improve the dynamic characteristics of the control loop, to suppress the effect from the sluggish mechanical inertia. In addition, a deep robustness analysis may be conducted taking into account unmodeled dynamics, or associated couplings and uncertainties. This can significantly improve the performance of methods based on the P&O principle.
Author Contributions: All authors have conceived and designed the system, performed the experiments, and analyzed the data. All authors contributed jointly to the writing and preparing revision of this manuscript. All authors have read and approved the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest and the founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the result.   Figure 11a