SVPWM调制的simulink仿真

1、SVPWM的生成

在FOC矢量控制中,Id、Iq 经过PID输出Vd、Vq。Vd、Vq经过反park变换成Vα、Vβ。再Vα、Vβ合成空间参考矢量Uref。那么怎么根据Vα、Vβ确定Uref所在扇区,然后确定所在扇区两个非零基础矢量的作用时间呢。

1.1、传统的计算方法

1.1.1、扇区的判断

把Uref转换为虚数形式Uref=Va+jVβ。那么其所在扇区 sector = floor(angle(Uref)/(pi/3)) + 1。

1.1.2、两个非零基础矢量的作用时间

那么其所在扇区的两个基础非零矢量的作用时间呢?推导下,假设已确定在第一扇区。
在这里插入图片描述
在两相静止参考坐标系(α,β)中,令 Uref 和 U4 间的夹角是 θ,由正弦定理可得:
在这里插入图片描述
因为 U4 = 2*Udc / 3,所以可以得到两个非零基础矢量的作用时间为
在这里插入图片描述
式中 m 为 SVPWM 调制系数(调制比)注意超过1就失真,m=sqrt(3)|Uref|/Udc。可以看出计算两个相邻基础向量的时间需要解三角函数,对单片机来说计算有点费时间。本来一个核心的FOC算法要求在十几二十个微秒内执行完。那么下面介绍一种不需要解三角函数的方法。

1.2、第二种方法

1.2.1、扇区的判断
在这里插入图片描述
由Vα、Vβ决定的空间电压矢量所处的扇区。假定合成的电压矢量落在第 I 扇区,可知其等价条件如下:
0° < arctan( Vβ/Vα ) < 60° 以上等价条件再结合矢量图几何关系分析,可以判断出合成电压矢量 Uref
落在第 I 扇区的充分必要条件为:Vα>0、Vβ>0且Vβ/Vα<sqrt(3)。同理可得:

在这里插入图片描述
若进一步分析以上的条件,有可看出参考电压矢量 Uref 所在的扇区完全由 U β U_{\beta } Uβ 3 U α − U β \sqrt{3}U_{\alpha }-U_{\beta} 3 UαUβ − 3 U α − U β -\sqrt{3}U_{\alpha }-U_{\beta} 3 UαUβ 三式决定,因此令:
{ U 1 = U β U 2 = 3 2 U α − 1 2 U β U 3 = − 3 2 U α − 1 2 U β \begin{matrix} \left\{ \begin{aligned} &U_{1}=U_{\beta}\\ &U_{2}=\frac{\sqrt{3}}{2}U_{\alpha}-\frac{1}{2}U_{\beta}\\ &U_{3}=-\frac{\sqrt{3}}{2}U_{\alpha}-\frac{1}{2}U_{\beta}\\ \end{aligned} \right. \end{matrix} U1=UβU2=23 Uα21UβU3=23 Uα21Uβ
再定义,若U1 > 0 ,则 A=1,否则 A=0;若U 2 > 0 ,则 B=1,否则 B=0;若U3 > 0,则 C=1,否则 C=0。
可以看出 A,B,C 之间共有八种组合,但由判断扇区的公式可知 A,B,C不会同时为 1 或同时为 0,所以实际的组合是六种,A,B,C 组合取不同的值对应着不同的扇区,并且是一一对应的,因此完全可以由 A,B,C 的组合判断所在的扇区。为区别六种状态,令 N=4C+2B+A,则可以通过下表计算参考电压矢量 Uref 所在的扇区。

N315462
扇区号

所以只要简单的逻辑运算就可以确定所在扇区,对于提高系统的响应速度和进行仿真都是很有意义的。

1.2.2、两个非零基础矢量的作用时间

根据伏秒平衡原理,以 Uref 处在第Ⅰ扇区时进行分析,根据图 2-9 有:

假设,电压参考矢量Uref在Uα、Uβ坐标轴上的投影为Uα、Uβ,Uref在两个相邻的非零基础矢量上的投影为 2 3 U d c T 4 T s \frac{2}{3}U_{dc} \frac{T4}{Ts} 32UdcTsT4 2 3 U d c T 6 T s \frac{2}{3}U_{dc} \frac{T6}{Ts} 32UdcTsT6,如果Uref在U4、U6上的投影再投影到Uα、Uβ坐标系上,则:

在这里插入图片描述
经过整理后得出:
在这里插入图片描述
表中两个非零矢量作用时间的比例系数为 K = 3 T s / U d c \sqrt{3}T_{s}/U_{dc} 3 Ts/Udc,注意Uref的幅值不能大于 U d c / 3 U_{dc}/\sqrt{3} Udc/3 防止Uref在正六边形内切圆外面失真。

如果Uref 处在第II扇区,根据图 2-10 有:
在这里插入图片描述
经过整理后得出:
在这里插入图片描述
同理可求得Uref在其它扇区中各矢量的作用时间,结果如表2-4所示。

扇区时间
T 4 = K U 2 、 T 6 = K U 1 T_{4}=KU_{2}、T_{6}=KU_{1} T4=KU2T6=KU1
T 2 = − K U 2 、 T 6 = − K U 3 T_{2}=-KU_{2}、T_{6}=-KU_{3} T2=KU2T6=KU3
T 2 = K U 1 、 T 3 = K U 3 T_{2}=KU_{1}、T_{3}=KU_{3} T2=KU1T3=KU3
T 1 = − K U 1 、 T 3 = − K U 2 T_{1}=-KU_{1}、T_{3}=-KU_{2} T1=KU1T3=KU2
T 1 = K U 3 、 T 5 = K U 2 T_{1}=KU_{3}、T_{5}=KU_{2} T1=KU3T5=KU2
T 4 = − K U 3 、 T 5 = − K U 1 T_{4}=-KU_{3}、T_{5}=-KU_{1} T4=KU3T5=KU1

1.1.3、三相占空比的计算

所在扇区的非零基础矢量作用时间已确定了,那么怎么计算各相占空比呢。
在这里插入图片描述
上图为在第一扇区时各相占空比。
{ T a = ( T s − T 4 − T 6 ) / 4 T b = T a + 1 2 T 4 T c = T b + 1 2 T 6 \begin{matrix} \left\{ \begin{aligned} &T_{a}=(T_{s}-T_{4}-T_{6})/4\\ &T_{b}=T_{a}+\frac{1}{2}T_{4}\\ &T_{c}=T_{b}+\frac{1}{2}T_{6} \end{aligned} \right. \end{matrix} Ta=TsT4T6/4Tb=Ta+21T4Tc=Tb+21T6
或者
{ T a = ( T s + T 4 + T 6 ) / 4 T b = T a − 1 2 T 4 T c = T b − 1 2 T 6 \begin{matrix} \left\{ \begin{aligned} &T_{a}=(T_{s}+T_{4}+T_{6})/4\\ &T_{b}=T_{a}-\frac{1}{2}T_{4}\\ &T_{c}=T_{b}-\frac{1}{2}T_{6} \end{aligned} \right. \end{matrix} Ta=Ts+T4+T6/4Tb=Ta21T4Tc=Tb21T6
在这里插入图片描述
上图为在第二扇区时各相占空比。
{ T b = ( T s − T 2 − T 6 ) / 4 T a = T b + 1 2 T 2 T c = T a + 1 2 T 6 \begin{matrix} \left\{ \begin{aligned} &T_{b}=(T_{s}-T_{2}-T_{6})/4\\ &T_{a}=T_{b}+\frac{1}{2}T_{2}\\ &T_{c}=T_{a}+\frac{1}{2}T_{6} \end{aligned} \right. \end{matrix} Tb=TsT2T6/4Ta=Tb+21T2Tc=Ta+21T6
或者
{ T b = ( T s + T 2 + T 6 ) / 4 T a = T b − 1 2 T 2 T c = T a − 1 2 T 6 \begin{matrix} \left\{ \begin{aligned} &T_{b}=(T_{s}+T_{2}+T_{6})/4\\ &T_{a}=T_{b}-\frac{1}{2}T_{2}\\ &T_{c}=T_{a}-\frac{1}{2}T_{6} \end{aligned} \right. \end{matrix} Tb=Ts+T2+T6/4Ta=Tb21T2Tc=Ta21T6
在这里插入图片描述
上图为在第三扇区时各相占空比。
{ T b = ( T s − T 2 − T 3 ) / 4 T c = T b + 1 2 T 2 T a = T c + 1 2 T 3 \begin{matrix} \left\{ \begin{aligned} &T_{b}=(T_{s}-T_{2}-T_{3})/4\\ &T_{c}=T_{b}+\frac{1}{2}T_{2}\\ &T_{a}=T_{c}+\frac{1}{2}T_{3} \end{aligned} \right. \end{matrix} Tb=TsT2T3/4Tc=Tb+21T2Ta=Tc+21T3
或者
{ T b = ( T s + T 2 + T 3 ) / 4 T c = T b − 1 2 T 2 T a = T c − 1 2 T 3 \begin{matrix} \left\{ \begin{aligned} &T_{b}=(T_{s}+T_{2}+T_{3})/4\\ &T_{c}=T_{b}-\frac{1}{2}T_{2}\\ &T_{a}=T_{c}-\frac{1}{2}T_{3} \end{aligned} \right. \end{matrix} Tb=Ts+T2+T3/4Tc=Tb21T2Ta=Tc21T3
在这里插入图片描述
上图为在第四扇区时各相占空比。
{ T c = ( T s − T 1 − T 3 ) / 4 T b = T c + 1 2 T 1 T a = T b + 1 2 T 3 \begin{matrix} \left\{ \begin{aligned} &T_{c}=(T_{s}-T_{1}-T_{3})/4\\ &T_{b}=T_{c}+\frac{1}{2}T_{1}\\ &T_{a}=T_{b}+\frac{1}{2}T_{3} \end{aligned} \right. \end{matrix} Tc=TsT1T3/4Tb=Tc+21T1Ta=Tb+21T3
或者
{ T c = ( T s + T 1 + T 3 ) / 4 T b = T c − 1 2 T 1 T a = T b − 1 2 T 3 \begin{matrix} \left\{ \begin{aligned} &T_{c}=(T_{s}+T_{1}+T_{3})/4\\ &T_{b}=T_{c}-\frac{1}{2}T_{1}\\ &T_{a}=T_{b}-\frac{1}{2}T_{3} \end{aligned} \right. \end{matrix} Tc=Ts+T1+T3/4Tb=Tc21T1Ta=Tb21T3
在这里插入图片描述
上图为在第五扇区时各相占空比。
{ T c = ( T s − T 1 − T 5 ) / 4 T a = T c + 1 2 T 1 T b = T a + 1 2 T 5 \begin{matrix} \left\{ \begin{aligned} &T_{c}=(T_{s}-T_{1}-T_{5})/4\\ &T_{a}=T_{c}+\frac{1}{2}T_{1}\\ &T_{b}=T_{a}+\frac{1}{2}T_{5} \end{aligned} \right. \end{matrix} Tc=TsT1T5/4Ta=Tc+21T1Tb=Ta+21T5
或者
{ T c = ( T s + T 1 + T 5 ) / 4 T a = T c − 1 2 T 1 T b = T a − 1 2 T 5 \begin{matrix} \left\{ \begin{aligned} &T_{c}=(T_{s}+T_{1}+T_{5})/4\\ &T_{a}=T_{c}-\frac{1}{2}T_{1}\\ &T_{b}=T_{a}-\frac{1}{2}T_{5} \end{aligned} \right. \end{matrix} Tc=Ts+T1+T5/4Ta=Tc21T1Tb=Ta21T5

在这里插入图片描述
上图为在第六扇区时各相占空比。
{ T a = ( T s − T 4 − T 5 ) / 4 T c = T a + 1 2 T 4 T b = T c + 1 2 T 5 \begin{matrix} \left\{ \begin{aligned} &T_{a}=(T_{s}-T_{4}-T_{5})/4\\ &T_{c}=T_{a}+\frac{1}{2}T_{4}\\ &T_{b}=T_{c}+\frac{1}{2}T_{5} \end{aligned} \right. \end{matrix} Ta=TsT4T5/4Tc=Ta+21T4Tb=Tc+21T5
或者
{ T a = ( T s + T 4 + T 5 ) / 4 T c = T a − 1 2 T 4 T b = T c − 1 2 T 5 \begin{matrix} \left\{ \begin{aligned} &T_{a}=(T_{s}+T_{4}+T_{5})/4\\ &T_{c}=T_{a}-\frac{1}{2}T_{4}\\ &T_{b}=T_{c}-\frac{1}{2}T_{5} \end{aligned} \right. \end{matrix} Ta=Ts+T4+T5/4Tc=Ta21T4Tb=Tc21T5

上式中 T a 、 T b 、 T c 、 T s T_{a}、T_{b}、T_{c}、T_{s} TaTbTcTs表示A、B,C三相的占空比,和PWM周期。如果把上式写成一个通用的式子。
{ t a = ( T s − T x − T y ) / 4 t b = ( t a + 1 2 T x ) t c = ( t b + 1 2 T y ) \begin{matrix} \left\{ \begin{aligned} &t_{a}=(T_{s}-T_{x}-T_{y})/4\\ &t_{b}=(t_{a}+\frac{1}{2}T_{x})\\ &t_{c}=(t_{b}+\frac{1}{2}T_{y})\\ \end{aligned} \right. \end{matrix} ta=TsTxTy/4tb=ta+21Tx)tc=tb+21Ty)
或者
{ t a = ( T s + T x + T y ) / 4 t b = ( t a − 1 2 T x ) t c = ( t b − 1 2 T y ) \begin{matrix} \left\{ \begin{aligned} &t_{a}=(T_{s}+T_{x}+T_{y})/4\\ &t_{b}=(t_{a}-\frac{1}{2}T_{x})\\ &t_{c}=(t_{b}-\frac{1}{2}T_{y})\\ \end{aligned} \right. \end{matrix} ta=Ts+Tx+Ty/4tb=ta21Tx)tc=tb21Ty)
其中下标 y 大于 x 。ta、tb 和 tc分别是相应的比较器的值,而不同扇区比较器的值分配如表下所示:

扇区123456
T a T_{a} Ta t a t_{a} ta t b t_{b} tb t b t_{b} tb t c t_{c} tc t c t_{c} tc t a t_{a} ta
T b T_{b} Tb t b t_{b} tb t a t_{a} ta t c t_{c} tc t b t_{b} tb t a t_{a} ta t c t_{c} tc
T c T_{c} Tc t c t_{c} tc t c t_{c} tc t a t_{a} ta t a t_{a} ta t b t_{b} tb t b t_{b} tb

2、SVPWM的simulink仿真

下面根据第二种方法来仿真:

2.1、simulink仿真波形图
在这里插入图片描述
在这里插入图片描述

2.2、simulink仿真程序

程序中去掉了Udc,然后也限制了过调制失真。

function [Ta,Tb,Tc] = fcn(ua,ub,udc,Ts)

u1 = ub;
u2 = sqrt(3)/2*ua - 1/2*ub;
u3 = -sqrt(3)/2*ua - 1/2*ub;

if u1 > 0
    a = 1;
else 
    a = 0;
end

if u2 > 0
    b = 1;
else 
    b = 0;
end

if u3 > 0
    c = 1;
else 
    c = 0;
end

N = 4*c + 2*b + a;

Ts = Ts*2; %因为PWM是中心对齐模式

k = sqrt(3)*Ts/udc; 

switch N
    case 1  %2扇区
        t1 = -k*u2;
        t2 = -k*u3;
        
        Tb = ((Ts)+ t1 + t2)/4;      
        Ta = Tb - t1/2;     
        Tc = Ta - t2/2; 
          
    case 2  %6扇区
        t1=-k*u3;
        t2=-k*u1;
        
        Ta = ((Ts)+ t1 + t2)/4;      
        Tc = Ta - t1/2;     
        Tb = Tc - t2/2; 
        
    case 3  %1扇区
        t1=k*u2;
        t2=k*u1;
               
        Ta = ((Ts)+ t1 + t2)/4;      
        Tb = Ta - t1/2;     
        Tc = Tb - t2/2; 
        
    case 4  %4扇区
        t1=-k*u1;
        t2=-k*u2;
              
        Tc = ((Ts)+ t1 + t2)/4;      
        Tb = Tc - t1/2;     
        Ta = Tb - t2/2; 
        
    case 5  %3扇区
        t1=k*u1;
        t2=k*u3;
                
        Tb = ((Ts)+ t1 + t2)/4;      
        Tc = Tb - t1/2;     
        Ta = Tc - t2/2; 
        
    case 6  %5扇区
        t1=k*u3;
        t2=k*u2;
        
        Tc = ((Ts) + t1 + t2)/4;      
        Ta = Tc - t1/2;     
        Tb = Ta - t2/2; 
        
    otherwise 
        Ta = Ts/4;      
        Tb = Ts/4;     
        Tc = Ts/4; 

end

end



参考模型:https://download.csdn.net/download/wanrenqi/35140354

  • 8
    点赞
  • 76
    收藏
    觉得还不错? 一键收藏
  • 6
    评论
三电平SVPWM(Space Vector Pulse Width Modulation)是一种广泛应用在电力电子变流器控制中的调制技术。其通过调整输出电压的脉冲宽度比例,实现对交流电源电压的调节。在Simulink仿真中,可以使用适当的模块和参数设置来实现三电平SVPWM仿真。 首先,我们需要在Simulink环境中建立一个模型。模型中包括一个电压源,用于提供三相交流电源;一个控制器,用于生成SVPWM调制信号;以及一个模拟电力电子变流器,用于将调制信号转换为输出电压。 在控制器模块中,可以使用三电平SVPWM算法生成相应的调制信号。SVPWM算法主要包括三个步骤:1)确定电机的空间矢量;2)通过空间矢量与三相相位角之间的关系,计算三相对应的调制信号的占空比;3)通过调制信号的占空比,生成相应的脉冲信号。 在模拟电力电子变流器模块中,可以使用合适的电路元件和参数设置,将调制信号转换为输出电压。模拟电力电子变流器一般包括逆变器和滤波电路。逆变器将调制信号转换为脉冲电压,而滤波电路则对脉冲电压进行滤波,得到平滑的输出电压。 在仿真过程中,可以通过调整控制器的参数、电力电子变流器的参数以及输入电压的波形等来观察输出电压的变化。可以对三电平SVPWM仿真进行多个场景的测试,例如改变输入电压的频率、改变负载阻抗等。通过仿真,可以评估三电平SVPWM在不同工况下的性能和稳定性。 总之,通过Simulink仿真,可以实现对三电平SVPWM的设计和调试,以及对其在不同工况下的性能评估。这种仿真方法有助于优化电力电子变流器的设计和控制策略,提高系统的效率和可靠性。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 6
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值