CN111799795B - 一种基于干扰观测的有源电力滤波器自适应滑模控制 - Google Patents

一种基于干扰观测的有源电力滤波器自适应滑模控制 Download PDF

Info

Publication number
CN111799795B
CN111799795B CN202010573671.2A CN202010573671A CN111799795B CN 111799795 B CN111799795 B CN 111799795B CN 202010573671 A CN202010573671 A CN 202010573671A CN 111799795 B CN111799795 B CN 111799795B
Authority
CN
China
Prior art keywords
hidden layer
sliding mode
vector
derivative
active power
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
CN202010573671.2A
Other languages
English (en)
Other versions
CN111799795A (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.)
Changzhou Campus of Hohai University
Original Assignee
Changzhou Campus of Hohai University
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 Changzhou Campus of Hohai University filed Critical Changzhou Campus of Hohai University
Priority to CN202010573671.2A priority Critical patent/CN111799795B/zh
Publication of CN111799795A publication Critical patent/CN111799795A/zh
Application granted granted Critical
Publication of CN111799795B publication Critical patent/CN111799795B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/01Arrangements for reducing harmonics or ripples
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/20Active power filtering [APF]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种基于干扰观测的有源电力滤波器自适应滑模控制,包括以下步骤:S1,定义开关函数H并根据开关函数H的定义以及基尔霍夫定理建立单相有源电力滤波器动力学方程,从而得到x的二阶导数
Figure DDA0002550277230000011
的通用数学模型;S2,定义有限时间干扰观测器的微分方程;S3,定义跟踪误差及其一阶导数、快速非奇异终端滑模面,通过滑模面求取等效控制项,通过干扰观测器微分方程求取扰动补偿项,得到理想的快速非奇异终端滑模控制律;S4,构造双隐层递归神经网络,通过网络输出逼近切换控制项的切换增益,最终实现自适应快速非奇异终端滑模控制律,并进行稳定证明。本发明能够快速实现无静差跟踪,并达到较低的电网电流畸变率。

Description

一种基于干扰观测的有源电力滤波器自适应滑模控制
技术领域
本发明涉及一种基于干扰观测的有源电力滤波器自适应滑模控制,属于智能控制技术领域。
背景技术
由于生产和生活中的绝大部分能量来源都是电网系统,因此保证电网系统的安全运行是有必要的。在实际的电网系统中,电网往往会受到用户负载侧的影响,也就是说如果负载侧存在大量非线性装置,那么不可避免的会产生谐波电流,而谐波电流一旦并入电网中,会产生诸如谐振、电路损耗、影响附近设备的通信质量等危害,严重时甚至危害电网系统的正常运作。目前,谐波补偿最经济的方式是在用户侧的并网点进行谐波补偿,这样做可以保证公共电网系统中依然是正弦电流,其他用户将不会受到公共电网中的谐波危害。
有源电力滤波器是目前更为可靠的一种谐波补偿装置,在实际应用中它不仅能补偿谐波电流还能补偿无功功率。有源电力滤波器在国外已经得到广泛的使用,然而在我国有源电力滤波器的工业产品仍然较少,这是因为有源电力滤波器的研究难点主要包含两个:第一,如何精确测量电网系统中的参考谐波电流含量;第二,如何设计电流控制器,以实现高精度的电流跟踪效果,即谐波补偿电流跟踪谐波参考电流。无功功率谐波检测算法是目前应用最多的谐波电流检测算法,它已经被证实拥有很好的谐波检测能力。然而电流控制器的设计却是多种多样的,目前没有公认的较好的电流控制算法,而在工业领域应用最多的滞环控制和PI控制控制算法尽管实现简单,但补偿效果不佳,无法发挥有源电力滤波器的优势。
发明内容
为解决传统控制算法的电流控制精度不高的问题,本发明提供一种基于干扰观测的有源电力滤波器自适应滑模控制,够快速实现无静差跟踪,并达到较低的电网电流畸变率。
本发明中主要采用的技术方案为:
一种基于干扰观测的有源电力滤波器自适应滑模控制,包括以下步骤:
S1,定义开关函数H来表示有源电力滤波器物理模型中开关管的导通情况,根据开关函数H的定义以及基尔霍夫定理建立单相有源电力滤波器动力学方程,然后根据动力学方程的状态变量i为x,得到x的二阶导数
Figure BDA0002550277210000021
的通用数学模型;
S2,为了补偿系统中的未知扰动影响,定义有限时间干扰观测器的微分方程;
S3,定义跟踪误差及其一阶导数,然后定义快速非奇异终端滑模面,通过滑模面求取等效控制项,通过干扰观测器微分方程求取扰动补偿项,并考虑切换控制项,将等效控制项、切换控制项和未知扰动补偿项相加得到理想的快速非奇异终端滑模控制律;
S4,构造双隐层递归神经网络,通过网络输出逼近切换控制项的切换增益,最终实现自适应快速非奇异终端滑模控制律,并进行稳定证明。
优选地,所述步骤S1的具体步骤如下:
S1-1:单相有源电力滤波器物理模型包含四个电力电子开关管,根据电路理论,开关管一共两种情况,因此开关函数H具体定义为:
Figure BDA0002550277210000031
S1-2:根据定义的开关函数H以及基尔霍夫定理,建立单相有源电力滤波器动力学方程如下所示:
Figure BDA0002550277210000032
其中,i表示谐波补偿电流,UMN=Udc*H为有源电力滤波器交流侧电压;
S1-3:定义状态变量x1=i,x2=di/dt,在考虑和时间相关的不确定性扰动函数g(t)有界的情况下,单相有源电力滤波器二阶数学模型如下所示:
Figure BDA0002550277210000033
其中,
Figure BDA0002550277210000034
是已知标称函数;在标称函数中,L是交流侧线路总电感,R是交流侧线路总电阻,Udc是直流侧电压,Us是电网电压,控制变量为u=H,不确定扰动界限满足
Figure BDA0002550277210000035
优选地,所述步骤S2中定义的有限时间干扰观测器的微分方程如下所示:
Figure BDA0002550277210000036
其中,Ψ=f(x)+bu是单相有源电力滤波器二阶数学模型中的标称函数部分,v0,v1是有限时间观测器的内部状态变量,
Figure BDA0002550277210000037
是有限时间干扰观测器估计的数学模型状态变量x2
Figure BDA0002550277210000041
是有限时间干扰观测器估计的未知扰动g(t),
Figure BDA0002550277210000042
是有限时间干扰观测器估计的未知扰动的一阶导数
Figure BDA00025502772100000413
有限时间干扰观测器的增益λ012,K是大于零的可调增益;
有限时间干扰观测器的观测误差定义为:
Figure BDA0002550277210000044
其中,L是未知扰动的二阶导数的上界,即
Figure BDA0002550277210000045
优选地,所述步骤S3的具体步骤如下:
S3-1:定义跟踪误差为e=x-r,跟踪误差的一阶导数为
Figure BDA0002550277210000046
其中r为参考电流信号;
S3-2:定义快速非奇异终端滑模面为
Figure BDA0002550277210000047
其中,可调滑模参数满足k1>0,k2>0,1<α1<2,α2>α1;sign(e)是关于跟踪误差的符号函数,
Figure BDA0002550277210000048
是关于跟踪误差一阶导数的符号函数;
因此,快速非奇异终端滑模面的一阶导数为:
Figure BDA0002550277210000049
S3-3,令
Figure BDA00025502772100000410
可得等效控制项ueq为:
Figure BDA00025502772100000411
S3-4,由于未知扰动g(t)的存在,定义一个扰动补偿项ud为:
Figure BDA00025502772100000412
其中,
Figure BDA0002550277210000051
是未知扰动g(t)的估计值,
Figure BDA0002550277210000052
由有限时间干扰观测器的输出
Figure BDA0002550277210000053
进行估计;
S3-5,为了保证控制系统的鲁棒性,定义切换控制项usw为:
Figure BDA0002550277210000054
其中,Kw>0是保证李雅普诺夫函数为半正定的任意可调参数,sign(s)是关于快速非奇异终端滑模面的符号函数;
S3-6,所述理想的快速非奇异终端滑模控制律由等效控制项,扰动补偿项和切换控制项构成:
Figure BDA0002550277210000055
优选地,所述步骤S4的具体步骤如下:
S4-1:构造双隐层递归神经网络,所述双隐层递归神经网络的结构包含4层结构,分别为输入层、第一隐含层、第二隐含层和输出层,同时输出层的结果将反馈给输入层;
输入层:输入层的第j个节点的输出θj表示为:
θj=inj·Wrj·exYj=1,2,...,m (11);
其中,inj为双隐层递归神经网络的第j个输入,exY为上一时刻神经网络的输出值,Wrj为第j个输入层节点的反馈权值,反馈权值向量定义为Wr=[Wr1 Wr2 … Wrj];
第一隐含层:第一隐含层的第j个节点输出结果φ1j为:
Figure BDA0002550277210000056
其中,第一隐含层输出向量为Φ1=[φ11 φ12 … φ1j],且φ1j表示第一隐含层第j个节点的输出,第一隐含层的高斯函数中心向量为C1=[c11,c12,…,c1n]T∈Rn×1,高斯函数基宽向量为B1=[b11,b12,…,b1n]T∈Rn×1,且c1n是第一隐含层的第n个节点中心向量,且b1n是第一隐含层的第n个节点中心向量,Rn×1表示实数域内n行1列的向量;
第二隐含层:第二隐含层第k个节点输出结果φ2k为:
Figure BDA0002550277210000061
其中,第二隐含层输出向量为Φ2=[φ21 φ22 ... φ2k],且φ2k表示第二隐含层第k个节点的输出,第二隐含层的高斯函数中心向量为C2=[c21 c22...c2l]T∈Rl×1,高斯函数基宽向量为B2=[b21 b22...b2l]T∈Rl×1,且c2l是第二隐含层的第l个节点中心向量,b2l是第二隐含层的第l个节点中心向量,Rl×1表示实数域内l行1列的向量;
输出层:双隐层递归神经网络的输出结果为:
Y=WTΦ2=W1φ21+W2φ22+...+Wlφ2l (14);
其中,W=[W1 W2 ... Wl]是双隐层递归神经网络的输出权值向量,Wl表示第二隐含层的第l个节点与输出值之间的权值向量;
S4-2:根据最优逼近理论,存在最优参数
Figure BDA0002550277210000062
使得
Figure BDA0002550277210000063
其中ε为最优逼近误差,用双隐层递归神经网络的输出(14)代替公式(9)中切换增益Kw以表示估计的切换增益
Figure BDA0002550277210000064
Figure BDA0002550277210000065
最终估计的切换控制项定义为
Figure BDA0002550277210000066
则基于有限时间扰动观测器的自适应滑模控制律为:
Figure BDA0002550277210000071
S4-3:当切换增益被双隐层递归神经网络逼近以后,逼近误差定义为:
Figure BDA0002550277210000072
其中
Figure BDA0002550277210000073
为了求取网络参数的自适应律,在
Figure BDA0002550277210000074
处对
Figure BDA0002550277210000075
进行泰勒展开可得:
Figure BDA0002550277210000076
其中,
Figure BDA0002550277210000077
是第一隐含层中心向量的估计误差,
Figure BDA0002550277210000078
是第一隐含层基宽向量的估计误差,
Figure BDA0002550277210000079
是第二隐含层中心向量的估计误差,
Figure BDA00025502772100000710
是第二隐含层基宽向量的估计误差,
Figure BDA00025502772100000711
是反馈权值的估计误差,Oh为泰勒展开的高阶项,
Figure BDA00025502772100000712
Figure BDA00025502772100000713
是第二隐含层的输出向量Φ2分别对最优参数
Figure BDA00025502772100000714
Figure BDA00025502772100000715
的导数;
S4-4,为了保证采用公式(15)所示的控制律的单相有源电力滤波器系统是稳定的,考虑如下自适应律:
Figure BDA00025502772100000716
Figure BDA00025502772100000717
Figure BDA00025502772100000718
Figure BDA0002550277210000081
Figure BDA0002550277210000082
Figure BDA0002550277210000083
其中,
Figure BDA0002550277210000084
分别为双隐层递归神经网络的权值、反馈权值、第一隐含层中心、第一隐含层基宽、第二隐含层中心和第二隐含层基宽参数的逼近误差的一阶导数,η123456是可调正常数。
优选地,为了证明控制系统稳定性,定义李雅普诺夫函数为:
Figure BDA0002550277210000085
对李雅普诺夫函数求一阶导数,并将等式(6)和自适应滑模控制律(15)代入李雅普诺夫函数一阶求导后的公式可得:
Figure BDA0002550277210000086
其中,
Figure BDA0002550277210000087
是定义的中间变量,无特殊含义;
进一步简化李雅普诺夫函数的一阶导数(25),将等式(16)—(23)代入公式(25)可得:
Figure BDA0002550277210000091
其中,
Figure BDA0002550277210000092
是神经网络逼近的高阶误差上界,
Figure BDA0002550277210000093
是观测器误差的上界,实际切换增益K'w与估计的切换增益
Figure BDA0002550277210000094
误差定义为
Figure BDA0002550277210000095
并且满足|σ|≤Γ,根据李雅普诺夫稳定性理论,只要保证
Figure BDA0002550277210000096
即可保证李雅普诺夫的一阶导数是半正定的,即
Figure BDA0002550277210000097
有益效果:本发明提供一种基于干扰观测的有源电力滤波器自适应滑模控制,实现简单,能提高动态补偿性能,并有效降低电网电流的谐波总畸变率,不仅能克服系统未知扰动带来的负面影响,同时能保证系统的跟踪误差在有限时间内收敛,此外,本发明采用双隐层递归神经网络对切换增益进行逼近,保证切换增益达到最优值,经过李雅普诺夫稳定性证明设计的自适应滑模控制律在单相有源电力滤波器中是稳定可行的,因此在实际有源电力滤波器中拥有较高的实用性。
附图说明
图1为本发明方法的基于干扰观测器的自适应滑模控制器原理图;
图2为本发明的单相有源电力滤波器的结构图;
图3为本发明的三相并联电源型有源电力滤波器的结构图;
图4为本发明的双隐层递归神经网络结构图;
图5为本发明的电网电流曲线图;
图6为本发明的谐波电流跟踪曲线图;
图7为本发明的跟踪误差图;
图8为本发明的双隐层递归神经网络逼近的切换增益曲线;
图9为本发明的稳态补偿下电网电流畸变率图;
图10为本发明的动态补偿下电网电流畸变率图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
图1为本发明的涉及原理图,图1表明利用传感器从负载电流中检测出谐波电流作为参考信号r。首先设计一种快速非奇异终端滑模面,并根据滑模面可以求取等效控制项ueq;其次考虑系统的未知扰动,利用状态变量x和控制律
Figure BDA0002550277210000101
等参数构造扰动补偿项ud;最后利用双隐层递归神经网络逼近切换增益,得到逼近后的切换控制项
Figure BDA0002550277210000102
最终所设计的控制律由切换控制、等效控制和扰动补偿项相加所得,自适应控制律经过有源电力滤波器系统后输出一个谐波补偿信号x,利用负反馈使系统的误差趋于零,最终实现谐波电流快速、无静差跟踪参考谐波电流的目的。
本发明的一种基于干扰观测的有源电力滤波器自适应滑模控制,包括以下步骤:
S1,建立单相有源电力滤波器数学模型
如图2所示,为单相有源电力滤波器的结构图,其中,Us是电网电压,iL是负载电流,is是电网电流,ic是补偿谐波电流,
Figure BDA0002550277210000111
是参考谐波电流,L是交流侧线路总电感,R是交流侧线路总电阻,Qi(i=1,2,3,4)是IGBT电力电子开关器件,Udc是直流侧电压。在图2中一共包含4个电力电子开关器件,根据电路理论,开关管一共两种情况,因此开关函数H具体定义为:
Figure BDA0002550277210000112
基于上述开关函数和基尔霍夫定理,建立单相有源电力滤波器动力学方程如下所示:
Figure BDA0002550277210000113
其中,UMN=Udc*H为有源电力滤波器交流侧电压,对上式简化后得到:
Figure BDA0002550277210000114
定义状态变量x1=i,x2=di/dt,其中i表示谐波补偿电流,在考虑不确定扰动g(t)是有界的情况下,在考虑不确定性扰动g(t)有界的情况下,单相有源电力滤波器二阶数学模型如下所示:
Figure BDA0002550277210000115
其中,
Figure BDA0002550277210000116
是已知标称函数,在标称函数中,L是交流侧线路总标称电感,R是交流侧线路总标称电阻,Udc是直流侧稳定的电压,Us是电网电压,控制变量为u=H,不确定扰动的界限满足
Figure BDA0002550277210000117
S2、定义有限时间干扰观测器的微分方程如下所示:
Figure BDA0002550277210000121
其中,Ψ=f(x)+bu是单相有源电力滤波器二阶数学模型中的标称函数部分,v0,v1是有限时间观测器的内部状态变量,
Figure BDA0002550277210000122
是有限时间干扰观测器估计的数学模型状态变量x2
Figure BDA0002550277210000123
是有限时间干扰观测器估计的未知扰动g(t),
Figure BDA0002550277210000124
是有限时间干扰观测器估计的未知扰动的一阶导数
Figure BDA0002550277210000129
有限时间干扰观测器的增益λ012,K是大于零的可调增益;
有限时间干扰观测器的观测误差定义为:
Figure BDA0002550277210000126
其中,L是未知扰动的二阶导数的上界,即
Figure BDA0002550277210000127
S3,定义跟踪误差及其一阶导数,利用跟踪误差和跟踪误差的导数来定义快速非奇异终端滑模面,根据滑模面求取得到等效控制项,然后定义切换控制项并加上扰动补偿项,将等效控制项、切换控制项和扰动补偿项相加得到的理想滑模控制律,具体步骤为:
S3-1:定义跟踪误差为e=x-r,跟踪误差的一阶导数为
Figure BDA0002550277210000128
其中r为参考电流信号;
S3-2:定义快速非奇异终端滑模面为
Figure BDA0002550277210000131
其中,可调滑模参数满足k1>0,k2>0,1<α1<2,α2>α1;sign(e)是关于误差的符号函数,
Figure BDA0002550277210000132
是关于误差一阶导数的符号函数;
因此,快速非奇异终端滑模面的一阶导数为:
Figure BDA0002550277210000133
S3-3,令
Figure BDA0002550277210000134
可得等效控制项ueq为:
Figure BDA0002550277210000135
S3-4,由于未知扰动g(t)的存在,定义一个扰动补偿项ud为:
Figure BDA0002550277210000136
其中,
Figure BDA0002550277210000137
是未知扰动g(t)的估计值,
Figure BDA0002550277210000138
由有限时间干扰观测器的输出
Figure BDA0002550277210000139
进行估计;
S3-5,为了保证控制系统的鲁棒性,定义切换控制项usw为:
Figure BDA00025502772100001310
其中,切换增益Kw>0是保证李雅普诺夫函数为半正定的任意可调参数,sign(s)是关于快速非奇异终端滑模面的符号函数;
S3-6,所述理想的快速非奇异终端滑模控制律由等效控制项,扰动补偿项和切换控制项构成:
Figure BDA00025502772100001311
S4,步骤S3中设计的理想滑模控制律能够保证单相有源电力滤波器系统的稳定,但上述滑模控制器的切换增益Kw选取是困难,因此本发明利用双隐层递归神经网络对切换增益进行逼近,使得增益能够达到最优值,具体包括以下步骤:
S4-1,构造双隐层递归神经网络,如图4所示,双隐层递归神经网络的结构包含4层结构,分别为输入层、第一隐含层、第二隐含层和输出层,同时输出层的结果将反馈给输入层;
输入层:输入层的第j个节点的输出θi表示为:
θj=inj·Wrj·exY,j=1,2,...,m (11);
其中,inj为双隐层递归神经网络的第j个输入,exY为上一时刻神经网络的输出值,Wrj为第j个输入层节点的反馈权值,反馈权值向量定义为Wr=[Wr1 Wr2 … Wrj];
第一隐含层:第一隐含层的第j个节点输出结果φ1j为:
Figure BDA0002550277210000141
其中,第一隐含层输出向量为Φ1=[φ11 φ12 ... φ1j],且φ1j表示第一隐含层第j个节点的输出,第一隐含层的高斯函数中心向量为C1=[c11,c12,…,c1n]T∈Rn×1,高斯函数基宽向量为B1=[b11,b12,…,b1n]T∈Rn×1,且c1n是第一隐含层的第n个节点中心向量,且b1n是第一隐含层的第n个节点中心向量,Rn×1表示实数域内n行1列的向量;
第二隐含层:第二隐含层第k个节点输出结果φ2k为:
Figure BDA0002550277210000142
其中,第二隐含层输出向量为Φ2=[φ21 φ22 ... φ2k],且φ2k表示第二隐含层第k个节点的输出,第二隐含层的高斯函数中心向量为C2=[c21 c22...c2l]T∈Rl×1,高斯函数基宽向量为B2=[b21 b22...b2l]T∈Rl×1,且c2l是第二隐含层的第l个节点中心向量,b2l是第二隐含层的第l个节点中心向量,Rl×1表示实数域内l行1列的向量;
输出层:双隐层递归神经网络的输出结果为:
Y=WTΦ2=W1φ21+W2φ22+...+Wlφ2l (14);
其中,W=[W1 W2 ... Wl]是双隐层递归神经网络的输出权值向量,Wl表示第二隐含层的第l个节点与输出值之间的权值向量;
S4-2,根据最优逼近理论,存在最优参数
Figure BDA0002550277210000151
W*,使得
Figure BDA0002550277210000152
其中ε为最优逼近误差,用双隐层递归神经网络的输出(14)代替公式(9)中切换增益Kw以表示估计的切换增益
Figure BDA0002550277210000153
Figure BDA0002550277210000154
最终的切换控制项定义为
Figure BDA0002550277210000155
则基于有限时间扰动观测器的自适应滑模控制律为:
Figure BDA0002550277210000156
S4-3:当切换增益被双隐层递归神经网络逼近以后,逼近误差定义为:
Figure BDA0002550277210000157
其中
Figure BDA0002550277210000158
为了求取网络参数的自适应律,在
Figure BDA0002550277210000161
处对
Figure BDA0002550277210000162
进行泰勒展开可得:
Figure BDA0002550277210000163
其中,
Figure BDA0002550277210000164
是第一隐含层中心向量的估计误差,
Figure BDA0002550277210000165
是第一隐含层基宽向量的估计误差,
Figure BDA0002550277210000166
是第二隐含层中心向量的估计误差,
Figure BDA0002550277210000167
是第二隐含层基宽向量的估计误差,
Figure BDA0002550277210000168
是反馈权值的估计误差,Oh为泰勒展开的高阶项,
Figure BDA0002550277210000169
Figure BDA00025502772100001610
是第二隐含层的输出向量Φ2分别对最优参数
Figure BDA00025502772100001611
Figure BDA00025502772100001612
的导数;
S4-4,为了保证采用公式(15)所示的控制律的单相有源电力滤波器系统是稳定的,考虑如下自适应律:
Figure BDA00025502772100001613
Figure BDA00025502772100001614
Figure BDA00025502772100001615
Figure BDA00025502772100001616
Figure BDA00025502772100001617
Figure BDA00025502772100001618
其中,
Figure BDA00025502772100001619
分别为双隐层递归神经网络的权值、反馈权值、第一隐含层中心、第一隐含层基宽、第二隐含层中心和第二隐含层基宽参数的逼近误差的一阶导数,η123456是可调正常数。
当前述步骤完成以后,为了证明控制系统稳定性,定义如下李雅普诺夫函数来证明系统的稳定性:
Figure BDA0002550277210000171
对李雅普诺夫函数求一阶导数,并将等式(6)和自适应滑模控制律(15)代入李雅普诺夫函数一阶求导后的公式可得:
Figure BDA0002550277210000172
其中,
Figure BDA0002550277210000173
是定义的中间变量,无特殊含义;
进一步简化李雅普诺夫函数的一阶导数(25),将等式(16)—(23)代入公式(25)可得:
Figure BDA0002550277210000174
其中,
Figure BDA0002550277210000175
是神经网络逼近的高阶误差上界,
Figure BDA0002550277210000176
是观测器误差的上界,实际切换增益与估计的切换增益误差定义为
Figure BDA0002550277210000177
并且满足σ≤Γ,同时不需要知道上界的值,根据李雅普诺夫稳定性理论,只要保证
Figure BDA0002550277210000178
即可保证李雅普诺夫的一阶导数是半正定的,即
Figure BDA0002550277210000179
当稳定性证明完成之后,在MATLAB/Simulink中搭建单相有源电力滤波器模型,仿真过程中采用以下参数:
电网电压为Us=24V,f=50Hz;非线性稳态负载的电阻R1=5Ω,R2=15Ω,电容C=1000uF,动态非线性负载的电阻为R1=5Ω,R2=15Ω,电容C=1000uF,有源电力滤波器主电路参数包括线路电感为Lc=10mH,电阻为Rc=0.1Ω;直流侧电压采用传统的PI控制方法,参考电压被设定为50V。实验结果图如图5、图6、图7、图8、图9、图10所示。
为了方便观察,本发明在仿真时,在0.1s时将单相有源电力滤波器并入电网,即此时有源电力滤波器开始工作进行谐波电流补偿,同时于0.4s接入一个非线性负载来模拟动态仿真情况。
图5是本发明的电网电流曲线图,可以看出,在0.1s之前,电源线路是严重畸变的,而当有源电力滤波器并网以后,电源电流很快趋于正弦波形,当负载变化时,电源电流依然是正弦波形。
图6是本发明的谐波电流跟踪曲线图,可以看到当单相有源电力滤波器开始工作以后,补偿谐波电流在很短的时间内就跟踪上了参考谐波电流,即使在0.4s负载变化时,依然能很快跟踪上。
为了详细比较跟踪性能,本发明绘制了如图7所示的跟踪误差图,从跟踪误差图可以看出,跟踪误差较为平滑,并且误差范围很小。
图8是本发明的双隐层递归神经网络逼近的切换增益曲线,可以看到,从仿真开始,增益就处于递增状态,在0.1s之后逼近值就趋于稳定。
为了直观反映本发明的有益效果,测量了0.3s时的总谐波畸变率如图9所示,此时的总谐波畸变率为4.24%,当负载在0.4s变化时,测量0.6s的总谐波畸变率如图10所示,此时的值是2.75%。这说明了本发明提出的基于干扰观测器的自适应滑模控制律拥有很强的稳态还是动态补偿性能,并且补偿后的电流畸变率能达到很低。
上述实施例主要针对单相有源电力滤波器设计,但实际上设计的控制器不仅适用于单相有源电力滤波器,针对如图3所示的三相三线制有源电力滤波器数学模型同样适用,为了进行说明,类似单相模型,利用电压和电流定理,建立如下三相动力学方程:
Figure BDA0002550277210000191
其中,i=[i1 i2 i3]T为补偿电流向量,i1 i2 i3分别对应a、b、c三相的电流,dk=[d1k d2k d3k]T为开关状态函数向量;d1k d2k d3k分别对应a、b、c三相的开关状态函数。因此第n相的开关状态函数定义为
Figure BDA0002550277210000192
ck为第k相开关函数,它的定义类似单相,具体为:
Figure BDA0002550277210000193
定义新的状态向量x'1=[i1 i2 i3]T
Figure BDA0002550277210000194
三相有源电力滤波器的动力学方程可简化为如下二阶通用模型:
Figure BDA0002550277210000195
其中
Figure BDA0002550277210000196
为三维的列向量,控制变量为u=dk,可以看到单相通用模型(3)和三相通用模型(29)的区别仅在于一个是一维标量,一个是三维列向量。由于三相数学模型的每一相之间可以认为是没有耦合的,因此针对三相模型的设计可以化简为分别对每一相进行设计,其设计思路与单相有源电力滤波器一致。
本发明所讲述的具体实施例所保护的内容是控制算法的设计,该算法并没有特别指要求适用的有源电力滤波器结构,本具体实施例仅仅为了方便讲述采用单相有源电力滤波器结构进行说明。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (2)

1.一种基于干扰观测的有源电力滤波器自适应滑模控制,其特征在于,包括以下步骤:
S1,定义开关函数H来表示有源电力滤波器物理模型中开关管的导通情况,根据开关函数H的定义以及基尔霍夫定理建立单相有源电力滤波器动力学方程,然后根据动力学方程的状态变量i为x,得到x的二阶导数
Figure FDA0003728190570000016
的通用数学模型;
所述步骤S1的具体步骤如下:
S1-1:单相有源电力滤波器物理模型包含四个电力电子开关管,根据电路理论,开关管一共两种情况,因此开关函数H具体定义为:
Figure FDA0003728190570000011
S1-2:根据定义的开关函数H以及基尔霍夫定理,建立单相有源电力滤波器动力学方程如下所示:
Figure FDA0003728190570000012
其中,i表示谐波补偿电流,UMN=Udc*H为有源电力滤波器交流侧电压;
S1-3:定义状态变量x1=i,x2=di/dt,在考虑和时间相关的不确定性扰动函数g(t)有界的情况下,单相有源电力滤波器二阶数学模型如下所示:
Figure FDA0003728190570000013
其中,
Figure FDA0003728190570000014
是已知标称函数;在标称函数中,L是交流侧线路总电感,R是交流侧线路总电阻,Udc是直流侧电压,Us是电网电压,控制变量为u=H,不确定扰动界限满足
Figure FDA0003728190570000015
S2,为了补偿系统中的未知扰动影响,定义有限时间干扰观测器的微分方程;
所述步骤S2中定义的有限时间干扰观测器的微分方程如下所示:
Figure FDA0003728190570000021
其中,Ψ=f(x)+bu是单相有源电力滤波器二阶数学模型中的标称函数部分,v0,v1是有限时间观测器的内部状态变量,
Figure FDA0003728190570000022
是有限时间干扰观测器估计的数学模型状态变量x2
Figure FDA0003728190570000023
是有限时间干扰观测器估计的未知扰动g(t),
Figure FDA0003728190570000024
是有限时间干扰观测器估计的未知扰动的一阶导数
Figure FDA0003728190570000025
有限时间干扰观测器的增益λ012,K是大于零的可调增益;
有限时间干扰观测器的观测误差定义为:
Figure FDA0003728190570000026
其中,L是未知扰动的二阶导数的上界,即
Figure FDA0003728190570000027
S3,定义跟踪误差及其一阶导数,然后定义快速非奇异终端滑模面,通过滑模面求取等效控制项,通过干扰观测器微分方程求取扰动补偿项,并考虑切换控制项,将等效控制项、切换控制项和未知扰动补偿项相加得到理想的快速非奇异终端滑模控制律;
所述步骤S3的具体步骤如下:
S3-1:定义跟踪误差为e=x-r,跟踪误差的一阶导数为
Figure FDA0003728190570000031
其中r为参考电流信号;
S3-2:定义快速非奇异终端滑模面为
Figure FDA0003728190570000032
其中,可调滑模参数满足k1>0,k2>0,1<α1<2,α2>α1;sign(e)是关于跟踪误差的符号函数,
Figure FDA0003728190570000033
是关于跟踪误差一阶导数的符号函数;
因此,快速非奇异终端滑模面的一阶导数为:
Figure FDA0003728190570000034
S3-3,令
Figure FDA0003728190570000035
可得等效控制项ueq为:
Figure FDA0003728190570000036
S3-4,由于未知扰动g(t)的存在,定义一个扰动补偿项ud为:
Figure FDA0003728190570000037
其中,
Figure FDA0003728190570000038
是未知扰动g(t)的估计值,
Figure FDA0003728190570000039
由有限时间干扰观测器的输出
Figure FDA00037281905700000310
进行估计;
S3-5,为了保证控制系统的鲁棒性,定义切换控制项usw为:
Figure FDA00037281905700000311
其中,Kw>0是保证李雅普诺夫函数为半正定的任意可调参数,sign(s)是关于快速非奇异终端滑模面的符号函数;
S3-6,所述理想的快速非奇异终端滑模控制律由等效控制项,扰动补偿项和切换控制项构成:
Figure FDA0003728190570000041
S4,构造双隐层递归神经网络,通过网络输出逼近切换控制项的切换增益,最终实现自适应快速非奇异终端滑模控制律,并进行稳定证明;
所述步骤S4的具体步骤如下:
S4-1:构造双隐层递归神经网络,所述双隐层递归神经网络的结构包含4层结构,分别为输入层、第一隐含层、第二隐含层和输出层,同时输出层的结果将反馈给输入层;
输入层:输入层的第j个节点的输出θj表示为:
θj=inj·Wrj·exY,j=1,2,...,m (11);
其中,inj为双隐层递归神经网络的第j个输入,exY为上一时刻神经网络的输出值,Wrj为第j个输入层节点的反馈权值,反馈权值向量定义为Wr=[Wr1 Wr2…Wrj];
第一隐含层:第一隐含层的第j个节点输出结果φ1j为:
Figure FDA0003728190570000042
其中,第一隐含层输出向量为Φ1=[φ11 φ12…φ1j],且φ1j表示第一隐含层第j个节点的输出,第一隐含层的高斯函数中心向量为C1=[c11,c12,…,c1n]T∈Rn×1,高斯函数基宽向量为B1=[b11,b12,…,b1n]T∈Rn×1,且c1n是第一隐含层的第n个节点中心向量,且b1n是第一隐含层的第n个节点中心向量,Rn×1表示实数域内n行1列的向量;
第二隐含层:第二隐含层第k个节点输出结果φ2k为:
Figure FDA0003728190570000051
其中,第二隐含层输出向量为Φ2=[φ21 φ22…φ2k],且φ2k表示第二隐含层第k个节点的输出,第二隐含层的高斯函数中心向量为C2=[c21 c22...c2l]T∈Rl×1,高斯函数基宽向量为B2=[b21 b22...b2l]T∈Rl×1,且c2l是第二隐含层的第l个节点中心向量,b2l是第二隐含层的第l个节点中心向量,Rl×1表示实数域内l行1列的向量;
输出层:双隐层递归神经网络的输出结果为:
Y=WTΦ2=W1φ21+W2φ22+...+Wlφ2l (14);
其中,W=[W1 W2...Wl]是双隐层递归神经网络的输出权值向量,Wl表示第二隐含层的第l个节点与输出值之间的权值向量;
S4-2:根据最优逼近理论,存在最优参数
Figure FDA0003728190570000052
使得
Figure FDA0003728190570000053
其中ε为最优逼近误差,用双隐层递归神经网络的输出(14)代替公式(9)中切换增益Kw以表示估计的切换增益
Figure FDA0003728190570000054
Figure FDA0003728190570000055
最终估计的切换控制项定义为
Figure FDA0003728190570000056
则基于有限时间扰动观测器的自适应滑模控制律为:
Figure FDA0003728190570000057
S4-3:当切换增益被双隐层递归神经网络逼近以后,逼近误差定义为:
Figure FDA0003728190570000061
其中
Figure FDA0003728190570000062
为了求取网络参数的自适应律,在
Figure FDA0003728190570000063
处对
Figure FDA0003728190570000064
进行泰勒展开可得:
Figure FDA0003728190570000065
其中,
Figure FDA0003728190570000066
是第一隐含层中心向量的估计误差,
Figure FDA0003728190570000067
是第一隐含层基宽向量的估计误差,
Figure FDA0003728190570000068
是第二隐含层中心向量的估计误差,
Figure FDA0003728190570000069
是第二隐含层基宽向量的估计误差,
Figure FDA00037281905700000610
是反馈权值的估计误差,Oh为泰勒展开的高阶项,
Figure FDA00037281905700000611
Figure FDA00037281905700000612
是第二隐含层的输出向量Φ2分别对最优参数
Figure FDA00037281905700000613
Figure FDA00037281905700000614
的导数;
S4-4,为了保证采用公式(15)所示的控制律的单相有源电力滤波器系统是稳定的,考虑如下自适应律:
Figure FDA00037281905700000615
Figure FDA00037281905700000616
Figure FDA00037281905700000617
Figure FDA00037281905700000618
Figure FDA00037281905700000619
Figure FDA00037281905700000620
其中,
Figure FDA0003728190570000071
分别为双隐层递归神经网络的权值、反馈权值、第一隐含层中心、第一隐含层基宽、第二隐含层中心和第二隐含层基宽参数的逼近误差的一阶导数,η123456是可调正常数。
2.根据权利要求1所述的一种基于干扰观测的有源电力滤波器自适应滑模控制,其特征在于:为了证明控制系统稳定性,定义李雅普诺夫函数为:
Figure FDA0003728190570000072
对李雅普诺夫函数求一阶导数,并将等式(6)和自适应滑模控制律(15)代入李雅普诺夫函数一阶求导后的公式可得:
Figure FDA0003728190570000073
其中,
Figure FDA0003728190570000074
是定义的中间变量,无特殊含义;
进一步简化李雅普诺夫函数的一阶导数(25),将等式(16)—(23)代入公式(25)可得:
Figure FDA0003728190570000081
其中,
Figure FDA0003728190570000082
是神经网络逼近的高阶误差上界,
Figure FDA0003728190570000083
是观测器误差的上界,实际切换增益K'w与估计的切换增益
Figure FDA0003728190570000084
误差定义为
Figure FDA0003728190570000085
并且满足|σ|≤Γ,根据李雅普诺夫稳定性理论,只要保证
Figure FDA0003728190570000086
即可保证李雅普诺夫的一阶导数是半正定的,即
Figure FDA0003728190570000087
CN202010573671.2A 2020-06-22 2020-06-22 一种基于干扰观测的有源电力滤波器自适应滑模控制 Active CN111799795B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010573671.2A CN111799795B (zh) 2020-06-22 2020-06-22 一种基于干扰观测的有源电力滤波器自适应滑模控制

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010573671.2A CN111799795B (zh) 2020-06-22 2020-06-22 一种基于干扰观测的有源电力滤波器自适应滑模控制

Publications (2)

Publication Number Publication Date
CN111799795A CN111799795A (zh) 2020-10-20
CN111799795B true CN111799795B (zh) 2022-08-19

Family

ID=72803710

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010573671.2A Active CN111799795B (zh) 2020-06-22 2020-06-22 一种基于干扰观测的有源电力滤波器自适应滑模控制

Country Status (1)

Country Link
CN (1) CN111799795B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112987569B (zh) * 2021-02-08 2021-11-02 南京工业大学 一种扰动上界自适应分数阶全局终端滑模无模型控制方法
CN112949180B (zh) * 2021-03-02 2022-09-06 南京工业大学 一种基于hgdob与rbf神经网络的盾构风电回转支承试验台液压加载控制方法
CN113568423B (zh) * 2021-08-01 2024-01-16 西北工业大学 一种考虑电机故障的四旋翼无人机智能容错控制方法
CN114371623B (zh) * 2022-01-10 2022-12-23 南京工业大学 一种结构热试验神经网络观测器分数阶全局滑模控制方法
CN114861435B (zh) * 2022-04-29 2024-06-11 南京航空航天大学 有限时间观测器下的多智能体系统一致性滑模控制算法
CN115818490B (zh) * 2022-11-22 2023-07-21 哈尔滨工程大学 一种rov释放回收的半主动升沉补偿装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108923430A (zh) * 2018-07-16 2018-11-30 河海大学常州校区 有源电力滤波器神经网络全局快速终端滑模控制方法及计算设备
CN109921422A (zh) * 2018-08-13 2019-06-21 河海大学常州校区 基于单反馈神经网络的有源电力滤波器非奇异终端滑模控制方法
CN110277790A (zh) * 2019-07-17 2019-09-24 河海大学常州校区 一种有源电力滤波器终端滑模控制方法
CN110286594A (zh) * 2019-07-17 2019-09-27 河海大学常州校区 一种有源电力滤波器自适应动态终端滑模控制方法
CN110350546A (zh) * 2019-07-11 2019-10-18 河海大学常州校区 一种单相有源电力滤波器控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108923430A (zh) * 2018-07-16 2018-11-30 河海大学常州校区 有源电力滤波器神经网络全局快速终端滑模控制方法及计算设备
CN109921422A (zh) * 2018-08-13 2019-06-21 河海大学常州校区 基于单反馈神经网络的有源电力滤波器非奇异终端滑模控制方法
CN110350546A (zh) * 2019-07-11 2019-10-18 河海大学常州校区 一种单相有源电力滤波器控制方法
CN110277790A (zh) * 2019-07-17 2019-09-24 河海大学常州校区 一种有源电力滤波器终端滑模控制方法
CN110286594A (zh) * 2019-07-17 2019-09-27 河海大学常州校区 一种有源电力滤波器自适应动态终端滑模控制方法

Also Published As

Publication number Publication date
CN111799795A (zh) 2020-10-20

Similar Documents

Publication Publication Date Title
CN111799795B (zh) 一种基于干扰观测的有源电力滤波器自适应滑模控制
CN111781838B (zh) 一种基于观测器的有源电力滤波器自适应动态滑模方法
CN107612409B (zh) 带偏磁控制的精简矩阵变换器模型预测控制方法
CN109950922B (zh) 一种适用于vsc-hvdc的多步模型预测控制方法
CN106230257A (zh) 一种双向直流变换器反馈线性化反步滑模控制方法
CN106374528A (zh) 一种基于改进下垂控制的微电网分散型控制策略模型
CN105978373B (zh) 实现微电网稳定的三相逆变器逆推滑模控制方法及系统
CN108336758A (zh) 一种基于纹波关联法的光伏组件mppt算法
CN110137971B (zh) 一种三相交流电力弹簧的电压平稳控制方法
CN104993494B (zh) 一种基于四象限电力电子变流器的电机模拟装置及方法
Rong et al. Output feedback control of single-phase UPQC based on a novel model
CN107546767A (zh) 一种光伏并网逆变器的控制结构及控制方法
CN111614282A (zh) 矢量扇区定位方法、局部寻优模型预测控制方法及装置
CN105141164A (zh) 单相光伏并网逆变器的反演全局快速终端滑模控制方法
CN103887798B (zh) 有源电力滤波器的反演全局快速终端滑模控制方法
CN104868773A (zh) 基于Lyapunov状态函数的单相并网逆变器控制装置
CN110138241B (zh) 一种基于雅可比理论的单相电压源的稳定性分析方法
CN110429835A (zh) 一种基于lcl滤波的rbfnn分段在线优化无源控制系统及方法
CN105305448A (zh) 基于模糊pi复合控制的有源滤波器自抗扰控制方法
CN113411002B (zh) 一种基于极限学习机的滑模变结构的单相逆变器控制系统及方法
CN110277790B (zh) 一种有源电力滤波器终端滑模控制方法
CN109921422A (zh) 基于单反馈神经网络的有源电力滤波器非奇异终端滑模控制方法
CN110518625B (zh) 一种可变学习率bp-pid控制的并网逆变器直流分量抑制方法
CN110286594B (zh) 一种有源电力滤波器自适应动态终端滑模控制方法
CN103293965B (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