CN105068426B - 基于干扰补偿的电液位置伺服***连续滑模控制方法 - Google Patents
基于干扰补偿的电液位置伺服***连续滑模控制方法 Download PDFInfo
- Publication number
- CN105068426B CN105068426B CN201510522508.2A CN201510522508A CN105068426B CN 105068426 B CN105068426 B CN 105068426B CN 201510522508 A CN201510522508 A CN 201510522508A CN 105068426 B CN105068426 B CN 105068426B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mover
- centerdot
- mtd
- 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
Links
Landscapes
- Feedback Control In General (AREA)
- Steroid Compounds (AREA)
Abstract
本发明公开了一种基于干扰补偿的电液位置伺服***连续滑模控制方法,方法包括:建立电液位置伺服***的数学模型;设计电液位置伺服***两个通道的干扰观测器;设计基于干扰观测器的连续滑模控制器。本发明消除滑模控制方法中的不连续项的同时保证了该方法的强鲁棒性,使得滑模控制器的输出连续化,彻底消除了滑模控制的抖振问题;在控制器中不使用***的加速度信号,削弱了测量噪声对跟踪性能的恶化,利于在工程实际中运用;在***同时存在匹配和不匹配干扰和强非线性的情况下使用滑模控制方法依旧获得了渐近跟踪的稳态性能,保证了电液位置伺服***良好的控制性能。
Description
技术领域
本发明属于电液伺服控制技术领域,特别是一种基于干扰补偿的电液位置伺服***连续滑模控制方法。
背景技术
液压伺服***具有功重比大、响应快及抗负载刚性强等突出优点,在众多重要领域内得到广泛运用。电液伺服***是一个典型的非线性***,包含许多非线性特性和建模不确定性干扰。随着电液伺服***向高精度、高频响发展时,***呈现的非线性特性对***性能的影响越显著,而且干扰的存在会使以***名义模型设计的控制器不稳定或降阶,因此电液伺服***非线性特性和干扰是限制***性能提升的重要因素。随着工业及国防领域技术水平的不断进步,以往基于传统线性理论设计的控制器已逐渐不能满足***的高性能需求,因此必须针对电液伺服***中的非线性特性和存在的干扰研究更加先进的非线性控制策略。
针对电液伺服***的干扰和非线性控制的问题,许多方法相继被提出。在液压位置伺服***控制器的设计中,针对电液伺服***存在的多种干扰,反馈线性化控制的基本思想是通过在控制器中对非线性函数进行精确补偿以使误差动态线性化。虽然理论上可以获得完美的渐近跟踪性能,但是实际***的模型是不可能精确已知的,总会存在建模误差,因此会恶化理论分析获得的跟踪性能。而由于***中存在不匹配干扰,传统的滑模控制方法的基本思路是通过增大控制器的鲁棒性来克服不匹配和匹配干扰从而到达滑模面,但是,即使到达滑模面后,***地跟踪误差在不匹配干扰的影响下仍然无法为零,只能得到一个和不匹配干扰上确界相关的一个一致有界的稳态跟踪误差。并且,通过增大不连续项增益的方法来增加控制器的鲁棒性,在实际运用中很可能激发***高频动态,使***失稳。同时,由于滑模控制方法在进行控制器设计中使用了符号函数,从而使得控制器的输出不连续,虽然有的滑模控制策略补偿了***存在的干扰使得不连续项增益大为减少,削弱滑模控制的抖振,但因为控制器表达式中始终存在不连续项,从而使得滑模控制方法始终存在抖振,无法根除抖振,因而传统的滑模控制方法具有很大的工程局限性。
发明内容
本发明的目的在于提供一种基于干扰补偿的电液位置伺服***连续滑模控制方法。
实现本发明目的的技术解决方案为:一种基于干扰补偿的电液位置伺服***连续滑模控制方法,包括以下步骤:
步骤1、建立电液位置伺服***的数学模型;
步骤2、设计电液位置伺服***两通道的干扰观测器;
步骤3、设计基于干扰观测器的连续滑模控制器。
本发明与现有技术相比,其显著优点是:
(1)本发明消除滑模控制策略中的不连续项的同时保证了该方法的强鲁棒性,使得滑模控制器的输出连续化,彻底消除了滑模控制的抖振问题;
(2)在***同时存在匹配和不匹配干扰和强非线性的情况下使用滑模控制方法依旧获得了渐近跟踪的稳态性能,保证了电液位置伺服***良好的控制性能;
(3)本发明设计干扰观测器观测了液压位置伺服***第二通道和第三通道的干扰,并在控制器中将其完全补偿,抑制了干扰对控制性能的影响;
(4)在控制器中不使用***的加速度信号,削弱了测量噪声对跟踪性能的恶化,利于在工程实际中运用。
附图说明
图1为本发明的基于干扰补偿的电液位置伺服***连续滑模控制方法流程图。
图2为本发明的液压位置伺服***的原理图。
图3为本发明的基于干扰补偿的电液位置伺服***连续滑模控制方法原理示意图。
图4为本发明实施例中***输出对期望指令的跟踪曲线图。
图5为本发明实施例中***的位置跟踪误差随时间变化的曲线图。
图6为本发明实施例中基于干扰补偿的电液位置伺服***连续滑模控制器(C‐SMC)作用下和不确定性补偿的滑模控制器(UC‐SMC)作用下及反馈线性化控制器(FLC)作用下的位置跟踪误差随时间变化曲线图。
图7为本发明实施例中基于干扰补偿的电液位置伺服***连续滑模控制器(C‐SMC)作用下***的控制输入随时间变化的曲线图。
图8为本发明实施例中不确定性补偿的滑模控制器(UC‐SMC)作用下***的控制输入随时间变化的曲线图。
图9为本发明实施例中基于干扰补偿的电液位置伺服***连续滑模控制器作用下***的第二通道干扰观测曲线图。
图10为本发明实施例中基于干扰补偿的电液位置伺服***连续滑模控制器作用下***的第二通道干扰观测误差随时间变化曲线图。
图11为本发明实施例中基于干扰补偿的电液位置伺服***连续滑模控制器作用下***的第三通道干扰观测曲线图。
图12为本发明实施例中基于干扰补偿的电液位置伺服***连续滑模控制器作用下***的第三通道干扰观测误差随时间变化曲线图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
结合图1,本发明的液压位置伺服***的不确定性补偿的滑模控制方法,包括以下步骤:
步骤1、建立液压位置伺服***的数学模型;
步骤1‐1、如图2右半部分所示,液压位置伺服***为通过伺服阀控制的液压马达驱动惯性负载的***构成,图2左半部分为伺服阀控制的液压马达的原理示意图;根据牛顿第二定律,惯性负载的运动方程为:
式(1)中m为惯性负载参数;PL为液压马达两腔压差;A为液压马达的排量;B为粘性摩擦系数;为建模误差,包括m、PL、B的名义值与真实值之间的偏差以及外负载干扰;y为惯性负载的位移;为惯性负载的速度,为惯性负载的加速度;t为时间变量;
液压马达左右两腔的压力动态方程为:
式(2)中P1和P2分别为液压马达两腔的压力,和分别为P1和P2的导数;V1=V01+Ay,V2=V02-Ay,V1和V2分别表示液压马达两腔的控制容积;V01和V02分别为液压马达两腔的初始容积;βe为有效油液弹性模量;Ct为内泄漏系数;q1(t)和q2(t)分别为P1和P2动态方程的建模误差;Q1和Q2分别为液压马达的进油腔流量和回油腔流量;Q1和Q2与伺服阀位移xv的关系为:
式(3)中s(xv)的定义为:
其中,kq为流量增益,Cd流量系数;ω为阀芯面积梯度;ρ为油液密度;Ps为供油压力,Pr为回油压力;液压马达两腔压力满足0<Pr<P1<Ps,0<Pr<P2<Ps,|PL|<<PS;
由于考虑伺服阀动态需要安装额外的位移传感器来获取伺服阀阀芯的位移,而且对于跟踪性能只有微小的提升;因此大量相关的研究都忽略伺服阀的动态,假设采用的是高响应的伺服阀,阀芯位移与控制输入近似为比例环节即xv=kiu,故式(3)可以写成
式(5)中kt=kqki代表总的流量增益,ki为伺服阀增益,u为液压位置伺服***输入,
步骤1‐2、定义状态变量:则***的状态方程为:
式(6)中m、B、βe、A、kt、V01、V02和Ct在观测器和控制器的设计中为名义值,其与真实值之间的偏差集中归类到***干扰中处理,在第二通中是d(x,t),在第三通道中是q(t);其中:
由于液压***参数m,B,βe,kt和Ct受各种因素(如温度、组件磨损程度等)影响变化很大,因此为了简化***状态方程,定义:
因为|PL|<<PS,从而g(x)≠0;第二通道干扰d1(x,t)和第三通道干扰d2(x,t)都是有界的,即:|d1(x,t)|≤D1,|d2(x,t)|≤D2,其中D1、D2分别为|d1(x,t)|和|d1(x,t)|上界,都是已知正数,并且d1(x,t)一阶导数存在;则液压位置伺服***模型为
步骤2,设计两个通道的干扰观测器,步骤如下:
步骤2‐1、设计第二通道干扰观测器:
定义第二通道干扰观测器滑模面s1为:
s1=z1-x2 (10)
其中,z1为第二通道干扰观测器内动态;
式(11)中,k1、β1、ε1、p1和q1均为第二通道干扰观测器参数;p1<q1,均为正奇数,k1、β1、ε1均为正数,β1≥D1;
由式(10)、(11)有:
故第二通道干扰d1(x,t)的观测值设计为其表达式如下:
定义第二通道干扰观测器的Lyapunov方程:
由于β1≥D1,故:
又因,若存在一正定函数V0(t)满足下面不等式:
则有,V0(t)在时间ts内收敛到平衡点,其中
式(18)中,α>0,λ>0,0<γ<1;
由式(17)、(18)式有,V1(t)将在有限时间内收敛到平衡点,即存在一个时刻t2,在经过有限时间t2后,s1为零,此时也将收敛到零,而d1(x,t)的观测误差为:
则干扰的观测误差也将在有限时间t2内为0;即经过时间t2后
综上,得到第二通道干扰观测器表达式为:
步骤2‐2、设计第三通道干扰观测器:
定义第三通道干扰观测器滑模面s2为:
s2=z2-x3 (20)
其中,z2为第三通道干扰观测器内部动态;
式(21)中,k2、β2、ε2、p2和q2均为第三通道干扰观测器参数;其中p2<q2,均为正奇数,k2、β2、ε2均为正数,β2≥D2;
由式(20)、(21)有:
故第三通道干扰d2(x,t)的观测值设计为其表达式如下:
定义第三通道干扰观测器的Lyapunov方程:
又因β2>D2,则有,
此时,由式(17)、(18)式有,V2(t)将在有限时间内为零,即存在一个时刻t3,在经过t3后,s2为零,此时也将收敛到零,又因d2(x,t)估计误差
则干扰的估计误差也将在有限时间t3内为零,即经过时间t3后有:
综上,得到第三通道干扰观测器表达式为:
步骤3、设计基于干扰观测器的连续滑模控制器,具体如下:
定义液压伺服***的位置跟踪误差e0(t)、速度的跟踪误差e1(t)、加速度的跟踪误差e2(t)、加速度导数的跟踪误差e3(t):
e0(t)=x1-xd(t) (27)
其中,xd(t)为***参考信号,xd(t)是三阶连续可导的,且***参考位置信号xd(t)、***参考速度信号***参考加速度信号及***参考加速度导数的信号都是有界的;
定义连续滑模控制器滑模面s为:
其中c1、c2、c3为滑模控制器参数,且均大于零,并且使得表达式是Hurwitz的,则有:
设计连续滑模控制器u的表达式为:
其中k3、k4、ζ为控制器参数,且k3>0、k4>0、0<ζ<1。
步骤4、连续滑模控制器的稳定性测试;具体如下:
步骤4‐1、将式(33)代入式(32)有:
定义滑模控制器Lyapunov方程:
则有:
故,由式(17)、(18)可知,V(t)将在有限时间内为零,即存在一个时刻t1,在经过有限时间t1后,s为零,即:
又因s1、s2也是有限时间内为0,t1为s为零的时刻,t2为s1为零的时刻,t3为s2为零的时刻,则存在t4=max{t1,t2,t3},经过t4时刻后有:
又因:
则:
e2(t)+c1e1(t)+c2e0(t)+c3∫e0(t)dt=d1(t) (40)即:
步骤4‐2、当时:
则有,当t→∞有:
故e0(t)=x1-xd(t)在时间趋于无穷的条件下趋于零;
当δ为一已知正数时:
则有,当t→∞有:
故e0(t)=x1-xd(t)在时间趋于无穷的条件下收敛于一个一致稳定界内。
综上可知,针对液压位置伺服***设计的连续滑模控制器可以使***得到全局渐近稳定的结果。调节观测器系数k1、β1、ε1、p1、q1、k2、β2、ε2、p2、q2可以使观测器的跟踪误差在有限时间内趋于零,调节增益c1、c2、c3、k3、k4、ζ可以使***的跟踪误差在时间趋于无穷的条件下趋于零。液压位置伺服***的基于干扰补偿的电液位置伺服***连续滑模控制方法原理示意图如图3所示。
下面结合具体实施例对本发明做进一步说明。
实施例1
结合图1‐图3,为考核所设计的控制器性能,在MATLAB/Simulink中搭建***仿真模型,在仿真中取如下参数,仿真步长设置为0.0002s:
负载转动惯量m=40kg·m2,马达排量A=2×10-4m3/rad,粘性摩擦系数B=80N·m·s/rad,供油压力Ps=7MPa,回油压力Pr=0,油液弹性模量βe=2×108Pa,马达两腔初始容积V01=V02=1×10-3m3,泄漏系数Ct=9×10-12m3/s/Pa,总流量增益压力动态建模误差q1(t)=q2(t)=6×10-6m3·rad/s,外负载干扰f(t)=200[1-exp(-0.1t3)](N·m),添加未建模摩擦项Fc(x2)=20/π[arctan(100x2)]N·m。
给定***的期望指令为xd(t)=arctan(sin(t))[1‐exp(‐0.01t3)]rad。
取如下的控制器以作对比:
连续滑模控制器(Continuous‐Sliding Mode Controller,C‐SMC):取第二通道干扰观测器参数k1=10000,β1=250,ε1=0.06,p1=3,q1=5,第三通道干扰观测器参数k2=10000,β2=60,ε2=0.05,p2=5,q2=7,控制器参数c1=24,c2=192,c3=512,k3=3、k4=5、ζ=0.5。
匹配和不匹配不确定性补偿的滑模控制器(Uncertainties Compensation‐Sliding Mode Controller,UC‐SMC):不匹配不确定性观测器参数k1=10000,β1=250,ε1=0.06,p1=3,q1=5,第三通道干扰观测器参数k2=10000,β2=60,ε2=0.05,p2=5,q2=7,控制器参数c1=24,c2=192,c3=512,k=2。
反馈线性化控制器(Feedback Linearization Controller,FLC):l1=20、l2=20、l3=20。
连续滑模控制器作用下电液位置伺服***的位置输出对期望指令的跟踪曲线如图4所示,跟踪误差如图5所示;由图4可知,***的位置输出曲线基本与期望指令曲线重合,说明连续滑模控制器具有较好的控制精度,在该控制器的作用下,***表现出良好的跟踪性能;由图5可知,在连续滑模控制器作用下,液压伺服***的位置输跟踪精度很高,稳态跟踪误差收敛到幅值为2×10-3(rad)的界内。
图6为连续滑模控制器(Continuous‐Sliding Mode Controller,C‐SMC)、匹配和不匹配不确定性补偿的滑模控制器(Uncertainties Compensation‐Sliding ModeController,UC‐SMC)、反馈线性化控制器(Feedback Linearization Controller,FLC)分别作用下***的跟踪误差曲线;由图6可知,在***存在未建模干扰时,反馈线性化控制器的控制性能比连续滑模控制器和不确定性补偿的滑模控制器差,因为电液位置伺服***的模型不可能精确已知,总会存在未建模的干扰,而存在的未建模干扰将严重恶化了反馈线性化控制器的控制性能。而连续滑模控制器和不确定性补偿的滑模控制器因为滑模面参数选取是一样的,故两者具有相同的误差动态曲线,同时说明连续滑模控制器和不确定性补偿的滑模控制器具有一样的强鲁棒性。
图7为连续滑模控制器的控制输入曲线,图8为不确定性补偿的滑模控制器的输入曲线;由图7,图8及其对应的局部放大图对比可知,连续滑模控制器的控制器输入曲线为一低频连续曲线,便于在工程实际中运用,而不确定性补偿的滑模控制器虽然控制输入曲线的抖振不大,但是依然存在高频抖振,这在实际的工程运用中,易激发***潜在的未建模高频动态,导致***发散失稳。同时结合图6可知,虽然连续滑模控制器的控制输入曲线不存在高频抖振,但是仍然和不确定性补偿的滑模控制器具有一样的强鲁棒性。
图9为第二通道干扰和第二通道干扰观测曲线,图10为第二通道干扰和第二通道干扰观测误差随时间变化曲线,由图9可知,实际存在的第二通道干扰曲线和其观测曲线基本重合,所设计的干扰观测器对液压位置伺服***的第二通道干扰估计非常准确,同时结合图10可知,第二通道干扰的估计误差在经过很短的时间后迅速收敛到零。
图11为第三通道干扰和第三通道干扰观测曲线,图12为第三通道干扰和第三通道干扰观测误差随时间变化曲线;在图11中,实际存在的第三通道干扰曲线和其估计曲线基本重合,而由图12可知,其估计误差在有限的时间内迅速收敛到零,说明所设计的干扰观测器对***存在的匹配干扰具有较好的估计能力。
Claims (2)
1.一种基于干扰补偿的电液位置伺服***连续滑模控制方法,其特征在于,包括以下步骤:
步骤1、建立电液位置伺服***的数学模型;具体如下:
步骤1-1、电液位置伺服***为通过伺服阀控制的液压马达驱动惯性负载的***;根据牛顿第二定律,惯性负载的运动方程为:
<mrow>
<mi>m</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>P</mi>
<mi>L</mi>
</msub>
<mi>A</mi>
<mo>-</mo>
<mi>B</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>,</mo>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式(1)中m为惯性负载参数;PL为液压马达两腔压差;A为液压马达的排量;B为粘性摩擦系数;为建模误差,包括m、PL、B的名义值与真实值之间的偏差以及外负载干扰;y为惯性负载的位移;为惯性负载的速度,为惯性负载的加速度;t为时间变量;
液压马达左右两腔的压力动态方程为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>P</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>&beta;</mi>
<mi>e</mi>
</msub>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>A</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>t</mi>
</msub>
<msub>
<mi>P</mi>
<mi>L</mi>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>P</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>&beta;</mi>
<mi>e</mi>
</msub>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<mi>A</mi>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<msub>
<mi>C</mi>
<mi>t</mi>
</msub>
<msub>
<mi>P</mi>
<mi>L</mi>
</msub>
<mo>-</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Q</mi>
<mn>2</mn>
</msub>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式(2)中P1和P2分别为液压马达两腔的压力,和分别为P1和P2的导数;V1=V01+Ay,V2=V02-Ay,V1和V2分别表示液压马达两腔的控制容积;V01和V02分别为液压马达两腔的初始容积;βe为有效油液弹性模量;Ct为内泄漏系数;q1(t)和q2(t)分别为P1和P2动态方程的建模误差;Q1和Q2分别为液压马达的进油腔流量和回油腔流量;Q1和Q2与伺服阀位移xv的关系为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>q</mi>
</msub>
<msub>
<mi>x</mi>
<mi>v</mi>
</msub>
<mo>&lsqb;</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>P</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
</mrow>
</msqrt>
<mo>+</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>r</mi>
</msub>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>q</mi>
</msub>
<msub>
<mi>x</mi>
<mi>v</mi>
</msub>
<mo>&lsqb;</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>r</mi>
</msub>
</mrow>
</msqrt>
<mo>+</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>P</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式(3)中s(xv)的定义为:
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
<mi> </mi>
<msub>
<mi>x</mi>
<mi>v</mi>
</msub>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
<mi> </mi>
<msub>
<mi>x</mi>
<mi>v</mi>
</msub>
<mo><</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,kq为流量增益,Cd流量系数;ω为阀芯面积梯度;ρ为油液密度;Ps为供油压力,Pr为回油压力;液压马达两腔压力满足0<Pr<P1<Ps,0<Pr<P2<Ps,|PL|<<PS;
阀芯位移与控制输入近似为比例环节即xv=kiu,故式(3)写成:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>t</mi>
</msub>
<mi>u</mi>
<mo>&lsqb;</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>P</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
</mrow>
</msqrt>
<mo>+</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>r</mi>
</msub>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Q</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>t</mi>
</msub>
<mi>u</mi>
<mo>&lsqb;</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>r</mi>
</msub>
</mrow>
</msqrt>
<mo>+</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>P</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
</mrow>
</msqrt>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
式(5)中kt=kqki代表总的流量增益,ki为伺服阀增益,u为液压位置伺服***输入,
步骤1-2、定义状态变量:则***的状态方程为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>3</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A&beta;</mi>
<mi>e</mi>
</msub>
<msub>
<mi>k</mi>
<mi>t</mi>
</msub>
</mrow>
<mi>m</mi>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mi>u</mi>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>A</mi>
<mn>2</mn>
</msup>
<msub>
<mi>&beta;</mi>
<mi>e</mi>
</msub>
</mrow>
<mi>m</mi>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&beta;</mi>
<mi>e</mi>
</msub>
<msub>
<mi>C</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>+</mo>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式(6)中m、B、A、βe、kt、V01、V02和Ct在观测器和控制器的设计中为名义值,其与真实值之间的偏差集中归类到***干扰中处理,在第二通道中是d(x,t),在第三通道中是q(t);其中:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>b</mi>
<mo>=</mo>
<mfrac>
<mi>B</mi>
<mi>m</mi>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>m</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>A&beta;</mi>
<mi>e</mi>
</msub>
</mrow>
<mi>m</mi>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>P</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
</mrow>
</msqrt>
<mo>+</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>r</mi>
</msub>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>r</mi>
</msub>
</mrow>
</msqrt>
<mo>+</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<msub>
<mi>P</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
为简化***状态方程,定义:
因为|PL|<<PS,则g(x)≠0;第二通道干扰d1(x,t)和第三通道干扰d2(x,t)都是有界的,即:|d1(x,t)|≤D1,|d2(x,t)|≤D2,其中D1、D2分别为|d1(x,t)|和|d1(x,t)|上界,都是已知正数,并且d1(x,t)一阶导数存在;则液压位置伺服***模型为:
步骤2、设计电液位置伺服***两通道的干扰观测器;步骤如下:
步骤2-1、设计第二通道干扰观测器:
定义第二通道干扰观测器滑模面s1为:
s1=z1-x2 (10)
其中,z1为第二通道干扰观测器内部动态;
<mrow>
<msub>
<mover>
<mi>z</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</msubsup>
<mo>-</mo>
<mo>|</mo>
<mrow>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>|</mo>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
式(11)中,k1、β1、ε1、p1和q1均为第二通道干扰观测器参数;p1<q1,均为正奇数,k1、β1、ε1均为正数,β1≥D1;
由式(10)、(11)有:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>s</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mover>
<mi>z</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</msubsup>
<mo>-</mo>
<mo>|</mo>
<mrow>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>|</mo>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
故第二通道干扰d1(x,t)的观测值设计为其表达式如下:
<mrow>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</msubsup>
<mo>-</mo>
<mo>|</mo>
<mrow>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>|</mo>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
定义第二通道干扰观测器的Lyapunov方程:
<mrow>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
由于β1≥D1,故:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>V</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<msub>
<mover>
<mi>s</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>&lsqb;</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</msubsup>
<mo>-</mo>
<mo>|</mo>
<mrow>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>|</mo>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>/</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</msubsup>
<mo>-</mo>
<mo>|</mo>
<mrow>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<mo>+</mo>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&le;</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>/</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</msubsup>
<mo>+</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>&le;</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>/</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</msubsup>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>V</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</msup>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>V</mi>
<mn>1</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>/</mo>
<mn>2</mn>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
又因,若存在一正定函数V0(t)满足下面不等式:
<mrow>
<msub>
<mover>
<mi>V</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&alpha;V</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>&lambda;V</mi>
<mn>0</mn>
<mi>&gamma;</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mn>0</mn>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>t</mi>
<mo>></mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
则有,V0(t)在时间ts内收敛到平衡点,其中,
<mrow>
<msub>
<mi>t</mi>
<mi>s</mi>
</msub>
<mo>&le;</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>+</mo>
<mi>&gamma;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mi>l</mi>
<mi>n</mi>
<mfrac>
<mrow>
<msubsup>
<mi>&alpha;V</mi>
<mn>0</mn>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&gamma;</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&lambda;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
式(18)中,α>0,λ>0,0<γ<1;
由式(17)、(18)式有,V1(t)将在有限时间内收敛到平衡点,即存在一个时刻t2,在经过有限时间t2后,s1为零,此时也将收敛到零,而d1(x,t)的观测误差为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>d</mi>
<mo>~</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</msubsup>
<mo>-</mo>
<mo>|</mo>
<mrow>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>|</mo>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</msubsup>
<mo>-</mo>
<mo>|</mo>
<mrow>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>|</mo>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</msubsup>
<mo>-</mo>
<mo>|</mo>
<mrow>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>|</mo>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msub>
<mover>
<mi>z</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msub>
<mover>
<mi>s</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>19</mn>
<mo>)</mo>
</mrow>
</mrow>
则干扰的观测误差也将在有限时间t2内为0;即经过时间t2后
综上,得到第二通道干扰观测器表达式为:
<mrow>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&beta;</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>s</mi>
<mn>1</mn>
<mrow>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>/</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</msubsup>
<mo>-</mo>
<mo>|</mo>
<mrow>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>|</mo>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
</mrow>
步骤2-2、设计第三通道干扰观测器:
定义第三通道干扰观测器滑模面s2为:
s2=z2-x3 (20)
其中,z2为第三通道干扰观测器内部动态;
式(21)中,k2、β2、ε2、p2和q2均为第三通道干扰观测器参数;其中p2<q2,均为正奇数,k2、β2、ε2均为正数,β2≥D2;
由式(20)、(21)有:
故第三通道干扰d2(x,t)的观测值设计为其表达式如下:
定义第三通道干扰观测器的Lyapunov方程:
<mrow>
<msub>
<mi>V</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msubsup>
<mi>s</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>24</mn>
<mo>)</mo>
</mrow>
</mrow>
又因β2>D2,则有,
由式(17)、(18)式有,V2(t)将在有限时间内为零,即存在一个时刻t3,在经过t3后,s2为零,此时也将收敛到零,又因d2(x,t)估计误差
则干扰的估计误差也将在有限时间t3内为零,即经过时间t3后有:
综上,得到第三通道干扰观测器表达式为:
步骤3、设计基于干扰观测器的连续滑模控制器,具体如下:
定义液压伺服***的位置跟踪误差e0(t)、速度的跟踪误差e1(t)、加速度的跟踪误差e2(t)、加速度导数的跟踪误差e3(t):
e0(t)=x1-xd(t) (27)
<mrow>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>28</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>e</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>29</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>e</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&CenterDot;&CenterDot;&CenterDot;</mo>
</mover>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>30</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,xd(t)为***参考信号,xd(t)是三阶连续可导的,且***参考位置信号xd(t)、***参考速度信号***参考加速度信号及***参考加速度导数的信号都是有界的。
定义连续滑模控制器滑模面s为:
<mrow>
<mi>s</mi>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<msub>
<mi>e</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<mo>&Integral;</mo>
<msub>
<mi>e</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>31</mn>
<mo>)</mo>
</mrow>
</mrow>
其中c1、c2、c3为滑模控制器参数,且均大于零,并且使得表达式是Hurwitz的,则有:
设计连续滑模控制器u的表达式为:
其中k3、k4、ζ为控制器参数,且k3>0、k4>0、0<ζ<1。
2.根据权利要求1所述的基于干扰补偿的电液位置伺服***连续滑模控制方法,其特征在于,对步骤3设计的连续滑模控制器进行稳定性测试,具体如下:
将式(33)代入式(32)有:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>s</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mover>
<mi>s</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>s</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>3</mn>
</msub>
<mi>s</mi>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>4</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>|</mo>
<mi>s</mi>
<mo>|</mo>
</mrow>
<mi>&zeta;</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>3</mn>
</msub>
<mi>s</mi>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>4</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>|</mo>
<mi>s</mi>
<mo>|</mo>
</mrow>
<mi>&zeta;</mi>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>34</mn>
<mo>)</mo>
</mrow>
</mrow>
定义滑模控制器Lyapunov方程:
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>35</mn>
<mo>)</mo>
</mrow>
</mrow>
则有:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>V</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>s</mi>
<mover>
<mi>s</mi>
<mo>&CenterDot;</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>4</mn>
</msub>
<mi>s</mi>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>5</mn>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<msup>
<mrow>
<mo>|</mo>
<mi>s</mi>
<mo>|</mo>
</mrow>
<mi>&zeta;</mi>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>k</mi>
<mn>4</mn>
</msub>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>&zeta;</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<msub>
<mi>k</mi>
<mn>5</mn>
</msub>
<msup>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>&zeta;</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>36</mn>
<mo>)</mo>
</mrow>
</mrow>
故,由式(17)、(18)可知,V(t)将在有限时间内为零,即存在一个时刻t1,在经过有限时间t1后,s为零,即:
<mrow>
<mi>s</mi>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<msub>
<mi>e</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<mo>&Integral;</mo>
<msub>
<mi>e</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mo>-</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>s</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>37</mn>
<mo>)</mo>
</mrow>
</mrow>
又因s1、s2也是有限时间内为0,t1为s为零的时刻,t2为s1为零的时刻,t3为s2为零的时刻,则存在t4=max{t1,t2,t3},经过t4时刻后有:
<mrow>
<mi>s</mi>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<msub>
<mi>e</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<mo>&Integral;</mo>
<msub>
<mi>e</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>38</mn>
<mo>)</mo>
</mrow>
</mrow>
又因:
<mrow>
<msub>
<mi>x</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>bx</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>e</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>39</mn>
<mo>)</mo>
</mrow>
</mrow>
则:
e2(t)+c1e1(t)+c2e0(t)+c3∫e0(t)dt=d1(t) (40)
即:
<mrow>
<msub>
<mover>
<mi>e</mi>
<mo>&CenterDot;&CenterDot;&CenterDot;</mo>
</mover>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<msub>
<mover>
<mi>e</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<msub>
<mover>
<mi>e</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<msub>
<mi>e</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>d</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>41</mn>
<mo>)</mo>
</mrow>
</mrow>
当时:
则有,当t→∞有:
<mrow>
<msub>
<mover>
<mi>e</mi>
<mo>&CenterDot;&CenterDot;&CenterDot;</mo>
</mover>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<msub>
<mover>
<mi>e</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<msub>
<mover>
<mi>e</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<msub>
<mi>e</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>42</mn>
<mo>)</mo>
</mrow>
</mrow>
故e0(t)=x1-xd(t)在时间趋于无穷的条件下趋于零;
当δ为一已知正数时:
则有,当t→∞有:
<mrow>
<msub>
<mover>
<mi>e</mi>
<mo>&CenterDot;&CenterDot;&CenterDot;</mo>
</mover>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<msub>
<mover>
<mi>e</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<msub>
<mover>
<mi>e</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>c</mi>
<mn>3</mn>
</msub>
<msub>
<mi>e</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mi>&delta;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>43</mn>
<mo>)</mo>
</mrow>
</mrow>
故e0(t)=x1-xd(t)在时间趋于无穷的条件下收敛于一个一致稳定的界内。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510522508.2A CN105068426B (zh) | 2015-08-24 | 2015-08-24 | 基于干扰补偿的电液位置伺服***连续滑模控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510522508.2A CN105068426B (zh) | 2015-08-24 | 2015-08-24 | 基于干扰补偿的电液位置伺服***连续滑模控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105068426A CN105068426A (zh) | 2015-11-18 |
CN105068426B true CN105068426B (zh) | 2018-03-09 |
Family
ID=54497817
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510522508.2A Active CN105068426B (zh) | 2015-08-24 | 2015-08-24 | 基于干扰补偿的电液位置伺服***连续滑模控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105068426B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107165892B (zh) * | 2017-06-14 | 2018-05-29 | 武汉科技大学 | 一种电液伺服***的滑模控制方法 |
CN109541945B (zh) * | 2019-01-10 | 2021-11-02 | 中国科学院光电技术研究所 | 一种基于复合型扰动观测器的扰动抑制方法 |
CN112068447A (zh) * | 2020-09-30 | 2020-12-11 | 长沙航空职业技术学院 | 大型装备电液***高精度鲁棒位置控制方法 |
CN112486021B (zh) * | 2020-12-07 | 2021-10-08 | 燕山大学 | 一种针对非对称伺服液压位置跟踪***的低复杂控制方法 |
CN115291507B (zh) * | 2022-01-06 | 2023-05-26 | 兰州理工大学 | 一种矿山充填料浆浓度滑模控制方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104238572A (zh) * | 2014-07-23 | 2014-12-24 | 南京理工大学 | 基于扰动补偿的电机伺服***无抖动滑模位置控制方法 |
CN104360635A (zh) * | 2014-10-21 | 2015-02-18 | 南京理工大学 | 一种电机位置伺服***的抗干扰控制方法 |
CN104614984A (zh) * | 2014-11-20 | 2015-05-13 | 南京理工大学 | 一种电机位置伺服***的高精度控制方法 |
CN104698844A (zh) * | 2015-02-09 | 2015-06-10 | 南京理工大学 | 液压位置伺服***的不确定性补偿的滑模控制方法 |
CN104808653A (zh) * | 2015-04-24 | 2015-07-29 | 南京理工大学 | 基于滑模的电机伺服***加性故障检测和容错控制方法 |
-
2015
- 2015-08-24 CN CN201510522508.2A patent/CN105068426B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104238572A (zh) * | 2014-07-23 | 2014-12-24 | 南京理工大学 | 基于扰动补偿的电机伺服***无抖动滑模位置控制方法 |
CN104360635A (zh) * | 2014-10-21 | 2015-02-18 | 南京理工大学 | 一种电机位置伺服***的抗干扰控制方法 |
CN104614984A (zh) * | 2014-11-20 | 2015-05-13 | 南京理工大学 | 一种电机位置伺服***的高精度控制方法 |
CN104698844A (zh) * | 2015-02-09 | 2015-06-10 | 南京理工大学 | 液压位置伺服***的不确定性补偿的滑模控制方法 |
CN104808653A (zh) * | 2015-04-24 | 2015-07-29 | 南京理工大学 | 基于滑模的电机伺服***加性故障检测和容错控制方法 |
Non-Patent Citations (2)
Title |
---|
基于干扰观测器的非线性不确定***自适应滑模控制;于靖,等;《控制理论与应用》;20140831;第31卷(第8期);第993-995页 * |
基于滑模状态观测器的电液位置伺服***控制;乔继红;《计算机仿真》;20141130;第31卷(第11期);第367-371页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105068426A (zh) | 2015-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105068426B (zh) | 基于干扰补偿的电液位置伺服***连续滑模控制方法 | |
CN104698844B (zh) | 液压位置伺服***的不确定性补偿的滑模控制方法 | |
CN104345639B (zh) | 一种电液位置伺服***鲁棒自适应控制方法 | |
CN104635490B (zh) | 一种单出杆液压缸位置伺服***的输出反馈控制方法 | |
US11454177B2 (en) | Method of aero-engine on-line optimization and multivariable control based on model prediction | |
CN104345638B (zh) | 一种液压马达位置伺服***的自抗扰自适应控制方法 | |
CN104111607B (zh) | 一种考虑输入时滞的电机位置伺服***的控制方法 | |
CN105093936B (zh) | 电液位置伺服***的有限时间连续滑模控制方法 | |
CN105700347B (zh) | 一种含磁滞补偿的液压马达预设性能跟踪控制方法 | |
CN106483844B (zh) | 基于非线性鲁棒的电液伺服***自适应位置控制器的实现方法 | |
CN104317198A (zh) | 具有时变输出约束的电液伺服***非线性鲁棒位置控制方法 | |
CN110308651B (zh) | 基于扩张状态观测器的电液伺服***全状态约束控制方法 | |
CN104360635A (zh) | 一种电机位置伺服***的抗干扰控制方法 | |
CN104932259A (zh) | 一种电液位置伺服***的增益自调节的超螺旋滑模控制方法 | |
CN104898428A (zh) | 一种基于干扰估计的电液伺服***自适应鲁棒控制方法 | |
CN107165892A (zh) | 一种电液伺服***的滑模控制方法 | |
CN108107728A (zh) | 一种基于干扰补偿的电液位置伺服***控制方法 | |
CN107340714A (zh) | 一种纳米伺服***的自适应逆控制方法 | |
CN114296346A (zh) | 基于模糊自抗扰控制的电液位置伺服***重构控制方法及装置 | |
Man et al. | Hydraulic impulse-testing system with pressure difference feedback from a second-stage valve | |
Zhao et al. | A novel method of parameter identification of nonlinear electrohydraulic servo systems | |
Liang-kuan et al. | Compound control strategy for MDF continuous hot pressing electrohydraulic servo system with uncertainties and input saturation | |
CN105739311A (zh) | 基于预设回声状态网络的机电伺服***受限控制方法 | |
CN108345268B (zh) | 考虑输入时滞约束的电液伺服***位置跟踪控制方法 | |
Chen | A backstepping sliding mode control algorithm of electro-hydraulic position servo system of valve-controlled symmetric cylinder based on extended state observer |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20151118 Assignee: Liheng Technology (Jiangsu) Co.,Ltd. Assignor: NANJING University OF SCIENCE AND TECHNOLOGY Contract record no.: X2022980020673 Denomination of invention: A Continuous Sliding Mode Control Method for Electro hydraulic Position Servo System Based on Interference Compensation Granted publication date: 20180309 License type: Exclusive License Record date: 20221103 |