一种电感电流断续模式分数阶开关变换器的符号分析方法技术

技术编号:14195522 阅读:117 留言:0更新日期:2016-12-15 15:36
本发明专利技术公开了一种电感电流断续模式分数阶开关变换器的符号分析方法,该方法结合了谐波平衡的原理,通过将变换器中关于状态变量分数阶次的微分运算转换为微分算子,并将所有微分算子合并为对角符号矩阵,从而将求解非整数阶微积分运算的过程转化为矩阵运算和线性方程(组)求解的过程,相比较已有的针对分数阶开关变换器常用的在Matlab/Simulink中建立Oustaloup滤波器近似模型的分析方法,本发明专利技术方法除了能够解析地分析变换器状态变量纹波峰峰值大小、储能元件阶次变化对变换器工作状态的影响,还能够快速地获得状态变量的稳态周期解析解,并可以用于分析状态变量的谐波成分。

【技术实现步骤摘要】

本专利技术涉及分数阶开关变换器的建模与分析领域,尤其是指一种电感电流断续模式(DCM,Discontinuous-conduction Mode)分数阶开关变换器的符号分析方法。
技术介绍
过去针对开关变换器常用的建模与分析方法有:基于状态空间平均法的模型、离散迭代映射模型、基于电路原理(KCL、KVL)的分段线性模型以及结合谐波平衡与扰动法的等效小参量法,这些方法的分析对象均是整数阶的开关变换器,即变换器中的电容、电感都是整数阶的元件,然而现有参考文献1“Westerlund S.Dead Matter Has Memory[M].Kalmar,Sweden:Causal Consulting,2002,Chap.7.”及参考文献2“Podlubny I.Fractional Differential Equations[M].San Diego:Academic Press,1999,Chap.2.”指出实际上电容和电感在本质上均是分数阶的,这就需要为变换器建立相应的分数阶模型。现有参考文献3“王发强,马西奎.基于分数阶微积分的电感电流断续模式下Boost变换器的建模与分析[J].中国科学:技术科学,2013,43(4).368-374”同时考虑电感与电容的分数阶特性,首先建立了电感电流断续模式(DCM,Discontinuous-conduction Mode)下分数阶Boost变换器的状态空间平均模型,并基于Oustaloup滤波器的非整数阶频域逼近微积分算法,在Matlab/Simulink环境下建立了仿真模型(如图1、2所示),对分数阶DC-DC变换器随阶次变化的工作特性进行了初步的分析。依据这一思路,现有技术(比如参考文献4“王发强,马西奎.电感电流连续模式下Boost变换器的分数阶建模与仿真分析[J].物理学报,2011,60(7).070506-1–070506-8”等)分 别研究了电感电流连续以及伪连续工作模式下的分数阶开关变换器,图3和图4是分别通过图1和图2中建立的模型获得的电容电压和电感电流波形,已有的技术是通过在Matlab/Simulink中建立模块化模型的方式展示分数阶开关变换器的工作特性,并通过仿真波形的方式展示稳态下变换器状态变量的纹波;这种方法不能得到状态变量的稳态周期解析解,难以解析地分析纹波峰峰值大小。
技术实现思路
本专利技术的目的在于克服现有技术的缺点与不足,提供一种电感电流断续模式(DCM,Discontinuous-conduction Mode)分数阶开关变换器的符号分析方法,能够快速获得工作在DCM状态下分数阶开关变换器状态变量稳态周期解析解。为实现上述目的,本专利技术所提供的技术方案为:一种电感电流断续模式分数阶开关变换器的符号分析方法,包括以下步骤:1)建立分数阶开关变换器的数学模型工作在电感电流断续模式DCM状态下分数阶DC-DC变换器的系统状态描述为:式中,x=[iLm vCm]T表示系统的状态变量,包括第m个电感上的电流iLm、第m个电容上的电压vCm,p表示相应的电感L、电容C元件上对应的系统状态变量的阶次,A0和B0分别表示不受开关函数影响的系数矩阵,A1B1和A3B3分别表示受开关函数影响的系数矩阵;开关函数δ(1)和δ(3)定义为:其中,变换器在开环工作时占空比D1和D3为固定值,同时,令非线性部分为:f(q)=δ(q)(Aqx+Bq)τ=ωt,其中然后将针对状态变量的微积分运算转换为关于微积分算子的代数运算,即 由于变换器中存在多个状态变量,故将每个状态变量相应的微积分算子合并为微分算子对角符号矩阵矩阵κ中α、β这些元素用于表示不同状态变量的分数阶微积分阶次,当L、C均为整数阶时,κ=-I,I为单位矩阵,其中的+/-号分别表示对状态变量求积分/微分;通过将微分运算转换为关于微分算子的代数运算,能够将分数阶开关变换器的数学模型如下所示:式(3)中,G0为所有包含微分算子对角符号矩阵κ的Gki组成的列矩阵,k∈Eir表示当前第i阶修正量中谐波次数k,i、k的定义后同, 从Gki的形式能够体现出分数阶次对状态变量解析解的影响;将状态变量x以及开关函数δ(q)均展开为主部与小量余项之和的形式:将上式代入f(q)=δ(q)(Aqx+Bq),合并相同阶次余项小量,得:其中:式中,用表示所述状态变量x第i阶修正量的主部,用表示所述状态变量x第i阶修正量的余项小量;根据谐波平衡原理,将所述状态变量x与开关函数δ(q)的展开式(4)中主部和第i阶余项小量做傅里叶展开如下:其中aki表示第i阶修正量的k次谐波成分的幅值,所述开关函数δ(q)展开式系数表达式为:其中依据谐波平衡原理,将系数展开式(8)代入傅里叶展开式(7),依次求解状态变量的主振荡分量和各阶修正量;2)求主振荡分量首先,求解状态变量的主振荡分量,通常主振荡中只含有直流量,故设为:x0=a00=[I00 V00]T (9)当k=0,即G0=G00=A0,x0代入式(6)中再代入(4)式中,得:G00·x0+b0(A1x0+B1)+c0(A3x0+B3)+B0=0 (10)由式(10)求得变换器状态变量的主振荡分量:3)求各阶修正量根据主振荡分量余项R1中含有的谐波成分,设状态变量的一阶修正量形式如下:其中,a11=[I11 V11]T,c.c表示共轭项,后同;由状态变量的一阶修正量中的 谐波成分可知k∈E1r,代入式(6)中f1,得到一阶修正量表达式:Gk1·x1+(b0A1x1+b1B1+b1A1x0+b1B1)+(c0A3x1+c1B3+c1A3x0+c1B3)+B0=0 (13)根据式(13)能够获得关于谐波幅值a01和ak1的线性方程组;将参数代入所得当前阶次修正量的表达式,若当前阶次修正量的各次谐波幅值相比较上一阶修正量小于一个数量级,则不需做更高阶的修正,反之,继续依据上述过程继续求更高阶次的修正量;4)将主振荡分量和各阶修正量相加,获得关于状态变量的稳态周期解析解表达式。本专利技术与现有技术相比,具有如下优点与有益效果:由本专利技术所提方法的求解公式可以看出,采用本方法求分数阶开关变换器状态变量的稳态周期解析解,相当于将求解非整数阶微积分运算的复杂过程转化为矩阵运算和求线性方程(组)的过程,只要根据电路原理建立如公式(3)矩阵形式的分数阶开关变换器状态方程,然后将系数表达式代入各阶修正量公式,通过简单的矩阵乘除加减运算和消元就可以得到关于分数阶变换器状态变量稳态解的表达式。相比较过去纯数学领域提出的各类分数阶微积分方程的求解方法,本专利技术所提方法的求解过程结合了开关变换器的特点,避开了针对分数阶微积分运算原理的深入讨论,所得的解具有明显的物理意义,根据采用本专利技术所提方法获得的稳态解的形式,可以清楚的看到状态变量所包含的谐波成分,有利于对该类变换器展开更深入的分析。附图说明图1为参考文献3中在Matlab/Simulink中基于Oustaloup滤波器法建立的开环分数阶Boost变换器仿真模型。图2为图1中封装的基于分数阶频域近似方法的Oustaloup滤波器子系统。图3为图1中开环分数阶Boost变换器仿真模型输出电容电压的逐周期仿真波形。其中,横坐标为仿真的时间,本文档来自技高网
...
一种电感电流断续模式分数阶开关变换器的符号分析方法

【技术保护点】
一种电感电流断续模式分数阶开关变换器的符号分析方法,其特征在于,包括以下步骤:1)建立分数阶开关变换器的数学模型工作在电感电流断续模式DCM状态下分数阶DC‑DC变换器的系统状态描述为:dpxdtp=(A0+B0)+δ(1)(A1x+B1)+δ(3)(A3x+B3)---(1)]]>式中,x=[iLm vCm]T表示系统的状态变量,包括第m个电感上的电流iLm、第m个电容上的电压vCm,p表示相应的电感L、电容C元件上对应的系统状态变量的阶次,A0和B0分别表示不受开关函数影响的系数矩阵,A1B1和A3B3分别表示受开关函数影响的系数矩阵;开关函数δ(1)和δ(3)定义为:δ(1)=1,0<t≤D1T0,D1T<t≤T---(2.a)]]>δ(3)=0,0<t≤(1-D3)T1,(1-D3)T<t≤T---(2.b)]]>其中,变换器在开环工作时占空比D1和D3为固定值,同时,令非线性部分为:f(q)=δ(q)(Aqx+Bq)τ=ωt,其中然后将针对状态变量的微积分运算转换为关于微积分算子的代数运算,即由于变换器中存在多个状态变量,故将每个状态变量相应的微积分算子合并为微分算子对角符号矩阵矩阵κ中α、β这些元素用于表示不同状态变量的分数阶微积分阶次,当L、C均为整数阶时,κ=‑I,I为单位矩阵,其中的+/‑号分别表示对状态变量求积分/微分;通过将微分运算转换为关于微分算子的代数运算,能够将分数阶开关变换器的数学模型如下所示:G0(κ)x+f(1)+f(3)+B0=0G0(κ)=A0-κ---(3)]]>式(3)中,G0为所有包含微分算子对角符号矩阵κ的Gki组成的列矩阵,k∈Eir表示当前第i阶修正量中谐波次数k,i、k的定义后同,从Gki的形式能够体现出分数阶次对状态变量解析解的影响;将状态变量x以及开关函数δ(q)均展开为主部与小量余项之和的形式:x=x0+Σi=1nϵixiδ(1)=δ0(1)+Σi=1nϵiδi(1)δ(3)=δ0(3)+Σi=1nϵiδi(3)---(4)]]>将上式代入f(q)=δ(q)(Aqx+Bq),合并相同阶次余项小量,得:f(q)=f0(q)+ϵf1(q)+ϵ2f2(q)+...+ϵnfn(q)---(5)]]>其中:f0(q)=δ0(q)(Aqx0+Bq)=f0m(q)+R1(q)f1(q)=δ0(q)x1+δ1(q)(Aqx0+Bq)=f1m(q)+R2(q)f2(q)=Aq(δ0(q)x2+δ1(q)x1)+δ2(q)(Aqx0+Bq)=f2m(q)+R3(q)...---(6)]]>式中,用表示所述状态变量x第i阶修正量的主部,用表示所述状态变量x第i阶修正量的余项小量;根据谐波平衡原理,将所述状态变量x与开关函数δ(q)的展开式(4)中主部和第i阶余项小量做傅里叶展开如下:xi=a0i+Σk∈Eir(akiejkτ+a‾kie-jkτ)δi(1)=b0+Σk∈Eir(bkiejkτ+b‾kie-jkτ)δi(3)=c0+Σk∈Eir(ckiejkτ+c‾kie-jkτ)---(7)]]>其中aki表示第i阶修正量的k次谐波成分的幅值,所述开关函数δ(q)展开式系数表达式为:b0=1T∫0Tδ(1)dt=D1bki=12(αki-jβki)c0=1T∫0Tδ(3)dt=D3cki=12(αki-jβki)---(8)]]>其中依据谐波平衡原理,将系数展开式(8)代入傅里叶展开式(7),依次求解状态变量的主振荡分量和各阶修正量;2)求主振荡分量首先,求解状态变量的主振荡分量,通常主振荡中只含有直流量,故设为:x0=a00=[I00 V00]T                                 (9)当k=0,即G0=G00=A0,x0代入式(6)中再代入(4)式中,得:G00·x0+b0(A1x0+B1)+c0(A3x0+B3)+B0=0               (10)由式(10)求得变...

【技术特征摘要】
1.一种电感电流断续模式分数阶开关变换器的符号分析方法,其特征在于,包括以下步骤:1)建立分数阶开关变换器的数学模型工作在电感电流断续模式DCM状态下分数阶DC-DC变换器的系统状态描述为: d p x dt p = ( A 0 + B 0 ) + δ ( 1 ) ( A 1 x + B 1 ) + δ ( 3 ) ( A 3 x + B 3 ) - - - ( 1 ) ]]>式中,x=[iLm vCm]T表示系统的状态变量,包括第m个电感上的电流iLm、第m个电容上的电压vCm,p表示相应的电感L、电容C元件上对应的系统状态变量的阶次,A0和B0分别表示不受开关函数影响的系数矩阵,A1B1和A3B3分别表示受开关函数影响的系数矩阵;开关函数δ(1)和δ(3)定义为: δ ( 1 ) = 1 , 0 < t ≤ D 1 T 0 , D 1 T < t ≤ T - - - ( 2. a ) ]]> δ ( 3 ) = 0 , 0 < t ≤ ( 1 - D 3 ) T 1 , ( 1 - D 3 ) T < t ≤ T - - - ( 2. b ) ]]>其中,变换器在开环工作时占空比D1和D3为固定值,同时,令非线性部分为:f(q)=δ(q)(Aqx+Bq)τ=ωt,其中然后将针对状态变量的微积分运算转换为关于微积分算子的代数运算,即由于变换器中存在多个状态变量,故将每个状态变量相应的微积分算子合并为微分算子对角符号矩阵矩阵κ中α、β这些元素用于表示不同状态变量的分数阶微积分阶次,当L、C均为整数阶时,κ=-I,I为单位矩阵,其中的+/-号分别表示对状态变量求积分/微分;通过将微分运算转换为关于微分算子的代数运算,能够将分数阶开关变换器的数学模型如下所示: G 0 ( κ ) x + f ( 1 ) + f ( 3 ) + B 0 = 0 G 0 ( κ ) = A 0 - κ - - - ( 3 ) ]]>式(3)中,G0为所有包含微分算子对角符号矩阵κ的Gki组成的列矩阵,k∈Eir表示当前第i阶修正量中谐波次数k,i、k的定义后同,从Gki的形式能够体现出分数阶次对状态变量解析解的影响;将状态变量x以及开关函数δ(q)均展开为主部与小量余项之和的形式: x = x 0 + Σ i = 1 n ϵ i x i δ ( 1 ) = δ 0 ( 1 ) + Σ i = 1 n ϵ i δ i ( 1 ) δ ( 3 ) = δ 0 ( 3 ) + Σ i = 1 n ϵ i δ i ( 3 ) - - - ( 4 ) ]]>将上式代入f(q)=δ(q)(Aqx+Bq),合并相同阶次余项小量,得: f ( q ) = f 0 ( q ) + ϵf 1 ( q ) + ϵ 2 f 2 ( q ) + ... + ϵ n f n ( q ) - - - ( 5 ) ]]>其中: f 0 ( q ) = δ 0 ( q ) ( A q x 0 + B q ) = f 0 m ( q ) + R 1 ( q ) f 1 ( q ) = δ 0 ( q ) x 1 + δ 1 ( q ) ( A q x 0 + B q ) = f 1 m ( q ) + R 2 ( q ) f 2 ( q ) = A q ( δ 0 ( q ) x 2 + δ 1 ( q ) x 1 ) + δ 2 ( q ) ( A q x 0 + B q ) = f 2 m ( q ) + R 3 ( q ) ... - - - ( 6 ) ]]>式中,用表示所述状态变量x第i阶修正量的主部,用表示所述状态变量x第i阶修正量的余项小量...

【专利技术属性】
技术研发人员:陈艳峰胡劼陈曦张波
申请(专利权)人:华南理工大学
类型:发明
国别省市:广东;44

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1