CN106909711A - 一种求分数阶ccm开关变换器瞬态解的方法 - Google Patents
一种求分数阶ccm开关变换器瞬态解的方法 Download PDFInfo
- Publication number
- CN106909711A CN106909711A CN201710019197.7A CN201710019197A CN106909711A CN 106909711 A CN106909711 A CN 106909711A CN 201710019197 A CN201710019197 A CN 201710019197A CN 106909711 A CN106909711 A CN 106909711A
- Authority
- CN
- China
- Prior art keywords
- fractional order
- omega
- component
- transient
- switch converters
- 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.)
- Granted
Links
- 230000001052 transient effect Effects 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 63
- 230000010355 oscillation Effects 0.000 claims abstract description 70
- 230000008569 process Effects 0.000 claims abstract description 10
- 230000006870 function Effects 0.000 claims description 32
- 230000001939 inductive effect Effects 0.000 claims description 23
- 238000012937 correction Methods 0.000 claims description 20
- 238000013178 mathematical model Methods 0.000 claims description 16
- 238000004458 analytical method Methods 0.000 claims description 12
- 239000000203 mixture Substances 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000011438 discrete method Methods 0.000 claims description 6
- 230000000737 periodic effect Effects 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 230000005611 electricity Effects 0.000 claims description 2
- 238000004088 simulation Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000012795 verification Methods 0.000 description 2
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Dc-Dc Converters (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种求分数阶CCM开关变换器瞬态解的方法,利用等效小参量法,将复杂的分数阶微分运算转换为微分算子的指数运算,从而将复杂的分数阶变换器状态变量稳态周期解的求解转化为直流分量和各次谐波分量幅度的求解,而直流分量和谐波分量幅度的求解则能够利用谐波平衡法;只需要分析分数阶开关变换器状态变量主振荡方程的瞬态解,再结合其稳态解即可得到分数阶开关变换器状态变量的瞬态解析解。本发明方法能够较好的分析分数阶CCM开关变换器的瞬态过程。
Description
技术领域
本发明涉及分数阶开关变换器的建模与分析领域,尤其是指一种求分数阶CCM开关变换器瞬态解的方法。
背景技术
过去针对开关变换器常用的建模与分析方法有:离散迭代映射模型、基于状态空间平均法的模型、基于电路原理(KCL、KVL)的分段线性模型[1]以及结合谐波平衡与扰动法的等效小参量法[2],这些方法的分析对象均是整数阶的开关变换器,即变换器中的电容、电感都是整数阶的元件,然而文献[3][4]指出实际上电容和电感在本质上均是分数阶的,这就需要为变换器建立相应的分数阶模型。
谐波平衡法是分析非线性系统周期解的常用方法之一,但是它只适用求系统的稳态解,文献[5]提出的通用平均法将谐波平衡法推广到瞬态分析,但是它不易求得高次谐波。文献[6]用结合谐波平衡与扰动法的等效小参量法分析了开关变换器的瞬态过程,但并未考虑电容和电感的分数阶特性。文献[7]考虑了电感电容的分数阶特性,用等效小参量法分析了分数阶开关变换器的稳态特性,但未分析分数阶开关变换器的瞬态过程。
参考文献
[1]罗晓曙.DC-DC变换器的非线性动力学行为与混沌控制[M].科学出版社,2012,pp.23-24.
[2]Qiu S S,FilanovskI M.Calculation of steady-state oscillations innon linear circuits[J],Int.J.Electronics,1989,67(3),pp.403-414.
[3]Westerlund S.Dead Matter Has Memory![M].Kalmar,Sweden:CausalConsulting,2002,Chap.7.
[4]Podlubny I.Fractional Differential Equations[M].San Diego:AcademicPress,1999,Chap.2.
[5]Sanders SR,Noworolski J,Liu X.et al.Generalized averaging methodfor power conversion circuits.IEEE Trans.on Power Electronics,1991,PE-6(2);251-259.
[6]林波涛,邱水生.适于瞬态分析的等效小参量法及其在PWM开关变换器中的应用[J].电子科学学刊,1997,19(5),pp.649-657.
[7]Xi Chen,Yanfeng Chen,Bo Zhang,Qiuyuan Dong.A Modeling and AnalysisMethod for Fractional-oder DC-DC Converters.IEEE Trans.on Power Electronics,2016,Vol.pp,no.99,pp.1-1.
发明内容
本发明的目的是在于克服现有人机交互方式的不足,提出了一种求分数阶CCM开关变换器瞬态解的方法,能够较好的分析分数阶CCM开关变换器的瞬态过程。
为实现上述目的,本发明所提供的技术方案为:一种求分数阶CCM开关变换器瞬态解的方法,利用等效小参量法,将复杂的分数阶微分运算转换为微分算子的指数运算,从而将复杂的分数阶变换器状态变量稳态周期解的求解转化为直流分量和各次谐波分量幅度的求解,而直流分量和谐波分量幅度的求解则能够利用谐波平衡法;只需要分析分数阶开关变换器状态变量主振荡方程的瞬态解,再结合其稳态解即可得到分数阶开关变换器状态变量的瞬态解析解;其包括以下步骤:
S1、建立以分数阶微分方程描述的分数阶开关变换器的非线性数学模型;
S2、利用等效小参量法得到分数阶开关变换器的等效数学模型
利用等效小参量法求解S1中的非线性数学模型,得到描述分数阶开关变换器的等效数学方程组,即分数阶开关变换器的等效数学模型;该等效数学方程组包含一个求系统状态变量主振荡分量的主振荡分数阶微分方程,和一系列求状态变量修正量的分数阶微分方程;
S3、利用谐波平衡法求分数阶开关变换器系统状态变量的稳态周期解
利用谐波平衡法逐步求解S2中等效数学方程组中的各个分数阶微分方程的稳态解,得到分数阶开关变换器状态变量稳态周期解的近似解表达式,所得到的近似周期解包含直流分量和纹波分量,其中纹波分量由基波和各次谐波组成;
S4、求分数阶开关变换器主振荡分量瞬态解
根据分数阶微积分定义求S2中主振荡分数阶微分方程瞬态解,得到分数阶开关变换器状态变量的主振荡分量的瞬态解;
S5、获取分数阶开关变换器系统状态变量瞬态解的解析表达式
把S3所得到纹波分量与S4所得到的主振荡分量的瞬态解叠加,即得到分数阶开关变换器的瞬态解的解析表达式。
在步骤S1中,所建立的以分数阶微分方程描述的分数阶开关变换器的非线性数学模型如下:
G0(pα,pβ)x+G1f(x)=U (1)
式中,x=[iL vC]T表示开关变换器系统的状态变量,其中iL为电感电流瞬时值,vC为电容电压瞬时值,上标“T”表示求矩阵的转置;p代表微分算子,其定义为p=d/dt,相应的分数阶微分运算dα/dtα和dβ/dtβ分别记为pα和pβ,其中α和β分别表示分数阶电感和分数阶电容的阶次;系数矩阵G0(pα,pβ)、G1由具体的开关变换器的电路参数所确定,其中G0(pα,pβ)同分数阶微分运算相关;f(x)=δx是一个非线性函矢量函数,其中δ=1为一个表征开关变换器中受控开关通断状态的开关函数,当该受控开关导通时δ=1,当该受控开关断开时δ=0;U为变换器输入电压向量;
在步骤S2中,所述的等效小参量法的具体步骤如下:
S21、将开关函数δ表示成主振荡分量和修正量之和的级数形式:其中δ0表示开关函数的主振荡分量,δi表示开关函数的第i阶修正量,它们能够根据具体开关函数的傅里叶级数来确定;
S22、将待求解的状态变量x也表示成级数形式:其中x0表示状态变量的主振荡分量,xi表示状态变量的第i阶修正量,它们在具体的求解过程逐步确定;
S23、将δ和x的级数表达式代入非线性矢量函数f(x)=δx中,得到其中f0表示非线性矢量函数的主振荡分量,fi表示非线性矢量函数的第i阶修正量;
S24、将f0表示为f0=f0m+εR1,将fi表示为fi=fim+εRi+1,其中f0m为f0的主项,包含f0中所有与x0具有相同频率成分的项,R1为f0的余项,包含f0中所有与x0具有不同频率成分的项;同理,fim为fi的主项,包含fi中所有与xi具有相同频率成分的项,Ri+1为fi的余项,包含fi中所有与xi具有不同频率成分的项;
S25、将f0=f0m+εR1和fi=fim+εRi+1代入中,从而能够将f(x)表示为
在上述步骤S21~S25中,上标或下标i是一个整数,i=1,2,……;ε是引入的一个小量标记,εixi表明xi是状态变量x的第i阶小量,并且有εi+1xi+1<εixi<x0,当在运算过程中需要具体数值时ε=1;
S26、将和代入到公式(1)中,并令等式两边具有相同εi项的项分别相等,能够得到描述分数阶CCM开关变换器的等效数学模型,如下式(2):
式(2)中的第1个分数阶微分方程用于求状态变量的主振荡分量x0,称为主振荡方程;第2~n个分数阶微分方程用于求状态变量的各阶修正量xi(i=1,2,……n),称为修正量方程;
在步骤S3中,以指数函数表示状态变量的稳态周期解xS的近似数学表达式如下:
式中,直流分量XDC=A0为分数开关变换器的状态变量稳态主振荡分量;xac为纹波分量,其中A1为基波的的幅度向量,为其共轭;Ai(i=2,3,……,n)为第i次波的的幅度向量,为其共轭;ω为开关变换器的角频率,t表示时间变量,j为虚数单位;公式(3)中的状态变量的稳态周期解的数学表达式也能够用三角函数的形式表示如下式(4):
xS=XDC+xac=A0+2Re(A1)cosωt-2Im(A1)sinωt+2Re(A2)cos2ωt-2Im(A2)sin2ωt
+2Re(A3)cos3ωt-2Im(A3)sin3ωt+...+2Re(An)cos nωt-2Im(An)sin nωt
(4)
公式(4)中的Re(Ai)(i=2,3,……,n)分别表示复数向量Ai的实部,Im(Ai)分别表示复数向量Ai的虚部;
在步骤S4中,根据Grnwald-Letnikov定义用离散方法求分数阶主振荡微分方程G0(pα,pβ)x0+G1f0m=U的瞬态数值解,再通过Matlab拟合得到瞬态解的近似表达式;
将主振荡微分方程写成如下关于系统状态变量的分数阶微分方程组的形式:
其中,a1,a2,a3,a4,b1,b2,b3为跟具体电路参数有关的常系数,α和β分别表示分数阶电感和分数阶电容的阶次;根据Grnwald-Letnikov定义用离散方法得到分数阶电容电压和电感电流主振荡分量瞬态近似数值解分别为式(6a)和(6b):
式中vC0-trans和iL0-trans分别表示分数阶电容电压和电感电流主振荡分量的瞬态解,h为步长,vt-kh和it-kh表示t-kh时刻电容电压和电感电流的瞬时值,式中能够由递推公式得出,k为自然数,k=1,2,....;将式(6a)和(6b)所得到的数值解用Matlab多项式拟合能够得到主振荡分量瞬态解的近似解析表达式x0-trans=[iL0-trans vC0-trans]T,其中vC0-trans和iL0-trans的近似解析表达式分别为式(7a)、(7b):
式中V0和I0分别表示分数阶电容电压和电感电流主振荡分量瞬态解的直流分量,ωv和ωi分别表示分数阶电容电压和电感电流主振荡分量瞬态解的基波角频率,avk和bvk分别表示分数阶电容电压主振荡分量瞬态解第k次谐波中余弦分量和正弦分量的幅值,aik和bik分别表示分数阶电感电流主振荡分量瞬态解第k次谐波中余弦分量和正弦分量的幅值,自然数k=1,2,3,……n;
在步骤S5中,把S3所得到纹波分量xac与S4所得到的主振荡分量瞬态解x0-trans叠加,即可得到分数阶开关变换器的瞬态解析解的近似表达式xtrans≈x0-trans+xac。
本发明与现有技术相比,具有如下优点与有益效果:
1、由本发明所提方法的求解公式可以看出,采用本方法求分数阶开关变换器状态变量的瞬态解,相当于将求解非整数阶微积分运算的复杂过程转化为矩阵运算和求线性方程(组)的过程,计算量小,结果直观,清晰。
2、通过与PSIM电路仿真波形的对比验证,该发明所提供的方法能够较好的分析分数阶CCM开关变换器的瞬态过程,便于电路的优化控制与设。
附图说明
图1为CCM分数阶Buck变换器的电路原理图。
图2a为本发明方法得到的电感电流波形与PSIM电路仿真波形的对比图。
图2b为本发明方法得到的电容电压波形与PSIM电路仿真波形的对比图。
具体实施方式
下面结合具体实施例对本发明作进一步说明。
本实施例所提供的求分数阶CCM开关变换器瞬态解的方法,具体包括以下步骤:
S1、建立以分数阶微分方程描述的分数阶开关变换器的非线性数学模型:
G0(pα,pβ)x+G1f(x)=U (1)
上式中x=[iL vC]T表示开关变换器系统的状态变量,其中iL为电感电流瞬时值,vC为电容电压瞬时值,上标“T”表示求矩阵的转置;p代表微分算子,其定义为p=d/dt,相应的分数阶微分运算dα/dtα和dβ/dtβ可分别记为pα和pβ,其中α和β分别表示分数阶电感和分数阶电容的阶次;系数矩阵G0(pα,pβ)、G1由具体的开关变换器的电路参数所确定,其中G0(pα,pβ)同分数阶微分运算相关;f(x)=δx是一个非线性函矢量函数,其中δ=1为一个表征开关变换器中受控开关通断状态的开关函数,当该受控开关导通时δ=1,当该受控开关断开时δ=0;U为变换器输入电压向量。
S2、利用等效小参量法得到分数阶开关变换器的等效数学模型
具体为:利用等效小参量法求解分数阶CCM开关变换器的非线性数学模型,得到描述该分数阶开关变换器的等效数学方程组,即等效数学模型。该等效数学方程组包含一个求分数阶CCM开关变换器状态变量主振荡分量的分数阶微分方程,和一系列求状态变量各阶修正量的分数阶微分方程。
所述的等效小参量法的具体步骤如下:
S21、将开关函数δ表示成主振荡分量和修正量之和的级数形式:其中δ0表示开关函数的主振荡分量,δi表示开关函数的第i阶修正量,它们可根据具体开关函数的傅里叶级数来确定;
S22、将待求解的状态变量x也表示成如下级数形式:其中x0表示状态变量的主振荡分量,xi表示状态变量的第i阶修正量,它们在具体的求解过程逐步确定;
S23、将δ和x的级数表达式代入非线性矢量函数f(x)=δx中,得到其中f0表示非线性矢量函数的主振荡分量,fi表示非线性矢量函数的第i阶修正量;
S24、将f0表示为f0=f0m+εR1,将fi表示为fi=fim+εRi+1,其中f0m为f0的主项,包含f0中所有与x0具有相同频率成分的项,R1为f0的余项,包含f0中所有与x0具有不同频率成分的项;类似的,fim为fi的主项,包含fi中所有与xi具有相同频率成分的项,Ri+1为fi的余项,包含fi中所有与xi具有不同频率成分的项;
S25、将f0=f0m+εR1和fi=fim+εRi+1代入中,从而可将f(x)表示为
在上述步骤S21~S25中,上标或下标i是一个整数,i=1,2,……;ε是本方法所引入的一个小量标记,例如εixi表明xi是状态变量x的第i阶小量,并且有εi+1xi+1<εixi<x0,当在运算过程中需要具体数值时ε=1;
S26、将和代入到公式(1)中,并令等式两边具有相同εi项的项分别相等,可得到描述分数阶CCM开关变换器的等效数学模型,如下式(2):
上式(2)中所述的第1个分数阶微分方程用于求状态变量的主振荡分量x0,称为主振荡方程;第2~n个分数阶微分方程用于求状态变量的各阶修正量xi(i=1,2,……n),称为修正量方程。
S3、利用谐波平衡法求分数阶开关变换器系统状态变量的稳态周期解
利用谐波平衡法逐步求解公式(2)所述的等效数学方程模型中的各个分数阶微分方程的稳态解,得到分数阶开关变换器状态变量稳态周期解的近似解表达式,所得到的近似周期解包含直流分量和纹波分量,其中纹波分量由基波和各次谐波组成。所得到的以指数函数表示状态变量的稳态周期解xS的近似数学表达式如下:
式中直流分量XDC=A0为分数开关变换器的状态变量稳态主振荡分量;xac为纹波分量,其中A1为基波的的幅度向量,为其共轭;Ai(i=2,3,……,n)为第i次波的的幅度向量,为其共轭;ω为开关变换器的角频率,t表示时间变量,j为虚数单位。公式(3)中的状态变量的稳态周期解的数学表达式也可以用三角函数的形式表示如下式(4):
xS=XDC+xac=A0+2Re(A1)cosωt-2Im(A1)sinωt+2Re(A2)cos2ωt-2Im(A2)sin2ωt
+2Re(A3)cos3ωt-2Im(A3)sin3ωt+...+2Re(An)cos nωt-2Im(An)sin nωt
(4)
公式(4)中的Re(Ai)(i=2,3,……,n)分别表示复数向量Ai的实部,Im(Ai)分别表示复数向量Ai的虚部。
S4、求分数阶开关变换器主振荡分量瞬态解
根据Grnwald-Letnikov定义用离散方法求分数阶主振荡微分方程G0(pα,pβ)x0+G1f0m=U的瞬态数值解,再通过Matlab拟合得到瞬态解的近似表达式。
将主振荡微分方程写成如下关于系统状态变量的分数阶微分方程组的形式:
其中a1,a2,a3,a4,b1,b2,b3为跟具体电路参数有关的常系数,α和β分别表示分数阶电感和分数阶电容的阶次。根据Grnwald-Letnikov定义用离散方法可得到分数阶电容电压和电感电流主振荡分量瞬态近似数值解分别为式(6a)和(6b):
式中,vC0-trans和iL0-trans分别表示分数阶电容电压和电感电流主振荡分量的瞬态解,h为步长,vt-kh和it-kh表示t-kh时刻电容电压和电感电流的瞬时值,式中可以由递推公式得出,k为自然数,k=1,2,....。将式(6a)和(6b)所得到的数值解用Matlab多项式拟合可以得到主振荡分量瞬态解的近似解析表达式x0-trans=[iL0-trans vC0-trans]T,其中vC0-trans和iL0-trans的近似解析表达式分别为式(7a)、(7b):
式中,V0和I0分别表示分数阶电容电压和电感电流主振荡分量瞬态解的直流分量,ωv和ωi分别表示分数阶电容电压和电感电流主振荡分量瞬态解的基波角频率,avk和bvk分别表示分数阶电容电压主振荡分量瞬态解第k次谐波中余弦分量和正弦分量的幅值,aik和bik分别表示分数阶电感电流主振荡分量瞬态解第k次谐波中余弦分量和正弦分量的幅值,自然数k=1,2,3,……n。
S5、获取分数阶开关变换器系统状态变量瞬态解
把S3所得到纹波分量xac与S4所得到的主振荡分量瞬态解x0-trans叠加,即得到分数阶开关变换器的瞬态解析解的近似表达式xtrans≈x0-trans+xac。
下面结合附图1、2a、2b对本实施例上述方法进行具体说明。
图1为CCM分数阶Buck变换器的电路原理图,其中,S表示开关,SD表示二极管,VDC表示直流电源,RL表示负载电阻,Lα与Cβ分别表示分数阶电容和分数阶电感。对于图1中的分数阶电感电流连续模式(CCM)Buck变换器,其电路参数为开关频率fs=2.5kHz(角频率为ω=5π×103rad/s),输入电压VDC=20V,分数阶电感L=20mH,阶次α=0.9,分数阶电容C=47μF,阶次β=0.9,负载电阻RL=22Ω,占空比D=0.6。取状态变量x=[iL vC]T,iL和vC分别为分数阶电感电流和电容电压。
根据上述方法步骤可求得分数阶Buck变换器状态变量稳态周期解纹波分量的近似解析表达式为xac=[iLac vCac]T,其中iLac和vCac分别为分数阶电感电流和电容电压纹波,它们分别由式(8a)和(8b)表示:
iLac≈-0.1028cosωt-0.0148sinωt-0.0076cos2ωt+0.0015sin 2ωt
-0.0056cos3ωt-0.0056sin3ωt (8a)
vCac=-0.0596cosωt-0.3516sinωt-0.0308cos2ωt-0.0072sin2ωt
+0.0054cos3ωt-0.0086sin3ωt (8b)
分数阶CCM-Buck变换器主振荡微分方程可表示为式(9):
根据步骤S4可得到主振荡分量近似瞬态解如式(7a)和(7b)所示,其中V0=-2.03V,ωv=218rad/s,I0=16.16A,ωi=158rad/s,其它各系数如下表1所示。
表1
根据步骤S5可得分数阶CCM-Buck变换器状态变量的瞬态解析解。
将本发明方法分别与PSIM电路仿真得到的瞬态电感电流、电容电压波形进行比较,如图2a、2b仿真结果对比验证图所示,图中实线为本发明得到的波形,虚线为PSIM电路仿真得到的波形。从图中可以看出两条曲线相符合,说明本发明所提出的方法是有效的。
以上所述实施例只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。
Claims (2)
1.一种求分数阶CCM开关变换器瞬态解的方法,其特征在于:利用等效小参量法,将复杂的分数阶微分运算转换为微分算子的指数运算,从而将复杂的分数阶变换器状态变量稳态周期解的求解转化为直流分量和各次谐波分量幅度的求解,而直流分量和谐波分量幅度的求解则能够利用谐波平衡法;只需要分析分数阶开关变换器状态变量主振荡方程的瞬态解,再结合其稳态解即可得到分数阶开关变换器状态变量的瞬态解析解;其包括以下步骤:
S1、建立以分数阶微分方程描述的分数阶开关变换器的非线性数学模型;
S2、利用等效小参量法得到分数阶开关变换器的等效数学模型
利用等效小参量法求解S1中的非线性数学模型,得到描述分数阶开关变换器的等效数学方程组,即分数阶开关变换器的等效数学模型;该等效数学方程组包含一个求系统状态变量主振荡分量的主振荡分数阶微分方程,和一系列求状态变量修正量的分数阶微分方程;
S3、利用谐波平衡法求分数阶开关变换器系统状态变量的稳态周期解
利用谐波平衡法逐步求解S2中等效数学方程组中的各个分数阶微分方程的稳态解,得到分数阶开关变换器状态变量稳态周期解的近似解表达式,所得到的近似周期解包含直流分量和纹波分量,其中纹波分量由基波和各次谐波组成;
S4、求分数阶开关变换器主振荡分量瞬态解
根据分数阶微积分定义求S2中主振荡分数阶微分方程瞬态解,得到分数阶开关变换器状态变量的主振荡分量的瞬态解;
S5、获取分数阶开关变换器系统状态变量瞬态解的解析表达式
把S3所得到纹波分量与S4所得到的主振荡分量的瞬态解叠加,即得到分数阶开关变换器的瞬态解的解析表达式。
2.根据权利要求1所述的一种求分数阶CCM开关变换器瞬态解的方法,其特征在于,在步骤S1中,所建立的以分数阶微分方程描述的分数阶开关变换器的非线性数学模型如下:
G0(pα,pβ)x+G1f(x)=U (1)
式中,x=[iL vC]T表示开关变换器系统的状态变量,其中iL为电感电流瞬时值,vC为电容电压瞬时值,上标“T”表示求矩阵的转置;p代表微分算子,其定义为p=d/dt,相应的分数阶微分运算dα/dtα和dβ/dtβ分别记为pα和pβ,其中α和β分别表示分数阶电感和分数阶电容的阶次;系数矩阵G0(pα,pβ)、G1由具体的开关变换器的电路参数所确定,其中G0(pα,pβ)同分数阶微分运算相关;f(x)=δx是一个非线性函矢量函数,其中δ=1为一个表征开关变换器中受控开关通断状态的开关函数,当该受控开关导通时δ=1,当该受控开关断开时δ=0;U为变换器输入电压向量;
在步骤S2中,所述的等效小参量法的具体步骤如下:
S21、将开关函数δ表示成主振荡分量和修正量之和的级数形式:其中δ0表示开关函数的主振荡分量,δi表示开关函数的第i阶修正量,它们能够根据具体开关函数的傅里叶级数来确定;
S22、将待求解的状态变量x也表示成级数形式:其中x0表示状态变量的主振荡分量,xi表示状态变量的第i阶修正量,它们在具体的求解过程逐步确定;
S23、将δ和x的级数表达式代入非线性矢量函数f(x)=δx中,得到其中f0表示非线性矢量函数的主振荡分量,fi表示非线性矢量函数的第i阶修正量;
S24、将f0表示为f0=f0m+εR1,将fi表示为fi=fim+εRi+1,其中f0m为f0的主项,包含f0中所有与x0具有相同频率成分的项,R1为f0的余项,包含f0中所有与x0具有不同频率成分的项;同理,fim为fi的主项,包含fi中所有与xi具有相同频率成分的项,Ri+1为fi的余项,包含fi中所有与xi具有不同频率成分的项;
S25、将f0=f0m+εR1和fi=fim+εRi+1代入中,从而能够将f(x)表示为
在上述步骤S21~S25中,上标或下标i是一个整数,i=1,2,……;ε是引入的一个小量标记,εixi表明xi是状态变量x的第i阶小量,并且有εi+1xi+1<εixi<x0,当在运算过程中需要具体数值时ε=1;
S26、将和代入到公式(1)中,并令等式两边具有相同εi项的项分别相等,能够得到描述分数阶CCM开关变换器的等效数学模型,如下式(2):
式(2)中的第1个分数阶微分方程用于求状态变量的主振荡分量x0,称为主振荡方程;第2~n个分数阶微分方程用于求状态变量的各阶修正量xi,i=1,2,……n,称为修正量方程;
在步骤S3中,以指数函数表示状态变量的稳态周期解xS的近似数学表达式如下:
式中,直流分量XDC=A0为分数开关变换器的状态变量稳态主振荡分量;xac为纹波分量,其中A1为基波的的幅度向量,为其共轭;Ai为第i次波的的幅度向量,i=2,3,……,n,为其共轭;ω为开关变换器的角频率,t表示时间变量,j为虚数单位;公式(3)中的状态变量的稳态周期解的数学表达式也能够用三角函数的形式表示如下式(4):
xS=XDC+xac=A0+2Re(A1)cosωt-2Im(A1)sinωt+2Re(A2)cos2ωt-2Im(A2)sin2ωt
+2Re(A3)cos3ωt-2Im(A3)sin3ωt+...+2Re(An)cosnωt-2Im(An)sinnωt
(4)
公式(4)中的Re(Ai)(i=2,3,……,n)分别表示复数向量Ai的实部,Im(Ai)分别表示复数向量Ai的虚部;
在步骤S4中,根据Grnwald-Letnikov定义用离散方法求分数阶主振荡微分方程G0(pα,pβ)x0+G1f0m=U的瞬态数值解,再通过Matlab拟合得到瞬态解的近似表达式;
将主振荡微分方程写成如下关于系统状态变量的分数阶微分方程组的形式:
其中,a1,a2,a3,a4,b1,b2,b3为跟具体电路参数有关的常系数,α和β分别表示分数阶电感和分数阶电容的阶次;根据Grnwald-Letnikov定义用离散方法得到分数阶电容电压和电感电流主振荡分量瞬态近似数值解分别为式(6a)和(6b):
式中vC0-trans和iL0-trans分别表示分数阶电容电压和电感电流主振荡分量的瞬态解,h为步长,vt-kh和it-kh表示t-kh时刻电容电压和电感电流的瞬时值,式中能够由递推公式得出,k为自然数,k=1,2,....;将式(6a)和(6b)所得到的数值解用Matlab多项式拟合能够得到主振荡分量瞬态解的近似解析表达式x0-trans=[iL0-trans vC0-trans]T,其中vC0-trans和iL0-trans的近似解析表达式分别为式(7a)、(7b):
式中V0和I0分别表示分数阶电容电压和电感电流主振荡分量瞬态解的直流分量,ωv和ωi分别表示分数阶电容电压和电感电流主振荡分量瞬态解的基波角频率,avk和bvk分别表示分数阶电容电压主振荡分量瞬态解第k次谐波中余弦分量和正弦分量的幅值,aik和bik分别表示分数阶电感电流主振荡分量瞬态解第k次谐波中余弦分量和正弦分量的幅值,自然数k=1,2,3,……n;
在步骤S5中,把S3所得到纹波分量xac与S4所得到的主振荡分量瞬态解x0-trans叠加,即可得到分数阶开关变换器的瞬态解析解的近似表达式xtrans≈x0-trans+xac。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710019197.7A CN106909711B (zh) | 2017-01-11 | 2017-01-11 | 一种求分数阶ccm开关变换器瞬态解的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710019197.7A CN106909711B (zh) | 2017-01-11 | 2017-01-11 | 一种求分数阶ccm开关变换器瞬态解的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106909711A true CN106909711A (zh) | 2017-06-30 |
CN106909711B CN106909711B (zh) | 2020-04-28 |
Family
ID=59206496
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710019197.7A Expired - Fee Related CN106909711B (zh) | 2017-01-11 | 2017-01-11 | 一种求分数阶ccm开关变换器瞬态解的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106909711B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110221645A (zh) * | 2019-06-11 | 2019-09-10 | 广州金升阳科技有限公司 | 一种纹波电流产生电路 |
CN110348057A (zh) * | 2019-06-12 | 2019-10-18 | 华南理工大学 | 快速获得电力电子dc/dc变换器动态和稳态特性的方法 |
CN110572053A (zh) * | 2019-10-28 | 2019-12-13 | 重庆三峡学院 | 一种分数阶d类半波整流器及其参数设计方法 |
CN110611420A (zh) * | 2019-10-28 | 2019-12-24 | 重庆三峡学院 | 一种单相分数阶buck型交流斩波器及其参数设计方法 |
CN112507643A (zh) * | 2020-12-21 | 2021-03-16 | 华南理工大学 | 融合卡尔曼滤波技术的超高频变换器分析方法 |
CN113408126A (zh) * | 2021-06-17 | 2021-09-17 | 华南理工大学 | 一种求分数阶甚高频谐振变换器瞬态解的解耦方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104008231A (zh) * | 2014-05-08 | 2014-08-27 | 华南理工大学 | 一种数字电流模控制的开关功率变换器的建模及分析方法 |
CN104978304A (zh) * | 2015-07-24 | 2015-10-14 | 华南理工大学 | 电流连续模式分数阶开关变换器的符号分析方法及装置 |
WO2016161363A1 (en) * | 2015-04-03 | 2016-10-06 | Avatekh, Inc. | Method and apparatus for regulated three-phase ac-to-dc conversions with high power factor and low harmonic distortions |
CN106227925A (zh) * | 2016-07-12 | 2016-12-14 | 华南理工大学 | 一种电感电流断续模式分数阶开关变换器的符号分析方法 |
-
2017
- 2017-01-11 CN CN201710019197.7A patent/CN106909711B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104008231A (zh) * | 2014-05-08 | 2014-08-27 | 华南理工大学 | 一种数字电流模控制的开关功率变换器的建模及分析方法 |
WO2016161363A1 (en) * | 2015-04-03 | 2016-10-06 | Avatekh, Inc. | Method and apparatus for regulated three-phase ac-to-dc conversions with high power factor and low harmonic distortions |
CN104978304A (zh) * | 2015-07-24 | 2015-10-14 | 华南理工大学 | 电流连续模式分数阶开关变换器的符号分析方法及装置 |
CN106227925A (zh) * | 2016-07-12 | 2016-12-14 | 华南理工大学 | 一种电感电流断续模式分数阶开关变换器的符号分析方法 |
Non-Patent Citations (2)
Title |
---|
杨亚泽: "基于等效小参量法的DC-DC开关变换器建模与仿真", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
林波涛,丘水生: "适于瞬态分析的等效小参量法及其在PWM开关变换器中的应用", 《电子科学学刊》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110221645A (zh) * | 2019-06-11 | 2019-09-10 | 广州金升阳科技有限公司 | 一种纹波电流产生电路 |
CN110348057A (zh) * | 2019-06-12 | 2019-10-18 | 华南理工大学 | 快速获得电力电子dc/dc变换器动态和稳态特性的方法 |
CN110348057B (zh) * | 2019-06-12 | 2020-12-22 | 华南理工大学 | 快速获得电力电子dc/dc变换器动态和稳态特性的方法 |
CN110572053A (zh) * | 2019-10-28 | 2019-12-13 | 重庆三峡学院 | 一种分数阶d类半波整流器及其参数设计方法 |
CN110611420A (zh) * | 2019-10-28 | 2019-12-24 | 重庆三峡学院 | 一种单相分数阶buck型交流斩波器及其参数设计方法 |
CN112507643A (zh) * | 2020-12-21 | 2021-03-16 | 华南理工大学 | 融合卡尔曼滤波技术的超高频变换器分析方法 |
CN113408126A (zh) * | 2021-06-17 | 2021-09-17 | 华南理工大学 | 一种求分数阶甚高频谐振变换器瞬态解的解耦方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106909711B (zh) | 2020-04-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106909711A (zh) | 一种求分数阶ccm开关变换器瞬态解的方法 | |
CN106874548A (zh) | 一种基于双重傅里叶变换分析逆变器的方法 | |
Gateau et al. | Multicell converters: active control and observation of flying-capacitor voltages | |
Merabet et al. | A comparative study of harmonic currents extraction by simulation and implementation | |
Hou et al. | Adaptive fuzzy backstepping control of three-phase active power filter | |
Seltzer et al. | Small signal phasor modeling of dual active bridge series resonant DC/DC converters with multi-angle phase shift modulation | |
CN104022662A (zh) | 基于模型预测控制的pwm整流器控制方法及装置 | |
Shukla et al. | Implementation of empirical mode decomposition based algorithm for shunt active filter | |
CN101242101B (zh) | 电力系统负荷谐波稳态模型的建立方法及仿真方法 | |
CN104978304A (zh) | 电流连续模式分数阶开关变换器的符号分析方法及装置 | |
Hwang et al. | A new modelling framework for power supply networks with converter based loads and generators-the Harmonic State-Space | |
CN106227925B (zh) | 一种电感电流断续模式分数阶开关变换器的符号分析方法 | |
CN103257270A (zh) | 基于傅里叶级数的锯齿载波混沌spwm频谱分析方法 | |
CN103532415A (zh) | 基于ghγ坐标系的四桥臂变流器空间矢量调制方案 | |
CN105322818A (zh) | 一种基于新型模型预测控制的三相pwm整流的控制方法 | |
CN102930093B (zh) | 一种两电平变流器的建模方法 | |
Hogan et al. | Load and source electronic emulation using resonant current control for testing in a microgrid laboratory | |
Li et al. | Improvement of model predictive control for single‐phase inverters by using sinusoidal signal observers | |
CN201374643Y (zh) | 应用于输入不平衡三相升压型交直流转换器的控制装置 | |
CN109787280B (zh) | 基于起始阻抗角的逆变器系统线路阻抗模拟方法 | |
Carrero et al. | Experimental evaluation of a cascade sliding mode-PI controller for a coupled-inductor boost converter | |
Komurcugil et al. | Control strategy for single-phase PWM ac/dc voltage-source converters based on Lyapunov's direct method | |
Vacheva et al. | Application of neural network for smart control of a buck DC/DC converter | |
Josefova et al. | Analysis of current drawn by the voltage-source active rectifier from the electricity network | |
Andriyanov et al. | Mathematical model of a three-phase transistor controlled rectifier with a PI voltage controller in regeneration mode |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200428 |