Open access peer-reviewed chapter

Harmonic Mitigation for VSI Using DP-Based PI Controller

Written By

Ragini Meshram, Monika Madhusoodan and Sushama Wagh

Submitted: April 2nd, 2017 Reviewed: January 29th, 2018 Published: November 5th, 2018

DOI: 10.5772/intechopen.74587

Chapter metrics overview

862 Chapter Downloads

View Full Metrics


The chapter focuses on the strengths of dynamic phasor (DP) based model over conventional time domain model and the controller designed using it for selected harmonic mitigation. To validate the effectiveness of the controller designed using DP-based model, the single-phase voltage source inverter (VSI) with various loads has been considered including effects of intermittent nature of renewable energy sources e.g. photovoltaic module. The DP technique offers distinct advantages in modelling, simulation, and control with respect to the time domain models. With the assets of DP modelling technique based on the measurement of harmonic coefficients, the PI controller is designed which eliminates selected voltage and current harmonics in VSI and results are compared with the repetitive control technique. It has been proved through simulations that as compared to conventional technique, the proposed DP-based PI controller eliminates multiple selected frequencies effectively.


  • dynamic phasor (DP)
  • harmonics
  • photovoltaic (PV) inverter
  • repetitive control (RC)
  • voltage source inverter (VSI)

1. Introduction

In modern power system, the growth of photovoltaic (PV) due to increase in efficiency, clean sources of energy and decrease in cost of solar technology promotes substantial power integration. According to their functionality and operating requirements, it is classified in standalone and grid connected [1]. In standalone system, remote area is supplied by dc or ac power with converters and energy storage devices, on the other hand, in grid connected system the generated power is supplied to the utility services without any energy storage equipments. Distributed generation (DG) systems making use of renewable energy sources are being designed and connected to grid. In distribution feeder, the high penetration of PV generators has considerable impact on the system behaviour. For better insight into these impacts development and fine tuning of existing power system and power converters in the simulation environment is essential. For fast computation of large distribution systems, the dynamical models for various components and various stages are essential. The chapter concentrates on development of the dynamic phasor (DP)-based model of pulse width modulated (PWM) converters that are capturing the transients of interest. The topology of PV system considered, consist of, three power stage circuit stages, such as dc/dc boost converter for maximum power point tracking (MPPT), low voltage single phase dc/ac inverter and filter inductance with grid.

In grid-connected PV inverter systems, the important power quality issue is current distortion [2, 3, 4]. According to IEEE Standard 1547–2003 and IEC standard 61,727, the total harmonic distortion (THD) for the grid current is lower than 5% to avoid unfavourable effects on the other equipments that are connected to grid. In addition, for odd harmonics from 3rd to 9th, the limitation is 4% and even harmonics are limited to 25% of the odd harmonic limits [5]. These acceptable current distortion levels are defined for grid-connected PV inverter systems in rated operational mode. Since, solar radiations are not steady in PV inverter systems, the output current of PV is less than its rated value.

Many filter topologies are used in literature for reduction of harmonics, the higher order power filter named as LLCL filter, which inserts a small inductor in the branch of the capacitor in the traditional LCL filter to create a series resonant circuit at the switching frequency [4, 6]. It can attenuate current ripple component of switching frequency much better than LCL filter [7], also decrease in total inductance which raise the characteristic resonance frequency for better operation of inverter system control [8, 9, 10, 11].

For analysis of harmonics in PWM systems and the harmonics generated by switching, the DP-based model is best available tool as explained in [18]. The DP-based model is developed from generalised averaging technique and is capable of converting periodic varying state variables into dc state variables which is a widely employed method for modelling of oscillatory systems [12, 13]. Contrary to different types of existing numerical methods available in literature, the DP-based technique [14, 15] focuses on the frequencies of interest and provides accurate simulations for larger time steps [16].

The chapter highlights the important advantages of DP-based models to capture the harmonic content in the signals. The harmonics of the DP model is represented by Fourier coefficients which are determined by harmonic balance equation which contains useful information according to the selected harmonics. This information can be used for better controller design as well as deeper analysis on system performance according to oscillating dynamics and its interaction within and with system. A DP-based model can also be formed using sequence components of alternators and its voltage controller can be designed as mentioned in [17].

Control strategies that are applied to the power electronics converter could support the power conditioning functionality of the current and grid-connected inverters. With respect to the current control, proportional resonant (PR) controller is the most popular controller for single-phase PV inverters. It is well known that the odd harmonics (3rd, 5th and 7th) are dominant in the output current of single-phase grid-connected inverters. Since PR controller cannot reject all harmonics appearing in the grid current, the PR controlled inverters may not be able to feed high quality currents into the grid. In order to eliminate the current harmonics distortion effectively, repetitive controller (RC) can be used with phase compensator to track and eliminate all the harmonics in the system.

The reduction in THD indicates the accuracy with which the periodic signal is tracked at the steady state. For achieving precise tracking of periodic signal, one of the most preferred approaches is the RC [19, 20, 21]. The fundamental operating principle of the RC is based on internal model principle [22], where it observes the systems’ periodic signals for one-cycle period and a corresponding compensating signal is generated in next cycle period to ensure precise tracking of the signal [23]. From the frequency point of view, the RC performs finer error cancellation for periodic signals by presenting a large magnitude of loop gain at the fundamental frequency and its integral multiples. This can be identified as a form of period-based integral control [24, 25, 26, 27, 28].


2. Background preliminaries

2.1. Concept of repetitive control

The main objective of RC is to asymptotically track the reference signal while rejecting disturbances. Any periodic signal with period T can be generated by a positive feedback system with the specified initial function as shown in Figure 1(a). According to Internal Model Principle [29], it is necessary to include model of Figure 1(a) into closed loop system in order to achieve perfect tracking or external disturbance rejection. In case of inverter control, when reference signal is the periodic signal, a pair of conjugate imaginary poles should be included at the frequency win the closed loop which can be provided by 1/(s2 + w2), as used in proportional resonant controller scheme. The inverter output voltage is full of harmonics and it is essential to reduce the level of concern harmonics to obtain low THD. One of the solutions is to include bundle of pairs of conjugate imaginary poles at different frequencies. Another solution could be so-called RC [30, 31, 32], which adopts infinite dimensional internal model to provide series of conjugate poles at all concerned harmonic frequencies. The stability and robustness of RC system is improved by introducing a low pass filter Q(s) in series with time delay element. Any low pass filter satisfying Qs<1can be used in minimum phase systems, but there are bandwidth restrictions on Q(s) for non-minimal phase plant [30].

Figure 1.

Repetitive control system (a) periodic signal generator (b) general structure of plug-in repetitive control (c) closed-loop control without RC.

The design and synthesis methods for modified RC system vary with different configurations as mentioned in [25, 26, 27, 28, 31]. The plug-in RC system depicted in Figure 1(b) is most commonly used structure. The design problem is mainly to choose and optimise the dynamic compensator Gf (s) and the low pass filter Q(s). The choice of controller parameters involves a trade-off between steady-state accuracy, robustness, and transient response of the system. The closed-loop control of single-phase VSI without RC is shown in Figure 1(c), where R(s) is the reference signal, D(s) is the disturbance, Gp(s) is the plant and Gc(s) is the closed-loop controller. The closed-loop transfer function H(s) of the control system without RC is represented in (1) and the tracking error E0(s) of the signal due to periodic input R(s) and disturbance D(s) is given by


To assure the effectiveness of RC for accurate tracking, the system tracking error E(s) with RC is written in terms of tracking error E0(s) of the original system without RC and is expressed as


To assure stability, the dynamic compensator must be chosen such that


Eq. (2) can be expressed as the sum of a geometric progression


where eTsimplies a one-cycle delay of fundamental period. Eq. (4) indicates that the first cycle error e0(t) = L−1[E0(s)] for t ∈ [0, T] is unaffected by RC. If


then the error can be simplified to Es=E0s1QseTs, which is the error of second cycle. If Q(s) can be achieved by


where w = 2π/T, then the steady-state tracking will converge to zero at each of the harmonics. However, because of physical nonideality of the system, the design requirements (5), (6) cannot be practically satisfied for all of the harmonics, rather it is more appropriate to adopt the conditions


thus, the task of designing the RC is to meet the harmonic elimination requirements of (7) while satisfying stability conditions (1) and (3).

2.2. The dynamic phasor-based model

The voltages and currents in power electronic converters and electrical drives are typically periodic in steady state, and mostly nonsinusoidal. The dynamics of interest for analysis and control are often those of deviations from periodic behaviour, for instance as manifested in deviations of the envelope of a quasi-sinusoidal waveform from its steady-state value. For analysis of the steady state, generally used methods are phasor, harmonic, and describing function [12, 13, 18]. With analytical approach reviewed, the main aim of this research is systematic development of phasor dynamics, from which the dynamic behaviour of the original waveform or its envelope can efficiently be deduced.

2.2.1. Definition and concepts used in DP model

The DP modelling technique is based on the generalised averaging method wherein the complex time domain signal is represented in the interval (t‐T, t) with the Fourier series representation. The time varying Fourier coefficients are termed as DPs, which is essentially a frequency-domain method where various harmonics developed over time can be decomposed. In power system steady-state analysis, these Fourier coefficients Xk (t) are also called as phasors. During transients, the system is not in a pure periodic state but nearly periodic state. The idea is now to extend this approach to nearly periodic signals [13] and to approximate x(τ) in the interval t∈ (t−T, t] with a Fourier series representation given as


In this representation, as the signal x(τ) is nearly periodic and since the interval under consideration slides as a function of time, the Fourier coefficients Xk (t) are time varying. This expression can also interpreted as an orthogonal signal expansion of the function x(t) with orthonormal basis ejkωs t. The kth phasor (Fourier coefficient) of X is defined as:


where ws = 2π/Ts. Note that the phasors are defined over a moving time-window and hence, time dependent. Additionally, if X is periodic with time Ts, then Xk (t) is constant. In the DP approach, only a few coefficients provide a good approximation of the original waveform. Some important properties of DPs are-

  1. The relation between the derivatives of x(τ) and the derivatives of Xk (t) is given by


  2. The product of two time-domain variables equals a discrete time convolution of the two DP sets of the variables, which is given by


  3. For a complex valued signal x, the relationship between Xk and X−k is given as Xk=Xk

2.3. Advantages and limitations of dynamic phasors

The DP approach offers a numbers of advantages over conventional methods.

  1. The oscillating waveforms of ac circuits become constant or slowly-varying in the DP domain and different frequency components can be handled separately with convenience.

  2. It approximates a periodically switched system with a continuous system, thereby converting periodic varying state variables into dc state variables.

  3. At steady state the DP Xk becomes constant.

  4. As the variations of DP Xk are much slower than the instantaneous quantities x(t), they can be used to compute the fast electromagnetic transients with larger step sizes, so that it makes simulation potentially faster than conventional time domain EMTP-like simulation.

  5. The selection of DP index-k gives a wider bandwidth in the frequency domain than traditional slow quasi-stationary models used in transient stability programs.

  6. The selection and variation of k also gives the possibility of showing coupling between various quantities and addressing particular problems at different frequencies.

The only disadvantage of the DP approach is that the number of variables and equations are higher than in the original equations.


3. DP-based modelling

3.1. Single-phase inverter with R-load

The typical structure of a single-phase inverter is presented in Figure 2 which includes H-bridge with switching function Si that is fed by LVDC link voltage (VdcL) and LC circuit is connected between inverter output and load resistance (RL) for filtering purposes. LC filter consists of filter inductor Lfi, filter resistance Rfi, and filter capacitor Cfi.

3.1.1. System modelling of inverter with R-load

Figure 2.

Single-phase inverter with LC output filter.

The dynamic model of inverter can be represented as:


where, Si (t) is a switching function that denotes the switching status as:


In this case, the switch is controlled using fixed frequency switching, i.e. time axis is divided into intervals [kTs, (k + 1)Ts], where Ts > 0 is the switch period and k ∈ N. The duty cycle dk ∈ [0, 1] is chosen at the beginning of each switch. The duty cycle determines the fraction of time the switch is active (1 & 0 represents the on-mode and off-mode respectively) and thus controls the system dynamics. The aim of controller is to regulate inverter output voltage (Vcfi) at desired value of voltage, which is a pure sinusoidal signal with angular frequency wand amplitude Vp.

3.1.2. DP model of inverter with R-load

The DP-based dynamic model for inverter topology shown in Figure 2 has been developed based on theory explained in Section 2.2. In order to obtain an accurate approximation, the foremost step is to fix Fourier coefficients for the DP model as specified in Table 1.

State variablesDP index kDiscontinuous variablesDP index k
ifi, Vcfi1, 3, 5Si1, 3, 5

Table 1.

Fourier coefficients of the system.

Converting (12) and (13) into DP form, the state variables, inductor current and capacitor voltage are modified to


Eqs. (15)(16) represented in real and imaginary coefficients as


(i) Modelling of switching function Si

Let the duty cycle di and switching period T are the variables associated with a switch and H-bridge of inverter consists of 4 switches Qi1−Qi4 (with reference to Figure 2). Representing state 1 in Si (t) for Qi1 and Qi2 ON while state −1 for Qi3 and Qi4 ON, the switching function (Si (t)) and its DP form Sikis defined in (21).


Substituting (21) into (17)(20), and organising dynamical equations into state-space representation will result in ddtX=AX+Bu, where


The state Eqs. (22)(25) implies that the dynamics of single-phase inverter represented using index-1, 3 and 5 of output filter inductor current (ifi) and index-1, 3 and 5 of output filter capacitor voltage (Vcfi) as state variables. If more indices-k were used, then the model would be more accurate, however, the resulting model would be too complex to provide insightful information for controller design.

3.2. Single-phase inverter with RL-load

The typical structure of a single-phase inverter is same as presented in Figure 2 except the load part, where R is replaced by combination of RL and inductance LL as a RL load. State variables for this case are output filter inductor current (ifi), output filter capacitor voltage (Vcfi), and load current (iL).

The dynamic model of inverter with RL load can be represented as:


The DP form of (26)(28) for the state variables inductor current, capacitor voltage and load current are:


Using the switching function model derived in Section 3.1 the detailed DP model can be summarised as:


The state space representation ddtX=AX+Bucan be obtained by substituting (21) into (32)(37), where


3.3. Single-phase grid-connected inverter for PV system

The configuration of system consist of three stages, the boost converter stage which is used for tracking of maximum power of photovoltaic array, the inverter stage which generates the behaviour of a sine square wave and the filtering step through L arrangement with grid is represented in Figure 3. The input voltage VPV is from photovoltaic array and the system employs signals at different frequencies:

  1. The boost converter with switch μhaving switching frequency (fc), which is of kHz.

  2. The switching function (S) constitutes inverter switching frequency (fi) in range of kHz.

  3. The main systems frequency (f) of 50 Hz.

Figure 3.

Single-phase grid-connected PV inverters power circuit topology.

3.3.1. System modelling of grid-connected PV inverter

In Figure 3, VPV is the input voltage to the system from photovoltaic renewable source where iL is a input current to boost converter, VC is the output voltage of the boost converter and also the input voltage to the inverter, μis the control signal of the boost converter. In inverter output side iLf is the current passing through the filter inductor while VAC = Esin(2πft) is the grid voltage and switching function for all inverter switches (S1 to S4) is S.

For the topology considered in Figure 3, different operating states can be obtained depending on the position of switches, μand S. The dynamical model of single-phase grid-connected PV system is obtained in (42)(44).


3.3.2. DP model of grid-connected PV inverter

Similar procedure is followed for the model development of PV inverter system which is shown in Figure 3. The considered Fourier coefficients for the system are specified in Table 2. Converting (42)(44) into DP domain, the revised model is represented in (45)(47).

State variablesDP index kDiscontinuous variablesDP index k
iL, VC0, 2, 4μ0, 2, 4
iLf1, 3, 5S1, 3, 5

Table 2.

Fourier coefficients of the system.

(i) Modelling of the boost converter switching function μ.

Let, the duty cycle d and switching period Ts be the variables associated with a boost converter switch and state 1 and 0 represents ON and OFF mode of the switch respectively. The switching function (μ) and its DP form μkis provided in (48).


(ii) Modelling of the inverter switching function S.

Let, the duty cycle di and switching period T be the variables associated with a inverter switch and H-bridge inverter consisting of 4 switches S1 to S4 (as shown in Figure 3). Using state 1 to represent S1 and S2 ON while state −1 for S 3 and S4 ON, the switching function (S(t)) and its DP notation Skcan be defined as (49).


Expanding (45)(47) for every coefficient of Table 2, a system dynamics can be obtained as:


Corresponding to the boost converter switching frequency, for both current (iL) and voltage (VC) signals, the DPs are chosen as the direct current (k = 0) and the second harmonic (k = 2). On the other hand, for inverter, the DP form of the filter current (iLf) is chosen as fundamental and the third harmonic frequencies (k = 1, 3). In addition, if switching harmonics needs to be accommodated in the model, the DP index should be extended to k = mfc for boost converter and k = mfi for inverter.


4. DP-based PI controller design

4.1. Development of PI controller using DP model

To design a PI controller using DP-based small signal model, for maintaining constant desired voltage at the output of the inverter, a nested controller is proposed as shown in Figure 4 with inner current and outer voltage control loops.

Figure 4.

PI control scheme for DP model of single-phase inverter.

Controller design and stability analysis for inverter requires the deviation of small-signal control-to-output transfer function, which is the dynamic response for small perturbation in control signal. In state-space model of single-phase inverter, when there is a small perturbation in di, all the state variables will deviate from their steady states. Assume that the input voltage VdcL0is constant, the capitalised variables represents steady state values, lower case variables are the actual states and Δ variables represents small-signal states. The deviation of state variables are defined as:


The small signal model of the single-phase inverter is given as: ddtΔXinv=AinvΔXinv+BinvΔdiwhere,


Eqs. (56)(58) can be used to calculate the control-to-output transfer function of single-phase inverter.

The real and imaginary references for the voltage loop are generated by phase-locked-loop (PLL) with the desired ac voltage. The outer voltage loop generates the current references for inner current loop whereas the inner current loop generates desired duty ratio.

4.1.1. Design of the inner current loop controllers

The inner current loop dynamics is described as


In (59), ΔifikIand ΔifikRare considered as perturbation which are further eliminated by the decoupling structure.

Describing transfer function for both ΔifikRand ΔifikI


the current control loop is built using a PI controller HPIC (s) and the closed loop transfer function is H0C (s) as given in (61)


To compute TiC in terms of T0C and hence, KpC for given fixed damping coefficient ξC, imposing a closed-loop transfer function of the form (61) leads to (62).


4.1.2. Design of the outer voltage loop controllers

The dynamics of outer voltage loop is described as


The transfer function from current component to voltage is


The voltage control loop is built using PI controller HPIV (s) and the voltage closed loop transfer function H0V(s) is given in (65)


By imposing closed loop transfer function of the form (65) for fixed damping coefficients ξV, the time constant TiV can be investigated in terms of T0V and hence, KpV is calculated as (66)


It should be noted that the outer voltage control loop time constant is larger than inner current control loop.


5. Performance evaluation of harmonic mitigation for VSI

Using simulation parameters as shown in Table 3, the models for single-phase inverter and grid-connected PV inverters are developed in MATLAB/Simulink and effectiveness of controller is validated for the single-phase inverter with R-load.

Single-phase inverter with R and RL loadGrid-connected PV inverter
Rfi = 0.07 WLfi = 1.07 mHCfi = 30mFL = 15 mHLf = 20 mHC = 20mF
RL = 35 WLL = 69.6 mHVPV = 200 VVAC = 230Vrms

Table 3.

Simulation parameters for harmonic mitigation of VSI.

5.1. Single-phase inverter with R and RL load

5.1.1. Response of VSI using conventional time domain model

Using input voltage to the inverter as 400 V with constant modulation of 0.8, the response obtained by conventional time domain model without controller is shown in Figure 5 and Figure 6 for only R and RL load respectively. It can be seen from Figures 5 and 6(a) and (b), irrespective of type of load R or RL, the inverter output filter current is more distorted as compared to output filter voltage which indicates the necessity of proper controller to eliminate unwanted harmonics in voltage or current. Figures 5 and 6(c) and (d) shows switching signal and discontinuous output voltage of inverter respectively.

Figure 5.

Performance analysis of single-phase inverter with R-load but without controller. (a) Output filter voltage of inverter (Vcfi) (b) output filter current of inverter (Ifi). (c) Switching signal (d) discontinuous output voltage.

Figure 6.

Performance analysis of single-phase inverter with RL-load but without controller. (a) Output filter voltage of inverter (Vcfi) (b) output filter current of inverter (Ifi). (c) Switching signal (d) discontinuous output voltage.

5.1.2. Response of VSI using DP model

In this section, open loop analysis of simulation results using DP model are built within the SimPowerSystem of the MATLAB/Simulink library with the change in load resistance from 35 Ω to 25 Ω at 0.2 sec. DP coefficient dynamics of state variables of inverter are presented in Figures 7 and 8. With DP index k = 1, 3, 5 inverter filter output voltage and current is shown in Figures 7 and 8(a) and (b) for R and RL load respectively.

Figure 7.

Response of single-phase inverter with R-load using DP-based model without controller(a) output filter voltage (Vcfi) (b) output filter current (Ifi) (c) re-coefficients of DP Vcfi (d) Im-coefficients of DP Vcfi (e) re-coefficients of DP Ifi (f) Im-coefficients of DP Ifi (g) re-coefficients of DP Si (h) Im-coefficients of DP Si.

Figure 8.

Response of single-phase inverter with RL-load using DP-based model without controller. (a) Output filter voltage (Vcfi) (b) output filter current (Ifi) (c) re-coefficients of DP Vcfi (d) Im-coefficients of DP Vcfi (e) re-coefficients of DP Ifi (f) Im-coefficients of DP Ifi (g) re-coefficients of DP iL (h) Im-coefficients of DP iL (i) re-coefficients of DP Si (j) Im-coefficients of DP Si.

The corresponding real and imaginary parts of harmonic coefficients are separately shown in Figures 7 and 8(c) and (d) for inverter filter voltage and Figures 7 and 8(e) and (f) for load current for R and RL load respectively. The switching function is also separated in harmonic coefficients that are shown in Figures 7 and 8(g) and (h).

It should be noted that as compared to conventional time domain model, DP-based model provides better insight in the features of output current and voltage waveforms. Using this valuable precise details, the controller could be designed to target selected harmonics.

5.2. Response comparison of single-phase grid-connected inverter for PV system using time domain and DP-based model

The single-phase grid-connected PV model is simulated in MATLAB/Simulink using input PV voltage as 200 V. The system consists of three stages that is boost stage, inverter stage, and filter stage. The variable of interest for harmonic analysis is the filter current in inverter stage and Figure 9 shows the inverter filter current with harmonics. Same single-phase grid-connected PV system response is observed using DP-based model using higher index by plotting harmonic coefficients separately for output filter current as shown in Figure 10.

Figure 9.

Response of inductor filter current (iLf) of grid-connected PV inverter.

Figure 10.

Response of inductor filter current (iLf) of grid-connected PV inverter using DP-based model. (a) Re-coefficients of DP iLf and (b) Im-coefficients of DP iLf.

5.3. Performance analysis of VSI with controller

Again it should be noted that DP-based model provides more accurate information about harmonic content as compared to time domain model. In both models, it is clearly seen that the controller is highly recommended for mitigation of dominant harmonic components. As a result, a conventional RC is designed and used with time domain model. On the other hand, a proposed DP-based PI controller is also designed to show its effectiveness as compared to RC in harmonic mitigation. The aim of the repetitive controller is to track the sinusoidal reference voltage of 230Vrms and eliminate all harmonics present in output voltage of inverter system. In Figure 11(a), the tracking of output voltage signal (Vcfi) to the desired reference signal with load resistance of 100 Ω is observed and shows the elimination of harmonics present in Figure 5(a). The corresponding current harmonics are also eliminated with this controller as shown in Figure 11(b) which were present in Figure 5(b). A DP-based PI controller is designed with nested form as shown in Figure 4 with outer voltage loop and the inner current loop. The controlled DP coefficients of inverter filter output voltage are shown in Figure 12(a) and (b) and the actual filter output voltage (Vcfi) is shown in Figure 12(c). Though, PI gains can be calculated by numerous techniques like loop shaping, pole placement, root locus, prediction-based method etc., in the proposed controller the gains are computed as explained in Section 4, specifically designed to eliminate 3rd and 5th harmonic coefficients.

Figure 11.

Performance analysis of single phase inverter with repetitive controller. (a) Inverter output filter voltage (Vcfi) and (b) inverter output filter current (Ifi).

Figure 12.

Performance analysis of DP single-phase inverter model with PI controller. (a) Re-coefficients of DP Vcfi (b) Im-coefficients of DP Vcfi (c) controlled filter output voltage (Vcfi).

Comparing response of VSI with RC- and DP-based PI controller Figures 11 and 12, it can be observed that proposed controller has achieved required quality of output voltage by eliminating selected harmonics. It was possible to target individual harmonic components because of higher index DP model which gave detailed precise information about the elements responsible for deviation from desired quality output. It should be noted that, irrespective of harmonics, the DP-based controller helps in providing better voltage regulation for various disturbances such as intermittent nature of renewable input sources and load variations [32].


6. Conclusions

In this chapter, the DP technique is used to model single-phase inverter with R-load, RL-load and grid-tied PV inverters with multiple frequencies. The overall design aims to obtain acceptable signal approximations with the selected DP coefficients. The harmonics coefficients (DP index-k) for each model are chosen after analysing the main harmonic from the original systems. These DP-based equations have been introduced to compute total, fundamental and distorting dc/ac voltage and currents. In view of this, the chapter has proposed two major contributions: strengths of DP-based model over time domain model and effectiveness of DP-based PI controller over RC. The simulation results comparing various circuit topologies in time domain and DP-based model has shown that the DP-based model provides more detailed information about harmonic contents in the output variables. Using this pinpoint knowledge of selected harmonics of concern, a DP-based controller with higher index can be designed to precisely target the required response which is otherwise not possible in general time domain model. As a result, the DP-based PI controller has shown expected quality improvements as compared to well-known RC controller in voltage regulation and selected harmonic elimination in both steady-state and transient state. It could be also useful in designing best possible operations by minimising losses in system.



The authors would like to acknowledge Dr. N. M. Singh, Department of Electrical Engineering, VJTI, Mumbai and Dr. A. M. Stanković, Department of Electrical and Computer Engineering, Tufts University, USA for their continuous guidance throughout this research.


Procedure to be adopted for three-phase case study:

The switching function of the switch Sjk is defined as

Sjk=1,switchSjkis closed0,switchSjkis openE67

where j represents the phase i.e. j ∈ {a, b, c} and variable k ∈ {p, n} is used to distinguish the upper and lower component of phase leg. The switching constraints for upper switch Sjp and the lower switch Sjn is Sjp + Sjn = 1.

The phase switching function defines if the output terminal of the phase leg (point j) is connected to the positive potential or to the negative potential of the dc link voltage by the single-pole double-throw switch Sj. The phase switching function (Sj) is defined as follows

Sj=1,Sjis connected to pointp1,Sjis connected to pointnE68

Assuming the case of three-phase inverter with R-load, the notations will be followed as given in Section 3.1. The dynamic equations are as:


where ifiab=ifiaifib, ifibc=ifibific, ifica=ificifia, Siab=SiaSibc, Sibc=SibSic, Sica=SicSia, Vcfiab=VcfiaVcfib, Vcfibc=VcfibVcficand Vcfica=VcficVcfia.

In Sine Pulse Width Modulation (SPWM) technique, Sij can be replaced by its fundamental component dij


where j = a, b, c, ϕa = ϕ0, ϕb = ϕ0 + 2π/3, ϕc = ϕ0 − 2π/3, ϕ0 is initial phase angle and m is amplitude modulation ratio.

The state variables in of (69) and (70) can be approximated in DP form as sum of the their DC components and index-k components


The input dc voltage can be considered as a constant VdcL, thus it has


where M=sin2πdijϕjand N=cos2πdijϕj1.


  1. 1. Baekhoej KS, Pedersen JK, Blaabjerg F. A review of single-phase grid-connected inverters for photovoltaic modules. IEEE Transactions on Industry Applications. 2005;41(5):1292-1306
  2. 2. Xu Jinming, Qiang Qian, Binfeng Zhang, Huizhen Wang, and Shaojun Xie. An adaptive algorithm for grid-connected inverter to suppress current harmonics and instabilities due to grid impedance and distortion. In: Applied Power Electronics Conference and Exposition (APEC), 2017 IEEE;2017. IEEE. pp. 2638-2643
  3. 3. Aleksandr R, Simoes MG, Al-Durra A, Muyeen SM. LCL filter design and performance analysis for grid-interconnected systems. IEEE Transactions on Industry Applications. 2014;50(2):1225-1232
  4. 4. Wu W, He Y, Blaabjerg F. An LLCL power filter for single-phase grid-tied inverter. IEEE Transactions on Power Electronics. 2012;27(2):782-789
  5. 5. Yang Y, Zhou K, Blaabjerg F. Current harmonics from single-phase grid-connected inverters examination and suppression. IEEE Journal of Emerging and Selected Topics in Power Electronics. 2016;4(1):221-233
  6. 6. Grahame HD, McGrath BP. Opportunities for harmonic cancellation with carrier-based PWM for a two-level and multilevel cascaded inverters. IEEE Transactions on Industry Applications. 2001;37(2):574-582
  7. 7. Erika T, Holmes DG. Grid current regulation of a three-phase voltage source inverter with an LCL input filter. IEEE Transactions on Power Electronics. 2003;18(3):888-895
  8. 8. Miguel C, Miret J, Matas J, Vicuna LG d, Guerrero JM. Control design guidelines for single-phase grid-connected photovoltaic inverters with damped resonant harmonic compensators. IEEE Transactions on Industrial Electronics. 2009;56(11):4492-4501
  9. 9. Abhijit K, John V. Mitigation of lower order harmonics in a grid-connected single-phase PV inverter. IEEE Transactions on Power Electronics. 2013;28(11):5024-5037
  10. 10. Marco L, Teodorescu R, Blaabjerg F. Multiple harmonics control for three-phase grid converter systems with the use of PI-RES current controller in a rotating frame. IEEE Transactions on Power Electronics. 2006;21(3):836-841
  11. 11. Cardoso MF, Spaduto RR, Freitas LCG d, Tavares CE, Macedo JR, Rezende PHO. Harmonic distortion analysis in a low voltage grid-connected photovoltaic system. IEEE Latin America Transactions. 2015;13(1):136-142
  12. 12. Sanders Seth R, Mark Noworolski J, Liu XZ, Verghese GC. Generalized averaging method for power conversion circuits. IEEE Transactions on Power Electronics. 1991;6(2):251-259
  13. 13. Stefan A, Mariéthoz S, Morari M. Dynamic phasor model predictive control of switched mode power converters. IEEE Transactions on Control Systems Technology. 2015;23(1):349-356
  14. 14. Seddik B, Munteanu I, Bratcu AI. Power electronic converters modeling and control. Advanced Textbooks in Control and Signal Processing. 2014;454:454
  15. 15. Alberto C-M, Pérez-Cisneros MA, Domínguez-Navarro JA, Osuna-Enciso V, Zúñiga-Grajeda V, Gurubel-Tun KJ. Dynamic phasors modeling for a single phase two stage inverter. Electric Power Systems Research. 2016;140:854-865
  16. 16. Parimi M, Monika M, Rane M, Wagh S, Stankovic A. Dynamic phasor-based small-signal stability analysis and control of solid state transformer. IEEE 6th International Conference on Power Systems (ICPS). IEEE. pp. 1-6;2016
  17. 17. Zhixin M, Piyasinghe L, Khazaei J, Fan L. Dynamic phasor-based modeling of unbalanced radial distribution systems. IEEE Transactions on Power Systems. 2015;30(6):3102-3109
  18. 18. Stefan A, Jönsson U. Harmonic analysis of pulse-width modulated systems. Automatica. 2009;45(4):851-862
  19. 19. Souvik D, Sahoo SK, Panda SK. Single-phase inverter control techniques for interfacing renewable energy sources with microgrid part I: Parallel-connected inverter topology with active and reactive power flow control along with grid current shaping. IEEE Transactions on Power Electronics. 2011;26(3):717-731
  20. 20. Fatih E. Plug-in repetitive control of single-phase grid-connected inverter for AC module applications. IET Power Electronics. 2017;10(1):47-58
  21. 21. Hitoshi H, Sayama K, Ohishi K, Shimizu T. Current control system based on repetitive control and disturbance observer for single-phase five-level inverter. In: Industrial Electronics Society, IECON. IEEE. pp. 004370-004375;2015
  22. 22. Paolo M, Marafao FP. Repetitive-based control for selective harmonic compensation in active power filters. IEEE Transactions on Industrial Electronics. 2004;51(5):1018-1024
  23. 23. Zhou Xuesong, Daichun Song, Youjie Ma, and Deshu Cheng. Grid-connected control and simulation of single-phase two-level photovoltaic power generation system based on repetitive control. International Conference on Measuring Technology and Mechatronics Automation (ICMTMA). Vol. 2. IEEE; 2010. pp. 366-369
  24. 24. Keliang Z, Yang Y, Blaabjerg F. Frequency adaptive repetitive control of grid-tied single-phase PV inverters. Energy Conversion Congress and Exposition (ECCE), 2015 IEEE; 2015. pp. 1689-1693
  25. 25. Weiss G. Repetitive control systems: Old and new ideas. In: Systems and Control in the twenty-First Century. Boston, MA: Birkhäuser; 1997. pp. 389-404
  26. 26. Sufen C, Lai YM, Tan S-C, Tse CK. Analysis and design of repetitive controller for harmonic elimination in PWM voltage source inverter systems. IET Power Electronics. 2008;1(4):497-506
  27. 27. Gerardo E, Martinez PR, Leyva-Ramos J. Analog circuits to implement repetitive controllers with feedforward for harmonic compensation. IEEE Transactions on Industrial Electronics. 2007;54(1):567-573
  28. 28. Keliang Z, Wang D, Low K-S. Periodic errors elimination in CVCF PWM DC/AC converter systems: Repetitive control approach. IEEE Proceedings-Control Theory and Applications. 2000;147(6):694-700
  29. 29. Francis Bruce A, Wonham WM. The internal model principle for linear multivariable regulators. Applied Mathematics and Optimization. 1975;2(2):170-194
  30. 30. Li C, Dongchun Z, Xianyi Z. Theory and applications of the repetitive control. SICE 2004 Annual Conference. 2004;1:27-34. IEEE
  31. 31. Wu-Sung Y, Tsai M-C. Analysis and estimation of tracking errors of plug-in type repetitive control systems. IEEE Transactions on Automatic Control. 2005;50(8):1190-1195
  32. 32. Meshram RV, Bhagwat M, Khade S, Wagh SR, Stankovi’c AM, Singh NM. Port-controlled Phasor Hamiltonian modeling and IDA-PBC control of solid-state transformer. IEEE Transactions on Control Systems Technology. Nov 2017

Written By

Ragini Meshram, Monika Madhusoodan and Sushama Wagh

Submitted: April 2nd, 2017 Reviewed: January 29th, 2018 Published: November 5th, 2018