CN108710296A - 微陀螺仪的分数阶自适应快速终端滑模控制方法 - Google Patents
微陀螺仪的分数阶自适应快速终端滑模控制方法 Download PDFInfo
- Publication number
- CN108710296A CN108710296A CN201810384011.2A CN201810384011A CN108710296A CN 108710296 A CN108710296 A CN 108710296A CN 201810384011 A CN201810384011 A CN 201810384011A CN 108710296 A CN108710296 A CN 108710296A
- Authority
- CN
- China
- Prior art keywords
- gyroscope
- fractional order
- parameter
- adaptive
- terminal sliding
- 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
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive 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/042—Adaptive 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)
- Gyroscopes (AREA)
Abstract
本发明公开了一种微陀螺仪的分数阶自适应快速终端滑模控制方法,包括:建立微陀螺仪***的无量纲动力学方程数学模型和参考轨迹模型;构建基于分数阶的自适应快速终端滑模控制器。本发明能够实现微陀螺仪实时跟踪的目标,使误差在有限时间内快速收敛到零,增强了***的鲁棒性,在有外界干扰的情况下仍然保持很好的性能;根据分数阶快速终端滑模面设计分数阶自适应律,基于Lyapunov稳定性判据设计一种自适应辨识方法,在线实时估计微陀螺仪的各项未知***参数,与整数阶相比,多了可调项,提高了控制效果和参数估计效果;采用自适应控制方法对***不确定性和外界干扰总和的上界的参数进行估计,达到对制造误差和外界扰动的自动跟踪。
Description
技术领域
本发明涉及微陀螺仪控制技术领域,具体是一种微陀螺仪的分数阶自适应快速终端滑模控制方法。
背景技术
微陀螺仪是测量惯性导航和惯性制导***角速度的传感器,因为它在任何环境下都能够自主导航,所以自出现以来就受到人们的广泛关注,在航天、航海、航空和军事等领域得到了普遍的应用。但是生产和制造过程中存在误差且易受温度影响,造成元件特性与设计之间的差异,从而导致微陀螺仪的性能降低。另外,微陀螺仪属于多输入多输出***并且***参数存在不确定性及易受外界环境的影响,使得陀螺仪追踪效果并不理想。传统的滑模控制方法中滑模面的设计都是采用误差的比例、微分或积分的线性组合,其中微分或积分的阶数都是整数,微陀螺仪的跟踪效果比较差,***参数和角速度估计效果也比较差,容易引起抖振。
发明内容
本发明为了避免上述不足之处,提出一种微陀螺仪的分数阶自适应快速终端滑模控制方法,在线实时估计陀螺仪参数矩阵,有效提高控制效果和参数估计效果,保证***的跟踪误差在有限时间内收敛到零,并且有效降低抖振。
本发明采用的技术方案是:
一种微陀螺仪的分数阶自适应快速终端滑模控制方法,包括如下步骤:
一、建立微陀螺仪***的无量纲动力学方程数学模型;
二、构建微陀螺仪***的参考轨迹模型;
三、构建基于分数阶的自适应快速终端滑模控制器,首先,基于Lyapunov稳定性理论设计分数阶快速终端滑模控制器,确保***状态可以在很短的有限时间内收敛到平衡点;然后,设计分数阶自适应快速终端滑模控制律,在线实时估计微陀螺仪的***参数;最后,为了避免控制***的抖振现象,采用自适应方法对***不确定性和外界干扰总和上界的参数进行估计,达到对制造误差和外界扰动的自动跟踪。
本发明的微陀螺仪的分数阶自适应快速终端滑模控制方法中,首先采用分数阶快速终端滑模面,根据Lyapunov函数设计带有分数阶的快速终端滑模控制律,快速终端滑模保证***的跟踪误差迅速在有限时间内收敛到零。再通过调节分数阶阶数,来提高***的跟踪性能。然后,采用自适应控制的方法对***不确定性和外界干扰总和的上界的参数进行估计,设计分数阶自适应快速终端滑模控制器,达到对制造误差和外界干扰的自动跟踪,从而减小***抖振。
具体地,Lyapunov(李雅普诺夫)函数V(x,t)是可以根据需要自由选取的,但是要符合李雅普诺夫函数的要求(即李雅普诺夫稳定性判据),即V(x,t)是正定的,当x≠0时, V(x,t)>0;当x=0时,V(0,t)=0,具有连续的偏导数;V(x,t)的一阶导是负半定的。
优选地,所述步骤一中:所述微陀螺仪的无量纲动力学方程数学模型为:
其中,q为微陀螺仪数学模型无量纲化后的位置向量,表示对位置向量q求一阶导数,表示对位置向量q求二阶导数;Db为无量纲化后阻尼矩阵,K为无量纲化后弹簧系数矩阵,Ω为无量纲化后角速度矩阵,u为微陀螺仪的控制输入,f为***不确定性和外界干扰总和,且有:
其中,ΔDb为惯性矩阵Db+2Ω的未知参数的不确定性,ΔK为惯性矩阵K的未知参数的不确定性,d(t)是外界干扰;
假设***不确定性和外界干扰总和满足:
其中,||f||表示求向量f的模,b0,b1,b2是***的不确定性和外界干扰的上界的参数,为正常数。
优选地,所述步骤二中:所述微陀螺仪***的参考轨迹模型为::
qr1=A1sin(ω1t),qr2=A2sin(ω2t) (4)
其中qr1、qr2是指x轴、y轴两轴的参考轨迹模型;A1、A2分别是微陀螺仪在x轴、y轴方向上的振幅;ω1、ω2分别是微陀螺仪在x轴、y轴方向上给定的振动频率;t是时间变量。
优选地,所述步骤三中,基于Lyapunov稳定性理论,设计分数阶快速终端滑模控制器,设计步骤包括:
定义微陀螺仪实际位置和参考轨迹之间的跟踪误差e为:
e=q-qr (5)
其中,qr为微陀螺仪***两轴的参考轨迹;
针对微陀螺仪设计分数阶终快速端滑模面s为:
其中,λ1,λ2,γ为正实数,0<γ<1,Dα-1e是指对e求α-1阶导,这是分数阶专用的写法,D 表示求导,α是指分数阶的阶数,|e|表示对e求绝对值,sign为符号函数,则微陀螺仪分数阶快速终端滑模面s的导数为:
其中,Dαe是对Dα-1e求导的结果;
不考虑***参数不确定性和外界干扰,令得到等效控制律ueq为:
不连续的控制项Δu为:
则基于分数阶的快速终端滑模控制律设计为:
设计Lyapunov函数V1为:
对上式两边求导并代入u',得到说明微陀螺仪***是稳定的。
优选地,所述步骤三中,基于分数阶的自适应快速终端滑模控制律设计步骤包括:
定义微陀螺仪无量纲化模型中估计误差以及***不确定性和外界干扰总和上界的参数估计误差分别为:
其中,分别是微陀螺仪参数矩阵Db,K,Ω的估计值,分别为参数矩阵Db,K,Ω的参数估计误差;是***不确定性和外界干扰总和上界的参数估计值,为相应的估计误差;
基于分数阶的自适应快速终端滑模控制律u为:
其中:
定义Lyapunov函数V为:
其中,M=MT>0,N=NT>0,P=PT>0为对称正定矩阵,tr{·}表示矩阵的求迹运算,xi为正常数;
为了保证Lyapunov函数的导数在线实时估计微陀螺仪的***参数以及***不确定性和外界干扰总和上界的参数,设计微陀螺仪参数矩阵Db,K,Ω的估计值和***不确定和外界干扰总和上界参数估计值的自适应律分别为:
为陀螺仪***阻尼矩阵Db的估计值的自适应律的转置,为陀螺仪***的弹簧系数矩阵 K的估计值的自适应律的转置,为陀螺仪角速度矩阵的估计值的自适应律的转置;当选取上述参数自适应律时,满足李雅普诺夫稳定性定理,保证微陀螺仪***全局稳定性。
与现有技术相比,本发明的有益效果在于:1)本发明设计快速终端滑模面,确保了***跟踪误差迅速在有限时间内收敛到零;在此基础上引入分数阶,即带有偏差的分数阶导数的快速终端滑模面,则多了可以调节的微分阶数项,增加了调节自由度,使得控制效果有所改善;当微陀螺仪***参数和角速度看作未知变量时,基于Lyapunov稳定性理论,设计参数的自适应律,实时更新微陀螺仪的角速度和其他***参数的估计值;同时,考虑到制造误差和外界扰动的上界未知的情况,采用自适应方法对***不确定性和外界干扰总和上界的参数进行估计,达到对制造误差和外界干扰的自动跟踪;2)本发明具有提高陀螺仪***控制性能和鲁棒性的特点,能够实现微陀螺仪较好的跟踪效果,使跟踪误差在有限时间内收敛到零并且降低***抖振,以及辨识出***参数和角速度,辨识***参数可以了解***的物理参数,如果某个参数距离给定的值很远,说明***内部出现问题,需要更换原件等。
附图说明
图1为本发明分数阶自适应快速终端滑模控制***的原理图;
图2为本发明采用自适应快速终端滑模控制方法的微陀螺仪X,Y轴跟踪效果图;
图3为本发明采用自适应快速终端滑模控制方法的跟踪误差图;
图4为本发明采用自适应快速终端滑模控制方法的两轴控制输入响应曲线图。
具体实施方式
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制本发明的保护范围。
如图1所示,一种微陀螺仪的分数阶自适应快速终端滑模控制方法,包括以下步骤:
一、建立微陀螺仪***的无量纲动力学方程数学模型
一般的微机械振动陀螺仪由三个部分组成:弹性材料所支撑悬挂的质量块、静电驱动装置和感测装置,将其简化为一个由质量块和弹簧构成的有阻尼振荡***。对z轴微陀螺仪而言,可以认为质量块被限制只能在x-y平面内运动,而不能沿z轴运动,只围绕z轴旋转。
根据旋转系中的牛顿定律,考虑加工制造误差,再对数学模型进行无量纲化处理及等效变换后,得到的微陀螺仪的数学模型如下:
其中:q为微陀螺仪数学模型无量纲化后的位置向量,q1,q2分别为微陀螺仪数学模型无量纲化后x轴、y轴的位置向量;为参考位移,Q1,Q2分别为微陀螺仪x轴、y轴的位置;表示对位置向量q求一阶导数,表示对位置向量q求二阶导数;为微陀螺仪的控制输入,ux,uy分别表示x轴和y轴的控制输入;为无量钢化后阻尼矩阵,其中,箭头→表示后者由前者计算得到是等于关系,Dxx,Dyy分别为微陀螺仪x轴、y轴的阻尼系数,Dxy为两轴的耦合阻尼系数,m为陀螺仪的质量,ω0为两轴的共振频率;为无量纲化后弹簧系数矩阵,其中,kxx,kyy分别为微陀螺仪x轴、y轴的弹簧系数,kxy为两轴耦合的弹簧系数;为无量纲化后角速度矩阵,其中Ω* z为z轴方向上的输入角速度,微陀螺仪是沿着x-y平面运动,但是同时它也在围绕z轴旋转,所以有z轴角速度;ΔDb为惯性矩阵Db+2Ω的未知参数的不确定性,ΔK为惯性矩阵K的未知参数的不确定性,d(t)是外界干扰。
从而式(1)可以改写为:
其中,f为***不确定性和外界干扰总和,且有:
假设***不确定性和外界干扰总和满足:
b0,b1,b2是***的不确定性和外界干扰的上界的参数,为正常数,||f||表示求向量的模值。
二、构建微陀螺仪***的参考轨迹模型
微陀螺仪的理想动态特性是一种无能量损耗,x轴、y轴两轴间无动态耦合的稳定正弦振荡,参考轨迹模型描述如下:
qr1=A1sin(ω1t),qr2=A2sin(ω2t) (5)
其中qr1、qr2是指x轴、y轴两轴的参考轨迹模型;A1、A2分别是微陀螺仪在x轴、y轴方向上的振幅;ω1、ω2分别是微陀螺仪在x轴、y轴方向上给定的振动频率;t是时间变量。
将(5)式转化成微分方程形式为:
其中,qr=[qr1 qr2]T为微陀螺仪***两轴的参考轨迹,表示对x轴、y轴的参考轨迹qr求二阶导,Kr=diag(ω1 2,ω2 2),diag()表示对角矩阵,矩阵Kr表示微陀螺仪在x轴、y轴方向上给定的振动频率的对角矩阵。
三、构建基于分数阶的自适应快速终端滑模控制器
基于分数阶的自适应快速终端滑模控制器设计步骤如下:
3.1)基于Lyapunov稳定性方法设计分数阶快速终端滑模控制器,确保***状态迅速在有限时间内收敛至平衡点。
定义微陀螺仪实际位置和参考轨迹之间的跟踪误差e为:
e=q-qr (7)
其中,qr为微陀螺仪***x、y轴的参考轨迹。
针对微陀螺仪设计分数阶终快速端滑模面s为:
其中,λ1,λ2,γ为正实数,0<γ<1,Dα-1e是指对e求α-1阶导,这是分数阶专用的写法,D 表示求导,α是指分数阶的阶数,|e|表示对e求绝对值,sign为符号函数,则微陀螺仪分数阶快速终端滑模面s的导数为:
其中,Dαe是对Dα-1e求导的结果。
不考虑***参数不确定性和外界干扰,令得到等效控制律ueq为:
不连续的控制项Δu为:
则基于分数阶的快速终端滑模控制律u'设计为:
稳定性证明:
设计Lyapunov函数V1为:
对上式两边求导并代入公式(4)、(9)、(12)可以得到:
根据Lyapunov稳定性判据可以证明微陀螺仪***是稳定的;
3.2)设计基于分数阶的自适应快速终端滑模控制器的控制律u:
由于微陀螺仪无量纲化模型中Db,K,Ω三个参数是未知的,并且在实际应用过程中,***存在不确定性和外界干扰,因此基于分数阶的快速终端滑模控制律u'公式(12)不能直接应用。基于自适应理论知识,用估计值代替(12)中的三个陀螺仪参数矩阵,设计三个参数估计值的自适应算法,在线实时更新估计;同时为了避免抖振问题,采用自适应控制的方法对***不确定性和外界干扰总和的上界的参数进行估计,用估计值代替(12) 中***不确定性和外界干扰总和上界的参数b0,b1,b2,设计分数阶自适应快速终端滑模控制器,达到对制造误差和外界扰动的自动跟踪。
定义微陀螺仪无量纲化模型中参数估计误差以及***不确定性和外界干扰总和上界的参数估计误差分别为:
基于分数阶的快速终端滑模控制律公式(12)可以转换为基于分数阶的自适应快速终端滑模控制律u:
其中:
稳定性证明:
将控制律式(16)代入公式(9)可得:
定义李雅普诺夫函数V为:
其中,M=MT>0,N=NT>0,P=PT>0为对称正定矩阵,称作自适应固定增益矩阵,它们是可调整的参数,是可人为选取的;xi为正常数;tr{·}表示矩阵的求迹运算。
对Lyapunov函数等式(19)两边求导并代入公式(18)得:
由于Db=Db T,K=KT,Ω=-ΩT,并且为标量,因此:
同理可得:
将公式(21)、(22)代入公式(20)得:
为保证设计参数自适应律分别为:
为陀螺仪***阻尼矩阵Db的估计值的自适应律的转置,为陀螺仪***的弹簧系数矩阵K的估计值的自适应律的转置,为陀螺仪角速度矩阵的估计值的自适应律的装置,是***不确定性和外界干扰总和上界参数的自适应律。
那么V对时间t的微分变为:
由于满足李雅普诺夫稳定性定理,说明跟踪轨迹在有限时间内到达了滑模面,并停留在滑模面上,从而验证本发明所提出的控制方法能够实现对微陀螺仪的分数阶自适应快速终端滑模的控制,具有鲁棒作用,提高了***控制性能,使得跟踪效果更好且可以在有限时间内收敛到零,并且降低抖振,同时在线辨识微陀螺仪***参数,从而了解微陀螺仪***的物理参数,如果参数辨识中某个估计值距离给定的值Db,K,Ω很远,说明***内部出现问题,需要更换原件等。
为了更加直观地显示本发明提出的分数阶自适应快速终端滑模控制方法的优越性,在 MATLAB/SIMULINK对本发明控制方法进行计算机仿真实验。
参考现有文献,选取微陀螺仪的参数为:
m=1.8×10-7kg,Dxx=1.8×10-6N s/m,Dyy=1.8×10-6Ns/m,
Dxy=3.6×10-7Ns/m,kxx=63.955N/m,kyy=95.92N/m,kxy=12.779N/m。
假定未知的输入角速度为Ωz *=100rad/s。首先微陀螺仪参数进行无量纲化处理,对于振动式微陀螺仪,基础质量块沿x轴,y轴振动幅值在亚微米级别,参考位移选取为q0=1μm较为合理。因为微陀螺仪的操作频率在kHz范围,所以选取固有频率为ω0=1kHz。得到微陀螺仪的无量纲化参数为:
ωx 2=355.3,ωy 2=532.9,ωxy=70.99,dxx=0.01
dyy=0.01,dxy=0.002,Ωz=0.1
仿真实验中,设***的初始条件为:q1(0)=0.5,q2(0)=0.5,q1(0),q2(0)是指微陀螺仪在x轴、y轴的初始位置,此刻t=0,微陀螺仪的两轴期望运行轨迹(参考轨迹) 为:qr1=sin(4.17t),qr2=1.2sin(5.11t),三个参数矩阵的估计初值分别为:
在滑模控制律中,滑模面参数取λ1=10,λ2=0.5,γ=0.5,取M=N=P=diag(150,150),x0=1.0,x1=1.6,x2=1.6,取分数阶阶数α=0.9。
当外界干扰取白噪声信号d=[0.5*randn(1,1);0.5*randn(1,1)]时,randn(1,1)表示产生随机数,仿真结果如图2至图4所示。
图2为微陀螺仪X、Y轴的跟踪曲线,从图中可以看出采用分数阶自适应快速终端滑模控制方法后微陀螺仪的X、Y轴轨迹能够很好追踪上参考轨迹。
图3为跟踪误差曲线图,从图中可以看出,X、Y轴上的跟踪误差e能在有限时间内收敛到零。
图4为两轴控制输入响应曲线。自适应方法对控制器中***不确定性和外界干扰总和的上界的参数进行估计,有效降低了抖振。
本发明应用于微陀螺仪的基于分数阶自适应快速终端滑模控制方法,该方法对于微陀螺仪进行有效、可靠的控制,在对***参数未知的情况下,能够有效估计出***的各项参数,并且保证***的全局稳定性。在滑模控制的设计中,选取分数阶快速终端滑模面,在保留滑模控制优势的基础上,又多了可以调节的微分阶数项,增加了调节自由度,改善了控制性能;快速终端滑模面保证了***误差迅速在有限时间内收敛到零,提高了跟踪性能;基于自适应机制,设计自适应律,实时估计陀螺仪***的参数和角速度;另外,采用自适应方法对***不确定性和外界干扰总和的上界的参数进行估计,达到对制造误差和外界扰动的自动跟踪,有效降低了抖振。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上大的限制,虽然本发明以较佳实例揭露如上,然而并非用以限制本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实例,但凡是为脱离本发明的技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与装饰,均仍属于本发明技术方案的范围内。
Claims (5)
1.一种微陀螺仪的分数阶自适应快速终端滑模控制方法,其特征在于:包括如下步骤:
一、建立微陀螺仪***的无量纲动力学方程数学模型;
二、构建微陀螺仪***的参考轨迹模型;
三、构建基于分数阶的自适应快速终端滑模控制器,首先,基于Lyapunov稳定性理论设计分数阶快速终端滑模控制器,确保***状态可以在很短的有限时间内收敛到平衡点;然后,设计分数阶自适应快速终端滑模控制律,在线实时估计微陀螺仪的***参数;最后,为了避免控制***的抖振现象,采用自适应方法对***不确定性和外界干扰总和上界的参数进行估计,达到对制造误差和外界扰动的自动跟踪。
2.根据权利要求1中所述的一种微陀螺仪的分数阶自适应快速终端滑模控制方法,其特征在于:所述步骤一中:所述微陀螺仪的无量纲动力学方程数学模型为:
其中,q为微陀螺仪数学模型无量纲化后的位置向量,表示对位置向量q求一阶导数,表示对位置向量q求二阶导数;Db为无量纲化后阻尼矩阵,K为无量纲化后弹簧系数矩阵,Ω为无量纲化后角速度矩阵,u为微陀螺仪的控制输入,f为***不确定性和外界干扰总和,且有:
其中,ΔDb为惯性矩阵Db+2Ω的未知参数的不确定性,ΔK为惯性矩阵K的未知参数的不确定性,d(t)是外界干扰;假设***不确定性和外界干扰总和满足:
其中,||f||表示求向量f的模,b0,b1,b2是***的不确定性和外界干扰的上界的参数,为正常数。
3.根据权利要求1中所述的一种微陀螺仪的分数阶自适应快速终端滑模控制方法,其特征在于:所述步骤二中:所述微陀螺仪***的参考轨迹模型为:
qr1=A1sin(ω1t),qr2=A2sin(ω2t) (4)
其中qr1、qr2是指x轴、y轴两轴的参考轨迹模型;A1、A2分别是微陀螺仪在x轴、y轴方向上的振幅;ω1、ω2分别是微陀螺仪在x轴、y轴方向上给定的振动频率;t是时间变量。
4.根据权利要求1中所述的一种微陀螺仪的分数阶自适应快速终端滑模控制方法,其特征在于:所述步骤三中,基于Lyapunov稳定性理论,设计分数阶快速终端滑模控制器,设计步骤包括:
定义微陀螺仪实际位置和参考轨迹之间的跟踪误差e为:
e=q-qr (5)
其中,qr为微陀螺仪***两轴的参考轨迹;
针对微陀螺仪设计分数阶终快速端滑模面s为:
其中,λ1,λ2,γ为正实数,0<γ<1,Dα-1e是指对e求α-1阶导,这是分数阶专用的写法,D表示求导,α是指分数阶的阶数,|e|表示对e求绝对值,sign为符号函数,则微陀螺仪分数阶快速终端滑模面s的导数为:
其中,Dαe是对Dα-1e求导的结果;
不考虑***参数不确定性和外界干扰,令得到等效控制律ueq为:
不连续的控制项Δu为:
则基于分数阶的快速终端滑模控制律设计为:
设计Lyapunov函数V1为:
对上式两边求导并代入u',得到说明微陀螺仪***是稳定的。
5.根据权利要求4中所述的一种微陀螺仪的分数阶自适应快速终端滑模控制方法,其特征在于:所述步骤三中,基于分数阶的自适应快速终端滑模控制律u设计步骤包括:
定义微陀螺仪无量纲化模型中估计误差以及***不确定性和外界干扰总和上界的参数估计误差分别为:
其中,分别是微陀螺仪参数矩阵Db,K,Ω的估计值,分别为参数矩阵Db,K,Ω的参数估计误差;是***不确定性和外界干扰总和上界的参数估计值,为相应的估计误差;
基于分数阶的自适应快速终端滑模控制律u为:
其中:
定义Lyapunov函数V为:
其中,M=MT>0,N=NT>0,P=PT>0为对称正定矩阵,tr{·}表示矩阵的求迹运算,xi为正常数;
为了保证Lyapunov函数的导数在线实时估计微陀螺仪的***参数以及***不确定性和外界干扰总和上界的参数,设计微陀螺仪参数矩阵Db,K,Ω的估计值和***不确定和外界干扰总和上界参数估计值的自适应律分别为:
为陀螺仪***阻尼矩阵Db的估计值的自适应律的转置,为陀螺仪***的弹簧系数矩阵K的估计值的自适应律的转置,为陀螺仪角速度矩阵的估计值的自适应律的转置,是***不确定性和外界干扰总和上界参数的自适应律;当选取上述参数自适应律时,满足李雅普诺夫稳定性定理,保证微陀螺仪***全局稳定性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810384011.2A CN108710296B (zh) | 2018-04-26 | 2018-04-26 | 微陀螺仪的分数阶自适应快速终端滑模控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810384011.2A CN108710296B (zh) | 2018-04-26 | 2018-04-26 | 微陀螺仪的分数阶自适应快速终端滑模控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108710296A true CN108710296A (zh) | 2018-10-26 |
CN108710296B CN108710296B (zh) | 2021-06-01 |
Family
ID=63867431
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810384011.2A Active CN108710296B (zh) | 2018-04-26 | 2018-04-26 | 微陀螺仪的分数阶自适应快速终端滑模控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108710296B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109557816A (zh) * | 2018-12-28 | 2019-04-02 | 武汉工程大学 | 一种压电陶瓷执行器迟滞特性的抑制方法、***及介质 |
CN110058525A (zh) * | 2019-05-17 | 2019-07-26 | 金陵科技学院 | 分数阶离心飞轮调速***有限时间镇定控制方法 |
CN110118609A (zh) * | 2019-05-16 | 2019-08-13 | 中国计量大学 | 基于自适应拟随机数搜索的热电偶分数阶***辨识方法 |
CN110262242A (zh) * | 2019-07-01 | 2019-09-20 | 河海大学常州校区 | 一种微陀螺仪自适应分数阶滑模控制方法 |
CN110647036A (zh) * | 2019-09-23 | 2020-01-03 | 河海大学常州校区 | 微陀螺仪自适应双反馈模糊神经网络分数阶滑模控制方法 |
CN111290015A (zh) * | 2020-03-09 | 2020-06-16 | 贵州大学 | 具有约束的分数阶自持式机电地震仪***加速稳定控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013074254A1 (en) * | 2011-11-14 | 2013-05-23 | Qualcomm Mems Technologies, Inc. | Combined resonators and passive circuit components for filter passband flattening |
CN103728882A (zh) * | 2014-01-07 | 2014-04-16 | 河海大学常州校区 | 微陀螺仪的自适应反演非奇异终端滑模控制方法 |
CN105204345A (zh) * | 2015-10-22 | 2015-12-30 | 河海大学常州校区 | 一种自适应分数阶滑模控制方法 |
CN107390517A (zh) * | 2017-07-21 | 2017-11-24 | 北京交通大学 | 用于列车ato***的鲁棒自适应非奇异终端滑模控制方法 |
CN107831655A (zh) * | 2017-10-23 | 2018-03-23 | 河海大学常州校区 | 微陀螺仪的分数阶自适应反演模糊滑模控制方法 |
-
2018
- 2018-04-26 CN CN201810384011.2A patent/CN108710296B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013074254A1 (en) * | 2011-11-14 | 2013-05-23 | Qualcomm Mems Technologies, Inc. | Combined resonators and passive circuit components for filter passband flattening |
CN103728882A (zh) * | 2014-01-07 | 2014-04-16 | 河海大学常州校区 | 微陀螺仪的自适应反演非奇异终端滑模控制方法 |
CN105204345A (zh) * | 2015-10-22 | 2015-12-30 | 河海大学常州校区 | 一种自适应分数阶滑模控制方法 |
CN107390517A (zh) * | 2017-07-21 | 2017-11-24 | 北京交通大学 | 用于列车ato***的鲁棒自适应非奇异终端滑模控制方法 |
CN107831655A (zh) * | 2017-10-23 | 2018-03-23 | 河海大学常州校区 | 微陀螺仪的分数阶自适应反演模糊滑模控制方法 |
Non-Patent Citations (1)
Title |
---|
黄滨: "高精度随动***的变结构控制", 《船电技术》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109557816A (zh) * | 2018-12-28 | 2019-04-02 | 武汉工程大学 | 一种压电陶瓷执行器迟滞特性的抑制方法、***及介质 |
CN109557816B (zh) * | 2018-12-28 | 2021-06-29 | 武汉工程大学 | 一种压电陶瓷执行器迟滞特性的抑制方法、***及介质 |
CN110118609A (zh) * | 2019-05-16 | 2019-08-13 | 中国计量大学 | 基于自适应拟随机数搜索的热电偶分数阶***辨识方法 |
CN110058525A (zh) * | 2019-05-17 | 2019-07-26 | 金陵科技学院 | 分数阶离心飞轮调速***有限时间镇定控制方法 |
CN110058525B (zh) * | 2019-05-17 | 2021-10-01 | 金陵科技学院 | 分数阶离心飞轮调速***有限时间镇定控制方法 |
CN110262242A (zh) * | 2019-07-01 | 2019-09-20 | 河海大学常州校区 | 一种微陀螺仪自适应分数阶滑模控制方法 |
CN110262242B (zh) * | 2019-07-01 | 2022-03-15 | 河海大学常州校区 | 一种微陀螺仪自适应分数阶滑模控制方法 |
CN110647036A (zh) * | 2019-09-23 | 2020-01-03 | 河海大学常州校区 | 微陀螺仪自适应双反馈模糊神经网络分数阶滑模控制方法 |
CN110647036B (zh) * | 2019-09-23 | 2022-04-22 | 河海大学常州校区 | 微陀螺仪自适应双反馈模糊神经网络分数阶滑模控制方法 |
CN111290015A (zh) * | 2020-03-09 | 2020-06-16 | 贵州大学 | 具有约束的分数阶自持式机电地震仪***加速稳定控制方法 |
CN111290015B (zh) * | 2020-03-09 | 2022-05-13 | 贵州大学 | 具有约束的分数阶自持式机电地震仪***加速稳定控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108710296B (zh) | 2021-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108710296A (zh) | 微陀螺仪的分数阶自适应快速终端滑模控制方法 | |
CN108227504A (zh) | 微陀螺分数阶自适应模糊神经反演终端滑模控制方法 | |
CN107831655B (zh) | 微陀螺仪的分数阶自适应反演模糊滑模控制方法 | |
CN103616818B (zh) | 微陀螺仪的自适应模糊神经全局快速终端滑模控制方法 | |
CN104281056B (zh) | 基于神经网络上界学习的微陀螺仪鲁棒自适应控制方法 | |
CN105045097B (zh) | 一种基于神经网络的微陀螺仪反演全局滑模模糊控制方法 | |
CN103728882B (zh) | 微陀螺仪的自适应反演非奇异终端滑模控制方法 | |
CN102508434B (zh) | 用于微陀螺仪的自适应模糊滑模控制器 | |
CN103116275B (zh) | 基于滑模补偿的微陀螺仪鲁棒神经网络控制***及方法 | |
CN103324087B (zh) | 基于神经网络的微陀螺仪的自适应反演控制***及方法 | |
CN102298322B (zh) | 基于模型参考自适应控制微陀螺仪的方法 | |
CN107831660A (zh) | 微陀螺仪自适应高阶超扭曲滑模控制方法 | |
CN105929694A (zh) | 一种微陀螺自适应神经网络非奇异终端滑模控制方法 | |
CN103279038B (zh) | 基于t-s模糊模型的微陀螺仪滑模自适应控制方法 | |
CN104503246B (zh) | 微陀螺仪***的间接自适应神经网络滑模控制方法 | |
CN104267604B (zh) | 微陀螺仪自适应神经网络全局滑模控制方法 | |
CN105278331A (zh) | 一种微陀螺的鲁棒自适应神经网络h无穷控制方法 | |
CN107807527A (zh) | 微陀螺仪增益可调的自适应超扭曲滑模控制方法 | |
CN109240083B (zh) | 微陀螺仪***的自适应模糊超扭曲滑模控制方法 | |
CN104155874B (zh) | 微陀螺仪的反演自适应模糊动态滑模控制方法 | |
CN107703757A (zh) | 微陀螺仪***的超扭曲滑模控制方法 | |
CN106338918A (zh) | 一种微陀螺自适应动态面双神经网络控制方法 | |
CN103345155B (zh) | 微陀螺仪的自适应反演控制***及方法 | |
CN103529701A (zh) | 微陀螺仪的神经网络全局滑模控制方法 | |
CN110389528B (zh) | 基于扰动观测的数据驱动mems陀螺仪驱动控制方法 |
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 |