This paper proposes a novel direct power control (DPC) strategy for neutralpointclamped (NPC) threelevel rectifiers, to directly control the active power, the reactive power and the neutral point potential of the rectifiers by referring to three precalculated vector influence tables and minimizing an objective function. In the three vector influence tables, the influences of different voltage vectors on the active power, the reactive power and the neutralpoint potential are shown explicitly. A conceptual description and control algorithm of the proposed controller are presented in this paper. Then, numerical simulations and experiments are carried out to validate the proposed method. Both the simulation and experimental results show that good performances during both the steadystate and transient operating conditions are achieved. As a result, the proposed strategy has been proven to be effective for NPC threelevel rectifiers.
I. INTRODUCTION
Neutralpointclamped (NPC) threelevel rectifiers have been widely used in medium power applications since Nabae
et al
. proposed the topology in 1981. This is due to the fact that it has a number of advantages such as low harmonic distortion of the input currents, low switching voltage stress and low d
_{u}
/d
_{t}
[1]

[6]
.
Although the threelevel NPC topology provides significant advantages over the conventional twolevel topology in highpower applications, it still has some inherent problems which have limited its practical application such as the voltage drifts and voltage ripples of the neutralpoint. Therefore, many hardware and software methods have been proposed for balancing the neutralpoint potential.
The highperformance control strategies of threelevel voltage source PWM rectifiers are mainly the voltage oriented control (VOC)
[7]
and direct power control (DPC)
[8]
, which are similar to the vector control (VC) and direct torque control (DTC)
[9]
for AC machines. The VOC guarantees high dynamic and steadystate performances via internal current control loops. This scheme decouples the converter currents into active and reactive power components. The control of the active and reactive powers is then achieved by controlling the decoupled converter currents using current controllers. One main drawback of such a system is that the performance is highly dependent on the applied current control strategy and the connected AC network conditions
[10]
. The DPC strategy is based on the evaluation of active and reactive instantaneous power error values to choose appropriate voltage vectors to control the active and reactive power. When compared with the VOC strategy, DPC does not need internal current control loops or separate PWM modules. It has a simple control algorithm, fast dynamic response and high efficiency. However, since an optimal switching table is essential to the conventional DPC system, most of the switching tables of threelevel rectifiers are established on the basis of qualitative analysis, without providing the quantitative influence of each space voltage vector on the active and reactive powers in each sector. Therefore, the performance of power control with the conventional switching table method is unsatisfactory.
There has been literature lot of research on how to overcome the shortcomings of the conventional DPC strategy. The three main methods are listed as follows:
1) Using multilevel hysteresis instead of twolevel hysteresis. Fourlevel
[11]
and fivelevel
[12]
hysteresis DPC are proposed to increase the flexibility in terms of selecting an appropriate voltage vector, but they are still relatively rough and use qualitative control instead of quantitative control.
2) Using SVM instead of a switching table. DPCSVM strategies usually calculate their output voltage references by means of closedloop schemes with PI controllers
[13]

[15]
, with errors of the active and reactive powers as the input of the controller. The power can be adjusted accurately and smoothly, and a constant switching frequency can be achieved. The main drawback of such a system is that the performance is highly dependent on the tuning of the PI controller.
3) Using predictive direct power control (PDPC). Rodriguez et al. proposed a new strategy that eliminates the hysteresis controllers and switching table
[16]
,
[17]
. In the proposed control strategy, the effects of each switching state on the input power are evaluated, and the switching state that minimizes a quality function is selected and applied during the next sampling period. However, since most PDPC methods require a large amount of online computation in every control cycle, how to reduce the computational burden has become a major issue for PDPC.
In this paper, a novel direct power control (DPC) strategy is proposed for NPC threelevel rectifiers. In order to control the active and reactive power and the neutral point potential, three precalculated vector influence tables are established and an optimized switching state is selected by an objective function. This paper is organized as follows. Section II introduces the threelevel PWM rectifier model and the neutralpoint potential. The influences of different voltage vectors on the active and reactive power and the neutralpoint potential are analyzed in Section III. The vector influence tables and the novel DPC algorithm are proposed in Section IV. Simulation and experimental studies are carried out in Section V and VI. Finally, some conclusions are summarized in Section VII.
II. THREELEVEL VOLTAGESOURCE RECTIFIER
 A. Model of a ThreeLevel PWM Rectifier
A typical NPC threelevel voltagesource rectifier is shown in
Fig. 1
.
e
_{a}
,
e
_{b}
,
e
_{c}
,
i
_{a}
,
i
_{b}
, and
i
_{c}
are the threephase grid voltages and currents, and
i
_{0}
is the neutralpoint current.
L
and
R
represent the ACside line inductance and equivalent resistance, respectively.
C
_{1}
and
C
_{2}
are the capacitances of the DCbus’s upper and lower capacitors, whose voltages are
U
_{dc1}
and
U
_{dc2}
.
Main circuit of a NPC threelevel voltagesource rectifier.
The switching function
S_{j}
(
S_{j}
=1, 0 and 1,
j
=a, b, c) is introduced to describe the states of the four switching devices of each phase. Taking phase a as an example, when T
_{a1}
and T
_{a2}
are on,
S
_{a}
=1, and the voltage between point a to point o is
U
_{dc}
/2. When T
_{a2}
and T
_{a3}
are on,
S
_{a}
=0, and the voltage between point a to point o is 0. When T
_{a3}
and T
_{a4}
are on,
S
_{a}
=1, and the voltage between point a to point o is 
U
_{dc}
/2. As a result, each phase voltage has three levels,
U
_{dc}
/2, 0 and 
U
_{dc}
/2, leading to 3
^{3}
=27 space voltage vectors for a threephase threelevel rectifier, as shown in
Fig. 2
. The space is divided into 12 zones according to the positions of the space voltage vectors. Based on their magnitudes, these vectors can be divided into four groups: zero vectors (
u_{Z}
), small vectors (
u_{S}
), medium vectors (
u_{M}
), and large vectors (
u_{L}
). The magnitudes of the different voltage vectors are listed in
Table I
.
Space voltage vectors of NPC threelevel rectifiers.
VOLTAGE VECTORS OF NPC THREELEVEL RECTIFIERS
VOLTAGE VECTORS OF NPC THREELEVEL RECTIFIERS
According to Kirchhoff's voltage law, the mathematical model in the dq rotating reference frame can be expressed by:
where
e
_{d}
,
e
_{q}
,
i
_{d}
,
i
_{q}
,
u
_{d}
and
u
_{q}
are the grid voltages, AC currents and rectifier input voltages of the PWM rectifier in the dq coordinates, respectively.
ω
is the grid voltage angular velocity.
The expressions of the instantaneous active power
p
and the reactive power
q
can be calculated as follows:
In the dq rotating coordinates, the vector of threephase grid voltages can be expressed as:
where
E
_{l}
is the rms of the grid line voltage.
The expression of the instantaneous active power
p
and the reactive power
q
can be simplified as follows:
By multiplying both sides of equation (1) by
e
_{d}
, the following expressions can be obtained:
By substituting (3) and (4) into (5), the DPC model of a threelevel NPC rectifier can be expressed as:
 B. NeutralPoint Potential of ThreeLevel PWM Rectifiers
In the high frequency mathematical model of threelevel rectifiers
[18]
, the relationship between the DCbus neutralpoint voltage variation and the switching state is obtained as:
where
C
is the sum of
C
_{1}
and
C
_{2}
, and Δ
U
_{dc}
is the difference between the two capacitorvoltage values.
III. INFLUENCE OF SPACE VOLTAGE VECTORS ON ACTIVE AND REACTIVE POWER AND NEUTRALPOINT POTENTIAL
 A. Influence of Space Voltage Vectors on Active and Reactive Power
Each voltage vector of a NPC threelevel rectifier can be decomposed into d and q components. Take the first switching state as an example. Its corresponding voltage vector
u
_{1}
projected on the dq frame is presented in
Fig. 3
(a). According to
Table I
, the magnitude of
u
_{1}
is
Therefore, the dq components of
u
_{1}
can be expressed as:
where
u
_{1}
_{d}
and
u
_{1}
_{q}
are projections on the daxis and qaxis of
u
_{1}
, respectively, and
θ
is the angular position of the grid voltage vector.
(a) Space voltage vector projected on dq frame. (b) AC current vector projected on abc frame.
According to (6), it can be noted that the change rates of the active and reactive power are related to the daxis and qaxis components of the voltage vector, respectively. The influence function of the active power
ξ
, the influence function of the reactive power
μ
, and the parameters
m
_{1}
and
m
_{2}
are defined as:
It can be noted that both
ξ
and
μ
are related to the voltage vector. By substituting (8) into (9) and (10), the influence functions of the active and reactive power of the first switching state are obtained as:
The
ξ
and
μ
functions of all of the voltage vectors of NPC rectifiers are calculated based on the aforementioned contents. The results of the
ξ
functions are shown in
Fig. 4
. Since the values of the
ξ
functions are very small, they are amplified 24 times in
Fig. 4
to make them easier to analyze. The
ξ
functions of the zero voltage vectors are not shown in
Fig. 4
because their values are always 0. It can be seen from
Fig. 4
that the
ξ
functions of the nonzero voltage vectors can be classified into three categories based on the magnitude of the space voltage vector and there are several characteristics of them:
Plot graphs of function ξ when the space voltage vector belongs to (a) Big voltage vectors. (b) Medium voltage vectors. (c) Small voltage vectors.
1) All the
ξ
functions are cosine functions. Therefore they are periodic functions.
2) In each category, the
ξ
functions of the different space voltage vectors have the same magnitude, but different phase angles. For example, in the big voltage vector category,
ξ
_{2}
lags
ξ
_{1}
by π/3, and
ξ
_{3}
lags
ξ
_{2}
by π/3, and so on.
The results of the
_{μ}
functions are very similar to those of the
ξ
functions except that all of the
_{μ}
functions are sine functions instead of cosine functions.
By substituting (9)(12) into (6), the relationships between the influence function
ξ
,
_{μ}
and the change rates of the active and reactive power can be expressed as:
It can be noted from (15) that the relationships are proportional and can be simplified as follows:
 B. Influence of Space Voltage Vectors on NeutralPoint Potential
According to (7), the neutral point current
i
_{0}
can be expressed by:
Since different values of
S_{j}
(
S_{j}
=1, 0 and 1,
j
=a, b, c) correspond to different space voltage vectors, by substituting different values of
S_{j}
into (17), the relationship between the space voltage vector and the neutralpoint current can be obtained, as shown in
Table II
. Notice that only the medium voltage vectors and the small voltage vectors are related to the unbalance of the neutralpoint potential. As a result, only the neutral point currents of the medium and small voltage vectors are shown in
Table II
.
THE RELATIONSHIP BETWEEN RECTIFIER VOLTAGE VECTOR AND NEUTRAL POINT CURRENT
THE RELATIONSHIP BETWEEN RECTIFIER VOLTAGE VECTOR AND NEUTRAL POINT CURRENT
It can be noted from (7) that the change rate of Δ
U
_{dc}
is related to the switching state and the AC currents. The neutralpoint potential influence function
δ
is defined as:
where 
i
 is the amplitude of the AC current vector.
As shown in
Fig. 3
(b), the AC current vector
i
can be expressed in the abc frame as:
where
α
is the angular position of the AC current vector.
According to (17)(19), the
δ
functions can be expressed as:
It can be noticed that all of the
δ
functions are cosine functions, and that they have characteristics similar to those of the
ξ
and
μ
functions.
According to (7), (17) and (18), the relationship between the influence function
δ
and the change rate of Δ
U
_{dc}
can be expressed as:
It can be noted from (21) that the relationship is proportional and can be simplified as follows:
IV. PROPOSED DIRECT POWER CONTROL
 A. Vector Influence Tables of the Active and Reactive Power
In order to show the values of
ξ
and
μ
explicitly, two offline vector influence tables are needed. In these tables, each value of
ξ
or
μ
is not an exact value, but the average value of the small zone it belongs to. A process of discretization of the angular location of the grid voltage vector
θ
is carried out to divide the whole plane evenly into 12 zones. Since
θ
∈[0,2π], each zone has π/6 of the angular scope of
θ
.
The average values of
ξ
and
μ
in each zone are calculated as:
where
n
_{θ}
={1,2,3…12}, round[] denotes the rounding of the number to the nearest integer, and
f
_{ξ}
and
f
_{μ}
are the influence factors of the active power and reactive power, respectively. In order to reduce the rounding errors, a coefficient
k
is added. When
k
is set to 24, the values of
f
_{ξ}
and
f
_{μ}
are integers between 24 and +24.
According to (16), (23) and (24), the approximate proportional relationship between
f
_{ξ}
,
f
_{μ}
and the change rates of the active and reactive power can be expressed as:
where avg( ) denotes the average value in each zone.
By substituting (8) into (23) and (24), the values of
f
_{ξ}
and
f
_{μ}
in each of the 12 zones corresponding to
u
_{1}
are calculated. In a similar way, the values of
f
_{ξ}
and
f
_{μ}
in each zone corresponding to other voltage vectors can also be calculated and the vector influence tables of the active and reactive power are established. As shown in
Table III
and
Table IV
, the values of
f
_{ξ}
and
f
_{μ}
are listed at the bottom of the vector influence tables. The values of
f
_{ξ}
and
f
_{μ}
of the zero voltage vectors are not shown in
Table III
or
Table IV
because they are always 0. The grayscale color of each small block is colored proportional to the value in the block. At the top of the table are the numbers of
n
_{θ}
arranged in different orders in accordance with different space voltage vectors. For example, in the big voltage vector category,
u
_{3}
lags
u
_{1}
by π/3. Therefore, by shifting the numbers of
n
_{θ}
in the first row where
u
_{1}
belongs to the left by two blocks, the numbers of the second row where
u
_{3}
belongs are obtained.
VECTOR INFLUENCE TABLE OF ACTIVE POWER
VECTOR INFLUENCE TABLE OF ACTIVE POWER
VECTOR INFLUENCE TABLE OF REACTIVE POWER
VECTOR INFLUENCE TABLE OF REACTIVE POWER
To obtain the values of
f
_{ξ}
and
f
_{μ}
in the tables, the following steps should be followed. For example, if the value of
f
_{ξ}
of vector
u
_{16}
when
n
_{θ}
=7 is to be found, the specific steps are presented as follows:

Step 1: Find the voltage vectoru16in the lefttop area ofTable III.

Step 2: In the row that vectoru16belongs to, find the number 7 ofnθand the column it belongs to.

Step 3: In that column, find the value in the row of the type of voltage vector at the bottom of the table. In the case of vectoru16, since it is a small voltage vector, the number 8 in the row of the small voltage vector is the value offξ.
 B. Vector Influence Table of the NeutralPoint Potential
The influence factor of the neutralpoint potential
f
_{δ}
is calculated as:
where
n
_{α}
={1,2,3…12}.
According to (22) and (27), the approximate proportional relationship between
f
_{δ}
and the change rate of Δ
U
_{dc}
can be expressed as:
The vector influence table of the neutralpoint potential can be established in a similar way, as shown in
Table V
. Notice that only the medium and small voltage vectors are related to the unbalance of the neutralpoint potential, and only the values of
f
_{δ}
of the medium and small voltage vectors are listed in
Table V
. In other words, the values of
f
_{δ}
of the large and zero voltage vectors are always 0. The steps to obtain the values of
f
_{δ}
in the table are similar to the steps illustrated in Section A.
VECTOR INFLUENCE TABLE OF NEUTRAL POINT POTENTIAL
VECTOR INFLUENCE TABLE OF NEUTRAL POINT POTENTIAL
 C. Proposed DPC Block Diagram and Objective Function
A block diagram of the proposed DPC is presented in
Fig. 5
. It can be seen that three vector influence tables and an objective function optimization part are used instead of the hysteresis comparators and switching table of the conventional method. The reference value of the active power
p
* is generated by a PI controller in the outer voltage loop while the reference value of the reactive power
q
* is set to 0.
Proposed DPC Block diagram.
The values of the instantaneous active power
p
and the reactive power
q
can be calculated as follows [28]:
According to (25), (26) and (28), the reference values of
f
_{ξ}
,
f
_{μ}
and
f
_{δ}
are defined as:
where
K
_{i}
and
f
_{δ}
are parameters which can be calculated by:
where
T
_{s}
is the sampling period.
In order to achieve power tracking and neutral point balance at the same time, the optimized space voltage vector should be selected so that the values of
f
_{ξ}
,
f
_{μ}
and
f
_{δ}
can approach as near as possible to their reference values of
f
_{ξ}
*,
f
_{μ}
* and
f
_{δ}
*. Therefore, the objective function can be described as:
where
λ
_{δ}
is the weighting factor. The objective function
g
is calculated for each set of values of
f
_{ξ}
,,
f
_{μ}
and
f
_{δ}
corresponding to the 27 space voltage vectors. The space voltage vector that corresponds to the minimal value of g will be selected and applied during the next control period. By giving
λ
_{δ}
different values, the priorities of the power control and the neutral point balance control can be made differently.
V. SIMULATION RESULTS
The performances of the conventional twolevel hysteresis DPC method and the proposed DPC method are given and compared by Matlab/Simulink. The simulation parameters are presented in
Table VI
.
SIMULATION AND EXPERIMENTAL PARAMETERS
SIMULATION AND EXPERIMENTAL PARAMETERS
The dynamic performances of the conventional and proposed methods are shown in
Fig. 6
, where an abrupt change of load (from 80Ω to 40Ω) is introduced at
t
=0.4s.
Simulation result of dynamic performance. (a) active and reactive power, neutralpoint potential and AC currents for the conventional method. (b) active and reactive power, neutralpoint potential and AC currents for the proposed method.
It can be seen that both the conventional and the proposed methods have a fast dynamic response. However, the conventional method causes larger power and neutral point potential ripples and a seriously distorted AC current, which are greatly improved by the proposed method. To further illustrate and validate the improvements of the proposed method, the standard deviations of the active and reactive power
σ
_{p}
and
σ
_{q}
and the neutral point potential
σ
_{npp}
are calculated. The standard deviation of a generic variable
x
is expressed as:
where
n
is the number of samples, and
According to the calculated results,
σ
_{p}
and
σ
_{q}
of the conventional method during the time period of 0.55s to 0.6s are 15.22W and 19.85Var, while
σ
_{p}
and
σ
_{q}
of the proposed method are only 4.86W and 6.06Var. The standard deviations of the active and reactive power are reduced by 68.07% and 69.47%, respectively. In addition, the values of
σ
_{npp}
of the conventional method and the proposed method are 292.9mV and 40.5mV, respectively. The ripple size of the neutral point potential is reduced by 86.17%.
In
Fig. 7
, the robustness against system parameter variations is tested for the proposed method. In the simulation, the dq inductances are increased to 15mH (150% of their nominal values), and the other simulation conditions are the same as those in
Fig. 6
. It can be seen that the variation of line inductance does not have a significant effect on system performance with the values of
σ
_{p}
and
σ
_{q}
being 5.02W and 6.16Var.
Simulation result of active and reactive power, neutralpoint potential and AC currents for the proposed method with 1.5L.
Fig. 8
shows the standard deviations of the active and reactive power, the total harmonic distortion (THD) of the AC currents and the standard deviation of the neutral point potential of the proposed method with different
K
_{i}
and
K
_{δ}
. According to (31) and the parameters listed in
Table VI
,
E
_{1}
is estimated to be 69V, 
i
 is estimated to be 4A, and the theoretical values of
K
_{i}
and
K
_{δ}
can be calculated as 1.4W and 1.5mV, respectively. According to (11) and (12),
m
_{1}
and
m
_{2}
are estimated to be 0.70 and 0.16, respectively. It can be seen from
Fig. 8
(a) that when
K
_{i}
changes within a range between 0.8W and 2W, the values of
σ
_{p}
and
σ
_{q}
remain within a range between 4W/Var and 8W/Var, and the THDs of the AC currents remain within a range between 1.2% and 2.1%. It can also be seen that both
σ
_{p}
and
σ
_{q}
and the THD of the AC currents are minimal when
K
_{i}
is approximately 1.4W, which is consistent with the theoretical value of
K
_{i}
. It can be seen from
Fig. 8
(b) that when
K
_{δ}
changes within a range between 0.6mV and 2.7mV, the values of
σ
_{npp}
are within a range between 35mV and 55mV, and
σ
_{npp}
is minimal when
K
_{δ}
is approximately 1.5mV, which is consistent with the theoretical value of
K
_{δ}
.
Simulation results with the variation of K_{i}、K_{δ}. (a) σ_{p}、σ_{q} and THD of i_{A} with the variation of K_{i}. (b) σ_{npp} with the variation of K_{δ}.
VI. EXPERIMENTAL RESULTS
To verify the feasibility and effectiveness of the proposed method, a threelevel rectifier prototype was implemented, as shown in
Fig. 9
. The system consists of a main circuit, a control circuit, a driver circuit and a sampling circuit. The complete system was connected to the utility grid through a threephase variac. The control system employed a TMS320F28335 digital signal processor for the control strategy. The sampling periods of the conventional and proposed methods are both 50μs. The parameters of the experimental system are presented in
Table VI
. According to section V,
K
_{i}
and
K
_{δ}
are chosen as 1.4W and 1.5mV, respectively, and
m
_{1}
and
m
_{2}
are chosen to be 0.70 and 0.16, respectively.
Photograph of experimental setup.
 A. SteadyState Performance
The steadystate performances of the two methods are presented in
Fig. 10
with the load being 40Ω and the ACside inductance being 10mH. It is clearly seen that the conventional DPC method results in significant current harmonic distortions as well as power and neutral point potential ripples, while the proposed DPC method has achieved a greatly improved performance with smaller ripples of the active and reactive powers, a wellsmoothed neutral point potential and more sinusoidal AC currents. To be specific, with the conventional method, as shown in
Fig. 10
(a),
σ
_{p}
and
σ
_{q}
are 30.37W and 26.78Var, while with the proposed method, as shown in
Fig. 10
(b),
σ
_{p}
and
σ
_{q}
are only 15.17W and 20.35Var. The percentages of the standard deviation reduction of the active and reactive power are 50.05% and 24.01%, respectively. The THD of the AC currents is further calculated. The result shows that the THDs of the AC currents of the conventional and proposed methods are 5.92% and 1.40%, respectively. In addition, the ripple size of the neutral point potential is reduced by 69.60% since the values of
σ
_{npp}
of the conventional method and proposed method are 0.556V and 0.169V, respectively.
Experimental waveforms for steadystate performance of active and reactive power, neutral point potential and AC currents with the load being 40Ω. (a) the conventional method. (b) the proposed method.
 B. Dynamic Performance
The dynamic performances of the two methods are presented in
Fig. 11
with load stepping from 80Ω to 40Ω. Once again, the improvements are evident in the case of the proposed DPC method. It can be seen that with an abrupt change of load, for both the conventional and the proposed methods, the active and reactive power and the neutral point potential respond rapidly and meet the reference values only after a short period. However, remarkable ripple reductions of the active and reactive power and the neutral point potential can be observed in the proposed method.
Experimental waveforms of active and reactive power, neutral point potential and AC currents for step change of load from 80Ω to 40Ω. (a) the conventional method. (b) the proposed method.
 C. Impact of the Variations of the Line Inductance
Further tests on the impact of line inductance variations on the steadystate performance with the proposed DPC strategy are carried out.
Fig. 12
shows the waveforms of the active and reactive power, the neutral point potential and the AC currents of the proposed method when there are errors in the ACside inductance parameter. As shown in
Fig. 12
(a), with the line inductance
L
being 8mH while the other experimental conditions are the same as those in
Fig. 10
(b),
σ
_{p}
and
σ
_{q}
of the proposed method become 25.81W and 23.35Var, and the THD of the AC currents becomes 1.63%. As shown in
Fig. 12
(b), with the line inductance
L
being 12mH while the other experimental conditions are the same as those in
Fig. 10
(b),
σ
_{p}
and
σ
_{q}
of the proposed method become 25.28W and 24.56Var, and the THD of the AC currents becomes 1.60%. In both cases, the ripple sizes of the neutral point potential remain small. It can be seen that such line inductance errors have no significant influence on system performance.
Experimental waveforms of active and reactive power, neutral point potential and AC currents for the proposed method with variations of the line inductance. (a) L=8mH. (b) L=12mH
VII. CONCLUSIONS
In this paper, a novel DPC algorithm has been proposed for NPC threelevel rectifiers. This paper defines the influence functions to analyze the relationships between the voltage vectors of the NPC rectifiers and the change rates of the active power, the reactive power and the neutral point potential. By discretizing and averaging the influence functions, the influence factors of the active power, the reactive power and the neutral point potential
f
_{ξ}
,
f
_{μ}
and
f
_{δ}
are obtained, and they have approximately proportional relationships with the change rates of the active power, the reactive power and the neutral point potential, respectively. Furthermore, three precalculated vector influence tables are established based on the influence factors, and by minimizing an objective function an appropriate vector can be selected. The performance of the novel DPC has been verified by simulation and experimental results. Compared to the conventional DPC method, remarkable ripple reductions of the active and reactive power and the neutral point potential can be achieved by the proposed method.
Acknowledgements
The authors are grateful to the National Key Basic Research Program of China (2013CB035602), and the Research Fund for the Doctoral Program of Higher Education of China (20110032120082) for supporting this work.
BIO
ChangLiang Xia was born in Tianjin, China, in 1968. He received his B.S. degree from Tianjin University, Tianjin, China, in 1990, and his M.S. and Ph.D. degrees from Zhejiang University, Hangzhou, China, in 1993 and 1995, respectively, all in Electrical Engineering. He is currently a Professor in the School of Electrical Engineering and Automation, Tianjin University, and also in the Tianjin Key Laboratory of Advanced Technology of Electrical Engineering and Energy, Tianjin Polytechnic University, Tianjin, China. In 2008, he became a “Yangtze Fund Scholar” Distinguished Professor and is currently supported by the National Science Fund for Distinguished Young Scholars. His current research interests include electrical machines and their control systems, power electronics, and the control of wind generators.
Zhe Xu was born in Handan, Hebei, China, in 1988. He received his B.S. and M.S. degrees in Electrical Engineering from Tianjin University, Tianjin, China, in 2011, 2014, respectively. His current research interests include control systems and power electronics.
JiaXin Zhao was born in Tianjin, China, in 1984. He received his B.S. degree from the China Agricultural University, Beijing, China, in 2006, and his M.S. and Ph.D. degrees from Tianjin University, Tianjin, China, in 2010 and 2014, respectively, all in Electrical Engineering. He is currently an Embedded Software Engineer for REnergy Electric Tianjin Ltd., Tianjin, China. His current research interests include electrical machines and motor drives, power electronics and wind power technology.
Xia C. L.
,
Gu X.
,
Shi T. N.
,
Yan Y.
2011
“Neutralpoint potential balancing of threelevel inverters in directdriven wind energy conversion system,”
IEEE Trans. Energy Conver.
26
(1)
18 
29
DOI : 10.1109/TEC.2010.2060487
Lee D. M.
,
Jung J. W.
,
Kwak S. S.
2011
“Simple space vector PWM scheme for 3level NPC inverters including the overmodulation region,”
Journal of Power Electronics
11
(5)
688 
696
DOI : 10.6113/JPE.2011.11.5.688
Ye Z. B.
,
Xu Y. M.
,
Li F.
,
Deng X. M.
,
Zhang Y. Z.
2014
“Simplified PWM strategy for neutralpointclamped (NPC) threelevel converter,”
Journal of Power Electronics
14
(3)
519 
530
DOI : 10.6113/JPE.2014.14.3.519
Li N.
,
Wang Y.
,
Lei W. J.
,
Niu R. G.
,
Wang Z. A.
2014
“Novel carrierbased pwm strategy of a threelevel NPC voltage source converter without lowfrequency voltage oscillation in the neutral point,”
Journal of Power Electronics
14
(3)
531 
540
DOI : 10.6113/JPE.2014.14.3.531
Zaragoza J.
,
Pou J.
,
Ceballos S.
,
Robles E.
,
Ibaez P.
,
Villate J. L.
2009
“A comprehensive study of a hybrid modulation technique for the neutralpointclamped converter,”
IEEE Trans. Ind. Electron.
56
(2)
294 
304
DOI : 10.1109/TIE.2008.2005132
Zhang Y. C.
,
Zhao Z. M.
,
Zhu J. G.
2011
“Hybrid PWM applied to highpower threelevel inverterfed inductionmotor drives,”
IEEE Trans. Ind. Electron.
58
(8)
3409 
3420
DOI : 10.1109/TIE.2010.2090836
Song W. S.
,
Feng X. Y.
,
Xiong C. L.
“a neutral point voltage regulation method with SVPWM control for singlephase threelevel NPC converters,”
IEEE Vehicle Power and Propulsion Conference
2008
1 
4
Li N.
,
Wang Y.
,
Li S.
,
Li Y. F.
,
Wang Z. A.
“Direct power control strategy used in threelevel NPC converters”
International Power Electronics and Motion Control Conference
2012
1675 
1679
Ortega C.
,
Arias A.
,
Caruana C.
,
Balcells J.
,
Asher G.
2010
“Improved waveform quality in the direct torque control of matrixconverterfed PMSM drives,”
IEEE Trans. Ind. Electron.
57
(6)
2101 
2110
DOI : 10.1109/TIE.2009.2033084
Kazmierkowski M. P.
,
Malesani L.
1998
“Current control techniques for threephase voltagesource PWM converters: A survey,”
IEEE Trans. Ind. Electron.
45
(5)
691 
703
DOI : 10.1109/41.720325
Chen W.
,
Zou Y. P.
,
Xu L. J.
“Direct power control for neutralpointclamped threelevel PWM rectifier,”
Proceedings of the IEEE International Conference on Industrial Technology
2008
1 
6
Serpa L. A.
,
Barbosa P. M.
,
Steimer P. K.
,
Kolar J. W.
“Fivelevel virtualflux direct power control for the active neutralpoint clamped multilevel inverter,”
IEEE Annual Power Electronics Specialists Conference
2008
1668 
1764
Ejlali A.
,
Arab Khaburi D.
“Power quality improvement using nonlinearload compensation capability of variable speed DFIG based on DPCSVM method,”
Power lectronics, Drive Systems and Technologies Conference
2014
280 
284
Xiao L.
,
Huang S. D.
,
Liu T.
,
Xu Q.
,
Huang K. Y.
“A novel SVMDPC control method for grid connected ACDC converters under asymmetrical fault,”
International Conference on Power Engineering, Energy and Electrical Drives
2011
1 
6
Xiao L.
,
Huang S. D.
,
Zhen L.
,
Xu Q.
,
Huang K. Y.
“Sliding mode SVMDPC for gridside converter of DPMSG under asymmetrical fault,”
International Conference on Electrical Machines and Systems
2011
1 
6
Rodriguez J.
,
Pontt J.
,
Correa J.
,
Lezana P.
,
Cortes P.
“Predictive power control of an AC/DC/AC converter,”
In Proceedings of IEEE Industry Applications Conference, 40th IAS Annual Meeting
2005
934 
939
Cortes P.
,
Rodriguez J.
,
Antoniewicz P.
,
Kazmierkowski M.
2008
“Direct power control of an AFE using predictive control,”
IEEE Trans Power Electron.
23
(5)
2516 
2523
DOI : 10.1109/TPEL.2008.2002065
Wu W. J.
,
Zhong Y. R.
,
Duan B.
2009
“A symmetrical threesection PWM for threelevel rectifier with medium vectors being eliminated,”
Transactions of China Electrotechnical Society
26
(4)
69 
76+94