我是靠谱客的博主 任性裙子,最近开发中收集的这篇文章主要介绍【电机控制算法】SVPWM算法的应用(基于Simulink与STM32外设场景的仿真)前言一、扇区判断二、各扇区基础矢量作用时间的计算三、STM32的外设场景四、互补输出且呈中心对称的PWM对五、载波和调制波六、观察相应扇区的3对PWM波形七、Simulink编程仿真,觉得挺不错的,现在分享给大家,希望可以做个参考。

概述

文章目录

  • 前言
  • 一、扇区判断
  • 二、各扇区基础矢量作用时间的计算
    • 1.若处于扇区I
    • 2.若处于扇区II
    • 3.若处于扇区III
    • 4.若处于扇区IV
    • 5.若处于扇区V
    • 6.若处于扇区VI
    • 7.In a word
  • 三、STM32的外设场景
    • 1. TIMx的计数器CNT
    • 2. TIMx的输入捕获寄存器CCRx
    • 3. TIMx的两种PWM模式
  • 四、互补输出且呈中心对称的PWM对
  • 五、载波和调制波
  • 六、观察相应扇区的3对PWM波形
    • 1.扇区I
    • 2.扇区II
    • 3.扇区III
    • 4.扇区IV
    • 5.扇区V
    • 6.扇区VI
    • 7.In a word
  • 七、Simulink编程仿真
    • 1. 输入与输出
    • 2. MATLAB编程
    • 3. 利用SVPWM算法的输出搭建中心对称且互补的PWM生成器
    • 4.仿真结果与波形


前言

授人以鱼不如授人以渔,笔者没有直接给出MCU的源码,而是希望读者能够触类旁通,能在仿真软件中搭建硬件的机制,又能根据硬件的特性编程出SVPWM。
笔者写下此文的目的,另一方面也是为了记录在MCU中实现七段式SVPWM最简单最高效率的计算方法,所以本文并不会把重点放在理解SVPWM上,而是适当推导,利用最终结论,并结合STM32的硬件场景,去编写SVPWM算法。
若您的重点是在于理解SVPWM的由来和始末,相信其他文章能够讲得比笔者更好。
在这里插入图片描述

一、扇区判断

我们需要知道电机转子当前所在的扇区,下一步才能生成与转子形成一定相位差的定子磁场,吸引转子以最大的力矩旋转。这时候需要利用到两相静止坐标轴(2s)的电压值 Vα和Vβ.
我们定义:

在这里插入图片描述

不难发现,U1+U2+U3=0。在芯片中编程时,可以使用这个关系对U2或者U3的表达式进行简化,从而避免更多的乘除法运算(因为更耗时)。
再定义:
若U1>0,A=1;否则A=0;
若U2>0,B=1;否则B=0;
若U3>0,C=1;否则C=0;
存在关系式: N=4C+2B+A
真值表:

N315462
扇区sectorIIIIIIIVVVI

二、各扇区基础矢量作用时间的计算

对于矢量,都是可以合成的。定义Uref为合成矢量,Ts为一个作用周期(也就是一个中心对称的互补PWM对的周期):

1.若处于扇区I

在这里插入图片描述
在这里插入图片描述

2.若处于扇区II

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

3.若处于扇区III

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

4.若处于扇区IV

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

5.若处于扇区V

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

6.若处于扇区VI

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

7.In a word

在这里插入图片描述

扇区sectorIIIIIIIVVVI
作用时间1T4=K*U2T2=-K*U2T2=K*U1T1=-K*U1T1=K*U3T4=-K*U3
作用时间2T6=K*U1T6=-K*U3T3=K*U3T3=-K*U2T5=K*U2T5=-K*U1

若Ts>作用时间1+作用时间2,那么,剩下的作用时间都平均分配给了两个零矢量。T(000)=T(111)=0.5 *(Ts-作用时间1-作用时间2)
当Ts=作用时间1+作用时间2,T(000)=T(111)=0;

三、STM32的外设场景

关于定时器PWM输出机制的理解上,还请翻阅:STM32定时器的单路PWM输出。

1. TIMx的计数器CNT

<1>非中心对称计数
在这里插入图片描述
<2>中心对称计数
在这里插入图片描述

2. TIMx的输入捕获寄存器CCRx

CCRx是指:当CNT计数到CCRx寄存器存储的计数值时,CHx的电平发生转换,有个图更形象:
在这里插入图片描述

3. TIMx的两种PWM模式

STM32主要有两种PWM发波模式。
PWM模式1:
CNT<CCR,CHx高电平,否则CHx低电平
在这里插入图片描述

PWM模式2:
CNT>CCR,CHx高电平,否则CHx低电平
在这里插入图片描述

四、互补输出且呈中心对称的PWM对

以CNT向上计数,PWM2模式为例。我们看看外设的应用场景。
单路非中心对称的PWM输出:
在这里插入图片描述
而互补输出且呈中心对称的PWM:
在这里插入图片描述
在这里插入图片描述
一般的MOS管,它的关断时间大于开通时间很多,反向恢复时间则更长,编程时往往选最大的时延考虑。

三相六臂全桥

三相六臂全桥

在这里插入图片描述
在这里插入图片描述

五、载波和调制波

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

六、观察相应扇区的3对PWM波形

以下均以CNT向上计数,PWM2模式为例.
在这里插入图片描述

波形由上到下依次是A,B,C三相。

1.扇区I

在这里插入图片描述
更正:Tx定义为零矢量后先发送的矢量的作用时间,Ty定义为其后发生的矢量的作用时间.

2.扇区II

在这里插入图片描述

3.扇区III

在这里插入图片描述

4.扇区IV

在这里插入图片描述

5.扇区V

在这里插入图片描述

6.扇区VI

在这里插入图片描述

7.In a word

在这里插入图片描述
对于CNT向上计数,PWM2模式:.

扇区sectorIIIIIIIVVVI
CCR1TaTbTcTcTbTa
CCR2TbTaTaTbTcTc
CCR3TcTcTbTaTaTb

【拓展】
CNT向上计数,PWM模式1?
在这里插入图片描述

CNT向下计数,PWM模式1?
在这里插入图片描述

CNT向下计数,PWM模式2?
在这里插入图片描述

这里留给读者自行推理,只要明白了例子,很快就能理解CCRx在不同模式下正确的赋值

七、Simulink编程仿真

总封装:
在这里插入图片描述
点进去里面,逐一分析。

1. 输入与输出

需要的输入
SVPWM的输入是反Park变换中得到的 Vα和Vβ,这两个是最关键的。
次之从我们的公共乘积因子K可知,我们需要输入母线的电压 Udc .
最后,输入我们希望生成的调制波频率 f_PWM .
我们得到的输出
三路输出比较寄存器CCRx的值
扇区
这样,Simulink中Svpwm模组中的函数function就已经确定了:
在这里插入图片描述
假如我们要求中心对称PWM对的频率是8K,那么他非中心对称的PWM(单路PWM)的频率就是2 * f_PWM,那么,中心对称的PWM到达中点的时间周期就是1/(2 * f_PWM);

2. MATLAB编程

采用载波向上计数,PWM模式2的思想进行搭建,编程。

%4个输入,4个输出
function [Tcm1,Tcm2,Tcm3,sector] = fcn( U_alpha,U_beta,U_dc,PWM_Freq)

%输出参数初始化
sector=0;
Tcm1=0;
Tcm2=0;
Tcm3=0;

%PWM采用中心对齐
%为了避免值过小,默认使用PWM频率的两倍作为载波值
%假如我们希望中心对称PWM的频率是8K,那么他对应的单路PWM频率应该是16K
%所以我们载波CNT的计数最大值是16K
Ts=PWM_Freq;
sqrt_3=sqrt(3);

%判断扇区部分(需要强制类型转换为浮点型)
U1=double(U_beta);
U2=double(sqrt_3*0.5*U_alpha-U_beta*0.5);
U3=double(-sqrt_3*0.5*U_alpha-U_beta*0.5);
%获得对应布尔值
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;
%根据N的公式结果判断扇区
switch N
    case 3
        sector=1;
    case 1
        sector=2;
    case 5
        sector=3;
    case 4
        sector=4;
    case 6
        sector=5;
    case 2
        sector=6;
end
%不同扇区不同矢量的作用时间计算
K=sqrt_3*Ts/U_dc;%公共乘积因子
switch (sector)
    case 1
        Tx=K*U2;
        Ty=K*U1;
    case 2
        Tx=-K*U2;
        Ty=-K*U3;
    case 3
        Tx=K*U1;
        Ty=K*U3;
    case 4
        Tx=-K*U1;
        Ty=-K*U2;
    case 5
        Tx=K*U3;
        Ty=K*U2;
    case 6
        Tx=-K*U3;
        Ty=-K*U1;
    otherwise
        Tx=0;
        Ty=0;
end
%过调制算法
if (Tx+Ty)>Ts
    Tx=Ts*Tx/(Tx+Ty);
    Ty=Ts*Ty/(Tx+Ty);
else
    Tx=Tx;
    Ty=Ty;
end
%零矢量总作用时间的二分之一,也就是T0和T7的值
T0=(Ts-Tx-Ty)*0.5;
%见六中的总结
Ta=T0*0.5;      %先跳转的PWM对的输出比较值
Tb=Ta+0.5*Tx;   %后跳转的PWM对的输出比较值
Tc=Tb+0.5*Ty;   %最后跳转的PWM对的输出比较值

%根据六中的扇区情况,赋值三路PWM对的输出比较值
switch sector
    case 1
        Tcm1=Ta;
        Tcm2=Tb;
        Tcm3=Tc;
    case 2
        Tcm1=Tb;
        Tcm2=Ta;
        Tcm3=Tc;
    case 3
        Tcm1=Tc;
        Tcm2=Ta;
        Tcm3=Tb;
    case 4
        Tcm1=Tc;
        Tcm2=Tb;
        Tcm3=Ta;
    case 5
        Tcm1=Tb;
        Tcm2=Tc;
        Tcm3=Ta;
    case 6
        Tcm1=Ta;
        Tcm2=Tc;
        Tcm3=Tb;
end
end

3. 利用SVPWM算法的输出搭建中心对称且互补的PWM生成器

总图:
在这里插入图片描述
需要的输入
三路输出比较寄存器CCRx的值
模拟CNT计数的三角波作为载波
获得的输出
三路互补且中心对称的PWM波,用于控制开关管

寻找载波模块:
在这里插入图片描述
拖出来,双击调整参数:
在这里插入图片描述

载波从0到1/(2 * f_PWM)时间,到达最大的计数值f_PWM。对应的是中心对称PWM周期的时间中点,即0.5Ts。
笔者这里使用数字量去代表时间,数字量拟用PWM频率代替,比如8K的PWM频率,它的Ts=2
PWM_Freq,那么它的时间中点就是f_PWM,这个也类似STM32的重装载值(最大65535);

点进PWM Generator
以下是笔者搭建的内容
在这里插入图片描述

4.仿真结果与波形

黄色表示A相
绿色表示B相
蓝色表示C相
CCRx的波形是马鞍波,而调制波直接输出到PMSM电机的三端上,所以端电压(端对地GND的电压,不是相电压)也是马鞍波!
在这里插入图片描述
中心对称的PWM对,作用在相上,会使用数字正弦波来等效模拟正弦波;
当然这个仿真可以通过CCRx的值作处理来获得,使用以下公式:

Un_A=double(Udc/3*(2*PWM_a-PWM_b-PWM_c));
Un_B=double(Udc/3*(2*PWM_b-PWM_a-PWM_c));
Un_C=double(Udc/3*(2*PWM_c-PWM_a-PWM_b));

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
这种数字正弦波,作用在电机上,感性负载的电流惯性,使电机产生与数字正弦波电压有一定相角差的正弦电流,同时三相的电流也是120°对称的。因为是数字正弦波,在正弦电流上我们可以看到一丝丝纹波,这是正常的。
作用在电机上的相电流:
在这里插入图片描述

最后

以上就是任性裙子为你收集整理的【电机控制算法】SVPWM算法的应用(基于Simulink与STM32外设场景的仿真)前言一、扇区判断二、各扇区基础矢量作用时间的计算三、STM32的外设场景四、互补输出且呈中心对称的PWM对五、载波和调制波六、观察相应扇区的3对PWM波形七、Simulink编程仿真的全部内容,希望文章能够帮你解决【电机控制算法】SVPWM算法的应用(基于Simulink与STM32外设场景的仿真)前言一、扇区判断二、各扇区基础矢量作用时间的计算三、STM32的外设场景四、互补输出且呈中心对称的PWM对五、载波和调制波六、观察相应扇区的3对PWM波形七、Simulink编程仿真所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(50)

评论列表共有 0 条评论

立即
投稿
返回
顶部