已经得到验证
function CMPR=svpwm(ualfa,ubeta)
global ts udc_ref
%在IF语句中出现的变量必须赋初值
persistent T1on T2on T3on t1 t2 sectorNO
if isempty(T1on)|isempty(T2on)|isempty(T3on)|isempty(t1)|isempty(t2)|isempty(sectorNO)
T1on=0
T2on=0
T3on=0
t1=0
t2=0
sectorNO=0
end
%#################################
%sectorNO 3(0 30')->1->5->4->6->2->3
%#################################
if ubeta>0
a=1
else a=0
end
if sqrt(3)*ualfa-ubeta>0
b=1
else b=0
end
if -sqrt(3)*ualfa-ubeta>0
c=1
else c=0
end
sectorNO=a+2*b+4*c
%#################################
%
%#################################
x=sqrt(3)*ubeta*ts/udc_ref
y=(1.5*ualfa+sqrt(3)/2*ubeta)*ts/udc_ref
z=(-1.5*ualfa+sqrt(3)/2*ubeta)*ts/udc_ref
%#################################
%voltage vector on time
%#################################
switch sectorNO
case 3
t1=-z
t2=x
case 1
t1=z
t2=y
case 5
t1=x
t2=-y
case 4
t1=-x
t2=z
case 6
t1=-y
t2=-z
case 2
t1=y
t2=-x
end
%#################################
%pwm timer value
%#################################
%饱和判断
if ts-t1-t2<0
t1=ts*t1/(t1+t2)
t2=ts-t1
end
ta=(ts-t1-t2)/4
tb=ta+t1/2
tc=tb+t2/2
switch sectorNO
case 3
T1on=ta
T2on=tb
T3on=tc
case 1
T1on=tb
T2on=ta
T3on=tc
case 5
T1on=tc
T2on=ta
T3on=tb
case 4
T1on=tc
T2on=tb
T3on=ta
case 6
T1on=tb
T2on=tc
T3on=ta
case 2
T1on=ta
T2on=tc
T3on=tb
end
CMPR=[T1on T2on T3on]
本文将介绍到：
2. 为什么最大矢量圆半径是
3. 为什么空间矢量的边界是六边形
4. 各扇区内基本矢量作用时间计算
在本文中，我们先做如下约定：
1. 线电压指电机三相绕组中任意两相绕组端点之间的电压，如 U_{UV} ；
2. 相电压指电机三相绕组中任意一相端点到电机绕组中心点O的电压，如 U_{UO} ；
3. 端电压指电机三相绕组中任意一相端点相对于GND的电压，如 U_{UG} ；
4. 中性点电压指电机三相绕组的中性点相对于GND的电压，如 U_{OG} ；
图1
大多数人最初接触电机控制时了解的可能是两两导通三相六状态的方波控制方式。这种控制方式根据电机的位置来选择导通不同桥臂上下各一个MOS管，理解起来比较容易。SVPWM与这种控制方式不同，SVPWM调制方式是将UVW三个桥臂分别定义为0,1两种状态，规定：
0代表下管开通上管关断；
1代表上管开通下管关断；
至于为什么这样，可能要去翻一番SVPWM的发展历史了，我们就姑且认为有人发现了这么一种控制方式，将桥臂按照如上所述定义。
三个桥臂的两种状态总共有八个组合，产生的结果如下：
这里为了便于记忆，我们将3位二进制与十进制的0-7对应起来。如下图，第一个桥臂是W相，它的数值为0，说明它的下MOS是导通的。第二个桥臂是V相，它的数值为1，说明它的上MOS是导通的。第三个桥臂是U相，它的数值是0，说明它的下MOS是导通的。这里与大家平时习惯的UVW(有些文献喜欢用ABC)顺序相反，但是后面所有的推导都是以U轴中心线和X轴正方向一致，V轴中心线在逆时针120度位置，W轴中心线在逆时针240度位置，所以结果是普遍适用的。采用此方法是方便记忆六扇区的矢量合成图分布。
图2
八个导通状态如下：
000 - v0 (zero vector)
001 - v1 (Phase +U)
010 - v2(Phase +V)
011 - v3 (Phase -W)
100 - v4(Phase +W)
101 - v5 (Phase -V)
110 - v6(Phase -U)
111 - v7 (zero vector)
下面通过电机控制的简图说明以上8种状态。这里用双向导通开关替代了上桥臂MOS和下桥臂MOS。当开关连接到上端，则上MOS导通。当开关连接到下端，则下MOS导通。
图3 000 - v0 (zero vector)
图4 001 - v1 (Phase +U)
图5 010 - v2(Phase +V)
图6 011 - v3 (Phase -W)
图7 100 - v4(Phase +W)
图8 101 - v5 (Phase -V)
图9 110 - v6(Phase -U)
图10 111 - v7 (zero vector)
熟悉了以上八种导通状态后，接下来分析各通电状态下，各相电压的大小以及合成电压矢量大小及方向。
如下图，U相桥臂的上管开通，下管关闭。V相和W相桥臂的上管关闭，下管导通。现规定电流从各相端口流入中性点为电流正方向，从中性点流出端口为负方向。则U相的电流为正，V/W相的电流为负。假设电机各相的阻抗相等，都为Z，则根据基尔霍夫电流定律，流过U相的电流是流过V/W相电流的2倍。计算各端口相对电机中性点的电压可以得到各相相电压如图所示，这里计算的中性点电压是相对于GND的，其他相电压都是端电压相对中性点的压降。根据右手定则，可得各相的磁场方向。同理，可以推广到电压矢量，如图右下角所示，可以计算得到三相电压合成矢量的数值为Udc，方向如图所示。为了便于理解相电压的计算，在本图中将右侧的星形电机内部绕组模型也与左边的MOS电路进行了连接。其实，左边已经有了电机绕组，这里大家清楚是为了表达而加的连接线即可。
图 11 001 - v1 (Phase +U)
关于其他导通状态与上图类似，这里直接给出结果，省略了文字描述。
图12 010 - v2(Phase +V)
图13 011 - v3 (Phase -W)
图14 100 - v4(Phase +W)
图15 101 - v5 (Phase -V)
图16 110 - v6(Phase -U)
以上六张图包含着大量的信息，我们来梳理一下。
线电压和相电压
这里提到的相电压是各端点相对于Y绕组中心的电压，线电压是两个端点之间的电压。可总结如下图：
图17
2. 8种通电状态下的基本矢量位置
如下图，6个导通状态下的电压矢量，加上000和111两个0矢量。为了快速记忆这张图各矢量的位置，可以将V1矢量与电机三相坐标系的U轴轴线正方向一致，则在120度位置的V轴与V2(010)重合，正好只有中间代表V桥臂的2进制为1(参见图2)。同理，在240度位置的W轴与V4（100）重合，因为100只有代表W相的2进制为1(参见图2)。两个相邻基本相量的间隔是60度，则剩余的空间矢量(相量)很好确定位置。比如，60度位置的相量就是V1和V2两个二进制的和，即：001 + 010 = 011 = V3。同理，可以在不参考任何资料的情况下快速将各基本矢量画出来。基本矢量的长度为Udc，上面图11-16中均有计算。
图18：
但是，如果按照这么来控制电机就类似于方波BLDC电机的6步换相控制了，每60电角度度切换一次导通状态，容易引起扭矩波动。尽量减少两个矢量切换之间的角度可以降低扭矩波动。为此，我们把60电角度的区域称之为一个扇区，如下图所示。
图19
如下图，在每个扇区内部，我们让矢量尽量均匀旋转，这样电机的扭矩波动就会降低。那好，怎么来实现60度内的矢量分布？
图20
如下图，任取第一扇区内的任一矢量，根据平行四边形原则，可以将其分解为沿V1和V2的量。也就是说我们可以根据相邻的两个基本矢量和一个或两个矢量(000或者111)计算它们所围成的扇区内的任意矢量的幅值和方向。这里，T代表的是采样周期，一般是PWM的周期，V1的作用时间是T1，V2的作用时间是T2，则0矢量的作用时间为T-T1-T2。
图21
到此，我们也就引出了SVPWM控制的基本思想。
通过第二节我们建立了UVW三相坐标下的空间矢量图(注意这里不是 alphabeta 坐标系)，其基本矢量的模长为Udc。接下来，我们再来思考一个问题，在每个扇区内，由两个基本矢量和零矢量合成的最大矢量圆半径是多少？为了推导这个问题，我们拿第一扇区为例进行说明。如下图，单独取出第一扇区并将其放大，第一扇区的两个基本矢量是V1和V3，零矢量可以为V000或者V111。
图22
可知，在30电角度时若要取得最大的电压合成矢量，则T1 = T2 = T/2(T0矢量作用时间为0)。根据以上信息，可以画出最大矢量圆，如下图绿色圆，所有的电压矢量要落在该矢量圆内部，以便进行线性调制。超出该圆后就成为过调制，有对应的处理办法。
图23
当超出该最大矢量圆半径后，在某些角度下，计算得到的T0作用时间为负值，如下图所示。
图24
有关空间矢量分析的论文和网络文章很多，但是极少有深入介绍为什么空间矢量的边界是六边形。大部分都是在推导出基本矢量后就直接默认给出这个六边形，并直接给出结论说合成矢量不能超越这个六边形边界，没有给出理论公式，也没有直观叙述。
在看文献时，这个问题困惑了我很久。现在就将自己的一点直观理解放在这里，作为抛砖引玉，希望有大神可以出来给大家更好的解释。
SVPWM是用伏秒平衡的原理来用基本矢量等效平均合成矢量。六边形出现在每60度电角度内不进行PWM调制(six-step mode)或者说出现在当整个扇区内T0都等于0的情况。
当不进行PWM调制时，每隔60度电角度合成矢量变换一次，类似于方波控制的六步换相。以第一扇区为例，本来合成矢量的大小是V1,并且持续大约60度电角度，不进行任何MOS开关的切换。当转过大约60度电角度时需要切换为合成矢量V2，这时存在一个切换过程，这个过程中V1逐渐减小，V2逐渐增加。它们的合成矢量便沿着六边形在第一扇区的边进行过渡。整个过程始终没有T0的参与。
图25
使用MATLAB编写程序，计算这个过程的合成矢量大小，便得到如图所示的六边形。
图26
接下来，换下脑子，讨论 另外一个问题：当我知道了合成电压矢量，如何计算各基本矢量的作用时间？这涉及到如何来设置PWM寄存器的数值。
关于各扇区相邻两基本矢量及零矢量的作用时间可以参考如下文章。里面有给出了每一步的计算公式，并配有图片方便理解。
详细推导SVPWM各扇区矢量作用时间
后续计划更新如下内容：
调制系数M
确定合成矢量所在扇区的方法
马鞍波形成探讨
SVPWM不同调制波形
SPWM注入零序矢量与SVPWM波形对比
等等。。。
大致弄懂了SVPWM的基本原理之后，下面我就使用MATLAB SIMULINK工具进行仿真
整体的仿真模块如下图所示：
整体的simulink仿真图如上图所示，虽然很多资料上都有，但是在实际按照资料搭建的过程中总是会遇到这样或者是那样的问题，在我经过多次失败的过程中，终于搭建出来正确的模型。下面我就具体介绍一下每个子模块的具体搭建方法以及参数的设定。
在本模型中PWM的开关周期为0.0002s,直流侧电压Udc=700V,仿真算法采用变步长ode23tb算法，且最大仿真步长（MAX STEP SIZE）为0.00001，其余值保持初始值不变。
PWM三角波可以用MATLAB SIMULINK仿真模块中的repeating sequence产生，具体设置如下：
其中对称的等腰三角波的幅值一般取频率的一般，使用示波器观察三角波如下：
需要注意的一点是，注意设定输入正弦波的采样时间，如果采样时间过大或者过小都会使得这个三角波失真。这里给出输入的两个正弦波的模块的设置：
接下来就分模块说说如何搭建：
首先是扇区N的判断模块：
扇区N的计算结果如下：
接下来是中间变量X,Y,Z的计算模块
T1和T2计算模块
切换时间Tcm1,Tcm2,Tcm3模块
切换时间的计算结果：
最后是计算相电压的模块，一般资料这个模块都没有把内部的结构详细给出
相电压Ua的计算结果：
就这样一个完整的SVPWM模块就完成了~
欢迎分享，转载请注明来源：内存溢出
微信扫一扫
支付宝扫一扫
评论列表（0条）