首页 / 专利库 / 信号处理 / 正弦波 / 一种求分数阶CCM开关变换器瞬态解的方法

一种求分数阶CCM开关变换器瞬态解的方法

阅读:1036发布:2020-05-11

专利汇可以提供一种求分数阶CCM开关变换器瞬态解的方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种求分数阶CCM 开关 变换器瞬态解的方法,利用等效小参量法,将复杂的分数阶微分运算转换为微分算子的指数运算,从而将复杂的分数阶变换器状态变量稳态周期解的求解转化为 直流分量 和各次谐波分量幅度的求解,而直流分量和谐波分量幅度的求解则能够利用谐波平衡法;只需要分析分数阶开关变换器状态变量主振荡方程的瞬态解,再结合其稳态解即可得到分数阶开关变换器状态变量的瞬态解析解。本发明方法能够较好的分析分数阶CCM开关变换器的瞬态过程。,下面是一种求分数阶CCM开关变换器瞬态解的方法专利的具体信息内容。

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,当该受控开关断开时δ=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)表示复数向量Ai的实部,i=2,3,……,n;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。

说明书全文

一种求分数阶CCM开关变换器瞬态解的方法

技术领域

[0001] 本发明涉及分数阶开关变换器的建模与分析领域,尤其是指一种求分数阶CCM开关变换器瞬态解的方法。

背景技术

[0002] 过去针对开关变换器常用的建模与分析方法有:离散迭代映射模型、基于状态空间平均法的模型、基于电路原理(KCL、KVL)的分段线性模型[1]以及结合谐波平衡与扰动法的等效小参量法[2],这些方法的分析对象均是整数阶的开关变换器,即变换器中的电容、电感都是整数阶的元件,然而文献[3][4]指出实际上电容和电感在本质上均是分数阶的,这就需要为变换器建立相应的分数阶模型。
[0003] 谐波平衡法是分析非线性系统周期解的常用方法之一,但是它只适用求系统的稳态解,文献[5]提出的通用平均法将谐波平衡法推广到瞬态分析,但是它不易求得高次谐波。文献[6]用结合谐波平衡与扰动法的等效小参量法分析了开关变换器的瞬态过程,但并未考虑电容和电感的分数阶特性。文献[7]考虑了电感电容的分数阶特性,用等效小参量法分析了分数阶开关变换器的稳态特性,但未分析分数阶开关变换器的瞬态过程。
[0004] 参考文献
[0005] [1]罗晓曙.DC-DC变换器的非线性动学行为与混沌控制[M].科学出版社,2012,pp.23-24.
[0006] [2]Qiu S S,FilanovskI M.Calculation of steady-state oscillations in non linear circuits[J],Int.J.Electronics,1989,67(3),pp.403-414.
[0007] [3]Westerlund S.Dead Matter Has Memory![M].Kalmar,Sweden:Causal Consulting,2002,Chap.7.
[0008] [4]Podlubny I.Fractional Differential Equations[M].San Diego:Academic Press,1999,Chap.2.
[0009] [5]Sanders SR,Noworolski J,Liu X.et al.Generalized averaging method for power conversion circuits.IEEE Trans.on Power Electronics,1991,PE-6(2);251-259.
[0010] [6]林波涛,邱生.适于瞬态分析的等效小参量法及其在PWM开关变换器中的应用[J].电子科学学刊,1997,19(5),pp.649-657.
[0011] [7]Xi Chen,Yanfeng Chen,Bo Zhang,Qiuyuan Dong.A Modeling and Analysis Method for Fractional-oder DC-DC Converters.IEEE Trans.on Power Electronics,2016,Vol.pp,no.99,pp.1-1.

发明内容

[0012] 本发明的目的是在于克服现有人机交互方式的不足,提出了一种求分数阶CCM开关变换器瞬态解的方法,能够较好的分析分数阶CCM开关变换器的瞬态过程。
[0013] 为实现上述目的,本发明所提供的技术方案为:一种求分数阶CCM开关变换器瞬态解的方法,利用等效小参量法,将复杂的分数阶微分运算转换为微分算子的指数运算,从而将复杂的分数阶变换器状态变量稳态周期解的求解转化为直流分量和各次谐波分量幅度的求解,而直流分量和谐波分量幅度的求解则能够利用谐波平衡法;只需要分析分数阶开关变换器状态变量主振荡方程的瞬态解,再结合其稳态解即可得到分数阶开关变换器状态变量的瞬态解析解;其包括以下步骤:
[0014] S1、建立以分数阶微分方程描述的分数阶开关变换器的非线性数学模型;
[0015] S2、利用等效小参量法得到分数阶开关变换器的等效数学模型
[0016] 利用等效小参量法求解S1中的非线性数学模型,得到描述分数阶开关变换器的等效数学方程组,即分数阶开关变换器的等效数学模型;该等效数学方程组包含一个求系统状态变量主振荡分量的主振荡分数阶微分方程,和一系列求状态变量修正量的分数阶微分方程;
[0017] S3、利用谐波平衡法求分数阶开关变换器系统状态变量的稳态周期解
[0018] 利用谐波平衡法逐步求解S2中等效数学方程组中的各个分数阶微分方程的稳态解,得到分数阶开关变换器状态变量稳态周期解的近似解表达式,所得到的近似周期解包含直流分量和纹波分量,其中纹波分量由基波和各次谐波组成;
[0019] S4、求分数阶开关变换器主振荡分量瞬态解
[0020] 根据分数阶微积分定义求S2中主振荡分数阶微分方程瞬态解,得到分数阶开关变换器状态变量的主振荡分量的瞬态解;
[0021] S5、获取分数阶开关变换器系统状态变量瞬态解的解析表达式
[0022] 把S3所得到纹波分量与S4所得到的主振荡分量的瞬态解叠加,即得到分数阶开关变换器的瞬态解的解析表达式。
[0023] 在步骤S1中,所建立的以分数阶微分方程描述的分数阶开关变换器的非线性数学模型如下:
[0024] G0(pα,pβ)x+G1f(x)=U  (1)
[0025] 式中,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,当该受控开关断开时δ=0;U为变换器输入电压向量;
[0026] 在步骤S2中,所述的等效小参量法的具体步骤如下:
[0027] S21、将开关函数δ表示成主振荡分量和修正量之和的级数形式:其中δ0表示开关函数的主振荡分量,δi表示开关函数的第i阶修正量,它们能够根据具体开关函数的傅里叶级数来确定;
[0028] S22、将待求解的状态变量x也表示成级数形式: 其中x0表示状态变量的主振荡分量,xi表示状态变量的第i阶修正量,它们在具体的求解过程逐步确定;
[0029] S23、将δ和x的级数表达式代入非线性矢量函数f(x)=δx中,得到其中f0表示非线性矢量函数的主振荡分量,fi表示非线性矢量函数的第i
阶修正量;
[0030] 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具有不同频率成分的项;
[0031] S25、将f0=f0m+εR1和fi=fim+εRi+1代入 中,从而能够将f(x)表示为
[0032] 在上述步骤S21~S25中,上标或下标i是一个整数,i=1,2,……;ε是引入的一个小量标记,εixi表明xi是状态变量x的第i阶小量,并且有εi+1xi+1<εixi<x0,当在运算过程中需要具体数值时ε=1;
[0033] S26、将 和 代入到公式(1)中,并令等式两边具有相同εi项的项分别相等,能够得到描述分数阶CCM开关变换器的等效数学模型,如下式(2):
[0034]
[0035] 式(2)中的第1个分数阶微分方程用于求状态变量的主振荡分量x0,称为主振荡方程;第2~n个分数阶微分方程用于求状态变量的各阶修正量xi(i=1,2,……n),称为修正量方程;
[0036] 在步骤S3中,以指数函数表示状态变量的稳态周期解xS的近似数学表达式如下:
[0037]
[0038] 式中,直流分量XDC=A0为分数开关变换器的状态变量稳态主振荡分量;xac为纹波分量,其中A1为基波的的幅度向量, 为其共轭;Ai(i=2,3,……,n)为第i次波的的幅度向量, 为其共轭;ω为开关变换器的频率,t表示时间变量,j为虚数单位;公式(3)中的状态变量的稳态周期解的数学表达式也能够用三角函数的形式表示如下式(4):
[0039] 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)[0040] 公式(4)中的Re(Ai)(i=2,3,……,n)表示复数向量Ai的实部,Im(Ai)表示复数向量Ai的虚部
[0041] 在步骤S4中,根据Grnwald-Letnikov定义用离散方法求分数阶主振荡微分方程G0(pα,pβ)x0+G1f0m=U的瞬态数值解,再通过Matlab拟合得到瞬态解的近似表达式;
[0042] 将主振荡微分方程写成如下关于系统状态变量的分数阶微分方程组的形式:
[0043]
[0044] 其中,a1,a2,a3,a4,b1,b2,b3为跟具体电路参数有关的常系数,α和β分别表示分数阶电感和分数阶电容的阶次;根据Grnwald-Letnikov定义用离散方法得到分数阶电容电压和电感电流主振荡分量瞬态近似数值解分别为式(6a)和(6b):
[0045]
[0046]
[0047] 式中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):
[0048]
[0049]
[0050] 式中V0和I0分别表示分数阶电容电压和电感电流主振荡分量瞬态解的直流分量,ωv和ωi分别表示分数阶电容电压和电感电流主振荡分量瞬态解的基波角频率,avk和bvk分别表示分数阶电容电压主振荡分量瞬态解第k次谐波中余弦分量和正弦分量的幅值,aik和bik分别表示分数阶电感电流主振荡分量瞬态解第k次谐波中余弦分量和正弦分量的幅值,自然数k=1,2,3,……n;
[0051] 在步骤S5中,把S3所得到纹波分量xac与S4所得到的主振荡分量瞬态解x0-trans叠加,即可得到分数阶开关变换器的瞬态解析解的近似表达式xtrans≈x0-trans+xac。
[0052] 本发明与现有技术相比,具有如下优点与有益效果:
[0053] 1、由本发明所提方法的求解公式可以看出,采用本方法求分数阶开关变换器状态变量的瞬态解,相当于将求解非整数阶微积分运算的复杂过程转化为矩阵运算和求线性方程(组)的过程,计算量小,结果直观,清晰。
[0054] 2、通过与PSIM电路仿真波形的对比验证,该发明所提供的方法能够较好的分析分数阶CCM开关变换器的瞬态过程,便于电路的优化控制与设。附图说明
[0055] 图1为CCM分数阶Buck变换器的电路原理图。
[0056] 图2a为本发明方法得到的电感电流波形与PSIM电路仿真波形的对比图。
[0057] 图2b为本发明方法得到的电容电压波形与PSIM电路仿真波形的对比图。

具体实施方式

[0058] 下面结合具体实施例对本发明作进一步说明。
[0059] 本实施例所提供的求分数阶CCM开关变换器瞬态解的方法,具体包括以下步骤:
[0060] S1、建立以分数阶微分方程描述的分数阶开关变换器的非线性数学模型:
[0061] G0(pα,pβ)x+G1f(x)=U  (1)
[0062] 上式中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,当该受控开关断开时δ=0;U为变换器输入电压向量。
[0063] S2、利用等效小参量法得到分数阶开关变换器的等效数学模型
[0064] 具体为:利用等效小参量法求解分数阶CCM开关变换器的非线性数学模型,得到描述该分数阶开关变换器的等效数学方程组,即等效数学模型。该等效数学方程组包含一个求分数阶CCM开关变换器状态变量主振荡分量的分数阶微分方程,和一系列求状态变量各阶修正量的分数阶微分方程。
[0065] 所述的等效小参量法的具体步骤如下:
[0066] S21、将开关函数δ表示成主振荡分量和修正量之和的级数形式:其中δ0表示开关函数的主振荡分量,δi表示开关函数的第i阶修正量,它们可根据具体开关函数的傅里叶级数来确定;
[0067] S22、将待求解的状态变量x也表示成如下级数形式: 其中x0表示状态变量的主振荡分量,xi表示状态变量的第i阶修正量,它们在具体的求解过程逐步确定;
[0068] S23、将δ和x的级数表达式代入非线性矢量函数f(x)=δx中,得到其中f0表示非线性矢量函数的主振荡分量,fi表示非线性矢量函数的第i
阶修正量;
[0069] 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具有不同频率成分的项;
[0070] S25、将f0=f0m+εR1和fi=fim+εRi+1代入 中,从而可将f(x)表示为
[0071] 在上述步骤S21~S25中,上标或下标i是一个整数,i=1,2,……;ε是本方法所引入的一个小量标记,例如εixi表明xi是状态变量x的第i阶小量,并且有εi+1xi+1<εixi<x0,当在运算过程中需要具体数值时ε=1;
[0072] S26、将 和 代入到公式(1)中,并令等式两边具有相同εi项的项分别相等,可得到描述分数阶CCM开关变换器的等效数学模型,如下式(2):
[0073]
[0074] 上式(2)中所述的第1个分数阶微分方程用于求状态变量的主振荡分量x0,称为主振荡方程;第2~n个分数阶微分方程用于求状态变量的各阶修正量xi(i=1,2,……n),称为修正量方程。
[0075] S3、利用谐波平衡法求分数阶开关变换器系统状态变量的稳态周期解
[0076] 利用谐波平衡法逐步求解公式(2)所述的等效数学方程模型中的各个分数阶微分方程的稳态解,得到分数阶开关变换器状态变量稳态周期解的近似解表达式,所得到的近似周期解包含直流分量和纹波分量,其中纹波分量由基波和各次谐波组成。所得到的以指数函数表示状态变量的稳态周期解xS的近似数学表达式如下:
[0077]
[0078] 式中直流分量XDC=A0为分数开关变换器的状态变量稳态主振荡分量;xac为纹波分量,其中A1为基波的的幅度向量, 为其共轭;Ai(i=2,3,……,n)为第i次波的的幅度向量, 为其共轭;ω为开关变换器的角频率,t表示时间变量,j为虚数单位。公式(3)中的状态变量的稳态周期解的数学表达式也可以用三角函数的形式表示如下式(4):
[0079] 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)[0080] 公式(4)中的Re(Ai)(i=2,3,……,n)表示复数向量Ai的实部,Im(Ai)表示复数向量Ai的虚部。
[0081] S4、求分数阶开关变换器主振荡分量瞬态解
[0082] 根据Grnwald-Letnikov定义用离散方法求分数阶主振荡微分方程G0(pα,pβ)x0+G1f0m=U的瞬态数值解,再通过Matlab拟合得到瞬态解的近似表达式。
[0083] 将主振荡微分方程写成如下关于系统状态变量的分数阶微分方程组的形式:
[0084]
[0085] 其中a1,a2,a3,a4,b1,b2,b3为跟具体电路参数有关的常系数,α和β分别表示分数阶电感和分数阶电容的阶次。根据Grnwald-Letnikov定义用离散方法可得到分数阶电容电压和电感电流主振荡分量瞬态近似数值解分别为式(6a)和(6b):
[0086]
[0087]
[0088] 式中,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):
[0089]
[0090]
[0091] 式中,V0和I0分别表示分数阶电容电压和电感电流主振荡分量瞬态解的直流分量,ωv和ωi分别表示分数阶电容电压和电感电流主振荡分量瞬态解的基波角频率,avk和bvk分别表示分数阶电容电压主振荡分量瞬态解第k次谐波中余弦分量和正弦分量的幅值,aik和bik分别表示分数阶电感电流主振荡分量瞬态解第k次谐波中余弦分量和正弦分量的幅值,自然数k=1,2,3,……n。
[0092] S5、获取分数阶开关变换器系统状态变量瞬态解
[0093] 把S3所得到纹波分量xac与S4所得到的主振荡分量瞬态解x0-trans叠加,即得到分数阶开关变换器的瞬态解析解的近似表达式xtrans≈x0-trans+xac。
[0094] 下面结合附图1、2a、2b对本实施例上述方法进行具体说明。
[0095] 图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分别为分数阶电感电流和电容电压。
[0096] 根据上述方法步骤可求得分数阶Buck变换器状态变量稳态周期解纹波分量的近似解析表达式为xac=[iLac vCac]T,其中iLac和vCac分别为分数阶电感电流和电容电压纹波,它们分别由式(8a)和(8b)表示:
[0097] iLac≈-0.1028cosωt-0.0148sinωt-0.0076cos2ωt+0.0015sin2ωt
[0098] -0.0056cos3ωt-0.0056sin3ωt  (8a)
[0099] vCac=-0.0596cosωt-0.3516sinωt-0.0308cos2ωt-0.0072sin2ωt
[0100] +0.0054cos3ωt-0.0086sin3ωt  (8b)
[0101] 分数阶CCM-Buck变换器主振荡微分方程可表示为式(9):
[0102]
[0103] 根据步骤S4可得到主振荡分量近似瞬态解如式(7a)和(7b)所示,其中V0=-2.03V,ωv=218rad/s,I0=16.16A,ωi=158rad/s,其它各系数如下表1所示。
[0104] 表1
[0105]
[0106]
[0107] 根据步骤S5可得分数阶CCM-Buck变换器状态变量的瞬态解析解。
[0108] 将本发明方法分别与PSIM电路仿真得到的瞬态电感电流、电容电压波形进行比较,如图2a、2b仿真结果对比验证图所示,图中实线为本发明得到的波形,虚线为PSIM电路仿真得到的波形。从图中可以看出两条曲线相符合,说明本发明所提出的方法是有效的。
[0109] 以上所述实施例只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈