CN114280944A - 一种具有输出约束的pmsm系统有限时间动态面控制方法 - Google Patents

一种具有输出约束的pmsm系统有限时间动态面控制方法 Download PDF

Info

Publication number
CN114280944A
CN114280944A CN202111678368.XA CN202111678368A CN114280944A CN 114280944 A CN114280944 A CN 114280944A CN 202111678368 A CN202111678368 A CN 202111678368A CN 114280944 A CN114280944 A CN 114280944A
Authority
CN
China
Prior art keywords
function
time
following
yields
substituting
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
Application number
CN202111678368.XA
Other languages
English (en)
Other versions
CN114280944B (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.)
Guizhou University
Original Assignee
Guizhou 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 Guizhou University filed Critical Guizhou University
Priority to CN202111678368.XA priority Critical patent/CN114280944B/zh
Publication of CN114280944A publication Critical patent/CN114280944A/zh
Application granted granted Critical
Publication of CN114280944B publication Critical patent/CN114280944B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Feedback Control In General (AREA)

Abstract

本发明公开了一种具有输出约束的PMSM系统有限时间动态面控制方法,引入合适的Lyapunov‑Krasovskii泛函来处理时滞问题,同时考虑非线性变换函数将输出约束问题转化为无约束问题。然后,在有限时间反演(backstepping)框架下,通过应用神经网络估计未知非线性函数和引入一阶滤波器来解决“复杂性爆炸”问题,设计了一种神经自适应有限时间动态面控制方法。此外,还证明了该系统的所有信号都是有限时间稳定的,并且在不违反输出约束的情况下,跟踪误差在有限时间内缩小到原点的一个小邻域处。

Description

一种具有输出约束的PMSM系统有限时间动态面控制方法
技术领域
本发明属于具有时滞和非对称时变输出约束的PMSM系统控制技术领域,涉及一种具有输出约束的PMSM系统有限时间动态面控制方法。
背景技术
由于永磁同步电机(PMSM)系统相对于步进电机和直流电机具有运行效率高、体积小、功率密度高等显著优点,因此在包括机器人、车辆和飞机在内的汽车工业中被广泛采用。然而,永磁同步电机系统是高度非线性的多变量耦合系统,对复杂的系统参数摄动和不断变化的外部干扰非常敏感。因此,迫切需要设计一种优秀的控制解决方案,有效地应对永磁同步电机系统所面临的上述挑战。在过去几十年中,永磁同步电机系统采用了多种控制方案,包括滑模控制(SMC)、自适应控制和反演控制。
在上述方法中,由于反演算法易于与模糊逻辑系统或神经网络等智能估计器集成,因此它在构造鲁棒或自适应控制器以在一定程度上克服不确定非线性系统中产生的参数摄动和外部干扰方面显示出其潜力。因此,在永磁同步电机系统中进行了大量的自适应反演控制研究。然而,传统反演控制导致的“复杂性爆炸”问题妨碍了前面的控制结果在实践中的广泛应用。为了避免这个问题,在文献“S.Ling,H.Wang,and P.X.Liu,“AdaptiveFuzzy Dynamic Surface Control of Flexible-Joint Robot Systems With InputSaturation,”IEEE-CAAJ.Autom.Sin.,vol.6,no.1,pp.97-106,2019”(Ling Song,WangHuanqing,Liu Peter X.,“具有输入饱和度的柔性关节机器人系统的自适应模糊动态表面控制,IEEE-CAA Journal of Automatica Sinica,卷:6,期:1,页码:97-106,2019)”中提出了一种动态表面控制解决方案,其中一阶滤波器在设计的每一步近似虚拟控制器的导数。随后,也有人通过引入跟踪微分器、一阶滤波器和指令滤波器等合适的估计器,在永磁同步电机系统中研究了许多改进的反演控制方法。遗憾的是,尽管上述控制方法可以在一定程度上提高系统性能,但它们没有涉及到能够有效提高永磁同步电机系统稳态和瞬态性能的另外两个研究前沿,即有限时间控制和约束控制。
从时间优化控制系统的角度来看,与渐近稳定方式相比,有限时间控制方案在原点附近表现出更好的收敛性和更强的鲁棒性。鉴于此,在过去的几年中,为提高永磁同步电机系统的响应和鲁棒性,已经进行了许多工作。例如,文献“L.Fang,L.Ma,S.Ding,andD.Zhao,“Robust finite-time stabilization of a class of high-order stochasticnonlinear systems subject to output constraint and disturbances,”Int.J.RobustNonlinear Control,vol.29,no.16,pp.5550-5573,2019(Fang Liand,Ma Li,Ding,ZhaoDean,“受输出约束和干扰影响的一类高阶随机非线性系统的鲁棒性有限时间稳定化,”International Journal of Robust and Nonlinear Control,卷:29,期:16,页码:5550-5573,2019)”中的工作通过分数功率积分器技术研究了永磁同步电机系统的有限时间稳定问题。在文献“J.Jiang,X.Zhou,W.Zhao,W.Li,and W.Zhang,“A fast integral slidingmode controller with an extended state observer for position control ofpermanent magnet synchronous motor servo systems,”Front.Inf.Technol.Electron.Eng.,vol.21,no.8,pp.1239-1250,2020(Jiang Jun-feng,Zhou Xiao-jun,Zhao Wei,Li Wei,Zhang Wen-dong,用于永磁同步电机伺服系统位置控制的带扩展状态观测器的快速积分滑动模式控制器,Frontiers of InformationTechnology&Electronic Engineering,卷:21,期:8,页码:1239-1250,2020)”中,永磁同步电机系统的有限时间控制器被构造成在有限时间内跟踪所需的输出信号,从而产生了优异的性能。文献“M.Sharma and I.Kar,“Finite Time Disturbance Observer BasedGeometric Control of Quadrotors,”.6thConference onAdvances inControlandOptimizationofDynamical Systems(ACODS),2020,pp.295-300(Sharma Manmohan,KarIndrani,基于有限时间扰动观测器的四旋翼飞机的几何控制(ACODS),6th Conference onAdvances in Control and Optimization of Dynamical Systems,卷:53,期:1,页码:295-300,2020)”中的研究将有限时间观测器考虑到反演设计中,以解决永磁同步电机系统的不确定性和动态问题。上述结果表明,考虑有限时间工具可以更快地实现跟踪误差,这是本文构建永磁同步电机系统有限时间控制方案的动机之一。
在研究约束控制的前沿,已经提出了许多重要的方法,如李雅普诺夫函数(BLF)、预测控制和规定的性能控制。其中,基于BLF的方法被视为解决状态和输出约束的有效方式。在文献“Q.Zhu and Y.Liu,“Neural network adaptive finite-time control ofstochastic nonlinear systems with full state constraints,”AsianJ.Control,vol.23,no.4,pp.1728-1739,2021(Zhu Qidan,Liu Yongchao,“具有完全状态约束的随机非线性系统的神经网络自适应有限时间控制,”Asian Journal of Control,卷:23,期:4,页码:1728-1739,2021)”中,对数型BLF出现在有限时间稳定控制中,以将状态变量限制在一定范围内。文献“L.Fang,L.Ma,S.Ding,and J.H.Park,“Finite-Time Stabilization ofHigh-Order Stochastic Nonlinear Systems With Asymmetric Output Constraints,”IEEE Trans.Syst.Man Cybern.,vol.51,no.11,pp.7201-7213,2021(Fang Liandi,Ma Li,Ding Shihong,ParkJu H.,“具有不对称输出约束的高阶随机非线性系统的有限时间稳定化,”IEEE Transactions On Systems Man Cybernetics-Systems,卷:51,期:11,页码:7201-7213,2021)”中给出的自适应反演控制采用BLF来处理不对称输出约束。尽管之前基于BLF的工作有助于将永磁同步电机系统的输出和状态限制在某些特定区域,但它们无法摆脱使用分段BLF来简化控制器设计和稳定性综合。为了避免使用分段非对称BLF,在文献“D.Li,H.Han,and J.Qiao,“Deterministic Learning-Based Adaptive Neural Controlfor Nonlinear Full-State Constrained Systems.,”IEEE Trans.neural networksLearn.Syst.,2021,doi:10.1109/TNNLS.2021.3126320(Li Dapeng,Han Honggui,QiaoJunfei,“非线性全状态受限系统的基于确定学习的自适应神经控制,”IEEETransactions on Neural Networks and Learning Systems,doi:10.1109/TNNLS.2021.3126320,2021)”中创建了坐标变换函数,将状态约束系统转换为无约束系统。最近,在文献“Z.Lv,Y.Ma,J.Liu,and J.Yu,“Full-State Constrained Adaptive FuzzyFinite-Time Dynamic Surface Control for PMSM Drive Systems,”Int.J.FuzzySyst.,vol.23,no.3,pp.804-815,2021(Lv Zhenxiang,Ma Yumei,Liu Jiapeng,YuJinpeng,“PMSM驱动系统的全状态约束的自适应模糊有限时间动态表面控制,”International Journal of Fuzzy Systems,卷:23,期:3,页码:804-815,2021)”中的工作将精细转换应用于具有全状态约束的永磁同步电机系统的有限时间控制设计。即使如此,他们也忽略了在设计的控制器中处理存在的时滞问题,这有损控制永磁同步电机系统的鲁棒性和有效性。
为了解决非线性控制系统的时滞问题,常用的工具是构造合适的Lyapunov-Krasovskii泛函来处理控制论中的时变扰动。文献“Z.Zhang,S.Chen,and Y.Zheng,“FullyDistributed Scaled Consensus Tracking of High-Order Multiagent Systems WithTime Delays and Disturbances,”IEEE Trans.Ind.Informatics,vol.18,no.1,pp.305-314,2022(Zhang Zheng,Chen Shiming,Zheng Yuanshi,具有时间延迟和干扰的高阶多代理系统的全分布式规模化共识跟踪,IEEE Transactions on Industrial Informatics,卷:18,期:1,页码:305-314,2022)”中的工作介绍了Lyapunov-Krasovskii泛函,它解决了时滞多智能体控制系统的时滞干扰问题。文献“S.Li,L.Ding,H.Gao,Y.-J.Liu,N.Li,andZ.Deng,“Reinforcement Learning Neural Network-Based Adaptive Control forState and Input Time-Delayed Wheeled Mobile Robots,”IEEE Trans.Syst.ManCybern.,vol.50,no.11,pp.4171-4182,2020(Li Shu,Ding Liang,Gao Haibo,Liu Yan-Jun,Li Nan,Deng Zongquan,基于强化学习神经网络的状态和输入延时轮式移动机器人的自适应控制,IEEE Transactions onSystems Man Cybernetics-Systems,卷:50,期:11,页码:4171-4182,2020)”中的工作将Lyapunov-Krasovskii泛函与径向基函数NNs(RBFNNs)相结合,提出了一种时滞系统的自适应神经反演控制方案。注意,控制结果是解决时滞非线性控制系统时滞效应的有效方案,但永磁同步电机系统的有限时间控制稳定设计不能直接应用。然后一个实际的问题是:如何结合Lyapunov-Krasovskii泛函技术来解决具有时滞和非对称时变输出约束的永磁同步电机系统的有限时间跟踪问题。
发明内容
本发明要解决的技术问题是:提供一种具有输出约束的PMSM系统有限时间动态面控制方法,以解决现有技术中存在的技术问题。
本发明采取的技术方案为:一种具有输出约束的PMSM系统有限时间动态面控制方法,该方法包括以下步骤:
(1)定义变量x1=θ,x2=ω,x3=iq,x4=id,对(d-q)坐标系下永磁同步电机的动力学模型进
行简化,得到如下式:
Figure BDA0003453101860000031
式(2)受以下输出约束:
x1∈Πx1:={x1∈R:yd(t)-F11(t)<x1(t)<yd(t)+F12(t)} (3)
其中,时变函数F11(t)>0和F12(t)>0表示已知约束边界,x1(t)表示输出变量,Δfi(x(t-τi)),i=1,...,4是时滞项,x=(x1,x2,x3,x4)T∈R4是式(2)的整体状态,τi,i=1,...,4是时间常数,
Figure BDA0003453101860000032
a2=3np(Ld-Lq)/2,b1=-Rs/Lq,b2=-npLd/Lq
Figure BDA0003453101860000033
b4=1/Lq,c1=-Rs/Ld,c2=npLq/Ld,c3=1/Ld
设1:参考信号yd(t)及其n阶导数
Figure BDA0003453101860000034
(n=0,...,4)有界且连续;约束函数F11(t),F12(t)及其k阶导数
Figure BDA0003453101860000035
(k=0,...,4)是有界的和连续的;
引理1:一个连续函数
Figure BDA00034531018600000315
由f(0,...,0)=0给出,其中
Figure BDA00034531018600000316
(i=1,2,...,n,mi>0),有光滑正函数
Figure BDA00034531018600000317
满足ωi(0)=0,这样
Figure BDA0003453101860000036
由引理1可知,式(2)的时滞项Δfi(x(t-τi)),i=1,...,4用
Figure BDA0003453101860000037
表示,然后,根据杨氏不等式,有:
Figure BDA0003453101860000038
引理2:考虑
Figure BDA0003453101860000039
Figure BDA00034531018600000310
是常数,然后,对于实变量y和z,以下不等式成立:
Figure BDA00034531018600000311
引理3:考虑实数r=1,....,m,和xr∈R,r=1,....,m,得到:
Figure BDA00034531018600000312
定义1:给出非线性系统
Figure BDA00034531018600000313
其中f(λ)表示系统状态为λ∈Rn.的光滑函数,如果对于所有初始条件λ(t0)=λ0,有ρ>0和稳定时间T(ρ,λ0)<∞,使所有t≥t0+T,稳定时间为||λ(t)||<ρ,,则平衡点为λ=0的非线性系统
Figure BDA00034531018600000314
称为有限时间半全局实际稳定;
引理4:对于非线性系统
Figure BDA0003453101860000041
如果存在光滑正定函数V(λ)和标量a>0,
Figure BDA0003453101860000042
b>0,σ>0和
Figure BDA00034531018600000426
则:
Figure BDA0003453101860000043
那么非线性系统
Figure BDA0003453101860000044
在有限时间内是稳定的,其中稳定时间可近似为:
Figure BDA0003453101860000045
其中t≥T,存在:
Figure BDA0003453101860000046
引理5:对于每个变量,
Figure BDA0003453101860000047
存在φ1<φ2,φ12是奇数整数。一个不等式适用于以下情况:
Figure BDA0003453101860000048
其中
Figure BDA0003453101860000049
和γ2=(2φ-1-2(1+φ)(φ-1))/(1+φ)>0;
引理6:对于
Figure BDA00034531018600000420
有一个集合Λ,它由
Figure BDA00034531018600000421
给出,然后,对于
Figure BDA00034531018600000427
满足不等式
Figure BDA00034531018600000422
(2)使用RBFNNs技术,未知非线性函数在闭合集中以任意精度估计,因此,有:
Figure BDA00034531018600000410
其中Z=[z1,z2,…,zn]T是输入向量,
Figure BDA00034531018600000423
是期望值,RBFNN的权重向量,l>1是节点数,δ(Z)是满足|δ(Z)|<δM,的估计误差,δM表示未知的有界参数,
Figure BDA00034531018600000411
是基函数向量,其中选择
Figure BDA00034531018600000412
作为应用的高斯函数:
Figure BDA00034531018600000413
其中μi=[μi1,...,μim]和
Figure BDA00034531018600000425
分别为域中心和高斯函数的宽度;
让理想权重向量
Figure BDA00034531018600000424
为:
Figure BDA00034531018600000414
其中
Figure BDA00034531018600000415
表示更新的权重向量;
因此:
Figure BDA00034531018600000416
其中
Figure BDA00034531018600000417
θi和||·||分别是未知变量和·的2-范式;
(3)设计有限时间动态表面控制
A、非线性误差相关转换函数:引入一个非线性转换函数,将输出约束系统式(3)转换为非约束系统;
定义2:非线性转换函数可构造为:
Figure BDA00034531018600000418
其中,跟踪误差s1=x1-yd,S1表示转换误差,F11(t)>0和F12(t)>0表示平滑时变函数,假设F11(t),F12(t)满足以下关系:|F11(t)|<C0,|F12(t)|<C1,其中常数C0>0和C1>0。
从式(15)中看出,对于满足F11(0)<s1(0)<F12(0)的每个初始值,当S1有界为t∈[0,+∞)时,确保了s1的有界性和约束性,F11(t)和F12(t)将缩写为F11和F12
区分S1给出:
Figure BDA00034531018600000419
其中
Figure BDA0003453101860000051
Figure BDA0003453101860000052
利用(3)和(16),得到输出无约束子系统,如下式所示:
Figure BDA0003453101860000053
B、设计神经自适应有限时间控制器:
设误差坐标变换,如下所示:
v1=S1,v2=x22c (20)
v3=x33c,v4=x4,
其中αic,i=2,3表示以下一阶滤波器的输出:
Figure BDA0003453101860000054
其中εi表示时间常数,虚拟信号αi在后面设计;
类似地,过滤器误差yi定义为:
yi=αici,i=2,3 (22)
集成(2)和(19),在(20)中vi,i=1,...,4求导得:
Figure BDA0003453101860000055
引入估算误差
Figure BDA0003453101860000056
如下所示:
Figure BDA0003453101860000057
其中变量
Figure BDA0003453101860000058
代表θi的估计值;
神经自适应有限时间控制器的设计步骤如下:
步骤1:选择Lyapunov函数V1为:
Figure BDA0003453101860000059
Lyapunov-Krasovskii函数VM为:
Figure BDA00034531018600000510
其中已知
Figure BDA00034531018600000513
为常数;
取VM在(26)的时间导数,得:
Figure BDA00034531018600000511
其中,参数Γ>0,和正函数ωik以消除时间延迟;
然后,将(26)中的V1导数与(24)结合,得到:
Figure BDA00034531018600000512
将(23)合并到(28)会得出:
Figure BDA0003453101860000061
根据(4),得出:
Figure BDA0003453101860000062
将(30)代入到(29)中会得出:
Figure BDA0003453101860000063
其中
Figure BDA0003453101860000064
Figure BDA0003453101860000065
因此式(31)带入
Figure BDA00034531018600000615
然后(31)变成
Figure BDA0003453101860000066
设计F1(X1)为:
Figure BDA0003453101860000067
其中
Figure BDA0003453101860000068
F1(X1)是未知的,因此,利用RBFNN估算F1(X1),如下所示:
Figure BDA0003453101860000069
其中常数为δM>0;
因此,(32)改写为:
Figure BDA00034531018600000610
根据杨氏不等式,得到:
Figure BDA00034531018600000611
其中,设计参数为d1>0;
将(36)带入(35)中得到:
Figure BDA00034531018600000612
将虚拟控制律α2和自适应律
Figure BDA00034531018600000613
设计为:
Figure BDA00034531018600000614
其中k11,k1211和α12是正常数β=β12,β12是两个奇整数,满足0<β1<β2
将(38)代入(37)得到:
Figure BDA0003453101860000071
对于(20)-(22)、(24)和(38),取y2的时间导数得:
Figure BDA0003453101860000072
其中
Figure BDA0003453101860000073
表示连续函数;
由于
Figure BDA0003453101860000074
在紧集中遵循给定初始条件的最大值,因此,存在这样的函数
Figure BDA0003453101860000075
Figure BDA0003453101860000076
其中
Figure BDA0003453101860000077
利用杨氏不等式,有:
Figure BDA0003453101860000078
将(42)代到(41)中有:
Figure BDA0003453101860000079
步骤2:考虑李雅普诺夫函数V2
Figure BDA00034531018600000710
其中已知常数为
Figure BDA00034531018600000715
将V2的微分与(24)结合,得:
Figure BDA00034531018600000711
将(23)和(43)组合成(45)将导出:
Figure BDA00034531018600000712
与(30)相似,有:
Figure BDA00034531018600000713
将(47)代入(46)会得出:
Figure BDA00034531018600000714
设计F2(X2)作为:
Figure BDA0003453101860000081
其中X2=[x1,...,x4,yd2c]T
将(49)代入(48)得到:
Figure BDA0003453101860000082
从(49)中知道F2(X2)也是未知的,因此,F2(X2)由以下RBFNN近似:
Figure BDA0003453101860000083
然后,(51)进一步重新表述为:
Figure BDA0003453101860000084
与(36)类似,以下不等式成立:
Figure BDA0003453101860000085
其中,设计参数为d2>0;
将(53)代入到(52)中得到:
Figure BDA0003453101860000086
与(38)类似,将虚拟控制律α3和自适应律
Figure BDA0003453101860000087
设计为:
Figure BDA0003453101860000088
其中k21,k2221和α22为正常数;
将(55)代入(54)得:
Figure BDA0003453101860000089
与(42)相似,得出:
Figure BDA00034531018600000810
其中函数
Figure BDA00034531018600000811
将(57)代入(56)得到:
Figure BDA0003453101860000091
步骤3:构造Lyapunov函数V3,如下式所示
Figure BDA0003453101860000092
其中已知常数为
Figure BDA00034531018600000912
V3的时间导数与(24)结合得到:
Figure BDA0003453101860000093
将(23)和(58)代入(57)得到:
Figure BDA0003453101860000094
与(30)类似,得到以下不等式:
Figure BDA0003453101860000095
将(62)与(61)结合,得:
Figure BDA0003453101860000096
将F3(X3)定义为:
F3(X3)=b1x3+b2x2x4+b3x2+3v3+v2 (64)
其中X3=[x2,x3,x42c3c]T
然后,(63)如下式所示:
Figure BDA0003453101860000097
F3(X3)是未知的;因此,存在这样一个
Figure BDA0003453101860000098
Figure BDA0003453101860000099
与(36)类似,它得到:
Figure BDA00034531018600000910
其中,设计参数为d3>0
然后,(65)表述为:
Figure BDA00034531018600000911
将实际控制器uq和自适应律
Figure BDA0003453101860000101
设计为:
Figure BDA0003453101860000102
其中k31,k3232和α32为正常数;
将(69)代入(68)得出:
Figure BDA0003453101860000103
步骤4:选择Lyapunov函数V4为:
Figure BDA0003453101860000104
其中已知常数为
Figure BDA00034531018600001013
V4与(24)的时间导数得:
Figure BDA0003453101860000105
将(23)和(70)代入(72)中得到:
Figure BDA0003453101860000106
与(30)类似,以下关系成立:
Figure BDA0003453101860000107
那么,(73)简化为:
Figure BDA0003453101860000108
设函数F4(X4)为:
F4(X4)=c1x4+c2x2x3+3v4 (76)
其中X2=[x2,x3,x4]T
那么,(75)可以构建为:
Figure BDA0003453101860000109
易知函数F4(X4)不确定,因此,存在
Figure BDA00034531018600001010
使得:
Figure BDA00034531018600001011
与(36)类似,得到:
Figure BDA00034531018600001012
其中,设计参数为d4>0;
将(79)代入(77)中得到:
Figure BDA0003453101860000111
将实际控制器ud和自适应律
Figure BDA0003453101860000112
设计为:
Figure BDA0003453101860000113
其中k41,k4241和α42为正常数;
将(81)代入(80)得到:
Figure BDA0003453101860000114
本发明的有益效果:与现有技术相比,本发明的效果如下:
1)通过引入跟踪误差的函数变换,将输出受限的永磁同步电机系统转化为一种新型的无约束永磁同步电机系统。与现有技术中基于分段非对称BLF的约束反演控制器相比,基于转换的方案对于构造具有非对称约束的非线性系统的反演控制器是方便而直接的;
2)与现有技术中的渐近控制结果不同,本发明通过引入一阶滤波器来规避“复杂性爆炸”,为无约束永磁同步电机系统设计了一种有限时间动态面控制方案与反演方法相结合,并在设计的控制器和自适应律中包含分数次幂项的方案,这种设计不仅保证了较高的跟踪精度和收敛速度,而且具有良好的抗干扰能力;
3)与现有技术中的有限时间控制结果不同,本发明通过引入精细的Lyapunov-Krasovskii泛函,研究了有限时间稳定控制的时滞效应。因此,所设计的有限时间控制方案适用于时滞同时存在的实际应用。
附图说明
图1为永磁同步电机系统控制原理结构示意图;
图2为输出信号x1和参考轨迹yd曲线;
图3为跟踪误差s1结果;
图4为x2的结果说明图;
图5为iq和id的结果说明;
图6为控制器uq和ud的轨迹。
具体实施方式
下面结合具体的实施例对本发明进行进一步介绍。
实施例1:如图1-6所示,一种具有输出约束的PMSM系统有限时间动态面控制方法,包括以下步骤:
A系统说明(d-q)坐标系下永磁同步电机的动力学模型可以表述为:
Figure BDA0003453101860000121
其中,ω为转子角速度(rad/s),θ为转子角度(°),iq为q-轴电流(A),idd-轴电流(A),uq为q-轴电压(V),ud为d-轴电压(V),J为转动惯量(kg·m2),B为摩擦系数(N/(rad/s)),
Figure BDA00034531018600001212
为永磁通量(Wb),Rs为定子线圈电阻(Ω),np为极对数,Lq为q-轴线圈电感(H),Ld为d-轴线圈电感(H),TL为负载力矩(N·m);
定义变量x1=θ,x2=ω,x3=iq,x4=id并考虑到时间延迟和非对称时变输出约束,则(1)可简化为:
Figure BDA0003453101860000122
受以下输出约束:
Figure BDA0003453101860000123
其中,时变函数F11(t)>0和F12(t)>0表示已知约束边界,x1(t)表示输出变量,Δfi(x(t-τi)),i=1,...,4是时滞项,x=(x1,x2,x3,x4)T∈R4是(2)的整个状态,τi,i=1,...,4是时间常数,及
Figure BDA00034531018600001213
a2=3np(Ld-Lq)/2,b1=-Rs/Lq,b2=-npLd/Lq
Figure BDA00034531018600001214
b4=1/Lq,c1=-Rs/Ld,c2=npLq/Ld,c3=1/Ld
对于永磁同步电机系统,首先涉及时间延迟Δfi(x(t-τi)),i=1,...,4,由(2)和(3)表示的不对称输出约束。与现有技术不同,包含的时间延迟项由整体状态构成,更符合实际情况。与常数约束和对称约束相比,(3)中关于输出变量的上下边界是时变的和不对称的,更具一般性。
本发明旨在设计一种具有有限时间特性的神经自适应动态表面控制器,以确保:
(a)跟踪误差s1=x1-yd在有限时间内缩小到原点的一个小邻域,并且闭环系统的整个信号是有界的;(b)输出信号x1(t)需要满足(3)所示的关系。为了实现这些目标,假设和引理如下:
假设1:参考信号yd(t)及其n阶导数
Figure BDA0003453101860000124
(n=0,...,4)有界且连续;约束函数F11(t),F12(t)及其k次导数
Figure BDA0003453101860000125
(k=0,...,4)是有界的和连续的。
引理1:一个连续函数
Figure BDA00034531018600001217
由f(0,...,0)=0给出,其中
Figure BDA00034531018600001216
(i=1,2,...,n,mi>0),有光滑正函数
Figure BDA00034531018600001215
满足ωi(0)=0,这样
Figure BDA0003453101860000126
由引理1可知,系统(2)的时滞项Δfi(x(t-τi)),i=1,...,4可用
Figure BDA0003453101860000127
表示。然后,根据杨氏不等式,有:
Figure BDA0003453101860000128
引理2:考虑
Figure BDA0003453101860000129
Figure BDA00034531018600001210
是常数。然后,对于实变量y和z,以下不等式成立:
Figure BDA00034531018600001211
引理3:考虑实数r=1,....,m,和xr∈R,r=1,....,m,,得到:
Figure BDA0003453101860000131
定义1:给出非线性系统
Figure BDA0003453101860000132
其中f(λ)表示系统状态为λ∈Rn.的光滑函数。如果对于所有初始条件λ(t0)=λ0,有ρ>0和稳定时间T(ρ,λ0)<∞,使所有t≥t0+T,稳定时间为||λ(t)||<ρ,,则平衡点为λ=0的非线性系统
Figure BDA0003453101860000133
称为有限时间半全局稳定。
引理4:对于非线性系统
Figure BDA0003453101860000134
如果存在光滑正定函数V(λ)和标量a>0,
Figure BDA0003453101860000135
b>0,σ>0和
Figure BDA00034531018600001327
则:
Figure BDA0003453101860000136
那么非线性系统
Figure BDA0003453101860000137
在有限时间内是稳定的,其中稳定时间可近似为:
Figure BDA0003453101860000138
其中t≥T,存在:
Figure BDA0003453101860000139
引理5:对于每个变量,
Figure BDA00034531018600001310
存在φ1<φ2,φ12是奇数整数。一个不等式适用于以下情况:
Figure BDA00034531018600001311
其中
Figure BDA00034531018600001312
和γ2=(2φ-1-2(1+φ)(φ-1))/(1+φ)>0;
引理6:对于
Figure BDA00034531018600001322
有一个集合Λ,它由
Figure BDA00034531018600001324
给出。然后,对于
Figure BDA00034531018600001329
满足不等式
Figure BDA00034531018600001323
此后,当上下文中没有混淆时,有时会省略函数变量。
神经网络系统和函数逼近:使用RBFNNs技术,能够将未知非线性函数在闭合集中以任意精度估计。因此,有:
Figure BDA00034531018600001313
其中Z=[z1,z2,…,zn]T是输入向量,
Figure BDA00034531018600001325
是期望值。RBFNN的权重向量,l>1是节点数,δ(Z)是满足|δ(Z)|<δM,的估计误差。δM表示未知的有界参数。
Figure BDA00034531018600001314
是基函数向量,其中选择
Figure BDA00034531018600001315
作为应用的高斯函数:
Figure BDA00034531018600001316
其中μi=[μi1,...,μim]和
Figure BDA00034531018600001328
分别为域中心和高斯函数的宽度。
让理想权重向量
Figure BDA00034531018600001326
为:
Figure BDA00034531018600001317
其中
Figure BDA00034531018600001318
表示更新的权重向量;因此:
Figure BDA00034531018600001319
其中
Figure BDA00034531018600001320
θi和||·||分别是未知变量和·的2-范式;
(3)设计有限时间动态表面控制
A、非线性误差相关转换函数:引入一个非线性转换函数,将输出约束系统式(3)转换为非约束系统;
定义2:非线性转换函数可构造为:
Figure BDA00034531018600001321
其中,跟踪误差s1=x1-yd,S1表示转换误差,F11(t)>0和F12(t)>0表示平滑时变函数,假设F11(t),F12(t)满足以下关系:|F11(t)|<C0,|F12(t)|<C1,其中常数C0>0和C1>0。
从式(15)中看出,对于满足F11(0)<s1(0)<F12(0)的每个初始值,当S1有界为t∈[0,+∞)时,确保了s1的有界性和约束性,F11(t)和F12(t)将缩写为F11和F12
区分S1给出:
Figure BDA0003453101860000141
其中
Figure BDA0003453101860000142
Figure BDA0003453101860000143
利用(3)和(16),得到输出无约束子系统,如下式所示:
Figure BDA0003453101860000144
B、设计神经自适应有限时间控制器:
设误差坐标变换,如下所示:
v1=S1,v2=x22c (20)
v3=x33c,v4=x4,
其中αic,i=2,3表示以下一阶滤波器的输出:
Figure BDA0003453101860000145
其中εi表示时间常数,虚拟信号αi在后面设计;
类似地,过滤器误差yi定义为:
yi=αici,i=2,3 (22)
集成(2)和(19),在(20)中vi,i=1,...,4求导得:
Figure BDA0003453101860000146
引入估算误差
Figure BDA0003453101860000147
如下所示:
Figure BDA0003453101860000148
其中变量
Figure BDA0003453101860000149
代表θi的估计值;
神经自适应有限时间控制器的设计步骤如下:
步骤1:选择Lyapunov函数V1为:
Figure BDA00034531018600001410
Lyapunov-Krasovskii函数VM为:
Figure BDA00034531018600001411
其中已知
Figure BDA00034531018600001413
为常数;
取VM在(26)的时间导数,可以得:
Figure BDA00034531018600001412
其中,参数Γ>0,和正函数ωik以消除时间延迟;
然后,将(26)中的V1导数与(24)结合,得到:
Figure BDA0003453101860000151
将(23)合并到(28)会得出:
Figure BDA0003453101860000152
根据(4),得出:
Figure BDA0003453101860000153
将(30)代入到(29)中会得出:
Figure BDA0003453101860000154
其中
Figure BDA0003453101860000155
Figure BDA0003453101860000156
因此式(31)带入
Figure BDA00034531018600001516
然后(31)变成
Figure BDA0003453101860000157
设计F1(X1)为:
Figure BDA0003453101860000158
其中
Figure BDA0003453101860000159
F1(X1)是未知的,因此,利用RBFNN估算F1(X1),如下所示:
Figure BDA00034531018600001510
其中常数为δM>0;
因此,(32)改写为:
Figure BDA00034531018600001511
根据杨氏不等式,得到:
Figure BDA00034531018600001512
其中,设计参数为d1>0;
将(36)带入(35)中得到:
Figure BDA00034531018600001513
将虚拟控制律α2和自适应律
Figure BDA00034531018600001514
设计为:
Figure BDA00034531018600001515
其中k11,k1211和α12是正常数β=β12,β12是两个奇整数,满足0<β1<β2
将(38)代入(37)得到:
Figure BDA0003453101860000161
对于(20)-(22)、(24)和(38),取y2的时间导数得:
Figure BDA0003453101860000162
其中
Figure BDA0003453101860000163
表示连续函数;
由于
Figure BDA0003453101860000164
在紧集中遵循给定初始条件的最大值,因此,存在这样的函数
Figure BDA0003453101860000165
Figure BDA0003453101860000166
其中
Figure BDA0003453101860000167
利用杨氏不等式,有:
Figure BDA0003453101860000168
将(42)代到(41)中有:
Figure BDA0003453101860000169
步骤2:考虑李雅普诺夫函数V2
Figure BDA00034531018600001610
其中已知常数为
Figure BDA00034531018600001614
将V2的微分与(24)结合,得:
Figure BDA00034531018600001611
将(23)和(43)组合成(45)将导出:
Figure BDA00034531018600001612
与(30)相似,有:
Figure BDA00034531018600001613
将(47)代入(46)会得出:
Figure BDA0003453101860000171
设计F2(X2)作为:
Figure BDA0003453101860000172
其中X2=[x1,...,x4,yd2c]T
将(49)代入(48)得到:
Figure BDA0003453101860000173
从(49)中知道F2(X2)也是未知的,因此,F2(X2)由以下RBFNN近似:
Figure BDA0003453101860000174
然后,(51)进一步重新表述为:
Figure BDA0003453101860000175
与(36)类似,以下不等式成立:
Figure BDA0003453101860000176
其中,设计参数为d2>0;
将(53)代入到(52)中得到:
Figure BDA0003453101860000177
与(38)类似,将虚拟控制律α3和自适应律
Figure BDA0003453101860000178
设计为:
Figure BDA0003453101860000179
其中k21,k2221和α22为正常数;
将(55)代入(54)得:
Figure BDA0003453101860000181
与(42)相似,得出:
Figure BDA0003453101860000182
其中函数
Figure BDA0003453101860000183
将(57)代入(56)得到:
Figure BDA0003453101860000184
步骤3:构造Lyapunov函数V3,如下式所示
Figure BDA0003453101860000185
其中已知常数为
Figure BDA00034531018600001813
V3的时间导数与(24)结合得到:
Figure BDA0003453101860000186
将(23)和(58)代入(57)得到:
Figure BDA0003453101860000187
与(30)类似,得到以下不等式:
Figure BDA0003453101860000188
将(62)与(61)结合,得:
Figure BDA0003453101860000189
将F3(X3)定义为:
F3(X3)=b1x3+b2x2x4+b3x2+3v3+v2 (64)
其中X3=[x2,x3,x42c3c]T
然后,(63)如下式所示:
Figure BDA00034531018600001810
F3(X3)是未知的;因此,存在这样一个
Figure BDA00034531018600001811
Figure BDA00034531018600001812
与(36)类似,它得到:
Figure BDA0003453101860000191
其中,设计参数为d3>0
然后,(65)表述为:
Figure BDA0003453101860000192
将实际控制器uq和自适应律
Figure BDA0003453101860000193
设计为:
Figure BDA0003453101860000194
其中k31,k3232和α32为正常数;
将(69)代入(68)得出:
Figure BDA0003453101860000195
步骤4:选择Lyapunov函数V4为:
Figure BDA0003453101860000196
其中已知常数为
Figure BDA00034531018600001911
V4与(24)的时间导数得:
Figure BDA0003453101860000197
将(23)和(70)代入(72)中得到:
Figure BDA0003453101860000198
与(30)类似,以下关系成立:
Figure BDA0003453101860000199
那么,(73)简化为:
Figure BDA00034531018600001910
设函数F4(X4)为:
F4(X4)=c1x4+c2x2x3+3v4 (76)
其中X2=[x2,x3,x4]T
那么,(75)可以构建为:
Figure BDA0003453101860000201
易知函数F4(X4)不确定,因此,存在
Figure BDA0003453101860000202
使得:
Figure BDA0003453101860000203
与(36)类似,得到:
Figure BDA0003453101860000204
其中,设计参数为d4>0
将(79)代入(77)中得到:
Figure BDA0003453101860000205
将实际控制器ud和自适应律
Figure BDA0003453101860000206
设计为:
Figure BDA0003453101860000207
其中k41,k4241和α42为正常数;
将(81)代入(80)得到:
Figure BDA0003453101860000208
至此,完成了永磁同步电机系统控制器的设计过程。为了更清楚地说明控制方案,方框图如图1所示。
稳定性分析:对于任何给定p>0,将紧集定义为:
Figure BDA0003453101860000209
定理1:在假设1下,针对永磁同步电机系统(2),设计了由控制律α23,uq,ud和自适应律
Figure BDA0003453101860000211
i=1,...,4组成的神经自适应有限时间动态表面控制方法。如果初始条件满足Ωi,i=1,...,4,F11(0)<s1(0)<F12(0)和yd∈(-d,d),则可确保所有控制目标。
证明:整体李雅普诺夫函数V被选择为:
Figure BDA0003453101860000212
结合(82),导出V(t)的导数:
Figure BDA0003453101860000213
利用杨氏不等式(10)和(24),得到:
Figure BDA0003453101860000214
其中
Figure BDA0003453101860000215
Figure BDA00034531018600002118
然后,(86)可以被重构为:
Figure BDA0003453101860000216
此外,在引理2中取
Figure BDA0003453101860000217
Figure BDA0003453101860000218
得出|z|β≤(1-β)g+|z|.,让z分别为ΓVM/2和
Figure BDA0003453101860000219
加εi≤2,i=2,3,有:
Figure BDA00034531018600002110
将(88)与(87)结合,得出:
Figure BDA00034531018600002111
其中b0=min{2ki1,(1/ε2-1/2),(1/ε3-1/2),αi1,Γ/21≤i≤4},
Figure BDA00034531018600002112
Figure BDA00034531018600002119
使a0=min{2βki2,(1/ε2-1/2),(1/ε3-1/2),
Figure BDA00034531018600002113
|1≤β。≤4},然后,通过使用引理3,(89)可以表示为:
Figure BDA00034531018600002114
需要选择合适的参数来保证a0>0和b0>0。对于(90)中的最后一个变量q0,需要通过引理2.7进行进一步讨论,如下所示:①对于
Figure BDA00034531018600002122
结合M(x)≥0,得出q0≤0;②对于v1∈Λ,我们可以得到
Figure BDA00034531018600002120
Figure BDA00034531018600002121
因此,v1q0都是有界的。此外,得到一个正常数
Figure BDA00034531018600002115
满足
Figure BDA00034531018600002116
那么(90)可以重写为:
Figure BDA00034531018600002117
通过使用引理4到(91),可以得到:
(I)非线性系统(2)在有限时间内是稳定的;
(II)存在一个常数T(称为设定时间),因此,对于任何
Figure BDA00034531018600002123
和所有t≥T,以下不等式成立:
Figure BDA0003453101860000221
其中,稳定时间T可近似为:
Figure BDA0003453101860000222
根据(84),得到:
Figure BDA0003453101860000223
这个不等式表示vi,i=1,...,4是有界的。同样,通过组合(15)可以得到y2,y3
Figure BDA0003453101860000224
i=1,...,4是有界的。利用(24),可以进一步得到
Figure BDA0003453101860000225
i=1,...,4是有界的。此外,从(15)中可以明显看出(F11(t)+s1(t))(F12(t)-s1(t))<(F11(t)+F12(t))2≤(C0+C1)2=Q。因此,结合(15)、(20)和(94),可以得出:
Figure BDA0003453101860000226
(95)表示s1有界。因为s1=x1-yd和参考信号yd是严格有界的,可以得出x1是有界的。因此,ξ1是有界的。然后,可以得出(38)中构造的α2及其导数
Figure BDA0003453101860000227
是有界的。结合y2=α2c2和(41),我们可以得出α2c
Figure BDA0003453101860000228
是有界的。考虑到v2=x22c,可以得出x2是有界的。同样,可以推断x3,x43,
Figure BDA0003453101860000229
α3c,、
Figure BDA00034531018600002210
uduq是有界的。总之,可以得出结论,由此产生的系统的所有信号都是有界的。
特别地,(95)揭示了通过充分调节设计参数,输出变量x1可以在有限时间后紧跟参考信号yd
此外,结合S1=s1/[(F11(t)+s1)(F12(t)-s1)],可以得出S1→±∞当且仅当s1→-F11(t)或s1→F12(t).。此外,基于v1=S1有界的事实,对于满足F11(0)<s1(0)<F12(0).的每个初始条件,它有-F11(t)<s1(t)<F12(t),
Figure BDA00034531018600002211
然后,考虑s1=x1-yd,它给出
Figure BDA00034531018600002212
对于满足yd(0)-F11(0)<x1(0)<yd(0)+F12(0).的任何初始条件。到目前为止,稳定性分析证明已经完成。
跟踪误差s1是控制能力的直观度量。根据公式(95),可以观察到跟踪误差s1的大小取决于变量Q,
Figure BDA00034531018600002213
b0。具体地说,从Q可以得出,跟踪误差的性能可以通过边界F11(t)F12(t)来实现。随后,使用
Figure BDA00034531018600002214
b0的定义,通过减少参数dii1i2,i=1,...,4和增加参数ki,
Figure BDA00034531018600002221
可以确保
Figure BDA00034531018600002215
的大小足够小,b0的值足够大。结果表明,跟踪误差s1可以调节得很小。然而,从(38)、(55)、(69)、(81)可以明显看出,增加ki,
Figure BDA00034531018600002222
和减少di可能导致控制信号的幅度大。因此,应考虑与控制努力相关的系统性能的权衡。
神经自适应有限时间动态表面控制,以解决时滞和非对称时变输出约束问题,而现有技术中并未涉及。此外,所设计的控制律和自适应律(包含分数次幂项αi2θ2β-1,i=1...,4)实现了永磁同步电机系统在有限时间内的快速稳定控制。因此,本发明设计的控制方案更具实用性。
仿真验证:为系统(2)提供了两个仿真案例,以验证控制方案的有效性。两个案例是:案例1:时滞项为Δfi=0,i=1,2,3,4,这意味着时滞对控制性能没有影响;
案例2:选择延时项如下:
Figure BDA00034531018600002216
Figure BDA00034531018600002217
Figure BDA00034531018600002218
Figure BDA00034531018600002219
选择永磁同步电机参数为J=0.003798Kg·m2,B=0.001158N·m/(rad/s),TL=1.5,
Figure BDA00034531018600002220
Ld=0.00285H,np=3,Lq=0.00315H,Rs=0.68Ω,时变函数为F11=0.8+2-0.3t,F12=0.7+2-0.2t。参考信号选择为yd=0.5(sin(t)+sin(2t));在仿真中,初始状态选择为xi(0)=0,i=1,...,4.。每个RBFNN包含11个节点,中心间隔为[-1111],宽度为10。设计控制参数选择为:α2c(0)=0,α3c(0)=0.5,
Figure BDA0003453101860000231
ε2=ε3=0.01,
Figure BDA0003453101860000232
ki1=10,ki2=30,αi1=38,di=0.1,αi2=0.05,1≤i≤4,β=97/101。
仿真结果如图2-6所示。图2表明,在情况1-2下,输出信号可以跟踪理想的信号曲线。图3给出了跟踪误差信号曲线。从图2-3可以看出,在实际操作中,系统输出和跟踪误差是有界的,并且不违反它们的约束。图4给出了实际控制信号和状态的响应。图5-6显示了实际控制信号和状态的轨迹。从以上结果可以看出,所设计的解决方案可以完美地执行,其性能趋于令人满意。
本发明基于神经自适应动态表面控制方法,研究了具有时滞和非对称时变输出约束的永磁同步电机系统的有限时间跟踪控制问题。首先,采用基于跟踪误差约束的非线性转换函数。这种策略不仅可以克服分段非对称BLF,而且输出约束系统的设计也可以简化为一般的无约束系统。
本发明将合适的Lyapunov-Krasovskii泛函和有限时间控制方法融合到通用的动态表面控制框架中,提出了一种新的具有分数次幂项自适应律的快速稳定控制方案,解决了有限时间内的“复杂性爆炸”和时滞问题。然后,给出了临界稳定性分析和仿真结果,验证了系统的所有信号都是有界的,以及所设计控制方案的有效性。我们未来的工作将这种设计扩展到通常会遇到物理限制和时间延迟的场景中,例如车辆、电动电梯、机器人和机床。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内,因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (1)

1.一种具有输出约束的PMSM系统有限时间动态面控制方法,其特征在于:该方法包括以下步骤:
(1)定义变量x1=θ,x2=ω,x3=iq,x4=id,对(d-q)坐标系下永磁同步电机的动力学模型进行简化,得到如下式:
Figure FDA0003453101850000011
式(2)受以下输出约束:
Figure FDA0003453101850000012
其中,时变函数F11(t)>0和F12(t)>0表示已知约束边界,x1(t)表示输出变量,Δfi(x(t-τi)),i=1,...,4是时滞项,x=(x1,x2,x3,x4)T∈R4是式(2)的整体状态,τi,i=1,...,4是时间常数,
Figure FDA0003453101850000013
a2=3np(Ld-Lq)/2,b1=-Rs/Lq,b2=-npLd/Lq
Figure FDA0003453101850000014
b4=1/Lq,c1=-Rs/Ld,c2=npLq/Ld,c3=1/Ld;ω为转子角速度,θ为转子角度,iq为q-轴电流,id为d-轴电流,uq为q-轴电压,ud为d-轴电压,J为转动惯量,B为摩擦系数,
Figure FDA0003453101850000015
为永磁通量,Rs为定子线圈电阻,np为极对数,Lq为q-轴线圈电感,Ld为d-轴线圈电感,TL为负载力矩;
设1:参考信号yd(t)及其n阶导数
Figure FDA0003453101850000016
有界且连续;约束函数F11(t),F12(t)及其k阶导数
Figure FDA0003453101850000017
是有界的和连续的;
引理1:一个连续函数
Figure FDA0003453101850000018
由f(0,...,0)=0给出,其中
Figure FDA0003453101850000019
有光滑正函数
Figure FDA00034531018500000110
满足ωi(0)=0,这样
Figure FDA00034531018500000111
由引理1可知,式(2)的时滞项Δfi(x(t-τi)),i=1,...,4用
Figure FDA00034531018500000112
表示,然后,根据杨氏不等式,有:
Figure FDA00034531018500000113
引理2:考虑
Figure FDA00034531018500000114
Figure FDA00034531018500000115
是常数,然后,对于实变量y和z,以下不等式成立:
Figure FDA00034531018500000116
引理3:考虑实数r=1,....,m,和xr∈R,r=1,....,m,得到:
Figure FDA00034531018500000117
定义1:给出非线性系统
Figure FDA00034531018500000118
其中f(λ)表示系统状态为λ∈Rn.的光滑函数,如果对于所有初始条件λ(t0)=λ0,有ρ>0和稳定时间T(ρ,λ0)<∞,使所有t≥t0+T,稳定时间为||λ(t)||<ρ,,则平衡点为λ=0的非线性系统
Figure FDA00034531018500000124
称为有限时间半全局实际稳定;
引理4:对于非线性系统
Figure FDA00034531018500000119
如果存在光滑正定函数V(λ)和标量a>0,
Figure FDA00034531018500000120
b>0,σ>0和
Figure FDA00034531018500000121
则:
Figure FDA00034531018500000122
那么非线性系统
Figure FDA00034531018500000123
在有限时间内是稳定的,其中稳定时间可近似为:
Figure FDA0003453101850000021
其中t≥T,存在:
Figure FDA0003453101850000022
引理5:对于每个变量,
Figure FDA0003453101850000023
φ=φ12,存在φ1<φ2,φ12是奇数整数。一个不等式适用于以下情况:
Figure FDA0003453101850000024
其中
Figure FDA0003453101850000025
Figure FDA0003453101850000026
引理6:对于
Figure FDA0003453101850000027
有一个集合Λ,它由
Figure FDA0003453101850000028
给出,然后,对于
Figure FDA0003453101850000029
满足不等式
Figure FDA00034531018500000210
(2)使用RBFNNs技术,对未知非线性函数在闭合集中以任意精度估计,因此,有:
Figure FDA00034531018500000211
其中Z=[z1,z2,…,zn]T是输入向量,
Figure FDA00034531018500000212
是期望值,RBFNN的权重向量,l>1是节点数,δ(Z)是满足|δ(Z)|<δM,的估计误差,δM表示未知的有界参数,
Figure FDA00034531018500000213
是基函数向量,其中选择
Figure FDA00034531018500000214
作为应用的高斯函数:
Figure FDA00034531018500000215
其中μi=[μi1,...,μim]和
Figure FDA00034531018500000216
分别为域中心和高斯函数的宽度;
让理想权重向量
Figure FDA00034531018500000224
为:
Figure FDA00034531018500000217
其中
Figure FDA00034531018500000218
表示更新的权重向量;
因此:
Figure FDA00034531018500000219
其中
Figure FDA00034531018500000220
θi和||·||分别是未知变量和·的2-范式;
(3)设计有限时间动态表面控制
A、非线性误差相关转换函数:引入一个非线性转换函数,将输出约束系统式(3)转换为非约束系统;
定义2:非线性转换函数可构造为:
Figure FDA00034531018500000221
其中,跟踪误差s1=x1-yd,S1表示转换误差,F11(t)>0和F12(t)>0表示平滑时变函数,假设F11(t),F12(t)满足以下关系:|F11(t)|<C0,|F12(t)|<C1,其中常数C0>0和C1>0。
从式(15)中看出,对于满足F11(0)<s1(0)<F12(0)的每个初始值,当S1有界为t∈[0,+∞)时,确保了s1的有界性和约束性,F11(t)和F12(t)将缩写为F11和F12
区分S1给出:
Figure FDA00034531018500000222
其中
Figure FDA00034531018500000223
Figure FDA0003453101850000031
利用(3)和(16),得到输出无约束子系统,如下式所示:
Figure FDA0003453101850000032
B、设计神经自适应有限时间控制器:
设误差坐标变换,如下所示:
Figure FDA0003453101850000033
其中αic,i=2,3表示以下一阶滤波器的输出:
Figure FDA0003453101850000034
其中εi表示时间常数,虚拟信号αi在后面设计;
类似地,过滤器误差yi定义为:
yi=αici,i=2,3 (22)
集成(2)和(19),在(20)中vi,i=1,...,4求导得:
Figure FDA0003453101850000035
引入估算误差
Figure FDA0003453101850000036
如下所示:
Figure FDA0003453101850000037
其中变量
Figure FDA0003453101850000038
代表θi的估计值;
神经自适应有限时间控制器的设计步骤如下:
步骤1:选择Lyapunov函数V1为:
Figure FDA0003453101850000039
Lyapunov-Krasovskii函数VM为:
Figure FDA00034531018500000310
其中已知l1>0为常数;
取VM在(26)的时间导数,得:
Figure FDA00034531018500000313
其中,参数Γ>0,和正函数ωik以消除时间延迟;
然后,将(26)中的V1导数与(24)结合,得到:
Figure FDA00034531018500000311
将(23)合并到(28)会得出:
Figure FDA00034531018500000312
根据(4),得出:
Figure FDA0003453101850000041
将(30)代入到(29)中会得出:
Figure FDA0003453101850000042
其中
Figure FDA0003453101850000043
Figure FDA0003453101850000044
因此式(31)带入
Figure FDA0003453101850000045
然后(31)变成
Figure FDA0003453101850000046
设计F1(X1)为:
Figure FDA0003453101850000047
其中
Figure FDA0003453101850000048
F1(X1)是未知的,因此,利用RBFNN估算F1(X1),如下所示:
Figure FDA0003453101850000049
其中常数为δM>0;
因此,(32)改写为:
Figure FDA00034531018500000410
根据杨氏不等式,得到:
Figure FDA00034531018500000411
其中,设计参数为d1>0;
将(36)带入(35)中得到:
Figure FDA00034531018500000412
将虚拟控制律α2和自适应律
Figure FDA00034531018500000413
设计为:
Figure FDA00034531018500000414
其中k11,k1211和α12是正常数β=β12,β12是两个奇整数,满足0<β1<β2
将(38)代入(37)得到:
Figure FDA00034531018500000415
对于(20)-(22)、(24)和(38),取y2的时间导数得:
Figure FDA0003453101850000051
其中
Figure FDA0003453101850000052
表示连续函数;
由于
Figure FDA0003453101850000053
在紧集中遵循给定初始条件的最大值,因此,存在这样的函数
Figure FDA0003453101850000054
Figure FDA0003453101850000055
其中
Figure FDA0003453101850000056
利用杨氏不等式,有:
Figure FDA0003453101850000057
将(42)代到(41)中有:
Figure FDA0003453101850000058
步骤2:考虑李雅普诺夫函数V2
Figure FDA0003453101850000059
其中已知常数为l2>0;
将V2的微分与(24)结合,得:
Figure FDA00034531018500000510
将(23)和(43)组合成(45)将导出:
Figure FDA00034531018500000511
与(30)相似,有:
Figure FDA00034531018500000512
将(47)代入(46)会得出:
Figure FDA00034531018500000513
设计F2(X2)作为:
Figure FDA00034531018500000514
其中X2=[x1,...,x4,yd2c]T
将(49)代入(48)得到:
Figure FDA0003453101850000061
从(49)中知道F2(X2)也是未知的,因此,F2(X2)由以下RBFNN近似:
Figure FDA0003453101850000062
然后,(51)进一步重新表述为:
Figure FDA0003453101850000063
与(36)类似,以下不等式成立:
Figure FDA0003453101850000064
其中,设计参数为d2>0;
将(53)代入到(52)中得到:
Figure FDA0003453101850000065
与(38)类似,将虚拟控制律α3和自适应律
Figure FDA0003453101850000066
设计为:
Figure FDA0003453101850000067
其中k21,k2221和α22为正常数;
将(55)代入(54)得:
Figure FDA0003453101850000068
与(42)相似,得出:
Figure FDA0003453101850000069
其中函数
Figure FDA00034531018500000610
将(57)代入(56)得到:
Figure FDA00034531018500000611
步骤3:构造Lyapunov函数V3,如下式所示
Figure FDA0003453101850000071
其中已知常数为l3>0;
V3的时间导数与(24)结合得到:
Figure FDA0003453101850000072
将(23)和(58)代入(57)得到:
Figure FDA0003453101850000073
与(30)类似,得到以下不等式:
Figure FDA0003453101850000074
将(62)与(61)结合,得:
Figure FDA0003453101850000075
将F3(X3)定义为:
F3(X3)=b1x3+b2x2x4+b3x2+3v3+v2 (64)
其中X3=[x2,x3,x42c3c]T
然后,(63)如下式所示:
Figure FDA0003453101850000076
F3(X3)是未知的;因此,存在这样一个
Figure FDA0003453101850000077
Figure FDA0003453101850000078
与(36)类似,它得到:
Figure FDA0003453101850000079
其中,设计参数为d3>0;
然后,(65)表述为:
Figure FDA00034531018500000710
将实际控制器uq和自适应律
Figure FDA00034531018500000711
设计为:
Figure FDA0003453101850000081
其中k31,k3232和α32为正常数;
将(69)代入(68)得出:
Figure FDA0003453101850000082
步骤4:选择Lyapunov函数V4为:
Figure FDA0003453101850000083
其中已知常数为l4>0;
取V4与(24)的时间导数得:
Figure FDA0003453101850000084
将(23)和(70)代入(72)中得到:
Figure FDA0003453101850000085
与(30)类似,以下关系成立:
Figure FDA0003453101850000086
那么,(73)简化为:
Figure FDA0003453101850000087
设函数F4(X4)为:
F4(X4)=c1x4+c2x2x3+3v4 (76)
其中X2=[x2,x3,x4]T
那么,(75)构建为:
Figure FDA0003453101850000088
易知函数F4(X4)不确定,因此,存在
Figure FDA0003453101850000089
使得:
Figure FDA00034531018500000810
与(36)类似,得到:
Figure FDA00034531018500000811
其中,设计参数为d4>0;
将(79)代入(77)中得到:
Figure FDA0003453101850000091
将实际控制器ud和自适应律
Figure FDA0003453101850000092
设计为:
Figure FDA0003453101850000093
其中k41,k4241和α42为正常数;
将(81)代入(80)得到:
Figure FDA0003453101850000094
CN202111678368.XA 2021-12-31 2021-12-31 一种具有输出约束的pmsm系统有限时间动态面控制方法 Active CN114280944B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111678368.XA CN114280944B (zh) 2021-12-31 2021-12-31 一种具有输出约束的pmsm系统有限时间动态面控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111678368.XA CN114280944B (zh) 2021-12-31 2021-12-31 一种具有输出约束的pmsm系统有限时间动态面控制方法

Publications (2)

Publication Number Publication Date
CN114280944A true CN114280944A (zh) 2022-04-05
CN114280944B CN114280944B (zh) 2024-02-13

Family

ID=80879850

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111678368.XA Active CN114280944B (zh) 2021-12-31 2021-12-31 一种具有输出约束的pmsm系统有限时间动态面控制方法

Country Status (1)

Country Link
CN (1) CN114280944B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114660946A (zh) * 2022-05-09 2022-06-24 电子科技大学 一种时滞成型过程系统的模糊自适应动态面控制方法
CN117215240A (zh) * 2023-09-28 2023-12-12 淮阴工学院 一种非线性机电伺服系统的预分配时间预定性能控制方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109873582A (zh) * 2019-01-29 2019-06-11 青岛大学 基于动态面的永磁同步电机有限时间位置跟踪控制方法
CN110336506A (zh) * 2019-08-20 2019-10-15 贵州大学 一种pmsm混沌系统神经网络反演控制方法
CN110347044A (zh) * 2019-07-15 2019-10-18 贵州大学 一种考虑输出约束的pmsm混沌系统神经网络动态面控制方法
US20200186069A1 (en) * 2018-12-07 2020-06-11 Trw Limited Method of characterising a synchronous motor
CN113031446A (zh) * 2021-03-15 2021-06-25 贵州大学 不确定时滞非线性系统非奇异性神经自适应跟踪控制方法
CN113064347A (zh) * 2021-03-15 2021-07-02 贵州大学 考虑非对称输入与输出约束的pmsm混沌系统自适应控制方法
GB202111727D0 (en) * 2021-05-13 2021-09-29 Univ Guizhou Adaptive event-triggered control method for permanent magnet synchronous generator coupled chaotic network system
CN113659895A (zh) * 2021-08-02 2021-11-16 青岛大学 基于指令滤波的永磁同步电动机全状态约束有限时间控制方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200186069A1 (en) * 2018-12-07 2020-06-11 Trw Limited Method of characterising a synchronous motor
CN109873582A (zh) * 2019-01-29 2019-06-11 青岛大学 基于动态面的永磁同步电机有限时间位置跟踪控制方法
CN110347044A (zh) * 2019-07-15 2019-10-18 贵州大学 一种考虑输出约束的pmsm混沌系统神经网络动态面控制方法
CN110336506A (zh) * 2019-08-20 2019-10-15 贵州大学 一种pmsm混沌系统神经网络反演控制方法
CN113031446A (zh) * 2021-03-15 2021-06-25 贵州大学 不确定时滞非线性系统非奇异性神经自适应跟踪控制方法
CN113064347A (zh) * 2021-03-15 2021-07-02 贵州大学 考虑非对称输入与输出约束的pmsm混沌系统自适应控制方法
GB202111727D0 (en) * 2021-05-13 2021-09-29 Univ Guizhou Adaptive event-triggered control method for permanent magnet synchronous generator coupled chaotic network system
CN113659895A (zh) * 2021-08-02 2021-11-16 青岛大学 基于指令滤波的永磁同步电动机全状态约束有限时间控制方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
孙建业;王志强;谷鑫;夏长亮;: "高速低载波比下永磁同步电机预测电流控制", 中国电机工程学报, no. 11 *
张钧星 等: "基于行为博弈进化算法的并联混合动力汽车控制策略参数优化", 《科学技术与工程》, vol. 17, no. 13 *
陈晓雷;林辉;吕帅帅;: "机载PMSM伺服作动系统约束反演控制", 电工技术学报, no. 21 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114660946A (zh) * 2022-05-09 2022-06-24 电子科技大学 一种时滞成型过程系统的模糊自适应动态面控制方法
CN117215240A (zh) * 2023-09-28 2023-12-12 淮阴工学院 一种非线性机电伺服系统的预分配时间预定性能控制方法

Also Published As

Publication number Publication date
CN114280944B (zh) 2024-02-13

Similar Documents

Publication Publication Date Title
Yu et al. Approximation-based discrete-time adaptive position tracking control for interior permanent magnet synchronous motors
CN111342720B (zh) 基于转矩观测的永磁同步电机自适应连续滑模控制方法
CN114280944A (zh) 一种具有输出约束的pmsm系统有限时间动态面控制方法
Luo et al. Chaos control of the permanent magnet synchronous motor with time-varying delay by using adaptive sliding mode control based on DSC
You et al. Adaptive neural network control using nonlinear information gain for permanent magnet synchronous motors
CN110112971B (zh) 一种基于有限时间动态面的异步电动机位置跟踪控制方法
CN106788046A (zh) 永磁同步电机命令滤波有限时间模糊控制方法
Liu et al. Second-order ESO-based current sensor fault-tolerant strategy for sensorless control of PMSM with B-phase current
Gabbi et al. Discrete-time sliding mode controller based on backstepping disturbance compensation for robust current control of PMSM drives
Yang et al. Speed sensorless control of a bearingless induction motor with combined neural network and fractional sliding mode
CN114706300B (zh) 具有扰动及输出约束的永磁同步电机系统有限时间控制法
Wang et al. An improved model-free active disturbance rejection deadbeat predictive current control method of pmsm based on data-driven
CN117691903A (zh) 基于rbf神经网络的永磁同步电机滑模控制方法
CN114567219A (zh) 一种永磁同步电机速度控制方法及装置
CN110378057B (zh) 一种内置式永磁同步电机抗干扰控制器及其设计方法
Gao et al. Neural network based dynamic surface integral nonsingular fast terminal sliding mode control for manipulators with disturbance rejection
CN111293941B (zh) 一种考虑铁损的永磁同步电动机有限时间动态面控制方法
Yingming et al. A learning observer-based control strategy for pmsm with position sensor fault in railway
CN112019111A (zh) 基于状态约束的永磁同步电动机随机系统模糊自适应控制方法
CN114519301B (zh) 具有时滞的非对称输出约束pmsm系统动态面跟踪控制方法
CN114172436B (zh) 基于观测器的异步电动机指令滤波离散控制方法
Lin et al. Improved Model-Free Sliding Mode Control of Linear Motor Based on Time-Varying Gain Model-Assisted Linear Extended State Observer
CN114244238B (zh) 双伺服电机系统转速同步容错控制方法
Liao et al. Low speed servo control for single Gimbal control moment gyros using back-stepping time-varying sliding model control
CN112039374B (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