Adaptive MPC Control for Wind Power Systems with VRB Storage Using SVR-Based Sensorless Estimation and FPNN-IPSO Optimization

1. Introduction

Fuzzy neural networks (FNNs), which combine the reasoning capabilities of fuzzy logic with the learning ability of neural networks, have gained increasing attention in control and modeling applications due to their ability to manage time-varying uncertainties [1,2,3]. Compared to traditional fuzzy systems or pure neural networks, FNNs offer better adaptability and robustness in nonlinear and uncertain environments. However, conventional FNNs are often inadequate for handling stochastic uncertainties, which are common in renewable energy systems [4,5,6].

To address these limitations, researchers have proposed integrating probabilistic neural networks (PNNs) with fuzzy structures. PNNs, based on Parzen window estimation and Bayes classification [7,8,9], are well-suited for managing probabilistic inputs and have been successfully applied in classification, fault detection, and nonlinear mapping tasks. A hybrid approach known as the fuzzy probabilistic neural network (FPNN) has emerged, aiming to combine the learning efficiency of FNNs and the uncertainty-handling capacity of probabilistic neural networks (PNNs) [10]. This architecture has demonstrated potential in various control scenarios involving stochastic disturbances [11,12].

In the context of wind power systems, maximum power capture (MPC) is essential for ensuring energy capture efficiency under fluctuating wind conditions. Recent studies have explored sensorless MPC techniques that eliminate the need for mechanical wind speed sensors, thereby enhancing system reliability and reducing costs. [13]. Among these techniques, three main categories have been identified. The first uses polynomial approximations of the power coefficient to infer wind speed [14,15], but the associated real-time root-solving process for high-order polynomials imposes a significant computational burden. The second approach utilizes power-mapping tables and lookup mechanisms [16], which suffer from memory constraints and limited resolution. The third employs turbine torque and rotor speed estimation [17], but requires an accurate torque observer, which complicates implementation.

Given the nonlinear and uncertain characteristics of wind energy systems, heuristic optimization methods have also been employed in controller design. Particle Swarm Optimization (PSO), introduced by Kennedy and Eberhart, is known for its fast convergence and ease of implementation. In particular, enhancements such as time-varying acceleration coefficients (TVAC) have been proposed to improve PSO's global search capabilities. Compared to genetic algorithms, PSO retains knowledge of previous optimal solutions, enabling more stable convergence in complex search spaces [18,19,20].

Support Vector Regression (SVR) has also been recognized as a powerful tool for nonlinear regression tasks, offering high generalization ability and robustness to local minima due to its convex optimization formulation [21]. SVR constructs a predictive model based on support vectors, making it particularly effective for wind speed estimation under sparse or noisy data conditions. When applied to wind energy systems, SVR can estimate real-time wind speed with high precision, serving as a valuable input for sensorless MPC control.

Motivated by the limitations of conventional control methods and recent advances in hybrid learning frameworks, this paper proposes a novel MPC control strategy for a permanent magnet synchronous generator (PMSG) wind turbine, based on an FPNN controller optimized by improved PSO (IPSO), and combined with SVR-based wind speed estimation. The objective is to improve power tracking accuracy and system stability without relying on mechanical wind sensors. Simulation results demonstrate the proposed method’s effectiveness compared to conventional PI and fuzzy control techniques.

2. Analytical Evaluation of the Wind Power Generation System

This study focuses on the control and modeling of a wind power generation system equipped with a PMSG and integrated with a vanadium redox battery (VRB) in a microgrid environment. As illustrated in Figure 1, the control system incorporates three main components: a wind speed estimator based on SVR, a speed controller for MPC, and current controllers in the synchronous reference frame. The SVR model is trained offline using historical data and applied online for real-time estimation of wind speed. This estimation serves as the basis for determining the optimal generator speed to maximize power output. The overall performance of the system is verified through a MATLAB/Simulink microgrid model, incorporating wind energy generation, VRB storage, and three-phase inverters. The dynamic behavior of the integrated system is analyzed under varying wind conditions.

2.1. Wind Turbine Characteristics and Modeling

To maximize energy extraction from wind, variable-speed wind turbines are commonly employed alongside power electronic interfaces, which decouple the generator speed from the fixed grid frequency. The aerodynamic power extracted by a wind turbine is mathematically represented as:

$$P_m={\displaystyle \frac{1}{2}}\rho A{C}_p(\lambda ,\beta ){\upsilon }^3$$

(1)

where ρ denotes the air density, A represents the swept area of the turbine blades, v is the wind speed, and C_p is the power coefficient. The coefficient C_p is a function of the tip-speed ratio λ and the blade pitch angle β. The tip-speed ratio is defined as:

where r denotes the blade radius and ω_r represents the rotor angular speed. The power coefficient C_p generally exhibits a nonlinear relationship with the tip-speed ratio λ and blade pitch angle β, and is often characterized using empirical models. Once the wind speed is estimated via the SVR module, the reference rotor speed ω_r^∗ required for maximum power point tracking (MPPT) can be determined as follows:

Upon estimation of the wind speed using the SVR module, the reference rotor speed ω_r^∗ required for MPPT can be calculated as:

$${\omega }_r^{*}={\displaystyle \frac{{\lambda }_{opt}}{r}}v$$

(4)

This control strategy ensures that the turbine operates at the optimal tip-speed ratio, thereby maximizing aerodynamic efficiency.

2.2. PMSG

The wind turbine is coupled to a three-phase PMSG, which converts the captured mechanical energy into electrical power. The relationship between the mechanical torque Tm and the electromagnetic torque Te of the generator can be described by the following equations: [22,23]:

$$T_m={\displaystyle \frac{P_m}{{\omega }_r}}$$

$$T_e={\displaystyle \frac{P_e}{{\omega }_e}}={\displaystyle \frac{2}{n_p}}{\displaystyle \frac{P_e}{{\omega }_r}}$$

(5)

where P_m​ and P_e​ are the mechanical and electrical powers, respectively, ω_e​ is the electrical angular frequency, and n_p is the number of poles. The mechanical dynamics of the PMSG can be described by:

$$J{\displaystyle \frac{d{\omega }_r}{dt}}=T_m-(n_p/2)T_e$$

(6)

where J is the moment of inertia. This dynamic model is crucial for the design and analysis of the speed control loop.

2.3. Principle of a Vanadium Redox Battery

Vanadium redox batteries (VRBs) are a type of flow battery that store and release energy by utilizing vanadium ions in different oxidation states as redox couples. Their unique electrochemical structure allows complete separation between energy storage capacity and power delivery capability, which makes them particularly advantageous in applications with highly variable power inputs such as wind energy systems.

The vanadium redox battery (VRB) consists of two vanadium-based electrolyte reservoirs, circulation pumps, and an electrochemical cell stack separated by an ion exchange membrane. During charging, VO₂⁺ in the positive electrolyte is reduced to VO²⁺, and V³⁺ in the negative electrolyte is reduced to V²⁺; the reverse oxidation reactions occur during discharge. The membrane permits H⁺ ion transfer to maintain charge neutrality while preventing vanadium species crossover.

From a control perspective, the dynamic behavior of a VRB is governed by a set of electrochemical and hydraulic equations, including charge transfer kinetics, diffusion-convection dynamics, and voltage-current characteristics of the cell stack. These dynamics are important to consider when integrating VRBs with wind generation systems, particularly in microgrids where real-time power balancing is critical.

In this study, the VRB functions as a buffer to absorb the mismatch between the intermittent wind generation and the instantaneous load demand. This helps maintain system voltage and frequency stability, especially during sudden changes in wind speed or during plug-and-play events of distributed generators. Moreover, the VRB is controlled to charge when excess power is available and to discharge when the generated power falls below demand, effectively flattening the power output profile of the wind turbine.

VRBs are preferred over other energy storage technologies owing to their advantageous characteristics, such as long cycle life and a high depth of discharge, and rapid dynamic response, and intrinsic safety due to the aqueous nature of its electrolytes. Furthermore, because the power and energy components are physically separated, the system can be flexibly scaled for either higher power or longer duration storage without redesigning the entire system. Given these advantages, the VRB is selected in this study as the energy storage component to enhance short-term power balance in the microgrid. Its integration strategy is discussed in the subsequent simulation model [24,25].

3. A SVR Approach for Wind Speed Estimation

Regression methods are algorithms used to estimate the relevance among system inputs and outputs based on sample or training data. Among these methods, SVR has been demonstrated to be particularly effective for power system applications [26].

Consider a set of training samples $\left\{(y_1,x_1),\dots (y_n,x_n)\right\}$, where $x_i$ and $y_i$ represent the input and output spaces, respectively, and n denotes the number of training instances. The objective of a regression issue is to identify a functional relationship that accurately estimates unknown parameters based on observed input-output data. [26,27].

The general form of the SVR estimation function can be formulated as [26]:

Let $\omega$ denote the weighting vector, b the bias term, $\phi (x)$ denotes a nonlinear mapping function that transforms input data from n dimensions to a higher-dimensional feature space, as illustrated in Figure 2. The dot product symbol express the inner product in this feature space. To determine the optimal weighting vector $\omega$, the cost function must be minimized [28,29,30]:

subject to

where $\epsilon$ denotes the permissible error margin and $\gamma$ serves as a user-specified regularization parameter that regulates the trade-off between the complexity of the model and the tolerance for training errors. To allow for deviations beyond the ε-insensitive zone in the training data, slack variables ${\zeta }_i$ and ${\zeta }_i^{*}$ are introduced.

The constraints incorporate $\epsilon$, which determines the tolerance margin for estimation errors within the $\epsilon$-insensitive zone. This parameter significantly impacts the quantity of support vectors involved in constructing the regression model. Specifically, as $\epsilon$ increases, the number of selected support vectors tends to decrease, resulting in a sparser model representation. With a properly selected value of $\gamma$, a trade-off can be implemented among minimizing the training error and decreasing the model complexity, as captured by the regularization term ${‖\omega ‖}^2$.

The key approach to satisfying these constraints is to formulate a Lagrangian function based on the objective function presented in equation (8) and maximizing the resulted function with respect to the dual variables ${\alpha }_i$ and ${\alpha }_i^{*}$. Applied to our optimization problem, the Lagrangian for (8) is

subjected to ${\alpha }_i,{\alpha }_i^{*},{\zeta }_i,{\zeta }_i^{*}>0$. At this stage, Equation (8) is subject to minimization with respect to the primal variables ($\omega$, b, ${\zeta }_i$, and ${\zeta }_i^{*}$), while simultaneously being maximized with respect to the Lagrange multipliers ${\alpha }_i$ and ${\alpha }_i^{*}$. So, by taking the partial derivatives of the Lagrangian function L with respect to the primal variables and equating each to zero, the conditions for optimality are obtained, the necessary conditions for optimality are obtained. [31,32].

Equation (9) can be converted into its corresponding dual formulation, and the solution to this dual problem is given by:

$$f(x)={\displaystyle \sum _{i=1}^n({\alpha }_i-{\alpha }_i^{*})\cdot K(x_i,x)}+b$$

(11)

which is subject to

$$0\le {\alpha }_i\le \gamma ,\begin{array}{cc}& 0\le {\alpha }_i^{*}\le \gamma \end{array}$$

In equation (12), the bias term b is formulated based on the Lagrange multipliers and the kernel function.

At this stage, all parameters in equation (11) have been computed offline. Therefore, equation (11) is employed online to estimate the output f(x) for any given input x (wind speed), as illustrated in Figure 2.

4. Proposed FPNN with IPSO Control System

The FPNN controller demonstrates robust performance in handling system uncertainties, owing to its inherent parallel processing structure and adaptive online learning mechanism. Furthermore, subsequent sections elaborate on the network architecture, the online learning algorithm, and the convergence analysis of the FPNN.

4.1. Fuzzy Probabilistic Neural Network (FPNN)

The proposed fuzzy probabilistic neural network (FPNN) consists of five layers, each designed to perform a specific signal transformation toward the final control output in Figure 3. In addition, the subsequent sections provide a detailed exposition of the signal propagation process and the principal role of each layer within the network architecture.

1) Layer1-Input layer: The inputs and outputs of the nodes are defined as follows:

where denotes the ith input to the input layer, and N indicates the Nth iteration. In this paper, the inputs of the FPNN are and , which represent the tracking error and its derivative, respectively.

2) Layer2-Membership layer: In the FPNN architecture, the receptive field is modeled using a Gaussian function. The equation of Gaussian function are defined by:

where is the output of the jth node corresponding to the ith input variable; denotes the center of the Gaussian function at the jth node associated with the ith input variable, and represents its corresponding width.

3) Layer3-Probabilistic layer: The receptive field function is typically represented by a Gaussian function, and its mathematical formulation is given by:

where denotes the output of the kth node corresponding to the jth input variable; is the center of the Gaussian function, and is the base of the Gaussian function.

4) Layer4-Rule layer: Each node corresponds to a specific rule within the knowledge base. In accordance with the Mamdani inference mechanism, each node utilizes a t-norm operation to generate an inference set based on its associated fuzzy rule, as expressed in equation (20). Subsequently, probabilistic information is handled using Bayes’ theorem [7,8], under the assumption that the fuzzy membership grades constitute independent variables, as shown in equation (21). Based on this formulation, the input and output of each node in this layer are defined as follows:

where and are the input of rule layer; and are set to one; is the output of the rule layer.

5) Layer5-Output layer: The node outputs are expressed as:

where y(N) denotes the output of the FPNN, which simultaneously functions as the control effort of the proposed controller, and w_l indicates the connective weight among the Layer 4 and the Layer 5.

4.2. Online Supervised Learning and Training Process

The learning mechanism of the FPNN is fundamentally based on the recursive computation of the gradient vector, wherein each component corresponds to the partial derivative of a predefined energy function with respect to an individual network parameter. To derive the online learning algorithm within a supervised gradient descent framework, the error function E is initially defined as follows:

Subsequently, the learning algorithm is detailed as follows:

1) Layer5: The propagated error term is defined as follows:

and the weight updates are computed as:

where η₁ denotes the learning rate, controlling the step size of the updates. The connective weight w_l is adjusted using the following update rule:

2) Layer4: The error term to be propagated is given:

3) Layer2: The following expression defines the error term that is propagated through the network:

By applying the chain rule, the update equations for the center and width of the membership function are derived as follows:

where the factors η₂ and η₃ are the learning rates. The parameter (η₁, η₂, η₃) will be optimized using the IPSO algorithm. The center m_j and base σ_j of the membership function are update according to the following formula:

Furthermore, the exact computation of the system Jacobian, , is challenging due to the unknown dynamics of the control system. To address this issue, a delta adaptation law is employed, as described below [33]:

4.3. Convergence Analysis

To analyze the stability of the control system, a candidate Lyapunov function $V\left(k\right)$—serving as an energy-like measure for the tracking error—is initially proposed. To ensure asymptotic stability, the weights of the neural network are adaptively adjusted via a backpropagation algorithm, propagating error signals from the output layer toward the input layer. This adaptive update scheme is designed such that the time derivative of the Lyapunov function remains negative semi-definite $\Delta V\left(k\right)=V(k)-V(k-1)<0$. According to Lyapunov’s direct method [34], if the Lyapunov function is positive definite $V\left(k\right)>0$ and its time derivative is negative semi-definite $\Delta V\left(k\right)<0$, then the function $V(k)$ qualifies as a valid Lyapunov candidate for the error dynamics. Consequently, the output tracking error $e\left(k\right)$ of the FPNN asymptotically converges to zero as time approaches infinity. This theoretical framework ensures both the stability of the closed-loop control system and the convergence of the neural network learning process.

4.4. Adjustment of Learning Rates Using IPSO

To further enhance the online learning capability of the FPNN, a hybrid optimization technique—IPSO, which incorporates elements of genetic algorithms and time-varying particle swarm optimization—is employed in this study to tune the learning rates η₁, η₂ and η₃. In the IPSO algorithm, each particle updates its position by leveraging both its individual experience and the collective experiences of its neighbors. This update process accounts for the particle’s current velocity and position, as well as the best positions previously identified by itself and its neighboring particles [35].

Two pseudo-random sequences, γ₁~U(0,1) and γ₂~U(0,1) , are employed to simulate the stochastic behavior of the algorithm. For each dimension d, $L_i^d$, ${\mathit{pbest}}_i^d$ represent the current position and the personal best position of the particle, respectively. The velocity update rule, incorporating these elements, is defined in equation (32). Besides, the inertia weight w=0 is set. IPSO can reduce a parameter settings. The acceleration coefficients f₁ and f₂ can be modified by the following equations (32) and (33) [36,37,38]. These parameters are referred to as time-varying acceleration coefficients.

Here, $V_i^d$ and $L_i^d$ represent the current velocity and position of the particle, respectively, while N_max denotes the maximum number of iterations. f_1i and f_2i are the initial parameters setting. f_1j and f_2j are the final parameters setting.

Step 1: Definition of Basic Conditions

With $L_i^d=[L_i^1,L_i^2,L_i^3]$ for learning rates (η₁, η₂, η₃), the population size is defined as P=15. Each particle is defined with a dimensionality of d=3. The parameters are optimized within predefined lower and upper bounds.

Step 2: Initialization of Position and Velocity

Each particle’s initial position $L_i^d\left(N\right)$ and velocity $V_i^d\left(N\right)$ are randomly assigned within the boundaries of the defined search domain. Each particle's initial pbest is assigned as its current position, and the gbest is determined by selecting the best pbest among the entire group. The $L_i^d\left(N\right)$ are randomly generated by

where $U[{\eta }_{\min }^d,{\eta }_{\max }^d]$ represents a value sampled from a uniform distribution within the lower and upper bounds of the learning rate, η_min and η_max, respectively.

Step 3: Formulation of the fitness function

For each vector $L_i^d$, a corresponding fitness value must be assigned and evaluated. In this study, an appropriate fitness function is formulated to evaluate the fitness value, as defined below:

where FIT is the fitness value, abs(·) is the absolute function; 0.1 is added in the dominant part to avoid approaching infinity.

Step 4: Selection of pbest and gbests

Each particle $L_i^d$ retains its own fitness value and identifies the highest value encountered thus far as ${\mathit{pbest}}_i^d$. And the maximum vector in the population ${\mathit{pbest}}_i^d=[{\mathit{pbest}}_1^d,{\mathit{pbest}}_2^d,...{\mathit{pbest}}_p^d]$ is reached. Furthermore, in the first iteration, each particle $L_i^d$ is initially assigned to ${\mathit{pbest}}_i^d$, and the particle with the highest fitness value among all pbest values is designated as the global best gbest.

Step 5: Check the gbest for updates

If the position of the gbest particle remains unchanged over a predefined number of iterations, a crossover operation is applied between gbest and a chromosome from the Genetic Algorithm (GA). The reorganization of position and velocity are

where f₃ are the acceleration factor; rand() represents a random number uniformly distributed in the interval [0, 1]. l_parent and l_child are parent and child generation of position. s_parent and s_child are parent and child generation of velocity. λ represents the interpolation factor between the parent and offspring generations, defined as a random number uniformly distributed in the range [0, 1].

Step 6: Update position and velocity

The updated velocity is then incorporated into the current position of the particle to compute its subsequent position, as described by equations (35) and (38).

Step 7: Convergence check

Steps 3 to 6 are iteratively executed until a significant improvement in the gbest fitness value is observed, or a predefined maximum number of generations is reached. The final and highest fitness value ${\mathit{gbest}}_i^d$ obtained is considered the optimal learning rate (η₁, η₂, η₃) for the FPNN.

5. Case Studies and Simulation Results

The performance of the proposed FPNN-IPSO controller is assessed and benchmarked against the TSK fuzzy controller [3] and the conventional PI controller. All three control strategies—FPNN with IPSO, TSK fuzzy, and PI—were implemented and evaluated through simulation. To validate the effectiveness of the proposed controller, three case studies were conducted. The comparative performance outcomes of the controllers are illustrated in Figure 4, Figure 5 and Figure 6 and are summarized in Table 1. Each case study was carried out under a wind profile characterized by a volatile sinusoidal waveform, with a sampling interval of 5 msec for wind velocity measurements. The demonstrations were carried out to check effective of a VRB system used for the stable supply of wind energy. When one of the inverters in the microgrid is controlled, the voltage and frequency can be maintained within the predefined acceptable ranges. The power output of the inverters can always be kept in the rating working point. In other words, the AC microgrid system is robustly stable with respect to plug and play operation of WTG. The average power of PI is compared with that of FPNN-IPSO and TSK fuzzy algorithm.

The wind turbine generator system employed in the simulation is characterized by the following parameters:

5.1. FPNN with IPSO Algorithm

The proposed FPNN-IPSO controller is evaluated in the configuration illustrated in Figure 1. As shown in Figure 4(a), the shaft speed successfully tracks the maximum power point at or below the rated generator speed, demonstrating effective performance. The validation of the MPPT capability is presented in Figure 4(b), which also includes the wind speed profile and the dynamic discrepancy between the turbine power and generator power. Figure 4(c) displays the power coefficient , which remains close to its optimal value throughout the entire wind speed profile, indicating efficient MPPT operation. The SVR-based speed tracking mechanism accurately follows the actual speed with minimal deviation across the wind profile. Figure 4(d) illustrates the speed tracking error, which is approximately 0.25 rad/s. The average power output achieved is 271 W, representing a 9.71% improvement over the conventional PI controller. These results confirm that the proposed controller delivers superior transient and steady-state performance, as evidenced by the data presented in Figure 4.

5.2. TSK Fuzzy-Based Algorithm

By replacing the FPNN-IPSO controller with a TSK fuzzy-based algorithm, the simulation results illustrated in Figure 5 indicate an average power output of 259 W over a 50 sec interval, corresponding to a 4.85% improvement compared to the conventional PI controller. The power coefficient remains approximately at 0.4412. As depicted in Figure 5(a), the shaft speed effectively tracks the maximum power point while operating below the rated generator speed. Figure 5(b) presents the validation of the maximum power tracking control, while Figure 5(c) displays the power coefficient, confirming the satisfactory performance of the TSK fuzzy-based controller in maintaining efficient power conversion. The results in Figure 5 show that magnitude of the rotation speed and power coefficient are relatively large. Compared to the PI controller, the TSK fuzzy controller offers improved transient and steady-state response, although it remains slightly less accurate than the proposed method.

5.3. PI Controller

When the FPNN-IPSO controller is replaced with the conventional PI controller, the simulation results are shown in Figure 6. The average power output achieved over the same 50-second period is 247 W. As illustrated in Figure 6(a), the shaft speed demonstrates the system's response under PI control, while Figure 6(b) provides the validation of the MPPT control. Figure 6(c) displays the power coefficient, which consistently remains close to the optimal value of approximately 0.4412. These results indicate that while the PI controller maintains reasonably effective tracking performance, its efficiency is comparatively lower than that of the FPNN-IPSO and TSK fuzzy-based controllers. The conventional PI controller achieves acceptable transient behavior; however, its closed-loop settling time is significantly longer compared to the FPNN with IPSO and TSK fuzzy methods.

5.4. Performance Comparison

The simulation-based performance comparison of various control strategies underscores the significance of MPPT across a wide range of wind speeds, as summarized in Table 1. This table presents a comparative analysis including average output power, the maximum error in the power coefficient, and the percentage increase in power achieved by each method. The proposed control algorithm demonstrates superior performance compared to existing methods reported in [17,25], and [39]. For instance, the maximum error in the power coefficient reaches approximately 23% in [17], while the maximum power deviation in [25] is around 7%. In contrast, the proposed method achieves a marked reduction in these errors, indicating enhanced tracking capability and control precision. Figure 7 illustrates the efficiency performance of the three control algorithms. With the implementation of the proposed strategy, the generator efficiency—defined from mechanical input to AC output—exceeds 86.11%. Additionally, when accounting for converter losses, the overall system efficiency—from mechanical input to DC output—remains above 81% under medium to high wind speed conditions. These results confirm that each system component can independently respond to dynamic conditions using a decentralized Plug-and-Play strategy, ensuring coordinated and stable operation of the microgrid.

6. Conclusion

This study proposes and investigates three distinct control algorithms specifically developed for a PMSG integrated within a variable-speed WECS. Under a varying wind speed profile, the proposed system is capable of performing MPPT. The speed control mechanism generates the generator torque reference, which is executed through an inner-loop current controller. Notably, all three MPPT algorithms proposed in this work are implemented without the use of wind speed sensors.

Among the methods evaluated, the conventional proportional–integral (PI) controller is shown to operate near the optimal point under certain conditions. However, the PI-based approach may exhibit limitations in systems characterized by strong nonlinearities or rapidly changing reference trajectories, resulting in suboptimal control performance.

Simulation results validate the efficacy of the proposed FPNN combined with IPSO control scheme in both trajectory forecasting and real-time power regulation. Experimental verification further demonstrates that the FPNN-IPSO-based MPPT strategy demonstrates high accuracy in both trajectory prediction and target tracking, ensuring precise alignment with the optimal operating point under varying wind conditions. Furthermore, the control strategy effectively regulates power output to meet load demands and enhance system stability. Additionally, it provides coordinated control across multiple components, contributing to microgrid optimization. The proposed method exhibits robust dynamic performance and maintains system stability even in the presence of parameter uncertainties, thereby ensuring reliable operation under diverse operating conditions.