SVPWM仿真详细说明important

附SVPWM的仿真实现1SVPWM的基本原理SPWM常用于变频调速控制系统，经典的SPWM控制主要目的是使变频器的输出电压尽量接近正弦波，并未关注输出的电流波形。而矢量控制的最终目的是得到圆形的旋转磁场，这样就要求变频器输出的电流波形接近正弦波。锁定得到圆形的旋转磁场这一目标，SVPWM控制技术利用逆变器各桥臂开关控制信号的不同组合，使逆变器的输出电压空间矢量的运行轨迹尽可能接近圆形。SVPWM是从电动机的角度出发，着眼于使电机获得幅值恒定的圆形磁场。图1所示为PWM逆变器的拓扑结构以及等效开关模型。ASBSCS4622dU2dU0'135ABC逆变器拓扑结构等效开关模型图1PWM逆变器电路电压源型逆变器常采用180导通型。用ABCSSS、、分别标记三个桥臂的状态，规定当上桥臂器件导通时桥臂状态为1，下桥臂导通时桥臂状态为0，当3个桥臂的功率开关管变化时，就会得到328种开关模式，每种开关模式对应一个电压矢量，矢量的幅值为23dU；有两种开关模式对应的电压矢量幅值为零，称为零矢量。例如：在某一时刻，设V1，V2，V3管处于开通状态，即10abcss,s，设为三相对称负载，各开关管的开通电阻均相等，则逆变器的等效电路为：图210abcss,s时逆变器的等效电路图这样，很容易就能得到该瞬时时刻的相电压：112333ANdBNdCNdvU,vU,vU(1)将其在静止坐标系中表示出来，如图3所示：图310abcss,s电压矢量图其中，U是合成的电压矢量，在两相静止坐标系（,坐标系）下，利用相电压合成电压矢量U的表达式：U2433jjANBNCNk(vveve)(2)其中，k为三相静止坐标系向两相静止坐标系转换的变换系数，变换分为基于等功率的坐标变换和基于等量的坐标变换，这里选择等量的坐标变换，则23k，式(2)即为：U243323jjANBNCN(vveve)(3)将式(1)的具体数值代入上式，则有：U1323jdUe(4)这样就得到了10abcss,s开关状态下的电压矢量，按照同样的方法分析另外7种开关状态，可以分别得到每种开关状态所对应的电压矢量，总结为表1所示。表1逆变器的不同开关状态对应的空间矢量表abcsss相电压矢量表达式矢量标号A相B相C相00000000u00113dU13dU23dU4323jdUe1u01013dU23dU13dU2323jdUe2u01123dU13dU13dU23jdUe3u10023dU13dU13dU023jdUe4u10113dU23dU13dU5323jdUe5u11013dU13dU23dU1323jdUe6u11100007u观察上表可知，三相VSR逆变器在不同的开关组合时的交流侧电压可用一个模为23dU的空间电压矢量在复平面上表示出来，这样就会得到8条空间矢量，如图4所示。图4电压矢量的空间分布与扇区分配显然触发电路每给逆变器发一组触发脉冲，就会在逆变器的交流侧得到一个电压矢量。SVPWM控制的最终目标是获得圆形的旋转电压矢量轨迹，在仅靠这8个电压矢量而不采取任何其它办法的情况下，就只能够得到轨迹为正六边形的旋转电压矢量。这与我们所追求的圆形旋转电压矢量相差甚远，必须引入多个中间矢量以逼近圆形的电压矢量轨迹，可以通过6个非零电压矢量和2个零电压矢量来合成我们所需要的中间矢量。虽然在同一时刻不可能存在两种开关状态，即不可能有两个电压矢量存在，但是若逆变器功率管的开关频率比其输出电压的频率高的多（100倍），每个电压矢量作用的时间极短，则就可以用基本的电压矢量来合成中间电压矢量，以逼近圆形的电压矢量轨迹。2SVPWM仿真模块的搭建上一节介绍了SVPWM控制技术的基本原理，本节的主要内容是介绍如何在Matlab/Simulink具体的实现这种技术。通过本节，要构建出一个可以实现这种SVPWM控制算法的模块，该模块的输入端为控制器发出的控制信号（***,,cbauuu），输出端应为6路触发脉冲。该模块主要包括以下子模块：扇区选择（SectorSelector）子模块；时间计算（TimeCalculating）子模块；时间配合（TimeMatching）子模块；触发脉冲产生（PulsesGenetator）子模块；2.1扇区的选择采用追踪电压型SVPWM控制技术的PWM整流器，其追踪的电压指令就是控制器发出的电压指令u***uju，**u,u分别是两相静止坐标系下,轴分量，它们均是时变的交流量，且相位相差90。***,,cbauuu分别为电压指令u*在三相旋转坐标系下的分量。所谓追踪电压型的SVPWM，就是利用8个基本的电压矢量去追踪给定电压矢量。六个长度不为零的矢量将一个周期分成了6个扇区，为了减少管子的开关次数以及增加系统的稳定性，合成目标矢量采用其所在扇区最近两个基本矢量和两个零矢量共同合成。如图4所示，例如当电压矢量指令u*出现在第Ⅰ扇区时，应当用0u、4u、6u、7u来合成中间电压矢量以追踪电压指令。表2基矢量选择表指令电压所在扇区ⅠⅡⅢⅣⅤⅥ选取的基电压矢量4u、6u0u、7u2u、6u0u、7u2u、3u0u、7u1u、3u0u、7u1u、5u0u、7u4u、5u0u、7u但是，我们还必须知道，以上仅是在已知指令电压矢量所在扇区下所进行的讨论，那么如何确定电压指令矢量所在的扇区？从图4可以看出，*u的正负可以决定矢量u*上半部分的三个扇区或者下半部分的三个扇区，剩下的任务就是判断u*在三个扇区中的哪一个，以区分Ⅰ、Ⅱ、Ⅲ为例，考虑临界情况如下页图5所示：图5临界扇区的判断由图（a）所示：033013302aausin()ucos()u...Nsin()ucos()u...N(1)由图（b）所示：033033302aausin()ucos()u...Nsin()ucos()u...N(2)式中N为扇区，为方便起见，令：1223333refrefarefauuusin()ucos()uusin()ucos()u(3)则可得到第Ⅰ扇区的判别条件为：122000refrefrefuuu(4)同理，其它各个扇区都可以通过这种方法列出判别条件，最后可得到参考电压123refrefrefu,u,u与电压指令u*所在关系如下表所示：表3扇区判断表3refu2refu1refuNM000111011001101010ⅡⅥⅠⅣⅢⅤ123456表中123refrefrefu,u,u大于零时取1，小于零时取0，N为扇区号。该算法可以很容易地判断电压指令u*所在扇区，且算法中部存在除法，因而不会有截断误差。由于在判断扇区过程中要用到电压指令在两相静止坐标系下,轴分量，而给定为三相旋转坐标系下的指令电压矢量，所以控制信号要先经过从三项旋转坐标系到亮相静止坐标系的变换。其变换关系矩阵为：/111222333022ABCC（5）MATLAB仿真模型为：2b1a0.866sqrt(3)/20.816sqrt(2/3)0.816sqrt(2/3)-0.866-sqrt(3)/2-0.5-1/2-0.5-1/23C2B1A图63/2变换模块最终生成的3/2变换模块为：ABCabC3s_2s2图7最终生成的3/2变换模块在MATLAB/Simulink环境下用来实现扇区的模块如下页图8所示。图中，Ref1、Ref2、Ref3是三个选通开关，当中间的输入信号大于零时，输出为1，小于零时输出为0。模块的总输出信号M是按照Ref3、Ref2、Ref1的顺序排列得到的二进制数值，并非实际中的扇区值，但是却与实际扇区间存在一一对应的关系，如表2所示。当然，也可以通过多路选通开关实现M到N的转换，但在实际中没有转换的必要，因为我们最终想得到的只是电压指令所在的空间位置，与各空间位置的编号没有关系。换言之，也可以按照M所在的位置安排扇区的编号，但出于习惯做法，各扇区仍按照图4进行分配。图8扇区选择的Matlab/Simulink实现最终生成的扇区选择子模块如下：图9扇区选择子模块当输入为图10所示的三项正弦信号时，输出波形为扇区序号波形，如图11所示。图10输入控制信号图11山区选择输出信号2.2时间计算在判定了指令电压矢量u*所在的扇区和所需要的基电压矢量后，接着计算两空间矢量的作用时间，仍以图5所示Ⅰ号扇区为例。设在一个开关周期（sT）内，120T,T,T分别为4u、6u和零矢量的作用时间，则由图4.4知：(6)将u***uju代入上式，并结合1223duuU，得：1212221206060*ss*ssTTuuucosTTTuusinTTTTT(7)1201233223ssddsdsTTTuuUUTTuUTTTT(8)这样就得到了电压指令u*在第Ⅰ扇区时，用来追踪电压指令的各基电压矢量作用的时间，同样的方法用于分析在其它扇区时的情况，可得u*在各个扇区的作用时间12T,T如下表所示：表4各个扇区中12T,T对应关系表ⅠⅡⅢⅣⅤⅥ1T-ZZX-X-YY2TXY-YZ-Z-X其中XYZ的值为：333223322sdssddssddTXuUTTYuuUUTTZuuUU(9)图12XYZ计算的Matlab/Simulink实现图1312T,T计算的Matlab/Simulink实现需要指出的是，在计算12T,T时有可能出现12sTTT的情况，因此，还必须进行12T,T的标准化:*1112*2212TTTTTTTTTT(10)即要对上述计算出来的电压矢量的作用时间进行调整，具体方法如式10所示，实现的模型如下：图1412T,T标准化的Matlab/Simulink实现图1512T,T计算子模块时间标准化后输出波形如图16所示。图16T1’波形2.3矢量合成方法研究与时间匹配用基电压矢量合成中间电压矢量追踪指令电压矢量，虽然在功率开关管的开关频率远大于输出电压频率时可近似认为它们同时存在，但是这毕竟是一种近似而实际中又不可能出现的情况，因此，有必要仔细研究基矢量的合成问题。仍以电压指令u*在第Ⅰ扇区时为例来说明常用的矢量合成方法。图14给出了三种常用的矢量合成方法：单三角形法，将零矢量（0u、7u）均匀地分布在指令电压矢量u*的起、终点上，然后依次由4u、6u按三角形方法合成。该方法的特点是：PWM谐波分量主要集中在开关频率swf及2swf上，在频率swf处谐波幅值较大。双三角形法，将零矢量（0u、7u）均匀地分布在指令电压矢量u*的起、终点上，但两空间矢量在中点相交而形成两个三角形，这种方法的开关函数波形对称。PWM谐波分量仍主要分布在开关频率的整数倍附近,谐波幅值比方法a）有所降低。改进的双三角形法，这种方法与b）相似,不同的是在矢量u*的中点处插入了零矢量7u，这样做的好处在于在频率swf处的谐波幅值明显降低。图17三种常用的矢量合成方法比较上述的三种方法，虽然法（c）开关频率较高且算法较复杂，但现代的IPM模块以及TI的DSP芯片完全能够满足要求，为了达到最佳的输出电压波形，本文采用该种方法。下面将详细介绍此法的合成过程。记0127,,,,sTTTTT分别为开关周期、0u、4u、6u、7u作用的时间，为了叙述的方便，引入123456,,,,,tttttt：001012012123470105276127;;;2222;22TTTTTTTTTttttTTTTtTTtTTT图18一个开关周期内基矢量变化图开关状态(abcsss)的变换过程为：000-100-110-111-110-100-000显然，每次变化只有一个功率开关管的状态