Finite Control Set Model Predictive Control for an LCL-Filtered Grid-Tied Inverter with Full Status Estimations under Unbalanced Grid Voltage

: This paper proposes a novel ﬁnite control set model predictive control (FCS-MPC) strategy with merely grid-injected current sensors for an inductance-capacitance-inductance (LCL)-ﬁltered grid-tied inverter, which can obtain a sinusoidal grid-injected current whether three-phase grid voltages are balanced or not. Compared with the conventional FCS-MPC method, four compositions are added in the proposed FCS-MPC algorithm, where the grid voltage observer (GVO) and Luenberger observer are combined together to achieve full status estimations (including grid voltage, capacitor voltage, inverter-side current, and grid-injected current), while the sequence extractor and the reference generator are applied to eliminate the double frequency ripples of the active or reactive power, or the negative sequence component (NSC) of the grid-injected current caused by the unbalanced grid voltage. Simulation model and experimental platform are established to verify the e ﬀ ectiveness of the proposed FCS-MPC strategy, with full status estimations under both balanced and unbalanced


Introduction
Grid-tied inverters have been widely utilized in distributed generation systems, since they are the interfaces between DC sources and power grids [1,2]. In regards to the control of grid-tied inverters, besides the classical linear control schemes [3,4], a large number of nonlinear control strategies, such as model predictive control (including continuous control set model predictive control (CCS-MPC) and finite control set model predictive control (FCS-MPC)), sliding mode control, passivity-based control [5][6][7], and so on, were proposed. Among them, the FCS-MPC attracted significant attentions in recent years, owing to the technique advantages, including no need of the modulator, straightforward handling of nonlinearities and constraints, quick dynamic responses, and simple implementation [8][9][10][11].
Recently, Falkowski et al. [12] proposed an FCS-MPC method for the inductance-capacitanceinductance (LCL)-filtered grid-tied ac-dc converter, which can damp the oscillations caused by the filter resonance and acquire the high performance of grid currents. However, this FCS-MPC method requires to measure the inverter-side current, the capacitor voltage grid voltage, and the grid-injected current by sensors, increasing the cost and complexity. Aiming at reducing the number of sensors, many sensorless control schemes, such as virtual flux [13,14], state observer [15,16], and so on, were widely adopted.
Then, by taking the variables i1, i2, and uc as the system states, the dynamic model of inverter system can be further described by: Additionally, x = [i1αβ i2αβ ucαβ] T is the state vector and y = i2αβ is the output of the system, which denotes that only the grid-injected current is measured in this paper. The inverter output voltage viαβ is obtained by the combinations of switching signals Sa, Sb, Sc, which are described in Table 1, and it can be expressed as a state space complex-vector in αβ stationary reference frame:

Discrete-Time Model of LCL Filter
By neglecting the parasitic resistances, the dynamic model of inverter system can be expressed as follows: where the variables of current and voltage are expressed in the form of complex-vectors in the αβ stationary coordinate system (i.e., i 1αβ = i 1α + ji 1β , i 2αβ = i 2α + ji 2β , v iαβ = v iα + jv iβ , u cαβ = u cα + ju cβ , v gαβ = v gα + jv gβ ). Then, by taking the variables i 1 , i 2 , and u c as the system states, the dynamic model of inverter system can be further described by: where Additionally, x = [i 1αβ i 2αβ u cαβ ] T is the state vector and y = i 2αβ is the output of the system, which denotes that only the grid-injected current is measured in this paper. The inverter output voltage v iαβ Energies 2019, 12, 2691 4 of 22 is obtained by the combinations of switching signals S a , S b , S c , which are described in Table 1, and it can be expressed as a state space complex-vector in αβ stationary reference frame: v iαβ (n) = "0", n = {0 , 7} 2 3 U dc e j π 3 (n−1) , n = {1 , 2 . . . , 6} . (3) For the digital implementation of the control algorithm, the continuous-time dynamic model of the LCL filter can be represented in discrete time (the sampling time is T s ) by adopting the zero-order-hold (ZOH) discretization method: where matrices A 1 , B 1 , B 2 are

Conventional FCS-MPC Scheme
With the improved calculation power, the new generation of the digital signal processor (DSP), the field-programmable gate array (FPGA), and the dSPACE platform were utilized to implement the computationally-complex control algorithm for power electronics and drives. Consequently, the FCS-MPC method attracted a lot of attention in recent years, since the problem of computational burden can be preliminary solved. Unlike other linear or nonlinear control methods applied in the grid-tied inverter, the FCS-MPC algorithm has an obvious advantage of no need of a modulation stage. For the FCS-MPC scheme, by adopting a traversal method (each inverter output voltage vector is used to evaluate a cost function), the inverter output voltage vector with minimum cost function was selected. Then, the driving signal could be deduced according to Table 1.
For the LCL-filtered grid-tied inverter, the conventional FCS-MPC algorithm was implemented in the following steps [12]:

1.
Measure i 1 , i 2 , u c , and v g by using the current and voltage sensors; 2.
Calculate u * c and i * 1 according to the given reference of grid-injected current i * 2 ; 3.
Deduce the references of three state variables at next step x * (k + 1) by utilizing the Lagrangian Extrapolation method; 4.
Obtain the predictions of three state variables in the next sampling instant x (k + 1) for all possible inverter output voltage vectors based on the discrete-time model of the LCL filter expressed in Equation (4); 5.
Construct the cost function and define the weighting factor; 6.
Select the optimal inverter output voltage vector by minimizing the cost function;

Theoretical Analysis of the Proposed FCS-MPC Scheme with Full Status Estimations under Unbalanced Grid Voltage Condition
For the conventional FCS-MPC algorithm, it requires two kinds of current sensors to measure the inductor currents as well as two types of voltage sensors to probe the capacitor voltages and the grid voltages, increasing the cost and control complexity. In order to reduce the number of sensors and enhance the control reliability under unbalanced grid voltage condition, a novel FCS-MPC strategy is proposed in this paper. As depicted in Figure 1, compared with the conventional FCS-MPC scheme described in Section 2, four compositions including the GVO, the Luenberger observer, the sequence extractor, and the references generator are added in the proposed FCS-MPC strategy, which can achieve full status estimations and eliminate the negative effects caused by the unbalanced grid voltages. The detailed diagram of the proposed strategy is depicted in Figure 2.

Theoretical Analysis of the Proposed FCS-MPC Scheme with Full Status Estimations under Unbalanced Grid Voltage Condition
For the conventional FCS-MPC algorithm, it requires two kinds of current sensors to measure the inductor currents as well as two types of voltage sensors to probe the capacitor voltages and the grid voltages, increasing the cost and control complexity. In order to reduce the number of sensors and enhance the control reliability under unbalanced grid voltage condition, a novel FCS-MPC strategy is proposed in this paper. As depicted in Figure 1, compared with the conventional FCS-MPC scheme described in Section 2, four compositions including the GVO, the Luenberger observer, the sequence extractor, and the references generator are added in the proposed FCS-MPC strategy, which can achieve full status estimations and eliminate the negative effects caused by the unbalanced grid voltages. The detailed diagram of the proposed strategy is depicted in Figure 2.  Figure 2. The block diagram of the proposed strategy under the unbalanced grid voltage condition.

Grid Voltage Observer
As shown in Figure 2, the adaptive filter consists of the SOGI, the filter gain coefficient k, and the output feedback. It is vital to know the value of the grid angular frequency when the adaptive filter is utilized. In practice, the grid frequency may deviate from the rated value. When assuming that the grid frequency is a constant, the effect of the adaptive filter will drop sharply when the grid frequency shifts, making it difficult to obtain the correct information of the grid voltage. Hence, an approach for detecting the grid frequency is required.
In a three-phase system, the phase lock method of SRF-PLL is usually adopted. Under the ideal grid voltage condition, the SRF-PLL yields great performance and tracks the variable grid frequency accurately. However, under the unbalanced grid voltage condition, the overall dynamic performance of the SRF-PLL would become unacceptably deteriorated, due to the negative effects caused by the unbalanced grid voltages. For handling this problem, a PLL based on the decoupled double synchronous reference frame is utilized [22], which isolates the PSC and NSC of the grid voltages, and then takes the PSC as the input signal of PLL. Additionally, the output grid angular frequency is taken as the input angular frequency of the adaptive filter, making the system achieve the frequencyadaptive under the unbalanced grid voltage condition.
In order to analyze the performance of the adaptive filter on tracking the grid frequency, the transfer functions are expressed as follows: where ωp and k set the center angular frequency and the damping factor of the adaptive filter, respectively. u is the input signal of the adaptive filter, while û and ˆ⊥ u are the output signals. The amplitude and phase response of the adaptive filter can be calculated as:

Grid Voltage Observer
As shown in Figure 2, the adaptive filter consists of the SOGI, the filter gain coefficient k, and the output feedback. It is vital to know the value of the grid angular frequency when the adaptive filter is utilized. In practice, the grid frequency may deviate from the rated value. When assuming that the grid frequency is a constant, the effect of the adaptive filter will drop sharply when the grid frequency shifts, making it difficult to obtain the correct information of the grid voltage. Hence, an approach for detecting the grid frequency is required.
In a three-phase system, the phase lock method of SRF-PLL is usually adopted. Under the ideal grid voltage condition, the SRF-PLL yields great performance and tracks the variable grid frequency accurately. However, under the unbalanced grid voltage condition, the overall dynamic performance of the SRF-PLL would become unacceptably deteriorated, due to the negative effects caused by the unbalanced grid voltages. For handling this problem, a PLL based on the decoupled double synchronous reference frame is utilized [22], which isolates the PSC and NSC of the grid voltages, and then takes the PSC as the input signal of PLL. Additionally, the output grid angular frequency is taken as the input angular frequency of the adaptive filter, making the system achieve the frequency-adaptive under the unbalanced grid voltage condition.
In order to analyze the performance of the adaptive filter on tracking the grid frequency, the transfer functions are expressed as follows: where ω p and k set the center angular frequency and the damping factor of the adaptive filter, respectively. u is the input signal of the adaptive filter, whileû andû ⊥ are the output signals. The amplitude and phase response of the adaptive filter can be calculated as: It can be deduced from Equations (7) and (8), if the angular frequency ω of the input signal u is equal to the center angular frequency ω p in the steady state, G 1 ( jω) = G 2 ( jω) = 1, ∠G 1 ( jω) = 0, ∠G 2 ( jω) = −π/2. Consequently, we can get the conclusion that in combination with the SRF-PLL, although the grid frequency is variable, the adaptive filter can track the input signals accurately without any error in the steady state. And, the output signals are a pair of orthogonal quantities, whereû and u are in same phase, butû ⊥ is 90 • lag respect to the input signal u.
For the LCL-filtered grid-tied inverter system, due to the high impedance of the filtering capacitor at low frequency, the current of this capacitor can be neglected. When neglecting the current of the filtering capacitor, the continue-time model of grid-tied inverter with LCL filter in αβ coordinate system can be expressed as follows: Since the differential of the sinusoidal signal can be transformed into the in-phase or inverted value of its quadrature signal, based on the adaptive filter, Equation (9) can be written in the following form: v And then, Equations (10) and (11) can be simplified as follows: wherev gαβ ,v ⊥ gαβ represent the outputs of the GVO (i.e., the estimated grid voltage and its quadrature signal, respectively).v ioptαβ ,v ⊥ ioptαβ andî 2αβ ,î ⊥ 2αβ are the outputs of the adaptive filter, whose input signals are v iopt and i 2 , respectively.

Luenberger Observer
To further reduce the number of sensors, the Luenberger observer is adopted to combine with the GVO, where the inverter-side current sensors and the capacitor voltage sensors can be saved. The state-space model of the Luenberger observer in the discrete-time domain can be described as follows: x where L = [l 1 l 2 l 3 ] T represents the observer gain vector and C c = [0 1 0] is the output vector, which denotes that the grid-injected current is measured in this paper. The grid voltagev g is observed by using GVO. By defining the estimation error of ∆x(k) = x(k) −x(k), based on Equation (13), the dynamics of state observation error is derived as: Hence, if the matrix of A 1 -LC c is Hurwitz, the error vector will converge to zero. The observability , which indicates that the system is observable.
Therefore, the eigenvalues of the observer can be assigned arbitrarily, and the characteristic polynomial of the observer can be set as: where p 1 , p 2 , and p 3 are the desired poles of the Luenberger observer. In order to obtain the values of L, p 1 , p 2 , and p 3 need to be ensured. It is usually easier to identify the poles first in the s-domain and then map them to the z-domain via z = exp(sT s ). In the s-domain, the closed-loop characteristic polynomial can be described as (s + α od )(s 2 + 2ζ or ω or s + ω 2 or ). Then, the poles p 1 , p 2 , and p 3 in the z-domain can be expressed as follows: where the pair of complex-conjugate poles, determined by ζ or and ω or , are set to decide the dominant dynamics of the estimation errors. The real pole α od is located at a higher frequency. ζ or is the damping ratio, usually set as 0.707. The range of the natural frequency ω or with respect to the resonance frequency of LCL filter is regarded from 0.5 to 1, and the value of α od is five to 10 times larger than the pair of complex-conjugate poles. Hence, the gain vector L can be deduced by solving Equation (15). A common method to solve Equation (15) is utilizing the MATLAB function (i.e., acker). When the observer gain vector L is deduced, three state variables can be estimated to be applied to the conventional FCS-MPC algorithm. Hence, the inverter-side current sensors and capacitor voltage sensors can be saved. However, since the negative effects introduced by the NSC of the unbalanced grid voltage are not considered into the control algorithm, the performance will be deteriorated when the grid voltage falls into unbalance.

Sequence Extractor
To eliminate the adverse effects caused by the NSC of unbalanced grid voltages, the sequence extractor of the grid voltages is required.
Under the unbalanced grid voltage condition, by neglecting the zero sequence components, the positive-and negative-components can be extracted from the grid voltage. Consequently, the grid voltages, which are estimated by utilizing the GVO in this paper, can be described as the sum of PSC and NSC in αβ reference frame: wherev g is the estimated grid voltage and V g is its amplitude.v p gα ,v p gβ andv n gα ,v n gβ are the PSC and NSC of the estimated grid voltage in αβ reference frame, respectively. ϕ p and ϕ n are the initial phase angles, and ω is the grid angular frequency.
By utilizing the delayed signal cancellation (DSC) method [23], when the delay time is T g /4, Equation (17) can be rewritten as: Hence, according to Equations (17) and (18), the PSC and NSC of the estimated grid voltage can be deduced in the following equations: v ⊥ gα andv ⊥ gβ are 90 • lag in respect tov gα andv gβ , respectively, thus they are equivalent to a delay of T g /4. Then, the Equation (19) can be written as follows based on complex-vectors in the αβ stationary coordinate system: Hence, based on the GVO, the PSC and NSC of the grid voltage can be deduced, which can be utilized for generating the references of grid-injected currents in the next part.

Reference Generator
For the unbalanced grid voltages and currents without a zero sequence, the grid voltage and grid-injected current also can be expressed as the sum of the PSC and NSC in dq reference frame: wherev g is the estimated grid voltage obtained by the GVO.v p gdq andv n gdq are the PSC and NSC of the estimated grid voltages in the dq rotating reference system, respectively. Similarly, i p 2dq and i n 2dq are the PSC and NSC of the grid-injected currents in the dq rotating reference system, respectively. ω is the grid angular frequency.
Based on Equation (21), the power of the grid side can be expressed as follows: According to Equation (22), the active power of p and reactive power of q can be described as: where p 0 and q 0 are the average values of the active power and reactive power, respectively. p c2 , p s2 and q c2 , q s2 are ripples of p and q, respectively. As derived in [24], the average value and ripples of p and q can be further expressed as follows: The relationship between PSC and NSC of the estimated grid voltage or grid-injected current in the dq rotating coordinate system and αβ stationary coordinate system is where x represents the estimated grid voltages or grid-injected currents. By utilizing Equation (25), Equation (24) can be expressed in αβ stationary coordinate system as: where It can be seen from Equations (26) and (27) that if the grid voltages are decided, then the inverter has four controllable freedoms (i p 2α , i p 2β , i n 2α , i n 2β ). This implies that only four control targets can be established [25]. Normally, the average values of the active power p 0 and reactive power q 0 are controlled to track their references. Consequently, only the remaining two control freedoms can be selected. According to desired control targets, Equations (26) and (27) can be simplified and solved in the following three cases also introduced in [17]: Case 1: To eliminate active power ripples In this case, the control target is to remove the double frequency ripples of active power (i.e., p c2 and p s2 in Equation (26) are set to zero (Hence, k 1 , k 2 are zero. q c2 and q s2 are uncontrolled and thus there are reactive power ripples under unbalanced grid voltage. The control objective is equivalent to handle the following equation: Based on Equations (26), (27), and (28), the reference of the grid-injected current can be calculated as follows: Case 2: To eliminate reactive power ripples In this case, the control target is to eliminate the reactive power ripples (i.e., q c2 and q s2 in Equation (26) are considered as zero). Hence, k 3 , k 4 are both zero. p c2 and p s2 are uncontrolled and thus active power ripples exist under unbalanced grid voltage conditions. The control objective is equivalent to solve the following equation: Similarly, based on Equations (26), (27), and (30), the reference of the grid-injected current can be obtained as: Case 3: To achieve balanced and sinusoidal grid-injected currents In this case, the control target is to achieve balanced and sinusoidal grid-injected currents, hence, i n 2α and i n 2β are regarded as zero. i p 2α and i p 2β should be used to satisfy the equations of p 0 = P * and q 0 = Q * . Since p c2 , p s2 and q c2 , q s2 are uncontrolled and thus both active power and reactive power ripples exist if grid voltages are unbalanced. The control objective is equivalent to deal with the following equation: According to Equations (26), (27), and (32), the reference of the grid-injected current can be deduced as: . A detailed comparison of the proposed FCS-MPC scheme with three different control targets under unbalanced grid voltage is presented in Table 2.

Detailed Implementation of the Proposed FCS-MPC Strategy
According to the analysis mentioned above, the reference of the grid-injected current can be deduced. Then, it is taken into the FCS-MPC algorithm to alleviate the adverse effects caused by the unbalanced grid voltage. Meanwhile, based on the GVO and Luenberger observer, at least six sensors are saved. Besides, since the practical implementation of FCS-MPC requires considering the negative effect of the computational delay, the discrete-time model of the LCL filter are shifted one step forward in order to eliminate this time delay (i.e., the values at the (k + 2) th instant rather than the (k + 1) th instant should be applied in the cost function) [26]. The flowchart of the proposed FCS-MPC algorithm is depicted in Figure 3, where the detailed implementation steps can be summarized as follows: 1.
Measure the grid-injected currents of phase a and b, and deduce the c-phase grid-injected current by 2.
Estimate the inverter-side inductor currentî 1 , the grid-injected currentî 2 , and the capacitor voltageû c by using the Luenberger observer, which is based on the observed grid voltagesv g and the measured grid-injected currents; 4.
Deduce the references of the inverter-side current i * 1αβ and the capacitor voltage u * cαβ in αβ stationary coordinate system based on i * 2αβ andv gαβ . 6.
Calculate the references of the three estimated state variables x * at the (k + 2) th instant by utilizing the Lagrangian extrapolation, which can be expressed as (x = [i 1 i 2 u c ] T ): 7.
Predict the three estimated state variables x at the (k + 1) th instant according to Equation (4), and calculate the estimated grid voltagev g at the (k + 1) th instant bŷ v g (k + 1) =v g (k)e jωT s . 8.
Predict the three estimated state variables x at the (k + 2) th instant, which is based on the discrete-time model of the LCL filter described in Equation (4); 9.
Construct the following cost function and then each voltage vector described in Table 1 is taken into the cost function. Consequently, seven different value of cost functions J can be obtained. where ε are the errors of controlled variables between reference value and estimated value. The weighting factor λ uc is set to achieve active damping [27], and λ i2 is tuned to obtain a high quality grid-injected current [12], and the weighting factors are designed according to [28]. 10. The inverter voltage vector with minimum cost function v iopt is selected. Correspondingly, the driving signals can be deduced according to Table 1.
uc ucα ucβ J = ε k + +ε k + +λ ε k + +ε k + +λ ε k + +ε k + (38) where ε are the errors of controlled variables between reference value and estimated value. The weighting factor λuc is set to achieve active damping [27], and λi2 is tuned to obtain a high quality grid-injected current [12], and the weighting factors are designed according to [28]. 10. The inverter voltage vector with minimum cost function viopt is selected. Correspondingly, the driving signals can be deduced according to Table 1.

Simulation Results
In order to verify the effectiveness of the proposed FCS-MPC strategy, the simulation test is carried out in MATLAB/Simulink. The system variables and parameters are listed as follows. The dcside voltage is 150 V, the three-phase grid voltage is 50 V (RMS) and the rated power is 750 W. When the three-phase grid voltage falls into unbalance, b-phase voltage dips to 20 V (RMS) in this paper. The sampling frequency/sampling time is 25 kHz/40 μs. The inverter-side inductance is 2.4 mH, the grid-side inductance is 1.2 mH, and the filtering capacitor is 6 μF. The simulation results are shown in Figures 4-8, respectively. Figure 4 shows the simulation results of the measured grid voltages, the estimated grid voltages, and their errors in the αβ stationary reference system, where the estimations are obtained by GVO. Figure 5 displays the simulation results of the three measured state variables and their estimations in b-phase, as well as the errors between measurements and estimations, where these estimations are calculated by the Luenberger observer. Additionally, the performance of the GVO and Luenberger observer is verified in cases 1, 2, and 3, respectively. It is worth noting that for all measured signals in this paper, including the grid voltage, the grid-injected current, the capacitor voltage, and the inverter-side current, only the grid-injected currents are applied for the GVO and Luenberger observer estimations. The other measured signals are only used for comparisons with the observation results.
It can be seen from Figure 4 that the observations of the grid voltage in α-axis and β-axis are both consistent with the measurements, and the errors of them under balanced and unbalanced grid voltages are both relatively small, which can be acceptable. Hence, the design of GVO is successful and effective. Since the reference of the grid-injected currents is variable for the different control

Simulation Results
In order to verify the effectiveness of the proposed FCS-MPC strategy, the simulation test is carried out in MATLAB/Simulink. The system variables and parameters are listed as follows. The dc-side voltage is 150 V, the three-phase grid voltage is 50 V (RMS) and the rated power is 750 W. When the three-phase grid voltage falls into unbalance, b-phase voltage dips to 20 V (RMS) in this paper. The sampling frequency/sampling time is 25 kHz/40 µs. The inverter-side inductance is 2.4 mH, the grid-side inductance is 1.2 mH, and the filtering capacitor is 6 µF. The simulation results are shown in Figures 4-8, respectively. Figure 4 shows the simulation results of the measured grid voltages, the estimated grid voltages, and their errors in the αβ stationary reference system, where the estimations are obtained by GVO. Figure 5 displays the simulation results of the three measured state variables and their estimations in b-phase, as well as the errors between measurements and estimations, where these estimations are calculated by the Luenberger observer. Additionally, the performance of the GVO and Luenberger observer is verified in cases 1, 2, and 3, respectively. It is worth noting that for all measured signals in this paper, including the grid voltage, the grid-injected current, the capacitor voltage, and the inverter-side current, only the grid-injected currents are applied for the GVO and Luenberger observer estimations. The other measured signals are only used for comparisons with the observation results.
It can be seen from Figure 4 that the observations of the grid voltage in α-axis and β-axis are both consistent with the measurements, and the errors of them under balanced and unbalanced grid voltages are both relatively small, which can be acceptable. Hence, the design of GVO is successful and effective. Since the reference of the grid-injected currents is variable for the different control targets under the unbalanced grid voltage, we can find that the amplitudes of the inverter-side currents and grid-injected currents vary with the change of three cases in Figure 5, where the dashed line expresses the amplitude of these two currents under the balanced grid voltages, and, the estimations of these two currents can track the measurements accurately. Additionally, the observed capacitor voltage also can follow the measured signal exactly, even under the unbalanced grid voltage condition. Consequently, the conclusion we can draw is that the estimations of the inverter-side currents, grid-injected currents, capacitor voltages, and grid voltages can replace the measurements and at least six sensors can be saved. Additionally, the effectiveness of the proposed FCS-MPC strategy for decreasing the number of sensors is preliminarily verified. targets under the unbalanced grid voltage, we can find that the amplitudes of the inverter-side currents and grid-injected currents vary with the change of three cases in Figure 5, where the dashed line expresses the amplitude of these two currents under the balanced grid voltages, and, the estimations of these two currents can track the measurements accurately. Additionally, the observed capacitor voltage also can follow the measured signal exactly, even under the unbalanced grid voltage condition. Consequently, the conclusion we can draw is that the estimations of the inverterside currents, grid-injected currents, capacitor voltages, and grid voltages can replace the measurements and at least six sensors can be saved. Additionally, the effectiveness of the proposed FCS-MPC strategy for decreasing the number of sensors is preliminarily verified.  To further verify the effectiveness of the proposed FCS-MPC scheme, the performance comparison results between the proposed control method and the FCS-MPC with full status measurements are displayed in Figure 6. It can be seen that. Compared with the FCS-MPC with full status measurements, the quality of the grid-injected current of the proposed control method declines a little; however, it still can meet the harmonic standard of IEEE (the THD of the grid-injected current  targets under the unbalanced grid voltage, we can find that the amplitudes of the inverter-side currents and grid-injected currents vary with the change of three cases in Figure 5, where the dashed line expresses the amplitude of these two currents under the balanced grid voltages, and, the estimations of these two currents can track the measurements accurately. Additionally, the observed capacitor voltage also can follow the measured signal exactly, even under the unbalanced grid voltage condition. Consequently, the conclusion we can draw is that the estimations of the inverterside currents, grid-injected currents, capacitor voltages, and grid voltages can replace the measurements and at least six sensors can be saved. Additionally, the effectiveness of the proposed FCS-MPC strategy for decreasing the number of sensors is preliminarily verified.  To further verify the effectiveness of the proposed FCS-MPC scheme, the performance comparison results between the proposed control method and the FCS-MPC with full status measurements are displayed in Figure 6. It can be seen that. Compared with the FCS-MPC with full status measurements, the quality of the grid-injected current of the proposed control method declines a little; however, it still can meet the harmonic standard of IEEE (the THD of the grid-injected current does not exceed 5%). Note that, compared with the FCS-MPC with full status measurements, at least six sensors can be saved for the proposed method.   Figure 7 depicts that the grid voltages vary from balanced to unbalanced at 0.25 s, and then back to balanced at 0.55 s. In order to achieve grid synchronization and ensure the maximum energy injected into the grid, Q* is set to zero in this paper. Hence, based on Equations (29), (31), and (33), the references of the grid-injected currents can be calculated. The bottom of Figure 7 reveals the performance of the grid-injected currents under the unbalanced grid voltages by utilizing the proposed FCS-MPC strategy based on three different control targets. Figure 8 represents the simulation results of the active power and reactive power under the balanced and unbalanced grid voltages. It can be seen from Figure 7 that the three-phase grid-injected currents are sinusoidal but unbalanced in case 1 and case 2; however, the grid-injected currents are balanced and sinusoidal in case 3. For Figure 8, we found that the double frequency ripples of active power are eliminated in case 1, the double frequency ripples of reactive power are eliminated in case 2, and the double frequency ripples of both the active and reactive power exist in case 3, which is in agreement with the theoretical analysis described in Table 2. Therefore, according to it can be verified that the proposed FCS-MPC strategy with merely grid-injected current sensors is effective, and it can obtain good performance of the grid-injected currents even under the unbalanced grid voltage condition.   Figure 7 shows the simulation results of the grid voltages and grid-injected currents under balanced and unbalanced grid voltages. The top of Figure 7 depicts that the grid voltages vary from balanced to unbalanced at 0.25 s, and then back to balanced at 0.55 s. In order to achieve grid synchronization and ensure the maximum energy injected into the grid, Q* is set to zero in this paper. Hence, based on Equations (29), (31), and (33), the references of the grid-injected currents can be calculated. The bottom of Figure 7 reveals the performance of the grid-injected currents under the unbalanced grid voltages by utilizing the proposed FCS-MPC strategy based on three different control targets. Figure 8 represents the simulation results of the active power and reactive power under the balanced and unbalanced grid voltages. It can be seen from Figure 7 that the three-phase grid-injected currents are sinusoidal but unbalanced in case 1 and case 2; however, the grid-injected currents are balanced and sinusoidal in case 3. For Figure 8, we found that the double frequency ripples of active power are eliminated in case 1, the double frequency ripples of reactive power are eliminated in case 2, and the double frequency ripples of both the active and reactive power exist in case 3, which is in agreement with the theoretical analysis described in Table 2. Therefore, according to Figures 4-8, it can be verified that the proposed FCS-MPC strategy with merely grid-injected current sensors is effective, and it can obtain good performance of the grid-injected currents even under the unbalanced grid voltage condition.

Experimental Results
To further verify the effectiveness of the proposed FCS-MPC strategy, the laboratory test-rig depicted in Figure 9 is established. The experimental parameters were the same as the simulation ones. The programmable ac source (Chroma 61830) is utilized to simulate the balanced and unbalanced grid voltage condition. The power stage consists of a two-level voltage-source inverter (Danfoss-FC320) with a dc-link voltage provided by Chroma 62150 H-600S DC power supply. The digital control algorithm is implemented in dSPACE 1202 platform, where a control desk project is established to regulate control parameters and reference values, as well as display the experimental To further verify the effectiveness of the proposed FCS-MPC scheme, the performance comparison results between the proposed control method and the FCS-MPC with full status measurements are displayed in Figure 6. It can be seen that. Compared with the FCS-MPC with full status measurements, the quality of the grid-injected current of the proposed control method declines a little; however, it still can meet the harmonic standard of IEEE (the THD of the grid-injected current does not exceed 5%). Note that, compared with the FCS-MPC with full status measurements, at least six sensors can be saved for the proposed method. Figure 7 shows the simulation results of the grid voltages and grid-injected currents under balanced and unbalanced grid voltages. The top of Figure 7 depicts that the grid voltages vary from balanced to unbalanced at 0.25 s, and then back to balanced at 0.55 s. In order to achieve grid synchronization and ensure the maximum energy injected into the grid, Q* is set to zero in this paper. Hence, based on Equations (29), (31), and (33), the references of the grid-injected currents can be calculated. The bottom of Figure 7 reveals the performance of the grid-injected currents under the unbalanced grid voltages by utilizing the proposed FCS-MPC strategy based on three different control targets. Figure 8 represents the simulation results of the active power and reactive power under the balanced and unbalanced grid voltages. It can be seen from Figure 7 that the three-phase grid-injected currents are sinusoidal but unbalanced in case 1 and case 2; however, the grid-injected currents are balanced and sinusoidal in case 3. For Figure 8, we found that the double frequency ripples of active power are eliminated in case 1, the double frequency ripples of reactive power are eliminated in case 2, and the double frequency ripples of both the active and reactive power exist in case 3, which is in agreement with the theoretical analysis described in Table 2. Therefore, according to Figures 4-8, it can be verified that the proposed FCS-MPC strategy with merely grid-injected current sensors is effective, and it can obtain good performance of the grid-injected currents even under the unbalanced grid voltage condition.

Experimental Results
To further verify the effectiveness of the proposed FCS-MPC strategy, the laboratory test-rig depicted in Figure 9 is established. The experimental parameters were the same as the simulation ones. The programmable ac source (Chroma 61830) is utilized to simulate the balanced and unbalanced grid voltage condition. The power stage consists of a two-level voltage-source inverter (Danfoss-FC320) with a dc-link voltage provided by Chroma 62150 H-600S DC power supply. The digital control algorithm is implemented in dSPACE 1202 platform, where a control desk project is established to regulate control parameters and reference values, as well as display the experimental results which cannot be probed by the Yokogawa DL 1640 digital oscilloscope assembled with two current probes HIOKI 3276, including the inverter-side currents, grid-injected currents, capacitor voltages, grid voltages, active power, and reactive power in this paper.
(Danfoss-FC320) with a dc-link voltage provided by Chroma 62150 H-600S DC power supply. The digital control algorithm is implemented in dSPACE 1202 platform, where a control desk project is established to regulate control parameters and reference values, as well as display the experimental results which cannot be probed by the Yokogawa DL 1640 digital oscilloscope assembled with two current probes HIOKI 3276, including the inverter-side currents, grid-injected currents, capacitor voltages, grid voltages, active power, and reactive power in this paper.  Figure 9. The experimental test setup. Figure 10 displays the experimental comparisons between measurements and estimations for the grid voltage in α-axis and β-axis, respectively, when grid voltage varies from unbalanced to balanced based on case 1. Figure 11 shows the experimental comparisons between measurements and estimations for the state variables in b-phase when grid voltage varies from unbalanced to balanced based on case 1. In Figures 10 and 11, the red line expresses the measurement, the green line  Figure 10 displays the experimental comparisons between measurements and estimations for the grid voltage in α-axis and β-axis, respectively, when grid voltage varies from unbalanced to balanced based on case 1. Figure 11 shows the experimental comparisons between measurements and estimations for the state variables in b-phase when grid voltage varies from unbalanced to balanced based on case 1. In Figures 10 and 11, the red line expresses the measurement, the green line denotes the estimation, and the blue line represents their error. Additionally, it can be found that the observations track the measurements accurately both in steady and dynamic states, reflecting the effectiveness of the GVO and Luenberger observer. Therefore, the sensors of inverter-side current, capacitor voltage, and grid voltage can be saved by taking the observations instead of measurements into the FCS-MPC algorithm. The experimental comparisons based on case 2 and case 3 also can be obtained, and the effects of track are similar to the case 1, thus it is not shown in this paper. It should be noticed that, for all measurements, except for the grid-injected currents, the other measured signals are merely utilized for comparisons with the estimated results.    Figure 12 depicts the experimental performance comparison results between the proposed control method and the FCS-MPC with full status measurements. Although the THD of the grid-injected current by adopting the proposed control method is higher than the one using the FCS-MPC with full status measurements, the difference between these two THD values of grid-injected current is not so large, where they both meet the harmonic standard of IEEE (the THD of grid-injected current does not exceed 5%). Therefore, the effectiveness of the proposed control method is verified. injected current by adopting the proposed control method is higher than the one using the FCS-MPC with full status measurements, the difference between these two THD values of grid-injected current is not so large, where they both meet the harmonic standard of IEEE (the THD of grid-injected current does not exceed 5%). Therefore, the effectiveness of the proposed control method is verified.   Figures 13-15 show the experimental results of grid-injected currents, active and reactive power by utilizing the proposed FCS-MPC scheme based on three different control targets, which are adopted to eliminate the negative effects caused by the unbalanced grid voltage, respectively. Due to the limitation of the experimental device, only a-phase and b-phase voltages and currents are measured. It can be observed from Figure 13 that the grid-injected currents are sinusoidal but unbalanced (the current in b-phase is largest, while the amplitude of currents in a-phase and c-phase are same) in the steady states, and the dynamic performance is good when the grid voltages vary  eliminate the negative effects caused by the unbalanced grid voltage, respectively. Due to the limitation of the experimental device, only a-phase and b-phase voltages and currents are measured. It can be observed from Figure 13 that the grid-injected currents are sinusoidal but unbalanced (the current in b-phase is largest, while the amplitude of currents in a-phase and c-phase are same) in the steady states, and the dynamic performance is good when the grid voltages vary from balanced to unbalanced and then back to balanced. Note that further optimization in the research is possible by increasing the sampling frequency, but it is not the subject of this paper, and the double frequency ripples of the active power are eliminated, while the ripples of the reactive power still exist under the unbalanced grid voltage. Hence, the experimental results are in agreement with the theoretical analysis of case 1.
(a) (b) Figure 12. The experimental performance comparison results between (a) the proposed control method and (b) the FCS-MPC with full status measurements. Figures 13-15 show the experimental results of grid-injected currents, active and reactive power by utilizing the proposed FCS-MPC scheme based on three different control targets, which are adopted to eliminate the negative effects caused by the unbalanced grid voltage, respectively. Due to the limitation of the experimental device, only a-phase and b-phase voltages and currents are measured. It can be observed from Figure 13 that the grid-injected currents are sinusoidal but unbalanced (the current in b-phase is largest, while the amplitude of currents in a-phase and c-phase are same) in the steady states, and the dynamic performance is good when the grid voltages vary from balanced to unbalanced and then back to balanced. Note that further optimization in the research is possible by increasing the sampling frequency, but it is not the subject of this paper, and the double frequency ripples of the active power are eliminated, while the ripples of the reactive power still exist under the unbalanced grid voltage. Hence, the experimental results are in agreement with the theoretical analysis of case 1.  As demonstrated in Figure 14, the performance of the grid-injected currents are acceptable both in steady (the current in b-phase is smallest, while the amplitude of currents in a-phase and c-phase are same) and dynamic states, and the double frequency ripples of the reactive power are eliminated, while the active power ripples still exist under the unbalanced grid voltage. Consequently, the experimental results verify the correctness of the theory of case 2.
Balanced Unbalanced Figure 13. Experimental results of (a) grid-injected current when grid voltage varies from balanced to unbalanced, (b) grid-injected current when grid voltage varies from unbalanced to balanced, and (c) active and reactive power when grid voltage varies from unbalanced to balanced in case 1.

Unbalanced Balanced
(c) Figure 14. Experimental results of (a) grid-injected current when grid voltage varies from balanced to unbalanced, (b) grid-injected current when grid voltage varies from unbalanced to balanced, and (c) active and reactive power when grid voltage varies from unbalanced to balanced in case 2. Figure 15 displays the sinusoidal grid-injected currents with different amplitudes under the balanced grid voltage and the unbalanced grid voltage. However, it can be found that the currents are balanced whether three-phase grid voltages are balanced or not, and the transient performance is good when the grid voltages change from balanced to unbalanced and then back to balanced. Furthermore, the double frequency ripples of the active and reactive power are not removed under the unbalanced grid voltages. These performances are also consistent with the theoretical analysis of case 3.
[70V/div] g v Balanced Unbalanced As demonstrated in Figure 14, the performance of the grid-injected currents are acceptable both in steady (the current in b-phase is smallest, while the amplitude of currents in a-phase and c-phase are same) and dynamic states, and the double frequency ripples of the reactive power are eliminated, while the active power ripples still exist under the unbalanced grid voltage. Consequently, the experimental results verify the correctness of the theory of case 2. Figure 15 displays the sinusoidal grid-injected currents with different amplitudes under the balanced grid voltage and the unbalanced grid voltage. However, it can be found that the currents are balanced whether three-phase grid voltages are balanced or not, and the transient performance is good when the grid voltages change from balanced to unbalanced and then back to balanced. Furthermore, unbalanced, (b) grid-injected current when grid voltage varies from unbalanced to balanced, and (c) active and reactive power when grid voltage varies from unbalanced to balanced in case 2. Figure 15 displays the sinusoidal grid-injected currents with different amplitudes under the balanced grid voltage and the unbalanced grid voltage. However, it can be found that the currents are balanced whether three-phase grid voltages are balanced or not, and the transient performance is good when the grid voltages change from balanced to unbalanced and then back to balanced. Furthermore, the double frequency ripples of the active and reactive power are not removed under the unbalanced grid voltages. These performances are also consistent with the theoretical analysis of case 3. Unbalanced Balanced (c) Figure 15. Experimental results of (a) grid-injected current when grid voltage varies from balanced to unbalanced, (b) grid-injected current when grid voltage varies from unbalanced to balanced, and (c) active and reactive power when grid voltage varies from unbalanced to balanced in case 3.
According to the simulation and experimental results shown above, it can be found that the proposed FCS-MPC strategy with merely grid-injected current sensors can work well whether the grid voltages are balanced or not.

Conclusions
In this paper, a novel FCS-MPC strategy with full status observations based on GVO and Luenberger observer is proposed for the LCL-filtered grid-tied inverter under the unbalanced grid voltage condition. A step-by-step design procedure of the proposed control scheme is described in detail. The simulation model and experimental platform are established to verify the performance of proposed control method, where the three different reference cases of the grid-injected currents in the proposed cost function are analyzed. By theoretical analysis, simulation, and experimental According to the simulation and experimental results shown above, it can be found that the proposed FCS-MPC strategy with merely grid-injected current sensors can work well whether the grid voltages are balanced or not.

Conclusions
In this paper, a novel FCS-MPC strategy with full status observations based on GVO and Luenberger observer is proposed for the LCL-filtered grid-tied inverter under the unbalanced grid voltage condition. A step-by-step design procedure of the proposed control scheme is described in detail. The simulation model and experimental platform are established to verify the performance of proposed control method, where the three different reference cases of the grid-injected currents in the proposed cost function are analyzed. By theoretical analysis, simulation, and experimental verification, the following conclusions can be drawn:

1.
By utilizing the proposed control algorithm, the errors between observations and measurements are small enough under both balanced and unbalanced grid voltages, and the observations have a good dynamic response. Therefore, these observations can replace the measurements used in the control.

2.
Under the unbalanced grid voltage condition, the proposed FCS-MPC strategy obtains the satisfactory performance in both the steady and dynamic states for three different control targets. Note that only the grid-injected currents are measured.