同步式SPWM两电平正弦脉宽调制逆变器（全波三角波）——正弦波为调制波（双重傅里叶分析）

SPWM正弦脉宽调制介绍

SPWM是调制波为正弦波、载波为三角波或锯齿波的一种脉宽调制法。

特点：
　原理简单，通用性强，控制和调节性能好，具有先出谐波、调节和稳定输出电压的多种作用，是一种比较好的波形改善法。
分类：
　分为两阶式和三阶式两种。（“阶”指的是PWM式逆变器输出电压在一个周期内的电压电平数）
　当阶数为3时，把输出电压基波半个周期内的脉冲数称为PWM逆变电路的脉冲数，在三相桥式逆变电路中，此脉冲数就是在基波一个周期内，同一个逆变开关的开通次数，他等于载波频率与基波频率（调制波频率）之比。

载波为全波三角波的单相二阶SPWM逆变器

工作原理：

电路与波形图：

谐波分析

从上述的波形图中，可以看出SPWM调制在调制波的各周期内所包含的脉冲模式是没有重复性的，因此，我们无法直接以调制波角频率 w s w_s ws​为基准，用傅里叶级数把他分解为调制波角频率倍数的谐波，而采用以载波的角频率 w c w_c wc​为基准，考察边频带谐波分布情况的方法是比较合适的，这就需要采用双重傅里叶级数的分析法来进行谐波分析，并且这种分析方法对于同步式也是正确的。
　双重傅里叶谐波分析表示：
　在载波 w c w_c wc​的一个周期内，将三角波（负峰值为起点）可以表示，他们的斜率为 + 2 U c π +\frac{2U_c}{\pi} +π2Uc​​和 − 2 U c π -\frac{2U_c}{\pi} −π2Uc​​，初始值为 + U c +U_c +Uc​和 − U c -U_c −Uc​，这样三角波的数学方程式就可以写成这样的形式了：
u c = { − ( w c t − 2 π k − π ) 2 U c π + U c 2 π k + π ≤ w c t ≤ 2 π k + 2 π ( w c t − 2 π k ) 2 U c π + U c 2 π k ≤ w c t ≤ 2 π k + π u_c=\begin{cases} -(w_ct-2\pi k-\pi)\frac{2U_c}{\pi}+U_c & 2\pi k + \pi \le w_ct \le 2\pi k +2\pi \\ (w_ct-2\pi k)\frac{2U_c}{\pi}+U_c & 2\pi k \le w_ct \le 2\pi k +\pi \end{cases} uc​={−(wc​t−2πk−π)π2Uc​​+Uc​(wc​t−2πk)π2Uc​​+Uc​​2πk+π≤wc​t≤2πk+2π2πk≤wc​t≤2πk+π​
而此时正弦调制波函数为：
u s = s i n ( w s t − φ ) u_s=sin(w_st-\varphi) us​=sin(ws​t−φ)
取调制度 U s U c = M ≤ 1 \frac{U_s}{U_c}= M \le 1 Uc​Us​​=M≤1,载波比 w c w s = N ≫ 1 \frac{w_c}{w_s} = N \gg 1 ws​wc​​=N≫1，N等于任意正整数。
　二阶SPWM波的采样点是正弦波与三角波的交点，在交点处存在数学关系 u s = u c u_s = u_c us​=uc​。
　在采样点a： U s s i n ( w s − φ ) = − ( w c t − 2 π k − π ) 2 U c π + U c U_s sin(w_s-\varphi ) = -(w_ct-2\pi k-\pi)\frac{2U_c}{\pi}+U_c Us​sin(ws​−φ)=−(wc​t−2πk−π)π2Uc​​+Uc​.
　取 X = w c t , Y = w s t − φ X = w_ct,Y = w_st-\varphi X=wc​t,Y=ws​t−φ则公式可以写为
　 X = 2 π ( k + 1 ) − π 2 ( 1 + M s i n Y ) X=2\pi (k+1)-\frac{\pi}{2}(1+MsinY) X=2π(k+1)−2π​(1+MsinY)
　在采样点b： X = 2 π k + π 2 ( 1 + M s i n Y ) X =2\pi k + \frac{\pi}{2}(1+MsinY) X=2πk+2π​(1+MsinY)
　同时，按照调制的方式，在 X = w c t X = w_ct X=wc​t的 2 π ( k − 1 2 ) 2\pi (k-\frac{1}{2}) 2π(k−21​)到 2 π ( k + 1 2 ) 2\pi (k+\frac{1}{2}) 2π(k+21​)的区间内，此时调制输出信号的表达式为： u L = { 　 E 2 π ( k + 1 ) − π 2 ( 1 + M s i n Y ) ≤ X ≤ 2 π k + π 2 ( 1 + M s i n Y ) 　 − E 2 π k + π 2 ( 1 + M s i n Y ) ≤ X ≤ 2 π ( k + 1 ) − π 2 ( 1 + M s i n Y ) 　 　 　 （ １ ） u_L = \begin{cases} 　E & 2\pi (k+1)-\frac{\pi}{2}(1+MsinY) \le X \le 2\pi k +\frac{\pi}{2}(1+MsinY)\\ 　-E & 2\pi k +\frac{\pi}{2}(1+MsinY) \le X \le 2\pi (k+1)-\frac{\pi}{2}(1+MsinY) 　\end{cases}　　（１） uL​={　E　−E​2π(k+1)−2π​(1+MsinY)≤X≤2πk+2π​(1+MsinY)2πk+2π​(1+MsinY)≤X≤2π(k+1)−2π​(1+MsinY)　​　　（１）
　其中 Y = w s − φ , X = w c t , w c w s = N , U s U c = M ≤ 1 , k = 0 , 1 , 2... Y = w_s - \varphi , X=w_ct , \frac{w_c}{w_s} = N, \frac {U_s}{U_c} = M \le 1,k=0,1,2... Y=ws​−φ,X=wc​t,ws​wc​​=N,Uc​Us​​=M≤1,k=0,1,2...
　假设m为相对于载波 w c w_c wc​的谐波次数，n为相对于调制波 w s w_s ws​的谐波次数，怎 u L u_L uL​的双重傅里叶级数表达式为
　 u L = 1 2 A 00 （ 直 流 分 量 ） + ∑ n = 1 ∝ ( A 0 n c o s n Y + B 0 n s i n n Y ) （ 调 制 波 谐 波 分 量 ） 　 + ∑ m = 1 ∝ ( A m 0 c o s m X + B m 0 s i n m X ) （ 载 波 谐 波 分 量 ） 　 + ∑ m = 1 ∝ ∑ − ∝ , n ≠ 0 + ∝ ( A m n c o s ( m X + n Y ) + B m n s i n ( m X + n Y ) ) （ 边 带 谐 波 分 量 ） 　 　 　 ( ２ ) u_L = \frac{1}{2} A_{00}（直流分量）+ \sum_{n=1}^{\propto}(A_{0n}cosnY+B_{0n}sinnY) （调制波谐波分量） \\ 　+ \sum_{m=1}^{\propto}(A_{m0}cosmX+B_{m0}sinmX) （载波谐波分量）\\ 　+\sum_{m=1}^{\propto} \sum_{-\propto ,n\ne 0}^{+\propto}(A_{mn}cos(mX+nY)+B_{mn}sin(mX+nY))（边带谐波分量）　　　(２) uL​=21​A00​（直流分量）+n=1∑∝​(A0n​cosnY+B0n​sinnY)（调制波谐波分量）　+m=1∑∝​(Am0​cosmX+Bm0​sinmX)（载波谐波分量）　+m=1∑∝​−∝,n​=0∑+∝​(Amn​cos(mX+nY)+Bmn​sin(mX+nY))（边带谐波分量）　　　(２)
　式中关系为：
　 A m n + j B m n = 2 ( 2 π ) 2 ∫ − π + π ∫ − π + π U L ⋅ e j ( m X + n Y ) d X d Y , m = 1 , 2 , 3... 　 　 （ ３ ） A_{mn}+jB_{mn} = \frac{2}{(2\pi)^{2}}\int_{-\pi}^{+\pi}\int_{-\pi}^{+\pi}U_L\cdot e^{j(mX+nY)}dXdY, m=1,2,3...　　（３） Amn​+jBmn​=(2π)22​∫−π+π​∫−π+π​UL​⋅ej(mX+nY)dXdY,m=1,2,3...　　（３）
　将在一个周期内的方程（1）代入（３）进行傅里叶积分的求解为：
　 A m n + j B m n = − 2 E ( 2 π ) 2 ∫ − π + π ∫ 2 π ( k − 1 2 ) 2 π ( k + 1 ) + π 2 ( 1 + M s i n Y ) e j ( m X + n Y ) d X d Y 　 + 2 E ( 2 π ) 2 ∫ − π + π ∫ 2 π ( k + 1 ) − π 2 ( 1 + M s i n Y ) 2 π k + π 2 ( 1 + M s i n Y ) e j ( m X + n Y ) d X d Y 　 − 2 E ( 2 π ) 2 ∫ − π + π ∫ 2 π k + π 2 ( 1 + M s i n Y ) 2 π ( k + 1 2 ) e i ( m X + n Y ) d X d Y 　 （ ４ ） A_{mn}+jB_{mn} = -\frac{2E}{(2\pi)^2}\int_{-\pi}^{+\pi}\int_{2\pi(k-\frac{1}{2})}^{2\pi(k+1)+\frac{\pi}{2}(1+MsinY)}e^{j(mX+nY)}dXdY \\ 　 +\frac{2E}{(2\pi)^2}\int_{-\pi}^{+\pi}\int_{2\pi(k+1)-\frac{\pi}{2}(1+MsinY)}^{2\pi k+\frac{\pi}{2}(1+MsinY)}e^{j(mX+nY)}dXdY \\ 　 -\frac{2E}{(2\pi)^2}\int_{-\pi}^{+\pi}\int_{2\pi k+\frac{\pi}{2}(1+MsinY)}^{2\pi(k+\frac{1}{2})}e^{i(mX+nY)}dXdY　（４） Amn​+jBmn​=−(2π)22E​∫−π+π​∫2π(k−21​)2π(k+1)+2π​(1+MsinY)​ej(mX+nY)dXdY　+(2π)22E​∫−π+π​∫2π(k+1)−2π​(1+MsinY)2πk+2π​(1+MsinY)​ej(mX+nY)dXdY　−(2π)22E​∫−π+π​∫2πk+2π​(1+MsinY)2π(k+21​)​ei(mX+nY)dXdY　（４）
　 因为 e j m ( 2 π ( k − 1 2 ) ) = e j m ( 2 π ( k − 1 2 ) ) = ( − 1 ) m e^{jm(2\pi(k-\frac{1}{2}))} = e^{jm(2\pi(k-\frac{1}{2}))}=(-1)^m ejm(2π(k−21​))=ejm(2π(k−21​))=(−1)m故公式（4）可以化简为 ( 4 ) = − E j m 2 ( π ) 2 ∫ − π + π [ 2 e j m ( 2 π k + π 2 ( 1 + M s i n Y ) − 2 e j m ( 2 π ( k + 1 ) − π 2 ( 1 + M s i n Y ) ⋅ e j n Y d X d Y ] (4)=-\frac{E}{jm2(\pi)^2}\int_{-\pi}^{+\pi}[2e^{jm(2\pi k+\frac{\pi}{2}(1+MsinY)}-2e^{jm(2\pi(k+1)-\frac{\pi}{2}(1+MsinY)}\cdot e^{jnY}dXdY] (4)=−jm2(π)2E​∫−π+π​[2ejm(2πk+2π​(1+MsinY)−2ejm(2π(k+1)−2π​(1+MsinY)⋅ejnYdXdY]
　 （计算易错点：进行比较大的计算的时候，一定要先搞清楚符号是前提，再就是公式套用的要准确，才能计算准确）
　 且存在关系有 e j m 2 π k = e i m 2 π ( k + 1 ) = 1 e^{jm2\pi k}=e^{im2\pi(k+1)}=1 ejm2πk=eim2π(k+1)=1则将此关系代入（4）式中，可以得到此时的化简关系为：
　 A m n + j B m n = − j E m ( π ) 2 ∫ − π + π [ e j ( m π 2 + m M π 2 s i n Y ) − e − j ( m π 2 + m M π 2 s i n Y ) ] e j n Y d X d Y A_{mn}+jB_{mn}=-\frac{jE}{m(\pi)^2}\int_{-\pi}^{+\pi}[e^{j(\frac{m\pi}{2}+\frac{mM\pi}{2}sinY)}-e^{-j(\frac{m\pi}{2}+\frac{mM\pi}{2}sinY)}]e^{jnY}dXdY Amn​+jBmn​=−m(π)2jE​∫−π+π​[ej(2mπ​+2mMπ​sinY)−e−j(2mπ​+2mMπ​sinY)]ejnYdXdY
　 由bessel理论以及 e j n π = ( − 1 ) n e^{jn\pi}=(-1)^n ejnπ=(−1)n,得：
　
代入化简可以得到目标的表达式化简为：

（在化简过程中需要注意到关系 s i n n π = 2 s i n n π 2 ⋅ cos ⁡ n π 2 = 0 sinn\pi=2sin\frac{n\pi}{2}\cdot \cos\frac{n\pi}{2}=0 sinnπ=2sin2nπ​⋅cos2nπ​=0的化简关系的使用 且记作这个公式为（5））

分类讨论谐波分析:

1、从结果的目标式可以知道，当m=0,n=时， s i n m + n 2 π = 0 sin\frac{m+n}{2}\pi=0 sin2m+n​π=0,因此恒定分量 A 00 = 0 A_{00}=0 A00​=0;
2、对于基波和基波的谐波分析，此时m=0，
此时，相关的关系为：

由傅里叶级数的关系： ∫ − π + π s i n Y c o s n Y ⋅ d Y = 0 \int_{-\pi}^{+\pi}sinYcosnY\cdot dY=0 ∫−π+π​sinYcosnY⋅dY=0,因此存在关系：
A 0 n + j B 0 n = j M E π ∫ − π π s i n Y ⋅ s i n n Y ⋅ d Y A_{0n}+jB_{0n}=\frac{jME}{\pi}\int_{-\pi}^{\pi}sinY\cdot sinnY\cdot dY A0n​+jB0n​=πjME​∫−ππ​sinY⋅sinnY⋅dY,由傅里叶级数关系，可以知道这个表达式只在基波，n=1的时候才存在数值即 B 01 = M E , A 01 = 0 B_{01}=ME,A_{01}=0 B01​=ME,A01​=0，即对于基波的谐波是不存在的。
3、对于载波及载波的谐波，n=0,由公式（5）得：
A m 0 + j B m 0 = 4 E m π J 0 ( m M π 2 ） s i n m π 2 A_{m0}+jB_{m0}=\frac{4E}{m\pi} J_0(\frac{mM\pi}{2}）sin\frac{m\pi}{2} Am0​+jBm0​=mπ4E​J0​(2mMπ​）sin2mπ​
即：
A m 0 = 4 E m π J 0 ( m M π 2 ) s i n m π 2 , B m 0 = 0 A_{m0}=\frac{4E}{m\pi}J_0(\frac{mM\pi}{2})sin\frac{m\pi}{2},B_{m0}=0 Am0​=mπ4E​J0​(2mMπ​)sin2mπ​,Bm0​=0他的表达式中只存在实数项（即输出对应的余弦项），且这个项的输出表达式可以看出，当m为偶数的时候 A m 0 = 0 A_{m0}=0 Am0​=0,无输出，而当m为奇数的时候才有对应的谐波输出，即输出的载波的谐波只有奇数次谐波，而不存在偶次谐波。
4、载波及载波m次谐波的上下变频分析：
此时的表达式为：
A m n = 4 E m π J n ( m M π 2 ) s i n ( m + n 2 π ) ⋅ c o s n π 2 A_{mn}=\frac{4E}{m\pi}J_n(\frac{mM\pi}{2})sin(\frac{m+n}{2}\pi)\cdot cos\frac{n\pi}{2} Amn​=mπ4E​Jn​(2mMπ​)sin(2m+n​π)⋅cos2nπ​
B m n = 4 E m π J n ( m M π 2 ) s i n ( m + n 2 π ) ⋅ s i n n π 2 B_{mn}=\frac{4E}{m\pi}J_n(\frac{mM\pi}{2})sin(\frac{m+n}{2}\pi)\cdot sin\frac{n\pi}{2} Bmn​=mπ4E​Jn​(2mMπ​)sin(2m+n​π)⋅sin2nπ​
对输出谐波表达式分析为：
当m为奇数，n为偶数 ( s i n n π 2 = 0 ) (sin\frac{n\pi}{2}=0) (sin2nπ​=0)此时的m+n为奇数时，输出只存在 A m n A_{mn} Amn​项；
当m为奇数，n为奇数 ( c o s n π 2 = 0 ) (cos\frac{n\pi}{2}=0) (cos2nπ​=0)此时的m+n为偶数 ( s i n ( m + n 2 π ) = 0 ) (sin(\frac{m+n}{2}\pi)=0) (sin(2m+n​π)=0).没有对应的谐波输出；
当m为偶数，n为偶数，此时的m+n为偶数 ( s i n ( m + n 2 π ) = 0 ) (sin(\frac{m+n}{2}\pi)=0) (sin(2m+n​π)=0).没有对应的谐波输出；
当m为偶数，n为奇数 ( c o s n π 2 = 0 ) (cos\frac{n\pi}{2}=0) (cos2nπ​=0)，此时的m+n为奇数时，输出只存在 B m n B_{mn} Bmn​项；
由于 w c w s = N \frac{w_c}{w_s}=N ws​wc​​=N是任意正整数，故将上面得到的系数结果代入方程中，便可以得到此时 u L u_L uL​的双重傅里叶级数表达式为：
u L = M E s i n ( w s t − φ ) + 4 E π ∑ m = 1 , 3 , 5 , 7... ∝ J 0 ( m M π 2 ) m s i n m π 2 c o s ( m N w s t ) + 4 E π ∑ m = 1 , 3 , 5 , 7... ∝ ∑ n = ± 2 , ± 4 , ± 6... ± ∝ J n ( m M π 2 ) m s i n ( m + n 2 π ) c o s n π 2 c o s [ ( m N + n ) w s t − n φ ] + 4 E π ∑ m = 2 , 4 , 6 , 8... ∝ ∑ n = ± 1 , ± 3 , ± 5... ± ∝ J n ( m M π 2 ) m s i n ( m + n 2 π ) s i n n π 2 s i n [ ( m N + n ) w s t − n φ ] 　 （ ６ ） u_L=MEsin(w_st-\varphi)+\frac{4E}{\pi}\sum_{m=1,3,5,7...}^{\propto}\frac{J_0(\frac{mM\pi}{2})}{m}sin\frac{m\pi}{2}cos(mNw_st)\\ +\frac{4E}{\pi}\sum_{m=1,3,5,7...}^{\propto}\sum_{n=\pm 2,\pm 4,\pm 6...}^{\pm \propto}\frac{J_n(\frac{mM\pi}{2})}{m}sin(\frac{m+n}{2}\pi)cos\frac{n\pi}{2}cos[(mN+n)w_st-n\varphi]\\ +\frac{4E}{\pi}\sum_{m=2,4,6,8...}^{\propto}\sum_{n=\pm 1,\pm 3,\pm 5...}^{\pm \propto}\frac{J_n(\frac{mM\pi}{2})}{m}sin(\frac{m+n}{2}\pi)sin\frac{n\pi}{2}sin[(mN+n)w_st-n\varphi]　（６） uL​=MEsin(ws​t−φ)+π4E​m=1,3,5,7...∑∝​mJ0​(2mMπ​)​sin2mπ​cos(mNws​t)+π4E​m=1,3,5,7...∑∝​n=±2,±4,±6...∑±∝​mJn​(2mMπ​)​sin(2m+n​π)cos2nπ​cos[(mN+n)ws​t−nφ]+π4E​m=2,4,6,8...∑∝​n=±1,±3,±5...∑±∝​mJn​(2mMπ​)​sin(2m+n​π)sin2nπ​sin[(mN+n)ws​t−nφ]　（６）
对公式进行化简，由三角函数的诱导公式为：

所以此时的 u L u_L uL​可以写成：
u L = M E s i n ( w s t − φ ) + 4 E π ∑ m = 1 , 3 , 5 , 7... ∝ J 0 ( m M π 2 ) m s i n m π 2 c o s ( m N w s t ) + 4 E π ∑ m = 1 , 3 , 5 , 7... ∝ ∑ n = ± 2 , ± 4 , ± 6... ± ∝ J n ( m M π 2 ) m s i n m π 2 c o s [ ( m N + n ) w s t − n φ ] + 4 E π ∑ m = 2 , 4 , 6 , 8... ∝ ∑ n = ± 1 , ± 3 , ± 5... ± ∝ J n ( m M π 2 ) m c o s m π 2 s i n [ ( m N + n ) w s t − n φ ] 　 （ 7 ） u_L=MEsin(w_st-\varphi)+\frac{4E}{\pi}\sum_{m=1,3,5,7...}^{\propto}\frac{J_0(\frac{mM\pi}{2})}{m}sin\frac{m\pi}{2}cos(mNw_st)\\ +\frac{4E}{\pi}\sum_{m=1,3,5,7...}^{\propto}\sum_{n=\pm 2,\pm 4,\pm 6...}^{\pm \propto}\frac{J_n(\frac{mM\pi}{2})}{m}sin\frac{m\pi}{2}cos[(mN+n)w_st-n\varphi]\\ +\frac{4E}{\pi}\sum_{m=2,4,6,8...}^{\propto}\sum_{n=\pm 1,\pm 3,\pm 5...}^{\pm \propto}\frac{J_n(\frac{mM\pi}{2})}{m}cos\frac{m\pi}{2}sin[(mN+n)w_st-n\varphi]　（7） uL​=MEsin(ws​t−φ)+π4E​m=1,3,5,7...∑∝​mJ0​(2mMπ​)​sin2mπ​cos(mNws​t)+π4E​m=1,3,5,7...∑∝​n=±2,±4,±6...∑±∝​mJn​(2mMπ​)​sin2mπ​cos[(mN+n)ws​t−nφ]+π4E​m=2,4,6,8...∑∝​n=±1,±3,±5...∑±∝​mJn​(2mMπ​)​cos2mπ​sin[(mN+n)ws​t−nφ]　（7）
或：
u L = M E s i n ( w s t − φ ) + 4 E π ∑ m = 1 , 3 , 5 , 7... ∝ J 0 ( m M π 2 ) m s i n m π 2 c o s ( m N w s t ) + 4 E π ∑ m = 1 , 2 , 3 , 4... ∝ ∑ n = ± 1 , ± 2 , ± 3... ± ∝ J n ( m M π 2 ) m s i n ( m + n 2 π ) ⋅ { c o s n π 2 c o s [ ( m N + n ) w s t − n φ ] + s i n n π 2 s i n [ ( m N + n ) w s t − n φ ] } 　 （ 8 ） u_L=MEsin(w_st-\varphi)+\frac{4E}{\pi}\sum_{m=1,3,5,7...}^{\propto}\frac{J_0(\frac{mM\pi}{2})}{m}sin\frac{m\pi}{2}cos(mNw_st)\\ +\frac{4E}{\pi}\sum_{m=1,2,3,4...}^{\propto}\sum_{n=\pm 1,\pm 2,\pm 3...}^{\pm \propto}\frac{J_n(\frac{mM\pi}{2})}{m}sin(\frac{m+n}{2}\pi)\cdot \\ \{cos\frac{n\pi}{2}cos[(mN+n)w_st-n\varphi]+sin\frac{n\pi}{2}sin[(mN+n)w_st-n\varphi]\}　（8） uL​=MEsin(ws​t−φ)+π4E​m=1,3,5,7...∑∝​mJ0​(2mMπ​)​sin2mπ​cos(mNws​t)+π4E​m=1,2,3,4...∑∝​n=±1,±2,±3...∑±∝​mJn​(2mMπ​)​sin(2m+n​π)⋅{cos2nπ​cos[(mN+n)ws​t−nφ]+sin2nπ​sin[(mN+n)ws​t−nφ]}　（8）
最终化简为：
u L = M E s i n ( w s t − φ ) + 4 E π ∑ m = 1 , 3 , 5 , 7... ∝ J 0 ( m M π 2 ) m s i n m π 2 c o s ( m N w s t ) + 4 E π ∑ m = 1 , 2 , 3 , 4... ∝ ∑ n = ± 1 , ± 2 , ± 3... ± ∝ J n ( m M π 2 ) m s i n ( m + n 2 π ) { c o s [ ( m N + n ) w s t − n φ − n π 2 ] } 　 （ 9 ） u_L=MEsin(w_st-\varphi)+\frac{4E}{\pi}\sum_{m=1,3,5,7...}^{\propto}\frac{J_0(\frac{mM\pi}{2})}{m}sin\frac{m\pi}{2}cos(mNw_st)\\ +\frac{4E}{\pi}\sum_{m=1,2,3,4...}^{\propto}\sum_{n=\pm 1,\pm 2,\pm 3...}^{\pm \propto}\frac{J_n(\frac{mM\pi}{2})}{m}sin(\frac{m+n}{2}\pi)\{cos[(mN+n)w_st-n\varphi -\frac{n\pi}{2}]\}　（9） uL​=MEsin(ws​t−φ)+π4E​m=1,3,5,7...∑∝​mJ0​(2mMπ​)​sin2mπ​cos(mNws​t)+π4E​m=1,2,3,4...∑∝​n=±1,±2,±3...∑±∝​mJn​(2mMπ​)​sin(2m+n​π){cos[(mN+n)ws​t−nφ−2nπ​]}　（9）

结论:当载波为全波三角波的单相二阶SPWM逆变器输出电压的谐波包含下面的成分：基波、载波、载波的m次谐波，载波及载波m次谐波的上下边频谐波，其中基波幅值与调制度成正比，故通过调制调制波（正弦波）幅值的大小就可以调节输出电压。当m为偶数时，载波的m次谐波也不存在；当m+n为偶数时，载波与载波的m次谐波的上下边频谐波也不存在。
　当M=0.2~1时，由方程所计算的频谱分布值如下表。由此表可以画出N等于任意正整数时的频谱分布，如下图所示。由图可以看出，频谱的分布不只与M的大小有关，也与载波比N有关；M可以改变谐波的幅值，而N可以改变谐波的频率。 N越大，谐波频率越高，滤波越容易。当N达到一定值时，只用几个微法的电容就可以达到比较好的滤波效果（可以将设备做的比较小），这就是提高载波频率的主要原因。
　
从中可以看出存在奇数次的m次的载波，而不存在偶次载波的谐波。

讨论各种情况下对于SPWM逆变的谐波分析

1、同步式SPWM在正弦调制波起点位置不同时，二阶SPWM波的双重傅里叶级数表达式的不同：

对于同步式SPWM，相对于载波三角波来说，正弦调制波起始点的位置共有四种，即起始点在三角波的正、负峰值处，上升沿与下降沿的零点处。
下表是根据不同的起点的载波三角波的相应波形，以及根据波形解出的二阶SPWM波的双重傅里叶级数表达式。

结论分析：

• 正弦波调制起点不同对方程的基波与谐波的幅值是没有影响的，对方程的谐波次数也没有影响。
• 他只影响到第二项与第三项前的正、负号，以及是正弦还是余弦函数。（起始点在三角波福峰值处，方程的第二项、第三项是正的余弦函数；在三角波的正峰值处，方程的第二项、第三项为负余弦函数；起始点在三角波下降沿零点处，方程的第二项、第三项是正的正弦函数；起始点在三角波的上升沿零点处时，方程的第二项、第三项是负的正弦函数）
• 即调制波起始点在三角波的正、负峰值处时，方程中含有正弦项和余弦项；而调制波起始点在三角波的上升或下降沿零点处时，方程中只含有正弦项，此时 u L ( t ) u_L(t) uL​(t)是奇函数，波形对称于原点。所以，就输出波形而论，正弦调制起始点位置在三角波的上升或下降沿零点处为好，他可以得到对称于原点的奇函数波形，这对于N值较小的同步调制十分有利的。

同步式SPWM N取奇数或偶数时对输出波形的影响

输出谐波方程：
u L = M E s i n ( w s t − φ ) + 4 E π ∑ m = 1 , 3 , 5 , 7... ∝ J 0 ( m M π 2 ) m s i n m π 2 c o s ( m N w s t ) + 4 E π ∑ m = 1 , 3 , 5 , 7... ∝ ∑ n = ± 2 , ± 4 , ± 6... ± ∝ J n ( m M π 2 ) m s i n m π 2 c o s [ ( m N + n ) w s t − n φ ] + 4 E π ∑ m = 2 , 4 , 6 , 8... ∝ ∑ n = ± 1 , ± 3 , ± 5... ± ∝ J n ( m M π 2 ) m c o s m π 2 s i n [ ( m N + n ) w s t − n φ ] u_L=MEsin(w_st-\varphi)+\frac{4E}{\pi}\sum_{m=1,3,5,7...}^{\propto}\frac{J_0(\frac{mM\pi}{2})}{m}sin\frac{m\pi}{2}cos(mNw_st)\\ +\frac{4E}{\pi}\sum_{m=1,3,5,7...}^{\propto}\sum_{n=\pm 2,\pm 4,\pm 6...}^{\pm \propto}\frac{J_n(\frac{mM\pi}{2})}{m}sin\frac{m\pi}{2}cos[(mN+n)w_st-n\varphi]\\ +\frac{4E}{\pi}\sum_{m=2,4,6,8...}^{\propto}\sum_{n=\pm 1,\pm 3,\pm 5...}^{\pm \propto}\frac{J_n(\frac{mM\pi}{2})}{m}cos\frac{m\pi}{2}sin[(mN+n)w_st-n\varphi] uL​=MEsin(ws​t−φ)+π4E​m=1,3,5,7...∑∝​mJ0​(2mMπ​)​sin2mπ​cos(mNws​t)+π4E​m=1,3,5,7...∑∝​n=±2,±4,±6...∑±∝​mJn​(2mMπ​)​sin2mπ​cos[(mN+n)ws​t−nφ]+π4E​m=2,4,6,8...∑∝​n=±1,±3,±5...∑±∝​mJn​(2mMπ​)​cos2mπ​sin[(mN+n)ws​t−nφ]
波形图：

结论：
在四种起始位置中，比较这8种波形可以知道：
当N为偶数时，方程第二项中mN为偶数，第三项中mN+n为偶数，第四项mN+n为奇数，所以方程中的第二、三项为偶次谐波；
当N为奇数时，方程第二项中mN为奇数，第三项mN+n为奇数，第四项中mN+n为奇数，所以方程中的所有项都是奇数次谐波，因此，就输出波形而言，N=奇数为好。
他可以让方程中的谐波全为奇次谐波，所得到的波形既对称于原点，又对称于纵轴（镜对称）所以为奇谐波函数。
当N=奇数时，比较正弦调制波起始点位置在三角波下降沿零点处和上升沿零点处两种波形可知，起始点在上升沿零点处，所得第一个脉冲为负脉冲，而当N值很小（如：N=3）时波形很差，不如起始点在下降沿处的波形好。
总结：在选择调制波起始点位置和选择N为奇数或偶数时，应选取起始点在下降沿零点处，N=奇数为最好