Model Predictive Current Control Method with Improved Performances for Three-Phase Voltage Source Inverters

In this paper, the model predictive current control (MPCC) method using two vectors has been proposed to control output currents of three-phase voltage source inverters (VSIs) with small current errors and current ripples. Also, the proposed method can reduce switching losses by applying the vector pre-selection technique to the MPCC for the VSI. The VSI generates seven voltage vectors to control the output currents, but the proposed method uses four available voltage vectors with one switch, which are classified by the vector pre-selection method clamping one leg and conducting the largest output current among the three legs to reduce the switching losses. In the proposed method, selecting two future voltage vectors among the four voltage vectors and dividing them in a future sampling period are determined by an optimization process. The proposed method results in the lower total loss, better total harmonic distortion (THD), and smaller current errors than the conventional method with half the sampling period of the proposed method due to the optimal process. Simulation and experimental results of the three-phase VSIs are presented in order to verify the effectiveness of the proposed method.


Introduction
In the last few decades, the control schemes for voltage source inverters (VSIs) have been presented and widely studied [1,2].The typical control schemes for voltage source inverter are hysteresis control and linear control with pulse width modulation (PWM). New control schemes, such as fuzzy logic control, neural networks, and predictive control, have been proposed due to the improvement of fast and powerful microprocessors. Among the new control schemes, model predictive control (MPC) using a finite-control-set is simple and effective for VSIs because it is simple and flexible without any PWM blocks [3][4][5]. Thanks to its simplicity and control flexibility, the MPC has been used in power electronics for VSIs, including active front ends [6][7][8], multi-level inverters [9][10][11][12], and matrix converters [13][14][15]. The VSIs can be controlled by using the model predictive control, which predicts the finite number of the future output current behaviors based on the load dynamic model. A cost function with respect to the error term between the predicted current and reference current is used to evaluate all predicted current values obtained by possible switching states to select one optimal switching state with the smallest cost value [16][17][18][19]. The optimal switching state is applied to the VSIs in one sampling period [20,21]. The conception of using two vectors has also recently been applied in the MPC for the torque control and flux control to reduce torque and flux ripples and the PWM rectifier to reduce power ripples [22]. This paper proposes a switching strategy, that is, the vector pre-selection technique to obtain high efficiency [23], and the model predictive current control method based on period to reduce current errors and ripples [24,25]. In order to obtain high efficiency, the switching strategy classifies the four available voltage vectors with one switch that clamps one leg carrying the largest current among the three legs in the VSIs. The proposed method selects two vectors among four available vectors classified by pre-selection method by using a cost function defined at the changing point of the two voltage vectors, as well as at the next sampling instants to evaluate the predicted current values. Therefore, the newly defined cost function takes into account current errors at two points, and the proposed method distributes two selected vectors within the sampling period in a manner that minimizes the squared current errors in the cost function. The proposed method leads to lower total loss [26], better current total harmonic distortion (THD), and smaller current errors compared to the model predictive current control (MPCC) method using one vector during one sampling period. Simulation and experimental results of the three-phase VSIs are presented in order to verify the effectiveness of the proposed method.

Conventional Model Predictive Current Control Method for VSIs
The basic structure of the three-phase VSI is shown in the following Figure 1. The three-phase VSI has three legs and each of which consists of two switches. The two switches in each leg operate complementary to each other. The output voltage vectors generated by the switches can be expressed as where , , and are the phase-to-neutral voltages. The three-phase VSI composed of six switches has a total of eight load voltage vectors divided into six active vectors and two zero vectors, as shown in Figure 2, in accordance with a combination of switching states. Since the two zero vectors among eight voltage vectors produce the same output voltage vectors, only seven voltage vectors are considered as control elements. The output current can be also expressed as  The three-phase VSI has three legs and each of which consists of two switches. The two switches in each leg operate complementary to each other. The output voltage vectors generated by the switches can be expressed as v = 2 3 v an + v bn e j(2π/3) + v cn e j(4π/3) where v an , v bn , and v cn are the phase-to-neutral voltages. The three-phase VSI composed of six switches has a total of eight load voltage vectors divided into six active vectors and two zero vectors, as shown in Figure 2, in accordance with a combination of switching states. Since the two zero vectors among eight voltage vectors produce the same output voltage vectors, only seven voltage vectors are considered as control elements. The output current can be also expressed as where i a , i b, and i c are the three-phase output currents. The output current dynamics can be described by the differential equation as where v, i, R, L, and e are the voltage vectors produced by the VSI, the output current vector, the load resistance, inductance, and load back-emf vector, respectively. The output current derivative di/dt in (3) can be replaced by the Forward Euler approximation with a sampling period T s as of switching states. Since the two zero vectors among eight voltage vectors produce the same output voltage vectors, only seven voltage vectors are considered as control elements. The output current can be also expressed as  Then, the output current at the next sampling period can be expressed as Seven future output currents generated by the seven voltage vectors v k in the k th sampling period can be predicted by using (5). A predefined cost function composed of errors between the future output currents and the future reference currents is used to select the optimal vector among seven available voltage vectors minimizing the errors at every sampling period. The two-step future output current at the (k + 2) th instant is obtained by shifting (5) one step forward as in [2] to prevent the unavoidable calculation delay, Thus, the cost function considering calculation delay can be defined as the second order errors as where i and i * are the output current and the reference output current in the αβ frame, respectively. The future reference current value required in the cost function can be obtained by the Lagrange extrapolation formula as The back-emf vector in (6) can be estimated by assuming that the future back-emf vector is equal to the present back-emf vector because the back-emf vector varies at a much lower frequency compared with the fast sampling frequency. The present back-emf vector can be calculated from (5) as Only one voltage vector selected from the seven possible voltage vectors is applied for one full sampling period. Thus, the conventional method uses the optimal voltage vector that minimizes the current error, and the vector is maintained for at least one sampling period.

Proposed Model Predictive Current Control Method to Improve Performance in Terms of Efficiency and Current Harmonics
The instantaneous switching losses of the VSI depend on the amplitude of the output currents through the switches at the moment of switching. Thus, the vector pre-selection method sorts the voltage vectors clamping one leg and carrying the largest current among the three legs to reduce the switching losses. In order to sort the voltage vectors among the seven voltage vectors of the VSI, In the inverse model equation, the output currents are replaced with the reference currents to reduce ripples of the future reference voltage vector. The vector pre-selection method selects the voltage vectors in accordance with the state of the future output currents and phase voltages at every sampling period. The three-phase future voltages obtained by using (10) are classified as v k+1 max , v k+1 mid , and v k+1 min according to the magnitude of each phase future voltage and are assigned as The prohibitive phase that should not be allowed to be clamped corresponds with the phase with the future voltage v k+1 mid , that is, the medium value among the three-phase future voltage to assure linear modulation ranges. In the two voltages classified as v k+1 max and v k+1 min , the future reference currents determine that one phase with a higher absolute output current is clamped to the positive or negative DC-link. Once it is determined that one leg of the VSIs is clamped, the upper switch or the lower switch is compelled to turn on according to the sign of the corresponding reference voltage. Figure 3 shows pre-selected vectors and the clamped switches according to the output currents and phase voltages. Four voltage vectors among seven voltage vectors can clamp one phase exposed to the higher output current between two voltages classified as v k+1 max and v k+1 min , so the four available voltage vectors can be pre-selected. In the inverse model equation, the output currents are replaced with the reference currents to reduce ripples of the future reference voltage vector. The vector pre-selection method selects the voltage vectors in accordance with the state of the future output currents and phase voltages at every sampling period. The three-phase future voltages obtained by using (10) are classified as , , and according to the magnitude of each phase future voltage and are assigned as The prohibitive phase that should not be allowed to be clamped corresponds with the phase with the future voltage , that is, the medium value among the three-phase future voltage to assure linear modulation ranges. In the two voltages classified as and , the future reference currents determine that one phase with a higher absolute output current is clamped to the positive or negative DC-link. Once it is determined that one leg of the VSIs is clamped, the upper switch or the lower switch is compelled to turn on according to the sign of the corresponding reference voltage. Figure 3 shows pre-selected vectors and the clamped switches according to the output currents and phase voltages. Four voltage vectors among seven voltage vectors can clamp one phase exposed to the higher output current between two voltages classified as and , so the four available voltage vectors can be pre-selected.    Table 1 shows the clamped switches, pre-selected vectors, and the conditions of the output currents and phase voltages, making switches clamping. The proposed method uses two vectors among the four available voltage vectors for reducing the output current errors and ripples. Optimal two voltage vectors and their optimal durations have to be selected to generate the output current closely tracking the reference current because of using the two vectors in one sampling period. In this paper, the proposed method using two vectors is proposed to reduce switching losses, current errors, and current ripples. Table 1. Conditions of the vector pre-selection.

Conditions Pre-Selected Vectors Clamped Switch Phase Voltages
Output Currents In order to reduce output current errors and ripples, the proposed method uses two vectors defined as the initial and the next voltage vector at every sampling point. The two vectors are divided by the optimization process in a sampling period. The selected voltage vectors in the proposed method are applied for the calculation of the cost function minimizing the output current errors at every sampling point. One sampling period T s is divided into the time intervals T k 1 and T k 2 depending on selected voltage vectors at every sampling point, as Each of the durations T k 1 and T k 2 has a positive value from 0 to T s because T k 1 and T k 2 are the time intervals for the initial voltage vector v k 1 and the next voltage vector v k 2 , respectively. In the proposed method, the predicted output current in (5) can be modified as In the one sampling period, the back-emf vector e k + T k 1 does not change considerably, so it can be estimated as The present back-emf is assumed by using an extrapolation of the past value by shifting (10) one-step backward as To prevent the control delay, the two-step future output current can be obtained by shifting (13) one step forward as In case that the proposed scheme uses the cost function (7), which considers the current errors at next sampling instant, the current errors at changing point of the two vectors could be higher than at the next sampling instant. So, a cost function of the proposed scheme has to contain the output currents at changing point of the two vectors. Therefore, the cost function for the proposed scheme is considered at two instants, as The reference output currents at the changing point t = (k + 1) + T k+1 1 of the two vectors can be obtained as Also, the actual output currents at the same point t = (k + 1) + T k+1 1 can be obtained by The cost function G considering two points of the current errors depends on T k+1 1 , so the optimal duration T k+1 1 for the vector v k+1 1 can be obtained by differentiating the cost function G with respect to T k+1 By reflecting (8), (16), (17), and (18), (19) to (20), the optimal duration T k+1 1 is determined as The duration T k+1 2 for the vector v k+1 2 during one sampling period can be determined as 11. Determining one optimal vector sets with and along with their optimal durations and .
12. Storing , , , and for the next application at the (k + 1) th instant. The proposed and conventional methods are compared in terms of current errors, THD values, the total losses, and the number of switchings. The current errors and the THD percentage obtained from the proposed method are shown in Figure 5 a, b, according to the switching frequency because the performance of the VSIs is very much affected by the switching frequency due to its operating principle. The current error between the reference currents and the output currents is decided as

One-step Lagrange extrapolation
Calculate future reference current at t = ((k+1)+ ) in (18) Calculate future current at t = ((k+1)+ ) in (19) Calculate and in (21) and (22) T Predict twostep future load current in (16) VSI R L e Select optimal vector set using cost function in (17) Predict onestep future load current in (13) Vector preselection algorithm The proposed and conventional methods are compared in terms of current errors, THD values, the total losses, and the number of switchings. The current errors and the THD percentage obtained from the proposed method are shown in Figure 5 a,b, according to the switching frequency because the performance of the VSIs is very much affected by the switching frequency due to its operating principle. The current error between the reference currents and the output currents is decided as where the value of N is set to 10,000, and the error is the average of 10,000 calculated values. As shown in Figure 5a, the proposed method shows the same results in the current error compared to the conventional method. The THD percentage of the output currents is defined as where i x1 ···i xn are the fundamental and n th -harmonic components of the output currents in phase x, respectively. The simulation was set up to the 8335 th -harmonic components. In Figure 5b, the THD of the proposed method is similar to the conventional method. Power losses of the VSI are calculated to evaluate the switching losses' reduced performance of the proposed method in the simulation. Power losses can be divided into conduction losses and switching losses. The conduction losses P cond in an Insulated Gate Bipolar Transistor (IGBT) can be calculated using the following expression: where V CE , I C , and T s are the collector-emitter voltage, collector current, and sampling period, respectively. In order to determine the constants in (25), a manufacturer's data sheet is used to apply exponential curve fitting on the on-state loss characteristics [8]. Similarly to the conduction losses, the switching losses P swit can be calculated using the following expression: where E on and E o f f are turn-on switching loss and turn-off switching loss, respectively. Turn-on switching loss and turn-off switching loss can be obtained by using the loss model with the switching energies obtained from the datasheet of the IGBT module (SKM50GB123D). As shown in Figure 5c, the total losses of the proposed method are lower than those of the conventional method under different switching frequency due to pre-selecting vectors clamping one leg and conducting the largest output current among the three legs. Figure 5d shows the average number of switchings. In Figure 5d, as the switching frequency increases, the number of switchings of the conventional method and the proposed method increases equally. where the value of N is set to 10,000, and the error is the average of 10,000 calculated values. As shown in Figure 5a, the proposed method shows the same results in the current error compared to the conventional method. The THD percentage of the output currents is defined as where ix1···ixn are the fundamental and n th -harmonic components of the output currents in phase x, respectively. The simulation was set up to the 8335 th -harmonic components. In Figure 5b, the THD of the proposed method is similar to the conventional method. Power losses of the VSI are calculated to evaluate the switching losses' reduced performance of the proposed method in the simulation. Power losses can be divided into conduction losses and switching losses. The conduction losses Pcond in an Insulated Gate Bipolar Transistor (IGBT) can be calculated using the following expression: where , , and are the collector-emitter voltage, collector current, and sampling period, respectively. In order to determine the constants in (25), a manufacturer's data sheet is used to apply exponential curve fitting on the on-state loss characteristics [8]. Similarly to the conduction losses, the switching losses Pswit can be calculated using the following expression: where and are turn-on switching loss and turn-off switching loss, respectively. Turn-on switching loss and turn-off switching loss can be obtained by using the loss model with the switching energies obtained from the datasheet of the IGBT module (SKM50GB123D). As shown in Figure 5c, the total losses of the proposed method are lower than those of the conventional method under different switching frequency due to pre-selecting vectors clamping one leg and conducting the largest output current among the three legs. Figure 5d shows the average number of switchings. In Figure 5d, as the switching frequency increases, the number of switchings of the conventional method and the proposed method increases equally.

4.1.Simulation and Experimental Results
The simulations were done for the following parameters Vdc = 260 V, the RLe load with R = 0.8 Ω, L = 12 mH, e = 20 V, and the amplitude of the reference current was 12 A. Figure 6 shows the output currents and switching signal S1 obtained from the conventional method using one vector with Ts = 250 μs, the conventional method using one vector with Ts = 125 μs, the conventional method using two vectors with Ts = 250 μs [19], and the proposed method using two vectors with vector pre-selection method with Ts = 250 μs. Since the proposed method uses two vectors, the sampling period of the conventional method is half of the sampling period of the proposed method in order to compare the two methods. In Figure 6, the output current ripples of the proposed scheme are lower than those of the conventional methods, and the switching signal S1, which is the upper switch of a-leg, stops operation around the peak output current. Also, the frequency spectrum in Figure 6 shows that the proposed scheme has lower THD values than conventional methods.

Simulation and Experimental Results
The simulations were done for the following parameters V dc = 260 V, the RLe load with R = 0.8 Ω, L = 12 mH, e = 20 V, and the amplitude of the reference current was 12 A. Figure 6 shows the output currents and switching signal S 1 obtained from the conventional method using one vector with T s = 250 µs, the conventional method using one vector with T s = 125 µs, the conventional method using two vectors with T s = 250 µs [19], and the proposed method using two vectors with vector pre-selection method with T s = 250 µs. Since the proposed method uses two vectors, the sampling period of the conventional method is half of the sampling period of the proposed method in order to compare the two methods. In Figure 6, the output current ripples of the proposed scheme are lower than those of the conventional methods, and the switching signal S 1 , which is the upper switch of a-leg, stops operation around the peak output current. Also, the frequency spectrum in Figure 6 shows that the proposed scheme has lower THD values than conventional methods.

4.1.Simulation and Experimental Results
The simulations were done for the following parameters Vdc = 260 V, the RLe load with R = 0.8 Ω, L = 12 mH, e = 20 V, and the amplitude of the reference current was 12 A. Figure 6 shows the output currents and switching signal S1 obtained from the conventional method using one vector with Ts = 250 μs, the conventional method using one vector with Ts = 125 μs, the conventional method using two vectors with Ts = 250 μs [19], and the proposed method using two vectors with vector pre-selection method with Ts = 250 μs. Since the proposed method uses two vectors, the sampling period of the conventional method is half of the sampling period of the proposed method in order to compare the two methods. In Figure 6, the output current ripples of the proposed scheme are lower than those of the conventional methods, and the switching signal S1, which is the upper switch of a-leg, stops operation around the peak output current. Also, the frequency spectrum in Figure 6 shows that the proposed scheme has lower THD values than conventional methods.   Figure 7 shows the output currents and switching signal S 1 obtained from the conventional method using one vector, the conventional method using two vectors [19], and the proposed method when the switching frequency f switching = 4 kHz instead of the sampling frequency. As shown in Figure 7, when operating at the same switching frequency, there was no significant difference in high-frequency components between the proposed method and the conventional methods. Also, it can be seen that the clamping was performed well in the proposed method.
Electronics 2019, 8,625 12 of 20 Figure 7 shows the output currents and switching signal S1 obtained from the conventional method using one vector, the conventional method using two vectors [19], and the proposed method when the switching frequency = 4 kHz instead of the sampling frequency. As shown in Figure 7, when operating at the same switching frequency, there was no significant difference in high-frequency components between the proposed method and the conventional methods. Also, it can be seen that the clamping was performed well in the proposed method. The dynamic responses for the conventional method using one vector with T s = 125 µs, the conventional method using two vectors with T s = 250 µs, and the proposed method with the T s = 250 µs are shown in Figure 8. The amplitude of the reference currents was changed from 6 A to 12 A, and the fundamental frequency of the reference currents was changed from 60 Hz to 90 Hz. In both cases, the output currents of the proposed method had a fast dynamic response like those of the conventional methods.  , and ic), a-leg upper switch S1, and the frequency spectrum of the a-phase current (ia) with |I * | = 12 A. (a) the conventional method using one vector, (b) the conventional method using two vectors, and (c) the proposed method with switching frequency = 4kHz. THD: total harmonic distortion.
The dynamic responses for the conventional method using one vector with Ts = 125 μs, the conventional method using two vectors with Ts = 250 μs, and the proposed method with the Ts = 250 μs are shown in Figure 8. The amplitude of the reference currents was changed from 6 A to 12 A, and the fundamental frequency of the reference currents was changed from 60 Hz to 90 Hz. In both cases, the output currents of the proposed method had a fast dynamic response like those of the conventional methods.  Figure 9 shows the dynamic responses for the conventional method using one vector, the conventional method using two vectors, and the proposed method with switching frequency = 4 kHz. The amplitude of the reference currents was changed from 6 A to 12 A, and the fundamental frequency of the reference currents was changed from 60 Hz to 90 Hz. In both cases, the output currents of the proposed method had a fast dynamic response like those of the  Figure 9 shows the dynamic responses for the conventional method using one vector, the conventional method using two vectors, and the proposed method with switching frequency f switching = 4 kHz. The amplitude of the reference currents was changed from 6 A to 12 A, and the fundamental frequency of the reference currents was changed from 60 Hz to 90 Hz. In both cases, the output currents of the proposed method had a fast dynamic response like those of the conventional methods.
(c) Figure 8. Simulation results of the three-phase output currents (ia, ib, and ic) and the a-leg upper switch S1 for the current magnitude step change from 6 A to 12 A and frequency step change from 60 Hz to 90 Hz (|I*| = 12 A). (a) the conventional method using one vector with Ts = 125 μs, (b) the conventional method using two vectors with Ts = 250 μs, and (c) the proposed method with Ts = 250 μs. Figure 9 shows the dynamic responses for the conventional method using one vector, the conventional method using two vectors, and the proposed method with switching frequency = 4 kHz. The amplitude of the reference currents was changed from 6 A to 12 A, and the fundamental frequency of the reference currents was changed from 60 Hz to 90 Hz. In both cases, the output currents of the proposed method had a fast dynamic response like those of the conventional methods. (c) Figure 9. Simulation results of the three-phase output currents (ia, ib, and ic) and the a-leg upper switch S1 for the current magnitude step change from 6 A to 12 A and frequency step change from 60 Hz to 90 Hz (|I*| = 12 A). (a) the conventional method using one vector, (b) the conventional method using two vectors, and (c) the proposed method with switching frequency = 4 kHz.
The MPCC method using two vectors with pre-selection method was implemented by configuring a prototype setup. The experimental setup was composed of the three-phase VSI with the IGBT module and RL load. A digital signal processor (DSP, TMS320F28335) was used to execute the switching algorithm for the VSI. Figure 10 shows the experimental waveforms of the output Figure 9. Simulation results of the three-phase output currents (i a , i b , and i c ) and the a-leg upper switch S 1 for the current magnitude step change from 6 A to 12 A and frequency step change from 60 Hz to 90 Hz (|I*| = 12 A). (a) the conventional method using one vector, (b) the conventional method using two vectors, and (c) the proposed method with switching frequency f switching = 4 kHz.
The MPCC method using two vectors with pre-selection method was implemented by configuring a prototype setup. The experimental setup was composed of the three-phase VSI with the IGBT module and RL load. A digital signal processor (DSP, TMS320F28335) was used to execute the switching algorithm for the VSI. Figure 10 shows the experimental waveforms of the output currents and the switching signal of the switch S 1 , which is the upper switch of a-leg. The proposed method was performed under the conditions with T s = 250 µs, V dc = 260 V, and the amplitude of the reference current |I * |= 12 A. In Figure 10, the output current ripples of the proposed method with T s = 250 µs are less than those of the conventional method with T s = 125 µs. Besides, the experimental results obtained from the proposed method and the conventional method were comparable to the simulation results, as shown in Figure 6. The switching signal S 1 stopped operation around the peak output current because of the vector pre-selection method. In Figure 10, the THD values were obtained from the digital oscilloscope MSO3054 set up to the 400 th harmonic components, which can be considered to be the maximum number of harmonics. The dynamic responses of the proposed method with the sampling period Ts = 250 μs and the conventional method with Ts = 250 μs are shown in Figure 11. The amplitude of the reference currents was changed from 6 A to 12 A, and the fundamental frequency of the reference currents was changed from 60 Hz to 90 Hz. In the step changes of magnitude and fundamental frequency of the reference currents, the output currents of the proposed method tracked the reference change as fast as those of the conventional method with half the sampling period.
(a) The dynamic responses of the proposed method with the sampling period T s = 250 µs and the conventional method with T s = 250 µs are shown in Figure 11. The amplitude of the reference currents was changed from 6 A to 12 A, and the fundamental frequency of the reference currents was changed from 60 Hz to 90 Hz. In the step changes of magnitude and fundamental frequency of the reference currents, the output currents of the proposed method tracked the reference change as fast as those of the conventional method with half the sampling period. The dynamic responses of the proposed method with the sampling period Ts = 250 μs and the conventional method with Ts = 250 μs are shown in Figure 11. The amplitude of the reference currents was changed from 6 A to 12 A, and the fundamental frequency of the reference currents was changed from 60 Hz to 90 Hz. In the step changes of magnitude and fundamental frequency of the reference currents, the output currents of the proposed method tracked the reference change as fast as those of the conventional method with half the sampling period.   Figure 11. Experimental results of the three-phase output currents (ia, ib, and ic) and the a-leg upper switch S1 for the current magnitude step change from 6 A to 12 A and frequency step change from 60 Figure 11. Experimental results of the three-phase output currents (i a , i b , and i c ) and the a-leg upper switch S 1 for the current magnitude step change from 6 A to 12 A and frequency step change from 60 Hz to 90 Hz(|I*| = 12 A). (a) the proposed method with T s = 250 µs, and (b) the conventional method with T s = 125 µs.

Performance Comparison
In order to evaluate the superiority of the proposed method in terms of efficiency, experimental data of loss and power efficiency were measured in the conventional and the proposed methods. The losses were measured with a power meter (Precision Power Analyzer 5500-Newtons 4th) in the experimental setup. The three conventional MPCC methods were compared: (1) Conv. (1 vector): MPCC method using a single optimal voltage vector [3] with the same sampling period with the proposed method. (2) Conv. (1 vector): MPCC method using a single optimal voltage vector [3] with half sampling period of the proposed method (3) Conv. (2 vectors): MPCC method using two optimal voltage vectors [19]. Figure 12 shows the power efficiency of the conventional methods and the proposed method measured by PPA. As shown in Figure 11, the power efficiency of the conventional method using one vector with T s = 125 µs showed the worst efficiency, and it is because the sampling frequency is twice as large as the other methods. When the sampling period was equal to 250 µs, the conventional method using two vectors showed the worst efficiency because this method uses two vectors during one sampling period. It can be clearly seen that the proposed method resulted in the highest efficiency because the clamping operation of switches exposed to highest current conduction can reduce the switching losses. Table 2 shows the performance comparison of the conventional methods and the proposed method. The conventional method using one vector during one sampling period showed high THD values of the output currents. To overcome this drawback, control methods using two vectors during one sampling period have been studied [19]. Although this method can improve the current quality by decreasing THD values of the output currents, it increases switching losses because of the usage of two vectors during one sampling period. The proposed method not only improved the THD values of the output currents by using two optimal vectors during one sampling period but also reduced Figure 11. Experimental results of the three-phase output currents (i a , i b , and i c ) and the a-leg upper switch S 1 for the current magnitude step change from 6 A to 12 A and frequency step change from 60 Hz to 90 Hz(|I*| = 12 A). (a) the proposed method with T s = 250 µs, and (b) the conventional method with T s = 125 µs.

Performance Comparison
In order to evaluate the superiority of the proposed method in terms of efficiency, experimental data of loss and power efficiency were measured in the conventional and the proposed methods. The losses were measured with a power meter (Precision Power Analyzer 5500-Newtons 4th) in the experimental setup. The three conventional MPCC methods were compared: (1) Conv. (1 vector): MPCC method using a single optimal voltage vector [3] with the same sampling period with the proposed method. (2) Conv. (1 vector): MPCC method using a single optimal voltage vector [3] with half sampling period of the proposed method (3) Conv. (2 vectors): MPCC method using two optimal voltage vectors [19]. Figure 12 shows the power efficiency of the conventional methods and the proposed method measured by PPA. As shown in Figure 11, the power efficiency of the conventional method using one vector with T s = 125 µs showed the worst efficiency, and it is because the sampling frequency is twice as large as the other methods. When the sampling period was equal to 250 µs, the conventional method using two vectors showed the worst efficiency because this method uses two vectors during one sampling period. It can be clearly seen that the proposed method resulted in the highest efficiency because the clamping operation of switches exposed to highest current conduction can reduce the switching losses. Table 2 shows the performance comparison of the conventional methods and the proposed method. The conventional method using one vector during one sampling period showed high THD values of the output currents. To overcome this drawback, control methods using two vectors during one sampling period have been studied [19]. Although this method can improve the current quality by decreasing THD values of the output currents, it increases switching losses because of the usage of two vectors during one sampling period. The proposed method not only improved the THD values of the output currents by using two optimal vectors during one sampling period but also reduced the switching loss by clamping switches exposed to highest current conduction. On the other hand, the execution time in the digital signal processor of the proposed method could be higher than the conventional methods because of the complexity of the proposed method, which needs the technique to pre-select the vectors and choose two vectors.
In addition, Tables 3 and 4 show performance comparisons between the conventional method and the proposed method in the simulation and experimental results, respectively. As shown in Tables 3  and 4, the THD values and the total losses of the proposed method were lower than those of the conventional method, with the half sampling frequency of the proposed method. Table 5 depicts the results of the proposed method obtained in simulation and experimental result, where data were almost similar. one vector with Ts = 125 μs showed the worst efficiency, and it is because the sampling frequency is twice as large as the other methods. When the sampling period was equal to 250 μs, the conventional method using two vectors showed the worst efficiency because this method uses two vectors during one sampling period. It can be clearly seen that the proposed method resulted in the highest efficiency because the clamping operation of switches exposed to highest current conduction can reduce the switching losses.  Table 2 shows the performance comparison of the conventional methods and the proposed method. The conventional method using one vector during one sampling period showed high THD Power efficiency

Conclusions
In this paper, the MPCC method for three-phase VSI based on the utilization of two vectors has been proposed to regulate output currents with small current errors and current ripples. Also, the vector pre-selection method is applied to the proposed method to reduce the VSIs switching losses by pre-selecting the four available voltage vectors. In the proposed method, selecting two future voltage vectors among the four available voltage vectors and dividing them in a future sampling period are determined by the newly defined cost function and the optimization process. The proposed method reduces the total loss versus the conventional MPCC method with the same switching frequency due to the optimal process involved in selecting two vectors and their time durations. In addition, the proposed method and the conventional method show the same performance in current error and THD. The effectiveness of the proposed method is verified from the study of the simulation and experiment of the three-phase VSIs.