A Fractional-Order Multi-Rate Repetitive Controller for Single-Phase Grid-Connected Inverters

Abstract

The multi-rate repetitive controller (MRC) can achieve zero steady-state error in tracking the reference current signal of grid-connected inverters, save the settling time effectively, and improve the running speed. However, when the grid frequency fluctuates, the harmonic suppression performance of MRC will degrade. Aiming at the problem of harmonic suppression performance degradation, a fractional-order MRC (FOMRC) based on the farrow structure fractional delay (FD) filter is proposed. Firstly, the equivalent digital model of MRC is established, and a Farrow structure fractional delay (FD) filter based on Taylor series expansion is selected as the internal model filter of MRC. The stability analysis and harmonic suppression characteristics of the FOMRC are analyzed. Then, the parameter design of FOMRC applied to an LCL single-phase grid-connected inverter control system is given. Finally, the simulation results show that the proposed method has better transient and steady-state performance than the CRC when the grid frequency fluctuates.

1. Introduction

With the proposal of a carbon peak and carbon neutral target in our country, the step of clean and low-carbon transition of the power system will be further accelerated. According to various policy orientations and the future direction of clean energy development, distributed energy will have a gradual growth and may become a major energy source in the future. Distributed energy is a kind of clean energy with strong adaptability to the environment. However, many non-linear power electronic switching devices are used in distributed generation, which will inject a large amount of harmonic current into the grid system and reduce the power quality [1]. The grid-connected inverter is the interface unit of power transmission from distributed generation system to the power grid, and its performance directly determines the quality of grid-connected current. According to the IEEE standard, the total harmonic distortion (THD) of grid-connected current must be less than 5%.

In order to reduce the harmonic content of the grid current of the inverter, L-type filter, and LCL-type filter are usually used to suppress the high-frequency harmonic at the switching frequency of the grid current. Under the same harmonic suppression requirements, the LCL filter is widely used because of its smaller actual volume and lower cost, which can better suppress grid current harmonics [2,3]. For the suppression of low-frequency harmonics in grid current, many controllers have been proposed. Among them, the commonly used include proportional integral (PI) controller [4], proportional resonant (PR) controller [5], proportional multi-resonant (PMR) controller, repetitive controller (RC), etc.

The repetitive controller (RC) [6,7,8,9] is widely used because it can realize the static error-free tracking and disturbance rejection of fundamental and harmonics of sinusoidal signals. However, when the conventional repetitive controller (CRC) is implemented digitally, the sampling rate of the repetitive controller and the switching frequency of the inverter power device are often 10 kHz or even higher, which makes the memory consumption of the controller too large and the calculation burden too heavy. For this reason, multiple rate repetitive controller (MRC) [10,11] has been proposed, in which the repetitive controller is set in a low sampling frequency environment, and the feedback part of the grid-connected current still adopts a high sampling rate. Compared with CRC, the memory consumption and calculation times of the MRC controller in each sampling period are greatly reduced [12]. Previous work has shown that MRC can maintain the same convergence rate and total harmonic distortion (THD) of grid-connected current as CRC, and CRC is a special case when MRC sampling factor is 1.

However, in the practical application of inverters, there is usually a small range fluctuation of the power grid frequency, which makes the internal mode of the RC periodic signal not match with the periodic signal to be tracked or eliminated, resulting in the attenuation of the resonant frequency gain [13,14]. The internal mode delay number $N=f_s/f_g$, that is, the ratio of the sampling frequency to the fundamental frequency. When grid frequency fluctuates, N and the reduced-sampled N may be a fraction, and fractional delay can not be achieved directly in digital systems, it will reduce the repetitive control gain and the harmonic suppression ability of the system. The variable sampling frequency method [15] can ensure that N is an integer to achieve appropriate harmonic suppression. However, the variable sampling rate means that the change of system dynamics, especially the change of the controlled object model, increases the difficulty of system stability analysis. The finite impulse response (FIR) filter based on the Lagrange interpolation method is widely used as the internal model filter of various repetitive controllers (odd RC [16], dual-mode RC [17], selective harmonic RC [18]) because of its simple structure and linear phase in the whole frequency band. The fractional delay link (FD) is approximated to obtain any numerical delay N, which improves the anti-frequency fluctuation ability of repetitive control. However, the FIR filter based on Lagrange interpolation needs to recalculate the coefficients of $M+1$ (M is the order of the FIR filter) sub-filters when the coefficients are updated, which greatly increases the amount of calculation and is slightly contrary to the original intention of using MRC in this paper.

Therefore, this paper proposes a fractional-order multi-rate repetitive controller (FOMRC) composed of a Farrow structure FD filter based on Taylor series expansion. FOMRC consists of an MRC with a low sampling rate and a relatively simple Farrow structure FD filter. Each sub-filter of the FD filter with a Farrow structure is designed offline in advance. For the new FD, only the fractional value of the fractional N is adjusted, which greatly reduces the computational burden of the controller. At present, the FOMRC remains to be further studied. This composite controller can not only reduce the computational burden of the digital system and retain the harmonic suppression ability of the CRC but can also accurately approximate the FD generated when the fundamental frequency of the power grid fluctuates. Improve the steady-state accuracy of multi-rate repetitive control.

2. Fractional Order Multi-Rate Repetitive Controller

2.1. MRC

Figure 1 shows the block diagram of a conventional plug-in MRC system. In the figure, the sampling period of the grid current feedback system is $T_s$, $E\left(z\right)$ is the error between the input signal $R\left(z\right)$ and the output signal $Y\left(z\right)$. In order to prevent the distortion of $E\left(z\right)$ signal outside the Nyquist frequency when entering the repetitive controller with a low sampling rate, the anti-spectral aliasing filter $F_1\left(z\right)$ is often used to intercept the signal in the $\left|\omega \right|>\pi /T_m$ part, and the signal $E\left(z_m\right)$ is obtained by downsampling as the input of RC. The output $U_r\left(z_m\right)$ of RC is up-sampled by zero-order holder (ZOH) to obtain $U_r\left(z\right)$, and then an anti-imaging filter $F_2\left(z\right)$ with cutoff frequency less than $\pi /T_m$ is connected to prevent distortion. $G_{pi}\left(z\right)$ is a proportional integral controller for stabilizing the closed-loop control system. $P\left(z\right)$ is the controlled object, and $D\left(z\right)$ is the disturbance source.

The sampling period of RC is $T_m$, which is m times of $T_s$, and m is the sampling factor. The relationship between the two sampling rates is

$$T_{\mathrm{m}}=m{T}_s;z=e^{s{T}_s};z_m=e^{s{T}_m}=z^m$$

(1)

In order to simplify the analysis, the multi-rate system is transformed into an equivalent single-rate system with a low sampling rate, and the structure is shown in Figure 2. $G_{pi}\left(z_m\right)$ and $P\left(z_m\right)$ are the equivalent transfer functions of PI controller and controlled object at low rates, respectively; $k_r$ is the internal mode gain; $Q\left(z_m\right)$ is used to improve the robustness of RC; $N_m=T_g/T_m$, $T_g$ is the fundamental frequency period of the reference signal and disturbance signal. $S\left(z_m\right)$ is a low-pass filter used to accelerate the high-frequency amplitude attenuation of the controlled object; $z_m^k$ is the phase-lead compensation, which is used to compensate the phase lag of the control system. The influence of $F_1\left(z\right)$, downsampling link, zero order holder and $F_2\left(z\right)$ on the system can also be compensated by the $z_m^k$.

The transfer function of MRC is

$$G_{mrc}=k_r\frac{z_m^{-N_m}Q\left(z_m\right)}{1-z_m^{-N_m}Q\left(z_m\right)}S\left(z_m\right)z_m^k$$

(2)

2.2. Farrow Structure FD Filter Design Based on Taylor Series

When $N_m$ is a fraction, $N_m$ can be divided into integer D and fractional d, $z_m^{-N_m}=z_m^{-D}{z}_m^{-d}$. The FD transfer function $z_m^{-d}$ can be expanded to a polynomial about d by the Taylor series [19]:

$$a_n\left(d\right)=z_m^{-d}=e^{-j\omega md}=\sum _{k=0}^M\frac{{(-j\omega m)}^k}{k!}{\left(d\right)}^k$$

(3)

where, $0\le d<1$.

The FD filter can be expressed as

$$G_d\left(z_m\right)=\sum _{n=0}^S{a}_n\left(d\right)z_m^{-n}$$

(4)

where S is the order of the FD filter and $a_n\left(d\right)$ is the polynomial of d.

Substituting (3) into (4), the Farrow [20] FD filter $G_d\left(z_m\right)$ can be obtained

$$G_d\left(z_m\right)=\sum _{k=0}^M\sum _{n=0}^S{a}_{nk}{z}_m^{-n}{d}^k=\sum _{k=0}^M{L}_k\left(z_m\right)d^k$$

(5)

$L_k\left(z_m\right)$ (k = 0,1,..., M) is the k subfilter in $G_d\left(z_m\right)$.

Usually choosing $S=M$, subfilter calculation formula based on Lagrange interpolation method is as follows [21,22]:

$$U=\left[\begin{array}{ccccc}{0}^0& 0^1& 0^2& \cdots & 0^S\\ 1^0& 1^1& 1^2& \cdots & 1^S\\ 2^0& 2^1& 2^2& \cdots & 2^S\\ ⋮& ⋮& ⋮& ⋮& ⋮\\ M^0& M^1& M^2& \cdots & M^S\end{array}\right]$$

(6)

$$z_{sub}={\left[\begin{array}{ccccc}1\hfill & z_m^{-1}\hfill & z_m^{-2}\hfill & \cdots \hfill & z_m^{-M}\hfill \end{array}\right]}^T$$

(7)

$$\begin{array}{ccc}\hfill f_{sub}& =U^{-1}{z}_{sub}\hfill & \hfill ={\left[\begin{array}{ccccc}{L}_0\left(z_m\right)\hfill & L_1\left(z_m\right)\hfill & L_2\left(z_m\right)\hfill & \cdots \hfill & L_M\left(z_m\right)\hfill \end{array}\right]}^T\end{array}$$

(8)

In the formula, $f_{sub}$ is a subfilter matrix, $z_{sub}$ is a delay operator matrix, and U is a Vandermon matrix. For different fractional delay $z_m^{-d}$, the amplitude response of the first and second-order FD filters based on the Lagrange interpolation method is shown in Figure 3. It can be seen that the bandwidth of the first-order FD filter is 50% of the Nyquist frequency, and the bandwidth of the second-order FD filter is 63.5% of the Nyquist frequency. It shows that the FD filter based on the Taylor series expansion Farrow structure can effectively approximate the fractional delay $z_m^{-d}$. Generally, the first-order FD is sufficient to provide sufficient bandwidth to compensate for harmonic distortion, but higher-order FD can further reduce the error. In this paper, the order M = 2.

Thus, the transfer function of FOMRC consisting of the Farrow structured FD filter based on the Taylor series is

$$G_{fomrc}=k_r\frac{z_m^{-D}{G}_d\left(z_m\right)Q\left(z_m\right)}{1-z_m^{-D}{G}_d\left(z_m\right)Q\left(z_m\right)}S\left(z_m\right)z_m^k$$

(9)

2.3. Stability Analysis of FOMRC

The steady-state error of the FOMRC system is obtained from Figure 2.

$$\begin{array}{cc}\hfill E\left(z_{\mathrm{m}}\right)=& \frac{\left({\mathrm{i}}_{ref}\left(z_m\right)-u_g\left(z_m\right)\right)}{1+\left(G_{rc}\left(z_m\right)+1\right)G_{op}\left(z_m\right)}\hfill \\ \hfill =& \frac{\left({\mathrm{i}}_{ref}\left(z_m\right)-u_g\left(z_m\right)\right){\left(1+G_{pi}\left(z_m\right)P\left(z_m\right)\right)}^{-1}\left(1-z_m^{-D}{G}_d\left(z_m\right)Q\left(z_m\right)\right)}{\left(1-z_m^{-D}{G}_d\left(z_m\right)Q\left(z_m\right)\left(1-k_r{z}_m^{k}S\left(z\right)G_{cl}\left(z_m\right)\right)\right)}\hfill \end{array}$$

(10)

where, $G_{cl}\left(z_m\right)=G_{pi}\left(z_m\right)P\left(z_m\right)/(1+G_{pi}\left(z_m\right)P\left(z_m\right))$, is the closed-loop transfer function of the system before RC is inserted. The FOMRC system is stable when the following two conditions are satisfied

➀ The roots of the $1+G_{\mathrm{pi}}\left(z_m\right)P\left(z_m\right)=0$ are in the unit circle.;

➁ $\left|z_m^{-D}{G}_d\left(z_m\right)Q\left(z_m\right)\left(1-k_r{z}_m^{k}S\left(z_m\right)G_{cl}\left(z_m\right)\right)\right|<1,\forall z=e^{j\omega T_m},0<\omega <\pi /T_m$.

It is known that stability condition 1 has nothing to do with FD. For condition 2,

$$\begin{array}{c}\hfill \left|Q\left(z_m\right)\left(1-k_r{z}_m^{k}S\left(z_m\right)G_{cl}\left(z_m\right)\right)\right|<{\left|z_m^{-D}{G}_d\left(z_m\right)\right|}^{-1},\forall z=e^{j\omega T_m},0<\omega <\pi /T_m\end{array}$$

(11)

Within the bandwidth of the FD filter, ${\left|z_m^{-D}{G}_d\left(z_m\right)\right|}^{-1}\to 1$, the stability of the system is independent of the FD filter.

2.4. Analysis of Harmonic Suppression Characteristics of FOMRC

When the sampling frequency is fixed and the grid frequency changes, $N_m$ may be a fraction. Taking the sampling factor m = 2 as an example, when the fundamental frequency $f_g$ changes from 50 Hz to 49.6 Hz and 50.4 Hz, the internal mode frequency characteristics of MRC ($N_m$ = 100) and FOMRC ($N_m$ = 100.8, 99.2) at the 8th harmonic frequency are shown in Figure 4. The open-loop gain of MRC decreased from 36.07 dB to 13.85 dB when $f_g$ changed from 50 Hz to 49.6 Hz, while FOMRC still provided higher open-loop gain, which was not affected by frequency fluctuation and had better harmonic suppression effect.

3. Modeling and Parameter Design of LCL Single-Phase Grid-Connected Inverter

As shown in Figure 5, the LCL single-phase grid-connected inverter is used as the control object for parameter design. $E_d$ is dc bus voltage, $u_i$ is inverter output voltage. $u_g$ is the grid voltage, $i_g$ is the grid current, and $i_{ref}$ is the reference current. $L_1$ and $L_2$ are filter inductance, C is the filter capacitor, and R is the series resistance of the capacitor to suppress the resonance peak generated by LCL filter.

Regardless of the equivalent resistance of the inductor, the transfer function from the input voltage $u_i$ to the grid current $i_g$ is

$$P\left(s\right)=\frac{CRs+1}{C{L}_1{L}_2{s}^3+C\left(L_1+L_2\right)R{s}^2+\left(L_1+L_2\right)s}$$

(12)

Table 1. Parameters of the LCL-type single-phase grid-connected inverter.

Parameters	Value
Inverter side inductance: $L_1$	3.8 $\mathrm{mH}$
Grid side inductance: $L_2$	2.2 $\mathrm{mH}$
Filter capacitance: C	10 $\mathrm{uF}$
Passive damping resistor: R	10 $\mathsf{\Omega }$
Dc bus voltage: $E_d$	380 $\mathrm{V}$
Grid frequency: $f_g$	50 $\mathrm{Hz}$
Sampling frequency: $f_s$	10 $\mathrm{kHz}$
Switching frequency: $f_{sw}$	10 $\mathrm{kHz}$
Switch dead time:	3 $\mathrm{us}$
Sampling ratio:	2

shows the parameters of an LCL-type single-phase grid-connected inverter. With parameters, (12) can be obtained by discretization

$$P\left(z\right)=\frac{0.006135z^2+0.004307z-0.002401}{z^3-2.005z^2+1.493z-0.4879}$$

(13)

Considering that the sampling factor m is too large, the sampling frequency of the repetitive controller will be too low, which will result in insufficient or over-compensation of phase, and affect the system stability and harmonic suppression accuracy. Therefore, the sampling factor m is set to 2. Select filter $F_1\left(z\right)=F_2\left(z\right)=0.15z^{-1}+0.7+0.15z$, the cutoff frequency is $2{\pi }^{*}2460<2{\pi }^{*}2500=\pi /T_m$.

According to [12], the parameters of the PI controller are designed and calculated. $k_p=10,k_i=1300$. $G_{cl}\left(z_m\right)$ has a amplitude margin of $9.95\mathrm{dB}$ and a phase margin of ${154}^{\circ }.Q\left(z_m\right)$ usually selects the low-pass filter with a gain of less than 1 or a constant of less than 1 to improve the stability of RC. In this paper, $Q\left(z_m\right)=0.25z_m+0.5+0.25z_m^{-1}$. In order to make $G_{cl}\left(z_m\right)$ have faster amplitude attenuation at high frequencies, select the fourth-order Butterworth low-pass filter $S\left(z_m\right)$ with a cut-off frequency of $1\mathrm{kHz}$. The expression of $S\left(z_m\right)$ is as follows

$$S\left(z_m\right)=\frac{0.04658z_m^4+0.1863z_m^3+0.2795z_m^2+0.1863z_m+0.04658}{z_m^4-0.7812z_m^3+0.68z_m^2-0.1827z_m+0.03012}$$

(14)

In order to determine the phase lead order k, the Bode diagram of $z_m^{k}S\left(z_m\right)G_{cl}\left(z_m\right)$ at $3,4,5,6$ is given in Figure 6. It can be seen that when $k=4$, the phase of $S\left(z_m\right)G_{cl}\left(z_m\right)$ is close to $0^{\circ }$ within $1\mathrm{kHz}$, and the compensation effect is the best.

Finally, select $k_r$ according to condition 2. According to (11)

$$\begin{array}{c}\hfill \left|Q\left(z_m\right)\left(1-k_r{z}_m^{k}S\left(z_m\right)G_{cl}\left(z_m\right)\right)\right|<1,\forall z=e^{j\omega T_m},0<\omega <\pi /T_m\end{array}$$

(15)

Define $H\left(z_{\mathrm{m}}\right)=Q\left(z_m\right)\left[1-k_r{z}_m^{k}S\left(z_m\right)G_{cl}\left(z_m\right)\right],z_m=e^{j\omega T_m}$. Then when the locus of $H\left(e^{j\omega T_m}\right)$ is in the unit circle, the system is stable. The internal model gain $k_r$ directly determines the error attenuation speed of the repetitive control system, and also affects the Nyquist curve of $H\left(e^{j\omega T_m}\right)$. It can be seen from Figure 7 that when $k_r\le 1.5$, the trajectory of $H\left(e^{j\omega T_m}\right)$ is always in the unit circle, and there is still a certain distance between the trajectory and the boundary of unit circle, indicating that the system has sufficient stability margin. The direction of the arrow represents the trajectory direction of the amplitude curve of $H\left(e^{j\omega T_m}\right)$ from 0 Hz to Nyquist frequency. When $k_r$ is 0.5, 1, 1.5, and 2, respectively, the corresponding points of the amplitude of $H\left(e^{j\omega T_m}\right)$ at 150 Hz and 1000 Hz are shown in the figure. It can be seen that when $k_r=1$, the low frequency region of $H\left(e^{j\omega T_m}\right)$ is closer to the center of the unit circle, and the controller has better low-frequency signal tracking ability in the low frequency region.

4. Simulation Analysis

In order to verify the performance of the FOMRC system, the LCL-type single-phase grid-connected inverter model is built in MATLAB / Simulink simulation environment, and the steady-state and dynamic performance are compared with CRC. The reference current is $10\mathrm{A}$. The CRC system, that is, sampling factor $m=1$, sampling frequency $f_s=10\mathrm{kHz}$, delay beat $N=f_s/f_g$, phase lead order $k=8$, other parameters remain unchanged.

4.1. Steady-State Performance

As can be seen from Figure 8, when $f_g=50\mathrm{Hz},N=T_g/T_s=200,N_m=T_g/T_m=100$, the output voltage and current waveforms of CRC and FOMRC are almost sinusoidal, the THD of CRC and FOMRC are $0.97\%$ and $1.33\%$, respectively, both of them can effectively suppress the harmonic in the power grid.

Figure 9 and Figure 10 show that when $f_g=49.6\mathrm{Hz},N=T_g/T_s=201.6,N_m=T_g/T_m=100.8;$ when $f_g=$$50.4\mathrm{Hz},N=T_g/T_s=198.4,N_m=T_g/T_m=99.2$. The voltage and current waveform and output current spectrum of CRC and FOMRC control systems under fractional delay are obtained. When $f_g=49.6\mathrm{Hz}$, the THD of CRC and FOMRC are $2.71\%$ and $1.26\%$, respectively. When $f_{\mathrm{g}}=50.4\mathrm{Hz}$, the THD of CRC and FOMRC are $3.33\%$ and $1.19\%$, respectively. It can be seen that when the fractional delay occurs, the output current waveform of CRC is seriously distorted and the harmonic content is high. The output current waveform of FOMRC is still a smooth sine curve, and the harmonic suppression ability is significantly improved compared with CRC.

Currently, a number of scholars have used a third-order FIR filter based on Lagrange interpolation to approximate the fractional delay [16,17,18] and achieved good harmonic suppression results. To verify that the FOMRC proposed in this paper has similar performance, Figure 11 gives a comparison plot of the THD results of CRC, MRC, and MRC using FIR to approximate the fractional delay, and the FOMRC proposed in this paper when the grid frequency fluctuates in the range of $50\pm 0.4\mathrm{Hz}$. It can be seen that the THD of MRC is slightly reduced compared with CRC due to the halving of the sampling rate when the grid frequency is 50 Hz. The larger the grid frequency offset, the more serious the degradation of the harmonic rejection ability of MRC. The proposed method can obtain smaller THD values, better frequency adaptivity, and the same anti-frequency fluctuation capability as the method using FIR filter with approximate fractional delay.

4.2. Dynamic Performance

Figure 12 and Figure 13 show the dynamic performance of CRC and FOMRC control systems with reference current amplitude reduced from 10 A to 6 A at t = 0.25 s. When the fundamental frequency is 49.6 Hz, the current error of the CRC system is ±1 A, and that of the FOMRC system is ±0.3 A. When the fundamental frequency is 50.4 Hz, the current error of the CRC system is ±1.1 A, and that of the FOMRC system is ±0.3 A. The errors of the two control systems can converge again in two fundamental periods when the current changes abruptly at different frequencies. It can be seen that FOMRC and CRC have similar error convergence speeds when the power grid frequency fluctuates, but can effectively reduce the error.

5. Conclusions

When the grid frequency fluctuates, the performance of the grid-connected inverter based on MRC will decrease. Therefore, this paper proposes FOMRC based on the Farrow structure FD filter to ensure that the resonant frequency of MRC can track the grid frequency in real time to obtain frequency adaptability.

The FOMRC system maintains a high open-loop gain at the fundamental and harmonic frequencies with variable grid frequency, allowing it to provide higher tracking progress and harmonic suppression performance, thus the system has a higher tracking gain and a lower THD.

What’s more, by setting the sampling factor to 2, the sampling frequency of the repetitive controller part is reduced to half of the original, which reduces the calculation. The simulation results verified the effectiveness of the proposed FOMRC scheme in terms of fast transient response, high tracking accuracy, and immunity to frequency variations.