CN113241973B - S型滤波器迭代学习控制直线电机轨迹跟踪控制方法 - Google Patents

S型滤波器迭代学习控制直线电机轨迹跟踪控制方法 Download PDF

Info

Publication number
CN113241973B
CN113241973B CN202110674284.2A CN202110674284A CN113241973B CN 113241973 B CN113241973 B CN 113241973B CN 202110674284 A CN202110674284 A CN 202110674284A CN 113241973 B CN113241973 B CN 113241973B
Authority
CN
China
Prior art keywords
controller
control
linear motor
iterative learning
follows
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110674284.2A
Other languages
English (en)
Other versions
CN113241973A (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.)
Changchun National Extreme Precision Optics Co ltd
Jilin University
Original Assignee
Changchun National Extreme Precision Optics Co ltd
Jilin 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 Changchun National Extreme Precision Optics Co ltd, Jilin University filed Critical Changchun National Extreme Precision Optics Co ltd
Priority to CN202110674284.2A priority Critical patent/CN113241973B/zh
Publication of CN113241973A publication Critical patent/CN113241973A/zh
Application granted granted Critical
Publication of CN113241973B publication Critical patent/CN113241973B/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
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/006Controlling linear motors
    • 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
    • H02P23/00Arrangements or methods for the control of AC motors characterised by a control method other than vector control
    • H02P23/0004Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
    • H02P23/0022Model reference adaptation, e.g. MRAS or MRAC, useful for control or parameter estimation
    • 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
    • H02P23/00Arrangements or methods for the control of AC motors characterised by a control method other than vector control
    • H02P23/14Estimation or adaptation of motor parameters, e.g. rotor time constant, flux, speed, current or voltage
    • 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
    • H02P25/064Linear motors of the synchronous type
    • 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
    • H02P6/00Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
    • H02P6/34Modelling or simulation for control purposes

Landscapes

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

Abstract

一种S型滤波器迭代学习控制直线电机轨迹跟踪控制方法,属于控制工程技术领域。本发明的目的是针对直线电机非线性干扰影响系统控制性能的问题,设计了超前校正的线性自抗扰控制器,然后根据设计好的自抗扰控制器结合迭代学习控制特点的S型滤波器迭代学习控制直线电机轨迹跟踪控制方法。本发明的步骤是:基于系统的控制模型设计超前校正的自抗扰控制器;设计S型结构的滤波器迭代学习控制器与控制器复合的直线电机运动平台轨迹跟踪控制器。本发明极大的提高了永磁同步直线电机平台的控制性能,能够有效提升直线电机的跟踪精度和抗扰能力。

Description

S型滤波器迭代学习控制直线电机轨迹跟踪控制方法
技术领域
本发明属于控制工程技术领域。
背景技术
在现代科技快速发展的前提条件下,从1840年第一台直线电机发明现世至今,直线电机越来越被广泛应用于现代工业领域。得益于其直接传动的性能,在取消了中间传动机构后,大大减小了机械传动时电机系统机械损耗及摩擦力扰动等的影响,但同时也降低了机械结构对系统抵抗外部扰动的能力。因此,为了能够获得更理想的运动控制精度,可以从结构设计和系统控制两方面进行考虑来减小上述问题带来的影响。
发明内容
本发明的目的是针对直线电机非线性干扰影响系统控制性能的问题,设计了超前校正的线性自抗扰控制器,然后根据设计好的自抗扰控制器结合迭代学习控制特点的S型滤波器迭代学习控制直线电机轨迹跟踪控制方法。
本发明的步骤是:
步骤1:基于系统的控制模型设计超前校正的自抗扰控制器;
永磁直线电机的数学模型可简化为
Figure BDA0003120072380000011
式中,系统时间常数T=1/(2.206*2*π),阻尼比ξ=1,s是拉普拉斯变换因子,比例系数KH=12.8;令
Figure BDA0003120072380000012
且设系统中所有的非线性动态和不可测扰动为d,则永磁同步直线电机平台的数学模型时域表达为:
Figure BDA0003120072380000013
Figure BDA0003120072380000014
式中,u为控制输入信号,y为系统输出信号,其中b0是已知输入的控制增益,f(·)为包含系统线性、非线性干扰,内部不确定性和外部扰动的总扰动;
Figure BDA0003120072380000015
Figure BDA0003120072380000016
则系统的状态方程为:
Figure BDA0003120072380000021
其中
Figure BDA0003120072380000022
C=[100],f(·)通过以下的LESO可以进行实时估计:
Figure BDA0003120072380000023
其中,
Figure BDA0003120072380000024
分别是是
Figure BDA0003120072380000025
的估计值,L=[β123]T为其可调增益;
为调参方便,将观测器的极点配置在-w0处,其特征方程为:
λ0(s)=s31s22s+β3=(s+ω0)3 (6)
式中,ω0为观测带宽,单位为rad/s,β1=3ω0,β2=3ω0 2,β3=ω0 3
根据自动控制原理相关知识可知,前向通道上的反馈控制器可以改善系统带宽,从而有效的减小跟踪误差;又知,系统带宽影响了输出跟随输入指令的能力强弱,带宽越宽,系统的动态性能就越好,其对非线性干扰的抑制能力就越强;超前校正能够很大程度上提高系统的相角裕度和剪切频率,但针对低频段的增益较低,难以很准确的跟踪低频信号和抑制干扰;合适的比例积分(PI)环节可以改善超前校正环节的性能,同时提高整个系统的低频段增益以增强其应对低频推力波动的能力,基于此考虑设计PI+超前校正的反馈控制器来对改善系统性能;
设计的PI+超前校正的反馈控制器表达式为:
Figure BDA0003120072380000026
式中的第一项为PI环节,第二项为超前环节;TI=1/ωI,T=1/ωT,τ=1/ωτ;Kc为比例系数,wI,wT,wτ均为转折频率;
控制器的控制律设计为:
Figure BDA0003120072380000027
当观测状态z3较好地估计总扰动f(·)时,系统数学模型表达式可改写为:
Figure BDA0003120072380000031
因此,校正之后的系统开环传递函数为:
Figure BDA0003120072380000032
根据自动控制原理中有关超前滞后控制器的设计原理,可知控制器的相关转折频率满足:
Figure BDA0003120072380000033
Figure BDA0003120072380000034
其中,h为中频段宽度,ωc为剪切频率;
为了更好的系统性能,选取中频段h=9,剪切频率ωc=60Hz,PI控制器的转折频率ωI=0.1Hz,比例系数可以根据频率点之间的关系按下式计算:
Kc=ωI·ωT·ωτ/b0 (13)
设计出的控制器传递函数为:
Figure BDA0003120072380000035
步骤2:设计S型结构的滤波器迭代学习控制器与步骤1的控制器复合的直线电机运动平台轨迹跟踪控制器;
“S”型迭代学习控制前馈量的更新方式可以表示为:
Figure BDA0003120072380000036
其中,LS为迭代学习控制律;QS为滤波器;ek=rk-yk为第k个迭代周期的误差信号;
Figure BDA0003120072380000037
为第k个迭代周期的迭代学习前馈控制信号,第k个迭代周期的控制输出为:
Figure BDA0003120072380000038
式中C为位置环PI+超前校正控制器,P为经扩张状态观测器补偿后的双积分串联型线性系统;
假设迭代过程中初始状态均为零,且
Figure BDA0003120072380000039
则可得出第k个迭代周期的误差传递函数为
Figure BDA00031200723800000310
从而,同样可得出第k+1个迭代周期的误差传递函数为
Figure BDA0003120072380000041
一般而言,位置输入rk不变,即rk+1=rk,根据式(17)和式(18)可知
Figure BDA0003120072380000042
因此,若要使系统收敛,即
Figure BDA0003120072380000043
结合不动点理论,故有
Figure BDA0003120072380000044
由上可知,迭代学习的主要任务是找到满足式(20)的迭代学习律LS,其中滤波器QS可以增强系统的鲁棒性;当
Figure BDA0003120072380000045
时,滤波器的存在能够保证其满足收敛条件,但是增强系统鲁棒性导致跟踪误差不能够趋近于零;
从收敛条件可以推出一种理想学习律,将控制律设计为灵敏度函数
Figure BDA0003120072380000046
的逆,即
Figure BDA0003120072380000047
可以实现对输入信号的完美跟踪;由于许多系统会存在非最小相位零点,直接求逆会导致系统不稳定,因而需要一种稳定的近似模型反演技术;一般而言系统的灵敏度函数为高阶函数,因此可以使用零相位跟踪控制技术来求系统的近似逆,过程如下:
1)通过域间转换将Tu(s)离散化,得到灵敏度函数的z域表达式Tu(z);
2)将系统动力学表达式按式(21)改写,其中B(z)划分为包含稳定(可逆)零点的多项式Bs(z)和包含不稳定(不可逆)Bu(z)零点的多项式:
Figure BDA0003120072380000048
多项式A(z)包含闭环系统的所有极点,Bu(z)可以写成n阶多项式的形式:
Bu(z)=bunzn+bu(n-1)zn-1+...+bu0 (22)
根据收敛条件得到学习律Ls的表达式如下:
Figure BDA0003120072380000051
又由于上式中分子的阶次并不总是等于分母的阶次,因此,可能需要若干个延迟单位来确保学习律L实现;因此
Tu(z)-1=z-pL (24)
3)再将多项式系数按照式的反转,并计算补偿增益:
Figure BDA0003120072380000052
其中
Figure BDA0003120072380000053
4)将式(25)和式(26)代入式(24)中可得学习律Ls的最终表达式;
综上所述,可推导S型滤波器迭代学习补偿方案的学习律表达式为:
Figure BDA0003120072380000054
本发明解决了带宽对系统动态性能的影响,考虑到实际系统中存在的自身非线性及外部扰动难以精确描述,选择了对系统模型依赖较小的自抗扰控制方案,考虑到非线性自抗扰控制调参困难不易用于实际系统,线性自抗扰控制在估计时变扰动时存在相位滞后的问题,进一步考虑超前校正具有提升系统带宽和控制精度的能力,从而设计了基于超前校正的线性自抗扰控制方案,最后,针对直线电机平台具有大量重复性扰动的现象,考虑迭代学习控制的学习能力对抑制重复性干扰的优异性能,设计了S型滤波器迭代学习控制进行前馈补偿。本发明极大的提高了永磁同步直线电机平台的控制性能,能够有效提升直线电机的跟踪精度和抗扰能力。
附图说明:
图1为本发明的系统框架图;
图2为本发明的直线电机平台实物图;
图3为本发明的基于超前校正的自抗扰控制器作用下位移跟踪曲线;
图4为本发明的基于超前校正的自抗扰控制器作用下轨迹跟踪误差曲线;
图5为本发明的加入S型滤波器迭代学习控制器后的位移跟踪曲线;
图6为本发明的加入S型滤波器迭代学习控制器后的轨迹跟踪误差曲线。
具体实施方式:
影响直线电机动态性能的最主要因素是推力波动,而推力波动是一种受电磁结构影响,也与电机状态息息相关的非线性扰动,因而始终存在于直线电机系统中。在对系统进行动力学分析和扫频分析建立直线电机的数学模型后,考虑到系统带宽影响了输出跟随输入指令的能力强弱,带宽越宽,系统的动态性能就越好,其对非线性干扰的抑制能力就越强,设计了基于超前校正的线性自抗扰控制器,然后通过分析系统重复运动时干扰的周期性,结合迭代学习控制的特点,设计了S型滤波器迭代学习控制器对重复性干扰进行补偿。考虑到直线电机在超精密加工领域的巨大优势和广泛应用前景,与其固有推力波动特性对运动控制精度影响之间的矛盾,如何解决这一问题至关重要。
以下结合附图对本发明作进一步描述:
本发明公开一种基于S型滤波器迭代学习控制直线电机轨迹跟踪控制方法,结合系统整体控制框图图1。
具体设计步骤如下:
步骤1:基于直线电机模型设计超前校正的线性自抗扰控制器:
直线电机的数学模型可简化为
Figure BDA0003120072380000061
式中,系统时间常数T=1/(2.206*2*π),阻尼比ξ=1,s是拉普拉斯变换因子,比例系数KH=12.8。
Figure BDA0003120072380000062
且设系统中所有的非线性动态和不可测扰动为d,则直线电机的数学模型时域表达为:
Figure BDA0003120072380000063
Figure BDA0003120072380000064
式中,u为控制输入信号,y为系统输出信号,其中b0是已知输入的控制增益,f(·)为包含系统线性、非线性干扰,内部不确定性和外部扰动的总扰动。
Figure BDA0003120072380000065
则系统的状态方程为:
Figure BDA0003120072380000066
其中
Figure BDA0003120072380000071
C=[1 0 0],f(·)通过以下的LESO可以进行实时估计:
Figure BDA0003120072380000072
其中,
Figure BDA0003120072380000073
分别是是
Figure BDA0003120072380000074
的估计值,L=[β123]T为其可调增益。
为调参方便,将观测器的极点配置在-w0处,其特征方程为:
λ0(s)=s31s22s+β3=(s+ω0)3 (6)
式中,ω0为观测带宽,单位为rad/s,β1=3ω0,β2=3ω0 2,β3=ω0 3
根据自动控制原理相关知识可知,前向通道上的反馈控制器可以改善系统带宽,从而有效的减小跟踪误差。又知,系统带宽影响了输出跟随输入指令的能力强弱,带宽越宽,系统的动态性能就越好,其对非线性干扰的抑制能力就越强。超前校正能够很大程度上提高系统的相角裕度和剪切频率,但针对低频段的增益较低,难以很准确的跟踪低频信号和抑制干扰。合适的比例积分(PI)环节可以改善超前校正环节的性能,同时提高整个系统的低频段增益以增强其应对低频推力波动的能力,基于此考虑设计PI+超前校正的反馈控制器来对改善系统性能。
设计的PI+超前校正的反馈控制器表达式为:
Figure BDA0003120072380000075
式中的第一项为PI环节,第二项为超前环节。TI=1/ωI,T=1/ωT,τ=1/ωτ。Kc为比例系数,wI,wT,wτ均为转折频率。
控制器的控制律设计为:
Figure BDA0003120072380000076
当观测状态z3较好地估计总扰动f(·)时,系统数学模型表达式可改写为:
Figure BDA0003120072380000077
因此,校正之后的系统开环传递函数为:
Figure BDA0003120072380000081
根据自动控制原理中有关超前滞后控制器的设计原理,可知控制器的相关转折频率满足:
Figure BDA0003120072380000082
Figure BDA0003120072380000083
其中,h为中频段宽度,ωc为剪切频率。
为了更好的系统性能,选取中频段h=9,剪切频率ωc=60Hz,PI控制器的转折频率ωI=0.1Hz,比例系数可以根据频率点之间的关系按下式计算:
Kc=ωI·ωT·ωτ/b0 (13)
设计出的控制器传递函数为:
Figure BDA0003120072380000084
步骤二:在步骤一的基础上设计S型滤波器迭代学习控制器:
“S”型迭代学习控制前馈量的更新方式可以表示为:
Figure BDA0003120072380000085
其中,LS为迭代学习控制律;QS为滤波器;ek=rk-yk为第k个迭代周期的误差信号;
Figure BDA0003120072380000086
为第k个迭代周期的迭代学习前馈控制信号,第k个迭代周期的控制输出为:
Figure BDA0003120072380000087
式中C为位置环PI+超前校正控制器,P为经扩张状态观测器补偿后的双积分串联型线性系统。
假设迭代过程中初始状态均为零,且
Figure BDA0003120072380000088
则可得出第k个迭代周期的误差传递函数为
Figure BDA0003120072380000089
从而,同样可得出第k+1个迭代周期的误差传递函数为
Figure BDA0003120072380000091
一般而言,位置输入rk不变,即rk+1=rk,根据式(17)和式(18)可知
Figure BDA0003120072380000092
因此,若要使系统收敛,即
Figure BDA0003120072380000093
结合不动点理论,故有
Figure BDA0003120072380000094
由上可知,迭代学习的主要任务是找到满足式(20)的迭代学习律LS,其中滤波器QS可以增强系统的鲁棒性。当
Figure BDA0003120072380000095
时,滤波器的存在能够保证其满足收敛条件,但是增强系统鲁棒性导致跟踪误差不能够趋近于零。
从收敛条件可以推出一种理想学习律,将控制律设计为灵敏度函数
Figure BDA0003120072380000096
的逆,即
Figure BDA0003120072380000097
可以实现对输入信号的完美跟踪。由于许多系统会存在非最小相位零点,直接求逆会导致系统不稳定,因而需要一种稳定的近似模型反演技术。
一般而言系统的灵敏度函数为高阶函数,因此可以使用零相位跟踪控制技术来求系统的近似逆,过程如下:
1)通过域间转换将Tu(s)离散化,得到灵敏度函数的z域表达式Tu(z)。
2)将系统动力学表达式按式(21)改写,其中B(z)划分为包含稳定(可逆)零点的多项式Bs(z)和包含不稳定(不可逆)Bu(z)零点的多项式:
Figure BDA0003120072380000098
多项式A(z)包含闭环系统的所有极点,Bu(z)可以写成n阶多项式的形式:
Bu(z)=bunzn+bu(n-1)zn-1+...+bu0 (22)
根据收敛条件得到学习律Ls的表达式如下:
Figure BDA0003120072380000101
又由于上式中分子的阶次并不总是等于分母的阶次,因此,可能需要若干个延迟单位来确保学习律L实现。因此
Tu(z)-1=z-pL (24)
3)再将多项式系数按照式的反转,并计算补偿增益:
Figure BDA0003120072380000102
其中
Figure BDA0003120072380000103
4)将式(25)和式(26)代入式(24)中可得学习律Ls的最终表达式。
综上所述,可推导S型滤波器迭代学习补偿方案的学习律表达式为:
Figure BDA0003120072380000104
试验:
搭建如图2所示的直线电机平台进行实验,验证基于S型滤波器迭代学习直线电机轨迹跟踪控制方法的有效性,得到最终实验结果。
直线电机的各项理论参数如表1所示:
表1直线电机各项理论参数
Figure BDA0003120072380000105
直线电机平台的运动部分主要组成为ELMO驱动器驱动的一维直线电机平台,其行程为100mm。位移传感器为分辨率为50nm的增量式光栅。首先通过光栅采集直线电机位置信号,将其与期望的位置信号比较得到误差信号,将仿真验证过的算法改写成PMAC控制器的上位机软件可识别的C语言程序,然后调用该算法,最后利用PMAC控制器中的D/A转化模块输出控制电压信号到驱动器中驱动电机进行直线运动。由于实际工况中的直线运动多为由多个S型曲线构成的步进轨迹,因此本课题的参考曲线选为三阶S曲线。同时,为了保证各控制方案的对比性,选择设计S曲线参数固定,其参数设置如下:加加速度为J=5m/s3,加速度为A=0.4m/s2,速度为V=0.015m/s,位移为X=10mm,非线性干扰用幅值范围为±50μm的白噪声替代。
下面验证所提的S型滤波器迭代学习前馈补偿控制的效果,分别采取单独的基于超前校正的线性自抗扰控制器和具有S型滤波器迭代学习前馈补偿的超前校正的线性自抗扰复合控制器进行比较。
基于超前校正的线性自抗扰控制器的各参数取值为b0=2000,w0=150(rad·s-1),wc=60Hz,wI=0.001Hz,h=9。
滤波器选用带宽为100的低通滤波器
Figure BDA0003120072380000111
其离散表达式为
Figure BDA0003120072380000112
图3和图4分别表示的是单独的基于超前校正的线性自抗扰控制器跟踪S型曲线的位移跟踪曲线和轨迹跟踪误差曲线,图5和图6分别表示的是具有S型滤波器迭代学习前馈补偿的复合控制器在迭代次数为1、2、5、10次的跟踪S型曲线的位移跟踪曲线和轨迹跟踪误差曲线。
表2为上述情况下的最大误差(maximum error,MAX)、均方根误差(root meansquare error,RMSE)、平均绝对误差(mean absolute error,MAE)的性能指标对比结果。
MAX=max|yi-yr| (28)
Figure BDA0003120072380000113
Figure BDA0003120072380000114
其中yr为参考信号,yi为系统输出,N为数据总数。
表2两种控制方法的性能指标对比
Figure BDA0003120072380000121

Claims (1)

1.一种S型滤波器迭代学习控制直线电机轨迹跟踪控制方法,其特征在于:其步骤是:
步骤1:基于系统的控制模型设计超前校正的自抗扰控制器;
永磁直线电机的数学模型可简化为
Figure FDA0003120072370000011
式中,系统时间常数T=1/(2.206*2*π),阻尼比ξ=1,s是拉普拉斯变换因子,比例系数KH=12.8;令
Figure FDA0003120072370000012
且设系统中所有的非线性动态和不可测扰动为d,则永磁同步直线电机平台的数学模型时域表达为:
Figure FDA0003120072370000013
Figure FDA0003120072370000014
式中,u为控制输入信号,y为系统输出信号,其中b0是已知输入的控制增益,f(·)为包含系统线性、非线性干扰,内部不确定性和外部扰动的总扰动;
Figure FDA0003120072370000015
则系统的状态方程为:
Figure FDA0003120072370000016
其中
Figure FDA0003120072370000017
C=[1 0 0],f(·)通过以下的LESO可以进行实时估计:
Figure FDA0003120072370000018
其中,
Figure FDA0003120072370000019
Figure FDA00031200723700000110
分别是y,
Figure FDA00031200723700000111
f(·)的估计值,L=[β123]T为其可调增益;
为调参方便,将观测器的极点配置在-w0处,其特征方程为:
λ0(s)=s31s22s+β3=(s+ω0)3 (6)
式中,ω0为观测带宽,单位为rad/s,β1=3ω0,β2=3ω0 2,β3=ω0 3
根据自动控制原理相关知识可知,前向通道上的反馈控制器可以改善系统带宽,从而有效的减小跟踪误差;又知,系统带宽影响了输出跟随输入指令的能力强弱,带宽越宽,系统的动态性能就越好,其对非线性干扰的抑制能力就越强;超前校正能够很大程度上提高系统的相角裕度和剪切频率,但针对低频段的增益较低,难以很准确的跟踪低频信号和抑制干扰;合适的比例积分(PI)环节可以改善超前校正环节的性能,同时提高整个系统的低频段增益以增强其应对低频推力波动的能力,基于此考虑设计PI+超前校正的反馈控制器来对改善系统性能;
设计的PI+超前校正的反馈控制器表达式为:
Figure FDA0003120072370000021
式中的第一项为PI环节,第二项为超前环节;TI=1/ωI,T=1/ωT,τ=1/ωτ;Kc为比例系数,wI,wT,wτ均为转折频率;
控制器的控制律设计为:
Figure FDA0003120072370000022
当观测状态z3较好地估计总扰动f(·)时,系统数学模型表达式可改写为:
Figure FDA0003120072370000023
因此,校正之后的系统开环传递函数为:
Figure FDA0003120072370000024
根据自动控制原理中有关超前滞后控制器的设计原理,可知控制器的相关转折频率满足:
Figure FDA0003120072370000025
Figure FDA0003120072370000026
其中,h为中频段宽度,ωc为剪切频率;
为了更好的系统性能,选取中频段h=9,剪切频率ωc=60Hz,PI控制器的转折频率ωI=0.1Hz,比例系数可以根据频率点之间的关系按下式计算:
Kc=ωI·ωT·ωτ/b0 (13)
设计出的控制器传递函数为:
Figure FDA0003120072370000031
步骤2:设计S型结构的滤波器迭代学习控制器与步骤1的控制器复合的直线电机运动平台轨迹跟踪控制器;
“S”型迭代学习控制前馈量的更新方式可以表示为:
Figure FDA0003120072370000032
其中,LS为迭代学习控制律;QS为滤波器;ek=rk-yk为第k个迭代周期的误差信号;
Figure FDA00031200723700000311
为第k个迭代周期的迭代学习前馈控制信号,第k个迭代周期的控制输出为:
Figure FDA0003120072370000033
式中C为位置环PI+超前校正控制器,P为经扩张状态观测器补偿后的双积分串联型线性系统;
假设迭代过程中初始状态均为零,且
Figure FDA0003120072370000034
则可得出第k个迭代周期的误差传递函数为
Figure FDA0003120072370000035
从而,同样可得出第k+1个迭代周期的误差传递函数为
Figure FDA0003120072370000036
一般而言,位置输入rk不变,即rk+1=rk,根据式(17)和式(18)可知
Figure FDA0003120072370000037
因此,若要使系统收敛,即
Figure FDA0003120072370000038
结合不动点理论,故有
Figure FDA0003120072370000039
由上可知,迭代学习的主要任务是找到满足式(20)的迭代学习律LS,其中滤波器QS可以增强系统的鲁棒性;当
Figure FDA00031200723700000310
时,滤波器的存在能够保证其满足收敛条件,但是增强系统鲁棒性导致跟踪误差不能够趋近于零;
从收敛条件可以推出一种理想学习律,将控制律设计为灵敏度函数
Figure FDA0003120072370000041
的逆,即
Figure FDA0003120072370000042
可以实现对输入信号的完美跟踪;由于许多系统会存在非最小相位零点,直接求逆会导致系统不稳定,因而需要一种稳定的近似模型反演技术;一般而言系统的灵敏度函数为高阶函数,因此可以使用零相位跟踪控制技术来求系统的近似逆,过程如下:
1)通过域间转换将Tu(s)离散化,得到灵敏度函数的z域表达式Tu(z);
2)将系统动力学表达式按式(21)改写,其中B(z)划分为包含稳定(可逆)零点的多项式Bs(z)和包含不稳定(不可逆)Bu(z)零点的多项式:
Figure FDA0003120072370000043
多项式A(z)包含闭环系统的所有极点,Bu(z)可以写成n阶多项式的形式:
Bu(z)=bunzn+bu(n-1)zn-1+...+bu0 (22)
根据收敛条件得到学习律Ls的表达式如下:
Figure FDA0003120072370000044
又由于上式中分子的阶次并不总是等于分母的阶次,因此,可能需要若干个延迟单位来确保学习律L实现;因此
Tu(z)-1=z-pL (24)
3)再将多项式系数按照式的反转,并计算补偿增益:
Figure FDA0003120072370000045
其中
Figure FDA0003120072370000046
4)将式(25)和式(26)代入式(24)中可得学习律Ls的最终表达式;
综上所述,可推导S型滤波器迭代学习补偿方案的学习律表达式为:
Figure FDA0003120072370000051
CN202110674284.2A 2021-06-17 2021-06-17 S型滤波器迭代学习控制直线电机轨迹跟踪控制方法 Active CN113241973B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110674284.2A CN113241973B (zh) 2021-06-17 2021-06-17 S型滤波器迭代学习控制直线电机轨迹跟踪控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110674284.2A CN113241973B (zh) 2021-06-17 2021-06-17 S型滤波器迭代学习控制直线电机轨迹跟踪控制方法

Publications (2)

Publication Number Publication Date
CN113241973A CN113241973A (zh) 2021-08-10
CN113241973B true CN113241973B (zh) 2022-07-05

Family

ID=77140256

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110674284.2A Active CN113241973B (zh) 2021-06-17 2021-06-17 S型滤波器迭代学习控制直线电机轨迹跟踪控制方法

Country Status (1)

Country Link
CN (1) CN113241973B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113934140A (zh) * 2020-07-14 2022-01-14 无锡市电子仪表工业有限公司 一种基于分布式滤波器的迭代学习目标轨迹跟踪方法
CN113759721B (zh) * 2021-09-02 2023-07-14 大连交通大学 一种永磁同步电动机l1自适应迭代学习控制方法
CN114185274B (zh) * 2021-12-06 2023-07-04 东北大学 基于迭代学习的钢铁生产过程重复性误差补偿控制方法
CN114265314B (zh) * 2021-12-23 2022-06-24 哈尔滨工业大学 一种基于fir滤波的鲁棒逆模型学习增益设计方法
CN115562017B (zh) * 2022-10-07 2024-08-23 中国计量科学研究院 一种传感器低频振动校准系统的迭代学习控制方法
CN117031967B (zh) * 2023-10-10 2024-01-23 华南理工大学 一种迭代学习自抗扰控制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109015661A (zh) * 2018-09-29 2018-12-18 重庆固高科技长江研究院有限公司 工业机器人迭代学习修正轨迹误差的方法
CN111817638A (zh) * 2020-06-10 2020-10-23 吉林大学 永磁同步直线电机平台的相位超前线性自抗扰控制器
CN113759721A (zh) * 2021-09-02 2021-12-07 大连交通大学 一种永磁同步电动机l1自适应迭代学习控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109015661A (zh) * 2018-09-29 2018-12-18 重庆固高科技长江研究院有限公司 工业机器人迭代学习修正轨迹误差的方法
CN111817638A (zh) * 2020-06-10 2020-10-23 吉林大学 永磁同步直线电机平台的相位超前线性自抗扰控制器
CN113759721A (zh) * 2021-09-02 2021-12-07 大连交通大学 一种永磁同步电动机l1自适应迭代学习控制方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
High-Procision Motion Control for a Linear Permanent Magnet Iron Core Synchronous Motor Drive in Position Platform;Mei-Yung Chen and Jian-Shiun Lu;《IEEE TRANSACTIONS ON INDUSTRIAL INFORMATICS》;20140228;第10卷(第1期);全文 *
基于Smith预估和性能加权函数的永磁直线同步电机鲁棒迭代学习控制;赵希梅等;《电工技术学报》;20161010(第19期);全文 *
基于迭代学习与FIR滤波器的PMLSM高精密控制;赵希梅等;《电工技术学报》;20170510;第32卷(第09期);全文 *

Also Published As

Publication number Publication date
CN113241973A (zh) 2021-08-10

Similar Documents

Publication Publication Date Title
CN113241973B (zh) S型滤波器迭代学习控制直线电机轨迹跟踪控制方法
CN108303885B (zh) 一种基于干扰观测器的电机位置伺服系统自适应控制方法
CN107544241B (zh) 压电陶瓷执行器迟滞的非线性pid逆补偿控制方法
Godbole et al. Performance analysis of generalized extended state observer in tackling sinusoidal disturbances
Lan et al. Optimal preview repetitive control with application to permanent magnet synchronous motor drive system
Jakovljević et al. Control of double-loop permanent magnet synchronous motor drives by optimized fractional and distributed-order PID controllers
CN107577146B (zh) 基于摩擦整体逼近的伺服系统的神经网络自适应控制方法
CN104199294A (zh) 电机伺服系统双神经网络摩擦补偿和有限时间协同控制方法
CN113517832B (zh) 一种低压伺服离散线性自抗扰控制方法
CN111708276B (zh) 基于线性状态观测器观测误差补偿的自适应鲁棒控制方法
CN110649845B (zh) 基于鲁棒广义预测控制的光电转台位置跟踪控制方法
CN111142404A (zh) 一种基于压电陶瓷驱动的微定位平台及其建模与控制方法
CN109634116B (zh) 一种分数阶机械式离心调速器系统的加速自适应稳定方法
CN110442026B (zh) 基于误差修正的扩张状态观测器、抗扰控制系统及设计方法
Hu et al. High-order sliding-mode observer based output feedback adaptive robust control of a launching platform with backstepping
CN113283090B (zh) 基于蜂群算法的飞机舵机电动伺服系统摩擦补偿方法
Huo et al. A compound control method for the rejection of spatially periodic and uncertain disturbances of rotary machines and its implementation under uniform time sampling
Youcef-Toumi et al. Control of robot manipulators using time delay
Yu et al. Observer-based data-driven constrained norm optimal iterative learning control for unknown non-affine non-linear systems with both available and unavailable system states
CN107544250B (zh) 一种少保守性的永磁球形电机轨迹跟踪控制方法
CN109995278B (zh) 一种考虑输入受限的电机伺服系统自调节控制方法
CN110471282A (zh) 一种同时估计系统状态和扰动输入的输出反馈控制方法
CN109194244B (zh) 一种面向电动伺服系统的控制方法及系统
Zhou et al. Preisach modeling of hysteresis and tracking control of a Thunder actuator system
Ma et al. Adaptive actuator compensation control with feedback linearization

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