CN114706300A - 具有扰动及输出约束的永磁同步电机系统有限时间控制法 - Google Patents

具有扰动及输出约束的永磁同步电机系统有限时间控制法 Download PDF

Info

Publication number
CN114706300A
CN114706300A CN202210096210.XA CN202210096210A CN114706300A CN 114706300 A CN114706300 A CN 114706300A CN 202210096210 A CN202210096210 A CN 202210096210A CN 114706300 A CN114706300 A CN 114706300A
Authority
CN
China
Prior art keywords
following
function
substituting
steps
time
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
CN202210096210.XA
Other languages
English (en)
Other versions
CN114706300B (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 CN202210096210.XA priority Critical patent/CN114706300B/zh
Publication of CN114706300A publication Critical patent/CN114706300A/zh
Application granted granted Critical
Publication of CN114706300B publication Critical patent/CN114706300B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了具有扰动及输出约束的永磁同步电机系统有限时间控制法,采用非线性变换函数将输出约束动态模型转换为无约束的动态模型;提出一种基于扰动观测器的有限时间控制策略估计不匹配的外界扰动,然后,通过结合神经网络来估计非线性函数和一阶滤波器来处理“复杂性爆炸”,构思出一种对永磁同步电机具有自适应学习能力的动态面控制方法。本发明有效解决有限时间反步框架下系统具有时变特性和不匹配的外部扰动的非对称输出约束的问题。

Description

具有扰动及输出约束的永磁同步电机系统有限时间控制法
技术领域
本发明属于具有时滞和非对称时变输出约束的PMSM系统控制技术领域,涉及具有扰动及输出约束的永磁同步电机系统有限时间控制法。
背景技术
与传统步进电机相比,由于永磁同步电机具有明显的优点,例如更高的运行效率,更好的功率密度,它们通常用于车辆,飞机,石油和化学工业。尽管如此,永磁同步电机是具有多变量耦合的非线性系统,并且容易受到某些不稳定性的影响,包括参数扰动和外在干扰。因此,设计一种有效的控制器来克服系统中存在的非线性,对于永磁同步电机的应用具有更深远的意义。在过去的几十年中,研究人员提出了许多关键的控制方法,例如滑模控制、自适应控制、反步控制等,以满足永磁同步电机控制系统的更高要求。
上述方法中,能够很好地将自适应神经网络和模糊逻辑系统控制结合的反步控制方法,对于具有任意不确定扰动的永磁同步电机构建高性能控制器十分有效。因此,在传统的反步方案中产生的虚拟控制律的多重微分导致“复杂性爆炸”,这增加了控制系统的计算负担。为了应对这一挑战,通过将卡尔曼滤波器、跟踪微分器、命令滤波器和一阶滤波器集成到经典的反步框架中,提出了许多关键方案。例如,提出扩展卡尔曼滤波器来评估周期性扰动及其导数。张钧星(张钧星等,具有外部扰动和约束输出的混沌永磁同步电机系统自适应神经动态表面控制[J].电气与电子工程的最新进展(以前称为电气与电子工程的最新专利),2020年)将从动态面控制方法衍生的命令扩展到永磁同步电机系统,以解决繁重的计算负担。尽管上文所提出的方法在一定程度上提高了系统的稳定性和鲁棒性,但它们仍然没有考虑到当前研究的前沿,如有限时间稳定性理论、扰动观测器和约束控制来加强系统性能。
由于有限时间控制方法能够从时间优化控制方案的角度提供出色的收敛性和优异的鲁棒性,因此设计了许多基于有限稳定性理论的有效策略来获得永磁同步电机的高收敛速度和响应。例如,有限时间控制方法被用来在有限时间内实现跟踪所需的电机位置。为了实现动态相应和抗干扰能力,一个有限时间控制器被设计以便系统输出在有限的时间内跟踪所需的速度信号。一种自适应有限时间控制方法被提出来保证复杂系统中的信号和所有状态误差都能有限时间收敛到有界区域。上述工作表明,通过有限时间方法可以实现更快速的有限时间跟踪误差收敛,这是在这项工作中构建针对永磁同步电机系统的有限时间控制方案的动机之一。
关于研究受约束的永磁同步电机控制系统,已经采用了许多有效的方法来处理它们,其中障碍Lyapunov函数通常被认为是将跟踪误差收敛到预设的起始紧集的有效的方法。利用两种类型的非对称和对称障碍Lyapunov函数来确保跟踪误差可以在有限的时间内被约束在起始点的小邻域内。一个积分型障碍Lyapunov函数被提出来处理基于输出约束的保守可行性条件。为了消除利用障碍Lyapunov函数引起的分段不对称性,提出一种非线性坐标转换函数,将具有状态约束的模型转换为无约束模型。利用非线性转换函数将全状态约束系统转变为无约束系统。针对永磁同步电机中的物理约束,提出一种有限时间H无穷大控制方法,以消除约束限制并得到出H无穷大性能。虽然前面的文献对约束系统做出了巨大贡献,但并没有考虑到损害受控系统稳定运行的外部干扰。
为了增强电机系统的抗干扰能力,许多研究人员已经考虑为增加扰动观测器。例如,集成了滑模扰动观测器来估计外部负载扭矩,以增强永磁同步电机的抗干扰和动态响应能力。一种基于跟踪微分器的扰动观测器被设计来增加开关增益,这可以抑制突然变化的干扰。为了增强对扰动的抵抗力,可以通过分数阶来估计不匹配的扰动。为了实现静止和动态跟踪性能,设计了一种扰动观测器来观测永磁同步电机的负载扰动和负载变化。基于扰动观测器的补偿方法被用来解决标称系统的外部干扰。最重要的是,基于扰动观测器的方法可以有效地消除未知干扰,从而在一定程度上获得抗干扰能力和鲁棒性。
发明内容
本发明要解决的技术问题是:提供具有扰动及输出约束的永磁同步电机系统有限时间控制法,以解决现有技术中存在的技术问题。
本发明采取的技术方案为:具有扰动及输出约束的永磁同步电机系统有限时间控制法,该方法包括以下步骤:
(1)定义变量x1=θ,x2=ω,x3=iq,x4=id,对(d-q)坐标系下永磁同步电机的动力学模型重新构建,得到如下式:
Figure BDA0003490846180000021
其中x1受限于:
Figure BDA0003490846180000022
其中已知边界约束条件是时变函数
Figure BDA0003490846180000023
Figure BDA0003490846180000024
x1(t)表示输出变量,xd为参考信号,Δd为外部干扰,x=(x1,...,x4)T∈R4为状态向量,
Figure BDA0003490846180000025
a2=3np(Ld-Lq)/2,b1=-Rs/Lq,b2=-npLd/Lq
Figure BDA0003490846180000026
b4=1/Lq,c1=-Rs/Ld,c2=npLq/Ld,c3=1/Ld;ω为转子角速度,θ为转子角度,iq为q-轴电流,id为d-轴电流,uq为q-轴电压,ud为d-轴电压,J为转动惯量,B为摩擦系数,
Figure BDA0003490846180000027
为永磁通量,Rs为定子线圈电阻,np为极对数,Lq为q-轴线圈电感,Ld为d-轴线圈电感,TL为负载力矩;
设1:期望轨迹xd(t)及其i阶导数
Figure BDA0003490846180000028
是光滑有界的,状态约束函数
Figure BDA0003490846180000029
及其j阶导数
Figure BDA00034908461800000210
是连续有界的;
引理1:对于任意实变量p和q,以及常量ki>0,i=1,2,3,满足如下条件:
Figure BDA00034908461800000211
引理2:对于任意变量
Figure BDA00034908461800000218
和实数
Figure BDA00034908461800000212
得到:
Figure BDA00034908461800000213
定义1:对于状态向量Λ∈Rn,将非线性正定系统定义为
Figure BDA00034908461800000214
其中f(Λ)表示正定函数;对任意初始条件Λ(t0)=Λ0,存在实数ρ>0和稳定时间T(ρ,Λ0)<∞满足||Λ(t)||<ρ;对任意的t≥t0+T,存在一个平衡点λ=0满足
Figure BDA00034908461800000215
那么非线性系统在有限时间内是半全局稳定的;
引理3:对于非线性系统
Figure BDA00034908461800000216
如果存在实数
Figure BDA00034908461800000217
和光滑函数V(Λ)>0满足:
Figure BDA0003490846180000031
那么复杂系统
Figure BDA00034908461800000315
能够实现有限时间稳定性,并且稳定时间通过以下公式进行估计:
Figure BDA0003490846180000032
如果t≥T,则如下公式成立:
Figure BDA0003490846180000033
引理4:对于每个变量
Figure BDA0003490846180000034
其中φ12,且φ1和φ2为奇数,如下不等式成立如下:
Figure BDA0003490846180000035
其中实数β1>0,β2>0,
Figure BDA0003490846180000036
引理5:引入以下系统:
Figure BDA0003490846180000037
其中Δi,(i=1,…,n)是第i个状态变量,
Figure BDA00034908461800000316
为设计常量,并且扰动项f(t)满足|f(t)|≤L0;然后系统将在有限的时间内收敛到原点,从而非线性系统实现有限时间稳定性;
(2)A、使用径向基神经网络估计未知的非线性函数,对未知非线性函数在具有任意精度的闭合集合中估计,因此,得:
Figure BDA0003490846180000038
其中Y=[y1,y2,…,yn]T为输入变量,
Figure BDA0003490846180000039
为径向基神经网络的期望权重向量,l>1表示节点数,估计误差Ξ(Y)满足|Ξ(Y)|<ΞM,其中ΞM表示为不确定的有界变量;W(Y)=[w1(Y),w2(Y),...,wl(Y)]T是基函数向量,选择第i个通用高斯函数wi(Y)为
Figure BDA00034908461800000310
其中δi=[δi1,...,δim]表示第i个高斯基函数的中心,
Figure BDA00034908461800000311
表示第i个高斯函数的宽度;
考虑期望权重向量
Figure BDA00034908461800000312
为:
Figure BDA00034908461800000313
其中
Figure BDA00034908461800000314
表示更新权重向量;
利用2范数来估计权重减轻神经网络的计算负担,因此,得到:
Figure BDA0003490846180000041
其中
Figure BDA0003490846180000042
表示第i个不确定变量,||·||代表·的2范数;
(3)设计基于有限时间扰动观测器的动态面控制
非线性坐标转换函数被用来将具有输出约束系统(2)转换为非约束系统,然后输出变量被约束在具有不对称时变边界的紧集中;
定义2:设计非线性转换函数为:
Figure BDA0003490846180000043
其中s1=x1-xd表示跟踪误差,ζ1表示转换误差,
Figure BDA0003490846180000044
是连续时变函数,假设存在两个正常量C0和C1,则变量
Figure BDA0003490846180000045
满足以下条件:
Figure BDA0003490846180000046
已知
Figure BDA0003490846180000047
Figure BDA0003490846180000048
是有界的,由此推断,从公式(15)可以推断出ζ1的有界性取决于跟踪误差s1;对于每个初始条件
Figure BDA0003490846180000049
当ζ1在t∈[0,+∞)内有界时,推断出s1(t)满足如下条件:
Figure BDA00034908461800000410
为了简化下文的表达式,分别将
Figure BDA00034908461800000411
Figure BDA00034908461800000412
缩写为
Figure BDA00034908461800000413
Figure BDA00034908461800000414
对ζ1求导,得到:
Figure BDA00034908461800000415
其中
Figure BDA00034908461800000416
Figure BDA00034908461800000417
基于(3)和(16),得到输出无约束子系统,如下式所示:
Figure BDA00034908461800000418
基于(15),得到无约束变量ζ1的值域是所有实数集,结合η11和假设1,可以推断出
Figure BDA00034908461800000419
此外,结合洛必达法则,得到:
Figure BDA00034908461800000420
当s1不受约束时,转换误差ζ1接近跟踪误差s1
B、基于结合辅助变量d,D,设计有限时间扰动观测器为:
Figure BDA00034908461800000421
其中
Figure BDA00034908461800000422
Figure BDA00034908461800000423
分别代表Δd,x2和D的估计值,定义
Figure BDA00034908461800000424
λ1和λ2并表示扰动观测器的正设计系数,κ1为适当的正常数;
定义
Figure BDA00034908461800000425
然后对
Figure BDA00034908461800000426
Figure BDA00034908461800000427
对时间t求导,得到:
Figure BDA0003490846180000051
假设不匹配扰动的变化可以忽略不计,则Δd有界,并且估计值
Figure BDA0003490846180000052
Figure BDA0003490846180000053
将分别在有限的时间内收敛到实际值Δd,x2和D;
C、设计基于扰动观测器的自适应有限时间控制器:
对跟踪误差采用非线性坐标变换:
Figure BDA0003490846180000054
其中
Figure BDA00034908461800000517
是下列第i个一阶滤波器的输出:
Figure BDA0003490846180000055
其中εi代表第i个滤波器的时间常数,第i个虚拟控制信号
Figure BDA00034908461800000518
将在之后设计;
类似地,第i个滤波器误差zi,i=2,3被定义为:
Figure BDA00034908461800000519
将(2)和(19)代入到(23)中,对(23)中的
Figure BDA00034908461800000516
求导:
Figure BDA0003490846180000056
定义第i个变量估计误差
Figure BDA0003490846180000057
为:
Figure BDA0003490846180000058
其中第i个变量
Figure BDA0003490846180000059
的估计值为
Figure BDA00034908461800000510
基于扰动观测器的自适应有限时间控制器的设计步骤如下:
步骤1、选择的第一个子Lyapunov函数V1为:
Figure BDA00034908461800000511
其中设计常数τ1>0;
结合(27),对(28)中的V1求导:
Figure BDA00034908461800000512
将(26)代入(29),得到:
Figure BDA00034908461800000513
定义F1(X1)为:
Figure BDA00034908461800000514
其中
Figure BDA00034908461800000515
采用径向基神经网络估计未知函数F1(X1):
Figure BDA0003490846180000061
其中设计常数ΞM>0.
因此,(30)变为:
Figure BDA0003490846180000062
根据杨氏不等式,得到:
Figure BDA0003490846180000063
其中设计常数p1>0;
将(34)代入(33)得到:
Figure BDA0003490846180000064
设计虚构控制律
Figure BDA00034908461800000616
和自适应律
Figure BDA0003490846180000065
为:
Figure BDA0003490846180000066
其中k11,k12,
Figure BDA00034908461800000617
表示正常数,α=α12和奇数α12满足如下条件:0<α12
将(36)代入到(35),得到
Figure BDA0003490846180000067
结合(24)-(26),(27)和(36),对z2求导:
Figure BDA0003490846180000068
其中
Figure BDA0003490846180000069
为一个连续函数;
推导出
Figure BDA00034908461800000610
Figure BDA00034908461800000611
满足在紧集的预设初始条件下服从的最大值,因此:
Figure BDA00034908461800000612
其中
Figure BDA00034908461800000613
利用杨氏不等式,得到:
Figure BDA00034908461800000614
将(40)代入到(37)中,得到:
Figure BDA00034908461800000615
步骤2、选择第二个子Lyapunov函数V2为:
Figure BDA0003490846180000071
其中τ2>0为已知常量;
结合(27),对V2进行求导,得到:
Figure BDA0003490846180000072
将(26)和(41)代入到(43)中,得到:
Figure BDA0003490846180000073
基于引理5,假设
Figure BDA0003490846180000074
其中q1表示正常数;然后将(44)简化为:
Figure BDA0003490846180000075
构造F2(X2)
Figure BDA0003490846180000076
其中
Figure BDA00034908461800000711
结合(46),然后(45)重新表述为:
Figure BDA0003490846180000077
与(32)类似,得到(47)中的F2(X2)未知,因此,提出一个径向基神经网络来估计F2(X2):
Figure BDA0003490846180000078
因此,(47)被重新表达为:
Figure BDA0003490846180000079
类似于(34),得到如下不等式:
Figure BDA00034908461800000710
其中设计常数p2>0;
将(50)代入(49),得到:
Figure BDA0003490846180000081
类似于(36),选择虚拟控制器
Figure BDA00034908461800000812
和自适应律
Figure BDA0003490846180000082
为:
Figure BDA0003490846180000083
其中设计常数k21,k22,
Figure BDA00034908461800000813
为正数;
将(52)代入到(51)中,得到:
Figure BDA0003490846180000084
类似于(40),得到:
Figure BDA0003490846180000085
其中函数
Figure BDA0003490846180000086
将(54)代入(53),得到
Figure BDA0003490846180000087
步骤3、选择第三个子Lyapunov函数V3为:
Figure BDA0003490846180000088
其中已知常量t3>0;
根据(27),将V3对时间t求导:
Figure BDA0003490846180000089
将(26)和(55)代入(57),得到:
Figure BDA00034908461800000810
设计F3(X3)为:
Figure BDA00034908461800000811
其中
Figure BDA00034908461800000814
因此,(58)重新表示为:
Figure BDA0003490846180000091
同样,F3(X3)是未知的,然后,应用径向基神经网络来估计
Figure BDA0003490846180000092
Figure BDA0003490846180000093
类似于(34),得到:
Figure BDA0003490846180000094
其中设计参数p3>0;
因此,(60)被重新表达为:
Figure BDA0003490846180000095
设计实际控制器uq和自适应律
Figure BDA0003490846180000096
为:
Figure BDA0003490846180000097
其中设计常数k31,k32,
Figure BDA00034908461800000913
为正数;
将(64)代入(63),得到:
Figure BDA0003490846180000098
步骤4、选择第四个子Lyapunov函数V4为:
Figure BDA0003490846180000099
其中设计常数τ4>0;
结合(27),对(66)中的V4求导,得到:
Figure BDA00034908461800000910
结合(26),并将(65)代入(67)得到:
Figure BDA00034908461800000911
构造F4(X4)为:
Figure BDA00034908461800000912
其中
Figure BDA0003490846180000109
然后,(68)被简化为:
Figure BDA0003490846180000101
由于F4(X4)是一个不确定的函数,应用径向基神经网络来估计它,得到:
Figure BDA0003490846180000102
类似于(34),得到:
Figure BDA0003490846180000103
其中设计常数p4>0;
将(72)代入(70),得到:
Figure BDA0003490846180000104
设计实际控制器ud和虚拟控制律
Figure BDA0003490846180000105
Figure BDA0003490846180000106
其中k41,k42,
Figure BDA0003490846180000108
为正常数;
将(74)代入(73),得到:
Figure BDA0003490846180000107
本发明的有益效果:与现有技术相比,本发明的效果如下:
1)利用合适的非线性转换函数,将具有不匹配外部干扰的输出约束永磁同步电机系统转换为新颖的无约束系统,与通常使用的基于分段障碍Lyapunov函数的反步法来解决非对称约束相比,这种设计的两个优点是,有约束和无约束系统都可以在不变的控制器下直接方便地求解,并且消除了对连续虚拟控制器的要求;
2)与基于扰动观测器渐近反步法不同,本发明提出了一种有限时间扰动观测器来估计无约束系统的不匹配外部干扰,从而增强了控制系统的抗干扰能力和动态响应能力;
3)与现有技术中的渐近稳定和有限时间稳定下的反步法不同,本发明的控制方案采用快速有限时间稳定理论将坐标变换和有限时间扰动观测器集成到动态面控制框架中,这项工作具有设计要求少、用户在实际应用中容易使用等优点,使其不受干扰,收敛速度快,对抑制输出约束和不匹配的外部干扰具有出色的鲁棒性。
附图说明
图1为永磁同步电机系统控制原理结构示意图;
图2为有限时间动态面控制、比例积分微分控制和神经动态面控制方案中输出信号x1和期望信号xd的曲线图,其中(a)为情况1中的系统,(b)为情况2中的系统;
图3为有限时间动态面控制、比例积分微分控制和神经动态面控制方案中跟踪误差s1的响应曲线图,其中(a)为情况1中的系统;(b)为情况2中的系统;
图4为状态变量x2的响应曲线图;其中(a)为情况1中的系统。(b)为情况2中的系统;
图5为状态变量id和iq的响应曲线图;(a)为情况1中的系统,(b)为情况2中的系统;
图6为控制器ud和uq的轨迹图;其中(a)为情况1中的系统。(b)为情况2中的系统。
具体实施方式
下面结合具体的实施例对本发明进行进一步介绍。
实施例1:如图1-6所示,具有扰动及输出约束的永磁同步电机系统有限时间控制法,包括以下步骤:
A系统说明(d-q)坐标系下永磁同步电机的动力学模型可以表述为:
Figure BDA0003490846180000111
其中,ω为转子角速度(rad/s),θ为转子角度(°),iq为q-轴电流(A),id为d-轴电流(A),uq为q-轴电压(V),ud为d-轴电压(V),J为转动惯量(kg·m2),B为摩擦系数(N/(rad/s)),
Figure BDA0003490846180000112
为永磁通量(Wb),Rs为定子线圈电阻(Ω),np为极对数,Lq为q-轴线圈电感(H),Ld为d-轴线圈电感(H),TL为负载力矩(N·m);
定义变量x1=θ,x2=ω,x3=iq,x4=id,对(d-q)坐标系下永磁同步电机的动力学模型重新构建,得到如下式:
Figure BDA0003490846180000113
其中x1受限于:
Figure BDA0003490846180000114
其中已知边界约束条件是时变函数
Figure BDA0003490846180000115
Figure BDA0003490846180000116
x1(t)表示输出变量,xd为参考信号,Δd为外部干扰,x=(x1,...,x4)T∈R4为状态向量,
Figure BDA0003490846180000117
a2=3np(Ld-Lq)/2,b1=-Rs/Lq,b2=-npLd/Lq
Figure BDA0003490846180000118
b4=1/Lq,c1=-Rs/Ld,c2=npLq/Ld,c3=1/Ld;ω为转子角速度,θ为转子角度,iq为q-轴电流,id为d-轴电流,uq为q-轴电压,ud为d-轴电压,J为转动惯量,B为摩擦系数,
Figure BDA0003490846180000119
为永磁通量,Rs为定子线圈电阻,np为极对数,Lq为q-轴线圈电感,Ld为d-轴线圈电感,TL为负载力矩;
首先处理了永磁同步电机系统中提出的不匹配干扰和方程非对称输出约束。与常数对称约束相比,方程采用的非对称时变约束上限和下限更具有一般性。
本发明的目标是设计一种基于扰动观测器的自适应有限时间动态面控制方案,该方案具有有限时间特性:
a)跟踪误差在有限的时间内收敛到原点的小邻域内,并且闭环系统中的所有信号都有界。
b)输出信号满足中的约束条件。
为了实现这些目标,给出了下面的假设和引理:
设1:期望轨迹xd(t)及其i阶导数
Figure BDA0003490846180000121
是光滑有界的,状态约束函数
Figure BDA0003490846180000122
及其j阶导数
Figure BDA0003490846180000123
是连续有界的;
引理1:对于任意实变量p和q,以及常量ki>0,i=1,2,3,满足如下条件:
Figure BDA0003490846180000124
引理2:对于任意变量
Figure BDA00034908461800001217
和实数
Figure BDA0003490846180000125
得到:
Figure BDA0003490846180000126
定义1:对于状态向量Λ∈Rn,将非线性正定系统定义为
Figure BDA0003490846180000127
其中f(Λ)表示正定函数;对任意初始条件Λ(t0)=Λ0,存在实数ρ>0和稳定时间T(ρ,Λ0)<∞满足||Λ(t)||<ρ;对任意的t≥t0+T,存在一个平衡点λ=0满足
Figure BDA0003490846180000128
那么非线性系统在有限时间内是半全局稳定的;
引理3:对于非线性系统
Figure BDA0003490846180000129
如果存在实数
Figure BDA00034908461800001210
和光滑函数V(Λ)>0满足:
Figure BDA00034908461800001211
那么复杂系统
Figure BDA00034908461800001212
能够实现有限时间稳定性,并且稳定时间通过以下公式进行估计:
Figure BDA00034908461800001213
如果t≥T,则如下公式成立:
Figure BDA00034908461800001214
引理4:对于每个变量
Figure BDA00034908461800001215
其中φ12,且φ1和φ2为奇数,如下不等式成立如下:
Figure BDA00034908461800001216
其中实数β1>0,β2>0,
Figure BDA0003490846180000131
引理5:引入以下系统:
Figure BDA0003490846180000132
其中Δi,(i=1,…,n)是第i个状态变量,
Figure BDA00034908461800001319
为设计常量,并且扰动项f(t)满足|f(t)|≤L0;然后系统将在有限的时间内收敛到原点,从而非线性系统实现有限时间稳定性;
为了简化表达,在不发生混淆的前提下,之后的陈述中有些函数将会被简写。
(2)神经网络系统和函数估计
A、使用径向基神经网络估计未知的非线性函数,对未知非线性函数在具有任意精度的闭合集合中估计,因此,得:
Figure BDA0003490846180000133
其中Y=[y1,y2,…,yn]T为输入变量,
Figure BDA0003490846180000134
为径向基神经网络的期望权重向量,l>1表示节点数,估计误差Ξ(Y)满足|Ξ(Y)|<ΞM,其中ΞM表示为不确定的有界变量;W(Y)=[w1(Y),w2(Y),...,wl(Y)]T是基函数向量,选择第i个通用高斯函数wi(Y)为
Figure BDA0003490846180000135
其中δi=[δi1,...,δim]表示第i个高斯基函数的中心,
Figure BDA0003490846180000136
表示第i个高斯函数的宽度;
考虑期望权重向量
Figure BDA0003490846180000137
为:
Figure BDA0003490846180000138
其中
Figure BDA0003490846180000139
表示更新权重向量;
利用2范数来估计权重减轻神经网络的计算负担,因此,得到:
Figure BDA00034908461800001310
其中
Figure BDA00034908461800001311
表示不确定变量,||·||代表·的2范数;
(3)设计基于有限时间扰动观测器的动态面控制
非线性坐标转换函数被用来将具有输出约束系统(2)转换为非约束系统,然后输出变量被约束在具有不对称时变边界的紧集中;
定义2:设计非线性转换函数为:
Figure BDA00034908461800001312
其中s1=x1-xd表示跟踪误差,
Figure BDA00034908461800001313
表示转换误差,
Figure BDA00034908461800001314
是连续时变函数,假设存在两个正常量C0和C1,则变量
Figure BDA00034908461800001315
满足以下条件:
Figure BDA00034908461800001316
已知
Figure BDA00034908461800001317
Figure BDA00034908461800001318
是有界的,由此推断,从公式(15)可以推断出ζ1的有界性取决于跟踪误差s1;对于每个初始条件
Figure BDA0003490846180000141
当ζ1在t∈[0,+∞)内有界时,推断出s1(t)满足如下条件:
Figure BDA0003490846180000142
为了简化表达,分别将
Figure BDA0003490846180000143
Figure BDA0003490846180000144
缩写为
Figure BDA0003490846180000145
Figure BDA0003490846180000146
对ζ1求导,得到:
Figure BDA0003490846180000147
其中
Figure BDA0003490846180000148
Figure BDA0003490846180000149
基于(3)和(16),得到输出无约束子系统,如下式所示:
Figure BDA00034908461800001410
基于(15),得到无约束变量ζ1的值域是所有实数集,结合η11和假设1,推断出
Figure BDA00034908461800001411
此外,结合洛必达法则,得到:
Figure BDA00034908461800001412
当s1不受约束时,转换误差ζ1接近跟踪误差s1;因此,可以得出结论,公式中的非线性转换函数可以同时解决有约束和非约束系统,这补偿了通过采用
Figure BDA00034908461800001425
来处理无约束系统而造成的特殊转换的不足。综上所述,非线性坐标变换函数不仅可以解决同一控制器下有状态约束和没有状态约束的系统,而且能够成功地规避基于分段不对称障碍Lyapunov方法的虚拟控制器的连续要求。因此,本发明提出的坐标变换方法更适合于设计与有限时间稳定法相结合的期望的反步控制器。
B、有限时间扰动观测器:
基于结合辅助变量d,D,设计有限时间扰动观测器为:
Figure BDA00034908461800001413
其中
Figure BDA00034908461800001414
Figure BDA00034908461800001415
分别代表Δd,x2和D的估计值,定义
Figure BDA00034908461800001416
λ1和λ2并表示扰动观测器的正设计系数,κ1为适当的正常数;
定义
Figure BDA00034908461800001417
然后对
Figure BDA00034908461800001418
Figure BDA00034908461800001419
对时间t求导,得到:
Figure BDA00034908461800001420
假设不匹配扰动的变化可以忽略不计,则Δd有界,并且估计值
Figure BDA00034908461800001421
Figure BDA00034908461800001422
将分别在有限的时间内收敛到实际值Δd,x2和D;
很容易地知道,通过选择合适的参数λ12和κ1,估计误差
Figure BDA00034908461800001423
Figure BDA00034908461800001424
可以在有限的时间内收敛到原点的一个小邻域内。同时,有限时间扰动观测器可以处理不匹配的扰动,并且避免了滑模控制带来的颤振现象。因此,可以通过选择适当的值来增强永磁同步电机系统的鲁棒性和收敛性。
C、设计基于扰动观测器的自适应有限时间控制器:
对跟踪误差采用非线性坐标变换:
Figure BDA0003490846180000151
其中
Figure BDA00034908461800001514
是下列第i个一阶滤波器的输出:
Figure BDA0003490846180000152
其中εi代表第i个滤波器的时间常数,第i个虚拟控制信号
Figure BDA00034908461800001515
将在之后设计;
使用一阶滤波器来估计
Figure BDA00034908461800001516
的时间导数,处理“复杂性的爆炸”。因此,通过应用该方法可以实现对系统(2)的基于扰动观测器的自适应有限时间动态面控制,有效地减轻计算负担。
类似地,第i个滤波器误差zi,i=2,3被定义为:
Figure BDA00034908461800001517
将(2)和(19)代入到(23)中,对(23)中的第i个误差转换变量
Figure BDA00034908461800001513
求导:
Figure BDA0003490846180000153
定义第i个变量的估计误差
Figure BDA0003490846180000154
为:
Figure BDA0003490846180000155
其中第i个变量
Figure BDA0003490846180000156
的估计值为
Figure BDA0003490846180000157
基于扰动观测器的自适应有限时间控制器的设计步骤如下:
步骤1、选择的第一个子Lyapunov函数V1为:
Figure BDA0003490846180000158
其中设计常数τ1>0;
结合(27),对(28)中的V1求导:
Figure BDA0003490846180000159
将(26)代入(29),得到:
Figure BDA00034908461800001510
定义F1(X1)为:
Figure BDA00034908461800001511
其中
Figure BDA00034908461800001512
采用径向基神经网络估计未知函数F1(X1):
Figure BDA0003490846180000161
其中设计常数ΞM>0.
因此,(30)变为:
Figure BDA0003490846180000162
根据杨氏不等式,得到:
Figure BDA0003490846180000163
其中设计常数p1>0;
将(34)代入(33)得到:
Figure BDA0003490846180000164
设计虚构控制律
Figure BDA00034908461800001616
和自适应律
Figure BDA0003490846180000165
为:
Figure BDA0003490846180000166
其中k11,k12,
Figure BDA00034908461800001617
表示正常数,α=α12和奇数α12满足如下条件:0<α12
将(36)代入到(35),得到
Figure BDA0003490846180000167
结合(24)-(26),(27)和(36),对z2求导:
Figure BDA0003490846180000168
其中
Figure BDA0003490846180000169
为一个连续函数;
推导出
Figure BDA00034908461800001610
Figure BDA00034908461800001611
满足在紧集的预设初始条件下服从的最大值,因此:
Figure BDA00034908461800001612
其中
Figure BDA00034908461800001613
利用杨氏不等式,得到:
Figure BDA00034908461800001614
将(40)代入到(37)中,得到:
Figure BDA00034908461800001615
步骤2、选择第二个子Lyapunov函数V2为:
Figure BDA0003490846180000171
其中τ2>0为已知常量;
结合(27),对V2进行求导,得到:
Figure BDA0003490846180000172
将(26)和(41)代入到(43)中,得到:
Figure BDA0003490846180000173
基于引理5,假设
Figure BDA0003490846180000174
其中q1表示正常数;然后将(44)简化为:
Figure BDA0003490846180000175
构造F2(X2)
Figure BDA0003490846180000176
其中
Figure BDA00034908461800001711
结合(46),然后(45)重新表述为:
Figure BDA0003490846180000177
与(32)类似,得到(47)中的F2(X2)未知,因此,提出一个径向基神经网络来估计F2(X2):
Figure BDA0003490846180000178
因此,(47)被重新表达为:
Figure BDA0003490846180000179
类似于(34),得到如下不等式:
Figure BDA00034908461800001710
其中设计常数p2>0;
将(50)代入(49),得到:
Figure BDA0003490846180000181
类似于(36),选择虚拟控制器
Figure BDA00034908461800001812
和自适应律
Figure BDA0003490846180000182
为:
Figure BDA0003490846180000183
其中设计常数k21,k22,
Figure BDA00034908461800001813
为正数;
通过有限时间扰动观测器(21)来估计不匹配的扰动,可以补偿虚拟控制器(52)中的扰动Δd,从而增强复杂系统的鲁棒性。
将(52)代入到(51)中,得到:
Figure BDA0003490846180000184
类似于(40),得到:
Figure BDA0003490846180000185
其中函数
Figure BDA0003490846180000186
将(54)代入(53),得到
Figure BDA0003490846180000187
步骤3、选择第三个子Lyapunov函数V3为:
Figure BDA0003490846180000188
其中已知常量τ3>0;
根据(27),将V3对时间t求导:
Figure BDA0003490846180000189
将(26)和(55)代入(57),得到:
Figure BDA00034908461800001810
设计F3(X3)为:
Figure BDA00034908461800001811
其中
Figure BDA00034908461800001912
因此,(58)重新表示为:
Figure BDA0003490846180000191
同样,F3(X3)是未知的,然后,应用径向基神经网络来估计
Figure BDA0003490846180000192
Figure BDA0003490846180000193
类似于(34),得到:
Figure BDA0003490846180000194
其中设计参数p3>0;
因此,(60)被重新表达为:
Figure BDA0003490846180000195
设计实际控制器uq和自适应律
Figure BDA0003490846180000196
为:
Figure BDA0003490846180000197
其中设计常数k31,k32,
Figure BDA00034908461800001913
为正数;
将(64)代入(63),得到:
Figure BDA0003490846180000198
步骤4、选择第四个子Lyapunov函数V4为:
Figure BDA0003490846180000199
其中设计常数τ4>0;
结合(27),对(66)中的V4求导,得到:
Figure BDA00034908461800001910
结合(26),并将(65)代入(67)得到:
Figure BDA00034908461800001911
构造F4(X4)为:
Figure BDA0003490846180000208
其中
Figure BDA0003490846180000209
然后,(68)被简化为:
Figure BDA0003490846180000201
由于F4(X4)是一个不确定的函数,应用径向基神经网络来估计它,得到:
Figure BDA0003490846180000202
类似于(34),得到:
Figure BDA0003490846180000203
其中设计常数p4>0;
将(72)代入(70),得到:
Figure BDA0003490846180000204
设计实际控制器ud和虚拟控制律
Figure BDA0003490846180000205
Figure BDA0003490846180000206
其中k41,k42,
Figure BDA00034908461800002010
为正常数;
将(74)代入(73),得到:
Figure BDA0003490846180000207
目前为止,控制器设计的整个过程已经完成。详细的控制结构演示如图1所示。稳定性分析:对于任何给定e>0,将紧集定义为:
Figure BDA0003490846180000211
定理1:基于假设1,针对永磁同步电机系统(2)所提出的基于扰动观测器的自适应有限时间动态面控制方法包括控制律
Figure BDA00034908461800002112
uq,ud和自适应律
Figure BDA0003490846180000212
如果满足初始条件Ωi,i=1,...,4,
Figure BDA0003490846180000213
和xd∈(-d,d),那么所有控制目标就会实现。
证明:选择整体Lyapunov函数为:
Figure BDA0003490846180000214
对(77)中的V求导,可以得到:
Figure BDA0003490846180000215
根据(9)、(27)和杨氏不等式,可以得到:
Figure BDA0003490846180000216
其中
Figure BDA0003490846180000217
因此,(78)可以被重写为:
Figure BDA0003490846180000218
此外,分别设计引理2中的参数为p=1,μ1=1-α,μ2=α,μ3=αα/(1-α),可以得到|q|α≤(1-α)g+|q|。分别设计参数q为
Figure BDA0003490846180000219
Figure BDA00034908461800002110
因此可以得到:
Figure BDA00034908461800002111
将(81)代入(80),可以得到:
Figure BDA0003490846180000221
其中
Figure BDA00034908461800002220
Figure BDA0003490846180000222
定义
Figure BDA0003490846180000223
基于引理2,(82)可以被整理为:
Figure BDA0003490846180000224
基于引理3,可以推导出以下特征:
(Ⅰ)非线性系统(2)是有限时间稳定的;
(II)对于任意
Figure BDA00034908461800002218
存在一个有限的时间T和所有t≥T,可以得到:
Figure BDA0003490846180000225
其中有限时间T通过如下公式得到:
Figure BDA0003490846180000226
从(85)中,可以推断出:
Figure BDA00034908461800002217
从(86)中可以推断出第i个误差转换变量
Figure BDA00034908461800002219
是有界的。与此类似,可以确保
Figure BDA0003490846180000227
Figure BDA0003490846180000228
的有界性。因此,根据(27),可以知道
Figure BDA0003490846180000229
有界。根据
Figure BDA00034908461800002210
可以得到
Figure BDA00034908461800002211
是有界的。对于
Figure BDA00034908461800002212
推导出
Figure BDA00034908461800002213
Figure BDA00034908461800002214
根据(87),可以推断出s1有界的。因为xd严格有界,根据s1=x1-xd可以得到x1有界。因为s1有界,可以得到η1有界,因此可以保证的有界性。根据(36),可以得出
Figure BDA00034908461800002221
有界。因此,从(24)可以得到
Figure BDA00034908461800002222
Figure BDA00034908461800002215
有界。基于此,结合(23),可以确保x2的有界性。类似地,可以得出结论:
Figure BDA00034908461800002216
有界。综上所述,可以得出闭环系统中所有信号是有界的结论。
特别是,从(85)中可以得出通过选择合适的设计参数
Figure BDA0003490846180000231
可以实现变量x1在有限时间内跟踪参考信号xd
此外,结合(15)和(23),当且仅当
Figure BDA0003490846180000232
ζ1的数值会接近无穷大,可以得到:对任意的t>0,有
Figure BDA0003490846180000233
对于初始条件,满足
Figure BDA0003490846180000234
因为s1=x1-xd,对s1的任意初始值,满足
Figure BDA0003490846180000235
对任意t>0,满足
Figure BDA0003490846180000236
基于上述讨论,可以证明所设计的控制器展示出了对被控系统优越的瞬态性能和稳定性。
跟踪误差s1的值可以反映控制效果的好坏。从(87)可以看出s1的值取决于变量σ0和b0的数值,尤其当σ0减少和b0增加时,s1的值接近无限大。而σ0和b0的值取决于变量
Figure BDA0003490846180000237
和设计参数
Figure BDA0003490846180000238
进一步,通过增加
Figure BDA00034908461800002313
和减少
Figure BDA0003490846180000239
的数值,s1的值可以被约束到预设的边界中。因此,应考虑与控制工作相关的系统性能权衡。
与现有技术相比,本发明的控制策略考虑了有限时间扰动观测器和动态面控制来消除外部干扰并确保输出约束不被违反。此外,设计的控制律以及自适应律中包含(2α-1)阶的
Figure BDA00034908461800002310
实现了对永磁同步电机系统的有限时间稳定控制。因此,本发明设计的控制方法在现实中更具功能性和普遍性。
为了验证本发明的效果,进行仿真实验和进行结果比较
1.1控制器的设计在本部分中,提供了两个模拟案例,用于验证对系统(2)设计的控制策略的有效性。为了更好地分析本发明设计的控制方案的普遍性和优缺点,针对两种情况设计不同的时变边界函数,并在每种方案中将比例积分微分控制和神经动态面控制作为对比方案进行仿真实验,同时考虑了存在外部干扰和没有外部干扰(Δd=0)的情况。
选择永磁同步电机参数作为:J=0.003798Kg·m2,B=0.001158N·m/(rad/s),TL=1.5,
Figure BDA00034908461800002311
Ld=0.00285H,np=3,Lq=0.00315H,Rs=0.68Ω。分别选择参考信号和扰动函数为xd=0.49(sin(t)+sin(2t))和Δd=20x2sin(2t)。
(1)神经动态面控制:考虑神经动态面控制的控制器设计,选择以下变量:
Figure BDA00034908461800002312
Figure BDA0003490846180000241
其中神经动态面控制方法的设计参数为k11=k21=k31=k41=80,
Figure BDA0003490846180000242
Figure BDA00034908461800002412
ε2=ε3=0.01。
(2)比例积分微分:考虑ud=0,选择比例积分微分控制器为:
Figure BDA0003490846180000243
其中设计参数kp=20,ki=0.05,kd=1.5。
同时,将本发明的方案与神经动态面控制和比例积分微分方法进行比较,并以如下3个定量指标作为衡量标准,详细结果列于表2中。
a)误差绝对值的积分:
Figure BDA0003490846180000244
b)时间和误差的绝对值乘积的积分:
Figure BDA0003490846180000245
c)误差平方的积分:
Figure BDA0003490846180000246
基于上述讨论,设计如下两种情况:
情况1:将时变边界设计为
Figure BDA0003490846180000247
并选择状态变量的起始值为xi(0)=0,i=1,...,4。每个径向基神经网络由11个节点组成,中心位于区间[-11,11]内,宽度为10。控制器的设计参数选择如下:
Figure BDA00034908461800002413
Figure BDA0003490846180000249
ε2=0.1,ε3=0.01,τ1=τ2=τ3=0.55,τ4=60,ki1=20,ki2=30,p1=0.06,p2=p3=p4=0.15,
Figure BDA00034908461800002414
1≤i≤4,α=17/19。
情况2:设计时变边界为
Figure BDA00034908461800002410
选择控制器的设计参数为
Figure BDA00034908461800002415
ε2=0.1,ε3=0.01,τ1=τ2=τ3=0.65,τ4=60,ki1=40,ki2=65,p1=0.06,p2=p3=p4=0.15,
Figure BDA0003490846180000252
α=97/101。
1.2仿真实验和比较结果
图2-6展示了直观的仿真实验和比较结果。每张图片包含两个子图片,子图片(a)和(b)分别代表情况1和情况2边界下三种控制策略(有限时间动态面控制,比例积分微分和神经动态面控制)的性能。图2为输出信号x1和期望轨迹xd的曲线。图3为跟踪误差s1的响应曲线。从图2-3可以看出,与比例积分微分和神经动态面控制方法相比,设计的控制器在外部扰动下输出信号x1对期望轨迹xd表现出优异的跟踪性能,并且跟踪误差不超出时变的边界。图4展示了有限时间动态面控制方法中状态变量x2的响应曲线。图5为状态变量id和iq的响应,图6为控制器ud和uq的轨迹。从图2-6可以明显看出,与扰动的神经动态面控制和比例积分微分控制相比,本发明设计的有限时间动态面控制策略可以在扰动情况下确保状态变量具有的极快响应能力和有限时间收敛速率。
为了简化表达式,给出了以下缩写:C-D和W-D分别表示包含干扰和无干扰的系统。
表2性能指标的比较结果
Figure BDA0003490846180000251
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内,因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (1)

1.具有扰动及输出约束的永磁同步电机系统有限时间控制法,其特征在于:该方法包括以下步骤:
(1)定义变量x1=θ,x2=ω,x3=iq,x4=id,对(d-q)坐标系下永磁同步电机的动力学模型重新构建,得到如下式:
Figure FDA0003490846170000011
其中x1受限于:
Figure FDA0003490846170000012
其中已知边界约束条件是时变函数
Figure FDA0003490846170000013
Figure FDA0003490846170000014
x1(t)表示输出变量,xd为参考信号,Δd为外部干扰,x=(x1,...x4)T∈R4为状态向量,
Figure FDA0003490846170000015
a2=3np(Ld-Lq)/2,b1=-Rs/Lq,b2=-npLd/Lq
Figure FDA0003490846170000016
b4=1/Lq,c1=-Rs/Ld,c2=npLq/Ld,c3=1/Ld;ω为转子角速度,θ为转子角度,iq为q-轴电流,id为d-轴电流,uq为q-轴电压,ud为d-轴电压,J为转动惯量,B为摩擦系数,
Figure FDA0003490846170000017
为永磁通量,Rs为定子线圈电阻,np为极对数,Lq为q-轴线圈电感,Ld为d-轴线圈电感,TL为负载力矩;
设1:期望轨迹xd(t)及其i阶导数
Figure FDA0003490846170000018
是光滑有界的,状态约束函数
Figure FDA0003490846170000019
及其j阶导数
Figure FDA00034908461700000110
是连续有界的;
引理1:对于任意实变量p和q,以及常量ki>0,i=1,2,3,满足如下条件:
Figure FDA00034908461700000111
引理2:对于任意变量
Figure FDA00034908461700000112
m和实数
Figure FDA00034908461700000120
得到:
Figure FDA00034908461700000113
定义1:对于状态向量Λ∈Rn,将非线性正定系统定义为
Figure FDA00034908461700000114
其中f(Λ)表示正定函数;对任意初始条件Λ(t0)=Λ0,存在实数ρ>0和稳定时间T(ρ,Λ0)<∞满足||Λ(t)||<ρ;对任意的t≥t0+T,存在一个平衡点Λ=0满足
Figure FDA00034908461700000115
那么非线性系统在有限时间内是半全局稳定的;
引理3:对于非线性系统
Figure FDA00034908461700000116
如果存在实数a>0,
Figure FDA00034908461700000117
b>0,σ>0,
Figure FDA00034908461700000118
和光滑函数V(Λ)>0满足:
Figure FDA00034908461700000119
那么复杂系统
Figure FDA0003490846170000021
能够实现有限时间稳定性,并且稳定时间通过以下公式进行估计:
Figure FDA0003490846170000022
如果t ≥ T,则如下公式成立:
Figure FDA0003490846170000023
引理4:对于每个变量
Figure FDA0003490846170000024
φ=φ12其中φ12,且φ1和φ2为奇数,如下不等式成立如下:
Figure FDA0003490846170000025
其中实数β1>0,β2>0,
Figure FDA0003490846170000026
γ2=(2φ-1-2(1+φ)(φ-1))/(1+φ)>0;
引理5:引入以下系统:
Figure FDA0003490846170000027
其中Δi,(i=1,…,n)是第i个状态变量,L0i,
Figure FDA0003490846170000028
为设计常量,并且扰动项f(t)满足|f(t)|≤L0;然后系统将在有限的时间内收敛到原点,从而非线性系统实现有限时间稳定性;
(2)A、使用径向基神经网络估计未知的非线性函数,对未知非线性函数在具有任意精度的闭合集合中估计,因此,得:
Figure FDA0003490846170000029
其中Y=[y1,y2,…,yn]T为输入变量,
Figure FDA00034908461700000210
为径向基神经网络的期望权重向量,l>1表示节点数,估计误差Ξ(Y)满足|Ξ(Y)|<ΞM,其中ΞM表示为不确定的有界变量;W(Y)=[w1(Y),w2(Y),...,wl(Y)]T是基函数向量,选择第i个通用高斯函数wi(Y)为
Figure FDA00034908461700000211
其中δi=[δi1,...,δim]表示第i个高斯基函数的中心,
Figure FDA00034908461700000212
表示高斯函数的宽度;
考虑期望权重向量
Figure FDA00034908461700000213
为:
Figure FDA00034908461700000214
其中
Figure FDA00034908461700000215
表示更新权重向量;
利用2范数来估计权重减轻神经网络的计算负担,因此,得到:
Figure FDA0003490846170000031
其中
Figure FDA0003490846170000032
Figure FDA0003490846170000033
表示第i个不确定变量,||·||代表·的2范数;
(3)设计基于有限时间扰动观测器的动态面控制
非线性坐标转换函数被用来将具有输出约束系统(2)转换为非约束系统,然后输出变量被约束在具有不对称时变边界的紧集中;
定义2:设计非线性转换函数为:
Figure FDA0003490846170000034
其中s1=x1-xd表示跟踪误差,ξ1表示转换误差,
Figure FDA0003490846170000035
是连续时变函数,假设存在两个正常量C0和C1,则变量
Figure FDA0003490846170000036
满足以下条件:
Figure FDA0003490846170000037
已知
Figure FDA0003490846170000038
Figure FDA0003490846170000039
是有界的,由此推断,从公式(15)可以推断出ξ1的有界性取决于跟踪误差s1;对于每个初始条件
Figure FDA00034908461700000310
当ξ1在t∈[0,+∞)内有界时,推断出s1(t)满足如下条件:
Figure FDA00034908461700000311
分别将
Figure FDA00034908461700000312
Figure FDA00034908461700000313
缩写为
Figure FDA00034908461700000314
Figure FDA00034908461700000315
对ξ1求导,得到:
Figure FDA00034908461700000316
其中
Figure FDA00034908461700000317
Figure FDA00034908461700000318
基于(3)和(16),得到输出无约束子系统,如下式所示:
Figure FDA00034908461700000319
基于(15),得到无约束变量ξ1的值域是所有实数集,结合η11和假设1,推断出αic,i=2,3,此外,结合洛必达法则,得到:
Figure FDA00034908461700000320
当s1不受约束时,转换误差ξ1接近跟踪误差s1
B、基于结合辅助变量d,D,设计有限时间扰动观测器为:
Figure FDA00034908461700000321
其中
Figure FDA00034908461700000322
Figure FDA00034908461700000323
分别代表Δd,x2和D的估计值,定义
Figure FDA00034908461700000324
λ1和λ2并表示扰动观测器的正设计系数,κ1为适当的正常数;
定义
Figure FDA00034908461700000325
然后对
Figure FDA00034908461700000326
Figure FDA00034908461700000327
对时间t求导,得到:
Figure FDA0003490846170000041
假设不匹配动的变化率忽略不计,则Δd有界,并且估计值
Figure FDA0003490846170000042
Figure FDA0003490846170000043
将分别在有限的时间内收敛到实际值Δd,x2和D;
C、设计基于扰动观测器的自适应有限时间控制器:
对跟踪误差采用非线性坐标变换:
Figure FDA0003490846170000044
其中
Figure FDA0003490846170000045
是下列第i个一阶滤波器的输出:
Figure FDA0003490846170000046
其中εi代表第i个滤波器的时间常数,第i个虚拟控制信号
Figure FDA0003490846170000047
将在之后设计;
类似地,滤波器误差zi,i=2,3被定义为:
Figure FDA0003490846170000048
将(2)和(19)代入到(23)中,对(23)中的
Figure FDA0003490846170000049
求导:
Figure FDA00034908461700000410
定义第i个变量的估计误差
Figure FDA00034908461700000411
为:
Figure FDA00034908461700000412
其中第i个变量
Figure FDA00034908461700000413
的估计值为
Figure FDA00034908461700000414
基于扰动观测器的自适应有限时间控制器的设计步骤如下:
步骤1、选择的第一个子Lyapunov函数V1为:
Figure FDA00034908461700000415
其中设计常数τ1>0;
结合(27),对(28)中的V1求导:
Figure FDA00034908461700000416
将(26)代入(29),得到:
Figure FDA00034908461700000417
定义F1(X1)为:
Figure FDA00034908461700000418
其中
Figure FDA00034908461700000419
采用径向基神经网络估计未知函数F1(X1):
Figure FDA0003490846170000051
其中设计常数ΞM>0.
因此,(30)变为:
Figure FDA0003490846170000052
根据杨氏不等式,得到:
Figure FDA0003490846170000053
其中设计常数p1>0;
将(34)代入(33)得到:
Figure FDA0003490846170000054
设计虚构控制律
Figure FDA0003490846170000055
和自适应律
Figure FDA0003490846170000056
为:
Figure FDA0003490846170000057
其中k11,k12,
Figure FDA0003490846170000058
表示正常数,α=α12和奇数α1,a2满足如下条件:0<a1<a2
将(36)代入到(35),得到
Figure FDA0003490846170000059
结合(24)-(26),(27)和(36),对z2求导:
Figure FDA00034908461700000510
其中
Figure FDA00034908461700000511
为一个连续函数;
推导出
Figure FDA00034908461700000512
Figure FDA00034908461700000513
满足在紧集的预设初始条件下服从的最大值,因此:
Figure FDA00034908461700000514
其中
Figure FDA00034908461700000515
利用杨氏不等式,得到:
Figure FDA00034908461700000516
将(40)代入到(37)中,得到:
Figure FDA00034908461700000517
步骤2、选择第二个子Lyapunov函数V2为:
Figure FDA0003490846170000061
其中τ2>0为已知常量;
结合(27),对V2进行求导,得到:
Figure FDA0003490846170000062
将(26)和(41)代入到(43)中,得到:
Figure FDA0003490846170000063
基于引理5,假设
Figure FDA0003490846170000064
其中q1表示正常数;然后将(44)简化为:
Figure FDA0003490846170000065
构造F2(X2)
Figure FDA0003490846170000066
其中
Figure FDA0003490846170000067
结合(46),然后(45)重新表述为:
Figure FDA0003490846170000068
与(32)类似,得到(47)中的F2(X2)未知,因此,提出一个径向基神经网络来估计F2(X2):
Figure FDA0003490846170000069
因此,(47)被重新表达为:
Figure FDA00034908461700000610
类似于(34),得到如下不等式:
Figure FDA00034908461700000611
其中设计常数p2>0;
将(50)代入(49),得到:
Figure FDA0003490846170000071
类似于(36),选择虚拟控制器
Figure FDA0003490846170000072
和自适应律
Figure FDA0003490846170000073
为:
Figure FDA0003490846170000074
其中设计常数k21,k22,
Figure FDA0003490846170000075
为正数;
将(52)代入到(51)中,得到:
Figure FDA0003490846170000076
类似于(40),得到:
Figure FDA0003490846170000077
其中函数
Figure FDA0003490846170000078
将(54)代入(53),得到
Figure FDA0003490846170000079
步骤3、选择第三个子Lyapunov函数V3为:
Figure FDA00034908461700000710
其中已知常量t3>0;
根据(27),将V3对时间t求导:
Figure FDA00034908461700000711
将(26)和(55)代入(57),得到:
Figure FDA00034908461700000712
设计F3(X3)为:
Figure FDA00034908461700000713
其中
Figure FDA00034908461700000714
因此,(58)重新表示为:
Figure FDA0003490846170000081
同样,F3(X3)是未知的,然后,应用径向基神经网络来估计
Figure FDA0003490846170000082
Figure FDA0003490846170000083
类似于(34),得到:
Figure FDA0003490846170000084
其中设计参数p3>0;
因此,(60)被重新表达为:
Figure FDA0003490846170000085
设计实际控制器uq和自适应律
Figure FDA0003490846170000086
为:
Figure FDA0003490846170000087
其中设计常数k31,k32,
Figure FDA0003490846170000088
为正数;
将(64)代入(63),得到:
Figure FDA0003490846170000089
步骤4、选择第四个子Lyapunov函数V4为:
Figure FDA00034908461700000810
其中设计常数τ4>0;
结合(27),对(66)中的V4求导,得到:
Figure FDA00034908461700000811
结合(26)并将(65)代入(67)得到:
Figure FDA00034908461700000812
构造F4(X4)为:
Figure FDA00034908461700000813
其中
Figure FDA0003490846170000091
然后,(68)被简化为:
Figure FDA0003490846170000092
由于F4(X4)是一个不确定的函数,应用径向基神经网络来估计它,得到:
Figure FDA0003490846170000093
类似于(34),得到:
Figure FDA0003490846170000094
其中设计常数p4>0;
将(72)代入(70),得到:
Figure FDA0003490846170000095
设计实际控制器ud和虚拟控制律
Figure FDA0003490846170000096
Figure FDA0003490846170000097
其中k41,k42,
Figure FDA0003490846170000098
为正常数;
将(74)代入(73),得到:
Figure FDA0003490846170000099
CN202210096210.XA 2022-01-26 2022-01-26 具有扰动及输出约束的永磁同步电机系统有限时间控制法 Active CN114706300B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210096210.XA CN114706300B (zh) 2022-01-26 2022-01-26 具有扰动及输出约束的永磁同步电机系统有限时间控制法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210096210.XA CN114706300B (zh) 2022-01-26 2022-01-26 具有扰动及输出约束的永磁同步电机系统有限时间控制法

Publications (2)

Publication Number Publication Date
CN114706300A true CN114706300A (zh) 2022-07-05
CN114706300B CN114706300B (zh) 2024-05-14

Family

ID=82167141

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210096210.XA Active CN114706300B (zh) 2022-01-26 2022-01-26 具有扰动及输出约束的永磁同步电机系统有限时间控制法

Country Status (1)

Country Link
CN (1) CN114706300B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116447078A (zh) * 2023-04-21 2023-07-18 曲阜师范大学 一种闭环信息非对称约束在线调整机舱悬浮暂态性能的控制方法
CN116743019A (zh) * 2023-04-21 2023-09-12 曲阜师范大学 一种基于约束边界动态调整对称型障碍Lyapunov函数有限时间机舱悬浮控制方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5341078A (en) * 1989-09-30 1994-08-23 Fanuc Ltd. Sliding mode control method for a machine having an operating section operatively coupled to a servomotor and wherein a switching variable is determined based on a torsion amount and a torsion speed
CN106788046A (zh) * 2017-02-20 2017-05-31 青岛大学 永磁同步电机命令滤波有限时间模糊控制方法
CN108303873A (zh) * 2017-12-28 2018-07-20 浙江工业大学 一种考虑控制量受限的永磁同步电机滑模控制器
CN108803325A (zh) * 2018-06-06 2018-11-13 黄山学院 永磁同步电机伺服系统鲁棒有限时间控制方法
CN109873582A (zh) * 2019-01-29 2019-06-11 青岛大学 基于动态面的永磁同步电机有限时间位置跟踪控制方法
CN110347044A (zh) * 2019-07-15 2019-10-18 贵州大学 一种考虑输出约束的pmsm混沌系统神经网络动态面控制方法
US10725440B1 (en) * 2016-02-26 2020-07-28 The Mathworks, Inc. Systems and method for parameter estimation for a permanent magnet synchronous machine

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5341078A (en) * 1989-09-30 1994-08-23 Fanuc Ltd. Sliding mode control method for a machine having an operating section operatively coupled to a servomotor and wherein a switching variable is determined based on a torsion amount and a torsion speed
US10725440B1 (en) * 2016-02-26 2020-07-28 The Mathworks, Inc. Systems and method for parameter estimation for a permanent magnet synchronous machine
CN106788046A (zh) * 2017-02-20 2017-05-31 青岛大学 永磁同步电机命令滤波有限时间模糊控制方法
CN108303873A (zh) * 2017-12-28 2018-07-20 浙江工业大学 一种考虑控制量受限的永磁同步电机滑模控制器
CN108803325A (zh) * 2018-06-06 2018-11-13 黄山学院 永磁同步电机伺服系统鲁棒有限时间控制方法
CN109873582A (zh) * 2019-01-29 2019-06-11 青岛大学 基于动态面的永磁同步电机有限时间位置跟踪控制方法
CN110347044A (zh) * 2019-07-15 2019-10-18 贵州大学 一种考虑输出约束的pmsm混沌系统神经网络动态面控制方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
杨雪婷;马玉梅;于金鹏;赵林;韩瑶;: "基于有限时间永磁同步电机命令滤波位置跟踪控制", 青岛大学学报(工程技术版), no. 03, 31 August 2017 (2017-08-31) *
王孟孟;贺思艳;李振勇;于海生;于金鹏;: "永磁同步电机的神经网络离散位置跟踪控制", 青岛大学学报(工程技术版), no. 03, 29 August 2018 (2018-08-29) *
韩吉霞;马飞越;佃松宜;罗连杰;胡怡;: "基于非线性干扰观测器不确定系统的终端滑模控制", 电光与控制, no. 02, 29 February 2020 (2020-02-29) *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116447078A (zh) * 2023-04-21 2023-07-18 曲阜师范大学 一种闭环信息非对称约束在线调整机舱悬浮暂态性能的控制方法
CN116743019A (zh) * 2023-04-21 2023-09-12 曲阜师范大学 一种基于约束边界动态调整对称型障碍Lyapunov函数有限时间机舱悬浮控制方法
CN116447078B (zh) * 2023-04-21 2024-01-09 曲阜师范大学 一种闭环信息非对称约束在线调整机舱悬浮暂态性能的控制方法
CN116743019B (zh) * 2023-04-21 2024-01-26 曲阜师范大学 含约束边界调整的对称障碍Lyapunov函数机舱悬浮控制方法

Also Published As

Publication number Publication date
CN114706300B (zh) 2024-05-14

Similar Documents

Publication Publication Date Title
Gao et al. An adaptive super twisting nonlinear fractional order PID sliding mode control of permanent magnet synchronous motor speed regulation system based on extended state observer
CN114706300A (zh) 具有扰动及输出约束的永磁同步电机系统有限时间控制法
CN113206623B (zh) 基于快速积分终端滑模及干扰估计的永磁同步电机有限时间调速控制方法
CN112910362B (zh) 永磁同步电机无模型预测电流控制方法
CN114280944B (zh) 一种具有输出约束的pmsm系统有限时间动态面控制方法
Lei et al. Discrete-time command filtered adaptive fuzzy fault-tolerant control for induction motors with unknown load disturbances
Zhang et al. Low‐complexity adaptive tracking control of MIMO nonlinear systems with unknown control directions
Jiang et al. Full state constraints-based adaptive fuzzy finite-time command filtered control for permanent magnet synchronous motor stochastic systems
CN114421835A (zh) 一种基于模糊偏差耦合反步滑模策略的多电机控制方法
Duff et al. Approximation of stability regions for large-scale time-delay systems using model reduction techniques
Zheng et al. A novel reduced-order load torque observer based discrete-time sliding mode control for PMSM speed servo system
Vo et al. Speed estimators using stator resistance adaptation for sensorless induction motor drive
Li et al. Disturbance‐observer‐based adaptive finite‐time dynamic surface control for PMSM with time‐varying asymmetric output constraint
Regaya et al. Adaptive proportional-integral fuzzy logic controller of electric motor drive
CN115313939A (zh) 永磁同步电机随机系统指令滤波模糊自适应控制方法
CN115102443A (zh) 一种永磁同步直线电机的控制方法、装置及存储介质
Gao et al. The neural network control approach for PMSM based on a high gain observer
Qu et al. Sliding-mode anti-disturbance speed control of permanent magnet synchronous motor based on an advanced reaching law
Liang et al. Dual-redundancy PMSM servo system: using single neuron PID controller
CN114499306B (zh) 具有时滞和扰动的永磁同步电机的自适应funnel动态面控制方法
Aguilar Mejía et al. Adaptive Speed Controller for a Permanent Magnet Synchronous Motor
CN111208728A (zh) 一种线性自抗扰控制方法、装置设备及存储介质
Zhao et al. Design of MRAC and Modified MRAC for the Turntable
Lee et al. Output-feedback linear parameter varying control of permanent magnet synchronous motors
CN113659894B (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