CN113644857A - 一种永磁直线电机有限级量化迭代学习控制方法 - Google Patents

一种永磁直线电机有限级量化迭代学习控制方法 Download PDF

Info

Publication number
CN113644857A
CN113644857A CN202110900720.3A CN202110900720A CN113644857A CN 113644857 A CN113644857 A CN 113644857A CN 202110900720 A CN202110900720 A CN 202110900720A CN 113644857 A CN113644857 A CN 113644857A
Authority
CN
China
Prior art keywords
quantization
input
equation
permanent magnet
output
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
CN202110900720.3A
Other languages
English (en)
Other versions
CN113644857B (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.)
Jiangnan University
Original Assignee
Jiangnan 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 Jiangnan University filed Critical Jiangnan University
Priority to CN202110900720.3A priority Critical patent/CN113644857B/zh
Publication of CN113644857A publication Critical patent/CN113644857A/zh
Application granted granted Critical
Publication of CN113644857B publication Critical patent/CN113644857B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P25/00Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
    • H02P25/06Linear motors
    • 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
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (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)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

本发明公开了一种永磁直线电机有限级量化迭代学习控制方法,涉及电机优化控制领域,该方法基于提升技术将重复运行的永磁直线电机转换为时间序列的输入输出矩阵模型,并将有限级对数量化器与编码解码机制相结合用以量化信号,从而减少传输数据量。针对采用编码解码方案量化信号的***,基于优化的思想设计迭代学习控制算法,根据性能指标函数得到迭代学习优化控制算法的前馈实现。基于压缩映射方法,证明了所设计的量化迭代学习优化算法在数学期望意义下的收敛性。该方法可以解决使用网络传输信号的永磁直线电机的跟踪控制问题,实现对期望轨迹的高精度跟踪。

Description

一种永磁直线电机有限级量化迭代学习控制方法
技术领域
本发明涉及电机优化控制领域,尤其是一种永磁直线电机有限级量化迭代学习控制方法。
背景技术
永磁直线电机是一种无中间传动机构的新型电机,可将电能直接转换为直线运动的机械能。由于其具有结构简单,推力体积比大,效率高,定位精确等优点,因而广泛应用于工业运输***、办公自动化及军事领域等。
对于执行重复性任务的永磁直线电机,迭代学习控制是一种能够有效提高轨迹跟踪精度的控制策略。而随着计算机与通信技术的发展,在永磁直线电机控制中采用网络控制***方案成为可能,其中永磁直线电机与迭代学习控制器位于不同站点并通过无线网络彼此通信,但受限于通信带宽有限及出于可靠性的考虑,需要减少数据传输量,量化信号作为一种解决该问题的常用技术手段,在实际应用中具有很大的优势,既可以实现***跟踪性能的要求,又能降低传输负担。因此对于采用网络通信的永磁直线电机,出现了量化数据背景下的迭代学习控制器设计问题。
在标称***上,许多迭代学习控制算法都取得了精确的跟踪性能和理想的收敛性能,并在实际应用中取得了良好的性能。然而,当传输信号存在量化误差时,常规的迭代学习控制算法的控制效果并不理想,由于实际应用中只能使用有限级量化器,而不是理想条件下的无限级量化器,这对控制器的设计提出了更高的要求。
发明内容
本发明人针对上述问题及技术需求,提出了一种永磁直线电机有限级量化迭代学习控制方法,将对数量化器与编码解码机制相结合设计了量化编码解码器,以迭代的方式实现信号精确传输,采用范数优化算法框架设计量化迭代学习控制算法,根据压缩映射法得到***的鲁棒有界收敛条件。
本发明的技术方案如下:
一种永磁直线电机有限级量化迭代学习控制方法,包括如下步骤:
第一步、建立永磁直线电机的动态模型:
动态模型采用下述动力学方程表示:
Figure BDA0003199692170000021
其中,R表示定子电阻,m表示电机动子部分的总质量,ψf表示定子永磁体励磁磁链,k1=π/τ,k2=1.5π/τ,τ表示极距,u(t)表示动子电压,p(t)表示电机位置,ω(t)表示电机动子速度;
第二步、构建永磁直线电机的离散状态空间方程:
选取满足香农采样定理采样周期Ts,使用欧拉法对连续***模型式(1)进行离散化,得到永磁直线电机的离散***模型如下:
Figure BDA0003199692170000022
将永磁直线电机的位置、动子速度定义为状态变量:x=[p ω]T,定义输入变量为动子电压u,输出为电机转速ω,则式(1)所示的永磁直线电机的动态模型描述为:
Figure BDA0003199692170000023
式中t和k分别代表采样时间和批次,批次过程的运行周期为T;在每个重复过程周期内,对于时间点t∈[0,T]取N个采样点;uk(t),yk(t)和xk(t)分别是该***第k批次t时刻的输入、输出和状态向量;A,B,C为式(3)中离散***各个参数矩阵,且满足CB≠0;并且假设***每个批次的初始状态一致,即xk(0)=x0
第三步、建立轨迹跟踪模型:
针对式(3)形式的线性离散***,将其状态空间表达式转换为时间序列的输入输出矩阵模型:
yk=Guk+dk (4)
其中:
Figure BDA0003199692170000031
dk=[CA CA2 CA3 … CAN]Txk(0)
uk=[uk(0),uk(1),...,uk(N-1)]T
yk=[yk(1),yk(2),...,yk(N)]T
G是时间序列上的输入输出传递矩阵;dk是***初始状态对输出的影响,假设xk(0)=0,
Figure BDA0003199692170000032
则dk=0;
第四步、设计量化编码解码器:
在采用网络传输信号的***中,由于网络带宽有限,需要在减少传输的信息量条件下实现精确跟踪轨迹的目标,因此可对信号进行量化处理,设计输入端量化编码解码器如下:
Figure BDA0003199692170000033
Figure BDA0003199692170000034
其中,0表示与***输入具有相同维度的零向量,uk(t)、
Figure BDA0003199692170000035
Figure BDA0003199692170000036
分别为编码器E1的输入、输出和内部状态;
Figure BDA0003199692170000037
是解码器D1的输出,即控制器输出uk(t)的估计值;q(·)是由式(7)定义的有限级对数量化器:
Figure BDA0003199692170000038
其中,v表示有限级对数量化器的输入;
Figure BDA0003199692170000039
μ为所选择的量化密度,量化水平如下:
Z={±zi|zi=μiz0,i=0,1,2,…,L-1}∪{0},0<μ<1,z0>0
其中,L表示正量化级数,z0>0表示最大允许级别,实际应用中可选择足够大的z0使得
Figure BDA0003199692170000041
避免量化器饱和以简化控制器设计;
设计输出端量化编码解码器如下:
Figure BDA0003199692170000042
Figure BDA0003199692170000043
其中,0表示与***输出具有相同维度的零向量,yk(t)、
Figure BDA0003199692170000044
Figure BDA0003199692170000045
分别为编码器E2的输入、输出和内部状态;
Figure BDA0003199692170000046
是解码器D2的输出,即***输出yk(t)的估计值;q(·)是由式(7)定义的有限级对数量化器,输入端有限级对数量化器的参数设置与输出端有限级对数量化器的参数设置可以不同;
第五步、建立量化前后信号的关系表达式:
有限级对数量化器输入v与输出q(v)存在量化误差Δv,为了便于控制器的设计,需要给出输入v与输出q(v)之间的关系;当
Figure BDA0003199692170000047
时,根据扇形有界得到q(v)=v+η·v=(1+η)v,|η|≤δ;当
Figure BDA0003199692170000048
时可得
Figure BDA0003199692170000049
因此得到有限级对数量化器输入v与输出q(v)的关系式如下:
q(v)=(1+η)v+d (10)
其中,η表示相对量化误差,满足|η|≤δ;d表示传输误差,满足
Figure BDA00031996921700000410
对于k批次t时刻的信号有:
q(vk(t))=(1+ηk(t))vk(t)+dk(t) (11)
其中,ηk(t)表示k批次t时刻信号的相对量化误差,dk(t)表示k批次t时刻信号的传输误差;
根据输入端量化编码解码器的定义与式(11),得到:
Figure BDA00031996921700000411
根据数学归纳法有
Figure BDA00031996921700000412
成立,于是
Figure BDA00031996921700000413
与uk+1(t)的关系式为:
Figure BDA0003199692170000051
将式(13)提升为向量形式,得到
Figure BDA0003199692170000052
与uk+1的向量关系式为:
Figure BDA0003199692170000053
其中:
Figure BDA0003199692170000054
根据输出端量化编码解码器的定义得到
Figure BDA0003199692170000055
与yk+1的向量关系式为:
Figure BDA0003199692170000056
其中:
Figure BDA0003199692170000057
相对量化误差独立于量化器输入,因此对于任意的k批次t时刻有E[vk(t)ηk(t)]=0;不同量化器产生的相对量化误差也是相互独立的,即
Figure BDA0003199692170000058
i≠j且i,j∈{1,2};而在相同的量化器中,相对量化误差ηk(t)在区间[-δ,δ]内均匀分布且对于任意k1,k2和t1,t2满足:
Figure BDA0003199692170000059
其中,
Figure BDA00031996921700000510
因此,通过取数学期望的方式化简含有
Figure BDA00031996921700000511
Figure BDA00031996921700000512
的表达式;
***中存在实际跟踪误差ek+1=yd-yk+1与辅助校正误差
Figure BDA00031996921700000513
真正体现***跟踪性能的是实际跟踪误差,控制器使用辅助校正误差修正当前批次输入信号;根据
Figure BDA00031996921700000514
得到:
Figure BDA00031996921700000515
因此,建立当前实际跟踪误差序列与前一批次辅助校正误差序列的关系式:
Figure BDA00031996921700000516
其中,定义
Figure BDA00031996921700000517
为传输误差,
Figure BDA00031996921700000518
β的取值与两侧量化器的参数选择有关,若输入输出端量化器选择的参数相同,则简化为
Figure BDA0003199692170000061
第六步、设计量化迭代学习控制轨迹跟踪优化算法:
考虑范数最优迭代学习控制框架,每批次的控制输入通过优化一个性能指标函数得到,性能指标函数的一般形式为:
Figure BDA0003199692170000062
性能指标函数包括实际跟踪误差和控制振荡,在优化过程中,分别用对称正定权重矩阵Q和R来表示其优先级,即Q=QT>0,R=RT>0;
取权重矩阵Q=qI,R=rI,并且诱导范数被定义为如下形式:
Figure BDA0003199692170000063
Figure BDA0003199692170000064
中不包含随机变量,因而其期望等于其本身;
定义
Figure BDA0003199692170000065
Ξ=σ2diag(γ1122,…,γNN),根据式(18)得到:
Figure BDA0003199692170000066
其中Ξ与σ2I都是对称正定矩阵;将式(20)代入性能指标函数式(19)得到:
Figure BDA0003199692170000067
有限级量化迭代学习控制更新律通过求解一个具有混合参数的优化问题得到,具体形式为:
Figure BDA0003199692170000068
为了解决式(22)的最小-最大化问题,利用拉格朗日对偶函数将最小-最大化问题转化为最小-最小化问题;因此,将最小-最大化问题重新表述为一个凸优化问题,从而保证全局最优解;
首先考虑关于传输误差wk+1的最大化问题,因为性能指标函数(21)只有第一项与wk+1有关,因此优化式(22)内部的最大化问题表示为:
Figure BDA0003199692170000069
最大化问题具有很强的对偶性;由于式(23)是带有约束的最大化问题,引入拉格朗日乘子λk+1,得到如下拉格朗日函数:
Figure BDA0003199692170000071
其中,拉格朗日乘子λk+1≥0,根据KKT条件,通过对wk+1和λk+1进行微分从而得到最优解,并用
Figure BDA0003199692170000072
Figure BDA0003199692170000073
表示最优解得到:
Figure BDA0003199692170000074
为了保证原性能指标函数的凸优化要求,拉格朗日函数(24)的海森矩阵应是负半定的,即
Figure BDA0003199692170000075
因此
Figure BDA0003199692170000076
可逆,用
Figure BDA0003199692170000077
表示其伪逆,则有:
Figure BDA0003199692170000078
Figure BDA0003199692170000079
是下列方程的解:
Figure BDA00031996921700000710
由式(26)和式(27)可以看出,最优解
Figure BDA00031996921700000711
Figure BDA00031996921700000712
的取值均依赖于uk+1,得到L(wk+1k+1)的对偶函数为:
Figure BDA00031996921700000713
将最大化问题(23)转化成对偶函数g(λk+1)的最小化问题,即:
Figure BDA00031996921700000714
其中g(λk+1)的具体表达式为:
Figure BDA00031996921700000715
将原始优化问题(22)外部的最小化问题与最小化问题(29)相结合,得到新的对偶函数,即原始优化问题(22)的最小化-最大化问题转化如下最小化-最小化问题:
Figure BDA00031996921700000716
其中Jdual(uk+1k+1)表示新的对偶性能指标函数且为凸函数,且有:
Figure BDA00031996921700000717
由于式(32)是凸函数,因此存在有效的数值算法来计算式(31)的全局最优解;其中最小化问题
Figure BDA0003199692170000081
是关于λk+1的凸优化问题,因此最优解
Figure BDA0003199692170000082
由对偶函数通过对λk+1微分得到,得到:
Figure BDA0003199692170000083
其中,
Figure BDA0003199692170000084
Figure BDA0003199692170000085
时,Qk+1是正定矩阵;
对uk+1进行微分得到:
Figure BDA0003199692170000086
由于(GTQk+1G+Ξ+R)可逆,将式(34)整理后得到有限级量化迭代学习控制更新律为:
Figure BDA0003199692170000087
其中:
Figure BDA0003199692170000088
增益矩阵Qk+1随批次改变,其值取决于优化问题(31)的解
Figure BDA0003199692170000089
第七步、分析量化迭代学习轨迹跟踪优化算法的收敛性:
对于期望跟踪轨迹yd,存在一个理想的输入ud满足yd=Gud,定义
Figure BDA00031996921700000810
则对于k+1批次有:
Figure BDA00031996921700000811
将控制更新律(35)代入式(37),得到:
Figure BDA00031996921700000812
根据输入端信号的关系式,得到辅助校正误差的等价形式为:
Figure BDA00031996921700000813
由于
Figure BDA00031996921700000814
是可逆矩阵,根据
Figure BDA00031996921700000815
有:
Figure BDA00031996921700000816
将式(39)、式(40)代入式(38)整理得到:
Figure BDA0003199692170000091
对式(41)两边取期望得到:
Figure BDA0003199692170000092
根据式(36)得到I-Ku-Kζ=0,令
Figure BDA0003199692170000093
对式(42)两边取范数,得到E(Δuk+1)的不等式为:
||E(Δuk+1)||≤||Ku-KeG+Kζ||||E(Δuk)||+b (43)
***经过k次迭代后,E(Δuk+1)的不等式转化为:
Figure BDA0003199692170000094
若选择的权重矩阵与量化密度使得约束条件成立:
||Ku-KeG+Kζ||≤ρ<1 (45)
根据压缩映射引理得到
Figure BDA0003199692170000095
因此式(43)简化为:
Figure BDA0003199692170000096
记||G||=c,根据E(ek)=GE(Δuk)得到:
Figure BDA0003199692170000097
即期望意义下的误差范数||E(ek)||收敛至一个有界值;
第八步、利用量化编码解码器的量化信号实现永磁直线电机的轨迹跟踪:
根据有限级量化迭代学习控制律确定永磁直线电机每一迭代批次的生成输入矢量,经量化编码解码器作用得到实际输入矢量,以实际输入矢量对永磁直线电机进行轨迹跟踪控制,永磁直线电机在实际输入矢量的控制作用下跟踪期望输出。
本发明的有益技术效果是:
本申请公开了针对永磁直线电机此类具有重复运动特征的线性***,将永磁直线电机作为被控对象,针对网络传输信号带宽有限的情况,设计量化编码解码器,将对数量化器与信号编码机制相结合,从而在迭代的过程中逐步减轻量化误差对跟踪性能的影响,以提高传输信号精度。本方法基于范数优化迭代学习框架,设计了一种量化迭代学习控制算法,保证了***跟踪误差在数学期望意义下的收敛性。
附图说明
图1是本申请提供的永磁直线电机的模型框图。
图2是本申请提供的永磁直线电机的实际输出曲线图。
图3是本申请提供的***实际跟踪误差2-范数收敛图。
图4是本申请提供的控制器计算得到的生成输入电压曲线图。
图5是本申请提供的永磁直线电机实际输入电压曲线图。
图6是本申请提供的权重矩阵Qk+1变化曲线图。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步说明。
结合图1-图6所示,请参考图1,其示出了本申请公开的永磁直线电机的模型框图。第k批次的控制器输出为uk,经编码器E1编码后通过网络传输,由解码器D1接收并解码得到实际控制矢量
Figure BDA0003199692170000101
作用于永磁直线电机可以得到***第k批次的实际输出yk,其与储存在期望轨迹存储器的设定期望值yd进行比较得到实际跟踪误差ek。将实际跟踪误差精度与设定的精度值做比较,若误差精度没有达到所设定的精度,则将实际输出yk经编码器E2编码后通过网络传输,由解码器D2接收并解码得到估计的输出值
Figure BDA0003199692170000102
其与储存在期望轨迹存储器的设定期望值进行比较得到辅助校正误差
Figure BDA0003199692170000103
将辅助校正误差
Figure BDA0003199692170000104
当前控制器输入uk与编码器E1的内部状态量
Figure BDA0003199692170000105
传递到优化迭代学习控制器生成下一批次的控制器输出uk+1,如此循环运行直至***实际输出与期望值之间的误差达到精度要求,则停止迭代,此时的控制器输入即为最优控制输入。
针对式(1)所示的永磁直线电机实际物理模型(也即动态模型),变量参数分别设定为:
R=8.6Ω,m=1.635kg,ψf=0.35Wb,τ=0.031m
***仿真时间设定为t=2s,采样时间设定为Ts=0.01s,则***的离散状态空间表达式的参数矩阵分别为:
Figure BDA0003199692170000106
本实施方式设定永磁直线电机的期望轨迹为:
Figure BDA0003199692170000111
单位为rad/s,同时令初始状态满足xk(0)=0。采用网络传输信号,由于带宽有限设计了有限级量化编码解码器,其中对数量化器的参数设置为μ=0.7。
选取权重矩阵Q=100I,R=0.1I,满足式(45)。当权重矩阵Q、R与量化密度μ确定时,量化迭代学习控制律中的Ku,Ke,Kζ也随之确定。本申请的上述量化迭代学习控制器基于STM32F103RCT6芯片实现,芯片的输入为电机控制电压u,并通过电压传感器采集得到。输入信号通过调理电路进入STM32F103RCT6芯片进行存储和计算,并构建迭代学习更新律,CPU计算后得到的信号为生成输入信号uk+1,经编码器E1编码后通过网络传输,由解码器D1接收并解码得到实际控制信号
Figure BDA0003199692170000112
实际控制信号作用于永磁直线电机,不断修正输出轨迹,直到跟踪上期望轨迹。永磁直线电机的动态模型(1)运行时,请参考图2,其示出了永磁直线电机应用量化迭代学习控制律(35)的轨迹跟踪效果图,经过一定批次k后,***的输出能准确跟踪到期望轨迹。图3表明量化迭代学习控制算法经过一定的迭代批次后依旧能够实现有界收敛,这也验证了算法的合理性以及有效性。
图4和图5分别为控制器计算得到的生成输入信号与***实际输入信号曲线图。对比图4与图5可知,随着迭代次数的增加,实际输入信号从初始的阶梯状逐渐变得平滑,经过多次迭代后实际输入信号与生成输入信号曲线基本一致,说明本申请所设计的量化编码解码器能够以迭代的方式实现信号的精确传输。图6表明随着迭代的进行权重矩阵Qk+1逐渐增大。
以上所述的仅是本申请的优选实施方式,本发明不限于以上实施例。可以理解,本领域技术人员在不脱离本发明的精神和构思的前提下直接导出或联想到的其他改进和变化,均应认为包含在本发明的保护范围之内。

Claims (1)

1.一种永磁直线电机有限级量化迭代学习控制方法,其特征在于,包括如下步骤:
第一步、建立永磁直线电机的动态模型:
所述动态模型采用下述动力学方程表示:
Figure FDA0003199692160000011
其中,R表示定子电阻,m表示电机动子部分的总质量,ψf表示定子永磁体励磁磁链,k1=π/τ,k2=1.5π/τ,τ表示极距,u(t)表示动子电压,p(t)表示电机位置,ω(t)表示电机动子速度;
第二步、构建永磁直线电机的离散状态空间方程:
选取满足香农采样定理采样周期Ts,使用欧拉法对连续***模型式(1)进行离散化,得到所述永磁直线电机的离散***模型如下:
Figure FDA0003199692160000012
将所述永磁直线电机的位置、动子速度定义为状态变量:x=[p ω]T,定义输入变量为动子电压u,输出为电机转速ω,则式(1)所示的永磁直线电机的动态模型描述为:
Figure FDA0003199692160000013
式中t和k分别代表采样时间和批次,批次过程的运行周期为T;在每个重复过程周期内,对于时间点t∈[0,T]取N个采样点;uk(t),yk(t)和xk(t)分别是该***第k批次t时刻的输入、输出和状态向量;A,B,C为式(3)中离散***各个参数矩阵,且满足CB≠0;并且假设***每个批次的初始状态一致,即xk(0)=x0
第三步、建立轨迹跟踪模型:
针对式(3)形式的线性离散***,将其状态空间表达式转换为时间序列的输入输出矩阵模型:
yk=Guk+dk (4)
其中:
Figure FDA0003199692160000021
dk=[CA CA2 CA3…CAN]Txk(0)
uk=[uk(0),uk(1),...,uk(N-1)]T
yk=[yk(1),yk(2),...,yk(N)]T
G是时间序列上的输入输出传递矩阵;dk是***初始状态对输出的影响,假设
Figure FDA0003199692160000022
则dk=0;
第四步、设计量化编码解码器:
在采用网络传输信号的***中,对信号进行量化处理,设计输入端量化编码解码器如下:
Figure FDA0003199692160000023
Figure FDA0003199692160000024
其中,0表示与***输入具有相同维度的零向量,uk(t)、
Figure FDA0003199692160000025
Figure FDA0003199692160000026
分别为编码器E1的输入、输出和内部状态;
Figure FDA0003199692160000027
是解码器D1的输出,即控制器输出uk(t)的估计值;q(·)是由式(7)定义的有限级对数量化器:
Figure FDA0003199692160000028
其中,v表示所述有限级对数量化器的输入;
Figure FDA0003199692160000029
μ为所选择的量化密度,量化水平如下:
Z={±zi|zi=μiz0,i=0,1,2,…,L-1}∪{0},0<μ<1,z0>0
其中,L表示正量化级数,z0>0表示最大允许级别;
设计输出端量化编码解码器如下:
Figure FDA0003199692160000031
Figure FDA0003199692160000032
其中,0表示与***输出具有相同维度的零向量,yk(t)、
Figure FDA0003199692160000033
Figure FDA0003199692160000034
分别为编码器E2的输入、输出和内部状态;
Figure FDA0003199692160000035
是解码器D2的输出,即***输出yk(t)的估计值;q(·)是由式(7)定义的有限级对数量化器;输入端有限级对数量化器的参数设置与输出端有限级对数量化器的参数设置可以不同;
第五步、建立量化前后信号的关系表达式:
有限级对数量化器输入v与输出q(v)存在量化误差Δv,当
Figure FDA0003199692160000036
时,根据扇形有界得到q(v)=v+η·v=(1+η)v,|η|≤δ;当
Figure FDA0003199692160000037
时可得q(v)=0=v+d,
Figure FDA0003199692160000038
因此得到所述有限级对数量化器输入v与输出q(v)的关系式如下:
q(v)=(1+η)v+d (10)
其中,η表示相对量化误差,满足|η|≤δ;d表示传输误差,满足
Figure FDA0003199692160000039
对于k批次t时刻的信号有:
q(vk(t))=(1+ηk(t))vk(t)+dk(t) (11)
其中,ηk(t)表示k批次t时刻信号的相对量化误差,dk(t)表示k批次t时刻信号的传输误差;
根据所述输入端量化编码解码器的定义与式(11),得到:
Figure FDA00031996921600000310
根据数学归纳法有
Figure FDA00031996921600000311
成立,于是
Figure FDA00031996921600000312
与uk+1(t)的关系式为:
Figure FDA00031996921600000313
将式(13)提升为向量形式,得到
Figure FDA00031996921600000314
与uk+1的向量关系式为:
Figure FDA0003199692160000041
其中:
Figure FDA0003199692160000042
根据所述输出端量化编码解码器的定义得到
Figure FDA0003199692160000043
与yk+1的向量关系式为:
Figure FDA0003199692160000044
其中:
Figure FDA0003199692160000045
所述相对量化误差独立于量化器输入,因此对于任意的k批次t时刻有E[vk(t)ηk(t)]=0;不同量化器产生的相对量化误差也是相互独立的,即
Figure FDA0003199692160000046
i≠j且i,j∈{1,2};而在相同的量化器中,相对量化误差ηk(t)在区间[-δ,δ]内均匀分布且对于任意k1,k2和t1,t2满足:
Figure FDA0003199692160000047
其中,
Figure FDA0003199692160000048
因此,通过取数学期望的方式化简含有
Figure FDA0003199692160000049
Figure FDA00031996921600000410
的表达式;
***中存在实际跟踪误差ek+1=yd-yk+1与辅助校正误差
Figure FDA00031996921600000411
真正体现***跟踪性能的是所述实际跟踪误差,控制器使用所述辅助校正误差修正当前批次输入信号;根据
Figure FDA00031996921600000412
得到:
Figure FDA00031996921600000413
因此,建立当前实际跟踪误差序列与前一批次辅助校正误差序列的关系式:
Figure FDA00031996921600000414
其中,定义
Figure FDA00031996921600000415
为传输误差,
Figure FDA00031996921600000416
Figure FDA00031996921600000417
β的取值与两侧量化器的参数选择有关,若输入输出端量化器选择的参数相同,则简化为
Figure FDA0003199692160000051
第六步、设计量化迭代学习控制轨迹跟踪优化算法:
考虑范数最优迭代学习控制框架,每批次的控制输入通过优化一个性能指标函数得到,所述性能指标函数的一般形式为:
Figure FDA0003199692160000052
所述性能指标函数包括所述实际跟踪误差和控制振荡,在优化过程中,分别用对称正定权重矩阵Q和R来表示其优先级,即Q=QT>0,R=RT>0;
取权重矩阵Q=qI,R=rI,并且诱导范数被定义为如下形式:
Figure FDA0003199692160000053
Figure FDA0003199692160000054
中不包含随机变量,因而其期望等于其本身;
定义
Figure FDA0003199692160000055
Ξ=σ2diag(γ1122,…,γNN),根据式(18)得到:
Figure FDA0003199692160000056
其中Ξ与σ2I都是对称正定矩阵;将式(20)代入性能指标函数式(19)得到:
Figure FDA0003199692160000057
有限级量化迭代学习控制更新律通过求解一个具有混合参数的优化问题得到,具体形式为:
Figure FDA0003199692160000058
为了解决式(22)的最小-最大化问题,利用拉格朗日对偶函数将所述最小-最大化问题转化为最小-最小化问题;因此,将所述最小-最大化问题重新表述为一个凸优化问题,从而保证全局最优解;
首先考虑关于传输误差wk+1的最大化问题,因为性能指标函数(21)只有第一项与wk+1有关,因此优化式(22)内部的最大化问题表示为:
Figure FDA0003199692160000059
所述最大化问题具有很强的对偶性;由于式(23)是带有约束的最大化问题,引入拉格朗日乘子λk+1,得到如下拉格朗日函数:
Figure FDA0003199692160000061
其中,拉格朗日乘子λk+1≥0,根据KKT条件,通过对wk+1和λk+1进行微分从而得到最优解,并用
Figure FDA0003199692160000062
Figure FDA0003199692160000063
表示最优解得到:
Figure FDA0003199692160000064
为了保证原性能指标函数的凸优化要求,拉格朗日函数(24)的海森矩阵应是负半定的,即
Figure FDA0003199692160000065
因此
Figure FDA0003199692160000066
可逆,用
Figure FDA0003199692160000067
表示其伪逆,则有:
Figure FDA0003199692160000068
Figure FDA0003199692160000069
是下列方程的解:
Figure FDA00031996921600000610
由式(26)和式(27)可以看出,最优解
Figure FDA00031996921600000611
Figure FDA00031996921600000612
的取值均依赖于uk+1,得到L(wk+1k+1)的对偶函数为:
Figure FDA00031996921600000613
将最大化问题(23)转化成对偶函数g(λk+1)的最小化问题,即:
Figure FDA00031996921600000614
其中g(λk+1)的具体表达式为:
Figure FDA00031996921600000615
将原始优化问题(22)外部的最小化问题与最小化问题(29)相结合,得到新的对偶函数,即原始优化问题(22)的最小化-最大化问题转化如下最小化-最小化问题:
Figure FDA00031996921600000616
其中Jdual(uk+1k+1)表示新的对偶性能指标函数且为凸函数,且有:
Figure FDA00031996921600000617
由于式(32)是凸函数,因此存在有效的数值算法来计算式(31)的全局最优解;其中最小化问题
Figure FDA0003199692160000071
是关于λk+1的凸优化问题,因此最优解
Figure FDA0003199692160000072
由对偶函数通过对λk+1微分得到,得到:
Figure FDA0003199692160000073
其中,
Figure FDA0003199692160000074
Figure FDA0003199692160000075
时,Qk+1是正定矩阵;
对uk+1进行微分得到:
Figure FDA0003199692160000076
由于(GTQk+1G+Ξ+R)可逆,将式(34)整理后得到所述有限级量化迭代学习控制更新律为:
Figure FDA0003199692160000077
其中:
Figure FDA0003199692160000078
增益矩阵Qk+1随批次改变,其值取决于优化问题(31)的解
Figure FDA0003199692160000079
第七步、分析量化迭代学习轨迹跟踪优化算法的收敛性:
对于期望跟踪轨迹yd,存在一个理想的输入ud满足yd=Gud,定义
Figure FDA00031996921600000710
则对于k+1批次有:
Figure FDA00031996921600000711
将控制更新律(35)代入式(37),得到:
Figure FDA00031996921600000712
根据输入端信号的关系式,得到所述辅助校正误差的等价形式为:
Figure FDA00031996921600000713
由于
Figure FDA00031996921600000714
是可逆矩阵,根据
Figure FDA00031996921600000715
有:
Figure FDA00031996921600000716
将式(39)、式(40)代入式(38)整理得到:
Figure FDA0003199692160000081
对式(41)两边取期望得到:
Figure FDA0003199692160000082
根据式(36)得到I-Ku-Kζ=0,令
Figure FDA0003199692160000083
对式(42)两边取范数,得到E(Δuk+1)的不等式为:
||E(Δuk+1)||≤||Ku-KeG+Kζ||||E(Δuk)||+b (43)
***经过k次迭代后,E(Δuk+1)的不等式转化为:
Figure FDA0003199692160000084
若选择的权重矩阵与量化密度使得约束条件成立:
||Ku-KeG+Kζ||≤ρ<1 (45)
根据压缩映射引理得到
Figure FDA0003199692160000085
因此式(43)简化为:
Figure FDA0003199692160000086
记||G||=c,根据E(ek)=GE(Δuk)得到:
Figure FDA0003199692160000087
即期望意义下的误差范数||E(ek)||收敛至一个有界值;
第八步、利用所述量化编码解码器的量化信号实现所述永磁直线电机的轨迹跟踪:
根据所述有限级量化迭代学习控制律确定永磁直线电机每一迭代批次的生成输入矢量,经所述量化编码解码器作用得到实际输入矢量,以所述实际输入矢量对永磁直线电机进行轨迹跟踪控制,所述永磁直线电机在所述实际输入矢量的控制作用下跟踪期望输出。
CN202110900720.3A 2021-08-06 2021-08-06 一种永磁直线电机有限级量化迭代学习控制方法 Active CN113644857B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110900720.3A CN113644857B (zh) 2021-08-06 2021-08-06 一种永磁直线电机有限级量化迭代学习控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110900720.3A CN113644857B (zh) 2021-08-06 2021-08-06 一种永磁直线电机有限级量化迭代学习控制方法

Publications (2)

Publication Number Publication Date
CN113644857A true CN113644857A (zh) 2021-11-12
CN113644857B CN113644857B (zh) 2023-08-18

Family

ID=78419844

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110900720.3A Active CN113644857B (zh) 2021-08-06 2021-08-06 一种永磁直线电机有限级量化迭代学习控制方法

Country Status (1)

Country Link
CN (1) CN113644857B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114347020A (zh) * 2021-12-22 2022-04-15 深圳市优必选科技股份有限公司 一种运动控制方法、运动控制装置及机器人

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102710212A (zh) * 2012-06-11 2012-10-03 桂林电子科技大学 永磁同步直线电机改进的迭代学习控制方法与控制***
CN106953575A (zh) * 2017-03-14 2017-07-14 中南大学 一种迭代学习的直接转矩控制方法与装置
CN107070341A (zh) * 2017-03-24 2017-08-18 中国科学院长春光学精密机械与物理研究所 基于鲁棒迭代学习控制的永磁同步电机转矩脉动抑制方法
CN107272409A (zh) * 2017-07-04 2017-10-20 浙江理工大学 一种基于迭代学习的直线伺服***振动抑制方法
CN109709802A (zh) * 2018-12-14 2019-05-03 江南大学 基于迭代学习控制的有源电子梯形电路的控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102710212A (zh) * 2012-06-11 2012-10-03 桂林电子科技大学 永磁同步直线电机改进的迭代学习控制方法与控制***
CN106953575A (zh) * 2017-03-14 2017-07-14 中南大学 一种迭代学习的直接转矩控制方法与装置
CN107070341A (zh) * 2017-03-24 2017-08-18 中国科学院长春光学精密机械与物理研究所 基于鲁棒迭代学习控制的永磁同步电机转矩脉动抑制方法
CN107272409A (zh) * 2017-07-04 2017-10-20 浙江理工大学 一种基于迭代学习的直线伺服***振动抑制方法
CN109709802A (zh) * 2018-12-14 2019-05-03 江南大学 基于迭代学习控制的有源电子梯形电路的控制方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114347020A (zh) * 2021-12-22 2022-04-15 深圳市优必选科技股份有限公司 一种运动控制方法、运动控制装置及机器人
CN114347020B (zh) * 2021-12-22 2024-02-13 深圳市优必选科技股份有限公司 一种运动控制方法、运动控制装置及机器人

Also Published As

Publication number Publication date
CN113644857B (zh) 2023-08-18

Similar Documents

Publication Publication Date Title
CN112085193B (zh) 一种神经网络的量化参数确定方法及相关产品
Kostina et al. Rate-cost tradeoffs in control
CN113625563B (zh) 一种直流电动机量化迭代学习容错控制方法
Xiao et al. Joint optimization of communication rates and linear systems
Wu et al. Sliding mode control of switched hybrid systems with time‐varying delay
CN110815225B (zh) 电机驱动单机械臂***的点对点迭代学习优化控制方法
Shen et al. A probabilistically quantized learning control framework for networked linear systems
CN110989347B (zh) 一种基于事件触发机制的网络化控制***和控制方法
Yang et al. On the feedback capacity of power-constrained Gaussian noise channels with memory
CN113644857A (zh) 一种永磁直线电机有限级量化迭代学习控制方法
CN108803601B (zh) 一种具有通信约束的移动机器人鲁棒预测跟踪控制方法
Wanigasekara et al. Performance of neural network based controllers and ΔΣ-based PID controllers for networked control systems: a comparative investigation
CN113159331B (zh) 一种网络化机器学习***的自适应稀疏度量化方法
CN113110063B (zh) 单轴进给***的鲁棒单调收敛点对点迭代学习控制方法
Eilbrecht et al. Distributed H2-optimized output feedback controller design using the ADMM
Gurt et al. Internal stability of dynamic quantised control for stochastic linear plants
CN115047763B (zh) 一种多无人机***的最小能量控制方法
Østergaard et al. Multiple descriptions for packetized predictive control
CN116893619B (zh) 一种工业机器人量化迭代学习控制方法
Fradkov et al. Hybrid quantised observer for multi-input-multi-output nonlinear systems
CN115107034A (zh) 一种单机械臂的量化迭代学习控制方法
CN115356927B (zh) 机器人的三闭环鲁棒预测函数控制方法
Zhang et al. Predictive compensation based quantization iterative learning control for nonlinear nonaffine discrete‐time systems
Xu et al. Quantized Distributed Nonconvex Optimization Algorithms with Linear Convergence under the Polyak--${\L} $ ojasiewicz Condition
CN116700000B (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