A Uniﬁed Active Frequency Regulating and Maximum Power Point Tracking Strategy for Photovoltaic Sources

: In order to optimize the extraction of solar energy, photovoltaic sources are commonly operated under the control of the so-called maximum power point (MPPT) strategy. However, as the rate of PV installations increases explosively, traditional MPPT algorithms may cause problems such as frequency deviation and power ﬂuctuations, making system frequency stability a challenge due to the inherent intermittent and stochastic nature of PVs. Consequently, in order to reduce the investment and maintenance costs of storage systems, innovative control is expected for PV sources to provide ancillary services for the system, especially for weak systems such as microgrids. In this paper, a novel active power control (APC) strategy, based on characteristic curve ﬁtting, is proposed to ﬂexibly regulate the PV output power. The transient process performance and robustness of the system are improved with the proposed APC strategy. In conjunction, an f – P droop mechanism is designed to provide a frequency regulating (FR) service for the AC microgrid. The comprehensive control strategy uniﬁes the FR function with the traditional MPPT function in a single control structure, allowing the PV source to operate either in the MPPT mode when the system frequency is nominal or in FR mode when the frequency exceeds it. The transition between MPPT and FR is autonomous and fully decentralized, which improves the PV generation efﬁciency as well as ensuring generation fairness among different parallel PV sources. Importantly, the proposed control strategy does not require any internal bundled energy within the PV generation system to achieve FR capability, but it effectively collaborates with the system-level energy storage system, thus reducing the necessary battery capacity. A detailed dynamic model of a PV generation system is constructed to validate the feasibility and effectiveness of the proposed control strategy.


Introduction
Solar energy, due to its relative low installation and maintenance costs and widespread resource distribution [1], has emerged as the predominant resource within renewable power systems. As more and more PV sources are being installed, intermittent and stochastic characteristics of irradiation are increasingly challenging and threatening to system stability control and safe operation. This problem is especially pronounced in weak systems such as AC microgrids and remote insular or mountainous power systems [2,3]. Traditionally, we use the energy storage system for frequency regulation and leave the maximum power point tracker to extract the highest available power from the PV array. As a result, PV sources do not inherently contribute to the system frequency regulation. Hence, from the perspective of the large power grid system, the projected growth of PV systems will put stress on other conventional generators and further weaken the power system stability [4]. Moreover, in future, PV systems may replace a significant part of rotational generation puts forward a universal controller consisting of a fast MPPT (FMPPT) controller and a slow MPPT (SMPPT) controller. The transition between controllers may not be seamless, potentially affecting system stability and inducing unwanted oscillations. A novel Newton quadratic interpolation-based APC algorithm is presented in Ref. [4] and utilized for system frequency regulation in Ref. [21]. The algorithm has relatively poor robustness and might cause transient problems when the irradiation or loads change rapidly and dramatically since the quadratic curve is quite different from the PV P-U curve.
To improve on methods presented in previous literature, a novel FR strategy for PV sources in AC microgrids with a traditional hierarchical control framework is introduced in this paper. With the proposed strategy, PV sources can either provide FR service when necessary or operate in MPPT mode to absorb as much of the solar resource as possible. The PV sources adaptively deliver primary FR control based on the system frequency deviation, and its FR capability naturally adjusts with irradiation variations. The proposed scheme unifies MPPT and FR within a single controller, eliminating the transition between different operating modes. Furthermore, similar to the commonly used droop control methods, it is fully decentralized and suitable for both distributed PV sources and large-scale PV plants. Importantly, the proposed control strategy does not require any internal bundled energy within the PV generation system to achieve FR capability, but it effectively collaborates with the system-level energy storage system, thus reducing the necessary battery capacity. A comparison between the proposed control strategy and the existing frequency regulation methods for PV sources is provided in Table 1. Ref. [8] Delta power control low medium √ Refs. [9][10][11] DC-link capacitor control high small √ Refs. [12][13][14] Deloading control low large √ Refs. [15,16] dp/dv control low √ √ large Ref. [18] Reference power tracking high large √ Ref. [18] FMPPT and SMPPT high large √ Refs. [19,20] Quadratic interpolation high √ √ large √ The proposed method Characteristic curve fitting high √ √ large √ √ " √ " in the table indicates that the method has the corresponding feature.
The remaining sections of this paper are organized as follows: Section 2 offers a general introduction to the proposed strategy and control structure. Section 3 elaborates on the specialized APC algorithm based on characteristic curve fitting. Section 4 examines the P-f droop control for PV sources and discusses its adaptability. Section 5 presents simulation tests, demonstrating the feasibility and effectiveness of the proposed strategy through the results. Finally, Section 6 concludes the paper.

Principles of the Proposed Control Strategy
The main features of the proposed PV control strategy are as follows: (1) The PV works in the FR mode to participate in the system frequency regulation when the system frequency surpasses its nominal value and the irradiance is sufficient. This control strategy aims to deliver frequency regulating services for the connected system, to cooperate with other synchronous sources, but not to become the sole grid-forming generator. (2) The PV works in the MPPT mode when the system frequency is at its nominal value or the irradiance is insufficient such as shadowing time. The insufficient irradiance may lead to possible loss of FR capability of PV sources but the system is still governed by The concept of droop is widely utilized for FR and multiple-source cooperation in an AC microgrid, irrespective of whether it is grid-connected or islanded. Typically, the role of P-f type source is assumed by the energy storage and grid-connected interlink converters. Since distributed generators such as PV sources cannot participate in the primary FR such as P-f type sources, load changes and power fluctuations may induce significant frequency deviations. This situation becomes more problematic for weak systems such as AC microgrids that have a high portion of PV sources.
The mathematical model of the PV array current at nominal irradiance and temperature is given by [17]: where I PV is the PV output current. N P and N s are the number of parallel-and seriesconnected PV modules and each module is composed of N PV cells. I sc,n and V oc,n are the short-circuit current and open-circuit voltage of the PV module at nominal condition (1000 W/m 2 , 25 • C). V PV is the operating voltage of the PV array and V t = NkT/q is the thermal voltage of the PV array with k being the Boltzmann constant, q being the electron charge and T being the temperature (Kelvin) of the PV cells. a is the ideality constant of the equivalent diode. The detailed explanations on PV modeling may be found in Refs. [22,23]. Then the output power of the PV array is derived from (1): The nonlinear characteristic between the PV output power, P PV , and its terminal voltage, V PV , is illustrated in Figure 1. As shown, the P-U curve can be divided into two segments: (1) the uphill segment and the downhill segment, determined by the sign of dP PV /dV PV . Operating within these segments reveals contrasting control characteristics: the PV output increases when the PV voltage rises in the uphill segment, while it decreases in the downhill segment when the PV voltage increases. This contrast complicates the design of frequency/voltage regulators when relying solely on conventional controller adaptations. Consider a typical two-level control structure for frequency regulation, as shown in Figure 2. If an abrupt change in irradiance or system load shifts the system frequency away from its reference, and the PV sources need to reduce their output to mitigate the frequency deviation, then level 2 should modulate the PV output power according to the P-U curve in Figure 2. Yet, Figure 2 also indicates that the V PV adjustment direction varies between the uphill and downhill segments. To decrease PV output in the uphill segment, a reduced V PV is required, while the opposite is true for the downhill segment. For effective frequency regulation with linear controllers, the PV operating range must be constrained to either the uphill or downhill segments. Hence, the maximum power point, particularly the value of V MPP , becomes critical for advanced PV applications such as frequency regulation.  To solve the aforementioned issues, a three-level hierarchical control strategy is designed for PV sources to enable them to contribute to primary FR. Figure 3 presents the control diagram of the proposed FR strategy with a typical single-stage PV generation system. As depicted, the overall strategy comprises three control levels: PV voltage regulation (level 1), active power control (level 2) and frequency droop control (level 3).
Level 3, the outermost control loop of the proposed control strategy, primarily focuses on the FR, forming the central part of the proposed cascaded controller. A frequency droop curve is employed for PV sources to modulate their output power based on the system frequency deviation. If the system frequency is maintained within an acceptable range, PV sources will autonomously switch to the MPPT mode from FR mode. This transition is governed by the system frequency and irradiation levels, which define the FR capability of PV sources. The input to level 3 is the frequency deviation, ∆ (the difference between its actual value and nominal value), and its output is the PV power reference. Essentially, level 3 determines the PV output power reference according to the system frequency deviation to provide FR service, while the inner control loops track the power references.
Level 2 identifies the accurate PV operating voltage, , that corresponds with the power reference, , provided by level 3, as controlling PV operating voltage is a  To solve the aforementioned issues, a three-level hierarchical control strategy is designed for PV sources to enable them to contribute to primary FR. Figure 3 presents the control diagram of the proposed FR strategy with a typical single-stage PV generation system. As depicted, the overall strategy comprises three control levels: PV voltage regulation (level 1), active power control (level 2) and frequency droop control (level 3).
Level 3, the outermost control loop of the proposed control strategy, primarily focuses on the FR, forming the central part of the proposed cascaded controller. A frequency droop curve is employed for PV sources to modulate their output power based on the system frequency deviation. If the system frequency is maintained within an acceptable range, PV sources will autonomously switch to the MPPT mode from FR mode. This transition is governed by the system frequency and irradiation levels, which define the FR capability of PV sources. The input to level 3 is the frequency deviation, ∆ (the difference between its actual value and nominal value), and its output is the PV power reference. Essentially, level 3 determines the PV output power reference according to the system frequency deviation to provide FR service, while the inner control loops track the power references.
Level 2 identifies the accurate PV operating voltage, , that corresponds with the power reference, , provided by level 3, as controlling PV operating voltage is a To solve the aforementioned issues, a three-level hierarchical control strategy is designed for PV sources to enable them to contribute to primary FR. Figure 3 presents the control diagram of the proposed FR strategy with a typical single-stage PV generation system. As depicted, the overall strategy comprises three control levels: PV voltage regulation (level 1), active power control (level 2) and frequency droop control (level 3).
Level 3, the outermost control loop of the proposed control strategy, primarily focuses on the FR, forming the central part of the proposed cascaded controller. A frequency droop curve is employed for PV sources to modulate their output power based on the system frequency deviation. If the system frequency is maintained within an acceptable range, PV sources will autonomously switch to the MPPT mode from FR mode. This transition is governed by the system frequency and irradiation levels, which define the FR capability of PV sources. The input to level 3 is the frequency deviation, ∆ f (the difference between its actual value and nominal value), and its output is the PV power reference. Essentially, level 3 determines the PV output power reference according to the system frequency deviation to provide FR service, while the inner control loops track the power references.
Level 2 identifies the accurate PV operating voltage, V re f PV , that corresponds with the power reference, P re f PV , provided by level 3, as controlling PV operating voltage is a common method for managing PV sources' output power. When P re f PV is less than the maximum PV , provided by level 2, the PV output power matches the required value dictated by level 3. Notably, level 3 generates suitable power generation references according to the system frequency deviation to achieve FR capability rather than directly regulating the inverter's output frequency.
Electronics 2023, 12, x FOR PEER REVIEW 6 of 18 common method for managing PV sources' output power. When is less than the maximum available power point (MAPP), a that corresponds to an operating point below the MAPP needs to be found. When exceeds the MAPP, (the PV operating voltage that corresponds to the MAPP) itself should be given out. The inherent nonlinearity in the PV's P-U curve makes the direct determination of challenging. Generally, there are three categories of solutions: (1) indirect measurement methods reliant on irradiation and temperature sensors; (2) searching methods such as perturb and observe (P&O) schemes and (3) estimation methods. Considering the drawbacks of additional measurement sensors and the relatively low convergence rate of P&O methods, this paper enhances conventional estimation methods by introducing a novel characteristic curve. The details of the proposed active power control algorithm are discussed in the subsequent section.
Level 1, the innermost control loop, employs a traditional dual-loop control structure to regulate the terminal voltage of the PV array. The outer loop at level 1 controls the PV array voltage, , and the reactive power, while the inner loop manages the current. Under the control of level 1, when the PV array voltage, , equals , provided by level 2, the PV output power matches the required value dictated by level 3. Notably, level 3 generates suitable power generation references according to the system frequency deviation to achieve FR capability rather than directly regulating the inverter's output frequency.  Figure 3. Structure of the PV source with proposed control system.

The Proposed APC Algorithm
Two complex issues arise when considering flexible active power control of PV sources: (1) The P-U characteristic alters with changing ambient conditions (irradiance and temperature) changes, implying that the PV operating point ( , ) varies under different environmental conditions; (2) a single photovoltaic output power corresponds to two distinct operating voltages, except for the MPP. The key point for the PV power control loop can be expressed as follows: For an arbitrary power dispatching reference, : If is lower than the MPP, then find the appropriate operating voltage (whether lower than the MPP voltage or higher, but with consistent maintenance) to equate with . While if exceeds MPP, then the MPP voltage should be found. In summary, the main task for PV power control is to find a PV operating voltage, * , that satisfies:

The Proposed APC Algorithm
Two complex issues arise when considering flexible active power control of PV sources: (1) The P-U characteristic alters with changing ambient conditions (irradiance and temperature) changes, implying that the PV operating point (U PV , P PV ) varies under different environmental conditions; (2) a single photovoltaic output power corresponds to two distinct operating voltages, except for the MPP. The key point for the PV power control loop can be expressed as follows: For an arbitrary power dispatching reference, P re f : If P re f is lower than the MPP, then find the appropriate operating voltage (whether lower than the MPP voltage or higher, but with consistent maintenance) to equate P PV with P re f . While if P re f exceeds MPP, then the MPP voltage should be found. In summary, the main task for PV power control is to find a PV operating voltage, U * , that satisfies: where P MPP PV is the real-time maximum PV power, P re f is the power dispatching reference, G is the irradiance and T is the temperature.
As stated before, it is hard to obtain the function f pv to solve the f pv U * , G, T = P re f for U * . The curve fitting and interpolation iteration methods are usually used to solve these problems in numerical calculations, for example, the quadratic interpolation used in [4]. Regarding previous power curtailment algorithms in literature, the control performance is usually unsatisfactory due to sensor dependance, computation complexity or algorithm robustness. Consequently, a characteristic curve fitting-based algorithm is introduced in this section. The fundamental concept involves utilizing a curve resembling the PV's characteristic to fit the P-U curve, and then iterating until it converges on a U PV value that corresponds with the given P re f PV . By enhancing the similarity between the real PV characteristic and the characteristic fitting curve, the robustness and convergence rate of the algorithm can be improved.
Considering two intersections with U-axis and the maximum point in the P-U curve of a PV array, the proposed characteristic fitting curve is given as: where a 1 , a 0 , b 1 and b 0 are the fitting coefficients and f is the proposed characteristic fitting curve function.
In the iterative process, only one point is sampled per iteration. Consider step k in the iterative process and assume that at least four points on the P-U curve have been sampled.
The four sampling points are denoted as The primary strategy of this method employes the fitted curve, which intersects the four sampled points, to emulate the real P-U curve and by continuous iteration, to make it converge locally to the target operating point, which is either the MPP or the power dispatch point. Once at least four sampling points are ready, the following equations are derived from Equation (4): where V k PV and P k PV with superscript k represent the PV voltage and output power in the kth iterative step, respectively. Equation (5) in then rewritten in the matrix form: where The fitting coefficients can then be solved by s = A −1 ·B or by the least square method if more than four sampling points are kept.
Once the valuables in s are obtained, substitute P re f into Equation (4) and rearrange it to obtain: Electronics 2023, 12, 3467 8 of 18 which can be regarded as a quadratic equation with V PV being the argument. This is defined as: where ∆ < 0 means that P re f has no intersections with the fitting curve, which indicates that the power dispatch reference, P re f , is higher than the MPP. Under this circumstance, the PV source should operate in MPPT mode as previously discussed and the maximum point of the fitting curve is employed for the next iteration step. The maximum point can be obtained by solving the zero point of the derivative of Equation (4). The MPP voltage can be approximated as: where V k pv_re f with superscript k represents the PV voltage reference for the kth iterative step. ∆ ≥ 0 means that P re f is less than or equal to the MPP and the PV source should now operate in the APC mode by regulating its output power to P re f . There are two solutions to Equation (7). The uphill segment solution (denoted by UHSS) is derived as: whereas the downhill segment solution (denoted by DHSS) is derived as: In [24,25], the DHSS is regarded as a more preferrable solution for PV voltage regulation for its better control performance in many scenarios. The focus of this paper is on the versatility and simplicity of the proposed APC algorithm and the DHSS is chosen as an example. It is still worthy of noting that the proposed algorithm applies to both DHSS and UHSS and that either can be selected.

Discuss on the Properties under Large Disturbances
As previously mentioned, two potential solutions exist for when the PV system operates under the MPP: one in the uphill segment and the other in the downhill segment. Ref. [26] provides a comprehensive analysis and modeling to reveal the intrinsic characteristics of these two different operating points of a PV generation system. To ease the analysis, the operating range was divided into four regions, namely: (1) the current-source region, (2) power region I, (3) power region II and (4) the voltage-source region, as shown in Figure 4, based on their respective output characteristics. The conclusion is that power regions I and II have demonstrated superior stability and controllable performance compared to the current-source and voltage-source regions. In the current-source and voltage-source regions, the phase margin decreases significantly, which needs to be avoided by controllers to prevent system instability.
Moreover, when the reference power output is exceptionally low, the corresponding PV array voltage of the LHSS may descend below the peak voltage of the AC side of the inverter, thereby leading to inverter failure [27]. For the RHSS, such a problem will never occur. Nevertheless, it is crucial to ensure that the PV operating point does not fall into the voltage-source region. When the reference power changes abruptly, the RHSS may have faster convergence rate than that of the LHSS. This can be explained intuitively in that the downhill segment of the P-U curve is steeper than the uphill segment. To mitigate the instability caused by large disturbances, the following saturation function can be used to where V   Moreover, when the reference power output is exceptionally low, the corresponding PV array voltage of the LHSS may descend below the peak voltage of the AC side of the inverter, thereby leading to inverter failure [27]. For the RHSS, such a problem will never occur. Nevertheless, it is crucial to ensure that the PV operating point does not fall into the voltage-source region. When the reference power changes abruptly, the RHSS may have faster convergence rate than that of the LHSS. This can be explained intuitively in that the downhill segment of the P-U curve is steeper than the uphill segment. To mitigate the instability caused by large disturbances, the following saturation function can be used to keep within an accepted range [ , ] (in the aforementioned power region I and power region II): where and are the predefined constants. The can be set as the maximum between the peak voltage of the AC side of the inverter and the upper boundary of the current-source region with a proper margin, whereas the can be set as the lower boundary of the voltage-source region to ensure the stable operation of the PV source.

PV Frequency Droop Control
Frequency droop control is commonly employed for primary FR and multi-source coordination (power sharing) in AC microgrids. For a traditional frequency droop type source in AC microgrids, its output current is determined based on the system frequency deviation. Traditionally, these frequency droop concepts cannot be applied to PV sources due to their lack of flexible power adjustment techniques. However, with the implementation of the APC algorithm presented in Section 3, this limitation is surmounted. A P-f type droop function, designed for PV frequency droop control, generates the reference PV output power in accordance with the deviation of the system frequency from its nominal

PV Frequency Droop Control
Frequency droop control is commonly employed for primary FR and multi-source coordination (power sharing) in AC microgrids. For a traditional frequency droop type source in AC microgrids, its output current is determined based on the system frequency deviation. Traditionally, these frequency droop concepts cannot be applied to PV sources due to their lack of flexible power adjustment techniques. However, with the implementation of the APC algorithm presented in Section 3, this limitation is surmounted. A P-f type droop function, designed for PV frequency droop control, generates the reference PV output power in accordance with the deviation of the system frequency from its nominal value [21]: P re f where P re f PV is the PV generation reference, (P nom , f nom ) is the nominal operating point, f sys is the system frequency and m is the droop coefficient, which indicates the change in PV output power for per unit voltage deviation, and can be designed as: According to Equation (13), P re f PV is calculated without consideration of the limit of maximum available power output, P MPP PV . When P re f PV reaches P MPP PV , it takes the value of P MPP PV and consequently loses the FR capacity. Figure 5 shows the P-f droop characteristics of the proposed PV frequency droop strategy. The limiting frequency, f lim , is defined as the frequency at which P nom − m· f sys − f nom is equal to the real-time P MPP PV : According to Equation (13), is calculated without consideration of the limit of maximum available power output, . When reaches , it takes the value of and consequently loses the FR capacity. Figure 5 shows the P-f droop characteristics of the proposed PV frequency droop strategy. The limiting frequency, , is defined as the frequency at which • is equal to the real-time : , It is apparent that the of the droop characteristic corresponds with and thus correlates with the irradiation and temperature. Specifically, the PV source generates if f while it reduces if f . The proposed strategy thus adapts to various irradiation and temperature conditions. If the irradiation levels are relatively high, the FR capability of the PV source is enhanced while it is diminished under insufficient irradiance. For example, Figure 5 illustrates that at 800 W/m , the PV sources activate FR only when exceeds 50.425 Hz (which is the at 800 W/m ), compared to 50 Hz at 1000 W/m .
In FR mode, a designated portion of the output power from PV sources is adjusted in response to system frequency deviations, based on the established f-P droop curve. This power modulation is determined by the droop coefficient, m, and the nominal power reference, . The coefficient m dictates the degree of power adjustment per frequency variation, while indirectly establishes the activation frequency of the FR for varying irradiance levels. Taking an 85 kW (1000 W/m ,25 °∁) PV source as an example, the concrete FR characteristics can be found in Table 2.  It is apparent that the f lim of the droop characteristic corresponds with P MPP PV and thus correlates with the irradiation and temperature. Specifically, the PV source generates P MPP PV if f < f lim while it reduces P PV if f > f lim . The proposed strategy thus adapts to various irradiation and temperature conditions. If the irradiation levels are relatively high, the FR capability of the PV source is enhanced while it is diminished under insufficient irradiance. For example, Figure 5 illustrates that at 800 W/m 2 , the PV sources activate FR only when f sys exceeds 50.425 Hz (which is the f lim at 800 W/m 2 ), compared to 50 Hz at 1000 W/m 2 .
In FR mode, a designated portion of the output power from PV sources is adjusted in response to system frequency deviations, based on the established f -P droop curve. This power modulation is determined by the droop coefficient, m, and the nominal power reference, P nom . The coefficient m dictates the degree of power adjustment per frequency variation, while P nom indirectly establishes the activation frequency of the FR for varying irradiance levels. Taking an 85 kW (1000 W/m 2 , 25 • C) PV source as an example, the concrete FR characteristics can be found in Table 2. Theoretically, P nom can be chosen based on the system requirements. When P nom is high, PV sources are unable to offer frequency regulation support at "low" (relatively low, but still higher than the nominal value) system frequency. Conversely, when P nom is low, PV sources have augmented FR capability, potentially sacrificing some benefits under high irradiation. Hence, there is a trade-off between economic benefits and FR capability. For weak systems such as a small-scaled, islanded AC microgrid with high PV penetration, which inherently has a weak FR capability, more support from PV sources is necessary, thus P nom should take a low-level value.
Importantly, the proposed adaptive control scheme is for efficiency promotion. In fact, if the irradiance is not sufficient, such as during shadowing time, the PV sources will try their best to harvest from the solar resource (work in MPPT mode) to enhance the system stability. The insufficient irradiation may lead to possible loss of frequency regulation capability of PV sources but the system is still governed by system-level storage or other synchronous sources. The proposed control strategy for PV sources aims to deliver frequency regulating services for the connected system, to cooperate with other synchronous sources, but not to become the only grid-forming power sources. This does not need bundled storage [28,29] in the PV generation system, but the system-level storage or synchronous sources are still needed. The service provided by PV sources may contribute to the reduction in the energy storage capacity required by the system, which may consequently reduce the installation and maintenance costs while enhancing the system operational reliability.

Simulation Results
A comprehensive dynamic model of an AC microgrid, composed of two PV sources, a frequency droop type energy storage and two AC loads, is established using Matlab/Simulink R2021b for simulation tests. The system configuration of the studied AC microgrid is shown in Figure 6.
high, PV sources are unable to offer frequency regulation support at "low" (relatively low, but still higher than the nominal value) system frequency. Conversely, when is low, PV sources have augmented FR capability, potentially sacrificing some benefits under high irradiation. Hence, there is a trade-off between economic benefits and FR capability. For weak systems such as a small-scaled, islanded AC microgrid with high PV penetration, which inherently has a weak FR capability, more support from PV sources is necessary, thus should take a low-level value. Importantly, the proposed adaptive control scheme is for efficiency promotion. In fact, if the irradiance is not sufficient, such as during shadowing time, the PV sources will try their best to harvest from the solar resource (work in MPPT mode) to enhance the system stability. The insufficient irradiation may lead to possible loss of frequency regulation capability of PV sources but the system is still governed by system-level storage or other synchronous sources. The proposed control strategy for PV sources aims to deliver frequency regulating services for the connected system, to cooperate with other synchronous sources, but not to become the only grid-forming power sources. This does not need bundled storage [28,29] in the PV generation system, but the system-level storage or synchronous sources are still needed. The service provided by PV sources may contribute to the reduction in the energy storage capacity required by the system, which may consequently reduce the installation and maintenance costs while enhancing the system operational reliability.

Simulation Results
A comprehensive dynamic model of an AC microgrid, composed of two PV sources, a frequency droop type energy storage and two AC loads, is established using Matlab/Simulink R2021b for simulation tests. The system configuration of the studied AC microgrid is shown in Figure 6.  The ESS maintains the frequency and voltage of the AC microgrid. The converter of the ESS employs the conventional P-f droop control and the droop coefficient is 0.01 Hz/kW. The nominal frequency and active power are 50 Hz and 0 kW, respectively. A traditional dual loop control structure is used, with voltage and current regulator gains set to K p_v_bat = 1, K i_v_bat = 100 and K p_i_bat = 3, K i_i_bvat = 200, respectively.
Two individual PV sources exist within the studied AC microgrid to demonstrate the power sharing coordination under the proposed control strategy. Each PV source delivers a nominal output power of 85 kW at 1000 W/m 2 and 25 • C. If the PV sources use the proposed P-f droop control, the nominal frequency is set as 50 Hz with a droop coefficient of 40,000 W/Hz. However, if the PV sources use the MPPT mode, the reference power value is set at 100 kW. The sample interval for the proposed APC algorithm is 0.02 s, with the PV source's dual loop gains set to K p_v_pv = 0.1, K i_v_pv = 1 and K p_i_pv = 3, K i_i_pv = 200, respectively.
In Section 5.1, Load 1 = 85 kW and Load 2 = 85 kW while in Section 5.2, Load 1 = 70 kW and Load 2 = 85 kW. The electrical parameters of the simulation system are summarized in Table 3. Table 3. Electrical parameters of the simulation system.

Effect of PV Frequency Droop Control
This scenario is mainly to study the effect of PV frequency droop control. Load 2 (85 kW), accounting for approximately 50% of the system's total load, is offline at 7 s. Comparative simulation tests are conducted, comparing PV sources operating under the conventional MPPT control with those under FR control. The dynamics of the storage output powers, PV output powers and the system frequencies under different control strategies are depicted in Figures 7-9, respectively.
Prior to the disconnection of Load 2, both sets of PV sources, under different control schemes, operate in the MPPT mode for both cases. The irradiation is set as 1000 W/m 2 and the output power from each PV source is equal. After 15 s, under traditional MPPT control, all power variations in the system (Load 2) are balanced by the system storage, with charging power rising from 0 to 82 kW. The system frequency increases from 50 to over 50.8 Hz. In contrast, under the proposed frequency droop strategy, both PV sources share the power imbalance burden along with the system storage in response to a sudden load change. After 15 s, each PV source gradually reduces their output power by 17 kW and system storage only raises its charging power from 0 to 47 kW. The system frequency rises to approximately 50.5 Hz, displaying superior frequency deviation mitigation compared to traditional methodologies. Electronics 2023, 12, x FOR PEER REVIEW 13 of 18

Effect of PV Different Irradiance
This scenario primarily investigates the effect of different irradiation levels. In this scenario, Load 2 (85 kW), which is 55% of the total system load in Section 5.2, is disconnected at the beginning of the test and recovers at 15 s. Similar to Section 5.1, there are still two cases, one under the MPPT scheme and the other under the proposed control strategy. Both cases are tested with three different irradiation levels (1000 W/m 2 , 800 W/m 2 and 600 W/m 2 , respectively). The dynamics of the storage output powers, PV output powers and the system frequencies are depicted, respectively, in  Both cases are tested with three different irradiation levels (1000 W/m , 800 W/m and 600 W/m , respectively). The dynamics of the storage output powers, PV output powers and the system frequencies are depicted, respectively, in  During the initial 15 s, the system frequency exceeds for 1000 W/m and 800 W/m and the PV sources operate in the FR mode under these two irradiation levels However, for 600 W/m , as the system frequency falls below the for 600 W/m , the operation mode defaults to MPPT. Upon reconnection of Load 2, the system frequency decreases, prompting the PV sources to enhance their output power in accordance with the droop curve. As depicted in Figure 12, under low irradiation (G = 800 W/m and G = 600 W/m in this case), the FR capability of PV sources diminishes, causing an adaptive switch to the MPPT mode. Consequently, greater frequency drops are observed. The sim ulation results validate that PV sources under higher solar irradiation have stronger FR capability. If the system's original FR capability is insufficient, it is preferable to take a lower in order to enhance the PV FR capacity.   Both cases are tested with three different irradiation levels (1000 W/m , 800 W/m and 600 W/m , respectively). The dynamics of the storage output powers, PV output powers and the system frequencies are depicted, respectively, in Figures 10-12.
During the initial 15 s, the system frequency exceeds for 1000 W/m and 800 W/m and the PV sources operate in the FR mode under these two irradiation levels.
However, for 600 W/m , as the system frequency falls below the for 600 W/m , the operation mode defaults to MPPT. Upon reconnection of Load 2, the system frequency decreases, prompting the PV sources to enhance their output power in accordance with the droop curve. As depicted in Figure 12, under low irradiation (G = 800 W/m and G = 600 W/m in this case), the FR capability of PV sources diminishes, causing an adaptive switch to the MPPT mode. Consequently, greater frequency drops are observed. The simulation results validate that PV sources under higher solar irradiation have stronger FR capability. If the system's original FR capability is insufficient, it is preferable to take a lower in order to enhance the PV FR capacity.   During the initial 15 s, the system frequency exceeds f lim for 1000 W/m 2 and 800 W/m 2 and the PV sources operate in the FR mode under these two irradiation levels. However, for 600 W/m 2 , as the system frequency falls below the f lim for 600 W/m 2 , the operation mode defaults to MPPT. Upon reconnection of Load 2, the system frequency decreases, prompting the PV sources to enhance their output power in accordance with the droop curve. As depicted in Figure 12, under low irradiation (G = 800 W/m 2 and G = 600 W/m 2 in this case), the FR capability of PV sources diminishes, causing an adaptive switch to the MPPT mode. Consequently, greater frequency drops are observed. The simulation results validate that PV sources under higher solar irradiation have stronger FR capability. If the system's original FR capability is insufficient, it is preferable to take a lower P nom in order to enhance the PV FR capacity.

Test under Real Operational Scenario
In this scenario, the performance of the proposed control strategy is studied under a real operational scenario. The PV system is simulated over a 30 s period in a typical cloudy day (real irradiance data from [4,30] are used, as shown in Figure 13). The sharp variation of the solar irradiance may be because of a passing cloud, or the shade of a building. Figures 14 and 15 show the time responses of PV output power and the system frequency under both the MPPT mode and the proposed FR mode, respectively. It is proved that the PV power output can follow the change of irradiance immediately, although some frequency and power spikes can be observed when the irradiance changes too sharply. When the PV system operates in the FR mode, and if the irradiance is abundant, the FR service could significantly mitigate the frequency deviation, otherwise, the PV system adaptively switches to the MPPT mode.

Test under Real Operational Scenario
In this scenario, the performance of the proposed control strategy is studied under a real operational scenario. The PV system is simulated over a 30 s period in a typical cloudy day (real irradiance data from [4,30] are used, as shown in Figure 13). The sharp variation of the solar irradiance may be because of a passing cloud, or the shade of a building. Figures 14 and 15 show the time responses of PV output power and the system frequency under both the MPPT mode and the proposed FR mode, respectively. It is proved that the PV power output can follow the change of irradiance immediately, although some frequency and power spikes can be observed when the irradiance changes too sharply. When the PV system operates in the FR mode, and if the irradiance is abundant, the FR service could significantly mitigate the frequency deviation, otherwise, the PV system adaptively switches to the MPPT mode.

Test under Real Operational Scenario
In this scenario, the performance of the proposed control strategy is studied under a real operational scenario. The PV system is simulated over a 30 s period in a typical cloudy day (real irradiance data from [4,30] are used, as shown in Figure 13). The sharp variation of the solar irradiance may be because of a passing cloud, or the shade of a building. Figures 14 and 15 show the time responses of PV output power and the system frequency under both the MPPT mode and the proposed FR mode, respectively. It is proved that the PV power output can follow the change of irradiance immediately, although some frequency and power spikes can be observed when the irradiance changes too sharply. When the PV system operates in the FR mode, and if the irradiance is abundant, the FR service could significantly mitigate the frequency deviation, otherwise, the PV system adaptively switches to the MPPT mode.

Conclusions
This study presents a unified strategy for active frequency regulation and maximum power point tracking (MPPT) pertaining to photovoltaic sources. Initially, a characteristic curve-fitting method is proposed, aimed at precisely and flexibly regulating the PV output power. The employment of a PV characteristic curve, mirroring the PV's real P-U curve, augments control performance and robustness. Based on the proposed active power control method, an innovative P-f type droop control is devised, enabling PV sources to deliver primary frequency regulation services. Under this control strategy, PV sources can adaptively participate in system frequency regulation or sustain the MPPT mode, based on the frequency and irradiation levels. Provided the irradiance is sufficient and the system frequency exceeds its nominal value, the PV source will relinquish some of its maximum available power proportionate to the frequency deviations. The frequency regulation capability attenuates with declining irradiance. Case study results substantiate that this strategy can curtail frequency deviation by approximately 25% during a 50% load shedding, with potential for further improvement via parameter settings.
The proposed control strategy is straightforward to implement as it does not require additional sensors, communication networks or complex calculations, thereby

Conclusions
This study presents a unified strategy for active frequency regulation and maximum power point tracking (MPPT) pertaining to photovoltaic sources. Initially, a characteristic curve-fitting method is proposed, aimed at precisely and flexibly regulating the PV output power. The employment of a PV characteristic curve, mirroring the PV's real P-U curve, augments control performance and robustness. Based on the proposed active power control method, an innovative P-f type droop control is devised, enabling PV sources to deliver primary frequency regulation services. Under this control strategy, PV sources can adaptively participate in system frequency regulation or sustain the MPPT mode, based on the frequency and irradiation levels. Provided the irradiance is sufficient and the system frequency exceeds its nominal value, the PV source will relinquish some of its maximum available power proportionate to the frequency deviations. The frequency regulation capability attenuates with declining irradiance. Case study results substantiate that this strategy can curtail frequency deviation by approximately 25% during a 50% load shedding, with potential for further improvement via parameter settings.
The proposed control strategy is straightforward to implement as it does not require additional sensors, communication networks or complex calculations, thereby

Conclusions
This study presents a unified strategy for active frequency regulation and maximum power point tracking (MPPT) pertaining to photovoltaic sources. Initially, a characteristic curve-fitting method is proposed, aimed at precisely and flexibly regulating the PV output power. The employment of a PV characteristic curve, mirroring the PV's real P-U curve, augments control performance and robustness. Based on the proposed active power control method, an innovative P-f type droop control is devised, enabling PV sources to deliver primary frequency regulation services. Under this control strategy, PV sources can adaptively participate in system frequency regulation or sustain the MPPT mode, based on the frequency and irradiation levels. Provided the irradiance is sufficient and the system frequency exceeds its nominal value, the PV source will relinquish some of its maximum available power proportionate to the frequency deviations. The frequency regulation capability attenuates with declining irradiance. Case study results substantiate that this strategy can curtail frequency deviation by approximately 25% during a 50% load shedding, with potential for further improvement via parameter settings.
The proposed control strategy is straightforward to implement as it does not require additional sensors, communication networks or complex calculations, thereby significantly reducing installation and maintenance costs and enhancing practical applicability. Nonetheless, it is imperative to consider the setting of the nominal point for PV sources to strike a balance between solar generation benefits and frequency regulation capability. Therefore, a trade-off exists between economic benefits and frequency regulation capability. Taking these considerations into account, the proposed strategy could be further refined in terms of optimal nominal point selection to boost the integration and utilization of PV sources, applicable to both distributed sources and large-scale PV plants.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.