CN115473467A - 基于模糊观测器的柔性关节机械臂指令滤波反步控制方法 - Google Patents

基于模糊观测器的柔性关节机械臂指令滤波反步控制方法 Download PDF

Info

Publication number
CN115473467A
CN115473467A CN202211123023.2A CN202211123023A CN115473467A CN 115473467 A CN115473467 A CN 115473467A CN 202211123023 A CN202211123023 A CN 202211123023A CN 115473467 A CN115473467 A CN 115473467A
Authority
CN
China
Prior art keywords
formula
equation
mechanical arm
representing
flexible joint
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.)
Pending
Application number
CN202211123023.2A
Other languages
English (en)
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.)
Qingdao University
Original Assignee
Qingdao 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 Qingdao University filed Critical Qingdao University
Priority to CN202211123023.2A priority Critical patent/CN115473467A/zh
Publication of CN115473467A publication Critical patent/CN115473467A/zh
Pending legal-status Critical Current

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
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, 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
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/0003Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
    • H02P21/001Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control using fuzzy control
    • 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
    • H02P21/13Observer control, e.g. using Luenberger observers or Kalman filters
    • 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
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • H02P21/18Estimation of position or speed

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Fuzzy Systems (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明属于柔性关节机械臂位置跟踪控制技术领域,具体公开了一种基于模糊观测器的柔性关节机械臂指令滤波反步控制方法,该方法针对模型参数未知的柔性关节机械臂系统,利用指令滤波及模糊观测器技术提出了指令滤波模糊控制方法,实现了对该系统期望轨迹的跟踪控制。设计模糊观测器来估计柔性机械臂的连杆角速度和电机转角速度;利用模糊自适应技术来解决模型参数不同程度的不确定性;同时,应用指令滤波技术解决了控制器设计过程中的计算复杂性问题;并通过Lyapunov控制原理证明了该柔性关节机械臂控制系统中所有变量的收敛性。本发明设计的观测器具有良好的观测效果,且解决了控制器设计过程中的计算复杂性问题,轨迹跟踪误差更小,系统鲁棒性更强。

Description

基于模糊观测器的柔性关节机械臂指令滤波反步控制方法
技术领域
本发明属于柔性关节机械臂位置跟踪控制技术领域,特别涉及一种基于模糊观测器的柔性关节机械臂指令滤波反步控制方法。
背景技术
近年来,由于机器人技术的迅猛发展,柔性关节机械臂在航空航天、工件装配、国防军工和教育行业等方面广泛应用,结合发展迅速的机器视觉技术,其在制造和医疗等行业均具有较高的研究价值。传统的刚性机械臂难以解决大量繁杂的高精度控制问题,大量研究和实践表明,在设计过程中若忽略机械臂的关节柔性,此机械臂在执行高精度作业时将受到很大的限制。相较于传统的刚性机械臂,柔性机械臂在搭载了谐波减速器的基础上有较为明显的连续性和多自由度数目的特点,其在面对复杂多障碍的环境情况、精细化的操作要求时,动作灵活精准,关节扭转范围宽泛。而柔性机械臂系统是强耦合的高阶非线性系统,这些特点限制了机械臂的控制性能。因此,柔性机械臂的精确控制问题,已是当今机器人领域的研究热点。控制方法包括奇异摄动法、智能控制、变结构控制、反馈线性化法、反步控制等。
其中反步控制是一种基于递归李雅普诺夫函数的方案,是一个有效的非线性系统控制框架,以此方法为基础在处理非线性控制方面的问题上有独特的优越性,反步控制在柔性机械臂的控制领域应用广泛。然而,在对虚拟控制函数连续求导时,控制器设计过程中所用到的传统反步法会存在“计算爆炸”问题,进而增加了计算量,导致柔性机械臂反步控制策略的应用具有很大的局限性。后来提出了动态面控制方法来解决“计算爆炸”,但是动态面控制存在的滤波误差会影响系统的控制精度。进一步提出了指令滤波技术,并结合误差补偿技术降低了滤波误差对控制精度的干扰。另外在一些控制方法中,要想实现对运行中柔性机械臂系统的更精确控制,需要准确了解运行系统中的各个状态变量,但是有一部分状态变量是不可测的,当想要观测柔性机械臂的速度项时就需要用到传感器,但传感器测量时易受干扰。
综上,目前柔性关节机械臂系统中存在未知项,连杆速度及电机转角速度状态不可测,关节位置跟踪在反步设计中推导复杂以及误差过大引起安全事故等技术问题。
发明内容
针对现有技术中存在的上述技术问题,本发明提出了一种基于模糊观测器的柔性关节机械臂指令滤波反步控制方法,以便对柔性关节机械臂系统进行精确位置跟踪控制。
本发明为了实现上述目的,采用如下技术方案:
步骤1.建立考虑不确定性的柔性关节机械臂动力学系统数学模型;
步骤2.设计模糊观测器估计得到的柔性关节机械臂模型中不可测的连杆速度和电机转角速度,并用模糊逻辑系统处理系统中不同程度的不确定项;
步骤3.设计指令滤波反步控制器,构建了二阶滤波器和误差补偿信号,用来解决反步设计中虚拟控制律的计算复杂度问题;
步骤4.选取Lyapunov函数进行推导,进而证明对由步骤3设计的基于模糊观测器的柔性关节机械臂指令滤波控制方法所控制系统Lyapunov稳定。
本发明具有如下优点:
(1)本发明方法设计了观测器估计柔性机械臂的连杆速度和电机转角速度,解决了机械臂速度项不可测问题,提高了系统的可靠性。
(2)本发明方法结合了指令滤波技术和误差补偿机制,解决了在控制器设计过程中的计算复杂性问题,消除了滤波误差的影响,不用再考虑期望信号高阶可导的需求,获得了更好的跟踪效果,提高了系统的控制精度。
(3)本发明方法考虑了系统中的不确定因素,利用模糊逻辑系统处理柔性机械臂系统的不确定性,提高了系统鲁棒性。
(4)本发明方法设计的观测器避免了实际系统中传感器的使用安装且减少了维护要求,解决了速度传感器测量时易受干扰的问题,更有利于实际应用。
附图说明
图1为本发明实施例中系统总体控制器以及系统复合被控对象的整体示意图。
图2为采用本发明方法后关节角位置和期望角位置的跟踪响应曲线图。
图3为采用本发明方法后关节角位置和期望角位置跟踪误差的响应曲线图。
图4为采用本发明方法后模糊观测器估计连杆速度项的跟踪曲线图。
图5为采用本发明方法后模糊观测器估计电机转角速度项的跟踪曲线图。
图6为采用本发明方法后该系统控制输入响应曲线图。
图7为采用本发明方法后应对不同程度不确定性B时位置的跟踪误差响应曲线图。
图8为采用本发明方法后应对不同程度不确定性
Figure BDA0003847211180000021
时位置的跟踪误差响应曲线图。
具体实施方式
本发明实施例述及了一种基于模糊观测器的柔性关节机械臂指令滤波反步控制方法,该方法针对模型参数未知的柔性关节机械臂系统,利用指令滤波及模糊观测器技术,提出了指令滤波模糊控制策略,实现了对该系统期望轨迹的跟踪控制,其包括如下步骤:首先设计模糊观测器来估计柔性机械臂的连杆角速度和电机转角速度;其次利用模糊自适应技术来解决模型参数的不确定性;同时应用指令滤波技术解决了控制器设计过程中的计算复杂性问题;并通过Lyapunov控制原理证明了该柔性关节机械臂控制系统中所有变量的收敛性。
下面结合附图以及具体实施方式对本发明作进一步详细说明:
图1为系统总体控制器以及柔性关节机械臂系统复合被控对象的整体示意图,其中,图1中涉及到指令滤波反步控制器、模糊观测器部分以及滤波补偿机制。
本发明方法的控制过程如下:对于柔性机械臂系统的关节位置、速度以及电机侧的角位置、速度四个输出量,输入控制器中控制信号;由于速度项不可测,输入到模糊观测器中产生速度项的估计,再将速度项的估计送入到控制器中;与给定目标位置作比较后产生的差值送入控制器,并将关节角位置与电机转角位置的估计误差送入控制器;将指令滤波反步控制器的控制信号输入到柔性机械臂系统中进而达成对系统的控制。
基于模糊观测器的柔性关节机械臂指令滤波反步控制方法,具体包括如下步骤:
步骤1.建立考虑不确定项的柔性关节机械臂动力学系统数学模型。
首先给出柔性关节机械臂的动态方程,如公式(1)所示:
Figure BDA0003847211180000031
其中,q、
Figure BDA0003847211180000032
分别表示连杆的位置、速度和加速度,qm
Figure BDA0003847211180000033
分别表示电机的转角位置、速度和加速度,M代表连杆的质量,l代表连杆的长度,g为重力加速度,
Figure BDA0003847211180000034
代表摩擦项即未知的干扰项,K表示关节刚度系数,J和B分别表示执行器的惯性和自然阻尼项,u代表电机控制输入,y表示系统的输出;设定
Figure BDA0003847211180000035
和B中包含不确定性。
设状态变量x1=q,
Figure BDA0003847211180000036
x3=qm,
Figure BDA0003847211180000037
将动态方程即公式(1)改写为:
Figure BDA0003847211180000038
为简化动态方程即公式(1),接下来定义如下变量,如公式(3)所示:
Figure BDA0003847211180000039
基于公式(3),进一步将公式(2)改写成:
Figure BDA0003847211180000041
设h2和h4是已知正常数,f2和f4是含有未知项的连续函数,且有连续有界的一阶导数。步骤2.设计模糊观测器,估计柔性关节机械臂的连杆速度和电机的转角速度。
令公式(4)中
Figure BDA0003847211180000042
进一步得到公式(5):
Figure BDA0003847211180000043
其中,Z1=[x1,x2,x3,x4]T
Figure BDA0003847211180000044
由模糊逻辑系统定义得知:
Figure BDA0003847211180000045
其中,θi是最优参数向量,
Figure BDA0003847211180000046
表示θi的估计,
Figure BDA0003847211180000047
为xi的状态估计值。
Figure BDA0003847211180000048
表示
Figure BDA0003847211180000049
与θi同为函数变量,
Figure BDA00038472111800000410
表示
Figure BDA00038472111800000411
的估计,
Figure BDA00038472111800000412
表示变量为
Figure BDA00038472111800000413
的函数块,
Figure BDA00038472111800000414
表示
Figure BDA00038472111800000415
的估计;
Figure BDA00038472111800000416
Figure BDA00038472111800000417
表示状态变量
Figure BDA00038472111800000418
与θi共同简化成的函数,
Figure BDA00038472111800000419
表示
Figure BDA00038472111800000420
的估计。
针对一个任意给定的正数ε,由模糊自适应逼近定理得到:
Figure BDA00038472111800000421
其中,
Figure BDA00038472111800000422
表示函数中变量为Z1且系数为θi T的函数块,δi(Z1)≤|ε|。
基于公式(5)并综合上式
Figure BDA00038472111800000423
将公式(4)的观测器设计为:
Figure BDA00038472111800000424
式中,
Figure BDA0003847211180000051
Figure BDA0003847211180000052
表示
Figure BDA0003847211180000053
的估计;定义观测器的误差e=[e1,e2,e3,e4]T,其中,ej表示估计误差,j=1,2,3,4;估计误差
Figure BDA0003847211180000054
D1、D2、D3、D4为待设计的正数;得到
Figure BDA0003847211180000055
Figure BDA0003847211180000056
其中,
Figure BDA0003847211180000057
表示
Figure BDA0003847211180000058
的估计,
Figure BDA0003847211180000059
表示
Figure BDA00038472111800000510
的估计。
假定一个函数fi(·)满足以下不等式:
Figure BDA00038472111800000511
Figure BDA00038472111800000512
为正常数,i=2,4;得到
Figure BDA00038472111800000513
其中,
Figure BDA00038472111800000514
表示Z1的估计。
Figure BDA00038472111800000515
其中,
Figure BDA00038472111800000516
是不确定正常数;
Figure BDA00038472111800000517
表示
Figure BDA00038472111800000518
的估计。
Figure BDA00038472111800000519
得到
Figure BDA00038472111800000520
ε表示一个任意小的正数。
设有一对称矩阵Q=QT>0,那么会有一对称矩阵P=PT>0满足ATP+PA=-Q。
步骤3.根据指令滤波和反步法设计基于模糊观测器的柔性关节机械臂的指令滤波反步控制方法,具体过程如下:
定义指令滤波器的公式如下:
Figure BDA00038472111800000521
其中,zo,1表示指令滤波器的一阶变量,zo,2表示指令滤波器的二阶变量。
xo,c=zo,1
Figure BDA00038472111800000522
为指令滤波器的两个输出;αo为指令滤波器的输入,αo表示虚拟控制函数;指令滤波器的初始状态αo(0)=zo,1(0),zo,2(0)=0。
当式(4)系统开始运行后,若输入信号满足:
Figure BDA00038472111800000523
则任意μ>0,存在ζ∈(0,1]和ωn≥0,使得|xo,co|≤μ,o=1,2,3,
Figure BDA00038472111800000524
均有界。
ρ1表示虚拟控制函数一阶导的上界,ρ2表示虚拟控制函数二阶导的上界。
由反步法的原理定义误差变量:
Figure BDA0003847211180000061
其中,zj表示误差变量,j=1,2,3,4;xd为给定的期望信号,xo,c为指令滤波器的输出信号,
Figure BDA0003847211180000062
为速度项的估计值。
定义误差补偿信号为:
Figure BDA0003847211180000063
其中,ξj表示误差补偿信号,vj表示补偿误差跟踪信号。
虚拟控制律和滤波误差补偿信号的具体形式将在以下设计过程中给出。
步骤3.1.选取Lyapunov函数为V0=eTPe,求导得:
Figure BDA0003847211180000064
由杨氏不等式得到:
2eTPF≤||e||2+(m2 2+m4 2)||P2||e||2 (11)
Figure BDA0003847211180000065
其中,m2、m4表示任意小的正数;将公式(11)、公式(12)代入公式(10)得:
Figure BDA0003847211180000066
其中,λmin(Q)是Q最小的特征值。
步骤3.2.选取第一个子系统的Lyapunov函数
Figure BDA0003847211180000067
求导得:
Figure BDA0003847211180000068
由于e2未知,由杨氏不等式得到:
Figure BDA0003847211180000069
设计虚拟控制律α1和补偿信号
Figure BDA00038472111800000610
为:
Figure BDA0003847211180000071
其中,k1为正常数,将公式(15)、公式(16)代入公式(14)得:
Figure BDA0003847211180000072
步骤3.3.选取Lyapunov函数
Figure BDA0003847211180000073
求导得:
Figure BDA0003847211180000074
其中,
Figure BDA0003847211180000075
为θ2的估计,
Figure BDA0003847211180000076
表示自适应律,λ2为正常数,h2为正常数。由杨氏不等式得到:
Figure BDA0003847211180000077
将公式(19)代入公式(18)得:
Figure BDA0003847211180000078
设计虚拟控制律α2、补偿信号
Figure BDA0003847211180000079
和自适应律
Figure BDA00038472111800000710
为:
Figure BDA00038472111800000711
其中,k2为正常数,n2为正常数,将公式(21)代入公式(20)得:
Figure BDA00038472111800000712
步骤3.4.选取Lyapunov函数
Figure BDA00038472111800000713
求导得:
Figure BDA0003847211180000081
由于e4未知,由杨氏不等式得到:
Figure BDA0003847211180000082
设计虚拟控制律α3和补偿信号
Figure BDA0003847211180000083
为:
Figure BDA0003847211180000084
其中,k3为正常数,将公式(24)、公式(25)代入公式(23)得:
Figure BDA0003847211180000085
步骤3.5.选取Lyapunov函数
Figure BDA0003847211180000086
求导得:
Figure BDA0003847211180000087
其中,λ4为正常数;
Figure BDA0003847211180000088
Figure BDA0003847211180000089
为θ4的估计,
Figure BDA00038472111800000810
表示自适应律。由杨氏不等式得到:
Figure BDA00038472111800000811
将公式(28)代入公式(27)得:
Figure BDA00038472111800000812
设计真实控制律u、补偿信号
Figure BDA00038472111800000813
和自适应律
Figure BDA00038472111800000814
为:
Figure BDA0003847211180000091
其中,k4、n4均为正的常数,将公式(30)代入公式(29)得:
Figure BDA0003847211180000092
步骤4.对基于模糊观测器的柔性关节机械臂指令滤波反步控制方法进行稳定性分析。选取V=V4,由公式(31)得:
Figure BDA0003847211180000093
由杨氏不等式得:
Figure BDA0003847211180000094
Figure BDA0003847211180000095
将公式(13)、公式(33)、公式(34)代入(32)得到:
Figure BDA0003847211180000096
其中,λmin(Q)-3-(m2 2+m4 2)||P||2>0,
Figure BDA0003847211180000097
Figure BDA0003847211180000098
Figure BDA0003847211180000099
其中,λmax(P)表示向量P的最大特征值;进一步由公式(35)得到:
Figure BDA00038472111800000910
式(36)表明,vj
Figure BDA0003847211180000101
均属于紧集:
Figure BDA0003847211180000102
因此,机械臂系统的所有信号都有界,滤波器的误差补偿信号ξj满足
Figure BDA0003847211180000103
其中
Figure BDA0003847211180000104
因为z1=v11且ξ1有界,跟踪误差z1有界,进而得到:
Figure BDA0003847211180000105
其中,minj(kj)表示kj中的最小值。
通过以上分析得知,vj、ξj、zj都是有界的,因此,在给定参数kj、mi、ni后,保持a不变,选择足够大的λi和足够小的ε,能够保证跟踪误差z1足够小。
为了验证所提控制方法的有效性,柔性机械臂控制系统选取以下参数进行仿真:
M=0.25kg,L=0.45m,g=9.8m/s2,K=5N·m/rad,J=5×10-4m/s2,B=0.01,B1=0.01。
选取摩擦项为
Figure BDA0003847211180000106
选取参考信号xd=sin(t)。
选择控制器参数为:k1=k2=k3=k4=10,b2=b4=0.5,n2=n4=3。
选择观测器参数为:D1=D3=100,D2=D4=10000,选择滤波器参数ωn=1100,ζ=0.8。
同时为了证明所提运用模糊逻辑系统对处理系统中不确定项的有效性,在其余参数不变的前提下,分别改变具有不确定性的摩擦项
Figure BDA0003847211180000109
和B进行对比仿真。
当B=0.01,分别取
Figure BDA0003847211180000107
Figure BDA0003847211180000108
分别取B=0.01;0.05;0.1。仿真结果,如图2至图8所示:
图2反映了本发明下角位置x1对期望轨迹xd的跟踪效果,图3为位置跟踪误差x1-xd,由图2和图3看出,该关节位置能较好地跟踪给定目标轨迹,具有较高的控制精度。
图4和图5反映了本发明方法下连杆角速度状态和电机转角速度状态的估计效果,由图4和图5看出,本发明方法设计的观测器可有效估计实际状态,观测效果良好。
图6反映了本发明方法下该系统控制输入响应曲线图。由图6看出,控制律曲线响应平稳,没有出现较大波动,控制效果平稳。
图7和图8反映了本发明方法下系统应对不同程度不确定性B和
Figure BDA0003847211180000111
时位置的跟踪误差的响应曲线图。由图7和图8看出,本发明方法可以处理系统模型中不同程度的不确定性,而几乎不改变系统的跟踪效果。
针对具有不确定项的柔性关节机械臂系统,本发明方法可以很好的实现对目标轨迹的跟踪控制,跟踪误差小,设计的观测器具有良好的观测效果,所用的控制方法可以很好地处理系统模型中不同程度的不确定性,系统的鲁棒性强。
综上,不论从理论稳定性证明还是从仿真结果看,本发明方法均达到了在模型参数具有不确定性下提高控制效果、降低计算复杂性以及满足对速度项有效估计的预期目标。
当然,以上说明仅仅为本发明的较佳实施例,本发明并不限于列举上述实施例,应当说明的是,任何熟悉本领域的技术人员在本说明书的教导下,所做出的所有等同替代、明显变形形式,均落在本说明书的实质范围之内,理应受到本发明的保护。

Claims (1)

1.基于模糊观测器的柔性关节机械臂指令滤波反步控制方法,其特征在于,
包括如下步骤:
步骤1.建立考虑不确定项的柔性关节机械臂动力学系统数学模型;
首先给出柔性关节机械臂的动态方程,如公式(1)所示:
Figure FDA0003847211170000011
其中,q、
Figure FDA0003847211170000012
分别表示连杆的位置、速度和加速度,qm
Figure FDA0003847211170000013
分别表示电机的转角位置、速度和加速度,M代表连杆的质量,l代表连杆的长度,g为重力加速度,
Figure FDA0003847211170000014
代表摩擦项即未知的干扰项,K表示关节刚度系数,J和B分别表示执行器的惯性和自然阻尼项,u代表电机控制输入,y表示系统的输出;设定
Figure FDA0003847211170000015
和B中包含不确定性;
设状态变量x1=q,
Figure FDA0003847211170000016
x3=qm,
Figure FDA0003847211170000017
将动态方程即公式(1)改写为:
Figure FDA0003847211170000018
为简化动态方程即公式(1),接下来定义如下变量,如公式(3)所示:
Figure FDA0003847211170000019
基于公式(3),进一步将公式(2)改写成:
Figure FDA00038472111700000110
设h2和h4是已知正常数,f2和f4是含有未知项的连续函数,且有连续有界的一阶导数;
步骤2.设计模糊观测器,估计柔性关节机械臂的连杆速度和电机的转角速度;
令公式(4)中
Figure FDA00038472111700000111
进一步得到公式(5):
Figure FDA0003847211170000021
其中,Z1=[x1,x2,x3,x4]T
Figure FDA0003847211170000022
由模糊逻辑系统定义得知:
Figure FDA0003847211170000023
其中,θi是最优参数向量,
Figure FDA0003847211170000024
表示θi的估计,
Figure FDA0003847211170000025
为xi的状态估计值;
Figure FDA0003847211170000026
表示
Figure FDA00038472111700000227
与θi同为函数变量,
Figure FDA0003847211170000027
表示
Figure FDA0003847211170000028
的估计,
Figure FDA0003847211170000029
表示变量为
Figure FDA00038472111700000210
的函数块,
Figure FDA00038472111700000211
表示
Figure FDA00038472111700000212
的估计;
Figure FDA00038472111700000213
Figure FDA00038472111700000214
表示状态变量
Figure FDA00038472111700000215
与θi共同简化成的函数,
Figure FDA00038472111700000216
表示
Figure FDA00038472111700000217
的估计;
针对一个任意给定的正数ε,由模糊自适应逼近定理得到:
Figure FDA00038472111700000218
其中,
Figure FDA00038472111700000219
表示函数中变量为Z1且系数为θi T的函数块,δi(Z1)≤|ε|;
基于公式(5)并综合上式
Figure FDA00038472111700000220
将公式(4)的观测器设计为:
Figure FDA00038472111700000221
式中,
Figure FDA00038472111700000222
Figure FDA00038472111700000223
表示
Figure FDA00038472111700000224
的估计;定义观测器的误差e=[e1,e2,e3,e4]T,其中,ej表示估计误差,j=1,2,3,4;估计误差
Figure FDA00038472111700000225
D1、D2、D3、D4为待设计的正数;得到
Figure FDA00038472111700000226
Figure FDA0003847211170000031
其中,
Figure FDA0003847211170000032
表示
Figure FDA0003847211170000033
的估计,
Figure FDA0003847211170000034
表示
Figure FDA0003847211170000035
的估计;
假定一个函数fi(·)满足以下不等式:
Figure FDA0003847211170000036
Figure FDA0003847211170000037
为正常数,i=2,4;得到
Figure FDA0003847211170000038
其中,
Figure FDA0003847211170000039
表示Z1的估计;
Figure FDA00038472111700000310
其中,
Figure FDA00038472111700000311
是不确定正常数;
Figure FDA00038472111700000312
表示
Figure FDA00038472111700000313
的估计;
Figure FDA00038472111700000314
得到
Figure FDA00038472111700000315
ε表示一个任意小的正数;
设有一对称矩阵Q=QT>0,那么会有一对称矩阵P=PT>0满足ATP+PA=-Q;
步骤3.根据指令滤波和反步法设计基于模糊观测器的柔性关节机械臂的指令滤波反步控制方法,具体过程如下:
定义指令滤波器的公式如下:
Figure FDA00038472111700000316
其中,zo,1表示指令滤波器的一阶变量,zo,2表示指令滤波器的二阶变量;
xo,c=zo,1
Figure FDA00038472111700000317
为指令滤波器的两个输出;αo为指令滤波器的输入,αo表示虚拟控制函数;指令滤波器的初始状态αo(0)=zo,1(0),zo,2(0)=0;
当式(4)系统开始运行后,若输入信号满足:
Figure FDA00038472111700000318
则任意μ>0,存在ζ∈(0,1]和ωn≥0,使得|xo,co|≤μ,o=1,2,3,
Figure FDA00038472111700000319
均有界;
ρ1表示虚拟控制函数一阶导的上界,ρ2表示虚拟控制函数二阶导的上界;
由反步法的原理定义误差变量:
Figure FDA00038472111700000320
其中,zj表示误差变量,j=1,2,3,4;xd为给定的期望信号,xo,c为指令滤波器的输出信号,
Figure FDA0003847211170000041
为速度项的估计值;
定义误差补偿信号为:
Figure FDA0003847211170000042
其中,ξj表示误差补偿信号,vj表示补偿误差跟踪信号;
虚拟控制律和滤波误差补偿信号的具体形式将在以下设计过程中给出;
步骤3.1.选取Lyapunov函数为V0=eTPe,求导得:
Figure FDA0003847211170000043
由杨氏不等式得到:
2eTPF≤||e||2+(m2 2+m4 2)||P||2||e||2 (11)
Figure FDA0003847211170000044
其中,m2、m4表示任意小的正数;将公式(11)、公式(12)代入公式(10)得:
Figure FDA0003847211170000045
其中,λmin(Q)是Q最小的特征值;
步骤3.2.选取第一个子系统的Lyapunov函数
Figure FDA0003847211170000046
求导得:
Figure FDA0003847211170000047
由于e2未知,由杨氏不等式得到:
Figure FDA0003847211170000048
设计虚拟控制律α1和补偿信号
Figure FDA0003847211170000049
为:
Figure FDA00038472111700000410
其中,k1为正常数,将公式(15)、公式(16)代入公式(14)得:
Figure FDA00038472111700000411
步骤3.3.选取Lyapunov函数
Figure FDA0003847211170000051
求导得:
Figure FDA0003847211170000052
其中,
Figure FDA0003847211170000053
为θ2的估计,
Figure FDA0003847211170000054
表示自适应律,λ2为正常数,h2为正常数;
由杨氏不等式得到:
Figure FDA0003847211170000055
将公式(19)代入公式(18)得:
Figure FDA0003847211170000056
设计虚拟控制律α2、补偿信号
Figure FDA0003847211170000057
和自适应律
Figure FDA0003847211170000058
为:
Figure FDA0003847211170000059
其中,k2为正常数,n2为正常数,将公式(21)代入公式(20)得:
Figure FDA00038472111700000510
步骤3.4.选取Lyapunov函数
Figure FDA00038472111700000511
求导得:
Figure FDA00038472111700000512
由于e4未知,由杨氏不等式得到:
Figure FDA00038472111700000513
设计虚拟控制律α3和补偿信号
Figure FDA0003847211170000061
为:
Figure FDA0003847211170000062
其中,k3为正常数,将公式(24)、公式(25)代入公式(23)得:
Figure FDA0003847211170000063
步骤3.5.选取Lyapunov函数
Figure FDA0003847211170000064
求导得:
Figure FDA0003847211170000065
其中,λ4为正常数;
Figure FDA0003847211170000066
Figure FDA0003847211170000067
为θ4的估计,
Figure FDA0003847211170000068
表示自适应律;
由杨氏不等式得到:
Figure FDA0003847211170000069
将公式(28)代入公式(27)得:
Figure FDA00038472111700000610
设计真实控制律u、补偿信号
Figure FDA00038472111700000611
和自适应律
Figure FDA00038472111700000612
为:
Figure FDA00038472111700000613
其中,k4、n4均为正的常数,将公式(30)代入公式(29)得:
Figure FDA00038472111700000614
步骤4.对基于模糊观测器的柔性关节机械臂指令滤波反步控制方法进行稳定性分析;
选取V=V4,由公式(31)得:
Figure FDA0003847211170000071
由杨氏不等式得:
Figure FDA0003847211170000072
Figure FDA0003847211170000073
将公式(13)、公式(33)、公式(34)代入(32)得到:
Figure FDA0003847211170000074
其中,
Figure FDA0003847211170000075
Figure FDA0003847211170000076
Figure FDA0003847211170000077
其中,λmax(P)表示向量P的最大特征值;进一步由公式(35)得到:
Figure FDA0003847211170000078
式(36)表明,vj
Figure FDA0003847211170000079
均属于紧集:
Figure FDA00038472111700000710
因此,机械臂系统的所有信号都有界,滤波器的误差补偿信号ξj满足
Figure FDA00038472111700000711
其中
Figure FDA00038472111700000712
因为z1=v11且ξ1有界,跟踪误差z1有界,进而得到:
Figure FDA00038472111700000713
其中,minj(kj)表示kj中的最小值;
通过以上分析得知,vj、ξj、zj都是有界的,因此,在给定参数kj、mi、ni后,保持a不变,选择足够大的λi和足够小的ε,能够保证跟踪误差z1足够小。
CN202211123023.2A 2022-09-15 2022-09-15 基于模糊观测器的柔性关节机械臂指令滤波反步控制方法 Pending CN115473467A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211123023.2A CN115473467A (zh) 2022-09-15 2022-09-15 基于模糊观测器的柔性关节机械臂指令滤波反步控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211123023.2A CN115473467A (zh) 2022-09-15 2022-09-15 基于模糊观测器的柔性关节机械臂指令滤波反步控制方法

Publications (1)

Publication Number Publication Date
CN115473467A true CN115473467A (zh) 2022-12-13

Family

ID=84332597

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211123023.2A Pending CN115473467A (zh) 2022-09-15 2022-09-15 基于模糊观测器的柔性关节机械臂指令滤波反步控制方法

Country Status (1)

Country Link
CN (1) CN115473467A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117885103A (zh) * 2024-03-14 2024-04-16 山东大学 基于降阶扩张状态观测器的柔性机械臂控制方法及系统

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117885103A (zh) * 2024-03-14 2024-04-16 山东大学 基于降阶扩张状态观测器的柔性机械臂控制方法及系统
CN117885103B (zh) * 2024-03-14 2024-05-17 山东大学 基于降阶扩张状态观测器的柔性机械臂控制方法及系统

Similar Documents

Publication Publication Date Title
CN108942924B (zh) 基于多层神经网络的模型不确定性机械臂运动控制方法
CN112817231B (zh) 一种具有强鲁棒性的机械臂高精度跟踪控制方法
CN108388114B (zh) 一种基于输出重定义的柔性机械臂复合控制方法
CN112241124A (zh) 一种自适应反演积分非奇异快速终端滑模控制器设计方法
CN111965976B (zh) 基于神经网络观测器的机器人关节滑模控制方法及系统
CN107193211B (zh) 基于自抗扰和反演技术的单臂机械手控制器及其设计方法
CN112947067A (zh) 一种水下机器人三维轨迹精确跟踪控制方法
CN109240092B (zh) 基于多智能体可重构模块化柔性机械臂轨迹跟踪控制方法
CN112549030B (zh) 一种基于抗干扰滤波的空间机械臂柔性关节位姿估计方法
CN114169230A (zh) 一种机器人动力学参数的辨识方法
CN115157238A (zh) 一种多自由度机器人动力学建模和轨迹跟踪方法
CN115284284A (zh) 基于状态受限的柔性机械臂奇异摄动控制方法
CN115981162A (zh) 一种基于新型扰动观测器的机器人系统滑模控制轨迹跟踪方法
CN112269317A (zh) 一种基于扩展卡尔曼滤波器的双边遥操作控制方法
CN115473467A (zh) 基于模糊观测器的柔性关节机械臂指令滤波反步控制方法
CN107765548B (zh) 基于双观测器的发射平台高精度运动控制方法
CN115816453A (zh) 基于tde的自适应超螺旋多变量快速终端滑模控制方法
CN114942593A (zh) 一种基于干扰观测器补偿的机械臂自适应滑模控制方法
CN113377006B (zh) 一种基于不变流形观测器的全局快速终端滑模控制方法
CN112936277B (zh) 一种水下机器人-机械手系统固定时间轨迹跟踪方法
CN114578697A (zh) 一种电机驱动机械手的多约束自适应控制方法
CN112223276B (zh) 基于自适应神经网络滑模控制的多关节机器人控制方法
CN116339141B (zh) 一种机械臂全局固定时间轨迹跟踪滑模控制方法
CN109048995B (zh) 一种三自由度Delta并联机器人的非线性关节摩擦力补偿方法
CN114003002B (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