CN112803809B - 基于s-k快速确定基本矢量的多电平变换器svm方法 - Google Patents

基于s-k快速确定基本矢量的多电平变换器svm方法 Download PDF

Info

Publication number
CN112803809B
CN112803809B CN202110133184.9A CN202110133184A CN112803809B CN 112803809 B CN112803809 B CN 112803809B CN 202110133184 A CN202110133184 A CN 202110133184A CN 112803809 B CN112803809 B CN 112803809B
Authority
CN
China
Prior art keywords
vector
switch state
zero sequence
absolute value
sequence voltage
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110133184.9A
Other languages
English (en)
Other versions
CN112803809A (zh
Inventor
李常学
王翠
段万政
曾文军
王沄禾
欧阳俊铭
李泽军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanchang Institute of Technology
Original Assignee
Nanchang Institute of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanchang Institute of Technology filed Critical Nanchang Institute of Technology
Priority to CN202110133184.9A priority Critical patent/CN112803809B/zh
Publication of CN112803809A publication Critical patent/CN112803809A/zh
Application granted granted Critical
Publication of CN112803809B publication Critical patent/CN112803809B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02MAPPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M7/00Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
    • H02M7/42Conversion of dc power input into ac power output without possibility of reversal
    • H02M7/44Conversion of dc power input into ac power output without possibility of reversal by static converters
    • H02M7/48Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
    • H02M7/483Converters with outputs that each can have more than two voltages levels
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02MAPPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
    • H02M7/00Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
    • H02M7/42Conversion of dc power input into ac power output without possibility of reversal
    • H02M7/44Conversion of dc power input into ac power output without possibility of reversal by static converters
    • H02M7/48Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
    • H02M7/53Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal
    • H02M7/537Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters
    • H02M7/5387Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters in a bridge configuration

Abstract

本发明公开了一种基于S‑K快速确定基本矢量的多电平变换器SVM方法,通过确定参考矢量所在的大扇区和小扇区,再确定基本矢量以及开关状态,最终形成脉冲信号。此方法将坐标系分为6个大扇区,利用旋转变换将其他5个大扇区的参考矢量定位转化到第Ⅰ大扇区完成,计算量大大减少。与现有技术相比,本发明的方案降低了参考矢量定位的难度,减少了参考矢量定位的计算量;通过S‑K算法快速确定了合成参考矢量的最近三个基本矢量对应的开关状态,通过简单的算法找到了零序电压绝对值最小的开关状态,极大简化了SVM算法的工作量,提高了SVM算法的速度。

Description

基于S-K快速确定基本矢量的多电平变换器SVM方法
技术领域
本发明属于多电平变换器调制领域,涉及一种基于S-K快速确定基本矢量的多电平变换器SVM方法。
背景技术
空间矢量调制(Space Vector Modulation,简称SVM)算法具有直流电压利用率高、谐波畸变率小、易于数字化实现等优点,被广泛应用于多电平变换器的调制。然而随着多电平变换器输出电平数的增加,空间矢量调制算法实现过程中基本矢量数量以及开关状态的数量都大大增加,导致合成参考矢量的基本矢量定位困难,且由于冗余开关状态大量增加,导致基本矢量对应开关状态的选择难度更大,最终制约了空间矢量调制算法在电平数大于七的多电平变换器上的应用。
因此,有必要设计一种基于S-K快速确定基本矢量的多电平变换器SVM 方法。
发明内容
本发明所要解决的技术问题是提供一种基于S-K快速确定基本矢量的多电平变换器SVM方法,该方法的脉冲调制简单,易于实施。
发明的技术解决方案如下:
本发明的目的是针对多电平变换器空间矢量调制算法存在的参考矢量定位计算复杂、繁琐以及定位后基本矢量对应开关状态的选择难度更大的问题,提出了基于S-K快速确定基本矢量的多电平变换器SVM方法,此处所述多电平变换器特指H桥级联多电平变换器。此方法将坐标系分为6个大扇区,利用旋转变换将其他5个大扇区的参考矢量定位转化到第Ⅰ大扇区完成,计算量大大减少。在第Ⅰ大扇区,通过使用一组公式快速定位参考矢量,然后利用S-K算法确定合成参考矢量的基本矢量对应的开关状态。该方法极大减少了传统SVM算法的计算量,简化了级联多电平变换器SVM算法,使得SVM算法适用于n级多电平变换器。三相n级多电平变换器拓扑结构如图1所示,图中每一相由n个电压型单相全桥逆变器子模块串联组成(单相全桥逆变器子模块即H桥),每一个H桥由四个IGBT及反并联二极管组成,H桥并联的电容表明子模块为电压型逆变,且所有子模块的输入电压都相等,用E表示。
为达到上述目的,本发明提供如下技术方案:基于S-K快速确定基本矢量的多电平变换器SVM方法,本方法的流程图如图2所示,具体包括以下步骤:
步骤一:根据SVM算法原理,传统笛卡尔坐标系(简称αβ坐标系)上相电压参考信号ura、urb、urc对应参考矢量Vr(v,v)的表达式为:
Figure GDA0003419104190000021
式中,v和v表示参考矢量Vr对应的坐标分量,相电压参考信号ura、urb、 urc分别为:
Figure GDA0003419104190000022
式中n表示级联单元数,m表示调制系数,ω表示参考电压信号角频率。
αβ坐标系上基本矢量V(α,β)与开关状态U(a,b,c)的映射函数为:
Figure GDA0003419104190000023
式中,α和β表示基本矢量V对应的坐标分量,U表示开关状态名称,a、 b、c分别表示变换器输出的三相相电压电平数,对于每一相由n个H桥级联组成的多电平变换器,a,b,c∈[±n,±(n-1),±(n-2),...,±2,±1,0]。
αβ坐标系上根据相电压参考信号计算的参考矢量Vr的坐标分量也可表示为:
Figure GDA0003419104190000024
步骤二:参考矢量的扇区判断及其旋转变换:
①对参考电压信号ura、urb、urc进行采样,并根据公式(1)计算参考矢量 Vr的坐标分量v和v
②计算参考矢量Vr的相位角θ:
θ=ωt (5)
③计算参考矢量Vr所在扇区的编号N:
Figure GDA0003419104190000031
式中:N表示大扇区的编号,N取值为1、2、3、4、5、6,分别与图3所示第Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ扇区一一对应(图3所示为αβ坐标系上二级五电平变换器参考矢量轨迹和基本矢量的分布图,图中加粗线条组成的三角形表示大扇区,分别标注为Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ扇区)。
④如果参考矢量Vr所在扇区编号N≥2,将参考矢量Vr顺时针旋转 (N-1)π/3角度至第Ⅰ扇区,得矢量Vr',计算矢量Vr'的相位角θ1,如果参考矢量 Vr所在扇区编号N=1,Vr'与Vr为同一个矢量。
Figure GDA0003419104190000032
步骤三:定位矢量Vr'并确定Vr'所在小三角形的编号:如图4所示,分别作第Ⅰ大扇区三角形三条边的垂线OP、OQ、OL,计算参考矢量Vr'在OP、OQ、 OL上的投影并取整得:
Figure GDA0003419104190000033
Figure GDA0003419104190000034
l=floor(2mnsinθ1)+1 (10)
S=-p+q2-q+l+1 (11)
式中:p、q、l分别为矢量Vr'在垂线OP、OQ、OL上的投影并取整所得结果;S为矢量Vr'所在小三角形的编号,以二级五电平变换器为例,其第Ⅰ大扇区小三角形的编号顺序如图5所示,图中共16个小三角形,小三角形分为4组,第一组1个小三角形,第二组3个小三角形,第三组5个小三角形,第四组7 个小三角形,具体的编号顺序为从左下角开始,沿左下角对边的方向、从下向上编号,且依次从1-4组编号。
步骤四:通过S-K算法确定合成参考矢量Vr最近的三个基本矢量V1、V2、V3对应的开关状态U1(a1,b1,c1)、U2(a2,b2,c2)、U3(a3,b3,c3),如图4所示:
①如果
Figure GDA0003419104190000041
Figure GDA0003419104190000042
否则,令
Figure GDA0003419104190000043
②确定包含参考矢量Vr的小三角形其中一个顶点矢量V1,该矢量对应的开关状态U1(a1,b1,c1):
当参考矢量Vr位于第Ⅰ扇区时,即N=1,V1对应的开关状态U1(a1,b1,c1)为:
Figure GDA0003419104190000044
式中ceil(*)表示向上取整函数。
当参考矢量Vr位于第Ⅱ扇区时,即N=2,V1对应的开关状态U1(a1,b1,c1)为:
Figure GDA0003419104190000045
式中floor(*)表示向下取整函数。
当参考矢量Vr位于第Ⅲ扇区时,即N=3,V1对应的开关状态U1(a1,b1,c1)为:
Figure GDA0003419104190000046
当参考矢量Vr位于第Ⅳ扇区时,即N=4,V1对应的开关状态U1(a1,b1,c1)为:
Figure GDA0003419104190000047
当参考矢量Vr位于第Ⅴ扇区时,即N=5时,V1对应的开关状态U1(a1,b1,c1)为:
Figure GDA0003419104190000051
当参考矢量Vr位于第Ⅵ扇区时,即N=6时,V1对应的开关状态U1(a1,b1,c1)为:
Figure GDA0003419104190000052
③根据参考矢量Vr所在小三角形三个顶点矢量之间的关系,判断矢量Vr'所在小三角形的编号S和K+1的奇偶性是否相同,确定基本矢量V2对应的开关状态U2(a2,b2,c2):
N=1时,S和K+1的奇偶性相同:
Figure GDA0003419104190000053
S和K+1的奇偶性不相同:
Figure GDA0003419104190000054
N=2时,S和K+1的奇偶性相同:
Figure GDA0003419104190000055
S和K+1的奇偶性不相同:
Figure GDA0003419104190000056
N=3时,S和K+1的奇偶性相同:
Figure GDA0003419104190000061
S和K+1的奇偶性不相同:
Figure GDA0003419104190000062
N=4时,S和K+1的奇偶性相同:
Figure GDA0003419104190000063
S和K+1的奇偶性不相同:
Figure GDA0003419104190000064
N=5时,S和K+1的奇偶性相同:
Figure GDA0003419104190000065
S和K+1的奇偶性不相同:
Figure GDA0003419104190000066
N=6时,S和K+1的奇偶性相同:
Figure GDA0003419104190000067
S和K+1的奇偶性不相同:
Figure GDA0003419104190000068
④计算基本矢量V3对应的开关状态U3(a3,b3,c3):
当N=1或N=2时,V3对应的开关状态U3(a3,b3,c3):
Figure GDA0003419104190000071
当N=3或N=4时,V3对应的开关状态U3(a3,b3,c3):
Figure GDA0003419104190000072
当N=5或N=6时,V3对应的开关状态U3(a3,b3,c3):
Figure GDA0003419104190000073
步骤五:计算基本矢量V1、V2、V3的开关状态U1'(a1',b1',c1')、U2'(a2',b2',c2')、U3'(a3',b3',c3'),使其对应零序电压绝对值最小:前面计算所得开关状态 U1(a1,b1,c1)、U2(a2,b2,c2)、U3(a3,b3,c3)不一定是零序电压绝对值最小的开关状态,此处是为了计算零序电压绝对值最小的开关状态。
①如果U1(a1,b1,c1)是基本矢量V1唯一的开关状态,则其零序电压也是唯一的,此时V1零序电压绝对值最小的开关状态U1'(a1',b1',c1')=U1(a1,b1,c1),否则:
给定变量i=1,如果:
Figure GDA0003419104190000074
基本矢量V1零序电压绝对值最小的开关状态U1'(a1',b1',c1')=U1(a1,b1,c1),否则:
Figure GDA0003419104190000075
基本矢量V1零序电压绝对值最小的开关状态U1'(a1',b1',c1')=U1'(a1+1,b1+1,c1+1);
或者i∈[2,3,…,n-1],使得:
Figure GDA0003419104190000076
式中,a1+i,b1+i,c1+i,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
a1+i-1,b1+i-1,c1+i-1,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
a1+i+1,b1+i+1,c1+i+1,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
基本矢量V1零序电压绝对值最小的开关状态U1'(a1',b1',c1')=U1'(a1+i,b1+i,c1+i)
②如果U2(a2,b2,c2)是基本矢量V2唯一的开关状态,则其零序电压也是唯一的,此时V2零序电压绝对值最小的开关状态U2'(a2',b2',c2')=U2(a2,b2,c2),否则:
给定变量i=1,如果:
Figure GDA0003419104190000081
基本矢量V2零序电压绝对值最小的开关状态U2'(a2',b2',c2')=U2(a2,b2,c2),否则:
Figure GDA0003419104190000082
基本矢量V2零序电压绝对值最小的开关状态
U2'(a2',b2',c2')=U2'(a2+1,b2+1,c2+1);
或者i∈[2,3,…,n-1],使得:
Figure GDA0003419104190000083
式中,a2+i,b2+i,c2+i,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
a2+i-1,b2+i-1,c2+i-1,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
a2+i+1,b2+i+1,c2+i+1,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
基本矢量V2零序电压绝对值最小的开关状态
U2'(a2',b2',c2')=U2'(a2+i,b2+i,c2+i);
③如果U3(a3,b3,c3)是基本矢量V3唯一的开关状态,则其零序电压也是唯一的,此时V3零序电压绝对值最小的开关状态U3'(a3',b3',c3')=U3(a3,b3,c3),否则:
给定变量i=1,如果:
Figure GDA0003419104190000084
基本矢量V3零序电压绝对值最小的开关状态U3'(a3',b3',c3')=U3(a3,b3,c3),否则:
Figure GDA0003419104190000091
基本矢量V3零序电压绝对值最小的开关状态U3'(a3',b3',c3')=U3'(a3+1,b3+1,c3+1);
或者i∈[2,3,…,n-1],使得:
Figure GDA0003419104190000092
式中,a3+i,b3+i,c3+i,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
a3+i-1,b3+i-1,c3+i-1,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
a3+i+1,b3+i+1,c3+i+1,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
基本矢量V3零序电压绝对值最小的开关状态
U3'(a3',b3',c3')=U3'(a3+i,b3+i,c3+i);
步骤六:利用伏秒平衡原理计算基本矢量V1、V2、V3的作用时间t1、t2、t3
①根据公式(3)计算基本矢量V111)、V222)、V333)的坐标分量:
Figure GDA0003419104190000093
Figure GDA0003419104190000094
Figure GDA0003419104190000095
②根据伏秒平衡原理计算基本矢量V1、V2、V3的作用时间t1、t2、t3
Figure GDA0003419104190000096
式中,TS表示采样周期。
步骤七:利用五段算法,确定开关状态切换路径为
U1'→U2'→U3'→U2'→U1',开关状态的作用时间分别为
t1/2→t2/2→t3→t2/2→t1/2。
步骤八:跟踪参考矢量的采样,重复步骤二到步骤七,完成多电平空间矢量调制。
与现有技术相比,本发明的上述方案降低了参考矢量定位的难度,减少了参考矢量定位的计算量;通过S-K算法快速确定了合成参考矢量的最近三个基本矢量对应的开关状态,通过简单的算法找到了零序电压绝对值最小的开关状态,极大简化了SVM算法的工作量,提高了SVM算法的速度。
附图说明
图1:三相n级多电平变换器拓扑结构图;
图2:基于S-K快速确定基本矢量的级联多电平变换器SVM方法的流程图;
图3:αβ坐标系上二级五电平变换器参考矢量轨迹和基本矢量的分布图;
图4:二级五电平变换器参考矢量Vr顺时针旋转(N-1)π/3角度至第Ⅰ扇区得矢量Vr'和第Ⅰ扇区大三角形三条边的垂线OP、OQ、OL;
图5:二级五电平变换器第Ⅰ大扇区的小三角形S编号图;
图6:第Ⅰ扇区中合成参考矢量Vr的三个基本矢量对应开关状态U1'、U2'、 U3'切换时序图;
图7:第Ⅱ扇区中合成参考矢量Vr的三个基本矢量对应开关状态U1'、U2'、 U3'切换时序图。
具体实施方式
以下将结合附图和具体实施例对本发明做进一步详细说明:
实施例1:传统笛卡尔坐标系αβ上,参考矢量轨迹为一个圆周,圆周的半径随级联单元数n和调制系数m不同而改变。以二级五电平变换器为例,即n=2,设参考电压信号调制系数m=0.85,参考电压信号频率为50Hz,采样周期 Ts=0.2ms,参考矢量轨迹和基本矢量分布如图3所示,图中,小三角形的顶点即为基本矢量,圆周为参考矢量轨迹。t=1.2ms,对参考电压信号进行采样,根据公式(4)计算得v=1.825,v=0.7223,即参考矢量为Vr(1.825,0.7223),如图5所示。参考矢量Vr(1.825,0.7223)的相位角θ=ωt=2πft=21.6°。根据公式(6)判断参考矢量Vr所在大扇区编号N=1,Vr和Vr'是同一个矢量。
根据公式(8)、(9)、(10)、(11)确定Vr'所在小三角形的编号S,p=3、q=4、 l=2,所以S=-p+q2-q+l+1=-3+42-4+2+1=12,即Vr'位于第Ⅰ扇区第12号小三角形内。
通过S-K算法确定合成参考矢量Vr最近的三个基本矢量V1、V2、V3:由于参考矢量Vr所在扇区编号N=1且Vr'所在小三角形编号S=12,S=12不满足
Figure GDA0003419104190000111
Figure GDA0003419104190000112
即得K=3;根据公式(12)求出a1=3+1-2=2、
Figure GDA0003419104190000113
c1=-2,即可确定包含参考矢量Vr的小三角形其中一个顶点矢量V1对应的开关状态U1(2,0,-2);由于S=12、K+1=4同为偶数,根据公式(18)求出a2=a1=2、b2=b1-1=-1、c2=-2,即可确定包含参考矢量Vr的小三角形第二个顶点矢量V2对应的开关状态U2(2,-1,-2),最后根据公式(30) 求出a3=a1-1=1、b3=b1-1=-1、c3=-2,确定包含参考矢量Vr的小三角形第三个顶点矢量V3对应的开关状态U3(1,-1,-2)。
计算基本矢量V1、V2、V3零序电压绝对值最小的开关状态U1'、U2'、U3';在五电平变换器中,a,b,c∈[±2,±1,0],U1(2,0,-2)是基本矢量V1唯一的开关状态,即U1'(2,0,-2)=U1(2,0,-2);U2(2,-1,-2)是唯一的开关状态,即 U2'(2,-1,-2)=U2(2,-1,-2);V3对应两个开关状态(1,-1,-2)、(2,0,-1),根据公式(40) 确定基本矢量V3零序电压绝对值最小的开关状态为U3'(2,0,-1)。
根据公式(3)计算开关状态U1(2,0,-2)、U2(2,-1,-2)、U3(1,-1,-2)对应的基本矢量V111)、V222)、V333)的坐标分量:
Figure GDA0003419104190000114
Figure GDA0003419104190000115
根据公式(45)计算基本矢量V1、V2、V3的作用时间t1、t2、t3
Figure GDA0003419104190000121
计算得t1=0.05ms、t2=0.02ms、t3=0.13ms。
采用五段算法,确定开关状态切换路径为U1'→U2'→U3'→U2'→U1',开关状态的作用时间分别为t1/2→t2/2→t3→t2/2→t1/2,如图6所示。
为了说明参考矢量Vr不在第Ⅰ扇区的情况,以t=7.8ms为例,对参考电压信号进行采样,根据公式(4)计算得v=-1.511、v=1.253,即参考矢量为 Vr(-1.511,1.253),如图5所示。参考矢量Vr(-1.511,1.253)的相位角θ=ωt=2πft=140.4°。根据公式(6)判断参考矢量Vr所在大扇区编号N=3,将Vr顺时针旋转2π/3角度至Ⅰ扇区,得矢量Vr',根据公式(7)计算Vr'的相位角θ1=20.4°。
根据公式(8)、(9)、(10)、(11)确定Vr'所在小三角形的编号S,p=3、q=4、 l=2,所以S=-p+q2-q+l+1=-3+42-4+2+1=12,即Vr'位于第Ⅰ扇区第12号小三角形内。
通过S-K算法确定合成参考矢量Vr最近的三个基本矢量V1、V2、V3:由于参考矢量Vr所在扇区编号N=1且Vr'所在小三角形编号S=12,S=12不满足
Figure GDA0003419104190000122
Figure GDA0003419104190000123
即得K=3;根据公式(14)求出a1=-2、 b1=3+1-2=2、
Figure GDA0003419104190000124
即可确定包含参考矢量Vr的小三角形其中一个顶点矢量V1对应的开关状态U1(-2,2,0);由于S=12、K+1=4同为偶数,根据公式(22)求出a2=a1=-2、b2=b1=2、c2=c1=-1,即可确定包含参考矢量Vr的小三角形第二个顶点矢量V2对应的开关状态U2(-2,2,-1);最后根据公式(31) 求出a3=a1=-2、b3=b1-1=1、c3=c1-1=-1,确定包含参考矢量Vr的小三角形第三个顶点矢量V3对应的开关状态U3(-2,1,-1)。
计算基本矢量V1、V2、V3零序电压绝对值最小的开关状态U1'、U2'、U3';在五电平变换器中,a,b,c∈[±2,±1,0],U1(-2,2,0)是基本矢量V1唯一的开关状态,即U1'(-2,2,0)=U1(-2,2,0);U2(-2,2,-1)是基本矢量V2唯一的开关状态,即 U2'(-2,2,-1)=U2(-2,2,-1);V3对应两个开关状态(-2,1,-1)、(-1,2,0),根据公式(40) 确定基本矢量V3零序电压绝对值最小的开关状态为U3'(-1,2,0)。
根据公式(3)计算开关状态U1(-2,2,0)、U2(-2,2,-1)、U3(-2,1,-1)对应的基本矢量V111)、V222)、V333)的坐标分量:
Figure GDA0003419104190000131
Figure GDA0003419104190000132
根据公式(45)计算基本矢量V1、V2、V3的作用时间t1、t2、t3
Figure GDA0003419104190000133
计算得t1=0.04ms、t2=0.03ms、t3=0.13ms。
采用五段算法,确定开关状态切换路径为U2'→U1'→U3'→U1'→U2',开关状态的作用时间分别为t2/2→t1/2→t3→t1/2→t2/2,如图7所示。

Claims (10)

1.一种基于S-K快速确定基本矢量的多电平变换器SVM方法,其特征在于,针对级联多电平H桥变换器,将三相参考电压信号ura、urb、urc通过3/2转换映射到αβ直角坐标系中得参考矢量轨迹,将开关状态通过3/2转换映射到αβ直角坐标系中得到离散的基本矢量;
步骤1:对三相参考电压信号进行采样,根据采样值计算参考矢量Vr的坐标分量v、v以及相位角θ,确定参考矢量所在大扇区的编号N;
步骤2:如果N=1,矢量Vr'跟参考矢量Vr重合,矢量Vr'的相位角θ1即为参考矢量Vr的相位角θ;如果N≥2,将参考矢量Vr旋转(N-1)π/3角度至第Ⅰ大扇区,得矢量Vr',计算矢量Vr'的相位角θ1;分别计算矢量Vr'在第Ⅰ大扇区三角形三条边的垂线上的投影并取整,根据所得结果计算矢量Vr'所在小三角形的编号S;
步骤3:根据矢量Vr'所在小三角形的编号S,计算得到K值,通过S-K算法确定合成参考矢量Vr最近的三个基本矢量V1、V2、V3对应的开关状态U1、U2、U3
步骤4:根据U1、U2、U3计算基本矢量V1、V2、V3零序电压绝对值最小的开关状态U1'、U2'、U3';利用伏秒平衡原理计算基本矢量V1、V2、V3的作用时间t1、t2、t3;采用五段算法确定开关状态切换路径。
2.如权利要求1所述的基于S-K快速确定基本矢量的多电平变换器SVM方法,其特征在于,按照以下公式确定参考矢量Vr的坐标分量v、v及其相位角θ:
Figure FDA0003419104180000011
θ=ωt (2)。
3.如权利要求1所述的基于S-K快速确定基本矢量的多电平变换器SVM方法,其特征在于,采用以下公式确定参考矢量Vr所在大扇区的编号N:
Figure FDA0003419104180000012
4.如权利要求1所述的基于S-K快速确定基本矢量的多电平变换器SVM方法,其特征在于,如果参考矢量Vr所在大扇区编号N=1,Vr'与Vr为同一个矢量,θ1=θ;如果参考矢量Vr所在大扇区编号N≥2,将参考矢量Vr顺时针旋转(N-1)π/3角度至第Ⅰ大扇区,得矢量Vr',矢量Vr'的相位角θ1
Figure FDA0003419104180000021
5.如权利要求1所述的基于S-K快速确定基本矢量的多电平变换器SVM方法,其特征在于,第Ⅰ大扇区三角形三条边的垂线为OP、OQ、OL,计算参考矢量Vr'在OP、OQ、OL上的投影并取整得:
Figure FDA0003419104180000022
Figure FDA0003419104180000023
l=floor(2mnsinθ1)+1 (7)
S=-p+q2-q+l+1 (8)
式中:p、q、l分别为矢量Vr'在垂线OP、OQ、OL上的投影并取整所得结果;S为矢量Vr'所在小三角形的编号。
6.如权利要求1所述的基于S-K快速确定基本矢量的多电平变换器SVM方法,其特征在于,通过S-K算法确定合成参考矢量Vr最近的三个基本矢量V1、V2、V3对应的开关状态U1(a1,b1,c1)、U2(a2,b2,c2)、U3(a3,b3,c3):
①如果
Figure FDA0003419104180000024
Figure FDA0003419104180000025
否则,令
Figure FDA0003419104180000026
②确定包含参考矢量Vr的小三角形其中一个顶点矢量V1,该矢量对应的开关状态U1(a1,b1,c1):当参考矢量位于第Ⅰ大扇区时,即N=1,V1对应的开关状态U1(a1,b1,c1)为:
Figure FDA0003419104180000027
当参考矢量位于第Ⅱ大扇区时,即N=2,V1对应的开关状态U1(a1,b1,c1)为:
Figure FDA0003419104180000031
当参考矢量位于第Ⅲ大扇区时,即N=3,V1对应的开关状态U1(a1,b1,c1)为:
Figure FDA0003419104180000032
当参考矢量为第Ⅳ大扇区时,即N=4,V1对应的开关状态U1(a1,b1,c1)为:
Figure FDA0003419104180000033
当参考矢量位于第Ⅴ大扇区时,即N=5时,V1对应的开关状态U1(a1,b1,c1)为:
Figure FDA0003419104180000034
当参考矢量位于第Ⅵ大扇区时,即N=6时,V1对应的开关状态U1(a1,b1,c1)为:
Figure FDA0003419104180000035
③根据参考矢量所在小三角形三个顶点矢量之间的关系,判断矢量Vr'所在小三角形的编号S和K+1的奇偶性是否相同,确定基本矢量V2对应的开关状态U2(a2,b2,c2):
N=1时,S和K+1的奇偶性相同:
Figure FDA0003419104180000041
S和K+1的奇偶性不相同:
Figure FDA0003419104180000042
N=2时,S和K+1的奇偶性相同:
Figure FDA0003419104180000043
S和K+1的奇偶性不相同:
Figure FDA0003419104180000044
N=3时,S和K+1的奇偶性相同:
Figure FDA0003419104180000045
S和K+1的奇偶性不相同:
Figure FDA0003419104180000046
N=4时,S和K+1的奇偶性相同:
Figure FDA0003419104180000047
S和K+1的奇偶性不相同:
Figure FDA0003419104180000048
N=5时,S和K+1的奇偶性相同:
Figure FDA0003419104180000051
S和K+1的奇偶性不相同:
Figure FDA0003419104180000052
N=6时,S和K+1的奇偶性相同:
Figure FDA0003419104180000053
S和K+1的奇偶性不相同:
Figure FDA0003419104180000054
④V3对应的开关状态U3(a3,b3,c3):
当N=1或N=2时,V3对应的开关状态U3(a3,b3,c3):
Figure FDA0003419104180000055
当N=3或N=4时,V3对应的开关状态U3(a3,b3,c3):
Figure FDA0003419104180000056
当N=5或N=6时,V3对应的开关状态U3(a3,b3,c3):
Figure FDA0003419104180000057
7.如权利要求1所述的基于S-K快速确定基本矢量的多电平变换器SVM方法,其特征在于,计算基本矢量V1、V2、V3的开关状态U1'(a1',b1',c1')、U2'(a2',b2',c2')、U3'(a3',b3',c3'),使对应的零序电压绝对值最小:
①如果U1(a1,b1,c1)是基本矢量V1唯一的开关状态,则其零序电压也是唯一的,此时V1零序电压绝对值最小的开关状态U1'(a1',b1',c1')=U1(a1,b1,c1),否则:
给定变量i=1,如果:
Figure FDA0003419104180000061
基本矢量V1零序电压绝对值最小的开关状态U1'(a1',b1',c1')=U1(a1,b1,c1),否则:
Figure FDA0003419104180000062
基本矢量V1零序电压绝对值最小的开关状态U1'(a1',b1',c1')=U1'(a1+1,b1+1,c1+1);
或者i∈[2,3,…,n-1],使得:
Figure FDA0003419104180000063
式中,a1+i,b1+i,c1+i,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
a1+i-1,b1+i-1,c1+i-1,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
a1+i+1,b1+i+1,c1+i+1,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
基本矢量V1零序电压绝对值最小的开关状态U1'(a1',b1',c1')=U1'(a1+i,b1+i,c1+i);
②如果U2(a2,b2,c2)是基本矢量V2唯一的开关状态,则其零序电压也是唯一的,此时V2零序电压绝对值最小的开关状态U2'(a2',b2',c2')=U2(a2,b2,c2),否则:
给定变量i=1,如果:
Figure FDA0003419104180000064
基本矢量V2零序电压绝对值最小的开关状态U2'(a2',b2',c2')=U2(a2,b2,c2),否则:
Figure FDA0003419104180000065
基本矢量V2零序电压绝对值最小的开关状态U2'(a2',b2',c2')=U2'(a2+1,b2+1,c2+1);
或者i∈[2,3,…,n-1],使得:
Figure FDA0003419104180000071
式中,a2+i,b2+i,c2+i,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
a2+i-1,b2+i-1,c2+i-1,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
a2+i+1,b2+i+1,c2+i+1,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
基本矢量V2零序电压绝对值最小的开关状态U2'(a2',b2',c2')=U2'(a2+i,b2+i,c2+i);
③如果U3(a3,b3,c3)是基本矢量V3唯一的开关状态,则其零序电压也是唯一的,此时V3零序电压绝对值最小的开关状态U3'(a3',b3',c3')=U3(a3,b3,c3),否则:
给定变量i=1,如果:
Figure FDA0003419104180000072
基本矢量V3零序电压绝对值最小的开关状态U3'(a3',b3',c3')=U3(a3,b3,c3),否则:
Figure FDA0003419104180000073
基本矢量V3零序电压绝对值最小的开关状态U3'(a3',b3',c3')=U3'(a3+1,b3+1,c3+1);
或者i∈[2,3,…,n-1],使得:
Figure FDA0003419104180000074
式中,a3+i,b3+i,c3+i,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
a3+i-1,b3+i-1,c3+i-1,∈[±n,±(n-1),±(n-2),...,±2,±1,0],
a3+i+1,b3+i+1,c3+i+1,∈[±n,±(n-1),±(n-2),...±2,±1,0],
基本矢量V3零序电压绝对值最小的开关状态U3'(a3',b3',c3')=U3'(a3+i,b3+i,c3+i)。
8.如权利要求1所述的基于S-K快速确定基本矢量的多电平变换器SVM方法,其特征在于,利用伏秒平衡原理计算基本矢量V1、V2、V3的作用时间t1、t2、t3
①根据公式3,利用开关状态U1(a1,b1,c1)、U2(a2,b2,c2)、U3(a3,b3,c3)计算基本矢量V111)、V222)、V333)的坐标分量:
Figure FDA0003419104180000081
Figure FDA0003419104180000082
Figure FDA0003419104180000083
②根据伏秒平衡原理计算基本矢量V1、V2、V3合成参考矢量Vr的作用时间t1、t2、t3
Figure FDA0003419104180000084
式中,TS表示采样周期。
9.如权利要求1所述的基于S-K快速确定基本矢量的多电平变换器SVM方法,其特征在于,采用五段算法,确定开关状态切换路径为U1'→U2'→U3'→U2'→U1',开关状态的作用时间分别为t1/2→t2/2→t3→t2/2→t1/2。
10.如权利要求1所述的基于S-K快速确定基本矢量的多电平变换器SVM方法,其特征在于,进入下一个周期后,跟踪参考矢量的采样,重复步骤1~步骤4,完成多电平变换器空间矢量调制。
CN202110133184.9A 2021-01-29 2021-01-29 基于s-k快速确定基本矢量的多电平变换器svm方法 Active CN112803809B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110133184.9A CN112803809B (zh) 2021-01-29 2021-01-29 基于s-k快速确定基本矢量的多电平变换器svm方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110133184.9A CN112803809B (zh) 2021-01-29 2021-01-29 基于s-k快速确定基本矢量的多电平变换器svm方法

Publications (2)

Publication Number Publication Date
CN112803809A CN112803809A (zh) 2021-05-14
CN112803809B true CN112803809B (zh) 2022-03-29

Family

ID=75813209

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110133184.9A Active CN112803809B (zh) 2021-01-29 2021-01-29 基于s-k快速确定基本矢量的多电平变换器svm方法

Country Status (1)

Country Link
CN (1) CN112803809B (zh)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103684013B (zh) * 2013-12-17 2017-04-05 中国矿业大学 一种多电平逆变器svpwm控制方法
CN107302318B (zh) * 2017-08-02 2019-06-07 南昌工程学院 一种基于h桥级联多电平逆变器的移相空间矢量调制方法
CN108418461B (zh) * 2018-04-04 2019-10-01 南昌工程学院 一种三角形连接级联h桥逆变器的空间矢量调制方法
CN110071655B (zh) * 2019-05-21 2020-06-09 南昌工程学院 一种简化的多电平变换器空间矢量调制方法

Also Published As

Publication number Publication date
CN112803809A (zh) 2021-05-14

Similar Documents

Publication Publication Date Title
CN110867887B (zh) Npc三电平逆变器最优序列模型预测控制装置及方法
US6151227A (en) Controller for power converter
CN104796024A (zh) 一种适用于三相级联h桥型变换器的简化多电平空间矢量脉宽调制方法及其调制软核
Li et al. Capacitor voltage balance control of five-level modular composited converter with hybrid space vector modulation
CN105186898A (zh) 适用于任意电平数单相级联h桥型变换器的简化多电平空间矢量脉宽调制方法及其调制软核
CN112803809B (zh) 基于s-k快速确定基本矢量的多电平变换器svm方法
CN105099221A (zh) 任意电平数单相级联三电平桥型变换器的简化多电平空间矢量脉宽调制方法及其调制软核
CN109347350B (zh) 一种三相多电平变流器及其电池soc均衡控制方法
CN112865572B (zh) 一种基于s-k及svpwm的级联多电平变换器
CN109039124B (zh) 基于移相空间矢量调制的mmc电容电压均衡控制方法
Sunbul et al. Simplified SVPWM method for the Vienna rectifier
CN112701947B (zh) 基于零序电压最小空间矢量调制控制的级联多电平变换器
Kumar et al. An effective Space-Vector PWM Method for Multi-level inverter based on two-level inverter
CN109004814B (zh) 一种用于mmc的子模块电容电压均衡控制系统
Botao et al. New SVPWM control scheme for three-phase diode clamping multilevel inverter with balanced dc voltages
Wang et al. Fast FCS-MPC-based SVPWM method to reduce switching states of multilevel cascaded H-bridge STATCOMs
CN114938151A (zh) 一种并网逆变器的控制方法
CN114884383A (zh) 一种简化型三电平逆变器的广义虚拟矢量调制方法及系统
CN110071653B (zh) 五电平npc变换器低调制度直流侧电容电压自平衡方法
CN112910298B (zh) 一种级联多电平变换器零序电压最小空间矢量调制方法
Massoud et al. Systematic analytical-based generalised algorithm for multilevel space vector modulation with a fixed execution time
Ji et al. A simplified virtual space vector modulation strategy for four-level t-type nnpc converter
CN113489406A (zh) 一种应用于九相开端绕组电机的svpwm谐波抑制方法
Periyaazhagar et al. Asymmetrical cascaded multi level inverter using control freedom pulse width modulation techniques
CN110601578B (zh) 一种最近电平等效的空间矢量调制方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant