A Modulated Model Predictive Current Controller for Interior Permanent-Magnet Synchronous Motors

: Model predictive current controllers (MPCCs) are widely applied in motor drive control and operations. To date, however, the presence of large current errors in conventional predictive current control remains a signiﬁcant predicament, due to harmonic distortions and current ripples. Naturally, noticeable current estimation inaccuracies lead to poor performance. To improve the above situation, a modulated model predictive current controller (MMPCC) is proposed for interior permanent-magnet synchronous motors (IPMSMs) in this paper. Two successive voltage vectors will be applied in a sampling period to greatly boost the number of candidate switching modes from seven to thirteen. A cost function, which is deﬁned as the quadratic sum of current prediction errors, is employed to ﬁnd an optimal switching mode and an optimized duty ratio to be applied in the next sampling period, such that the cost value is minimal. The e ﬀ ectiveness of the proposed method is veriﬁed through eight experiments using a TMS320F28379D microcontroller, and performance comparisons are made against an existing MPCC. In terms of quantitative improvements made to the MPCC, the proposed MMPCC reduces its current ripple and total harmonic distortion (THD) by, on average, 27.17% and 21.84%, respectively.


Introduction
Due to its advantages of high reliability, high power density, and minimal noise, the interior permanent-magnet synchronous motor (IPMSM) is regarded as a superior solution in a lot of electrical drives and industry applications [1,2]. The motor is characterized by its unique salient-pole structure for precision operation at diverse speeds. However, the motor also requires a highly dynamic control platform to take full advantage of these merits. Among the massive exploration of control strategies brought about by the rapid development of fast and inexpensive digital signal processors (DSPs), the well-known model predictive control (MPC) [3,4] stands out among the existing methods. The MPC holds the edge in terms of its predictive mechanism, its cost function formulation [5], and its online optimization.
Over the years, various applications in power converters and electrical drives have been carried out using model predictive controllers. Such improvements have elevated the status quo of control strategies, including IPMSM applications. However, there still remain underlying drawbacks in the conventional MPC, including scarcely regulated actuations [6], and the large accumulation of current and torque ripples [7,8]. This is primarily because of the limitations of voltage vector selections from the voltage source inverter (VSI). On the other hand, current tracking performance [9,10] has been put forward to improve the selection process of voltage vectors using the prediction error correction and improved deadbeat predictive current regulation. However, low current spikes remain difficult to achieve due to current update stagnation.
The various frameworks of MPC have also gained a great deal of attention in the optimization of flux weakening [11] and enhancement of dc link utilization [12]. The latter aims to expand higher speed operation beyond the nominal range by controlling the field strength and reducing the back electromotive force (EMF). However, this results in other problems concerning system robustness and dynamic response, because of the sensitivity to parameter variations and voltage saturations.
To overcome the above difficulties, the duty cycle modulation (DCM) technique has been considered. This strategy involves the selection of two or more active voltage vectors and applies them in a sampling period. Compared to the conventional MPC, the method achieves significant improvements in ripple reduction and prediction error [13][14][15]. However, the dynamic response is largely affected by the heavy on-line computation. A simple modulation strategy based on an adaptive saturation mechanism was adopted in [16] to tackle the problem. The torque error was found to be significantly reduced with an impressively decreased computational requirement; but the harmonics in the stator current waveform remain largely unimproved.
Studies have shown that predictive current control possesses an impressive capability for steady-state current tracking, thanks to high-speed processors [17][18][19] that overcome the limitation imposed by the computational load, hence making a lot of progress in both theory and application. This advantage has been successfully implemented to circumvent the shortcomings of load parameter dependency in predictive algorithms [20,21]. In contrast to traditional model-based predictive controllers, these methods are able to mitigate the effect of imprecise mathematical modeling in the prediction process. Current changes from the applied voltage vectors are utilized to predict future currents. A similar technique is adopted in [22], where two active voltage vectors with equal durations were applied in a sampling period. However, its efficiency is uncertain to some extent, due to the large number of candidate voltage vectors to be optimized. In addition, the durations of the two successive voltage vectors remain fixed.
To date, modulation has been widely employed in model predictive control. Equipped with modern microcontrollers, MPC has implemented this technology in a variety of applications, including stand-alone voltage source converters [23], multilevel converters of multi-terminal DC systems [24], and active-front-end rectifiers [25], to name a few. These implementations resulted in significant improvements to converters and rectifiers, but not IPMSMs. In [26], a linear range modulation was applied to the RL load to optimize the overmodulation caused by ineffective current tracking within the modulation zone. The application of a multilevel AC/AC converter is introduced in [27], in which steady-state multi-objective current tracking performance was analyzed. However, given the limited candidate switching states, the complexity of selecting the optimal active voltage vector becomes a disadvantage.
To improve the MPC as described in [20], this paper proposes a modulated model predictive current controller (MMPCC) for IPMSMs, whose main contributions are highlighted as follows: 1 In a sampling period, the applications of two successive voltage vectors are allowed to operate together with two optimal duty ratios while achieving the current predictions. It is worth noting that these optimized duty ratios are calculated online. 2 More candidate switching modes are obtained via the linear combinations of two voltage vectors, increased from the typical seven switching modes described in [20] to thirteen when using the proposed control strategy, thus enhancing control efficiency while retaining the simplicity of the algorithm. 3 The duty ratio is taken as a variable; hence, it is incorporated into the cost function and optimized. 4 This is the first implementation of such an MMPCC on the IPMSM drive system.
The rest of the article is organized as follows. An MPC-based current control is presented in Section 2, followed by Section 3, which discusses the proposed scheme in detail. The experimental platform, setup and results are provided in Section 4, and finally, the conclusion is given in Section 5.

Current Control Prediction Based on Conventional MPC
The algorithm of a typical model predictive current control (MPCC) is primarily based on the discretization of the associated mathematical model to predict the currents. Based on the criterion of cost function minimization [3], data are collected from the motor operation and then utilized so as to identify the optimal actuation.
The schematic diagram of a typical three-phase six-switch inverter-fed interior permanent-magnet synchronous motor is shown in Figure 1. The voltage source inverter contains six switches controlled by signals S a , S b , and S c . Other components include: the dc-link voltage, V DC ; the stator voltages, v a , v b , v c ; the stator currents, i a , i b , i c ; the extended back EMF, e a , e b , e c ; the stator resistance, r s ; and the q-axis inductance, L q . The algorithm of a typical model predictive current control (MPCC) is primarily based on the discretization of the associated mathematical model to predict the currents. Based on the criterion of cost function minimization [3], data are collected from the motor operation and then utilized so as to identify the optimal actuation.
The schematic diagram of a typical three-phase six-switch inverter-fed interior permanentmagnet synchronous motor is shown in Figure 1. The voltage source inverter contains six switches controlled by signals Sa, Sb, and Sc. Other components include: the dc-link voltage, VDC; the stator voltages, va, vb, vc; the stator currents, ia, ib, ic; the extended back EMF, ea, eb, ec; the stator resistance, rs; and the q-axis inductance, Lq.

Mathematical Model of IPMSM
The stator voltage equation of IPMSM expressed in the three-phase stationary reference frame is shown in Equation (1). However, due to its highly complex characteristics, the differential equation of electrical dynamics is given below under an ideal condition wherein the current loss, the hysteresis loss, and the magnetic saturations are assumed to be negligible [28].
where refers to the differential operator, stands for the phase-wise stator current, and represents the phase-wise stator flux linkage produced by the self-linked flux of stator currents and rotor magnets. The cross-coupling effects between the inductances brought by make the IPMSM a highly nonlinear system, and hence complex in nature. The following method can be employed to simplify the above expression, leading to a two-axis model [20].
where is the rotor flux linkage and is the rotor angle of the motor. The effect of crosscoupling in the two-axis model is considered to be the self-inductances represented by , , and , and the mutual inductances by , , , , and . The respective values can be obtained through

Mathematical Model of IPMSM
The stator voltage equation of IPMSM expressed in the three-phase stationary reference frame is shown in Equation (1). However, due to its highly complex characteristics, the differential equation of electrical dynamics is given below under an ideal condition wherein the current loss, the hysteresis loss, and the magnetic saturations are assumed to be negligible [28].
where d dt refers to the differential operator, i x stands for the phase-wise stator current, and ψ z represents the phase-wise stator flux linkage produced by the self-linked flux of stator currents and rotor magnets. The cross-coupling effects between the inductances brought by ψ z make the IPMSM a highly nonlinear system, and hence complex in nature. The following method can be employed to simplify the above expression, leading to a two-axis model [20].
where ψ m is the rotor flux linkage and θ re is the rotor angle of the motor. The effect of cross-coupling in the two-axis model is considered to be the self-inductances represented by L aa , L bb , and L cc , and the mutual inductances by L ab , L ba , L bc , L cb , L ac and L ca . The respective values can be obtained through The variables L A , L B , and L ls represent the DC component of the inductance, the magnitude of the cosine function of the inductance, and the coil leakage inductance, respectively. Given Equations (2)-(4), Equation (1) can be expanded to yield the following three-phase stator voltage expressions [20]: The transformation of motor inductances between stator and rotor can be simplified through the equivalent quadrature value of inductance in the two-axis model. The q-axis inductance, L q , is taken to be constant in this article and is defined as in [21].
Given the definition as of Equation (8), the three-phase stator voltage can be simplified as [21]: where e z is the phase-wise extended back EMF. Moreover, the a − b − c three-phase reference frame of IPMSM is not orthogonal in nature and the asymmetrical components are difficult to control. It is impossible to vary one axis without affecting the others. To circumvent these constraints and simplify the analysis, the α − β transformation is employed, which is widely adopted. This will reduce the projections of the phase quantities into the two-axis model reference frame. Following the α − β transformation, the counterpart of Equation (9) can be written as:

Optimal Voltage Vector Selection
The field-oriented control (FOC) is commonly used in motor drive systems powered by three-phase six-switch voltage source inverters (VSIs). Given the six-switch topology, a VSI can generate eight basic voltage vectors in general. Together with linear controllers [3], the FOC provides the three-phase voltages required for the IPMSM drive systems. Unlike FOC, the MPC determines the future voltage vector based on the predefined selection of switching states. Ideally, the optimal selection of the future voltage vector should yield the lowest cost function value.
The eight switching voltage vectors are shown in Figure 2, containing two zero-voltage vectors, V 0 and V 7 , and six non-zero or the active voltage vectors, V 1 , V 2 , V 3 , V 4 , V 5 , and V 6 . Associated with them are the switching states, S 0 -S 7 . Given in Table 1 are the details of the relationships between the switching states, the three-phase voltage outputs, the equivalent stator voltages, and the switching function. switching states, the three-phase voltage outputs, the equivalent stator voltages, and the switching function.
It can also be observed from Table 1 that the output voltages of and are zero, and they are, in theory, considered identical. This reduces the number of total switching states from eight to seven, and they are , , , , , , and .

Current Prediction Strategy
Compared to FOC, the MPC possesses certain flexibility and extensibility features. In addition, thanks to the performance advancement of modern microcontrollers, the MPC is able to predict the stator currents under all candidate switching modes within a sampling period. For the MPC to be realized and its operation carried out on a microcontroller, all continuous variables are approximated by discretized forms. As an example, the time derivative of current is approximated by: assuming the sampling period is short enough. Given the above condition, Equations (10) and (11) can be expressed in an equivalent discretized form as: It can also be observed from Table 1 that the output voltages of V 0 and V 7 are zero, and they are, in theory, considered identical. This reduces the number of total switching states from eight to seven, and they are S 0 , S 1 , S 2 , S 3 , S 4 , S 5 , and S 6 .

Current Prediction Strategy
Compared to FOC, the MPC possesses certain flexibility and extensibility features. In addition, thanks to the performance advancement of modern microcontrollers, the MPC is able to predict the stator currents under all candidate switching modes within a sampling period. For the MPC to be realized and its operation carried out on a microcontroller, all continuous variables are approximated by discretized forms. As an example, the time derivative of current is approximated by: assuming the sampling period T s is short enough. Given the above condition, Equations (10) and (11) can be expressed in an equivalent discretized form as: where x ∈ α, β and k stands for the (k)th period. With time delay compensation taken into account [29], the above equations can be rearranged to yield the (k + 1)th and (k + 2)th predicted currents as follows: where the superscripts "p" and "ˆ" represent the predicted current and estimated extended back-EMF, respectively. To circumvent the difficulty of getting the value ofê x (k + 1), one may again assume that the sampling period is sufficiently short and use: Given Equation (15), Equation (16) can be rewritten as where The voltage v x (k + 1), appearing in Equation (17), can be found from Table 1. To achieve optimal switching, the predicted current should be as close as possible to the current command. To that end, a cost function g is defined as follows: where the superscript "*" refers the command or, interchangeably, the reference. The switching state yielding the minimum cost value will then be selected and applied in the next sampling period. Details of the method are shown in the schematic diagram of Figure 3. where ∈ { , } and stands for the ( )th period. With time delay compensation taken into account [29], the above equations can be rearranged to yield the ( + 1)th and ( + 2)th predicted currents as follows: where the superscripts "p" and "^" represent the predicted current and estimated extended back-EMF, respectively. To circumvent the difficulty of getting the value of ( + 1), one may again assume that the sampling period is sufficiently short and use: Given Equation (15), Equation (16) can be rewritten as . The voltage ( + 1), appearing in Equation (17), can be found from Table 1. To achieve optimal switching, the predicted current should be as close as possible to the current command. To that end, a cost function is defined as follows: where the superscript "*" refers the command or, interchangeably, the reference. The switching state yielding the minimum cost value will then be selected and applied in the next sampling period. Details of the method are shown in the schematic diagram of Figure 3.

Proposed Modulated Model Predictive Current Control
Conventional MPC typically applies a single voltage vector in each sampling period, as shown in Figure 4a, and produces large current ripples when it is applied to the IPMSM drive system. To improve this situation, a modulated model predictive current control (MMPCC) is proposed, as shown in Figure 4b. The basic idea behind this is to make the voltage application time adaptive and optimized by incorporating the mechanism of vector modulation into the MPC. It can be seen from Figure 4b that two voltage vectors and are applied successively for current prediction in a sampling period, with durations and , respectively. One gains apparent benefits by using two modulated voltage vectors whose application durations are adapted and optimized online. This

Proposed Modulated Model Predictive Current Control
Conventional MPC typically applies a single voltage vector in each sampling period, as shown in Figure 4a, and produces large current ripples when it is applied to the IPMSM drive system. To improve this situation, a modulated model predictive current control (MMPCC) is proposed, as shown in Figure 4b. The basic idea behind this is to make the voltage application time adaptive and optimized by incorporating the mechanism of vector modulation into the MPC. It can be seen from Figure 4b that two voltage vectors V c1 s and V c2 s are applied successively for current prediction in a sampling period, with durations T c1 k and T c2 k , respectively. One gains apparent benefits by using two modulated voltage vectors whose application durations are adapted and optimized online. This is a very important and appealing feature, in contrast to conventional MPCCs without modulation, whose voltage application time is fixed. To reduce the difficulty of implementation and the computation time, thirteen switching modes were considered for current predictions.
Energies. 2019, 12, x ; doi: FOR PEER REVIEW 7 of 20 is a very important and appealing feature, in contrast to conventional MPCCs without modulation, whose voltage application time is fixed. To reduce the difficulty of implementation and the computation time, thirteen switching modes were considered for current predictions.

Design of Adaptive Duty Ratio Modulation
Taking the (k)th sampling period, for instance, the two durations are represented by and , respectively, for the two consecutive voltage vectors, and . The superscript " " refers to the resulting synthesized value of modulation. Accordingly, the sampling period can be expressed as: The two duty ratios corresponding to and are designated and , respectively, and are defined as: As a result of Equations (20) and (21), the sum of the two duty ratios and satisfies the following identity: Further details of Equation (19) are discussed in Section 3.3. Due to practical reasons and the limitations of the motor drive system, oftentimes, lower and upper limits are imposed on the duty ratios and . That is, 0.2 < ≤ 0.8 and = − .

Synthesized Voltage Vectors and Current Prediction
To implement the proposed method and reduce the computational complexity, thirteen synthesized voltage vectors, labeled as -, are used, each of which is synthesized by two different or identical voltage vectors. As shown in Figure 5, thirteen synthesized voltage vectors are projected onto the α-β plane, including a zero-voltage vector at the center, six co-adjacent nonzero voltage vectors -, and six non-zero voltage vectors across the edges -. As can be seen from Figure 4b, the application times and are adaptive, meaning that the two duty ratios and are adaptive. As such, the value of the (k)th synthesized voltage ( ) can be expressed as

Design of Adaptive Duty Ratio Modulation
Taking the (k)th sampling period, for instance, the two durations are represented by T c1 k and T c2 k , respectively, for the two consecutive voltage vectors, V c1 s and V c2 s . The superscript "c" refers to the resulting synthesized value of modulation. Accordingly, the sampling period T s can be expressed as: The two duty ratios corresponding to V c1 s and V c2 s are designated D c1 k and D c1 k , respectively, and are defined as: As a result of Equations (20) and (21), the sum of the two duty ratios D c1 k and D c1 k satisfies the following identity: Further details of Equation (19) are discussed in Section 3.3. Due to practical reasons and the limitations of the motor drive system, oftentimes, lower and upper limits are imposed on the duty ratios D c1

Synthesized Voltage Vectors and Current Prediction
To implement the proposed method and reduce the computational complexity, thirteen synthesized voltage vectors, labeled as V c 01 -V c 12 , are used, each of which is synthesized by two different or identical voltage vectors. As shown in Figure 5, thirteen synthesized voltage vectors are projected onto the α-β plane, including a zero-voltage vector V c 0 at the center, six co-adjacent non-zero voltage vectors V c 1 -V c 6 , and six non-zero voltage vectors across the edges V c 7 -V c 12 . As can be seen from Figure 4b, the application times T c1 k and T c2 k are adaptive, meaning that the two duty ratios D c1 k and D c2 k are adaptive. As such, the value of the (k)th synthesized voltage v c x (k) can be expressed as where V c1 x (k) and V c2 x (k) are the stator voltages corresponding to the two consecutive voltage vectors applied in the (k)th sampling period. Equation (22) implies: Next, substitution of (24) into (23) yields the following equation: Equation (25) can be utilized to obtain the (k − 1)th and (k + 1)th synthesized voltage vectors v c Please note that the voltage vectors V 0 ∼ V 7 correspond to the switching states S 0 -S 7 , as shown in Table 1. The thirteen synthesized voltage vectors and the two applied switching states are listed in Table 2.

Cost Function Design and Optimized Modulation Ratios
The proposed MMPCC retains an important criterion of model predictive control-the cost function design. In reference to Equations (11) and (15) 1)). Ideally, the optimal switching mode is that which yields the minimum cost value. However, the above calculation cannot be implemented if the duty ratio is unknown. To proceed, the optimal duty ratio needs to be determined first. Assuming a partial derivative of the cost function with respect to the duty ratio exists, one can easily obtain the following expression: Let = 0, one then gets the following optimal duty ratio, which is denoted as With the duty ratio in Equation (31) replaced by * , the cost function now becomes: * = ( It is worth mentioning that different synthesized voltages yield different optimal duty ratios. As such, the thirteen optimal duty ratios corresponding to the thirteen switching modes need to be

Cost Function Design and Optimized Modulation Ratios
The proposed MMPCC retains an important criterion of model predictive control-the cost function design. In reference to Equations (11) and (15), a new cost function G applicable to the MMPCC is defined as: Using Equations (25)- (27) and substituting (28) and (29) into (30), one can rewrite (30) as: . Ideally, the optimal switching mode is that which yields the minimum cost value. However, the above calculation cannot be implemented if the duty ratio D c1 k+1 is unknown. To proceed, the optimal duty ratio D c1 k+1 needs to be determined first. Assuming a partial derivative of the cost function G with respect to the duty ratio D c1 k+1 exists, one can easily obtain the following expression: Let ∂G ∂D c1 k+1 = 0, one then gets the following optimal duty ratio, which is denoted as D c1 * k+1 : With the duty ratio D c1 k+1 in Equation (31) replaced by D c1 * k+1 , the cost function now becomes: It is worth mentioning that different synthesized voltages yield different optimal duty ratios. As such, the thirteen optimal duty ratios corresponding to the thirteen switching modes need to be calculated online first. These optimal values are then substituted back into the cost function for minimization. Given Equation (34), thirteen cost values can be obtained by means of the thirteen synthesized voltage vectors listed in Table 2. The vector is considered to be optimal so long as it gives the minimum cost value. Given Table 2, the two switching states can be found through table-lookup and applied for durations T c1 * k+1 and T c2 * k+1 , respectively, in the next sampling period. Specifically, the optimal durations T c1 * k+1 and T c2 * k+1 are defined as The fundamental principle of the proposed MMPCC is illustrated in detail by Figure 6. The novel design of the cost function calculation in Equation (34) allows thirteen unique predicted values of stator currents, of which the cost value with optimal durations, as described in Equations (35) and (36), will be applied in the next sampling period. As a result, the optimal switching mode will naturally lead to the smallest current error in the (k + 2)th sampling period.  Table 2. The vector is considered to be optimal so long as it gives the minimum cost value. Given Table 2, the two switching states can be found through table-lookup and applied for durations * and * , respectively, in the next sampling period. Specifically, the optimal durations * and * are defined as * = * The fundamental principle of the proposed MMPCC is illustrated in detail by Figure 6. The novel design of the cost function calculation in Equation (34) allows thirteen unique predicted values of stator currents, of which the cost value with optimal durations, as described in Equations (35) and (36), will be applied in the next sampling period. As a result, the optimal switching mode will naturally lead to the smallest current error in the (k + 2)th sampling period.

Experimental Test Stand and Results
To validate the proposed MMPCC and evaluate its effectiveness, several experiments were carried out through the microcontroller TMS320F28379D of Texas Instruments (Dallas, TX, United States), and performance comparisons were made against the MPCC without modulation. The hardware drive circuit is shown in Figure 7, containing a power meter, an IPMSM, and including several other important components: Part A is composed of two current detection boards used to measure the a and b phase currents, respectively, and each board contains a current sensor of LA25-NP type by LEM Co. Ltd. (Plan-les-Ouates, Geneva, Switzerland), a current/voltage conversion circuit, and a 16-bit A/D converter. The bandwidth of the LA25-NP is 150 kHz, with ±25 A as its range. Part B includes an overcurrent protection circuit, an encoder circuit, and several HCPL2631 optocouplers made by Fairchild Semiconductor. The overcurrent protection circuit protects against excessive current, disabling the system instantaneously to avoid damage. Part C is the TMS320F28379D dual-core microcontroller made by Texas Instruments. The microcontroller consists of two sets of TMS320C28x 32-bit CPUs capable of high-speed data processing at up to 200 MHz. Part D consists of a low switching dissipation insulated gate bipolar transistors (IGBT) power module inverter (SCM1246MF) made by Sanken Electric Co. Ltd. (Niiza-shi, Saitama-ken, Japan), with an output breakdown voltage of 600 V and output current of 30 A. The inverter includes six sets of IGBTs

Experimental Test Stand and Results
To validate the proposed MMPCC and evaluate its effectiveness, several experiments were carried out through the microcontroller TMS320F28379D of Texas Instruments (Dallas, TX, United States), and performance comparisons were made against the MPCC without modulation. The hardware drive circuit is shown in Figure 7, containing a power meter, an IPMSM, and including several other important components: Part A is composed of two current detection boards used to measure the a and b phase currents, respectively, and each board contains a current sensor of LA25-NP type by LEM Co. Ltd. (Plan-les-Ouates, Geneva, Switzerland), a current/voltage conversion circuit, and a 16-bit A/D converter. The bandwidth of the LA25-NP is 150 kHz, with ±25 A as its range. Part B includes an overcurrent protection circuit, an encoder circuit, and several HCPL2631 optocouplers made by Fairchild Semiconductor. The overcurrent protection circuit protects against excessive current, disabling the system instantaneously to avoid damage. Part C is the TMS320F28379D dual-core microcontroller made by Texas Instruments. The microcontroller consists of two sets of TMS320C28x 32-bit CPUs capable of high-speed data processing at up to 200 MHz. Part D consists of a low switching dissipation insulated gate bipolar transistors (IGBT) power module inverter (SCM1246MF) made by Sanken Electric Co. Ltd. (Niiza-shi, Saitama-ken, Japan), with an output breakdown voltage of 600 V and output current of 30 A. The inverter includes six sets of IGBTs and freewheeling diodes, and three sets of pre-driver monolithic integrated circuits. Part E contains DC power supplies, RT-50C, providing two DC output channels of 5 V and 15 V. The device is equipped with built-in protection against short-circuit, overload, overvoltage, and capable of withstanding 300 V AC surge without damage. The IPMSM parameters are given in Table 3 while the corresponding composite parameters appearing in Equation (17) are listed in Table 4. equipped with built-in protection against short-circuit, overload, overvoltage, and capable of withstanding 300 V AC surge without damage. The IPMSM parameters are given in Table 3 while the corresponding composite parameters appearing in Equation (17) are listed in Table 4. To achieve a fair comparison, the sampling periods of both PCCs were set to be the same, at 100 microseconds, within which the proposed MMPCC and the MPCC can be executed. Respectively, the proposed MMPCC takes 62 microseconds, while the MPCC takes 22 microseconds to complete the execution.

Steady-State Response
To illustrate the response of the proposed method, steady-state performance was evaluated in the α-β reference frame with stator current commands running at different frequencies.

Steady-State Response
To achieve a fair comparison, the sampling periods of both PCCs were set to be the same, at 100 microseconds, within which the proposed MMPCC and the MPCC can be executed.
Respectively, the proposed MMPCC takes 62 microseconds, while the MPCC takes 22 microseconds to complete the execution.
To illustrate the response of the proposed method, steady-state performance was evaluated in the α-β reference frame with stator current commands running at different frequencies. Figure 8 shows the α-β current waveforms operating at an amplitude of 4 A and frequency of 30 Hz. As can be seen from Figure 8a, significant ripples exist at current peaks. Specifically, the current ripple is 0.1077 A, and the THD is 0.0812%. Figure 8b shows that the α-axis and β-axis currents of the proposed MMPCC follow the current commands very well, with current ripple being equal to 0.0705 A, and the THD being equal to 0.0812%. In terms of the improvement made to the MPCC, the proposed MMPCC reduces the current ripple by 34.54% and decreases the THD by 41.89%, according to Figure 8. shows the α-β current waveforms operating at an amplitude of 4 A and frequency of 30 Hz. As can be seen from Figure 8a, significant ripples exist at current peaks. Specifically, the current ripple is 0.1077 A, and the THD is 0.0812%. Figure 8b shows that the α-axis and β-axis currents of the proposed MMPCC follow the current commands very well, with current ripple being equal to 0.0705 A, and the THD being equal to 0.0812%. In terms of the improvement made to the MPCC, the proposed MMPCC reduces the current ripple by 34.54% and decreases the THD by 41.89%, according to Figure  8.
The next experiment is for performance evaluation at a low operating frequency. The experimental results of tracking the α-β current commands with amplitude of 4 A and frequency of 10 Hz are shown in Figure 9. As can be observed, Figure 9a exhibits more pronounced current ripples at peaks compared to Figure 9b. Particularly, the current ripple and THD produced by the MMPCC are reduced by 32.92% and 11.54% in terms of the improvement made to the MPCC. These results clearly demonstrate that the proposed MMPCC has better tracking performance under steady-state operations.

Transient Response
Among the many factors that need to be considered in terms of the operations of converters and electrical drives, two crucial ones stand out: fast response, and rejection of ambient electromagnetic disturbances. However, here, the focus is only placed on the former, as the latter is beyond the scope of this study. The transient response of the motor driver's output current is highly dependent on the current controller. To demonstrate the performance in this regard, a simple test is conducted The next experiment is for performance evaluation at a low operating frequency. The experimental results of tracking the α-β current commands with amplitude of 4 A and frequency of 10 Hz are shown in Figure 9. As can be observed, Figure 9a exhibits more pronounced current ripples at peaks compared to Figure 9b. Particularly, the current ripple and THD produced by the MMPCC are reduced by 32.92% and 11.54% in terms of the improvement made to the MPCC. These results clearly demonstrate that the proposed MMPCC has better tracking performance under steady-state operations. shows the α-β current waveforms operating at an amplitude of 4 A and frequency of 30 Hz. As can be seen from Figure 8a, significant ripples exist at current peaks. Specifically, the current ripple is 0.1077 A, and the THD is 0.0812%. Figure 8b shows that the α-axis and β-axis currents of the proposed MMPCC follow the current commands very well, with current ripple being equal to 0.0705 A, and the THD being equal to 0.0812%. In terms of the improvement made to the MPCC, the proposed MMPCC reduces the current ripple by 34.54% and decreases the THD by 41.89%, according to Figure  8. The next experiment is for performance evaluation at a low operating frequency. The experimental results of tracking the α-β current commands with amplitude of 4 A and frequency of 10 Hz are shown in Figure 9. As can be observed, Figure 9a exhibits more pronounced current ripples at peaks compared to Figure 9b. Particularly, the current ripple and THD produced by the MMPCC are reduced by 32.92% and 11.54% in terms of the improvement made to the MPCC. These results clearly demonstrate that the proposed MMPCC has better tracking performance under steady-state operations.

Transient Response
Among the many factors that need to be considered in terms of the operations of converters and electrical drives, two crucial ones stand out: fast response, and rejection of ambient electromagnetic disturbances. However, here, the focus is only placed on the former, as the latter is beyond the scope of this study. The transient response of the motor driver's output current is highly dependent on the current controller. To demonstrate the performance in this regard, a simple test is conducted

Transient Response
Among the many factors that need to be considered in terms of the operations of converters and electrical drives, two crucial ones stand out: fast response, and rejection of ambient electromagnetic disturbances. However, here, the focus is only placed on the former, as the latter is beyond the scope of this study. The transient response of the motor driver's output current is highly dependent on the current controller. To demonstrate the performance in this regard, a simple test is conducted containing an abrupt change in the current command during operations. Figure 10 shows the instantaneous behaviors where the α-axis current command undergoes a sudden jump from −4 A to 4 A at 0.05 s. As can be seen from Figure 10, both methods exhibit good transient responses. In this test, MMPCC performed slightly better than did the MPCC without modulation. Specifically, compared with the latter, the former reduced current ripple and THD by only 1.76% and 0.15%, respectively. A possible explanation for this is that the duration of comparison was too short, at only 0.1s. Additionally, the transient current error due to the sudden change in current command is much larger than the current ripple of interest. Figure 11, on the other hand, illustrates the transient performance of α-βcurrent command jumping from 1 A to 4 A at 0.2 s. Similar to the previous transient test, it is clearly visible that better current tracking performance is achieved by the proposed MMPCC. Current ripple is impressively reduced by 8.98%, while THD is slightly decreased by 0.68%. containing an abrupt change in the current command during operations. Figure 10 shows the instantaneous behaviors where the α-axis current command undergoes a sudden jump from −4 A to 4 A at 0.05 s. As can be seen from Figure 10, both methods exhibit good transient responses. In this test, MMPCC performed slightly better than did the MPCC without modulation. Specifically, compared with the latter, the former reduced current ripple and THD by only 1.76% and 0.15%, respectively. A possible explanation for this is that the duration of comparison was too short, at only 0.1s. Additionally, the transient current error due to the sudden change in current command is much larger than the current ripple of interest. Figure 11, on the other hand, illustrates the transient performance of α-β current command jumping from 1 A to 4 A at 0.2 s. Similar to the previous transient test, it is clearly visible that better current tracking performance is achieved by the proposed MMPCC. Current ripple is impressively reduced by 8.98%, while THD is slightly decreased by 0.68%.

Load-Torque and Speed Varying Condition
Further evaluation is conducted at various load-torque and speed conditions. To better illustrate the difference between the two current controllers in terms of their performances under steady-state current tracking, the waveforms of α-β current errors are provided.
Given in Figures 12 and 13 are the results of IPMSM running at a constant speed of 500 rpm under different load-torques of 1 N-m and 2 N-m, respectively. Figure 12a shows the corresponding current waveforms in which the MPCC generates 0.2463 A of current ripples and 5.8826% of THD. In that regard, the proposed MMPCC in yields smaller ripples and THD with high percentage performance improvements up to 44.05% and 62.29%, respectively, as depicted in Figure 12b. In addition, the current errors as shown in Figure 12d,f are significantly smaller than those in Figure containing an abrupt change in the current command during operations. Figure 10 shows the instantaneous behaviors where the α-axis current command undergoes a sudden jump from −4 A to 4 A at 0.05 s. As can be seen from Figure 10, both methods exhibit good transient responses. In this test, MMPCC performed slightly better than did the MPCC without modulation. Specifically, compared with the latter, the former reduced current ripple and THD by only 1.76% and 0.15%, respectively. A possible explanation for this is that the duration of comparison was too short, at only 0.1s. Additionally, the transient current error due to the sudden change in current command is much larger than the current ripple of interest. Figure 11, on the other hand, illustrates the transient performance of α-β current command jumping from 1 A to 4 A at 0.2 s. Similar to the previous transient test, it is clearly visible that better current tracking performance is achieved by the proposed MMPCC. Current ripple is impressively reduced by 8.98%, while THD is slightly decreased by 0.68%.

Load-Torque and Speed Varying Condition
Further evaluation is conducted at various load-torque and speed conditions. To better illustrate the difference between the two current controllers in terms of their performances under steady-state current tracking, the waveforms of α-β current errors are provided.
Given in Figures 12 and 13 are the results of IPMSM running at a constant speed of 500 rpm under different load-torques of 1 N-m and 2 N-m, respectively. Figure 12a shows the corresponding current waveforms in which the MPCC generates 0.2463 A of current ripples and 5.8826% of THD. In that regard, the proposed MMPCC in yields smaller ripples and THD with high percentage performance improvements up to 44.05% and 62.29%, respectively, as depicted in Figure 12b. In addition, the current errors as shown in Figure 12d,f are significantly smaller than those in Figure   Figure 11. Transient response with current magnitude jumping from 1 A to 4 A at 0.2 s: (a) MPCC without modulation; (b) the proposed MMPCC.

Load-Torque and Speed Varying Condition
Further evaluation is conducted at various load-torque and speed conditions. To better illustrate the difference between the two current controllers in terms of their performances under steady-state current tracking, the waveforms of α-β current errors are provided.
Given in Figures 12 and 13 are the results of IPMSM running at a constant speed of 500 rpm under different load-torques of 1 N-m and 2 N-m, respectively. Figure 12a shows the corresponding current waveforms in which the MPCC generates 0.2463 A of current ripples and 5.8826% of THD. In that regard, the proposed MMPCC in yields smaller ripples and THD with high percentage performance improvements up to 44.05% and 62.29%, respectively, as depicted in Figure 12b. In addition, the current errors as shown in Figure 12d,f are significantly smaller than those in Figure 12c,e, respectively. Figure 13 illustrates the current responses of the two control schemes with the IPMSM running at 500 rpm with a load-torque of 2 N-m. The current ripple and THD of the MPCC are reduced by 34.84% and 3.24%, respectively, by the proposed MMPCC, according to the data. Figure 14 shows the results under a constant speed of 1000 rpm and 1 N-m of load-torque. As can be seen, the THD index of the proposed MMPCC is much more satisfactory. It is found that the current ripple and THD using the proposed controller are reduced by 21.49% and 1.97%, respectively, compared to the MPCC. Similar results can also be observed in Figure 15. Based on Figures 12-15, both current controllers can perform the current tracking tasks for the IPMSM operating at 200, 500, and 1000 rpm successfully. However, the proposed MMPCC has much better current-tracking performance in significantly reducing the current errors. 12c,e, respectively. Figure 13 illustrates the current responses of the two control schemes with the IPMSM running at 500 rpm with a load-torque of 2 N-m. The current ripple and THD of the MPCC are reduced by 34.84% and 3.24%, respectively, by the proposed MMPCC, according to the data. Figure 14 shows the results under a constant speed of 1000 rpm and 1 N-m of load-torque. As can be seen, the THD index of the proposed MMPCC is much more satisfactory. It is found that the current ripple and THD using the proposed controller are reduced by 21.49% and 1.97%, respectively, compared to the MPCC. Similar results can also be observed in Figure 15. Based on Figures 12-15, both current controllers can perform the current tracking tasks for the IPMSM operating at 200, 500, and 1000 rpm successfully. However, the proposed MMPCC has much better current-tracking performance in significantly reducing the current errors.

Anaylsis
The experimental results are summarized in Table 5. For the steady-state tests, Figures 8 and 9 depict that the current ripples are reduced by approximately 34.54% and 32.92%, respectively. On the other hand, both the MMPCC and the MPCC have good transient responses, as revealed by Figures  10 and 11. Figure 12 and Table 5 show an interesting finding that the current ripple and THD using the MMPCC are reduced by 44.05% and 62.29%, respectively, compared to the MPCC without modulation. In addition, in Figures 13-15, one can observe that the reductions in current ripples are up to 34.84%, 21.49% and 38.78%, respectively. Based on Table 5 one may conclude that the proposed MMPCC largely enhances the current-tracking performance of the IPMSM drive system.

Anaylsis
The experimental results are summarized in Table 5. For the steady-state tests, Figures 8 and 9 depict that the current ripples are reduced by approximately 34.54% and 32.92%, respectively. On the other hand, both the MMPCC and the MPCC have good transient responses, as revealed by Figures 10  and 11. Figure 12 and Table 5 show an interesting finding that the current ripple and THD using the MMPCC are reduced by 44.05% and 62.29%, respectively, compared to the MPCC without modulation. In addition, in Figures 13-15, one can observe that the reductions in current ripples are up to 34.84%, 21.49% and 38.78%, respectively. Based on Table 5 one may conclude that the proposed MMPCC largely enhances the current-tracking performance of the IPMSM drive system.

Conclusions
An MMPCC for IPMSM is proposed and successfully implemented in this work. The new method modulates two voltage vectors in a sampling period with adaptive application durations and achieves current predictions effectively. Compared to the conventional MPC using seven cost values, the proposed MMPCC enlarges the candidate pool by employing thirteen cost values, while retaining simplicity. Through the novel design of a cost function, an optimal switching mode is obtained whose application time is optimized and calculated online. The presented experimental results demonstrate that the proposed method has good current tracking capabilities under steady-state, transient, and various operation conditions. Quantitatively, the improvements it made to the MPCC include the current ripple reduction by 27.17% and THD reduction by 21.84%, respectively, by means of averaging all data from the eight experiments, among which the highest current ripple reduction rate achieved was up to 44.05% in the best case. Thanks to the modulation and online optimization technique, the proposed MMPCC largely contributes in the reduction of current ripples and improves current tracking errors, compared to the MPCC without modulation, hence achieving significant performance improvement.