Harmonic Analysis and Attenuation Strategy for a Two-Stage Matrix Converter Fed by Dual-Inverter Based on Pulse Barycenter Method

: Since the DC-link of the dual-inverter two-stage matrix converter (TSMC) has no energy storage element, its grid-side current waveform is closely related to the load-side condition. The conventional modulation strategy of the dual-inverter TSMC is only the stack of the double inverters, without considering the cooperation between both modulation strategies. Thus, the ﬂuctuations in the grid-side currents synthesized by both load-side currents will be obvious. To solve this problem, this paper analyzes the effect of the current pulse position on the grid-side current waveform quality, and relevant relationship between the pulse position and current harmonics is obtained according to Fourier transform. On this basis, the harmonic attenuation strategy for the dual-inverter TSMC is proposed. In the proposed strategy, the grid-side current harmonic is reduced based on the optimal duty-cycle combination, obtained from the harmonic optimization problem. The simulation and experimental results are provided to verify that the theoretical analysis of grid-side current harmonic content is consistent with the actual condition, and the proposed modulation strategy can effectively reduce the grid-side current harmonic.


Introduction
Two-stage matrix converters (TSMCs) have wide application prospects thanks to the advantages of good input and output performance, invertible energy transmission, compact structure, flexible and adjustable grid-side power factor, and so on. The topology of TSMCs can be divided into the rectifier stage and inverter stage. There is no energy storage element in TSMCs, compared with AC-DC-AC converters, which reduce the volume and weight, improve the reliability, and prolong the life of the system [1][2][3]. TSMCs have great superiority in the field of aerospace and military equipment, where the requirement of the converter volume is strict.
Since the DC-link has no energy storage element, the inverter modulation strategies will affect the current waveform in the inverter stage, and then have an effect on the rectifier stage through the DC-link, which means the grid-side currents will be influenced by the load directly. Aiming at the harmonic problem of TSMC, the method of Fourier transform analysis is usually utilized. Specifically, the mathematical model about load-side voltage is established by combining the pulse width modulation (PWM) principle with the transfer function firstly, then analytical expression of each harmonic component is obtained by use of mathematic derivation and Fourier transform [4][5][6][7][8]. For example, reference [9] constructs the constraint equation according to the conduction angle in each modulation period and obtains the formula of the amplitude of harmonics by use of Fourier transform. Reference [10] proposes a harmonic calculation method with triple Fourier series and applies this method to calculate the output voltage harmonic of the Direct Matrix Converter (DMC) when the space vector PWM (SVPWM) strategy is used. Taking the "RMS current ripple" as the evaluating indicator, another common harmonic analysis method takes the integral of the errors between the expected current and actual current, applying this for inverter output waveform analysis [11,12]. For example, reference [12] proposes a current optimization strategy for the three-level inverter, in which the "RMS current ripple" is calculated under different modulation strategies, so the switching sequences according to the minimum "RMS current ripple" principle can be determined.
At present, the random PWM technology is usually adopted to reduce the modulation harmonics, i.e., for references [13,14], the grid-side current harmonics are suppressed by changing the carrier frequency or setting the distribution of the PWM voltage pulse within a certain range. Similarly, the SVPWM technology can also be improved to suppress the harmonics by correcting the duty cycles, such as adjusting the synthesizing sequence of voltage vectors, redistributing the duty cycles of active vectors, and so on [15,16]. Moreover, scholars also utilize different filters to suppress the grid-side current harmonics of MC, such as parallel active filter, hybrid filter, etc. [17][18][19][20][21][22][23].
As for the existence of the DC-link structure in TSMCs, the rectifier stage can connect two inverter stages; thus, TSMCs are suitable for operating conditions with a high requirement of power density, such as aircraft power systems, etc. [24,25]. When a DC-bus is linked with two parallel inverters, different modulation strategies need to be employed in double inverters. The currents of dual-inverters will be stacked and affect the grid-side current through the DC-link together, which makes the harmonic analysis of the grid-side current more complex. The harmonics of the gird-side current in dual-inverter TSMCs are related to the input frequency, the output frequency of double inverters, and the carrier frequency, etc. Therefore, the multiple-Fourier transform is needed to explain the coupling relationship in TSMCs, leading to a complex calculation process. In addition, the analysis method based on "RMS current ripple" needs to assume that the DC-bus voltage is constant, hence, it is not suitable for the harmonic analysis in topologies, such as TSMC. Besides, since the asymmetric modulation is normally utilized in TSMCs, there are nonlinear factors in this modulation method and it will lead to the distortion of the grid-side current waveform. Thus, both Fourier transform and the "RMS current ripple"-based method cannot be applied for harmonic analysis in dual-inverter TSMCs.
The method of Fourier transform occupies a dominant position in static harmonic analysis. Generally, fast Fourier transform (FFT) [26] or multiple-Fourier transform [27] are used to obtain the magnitude of fundamental and harmonic components of output voltage. However, the dynamic waveform quality evaluation methods have also been widely studied in recent years. For example, a dynamic harmonic evaluation method based on "effective value of output current fluctuation" is proposed in [28]. Firstly, the geometric relationship between the "effective value of output current fluctuation" and the triangle center of gravity of the output error current vector trajectory of the converter is established, and then the conclusion is derived that if the distance from the center of gravity to the origin is the shortest, the harmonic distortion rate of the converter is the smallest. Reference [29] proposes a dynamic harmonic evaluation method for a converter based on "effective value of stator flux fluctuation". Firstly, the flux error vector is decomposed in the d-q coordinate system, and the d-axis component and q-axis component of the "stator flux fluctuation vector" are solved, respectively. Then, the two components are combined to obtain the "effective value of stator flux fluctuation" generated in the unit sampling period, which is used to dynamically evaluate the harmonics in the converter. However, the existing harmonic analysis methods are limited by the heavy online calculation and cannot achieve dynamic harmonic evaluation. The pulse barycenter method (PBM) comprehensively considers multiple factors influencing the current harmonics, such as the pulse width, amplitude, and distribution, etc. [30][31][32][33]. On this basis, PBM can be used to analyze the grid-side current harmonics of TSMCs according to the pulse barycenter movement [34][35][36]. Specifically, the grid-side harmonic-equivalent-value-based PBM contains the information on the pulse location and the relationship between the grid-side current harmonics and the pulse barycenter can be effectively reflected by this value [36][37][38][39]. For a dual-inverter TSMC system, this paper employs PBM to analyze the effect of the current pulse distribution on the grid-side current harmonics in detail, then uses Fourier transformation to obtain the change rule of the pulse pattern with current harmonics; thus, a modulation strategy for dual-inverter TSMCs based on zero-vector duty-cycle correction (DCC) is proposed. The simulation and experimental results verify that the proposed modulation strategy according to PBM can effectively reduce the grid-side current harmonics and improve the quality of the grid-side current waveform.

Modulation Strategy of Dual-Inverter TSMCs
The topology of the double-inverter TSMC system is shown in Figure 1. From Figure 1, it can be seen that the DC-link has no energy storage element, and the DC-bus voltage is not the constant. Therefore, the operating of the rectifier stage needs to cooperate with the inverter stage, so as to implement the modulation of the whole system. contains the information on the pulse location and the relationship between the grid-side current harmonics and the pulse barycenter can be effectively reflected by this value [36][37][38][39]. For a dual-inverter TSMC system, this paper employs PBM to analyze the effect of the current pulse distribution on the grid-side current harmonics in detail, then uses Fourier transformation to obtain the change rule of the pulse pattern with current harmonics; thus, a modulation strategy for dual-inverter TSMCs based on zero-vector duty-cycle correction (DCC) is proposed. The simulation and experimental results verify that the proposed modulation strategy according to PBM can effectively reduce the grid-side current harmonics and improve the quality of the grid-side current waveform.

Modulation Strategy of Dual-Inverter TSMCs
The topology of the double-inverter TSMC system is shown in Figure 1. From Figure  1, it can be seen that the DC-link has no energy storage element, and the DC-bus voltage is not the constant. Therefore, the operating of the rectifier stage needs to cooperate with the inverter stage, so as to implement the modulation of the whole system.

The Modulation Strategy of Rectifier Stage
The grid-side voltage of the TSMC can be expressed as where, Uim and ωr represent the amplitude of the grid-side phase voltage and angular frequency, respectively.

The Modulation Strategy of Rectifier Stage
The grid-side voltage of the TSMC can be expressed as where, U im and ω r represent the amplitude of the grid-side phase voltage and angular frequency, respectively. There are nine switching states in the rectifier stage, and the corresponding grid-side current can be converted to the space vector form. Then, six active vectors and three zero vectors are obtained, where the six active vectors are marked as i 1 ∼ i 6 . The current complex plane is divided into six sectors, as shown in Figure 2a. There are eight switching states in the inverter stage, and the corresponding output voltage can also be converted to the space vector form. Then, six active vectors and two zero vectors are obtained, where the six active vectors are marked as u 1 ∼ u 6 . The voltage complex plane is divided into six sectors, as shown in Figure 2a.
There are nine switching states in the rectifier stage, and the correspon current can be converted to the space vector form. Then, six active vectors vectors are obtained, where the six active vectors are marked as 1 6 i i   . Th plex plane is divided into six sectors, as shown in Figure 2a. There are ei states in the inverter stage, and the corresponding output voltage can also to the space vector form. Then, six active vectors and two zero vectors are ob the six active vectors are marked as 1 6 u u    . The voltage complex plane i six sectors, as shown in Figure 2a. The grid-side reference current vector can be expressed as where, Ir and θrec represent the amplitude of the grid-side reference curre sector angle, respectively.
For the rectifier stage, the prior applied active vector in each control p The modulation index of the rectifier is defined as where, dc I is the average value of the DC-link current. In order to achiev voltage utilization ratio and reduce the times of commutation, zero vectors not applied in the modulation of the rectifier, and the modulation index m stant. Therefore, the duty ratios of the two active current vectors  1 i and pressed as [40] θ u θ u The grid-side reference current vector can be expressed as where, I r and θ rec represent the amplitude of the grid-side reference current vector and sector angle, respectively. For the rectifier stage, the prior applied active vector in each control period is set to ' → I µ ', and its duty cycle is d µ ; the second applied active vector is set to ' → I ν ', and its duty cycle is d ν . Taking sector I for example, the reference vector → I re f is synthesized by the active vector → i 1 and → i 2 . The formula of the reference vector in each PWM period is The modulation index of the rectifier is defined as where, I dc is the average value of the DC-link current. In order to achieve the maximal voltage utilization ratio and reduce the times of commutation, zero vectors are generally not applied in the modulation of the rectifier, and the modulation index m rec is kept constant.
Therefore, the duty ratios of the two active current vectors → i 1 and → i 2 can be expressed as [40]  

The Modulation Strategy of Inverter Stage
From the above analysis, it can be seen that the output current vector in the rectifier stage is synthesized by two active vectors. Thus, the DC voltage in the rectifier stage during each modulation period will be separated into two parts and then the modulation for the inverter stage needs to be implemented twice successively under these two segments of voltage. Meanwhile, in order to improve the safety and reliability of the system and reduce switching loss, the DC-link current is normally maintained at zero at the commutation instant in the rectifier stage, i.e., the voltage vectors of double inverters are zero.
The voltage vectors of the single inverter are shown in Figure 2b. The output voltage vector is set to → V re f , which is synthesized by two adjacent active vectors in each sector.
Taking sector I, for example, → V re f is assumed to be synthesized by where, θ inv is the sector angle of the output voltage vector; ω c is the output frequency of the inverter. The modulation index of the inverter is defined as where, U lm represents the amplitude of the load-side line voltage; θ in is the angle of grid-side voltage.
The distribution of the duty cycles of the dual-inverter TSMC system is shown in Figure 3.
the inverter stage needs to be implemented twice successively under these two seg of voltage. Meanwhile, in order to improve the safety and reliability of the syste reduce switching loss, the DC-link current is normally maintained at zero at the co tation instant in the rectifier stage, i.e., the voltage vectors of double inverters are z The voltage vectors of the single inverter are shown in Figure 2b.
where, θinv is the sector angle of the output voltage vector; ωc is the output freque the inverter.
The modulation index of the inverter is defined as where, Ulm represents the amplitude of the load-side line voltage; θin is the angle o side voltage. The distribution of the duty cycles of the dual-inverter TSMC system is sho Figure 3.  The duty cycles for the active vector Since the inverters in the system are not coupled, the modulation strategies of inverters are independent. Hence, the duty cycles for inverters can be calculated independently based on (8) and (9), and multiple combination patterns for the modulations can be constructed.

The Relationship between Current Pulse Barycenter and Load-Side Current
The current pulse barycenter (PB) of the single inverter could represent the distribution of the active vector pulse during each modulation period. The current waveform of Inverter 1 in each modulation period is shown in Figure 4. In Figure 4, i µ1 is the current of Inverter 1 when → I µ is applied; i ν1 is the current of Inverter 1 when → I ν is applied. kT s and (k + 1)T s are the beginning and ending instant of the kth period, respectively.
The time difference between the instant t µ1 and the beginning instant of the modulation period kT s is defined as current PB e µ1 (kT s ), where the instant of the geometric barycenter of i µ1 is achieved. e µ1 (kT s ) can be calculated by the duty cycles and control period T s .
i(t) is set to the average value of the inverter current in each period According to the area equivalence principle, the current i µ1 of Inverter 1 in the kth period can be approximately replaced by the ideal pulse with the amplitude T s i(kT s ) and delay time kT s + e µ1 (kT s ).
where, i(kT s ) is the average value of the inverter current during the kth period. Laplace transform is used to transform both sides of (11), and the formula of the current i µ1 (t) in the complex frequency domain is The infinite sum formula on the right side of Formula (12) can be further converted into integral form, then the frequency-domain expression of i µ1 (t) can be obtained as Define the variable τ as the function of t Substituting (14) into (13), the result is Formula (15) is the definition formula in the form of Laplace transform. Neglecting the high-frequency components, the current i µ1 (t) in time domain is From (16), it is known that the harmonic characteristics of i µ1 (t) and i(t) are similar, when e µ1 (t) is constant; the harmonics will generate in the differential block of (16), when e µ1 (t) is not constant. In other words, the PB of e µ1 (t) leads to the generation of harmonics in the current.

The Relationship between Current Pulse Barycenter and Grid-Side Current
Since the DC-bus of the dual-inverter TSMCs has no energy storage element, the grid-side current pulse can be synthesized by the composition of two inverter current pulses. Moreover, the change rule of the grid-side current PB is complex for the asymmetric modulation, so the duty cycles might be calculated firstly to obtain the PB of double inverters, and then the change rule of the grid-side current PB can be concluded. To sum up, (16) can be used to analyze the effect of the modulation strategy of inverters on the grid-side current harmonics.
Assuming that the reference vector in the rectifier stage is located in Section 2, the reference vector of Inverter 1 is located in Section 2, and the reference vector of Inverter 2 is located in Section 3, the distribution of vectors, grid-side current, and inverter currents are shown in Figure 5. The distribution coefficient of the zero vector in Inverter 1 is defined as x tion coefficient of the zero vector in Inverter 2 is defined as y, where x and y are value between 0 and 1. The distribution of the zero vector will change with th x and y, and the barycenter of the active vector pulse will also change accordin In Figure 5, the zero vector of Inverter 1 in each modulation period can be  Figure 5. Waveforms of vector distribution, grid-side current, and inverters' currents.
From Figure 2, it can be seen that the output line voltages in the rectifier are u ac and u bc . The DC-bus voltage is u ac when → i 2 is applied; the DC-bus voltage is u bc when → i 3 is applied. Moreover, the a-phase of the rectifier stage is turned on when u ac is applied, then i a will exist only when u ac is applied (the active vector → i 2 is applied). At this moment, i a is the sum of the current i µ1 of Inverter 1 and i µ2 of Inverter 2.
The distribution coefficient of the zero vector in Inverter 1 is defined as x; the distribution coefficient of the zero vector in Inverter 2 is defined as y, where x and y are the constant value between 0 and 1. The distribution of the zero vector will change with the variation in x and y, and the barycenter of the active vector pulse will also change accordingly.
In Figure 5, the zero vector of Inverter 1 in each modulation period can be expressed as Similarly, the zero vector of Inverter 2 in each control period can be expressed as From the definition of the current PB in Section 3.1, it is known that the PB of the current i µ1 and i µ2 of the inverter can be expressed as As the pulse waveform of i a is the sum of i µ1 and i µ2 , the PB of i a will also be determined by i µ1 and i µ2 . Considering the amplitude of i µ1 and i µ2 , the PB of i a can be expressed as The PB of i a in different sectors are different, and can be expressed as where, N r is the number of the sectors in the rectifier stage.
From (21), it is known that the PB of the grid-side current changes periodically, and it is related to the duty cycles in the rectifier stage and the inverter stage. Because of the duty cycles are the function of the input frequency ω r , output frequency ω c1, and ω c2 , the harmonics related to the input frequency and output frequency will exist in the grid-side current due to the variation in the current PB.

Modulation Strategy for Dual-Inverter TSMCs Based on Duty-Cycle Correction
From the analysis in Section 3.2, it is known that the current PB of the inverter changes with the variation in the value of the distribution coefficients x and y, and thereby, the current harmonic is affected. Therefore, the current harmonics can be suppressed by changing x and y. In order to calculate the grid-side current harmonic content under different values of x and y, the grid-side current harmonic content can be quantitatively calculated by combining with Fourier transformation with (15). The values of x and y, minimizing the grid-side current harmonic content, are obtained through the above procedure, and then the duty cycles of zero vectors for two inverters are corrected.
The Fourier coefficient represents the amplitude function of a certain frequency, and the harmonic content in the grid-side current can be calculated accordingly. The principle of the Fourier transformation is that waveform f (t) can be expressed by the form of infinite series, consisting of the sinusoidal harmonics with different orders, i.e., The amplitude of each harmonic a m can be calculated as By use of the transfer function method, the analysis result concluded that the function of the grid-side current is an odd function, so the input current i a (t) can be expressed as According to the analysis results in Section 3.2 and (16), the approximate expression of the gird-side current is obtained.
The average value in each period can be calculated in real time according to the grid-side current pulse waveform. In actual calculation, the ideal grid-side current can be expressed as i(t) = I im sin(ω r t) Substituting (21), (25), and (26) into (23), the amplitude of the m-th harmonic that changes with x and y can be obtained From (27), it can be known that the amplitude of the grid-side current harmonic is related to the current PB, and the current PB is related to the input frequency, output frequency, and the distribution coefficients x and y. Hence, the grid-side current harmonic content can be reduced by changing the values of x and y, and then the waveform quality of the grid-side current can be improved.
In order to analyze the changing rule of the grid-side current harmonic content expressed by (27), the following parameters listed in Table 1 are employed. Considering that the rectifier stage operates under zero-current commutation, the values of x and y vary from 0.2 to 0.8. For ease of calculation, only the percent of the 3rd, 5th, 7th, 9th, 11th, and 13th harmonics are calculated, and the harmonic orders higher than the 15th are neglected. The results of (27) are obtained by use of Maple17, and the results are shown in graphic form, as shown in Figure 6a. From Figure 6a, it can be seen that the harmonic content changes as the law of oloid. Based on this law, the distribution of the PWM pulse for the double invert be optimized, and the aim of decreasing the grid-side current harmonic content achieved.
In order to verify the results of the theoretical calculation, the values of x an selected and the a-phase grid-side current is analyzed by use of fast Fourier analy harmonic content is shown in Figure 6b, and it can be seen that the change trends in 6b are similar to Figure 6a.
In order to quantitatively compare the theoretical result in Figure 6a with th lation result in Figure 6b, the difference between the theoretical and simulation r shown in Table 2. From Table 2, it can be seen that the simulation results are slightly larger th theoretical results; this is because the high harmonics are neglected in the theoret culation, so its harmonic content in the theoretical result is smaller. Moreover, onl points are used in the simulation, which leads to the discontinuity of the duty cycle it is different from the ideal mathematical model obtained by theoretical calcula the graphic of the simulation results is not smooth.
From Figure 6, it can be seen that the harmonic content is smaller when x is 0. is 0 .8 (or x is 0.8 and y is 0.2). Therefore, the value of x is set to 0.2 and y is set to 0 is set to 0.8 and y is set to 0.2). Substituting the values of x and y into (17) and ( From Figure 6a, it can be seen that the harmonic content changes as the law of paraboloid. Based on this law, the distribution of the PWM pulse for the double inverters can be optimized, and the aim of decreasing the grid-side current harmonic content can be achieved.
In order to verify the results of the theoretical calculation, the values of x and y are selected and the a-phase grid-side current is analyzed by use of fast Fourier analysis. The harmonic content is shown in Figure 6b, and it can be seen that the change trends in Figure 6b are similar to Figure 6a.
In order to quantitatively compare the theoretical result in Figure 6a with the simulation result in Figure 6b, the difference between the theoretical and simulation result is shown in Table 2. From Table 2, it can be seen that the simulation results are slightly larger than the theoretical results; this is because the high harmonics are neglected in the theoretical calculation, so its harmonic content in the theoretical result is smaller. Moreover, only finite points are used in the simulation, which leads to the discontinuity of the duty cycles. Thus, it is different from the ideal mathematical model obtained by theoretical calculation, so the graphic of the simulation results is not smooth.
From Figure 6, it can be seen that the harmonic content is smaller when x is 0.2 and y is 0.8 (or x is 0.8 and y is 0.2). Therefore, the value of x is set to 0.2 and y is set to 0.8 (or x is set to 0.8 and y is set to 0.2). Substituting the values of x and y into (17) and (18), the optimized duty cycles of zero vectors xd 0Inv.1 and yd 0Inv.2 are obtained. The corrected distributions of PWM duty ratios of double inverters are shown in Figure 7. From Figure 7, it can be seen that the width of the pulse is not changed, even distribution of the active vector pulse changes. In other words, both the duty-cycles active vector and the amplitude of load-side voltage do not change.
The block diagram of the duty-cycle correction strategy is shown in Figure (7). After obtaining the duty cycles, the current harmonic co under different x and y are compared by use of (27) and the optimal pulse is ob According to (17) and (18) and the pulse distribution coefficient, the corrected duty of zero vectors can be obtained. Finally, the turn-on times of IGBTs in different sect calculated according to the corrected duty cycles and (8) and (9), and the application TAp1~TCn2 are obtained. Comparing the effect time of each IGBT with the triangular c the PWM pulses in Inverter 1 and Inverter 2 can be obtained by logic judgment.

Simulation Results
In order to verify the theoretical results, MATLAB/Simulink was used to es the simulation model of the dual-system TSMC. The parameters of the system are in Table 1.
Employing the modulation strategy without zero vector for the rectifier sta From Figure 7, it can be seen that the width of the pulse is not changed, even if the distribution of the active vector pulse changes. In other words, both the duty-cycles of the active vector and the amplitude of load-side voltage do not change.
The block diagram of the duty-cycle correction strategy is shown in Figure 8. First, the reference voltage of the inverter is obtained, and the angles θ 1 * and θ 2 * of the reference voltage are obtained by PLL. The duty cycles for the inverter d αInv.1 , d βInv.1 , d αInv.2 , and d βInv.2, and the duty cycles of zero vectors d 0Inv.1 and d 0Inv.2 can be calculated according to the angle of reference voltage and (7). After obtaining the duty cycles, the current harmonic contents under different x and y are compared by use of (27) and the optimal pulse is obtained. According to (17) and (18) and the pulse distribution coefficient, the corrected duty cycles of zero vectors can be obtained. Finally, the turn-on times of IGBTs in different sectors are calculated according to the corrected duty cycles and (8) and (9), and the application times T Ap1~TCn2 are obtained. Comparing the effect time of each IGBT with the triangular carrier, the PWM pulses in Inverter 1 and Inverter 2 can be obtained by logic judgment. From Figure 7, it can be seen that the width of the pulse is not changed, even if the distribution of the active vector pulse changes. In other words, both the duty-cycles of the active vector and the amplitude of load-side voltage do not change.
The block diagram of the duty-cycle correction strategy is shown in Figure 8. First, the reference voltage of the inverter is obtained, and the angles θ1* and θ2* of the reference voltage are obtained by PLL. The duty cycles for the inverter dαInv.1, dβInv.1, dαInv.2, and dβInv.2, and the duty cycles of zero vectors d0Inv.1 and d0Inv.2 can be calculated according to the angle of reference voltage and (7). After obtaining the duty cycles, the current harmonic contents under different x and y are compared by use of (27) and the optimal pulse is obtained. According to (17) and (18) and the pulse distribution coefficient, the corrected duty cycles of zero vectors can be obtained. Finally, the turn-on times of IGBTs in different sectors are calculated according to the corrected duty cycles and (8) and (9), and the application times TAp1~TCn2 are obtained. Comparing the effect time of each IGBT with the triangular carrier, the PWM pulses in Inverter 1 and Inverter 2 can be obtained by logic judgment.

Simulation Results
In order to verify the theoretical results, MATLAB/Simulink was used to establish the simulation model of the dual-system TSMC. The parameters of the system are shown in Table 1.
Employing the modulation strategy without zero vector for the rectifier stage, the waveforms of a-phase voltage usa and current isa obtained in simulation are shown in Figure 9, respectively. The output frequency is set to 40 Hz. The load-side line voltage uAB and phase current iA are shown in Figure 10. From these figures, it can be seen that the load-side line voltage is an ideal PWM waveform and output line current is near sinusoidal.

Simulation Results
In order to verify the theoretical results, MATLAB/Simulink was used to establish the simulation model of the dual-system TSMC. The parameters of the system are shown in Table 1.
Employing the modulation strategy without zero vector for the rectifier stage, the waveforms of a-phase voltage u sa and current i sa obtained in simulation are shown in Figure 9, respectively. The output frequency is set to 40 Hz. The load-side line voltage u AB   Figures 11 and 12 show the simulation waveforms and Fourier analy grid-side current when the distribution coefficients x = 0.5, y = 0.5 and respectively. In Figure 11, isa1 is the fundamental waveform of the grid-si the n-th harmonic of the grid-side current. From Figure 11, it can be se side current harmonic content decreases with the variation in the values   Figures 11 and 12 show the simulation waveforms and Fourier analy grid-side current when the distribution coefficients x = 0.5, y = 0.5 and respectively. In Figure 11, isa1 is the fundamental waveform of the grid-si the n-th harmonic of the grid-side current. From Figure 11, it can be see side current harmonic content decreases with the variation in the values  Figures 11 and 12 show the simulation waveforms and Fourier analysis results of the grid-side current when the distribution coefficients x = 0.5, y = 0.5 and x = 0.2, y = 0.8, respectively. In Figure 11, i sa1 is the fundamental waveform of the grid-side current; i san is the n-th harmonic of the grid-side current. From Figure 11, it can be seen that the grid-side current harmonic content decreases with the variation in the values of x and y. Figures 11 and 12 show the simulation waveforms and Fourier analysis grid-side current when the distribution coefficients x = 0.5, y = 0.5 and x respectively. In Figure 11, isa1 is the fundamental waveform of the grid-side the n-th harmonic of the grid-side current. From Figure 11, it can be seen side current harmonic content decreases with the variation in the values of

Analysis of Experimental Results
The experiments were carried out on the platform (shown in Figure  parameters shown in Table 1. The sampling frequency is set to 10 kHz; the s ysis of the grid-side current is executed to verify the distribution law of the

Analysis of Experimental Results
The experiments were carried out on the platform (shown in Figure 13), with the parameters shown in Table 1. The sampling frequency is set to 10 kHz; the spectrum analysis of the grid-side current is executed to verify the distribution law of the harmonics.

Analysis of Experimental Results
The experiments were carried out on the platform (shown in Figure 13 parameters shown in Table 1. The sampling frequency is set to 10 kHz; the spe ysis of the grid-side current is executed to verify the distribution law of the ha  Figure 14 shows the experimental waveforms of th line voltage uAB and current iA of double inverters when the input frequenc From Figure 14, it can be seen that the load-side line voltages of double inve ideal PWM waveform and the output line current is near sinusoidal.  Figures 14 and 15 show the experimental waveforms of the grid-side phase voltage u sa and current i sa and load-side line voltage u AB and current i A , when the input frequency is 50 Hz and x = 0.5, y = 0.5. Figure 14 shows the experimental waveforms of the load-side line voltage u AB and current i A of double inverters when the input frequency is 40 Hz. From Figure 14, it can be seen that the load-side line voltages of double inverters are an ideal PWM waveform and the output line current is near sinusoidal.
The experiments were carried out on the platform (shown in Figure  parameters shown in Table 1. The sampling frequency is set to 10 kHz; the s ysis of the grid-side current is executed to verify the distribution law of the    Figures 16 and 17 show the experimental waveforms of the grid-side phase voltage u sa and current i sa and load-side line voltage u AB and current i A when the input frequency is 50 Hz and x = 0.2, y = 0.8. Figure 16 shows the experimental waveforms of the load-side line voltage u AB and current i A of double inverters when the input frequency is 40 Hz.  Figure 16 shows the experimental waveforms of the load-sid line voltage uAB and current iA of double inverters when the input frequency is 40 Hz. Figure 18 shows the experimental waveforms and the Fourier analysis result of th grid-side phase current isa when x = 0.5, y = 0.5 and x = 0.2, y = 0.8. It can be seen that th harmonics of the grid-side current are suppressed by use of the duty-cycle correction stra egy, and it is consistent with the theoretical analysis.    Figure 16 shows the experimental waveforms of the load line voltage uAB and current iA of double inverters when the input frequency is 40 Hz Figure 18 shows the experimental waveforms and the Fourier analysis result o grid-side phase current isa when x = 0.5, y = 0.5 and x = 0.2, y = 0.8. It can be seen tha harmonics of the grid-side current are suppressed by use of the duty-cycle correction s egy, and it is consistent with the theoretical analysis.    Figure 18 shows the experimental waveforms and the Fourier analysis result of the grid-side phase current i sa when x = 0.5, y = 0.5 and x = 0.2, y = 0.8. It can be seen that the harmonics of the grid-side current are suppressed by use of the duty-cycle correction strategy, and it is consistent with the theoretical analysis.

Conclusions
This paper analyzed the harmonic characteristics of the grid-side current of a dualinverter TSMC system according to PBM, and analyzed the law of the effect from the pulse distribution on the harmonic content of the grid-side current. Further, a modulation strategy to reduce the current harmonics was proposed. The following results were proved by the simulation and experiments.

Conclusions
This paper analyzed the harmonic characteristics of the grid-side current of a dualinverter TSMC system according to PBM, and analyzed the law of the effect from the pulse distribution on the harmonic content of the grid-side current. Further, a modulation strategy to reduce the current harmonics was proposed. The following results were proved by the simulation and experiments.
(1) The grid-side current harmonics in the dual-inverter TSMC are affected by the pulse width and distribution of inverters, and the harmonic content changes with the variation in the distribution coefficients of the pulses in the form of paraboloid.
(2) By optimizing the pulse width and distribution of inverters, the harmonic content of the grid-side currents can be effectively reduced. The grid-side current harmonic content can be effectively reduced when the distribution coefficient x is set to a value near 0.2 (0.8) and y is set to a value near 0.8 (0.2).