Have a personal or library account? Click to login
Model Predictive Control of a New Low Cost 31-Level Inverter Cover

Model Predictive Control of a New Low Cost 31-Level Inverter

Open Access
|Oct 2025

Full Article

1.
Introduction

The development of DC–AC power inverters has progressed through multiple stages. There are three main factors that contribute to the enhancement of inverter performance: first, the growing demand for advanced high-performance applications; second, the development of semiconductor devices capable of handling larger amounts of power while minimising switching time and power losses; and third, the availability of high-performance low-cost digital controllers which enables the implementation of computationally demanding control algorithms.

The ability of a multilevel inverter (MLI) to closely resemble the desired sinusoidal output voltage has stimulated the tendency to construct MLI with a large number of levels. However, increasing the number of levels beyond a certain point becomes impractical in the basic MLI topologies due to the excessively large number of elements, which in turn increases complexity and cost while reducing reliability (Adly et al., 2023).

Several innovative topologies have been proposed building on the three basic topologies: Neutral-Point-Clamped (NPC) (Guennegues et al., 2009), Flying Capacitor (FC) (McGrath and Holmes, 2009) and Cascaded H-bridge (CHB) (Rodriguez et al., 2007). One of the earliest improvements, the asymmetrical CHB, is capable of achieving up to 3n levels for n cascaded cells rather than the 1 + 2n levels from the basic design. The maximum number of levels is obtained when the cascaded cells are supplied by DC sources with a ratio of 3 (Abdul Kadir and Hussien, 2004). However, this approach requires multiple active supplies with bidirectional current capability. To reduce the number of costly and bulky DC supplies, FCs are employed to provide additional DC links. However, the capacitor voltage must be maintained close to its average value. To achieve this, most modulation approaches exploit the availability of redundant states that have opposite effects on the capacitor voltage (Chen et al., 2020; Ramu et al., 2025).

New types of MLIs have been created, including the 31-level inverter, which has better output voltage quality and lower total harmonic distortion (THD) (Rodriguez et al., 2007). Adding more voltage levels improves the waveforms and enhances the whole system’s efficiency. This is why high-level MLIs are attractive for high-power and precise applications.

Several designs that deviate from the three basic topologies have been proposed with the aim of achieving the maximum number of voltage levels for a given number of elements and consequently, lower circuit complexity. The concept of nested uneven bridges has been employed to maximise the number of levels per switching device (Nanda et al., 2022). Although the number of levels increases exponentially with the number of nested bridges, this design lacks modularity.

Advanced control strategies are necessary to fully exploit the benefits of MLIs. Model predictive control (MPC) has emerged as a promising method that reduces harmonic distortion, improves DC-link voltage balance and provides fast dynamic response (Riar et al., 2015). MPC improves control performance by predicting the system future behaviour and changing the control signals accordingly. This is especially helpful for complicated MLI systems (Harbi et al., 2023).

Fixed-frequency MPC has been applied to reduce voltage ripples in FCs. For example, in a 5-level FC converter, this approach improves performance by ensuring accurate regulation of both FC and output voltages (Yang et al., 2024). However, this study does not address innovative inverter designs or other control targets, such as reducing switching losses and common mode voltage.

Finite Control Set-Model Predictive Control (FCS-MPC) has also been applied to optimise DC capacitor voltages in grid-connected and stand-alone inverters. For example, a nine-level Packed U-Cell inverter uses MPC to balance capacitor voltages while minimising THD and reducing capacitor voltage ripple. This control method allows the use of smaller DC capacitors in a nine-level Packed U-Cell inverter, making it suitable for various industrial applications, particularly for single-phase grid-connected devices (Noghani et al., 2019).

Karamanakos et al. (2018) demonstrated that the choice of error norm in the cost function significantly influences the performance of closed-loop systems in MPC with reference tracking. This study emphasises that the ℓ1-norm, which has been used in some works to simplify the cost function, leads to degraded performance and potential closed-loop instability. By contrast, the ℓ2-norm overcomes these limitations, providing improved stability. A case study based on a three-level inverter drive system highlighted the advantages of adopting the ℓ2-norm.

Hardware-in-the-loop (HIL) studies have been conducted to validate the real-time performance of MPC strategies in grid-connected converters. For instance, a simple predictive direct power control method has been tested through HIL studies to show its efficiency in reducing power ripple and Total Harmonic Distortion (THD) (Quang et al., 2024). The method was applied to a five-level Active Neutral Point Clamped (ANPC) inverter; however, its extension to innovative topologies is not straightforward.

Due to the complexities and ambiguities associated with tuning the individual cost function weighting parameters in FCS-MPC, some studies have considered sequential model predictive control (S-MPC) as an alternative. In this approach, a sequence of single-objective cost functions, ordered by their relative importance, is evaluated. The primary cost function generates a shortlist of N viable switching states, which are then further assessed by the secondary cost function. Ultimately, a switching vector that satisfies the constraints of both cost functions is selected and applied to the inverter. Bonaldo et al. (2025) conducted a comprehensive evaluation comparing S-MPC and FCS-MPC for a grid-connected NPC inverter and recommended S-MPC as a more robust alternative.

Another strategy involves a two-step optimisation process to balance FCs and neutral point (NP) voltages. In this approach, the first step focuses on FC voltage balancing, while the second step addresses NP voltage control. This method not only mitigates voltage fluctuations across the entire frequency range but also reduces the computational burden by dividing the optimisation process into manageable steps (Noori et al., 2023). However, this study did not discuss implementation challenges, and the considered nested NPC inverter does not optimise the inverter size.

In the present study, a single-phase 31-level inverter topology with a reduced component count is proposed, based on the concept presented by Omer et al. (2020). The present study aims to develop an effective voltage control strategy for this circuit. While the original inverter requires four isolated DC supplies, the proposed design uses merely two supplies. The other two supplies are replaced by two capacitors. A FCS-MPC strategy has been developed with three targets: first, to track the target output current with minimum error; second, to maintain the capacitor voltages near their nominal values; and third, to minimise the switching losses, considering that the switching devices are subjected to different voltage levels.

The paper is organised as follows: Section 2 introduces the proposed inverter circuit and defines its switching variables and switching states. Section 3 outlines the FCS-MPC control strategies and the associated computational methods. Section 4 presents the inverter model and simulation results. Section 5 discusses and compares the results with other studies. Section 6 presents the conclusions.

2.
Proposed MLI Topology and Switching States

The proposed inverter design, which features a reduced component count, is shown in Figure 1. The existing equivalent 31-level inverter has four isolated DC voltage supplies related by the ratio 1:2:5:10 (Omer et al., 2020). The first contribution of this work is the replacement of two of the four isolated DC supplies with capacitors. The capacitor voltages are E and 2E, where E is the voltage step and the maximum inverter output voltage is 15E. The 10 switching devices forming the inverter circuit, depicted in Figure 1 as MOSFETs, are grouped into five complementary pairs, that is, the state of switch Si is the opposite of switch Si′, where i ∈ {1, 2, 3, 4, 5}. The switching variable xi determines the state of switch Si, with 0 and 1 representing the OFF and ON states, respectively. The five-digit binary number (x5 x4 x3 x2 x1) represents the inverter switching state.

Figure 1.

The proposed single-phase 31-level inverter circuit.

From Figure 1, the output voltage is expressed as follows: (1) vo=vAA+vAB+vBB {v_o} = {v_{AA'}} + {v_{A'B'}} + {v_{B'B}}

The relationships between the switching state and the output voltage, as well as between the inverter switching state and the capacitor voltage variation, can be deduced by examining Figure 2. In Figure 2(a), the production of the voltage VAA is explained. This voltage is determined by the switching variables x3 and x1. The relationship can be summarised as follows: (2) vAA=4x3+x1E {v_{AA'}} = \left( {4{x_3} + {x_1}} \right)E

Figure 2.

The operation modes of the inverter system illustrate how the switching variables control the voltages: (a) The effect of (x3, x1) combinations on VAA; (b) The effect of (x4, x2) combinations on VB′B; (c) The effect of x5 on VA′B′.

In Figure 2(b), the production of the voltage VBB′ is shown. This voltage is determined by the switching variables x4 and x2. The relationship can be summarised as follows: (3) vBB=8x4+2x2E {v_{B'B}} = \left( {8{x_4} + 2{x_2}} \right)E

Figure 2(c) shows that the switching variable x5 determines VA′B′ as follows: (4) vAB=15Ex51 {v_{A'B'}} = 15E\left( {{x_5} - 1} \right)

Substituting Eqs. (2)–(4) into Eq. (1), yields the output voltage as follows: (5) voE=i=152i1xi16+x5¯ {{{v_o}} \over E} = \left( {\sum\nolimits_{i = 1}^5 {{2^{i - 1}}{x_i}} } \right) - 16 + \overline {{x_5}}

Figure 2 also shows the variation of the two capacitors voltages. These voltages depend on the load current (io) and the inverter switching state. The variation of the voltages of C1 and C2 over one switching interval (Ts) can be obtained by examining Figures 2(a) and 2(b), respectively: (6) δvC1=x3x1ioTsC1 \delta {v_{{C_1}}} = \left( {{x_3} - {x_1}} \right){{{i_o}{T_s}} \over {{C_1}}} (7) δvC0=x4x2ioTsC2 \delta {v_{{C_0}}} = \left( {{x_4} - {x_2}} \right){{{i_o}{T_s}} \over {{C_2}}} where δvC1 and δvC2 are the voltage variations over one switching interval, C1 and C2 represent the capacitance of the two capacitors and io is the load current.

3.
Controlling the Inverter by FCS-MPC Algorithm

This section presents the development of the control algorithm to track a reference current, stabilise the FCs voltages and minimise switching losses. The MPC approach has been selected due to its suitability for multi-target control. FCS-MPC has been shown to be one of the effective tools in converter control, particularly optimising inverter performance with multiple performance objectives (Muhaisen and Abdul Kadir, 2021). MPC is a model-based control strategy that uses a discrete-time model of the system to predict future behaviour. The core idea is to select the optimal switching state that minimises a predefined cost function. For inverters, the cost function usually includes factors related to current tracking error, common mode voltage suppression and capacitor voltage balancing (Garcia et al., 2024; Le et al., 2024; Nitheesh et al., 2022).

The cornerstones of an FCS-MPC algorithm are as follows:

  • Discrete-time model: The system’s future behaviour is predicted using a discrete-time model. In MLIs, the model predicts the current for all possible voltage vectors. The conventional MPC predicts the behaviour for one step ahead, while multiple step prediction horizons enhance the performance at the cost of exponential growth of computational burden.

  • Cost function: The cost function is designed to prioritise control objectives. For instance, minimising the current error and common mode voltage (CMV) are common objectives in three-phase inverters.

  • Voltage vector selection: The algorithm evaluates all possible voltage vectors and selects the one that minimises the cost function. However, the large number of voltage vectors in MLIs makes this process computationally intensive.

The following three subsections describe the three components of the discrete time model, and the fourth subsection describes the control system and the voltage vector selection.

3.1.
Reference current tracking

Considering an RL-load, the next desired output current has been forecasted based on the past values of the reference current and using a cubic spline extrapolation function. In the literature, the most common method to predict the reference output current is the Euler-based linear extrapolation (Sanca et al., 2023). More accurate estimation has been achieved with the higher order Runge-Kutta-based extrapolation, as presented by Adjei-Saforo et al. (2022). In this work, it has been noticed that the cubic spline method leads to more accurate prediction. Figure 3 shows a comparison of the reference current waveform with the prediction by Euler-based linear extrapolation and cubic spline interpolation. At 50 Hz reference frequency and 0.25 ms sampling time, the RMS error obtained using cubic spline interpolation is <5% of the error resultant using linear extrapolation.

Figure 3.

Comparison of the reference current with currents predicted by linear and cubic spline extrapolation.

The predicted reference current is determined as a function of the current value and the three most recent reference currents: (8) io,refk+1=fio,refk,io,refk1,io,refk2,io,refk3 {i_{o,ref}^{k + 1}}' = f\left( {i_{o,ref}^k,i_{o,ref}^{k - 1},i_{o,ref}^{k - 2},i_{o,ref}^{k - 3}} \right) where jo,refk+1 {j_{o,ref}^{k + 1}}' is the predicted next reference current. The past four current values are used to determine the parameters of the third-order extrapolation polynomial that predicts the next current.

To reduce computational complexity, it has been decided to identify the voltage magnitude required to realise the reference current instead of comparing it with the current produced by each inverter’s next switching state. This approach reduces computational effort, since comparing voltage levels for all inverter states is more efficient than calculating the resulting current for every possible switching state.

The (k + 1)th reference voltage is determined using the measured kth and estimated (k + 1)th reference current values, as follows: (9) vo,refk+1=Lio,refk+1iokTs+io,refk+1R {v_{o,ref}}^{k + 1} = L{{\left[ {{i_{o,ref}}^{k + 1} - {i_o}^k} \right]} \over {Ts}} + {i_{o,ref}}^{k + 1}R where L and R are the load inductance and resistance, respectively, and iok i_o^k is the measured output current at the kth step. For extended prediction horizons beyond one step, the (k + 2)th reference voltage is determined by the following equation: (10) vo,refk+2=Lio,refk+2iok2Ts+io,refk+2R {v_{o,ref}}^{k + 2} = L{{\left[ {{i_{o,ref}}^{k + 2} - {i_o}^k} \right]} \over {2Ts}} + {i_{o,ref}}^{k + 2}R where vo,refk+2 and io,refk+2 are the reference voltage and current at the second step ahead, respectively.

The cost function term corresponding to the required voltage to track a given reference current for the next switching state (i) and subsequent state (j) is given in Eq. (11): (11) Cvi,j=Kv1Wssvo,refk+1voiE2+Wssvo,refk+2vojE2 {C_v}\left( {i,j} \right) = {K_v}\left( {\left( {1 - {W_{ss}}} \right){{\left( {{{{v_{o,ref}}^{k + 1} - {v_o}^i} \over E}} \right)}^2} + {W_{ss}}{{\left( {{{{v_{o,ref}}^{k + 2} - {v_o}^j} \over E}} \right)}^2}} \right) where Cv is the cost function term corresponding to the voltage error, Kv is a weighting factor for reference voltage with 0 < Kv < 1 and WSS is a factor that represents the significance of the step (k + 2) compared with the following (k + 1) step, where 0 ≤ WSS < 0.5; WSS = 0 for a one-step prediction horizon. i and j are next and subsequent states indices, respectively, 0 ≤ i, j ≤ 31.

It must be emphasised that Eqs (9) and (10) describe an RL load; other types of loads such as machines, non-linear loads and grid-connected inverters have different models, many of which have been reported in the literature.

3.2.
Capacitor voltage balancing

According to the effect of the switching state on the capacitor voltages, given in Eqs (6) and (7), the cost function term representing the capacitor for a one-step prediction horizon is given in Eq. (12): (12) Cci=KcEVC1k+δVC1i2+2EVC2k+δVC2i2 {C_c}\left( i \right) = {K_c}\left[ {{{\left( {E - {V_{C1}}^k + \delta {V_{C1}}\left( i \right)} \right)}^2} + {{\left( {2E - {V_{C2}}^k + \delta {V_{C2}}\left( i \right)} \right)}^2}} \right] where Cc is the term of the cost function corresponding to the capacitor voltage, and Kc is the weighting factor for capacitors balancing, with 0 < Kc < (1 − Kv).

Eq. (12) can also be extended to two-step prediction horizon as follows: (13) Cci,j=Kc1Wss×EVC1k+δVC1i2+2EVC2k+δVC2i2+WssEVC1k+δVC1j2+2EVC2k+δVC2j2 {C_c}\left( {i,j} \right) = {K_c}\left[ {\matrix{ {\left( {1 - {W_{ss}}} \right) \times \left( {\matrix{ {{{\left( {E - {V_{C1}}^k + \delta {V_{C1}}\left( i \right)} \right)}^2} + } \cr {{{\left( {2E - {V_{C2}}^k + \delta {V_{C2}}\left( i \right)} \right)}^2}} \cr } } \right)} \cr { + {W_{ss}}\left( {\matrix{ {{{\left( {E - {V_{C1}}^k + \delta {V_{C1}}\left( j \right)} \right)}^2} + } \cr {{{\left( {2E - {V_{C2}}^k + \delta {V_{C2}}\left( j \right)} \right)}^2}} \cr } } \right)} \cr } } \right]

3.3.
Switching losses minimisation

As shown in Figure 2, switches S1 and S1′ are subjected to an off-state voltage of E, while S2 and S2′ experience 2E. It can also be seen that S3 and S3′, S4 and S4′ and S5 and S5′ are subjected to 4E, 8E and 15E, respectively. Since the off-state voltage of a switch influences the switching losses, the weighting factors in the switching-loss cost term are included to account for the different voltage levels of the switches.

The switching losses are determined based on the variation of the switching state and the voltage across each switch, as follows: (14) Cswi=Kswx1kx1i2+2x2kx2i2+4x3kx3i2+8x4kx4i2+15x5kx5i2 {C_{sw}}\left( i \right) = {K_{sw}}\left[ {\matrix{ {{{\left( {x_1^k - x_1^i} \right)}^2} + 2{{\left( {x_2^k - x_2^i} \right)}^2} + 4{{\left( {x_3^k - x_3^i} \right)}^2} + } \cr {8{{\left( {x_4^k - x_4^i} \right)}^2} + 15{{\left( {x_5^k - x_5^i} \right)}^2}} \cr } } \right]

Csw is the cost function term corresponding to the switching losses, and xyi x_y^i is a binary digit that represents the state of switch y in state i. Ksw is the weighting factor corresponding to the switching losses given by the following equation: (15) Ksw=1KvKc {K_{sw}} = 1 - {K_v} - {K_c}

Each of the five terms in the brackets in Eq. (14) takes a value of 1 if the corresponding switching device changes its state during the transition from the present state (state k) to state i, and 0 if it maintains its current state. When these terms are multiplied by the integers 1, 2, 4, 8, and 15, the result represents the ratios of the off-state voltages across each pair of switches. Since the converter operates in hard-switching mode, the switching losses are directly proportional to the off-state voltages. Therefore, the switching cost term is represented accordingly.

Similar to Eqs (11) and (13), Eq. (16) gives the switching loss cost for two-step prediction horizon: (16) Cswi=Ksw1Wssx1kx1i2+2x2kx2i2+4x3kx3i2+8x4kx4i2+15x5kx5i2+Wssx1kx1j2+2x2kx2j2+4x3kx3j2+8x4kx4j2+15x5kx5j2 \matrix{ {{C_{sw}}\left( i \right) = {K_{sw}}\left[ {\left( {1 - {W_{ss}}} \right)\left[ {{{\left( {x_1^k - x_1^i} \right)}^2} + 2{{\left( {x_2^k - x_2^i} \right)}^2} + 4{{\left( {x_3^k - x_3^i} \right)}^2} + 8{{\left( {x_4^k - x_4^i} \right)}^2} + 15{{\left( {x_5^k - x_5^i} \right)}^2}} \right] + {W_{ss}}} \right]} \hfill \cr {\left[ {{{\left( {x_1^k - x_1^j} \right)}^2} + 2{{\left( {x_2^k - x_2^j} \right)}^2} + 4{{\left( {x_3^k - x_3^j} \right)}^2} + 8{{\left( {x_4^k - x_4^j} \right)}^2} + 15{{\left( {x_5^k - x_5^j} \right)}^2}} \right]} \hfill \cr }

3.4.
The FCS-MPC controller

The flowchart of the FCS-MPC one- and two-step prediction horizon function is shown in Figure 4. The function begins by calculating the predicted values of the next reference currents and then the amplitude of the reference voltage. The iteration loops are then started by assigning large minimum cost which will be updated when a smaller cost is found. Then the outer loop (index i) is associated with the next step while, for two-step horizon, the inner loop (index j) is associated with the second next step.

Figure 4.

Flow chart of the FCS-MPC control function. FCS-MPC, finite control set-model predictive control.

In each iteration, the cost function is calculated, the index is recorded and the minimum cost value is updated when the calculated cost is smaller than the recorded minimum.

4.
Simulation and Results

The inverter model has been constructed on the Simulink platform, considering the converter parameters given in Table 1. The simulation model layout is shown in Figure 5. The MPC algorithm is implemented using an m-code script within a ‘User-Defined Function’ block of Simulink at the bottom of the model.

Figure 5.

Simulation model of the FCS-MPC controlled 31-level inverter. FCS-MPC, finite control set-model predictive control.

Table 1.

Circuit and 31-level inverter parameters.

ParameterSymbolValue
Voltage stepE100 V
The loadR, L100 Ω, L = 0.2 H
Reference current io* i_o^* 12 sin (100 πt)
Inverter capacitorsC1, C2100 µF
MOSFETRon, Roff0.1 Ω, 100 kΩ
Sampling timeTs0.5 ms
4.1.
Model and parameters

An initial investigation has been conducted with a one-step prediction horizon. The effect of changing the weighting factor on various performance indicators is shown in Figure 6. Figures 6(a)–(d) show the average and peak-to-peak ripple voltages of the two capacitors. It can be seen that increasing the capacitor voltage weighting factor (Kc) brings the average capacitors voltage closer to their nominal values. Also, a higher Kc reduced the ripple of both capacitor voltages. However, increasing Kc leads to higher switching losses as shown in Figure 6(e), which also indicates that increasing the reference tracking weighting factor (Kv) has a similar effect on the switching losses. Lowering both Kv and Kc allows more weight for the switching loss term, since the sum of the three weights is set to 1; in that case lower switching losses can be achieved at the expense of higher current distortion and capacitor voltage ripple.

Figure 6.

Variation of capacitor voltage switching losses and current distortion with the weighting factors with one-step prediction horizon. THD, total harmonic distortion.

4.2.
One-step prediction horizon

The weighting factors have been tuned to optimise the converter performance considering the following performance indicator (PI): (17) PI=10×THD+1500ΣVswfsw+io,refio2dtT+0.1vC1100+0.05vC2200 PI = 10 \times THD + {1 \over {500}}\Sigma {V_{sw}}{f_{sw}} + \sqrt {{{\smallint {{\left( {{i_{o,ref}} - {i_o}} \right)}^2}dt} \over T}} + 0.1\left| {{v_{C1}} - 100} \right| + 0.05\left| {{v_{C2}} - 200} \right| where the constants 10, 1/500, 0.1 and 0.05 are used as normalisation factors and the sign Σ denotes the summation for all switching devices, Vsw and fsw are the switch voltage and switching frequency and T is the reference current period (20 ms for 50 Hz reference current). The third term on the right-hand side represents the RMS of the current error.

The optimum points of the weighting factors were obtained using the MATLAB function fmin, which searches for the values that minimise the cost function defined in Eq. (17). This approach ensures that the selected weighting factors provide the best trade-off between current tracking performance and capacitor voltage regulation. The results from fmin directly determine the optimal combination of weighting factors, which are then applied in the predictive control strategy to achieve improved system performance. The optimised operating point was achieved with Kv = 0.7 and Kc = 0.22. The results are summarised in Table 2.

Table 2.

Optimum weighting factors and the corresponding performance indicator with one-step prediction.

VariableKvKcTHD (io)Vsw fsw io,refio2dt/T \sqrt {\left( {\smallint {{\left( {{i_{o,ref}} - {i_o}} \right)}^2}dt} \right)/T} Vc1,dcVc2,dc
Value0.70.220.009652200.01643100.1 V201.2 V

THD, total harmonic distortion.

For the same values of Kv and Kc, the performance of the inverter when starting with zero initial capacitor voltages is depicted in Figure 5. The resultant current is evidently very close to the reference current with RMS deviation from reference value of 0.01643 A as shown in Table 2, the THD of the current is <1%. Figure 7 shows that the capacitor voltages reach the target voltage in about 3 ms, without the need for any dedicated starting protocol. This reflects the stable and fast dynamic response of the FCS-MPC algorithm. The average capacitor voltage deviation is <1% of the average value, and the occasional ripple in the capacitor voltage is <3% of the reference.

Figure 7.

The 31-level inverter operation under the one-step prediction horizon MPC algorithm. The waveforms shown are reference and output current (top), the two capacitor voltages (middle) and the switching state (bottom). MPC, model predictive control.

Table 2 reveals that the sum of all switching signal pulses multiplied by the normalised switch off-state voltage is 220. This value corresponds either to two switching devices subjected to voltage E and operating at 220 Hz, or to 15 transitions of the switches subjected to the higher voltage (15E). Thanks to the large number of levels, the algorithm tendency to shift the switching towards the switching devices is subjected to lower voltage stress.

4.3.
Two-step prediction horizon

The results corresponding to the two-step horizon are presented in Table 3, and the corresponding voltage waveforms are given in Figure 8. No significant improvement has been achieved by extending the prediction horizon to two steps despite the added computational complexity. The two-step prediction requires 31 × 31 evaluations of the cost function at every step compared with 31 calculations for a one-step horizon. On the contrary, the performance indicator defined in Eq. (15) has improved by <5%.

Figure 8.

The 31-level inverter operation under the two-step prediction horizon MPC algorithm. The waveforms shown are reference and output current (top), the two capacitor voltages (middle) and the switching state (bottom). MPC, model predictive control.

Table 3.

Optimum weighting factors and the corresponding performance with two-step prediction.

VariableWSSKvKcTHD (io)Vsw fsw io,refio2dt/T \sqrt {\left( {\smallint {{\left( {{i_{o,ref}} - {i_o}} \right)}^2}dt} \right)/T} Vc1,dcVc2,dc
Value0.250.70.220.009732230.0152100.2 V201.1 V

THD, total harmonic distortion.

5.
Discussion and Comparison to Related Work

To evaluate the advantages of the proposed FCS-MPC-controlled MLI system, a comparison is presented in Table 4 with other published research that introduces innovative MLI topologies of a comparable number of levels. The comparison takes into consideration several factors as discussed in this section.

Table 4.

Comparison of the proposed 31 level inverter to other comparable designs.

[Ref]Saforo et al. (2020)Arif et al. (2021)Panda et al. (2021)Memon et al. (2024)Khasim and Dhanamjayulu (2022)Radhakrishnan et al. (2024)Proposed
NLevels17171713332531
Ndc2141442
Ndc/NLevels0.1180.0590.2350.0770.12120.160.064
Ns12121012161310
NswNlevel {{{N_{sw}}} \over {{N_{level}}}} 0.7060.7060.5880.9230.4850.520.322
Ncap4303002
TSV: ∑ Voff,sw/Vo.max, p3.3754.5542.9-2.5
TCV: ∑ Vcap/Vo.max, p10.87500.833000.2
(Max no. of switches in current path)/levels4/17 = 0.2346/17 = 0.3535/17 = 0.2945/13 = 0.3857/33 = 0.2126/25 = 0.245/31 = 0.161

TSV, total switches voltage. TCV, total capacitors voltage.

The proposed design achieves the second-lowest number of DC supplies required, Ndc, and the ratio of number of supplies to the number of levels (Ndc/levels). The lowest ratio is achieved by the converter presented by Arif et al. (2021), which uses a switched capacitors arrangement to generate multiple dc voltage levels. However, the switched capacitor converters require high switching frequency, larger capacitors and deliver lower output power which may restrict their applications.

As shown in Table 4, the proposed design has the lowest number of switching devices, Nsw, per level compared with the other six designs. This ratio is approximately 50% lower than that of the second-best design reported by Khasim and Dhanamjayulu (2022).

To obtain a higher number of voltage levels, some MLI designs use capacitors as voltage dividers, FCs or switched capacitors. Capacitor-less MLIs require a larger number of the more expensive, bulky isolated DC supplies. Compared with other designs using capacitors, the proposed design requires the fewest number of capacitors despite its large number of levels.

Total switches voltage (TSV) is calculated as the sum of the voltage stress across all the switches. The proposed MLI considerably reduces the TSV compared with all six other designs. The total capacitors voltage (TCV) is 24% or less than TCV of the three other designs that use capacitors. This factor is important as it affects the converter cost and size.

Finally, the proposed MLI outperforms the other six designs in the number of switching devices in the conduction path per level. This parameter influences the conduction losses of the converter. The proposed design has the second-lowest number of switching devices in the current path after the 17 level inverter presented by Saforo et al. (2020).

6.
Conclusion

This paper introduces a single-phase 31-level inverter that achieves reduced cost compared with other MLIs offering a similar number of levels. The proposed design uses 2 DC sources, 10 switching devices and 2 low voltage capacitors. It offers low switching losses due to fewer switching events per cycle and the number of switching devices in the current path is the smallest per level among six other recent designs reported in the literature.

This study highlights the significance of closed-loop MPC as an advanced control strategy for MLIs. By applying FCS-MPC to the proposed 31-level inverter, predictive optimisation can improve system performance while overcoming the constraints of traditional control techniques. This demonstrates the potential of MPC to play a central role in next-generation power conversion systems, particularly for renewable energy and high-power industrial applications.

Performance analysis shows that the proposed inverter achieves very low THD and maintains excellent voltage balance of the FCs, even at relatively low switching frequencies. This outstanding performance is enabled by the FCS-MPC algorithm, which not only satisfies multiple control objectives but also provides the flexibility to dynamically prioritise application-specific goals.

The results in this paper are based on modelling and simulation. Experimental validation on a hardware prototype is planned as future work to confirm the practical applicability of the proposed approach.

DOI: https://doi.org/10.2478/pead-2025-0023 | Journal eISSN: 2543-4292 | Journal ISSN: 2451-0262
Language: English
Page range: 342 - 356
Submitted on: Jul 13, 2025
|
Accepted on: Sep 30, 2025
|
Published on: Oct 31, 2025
In partnership with: Paradigm Publishing Services
Publication frequency: 1 issue per year

© 2025 Mohamad N. Abdul Kadir, Harith Al-Badrani, Yasir M. Y. Ameen, published by Wroclaw University of Science and Technology
This work is licensed under the Creative Commons Attribution 4.0 License.