CN105388758B - 一种基于Delta算子的液位控制系统的自适应滑模控制方法 - Google Patents
一种基于Delta算子的液位控制系统的自适应滑模控制方法 Download PDFInfo
- Publication number
- CN105388758B CN105388758B CN201510776263.6A CN201510776263A CN105388758B CN 105388758 B CN105388758 B CN 105388758B CN 201510776263 A CN201510776263 A CN 201510776263A CN 105388758 B CN105388758 B CN 105388758B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- mover
- 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.)
- Expired - Fee Related
Links
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000000694 effects Effects 0.000 claims abstract description 9
- 238000009415 formwork Methods 0.000 claims abstract description 7
- 239000007788 liquid Substances 0.000 claims abstract description 4
- 238000013461 design Methods 0.000 claims description 26
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000013016 damping Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 238000004886 process control Methods 0.000 abstract description 4
- 230000008901 benefit Effects 0.000 abstract description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 14
- 230000008859 change Effects 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000012369 In process control Methods 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005183 dynamical system Methods 0.000 description 2
- 238000010965 in-process control Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 238000001467 acupuncture Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 210000001258 synovial membrane Anatomy 0.000 description 1
- 239000010331 tianhuang Substances 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于Delta算子的液位控制系统的自适应滑模控制方法。在液位控制这类过程控制系统中,存在不确定性和时滞,这些都会影响系统的稳定性,基于delta算子,结合自适应控制,提出一种滑模控制方法。设计非线性滑模面,以线性矩阵不等式的形式给出时滞不确定系统渐进稳定的充分条件,考虑到输入扰动,用自适应边界估计的方法估计上界,从而最终构成完整的自适应滑模控制器。本发明通过设计一种特殊的非线性滑模面,改善了系统的运动品质;结合自适应控制和滑模控制的优点,所设计的控制律能够在系统具有时滞和不确定性有良好的控制效果。本发明用于带有时变时滞和系统建模不确定性的液位控制系统。
Description
技术领域
本发明涉及一种基于Delta算子的液位控制系统的自适应滑模控制方法,属于过程控制领域。
背景技术
在现代工业中,自动控制技术被广泛的应用,不少控制系统以温度、压力、流量、液位和成分等工艺参数作为被控变量,实现对系统的稳定控制。另一方面,随着计算机普及,数字控制器越来越多的被应用工业系统控制中,因此对系统建模的准确性也有了越来越高的要求。然而,各类过程控制系统作为一个复杂的被控对象,具有多输入多输出以及非线性、强耦合、时滞等各种各样的复杂问题。因此控制器需要在系统存在时滞和不确定性的情况下仍然具备较强的控制能力。
近些年,针对系统及控制算法复杂程度的大幅提高,能够保证系统响应速度与稳定性能的高速采样控制器受到日益增多的关注。然而,传统移位算子在高速采样会使离散系统的控制性能变差,甚至出现不稳定的现象。为此,Goodwin教授提出了delta算子方法,可以有效解决传统移位算子存在的问题。delta算子系统是连续时间系统和离散时间系统的统一描述形式,非常适合处理复杂系统、网络系统中的高速采样控制问题。
由于滑模控制的滑动模态对系统参数摄动和外加干扰有完全的自适应性,因此非常适合处理四旋翼直升机飞控系统的控制问题。它的控制是不连续的,控制过程中,闭环系统的结构不停的变化,迫使系统状态沿着预先设计好的滑模面运动,渐渐“滑”向状态平衡点,即渐近稳定。其最主要的优点是一旦系统状态量到达滑模面,系统便不受参数变化和外界扰动的影响。而自适应控制可以处理未知的不确定性和扰动的上界,系统本身不断地检测运行参数,控制效果不断地改变,使系统处于最优工作状态。两者广泛用于飞控系统中,为飞控系统的控制提供了新思路。
为了有效处理控制系统中存在的时滞和不确定性,近些年,研究者提出了很多有效的方法。高存臣则研究了时滞离散系统的控制问题;张端金针对系统控制和信号处理研究了delta算子方法;张彩虹针对delta算子不确定系统提出了一种滑膜控制方法。但现有方法对于同时含有时滞和不确定性的液位控制系统却很难有很好的控制效果,因此本发明有很好的实用性。
发明内容
发明目的:针对上述现有技术,一种基于Delta算子的液位控制系统的自适应滑模控制方法,构造特殊的非线性滑模面,有效的改善了系统的运动品质,以线性矩阵不等式给出系统渐进稳定的充分条件,结合自适应控制,设计滑模控制器,有效克服不确定性和扰动带来的影响。
技术方案:一种基于Delta算子的液位控制系统的自适应滑模控制方法,其特征在于:考虑过程控制系统中存在时变时滞和建模的不确定性,结合自适应控制,提出一种自适应滑模控制方法,使得控制系统在发生执行器有输入扰动的情况下能够稳定运行,根据所获取的模型参数,设计一种特殊的非线性滑模面,求解出系统的滑动模态,通过求解线性矩阵不等式使得系统滑动模态渐进稳定,进而结合自适应边界估计和等效控制设计出滑模控制律,最终构成控制器,包括如下具体步骤:
步骤1)获取系统数学模型:
其中,为系统状态变量,u(k)为控制输入,y(k)为可测输出,τk是不确定时滞项,满足0≤τm≤τk≤τM,τm,τM都是正数,都是常值矩阵,可控,且列满秩,为系统的建模不精确性以及参数摄动,φ(k)∈Rm为系统输入扰动;
步骤2)对系统(1)进行变换,存在非奇异变换:
将系统变为:
其中,
Φ(k)=ΔAx(k)+ΔAdx(k-τk)+Bφ(k),满足||Φ(k)||<η1+η2||x(k)||+η3||x(k-τk)||,η1,η2,η3为待定的正常数;
步骤3)根据滑模控制器的设计方法,首先设计滑模面,为了提高复杂系统的运动品质,在传统滑模面中加入非线性,设计如下滑模面:
其中,ψ(y(k),r)=diag{ψ(y,r)1,...,ψ(y,r)m},ψ(y(k),r)i=αie|y(k-1)-r|表示非线性项,r是一个预先给定的参考信号,αi为一调节参数,用于决定最终闭环系统的阻尼比,从滑模面的形式可知,切换函数中非线性部分的存在,有助于增加对控制输入的有效值,从而获得良好的控制效果,在系统的输出值到达预先给定值的过程中,非线性函数的值逐渐变化,当输出值y(k)达到预先设定值r时,非线性项ψ(y(k),r)将逐渐到达设定值,闭环系统的阻尼比也将会被改变,满足上述条件的函数都可以作为非线性函数,因此该非线性切换函数具有一般性,系数F由下面的不等式(5)和(6)确定
β,γ1,γ2>0 (6)
其中,N,L是正定矩阵,Λ11=(T-2)X+τML+ΩAA,
解出J,X后,可得到F=JX-1;
步骤4)根据步骤3)中设计的滑模面,设计控制律;
步骤4.1)由滑动模态的特点,令δs(k)=0,利用等效控制的思想求解出使得系统状态量维持在滑模面上的等效控制律如下:
ux=-(C(k+1)B)-1{[C(k+1)A+(C(k+1)-C(k))/T]×x(k)+C(k+1)Adx(k-τk)} (7)
步骤4.2)还需要设计不连续切换控制,迫使系统状态量向滑模面上运动,由于切换控制律的设计需要一切标称系统参数以外的范数上界信息,根据假设可知,参数摄动以及输入不确定性均有已知上界,而扰动的上界是未知的,因此,根据自适应估计的方法,定义三个自适应变量和用以动态逼近η1、η2和η3,且 和分别为对应估计偏差,取自适应律如下:
其中,θi,εi是待取定的系数,i=1,2,3;
步骤4.3)根据步骤4.2)中不确定项和扰动的估计值,设计如下不连续切换控制律:
其中,ε是一个正常数;
步骤5)综合步骤4.1)和步骤4.3),得出完整的控制律:
步骤6)根据系统状态,选择合适的参数,完成对其的自适应滑模控制。
具有如下优点:
(1)通过构造特殊的非线性滑模面,有助于增加对控制输入的有效值,从而获得良好的控制效果,使得控制系统具有更好的鲁棒性;
(2)利用线性矩阵不等式给出保证系统渐进稳定的充分条件,充分考虑到过程控制在实际运行过程中可能存在的不确定性和时滞现象,使得控制器的设计具有更好的实用性;
(3)引入自适应边界估计的方法估计出输入不确定性和扰动的大小,滑模控制律不断地改变参数,使得系统保守性更小,控制效果更佳。
本发明所用方法作为一种基于Delta算子的液位控制系统的自适应滑模控制方法,具有一定的实际应用价值,易于实现,容错能力强,能够有效提高液位控制系统的安全性。该方法可操作性强,应用方便、可靠。
附图说明
图1是本发明方法自适应滑模控制所基于的硬件结构示意图;
图2是本发明的方法的流程图;
图3是三容水箱的结构示意图;
图4是状态量响应曲线;
图5是控制输入响应曲线;
图6是MATLAB工具搭建的仿真模型。
具体实施方式
下面结合附图对本发明做更进一步的解释。
如图1所示,一种基于Delta算子的液位控制系统的自适应滑模控制方法基于入选装置实现:作为上位机的计算机1,计算机中将会实现自适应滑模算法2,并发出控制信号3;数据采集卡5,采集电压信号,并传递给计算机;三容水箱7,水箱底部装有液压传感器6,输水管道装有比例阀4,阀门比例由计算机进行控制。
如图2所示,一种基于Delta算子的液位控制系统的自适应滑模控制方法,其特征在于:考虑过程控制系统中存在时变时滞和建模的不确定性,结合自适应控制,提出一种自适应滑模控制方法,使得控制系统在发生执行器有输入扰动的情况下能够稳定运行,根据所获取的模型参数,设计一种特殊的非线性滑模面,求解出系统的滑动模态,通过求解线性矩阵不等式使得系统滑动模态渐进稳定,进而结合自适应边界估计和等效控制设计出滑模控制律,最终构成控制器,包括如下具体步骤:
步骤1)获取系统数学模型:
其中,为系统状态变量,u(k)为控制输入,y(k)为可测输出,τk是不确定时滞项,满足0≤τm≤τk≤τM,τm,τM都是正数,都是常值矩阵,可控,且列满秩,为系统的建模不精确性以及参数摄动,φ(k)∈Rm为系统输入扰动;
步骤2)对系统(1)进行变换,存在非奇异变换:
将系统变为:
其中,
Φ(k)=ΔAx(k)+ΔAdx(k-τk)+Bφ(k),满足||Φ(k)||<η1+η2||x(k)||+η3||x(k-τk)||,η1,η2,η3为待定的正常数;
步骤3)根据滑模控制器的设计方法,首先设计滑模面,为了提高复杂系统的运动品质,在传统滑模面中加入非线性,设计如下滑模面:
其中,ψ(y(k),r)=diag{ψ(y,r)1,...,ψ(y,r)m},ψ(y(k),r)i=αie|y(k-1)-r|表示非线性项,r是一个预先给定的参考信号,αi为一调节参数,用于决定最终闭环系统的阻尼比,从滑模面的形式可知,切换函数中非线性部分的存在,有助于增加对控制输入的有效值,从而获得良好的控制效果,在系统的输出值到达预先给定值的过程中,非线性函数的值逐渐变化,当输出值y(k)达到预先设定值r时,非线性项ψ(y(k),r)将逐渐到达设定值,闭环系统的阻尼比也将会被改变,满足上述条件的函数都可以作为非线性函数,因此该非线性切换函数具有一般性,系数F由下面的不等式(5)和(6)确定
β,γ1,γ2>0 (6)
其中,N,L是正定矩阵,Λ11=(T-2)X+τML+ΩAA,
解出J,X后,可得到F=JX-1;
步骤4)根据步骤3)中设计的滑模面,设计控制律;
步骤4.1)由滑动模态的特点,令δs(k)=0,利用等效控制的思想求解出使得系统状态量维持在滑模面上的等效控制律如下:
ux=-(C(k+1)B)-1{[C(k+1)A+(C(k+1)-C(k))/T]×x(k)+C(k+1)Adx(k-τk)} (7)
步骤4.2)还需要设计不连续切换控制,迫使系统状态量向滑模面上运动,由于切换控制律的设计需要一切标称系统参数以外的范数上界信息,根据假设可知,参数摄动以及输入不确定性均有已知上界,而扰动的上界是未知的,因此,根据自适应估计的方法,定义三个自适应变量和用以动态逼近η1、η2和η3,且 和分别为对应估计偏差,取自适应律如下:
其中,θi,εi是待取定的系数,i=1,2,3;
步骤4.3)根据步骤4.2)中不确定项和扰动的估计值,设计如下不连续切换控制律:
其中,ε是一个正常数;
步骤5)综合步骤4.1)和步骤4.3),得出完整的控制律:
步骤6)根据系统状态,选择合适的参数,完成对其的自适应滑模控制。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围,
下面以实际案例仿真说明实施方案的有效性。
采用浙江天煌科技实业有限公司生产的THJS-1三容水箱仿真平台作为具体的算法实验仿真对象。图3是其基本结构示意图。天煌THJS-1三容水箱实验装置由左、中、右三个水箱、储水箱及相应管路(含管路上手动阀)组成,每个水箱都装有液位传感器;三个水箱由两个泵供水,两套动力系统管道支路上都装有压力表来指示和测量管道内压力变化。管路上任一个手动阀都可作为干扰源,用以产生干扰信号;两套动力系统,其中任意一套可作为调节通道,另一套作为干扰通道,也可产生阶跃干扰信号,整个被控对象组成了一个复杂的过程控制系统。
考虑到系统参数摄动以及外界干扰,结合信号传输可能产生的时滞,设定高速采样周期T=0.01秒,应用delta算子系统模型,得到三容水箱单输入单输出的状态空间模型如下:
其中,hi(i=1,2,3)代表着水箱的高度,并取如下的仿真参数:
φ(k)=0.01sin(0.1πk),并设时滞τM=0.36,τm=0秒。
根据步骤2)可以得到非奇异矩阵为:
进而求得各个矩阵系数为:
A21=[0 0.1492],A22=-0.0746,
B=[0 0 1]T,并设C=[1 0 0],可以得到新的系统模型(3)。
根据步骤3)到5),取各参数r=0,α=0.1,θi=1,εi=0.1(i=1,2,3),EΔ为适当维的单位矩阵,取初始时刻系统的状态量矢量为(此时模型中的变量转化成x(k)):
x0=[x1 x2 x3]=[1 1 1]
根据本发明方法,对三容水箱进行自适应滑模控制,图4-图5为自适应滑模控制结果。图4是状态的响应曲线,图5是控制输入的曲线。
由图4可知,当系统发生执行器故障后,在本发明的滑模控制下,水箱液位在较短的时间内趋于稳定,且响应速度快,超调小。由图5的控制曲线可知,即使系统存在时滞和不确定性,自适应滑模控制律变化幅度很小,即该自适应滑模控制律能够很好地保证水箱的控制精度和安全性。
Claims (1)
1.一种基于Delta算子的液位控制系统的自适应滑模控制方法,其特征在于:考虑液位控制系统中存在时变时滞和建模的不确定性,结合自适应控制,提出一种自适应滑模控制方法,使得控制系统在发生执行器有输入扰动的情况下能够稳定运行,根据所获取的模型参数,设计一种特殊的非线性滑模面,求解出系统的滑动模态,通过求解线性矩阵不等式使得系统滑动模态渐进稳定,进而结合自适应边界估计和等效控制设计出滑模控制律,最终构成控制器,包括如下具体步骤:
步骤1)获取系统数学模型:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&delta;</mi>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mover>
<mi>A</mi>
<mo>&OverBar;</mo>
</mover>
<mo>+</mo>
<mi>&Delta;</mi>
<mover>
<mi>A</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>A</mi>
<mo>&OverBar;</mo>
</mover>
<mi>d</mi>
</msub>
<mo>+</mo>
<mi>&Delta;</mi>
<msub>
<mover>
<mi>A</mi>
<mo>&OverBar;</mo>
</mover>
<mi>d</mi>
</msub>
<mo>)</mo>
</mrow>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<msub>
<mi>&tau;</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mover>
<mi>B</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>+</mo>
<mi>&phi;</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<mi>C</mi>
<mo>&OverBar;</mo>
</mover>
<mover>
<mi>x</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,为系统状态变量,u(k)为控制输入,y(k)为可测输出,τk是不确定时滞项,满足0≤τm≤τk≤τM,τm,τM都是正数,都是常值矩阵,可控,且列满秩,为系统的建模不精确性以及参数摄动,φ(k)∈Rm为系统输入扰动;
步骤2)对系统(1)进行变换,存在非奇异变换:
将系统变为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&delta;</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>A</mi>
<mo>)</mo>
</mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mi>d</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&Delta;A</mi>
<mi>d</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<msub>
<mi>&tau;</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>B</mi>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>+</mo>
<mi>&phi;</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>C</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,
Φ(k)=ΔAx(k)+ΔAdx(k-τk)+Bφ(k),满足||Φ(k)||<η1+η2||x(k)||+η3||x(k-τk)||,η1,η2,η3为待定的正常数;
步骤3)根据滑模控制器的设计方法,首先设计滑模面,为了提高复杂系统的运动品质,在传统滑模面中加入非线性,设计如下滑模面:
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<mtable>
<mtr>
<mtd>
<msub>
<mi>C</mi>
<mi>m</mi>
</msub>
</mtd>
<mtd>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
</mtd>
</mtr>
</mtable>
<mo>&rsqb;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mo>&lsqb;</mo>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>F</mi>
<mo>-</mo>
<mi>&psi;</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>,</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<msubsup>
<mi>A</mi>
<mn>12</mn>
<mi>T</mi>
</msubsup>
</mrow>
</mtd>
<mtd>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
</mtd>
</mtr>
</mtable>
<mo>&rsqb;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,ψ(y(k),r)=diag{ψ(y,r)1,...,ψ(y,r)m},ψ(y(k),r)i=αie|y(k-1)-r|表示非线性项,r是一个预先给定的参考信号,αi为一调节参数,用于决定最终闭环系统的阻尼比,从滑模面的形式可知,切换函数中非线性部分的存在,有助于增加对控制输入的有效值,从而获得良好的控制效果,在系统的输出值到达预先给定值的过程中,非线性函数的值逐渐变化,当输出值y(k)达到预先设定值r时,非线性项ψ(y(k),r)将逐渐到达设定值,闭环系统的阻尼比也将会被改变,满足上述条件的函数都可以作为非线性函数,因此该非线性切换函数具有一般性,系数F由下面的不等式(5)和(6)确定
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&Lambda;</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mn>11</mn>
</msub>
<mi>X</mi>
<mo>-</mo>
<msub>
<mi>A</mi>
<mn>12</mn>
</msub>
<mi>J</mi>
<mo>+</mo>
<msub>
<mi>&Omega;</mi>
<mrow>
<mi>A</mi>
<mi>A</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>d</mi>
<mn>11</mn>
</mrow>
</msub>
<mi>X</mi>
<mo>-</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>d</mi>
<mn>12</mn>
</mrow>
</msub>
<mi>J</mi>
</mrow>
</mtd>
<mtd>
<msub>
<mi>&Omega;</mi>
<mrow>
<mi>A</mi>
<mi>H</mi>
</mrow>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<msub>
<mi>&Lambda;</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
</mrow>
</msub>
<mi>X</mi>
<mo>-</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>d</mi>
<mn>12</mn>
</mrow>
</msub>
<mi>J</mi>
<mo>+</mo>
<mi>L</mi>
<mo>/</mo>
<msub>
<mi>&tau;</mi>
<mi>M</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>H</mi>
<mn>1</mn>
</msub>
<mi>X</mi>
<mo>-</mo>
<msub>
<mi>H</mi>
<mn>2</mn>
</msub>
<mi>J</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mo>+</mo>
<msub>
<mi>&Omega;</mi>
<mrow>
<mi>A</mi>
<mi>H</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>XA</mi>
<mn>12</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>N</mi>
<mo>-</mo>
<mi>L</mi>
<mo>/</mo>
<msub>
<mi>&tau;</mi>
<mi>M</mi>
</msub>
</mrow>
</mtd>
<mtd>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>H</mi>
<mn>3</mn>
</msub>
<mi>X</mi>
<mo>-</mo>
<msub>
<mi>H</mi>
<mn>4</mn>
</msub>
<mi>J</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>XA</mi>
<mn>12</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>&beta;</mi>
<mi>I</mi>
<mo>+</mo>
<msub>
<mi>&Omega;</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<mi>I</mi>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<mi>I</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo><</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
β,γ1,γ2>0 (6)
其中,N,L是正定矩阵,Λ11=(T-2)X+τML+ΩAA,
<mrow>
<msub>
<mi>&Lambda;</mi>
<mn>22</mn>
</msub>
<mo>=</mo>
<msubsup>
<mi>XA</mi>
<mn>11</mn>
<mi>T</mi>
</msubsup>
<mo>-</mo>
<msup>
<mi>J</mi>
<mi>T</mi>
</msup>
<msubsup>
<mi>A</mi>
<mn>12</mn>
<mi>T</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>11</mn>
</msub>
<mi>X</mi>
<mo>-</mo>
<msub>
<mi>A</mi>
<mn>12</mn>
</msub>
<mi>J</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&tau;</mi>
<mi>M</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&tau;</mi>
<mi>m</mi>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>N</mi>
<mo>-</mo>
<mi>L</mi>
<mo>/</mo>
<msub>
<mi>&tau;</mi>
<mi>M</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&Omega;</mi>
<mrow>
<mi>A</mi>
<mi>A</mi>
</mrow>
</msub>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mi>&Omega;</mi>
<mrow>
<mi>A</mi>
<mi>A</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<msup>
<mi>&alpha;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>A</mi>
<mn>12</mn>
</msub>
<msubsup>
<mi>A</mi>
<mn>12</mn>
<mi>T</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<msup>
<mi>&alpha;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>A</mi>
<mrow>
<mi>d</mi>
<mn>12</mn>
</mrow>
</msub>
<msubsup>
<mi>A</mi>
<mrow>
<mi>d</mi>
<mn>12</mn>
</mrow>
<mi>T</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&beta;E</mi>
<mi>&Delta;</mi>
</msub>
<msubsup>
<mi>E</mi>
<mi>&Delta;</mi>
<mi>T</mi>
</msubsup>
<mo>,</mo>
<msub>
<mi>&Omega;</mi>
<mrow>
<mi>A</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<msup>
<mi>&alpha;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>A</mi>
<mn>12</mn>
</msub>
<msubsup>
<mi>H</mi>
<mn>2</mn>
<mi>T</mi>
</msubsup>
<mo>+</mo>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<msup>
<mi>&alpha;</mi>
<mn>2</mn>
</msup>
<msub>
<mi>A</mi>
<mrow>
<mi>d</mi>
<mn>12</mn>
</mrow>
</msub>
<msubsup>
<mi>H</mi>
<mn>4</mn>
<mi>T</mi>
</msubsup>
<mo>,</mo>
</mrow>
解出J,X后,可得到F=JX-1;
步骤4)根据步骤3)中设计的滑模面,设计控制律;
步骤4.1)由滑动模态的特点,令δs(k)=0,利用等效控制的思想求解出使得系统状态量维持在滑模面上的等效控制律如下:
ux=-(C(k+1)B)-1{[C(k+1)A+(C(k+1)-C(k))/T]×x(k)+C(k+1)Adx(k-τk)} (7)
步骤4.2)还需要设计不连续切换控制,迫使系统状态量向滑模面上运动,由于切换控制律的设计需要一切标称系统参数以外的范数上界信息,根据假设可知,参数摄动以及输入不确定性均有已知上界,而扰动的上界是未知的,因此,根据自适应估计的方法,定义三个自适应变量和用以动态逼近η1、η2和η3,且 和分别为对应估计偏差,取自适应律如下:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&delta;</mi>
<msub>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mn>1</mn>
</msub>
<msub>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&xi;</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&delta;</mi>
<msub>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mn>2</mn>
</msub>
<msub>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&xi;</mi>
<mn>2</mn>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&delta;</mi>
<msub>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>&theta;</mi>
<mn>3</mn>
</msub>
<msub>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&xi;</mi>
<mn>3</mn>
</msub>
<mo>|</mo>
<mo>|</mo>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<msub>
<mi>&tau;</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,θi,εi是待取定的系数,i=1,2,3;
步骤4.3)根据步骤4.2)中不确定项和扰动的估计值,设计如下不连续切换控制律:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>u</mi>
<mi>&Phi;</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mi>B</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>(</mo>
<mi>&epsiv;</mi>
<mo>+</mo>
<msub>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
<mo>+</mo>
<msub>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mn>3</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<mi>C</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<msub>
<mi>&tau;</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<mo>)</mo>
<mi>sgn</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,ε是一个正常数;
步骤5)综合步骤4.1)和步骤4.3),得出完整的控制律:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>u</mi>
<mo>=</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mi>B</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>{</mo>
<mo>&lsqb;</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>A</mi>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mi>C</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>T</mi>
<mo>&rsqb;</mo>
<mo>&times;</mo>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>A</mi>
<mi>d</mi>
</msub>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>-</mo>
<msub>
<mi>&tau;</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mi>&epsiv;</mi>
<mo>+</mo>
<msub>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<mi>C</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<mo>+</mo>
<msub>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<mi>C</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<mo>+</mo>
<msub>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mn>3</mn>
</msub>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<mi>C</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mo>|</mo>
<mi>x</mi>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<msub>
<mi>&tau;</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<mo>)</mo>
<mi>s</mi>
<mi>g</mi>
<mi>n</mi>
<mo>(</mo>
<mi>s</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>)</mo>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤6)根据系统状态,选择合适的参数,完成对其的自适应滑模控制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510776263.6A CN105388758B (zh) | 2015-11-10 | 2015-11-10 | 一种基于Delta算子的液位控制系统的自适应滑模控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510776263.6A CN105388758B (zh) | 2015-11-10 | 2015-11-10 | 一种基于Delta算子的液位控制系统的自适应滑模控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105388758A CN105388758A (zh) | 2016-03-09 |
CN105388758B true CN105388758B (zh) | 2018-04-03 |
Family
ID=55421150
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510776263.6A Expired - Fee Related CN105388758B (zh) | 2015-11-10 | 2015-11-10 | 一种基于Delta算子的液位控制系统的自适应滑模控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105388758B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106527148B (zh) * | 2016-12-15 | 2019-04-05 | 南京航空航天大学 | 一种离散时滞不确定系统的滑模鲁棒控制方法 |
CN107797451B (zh) * | 2017-10-21 | 2021-01-05 | 辽宁石油化工大学 | 未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法 |
CN108919647A (zh) * | 2018-07-23 | 2018-11-30 | 哈尔滨理工大学 | 一种具有随机非线性扰动的滑模控制方法 |
CN110426957B (zh) * | 2019-07-31 | 2020-03-13 | 深圳信息职业技术学院 | 一种基于时延观测器的水厂投药系统自适应滑模控制方法 |
CN110647136B (zh) * | 2019-09-29 | 2021-01-05 | 华东交通大学 | 一种牵引电机驱动系统复合故障检测与分离方法 |
CN112180727B (zh) * | 2020-10-01 | 2023-03-14 | 天水师范学院 | 一种大时滞非线性系统的多率采样控制方法、系统及应用 |
CN112799374B (zh) * | 2020-12-24 | 2023-01-10 | 南京财经大学 | 一种Delta算子切换粮食管理系统的全阶故障估计观测器的设计方法 |
CN113467528B (zh) * | 2021-07-23 | 2022-12-13 | 上海电力大学 | 基于滑模变结构控制的核电站除氧器水位控制设计方法 |
CN113534777B (zh) * | 2021-07-30 | 2022-08-02 | 淮阴工学院 | 一种针对二级化学反应器变时滞系统的故障估计方法 |
CN113777932B (zh) * | 2021-11-15 | 2022-02-22 | 南京信息工程大学 | 一种基于Delta算子的四旋翼自适应滑模容错控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6185467B1 (en) * | 1998-09-02 | 2001-02-06 | Cirrus Logic, Inc. | Adaptive discrete-time sliding mode controller |
CN102566418A (zh) * | 2012-02-17 | 2012-07-11 | 北京国电蓝天节能科技开发有限公司 | 一种无需物理参数的倒立摆自适应滑模控制方法 |
CN103279038A (zh) * | 2013-06-19 | 2013-09-04 | 河海大学常州校区 | 基于t-s模糊模型的微陀螺仪滑模自适应控制方法 |
CN104932259A (zh) * | 2015-05-20 | 2015-09-23 | 南京理工大学 | 一种电液位置伺服系统的增益自调节的超螺旋滑模控制方法 |
CN104965410A (zh) * | 2015-04-17 | 2015-10-07 | 华南理工大学 | 一种蒸汽发生器水位系统的动态滑模变结构控制器与方法 |
-
2015
- 2015-11-10 CN CN201510776263.6A patent/CN105388758B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6185467B1 (en) * | 1998-09-02 | 2001-02-06 | Cirrus Logic, Inc. | Adaptive discrete-time sliding mode controller |
CN102566418A (zh) * | 2012-02-17 | 2012-07-11 | 北京国电蓝天节能科技开发有限公司 | 一种无需物理参数的倒立摆自适应滑模控制方法 |
CN103279038A (zh) * | 2013-06-19 | 2013-09-04 | 河海大学常州校区 | 基于t-s模糊模型的微陀螺仪滑模自适应控制方法 |
CN104965410A (zh) * | 2015-04-17 | 2015-10-07 | 华南理工大学 | 一种蒸汽发生器水位系统的动态滑模变结构控制器与方法 |
CN104932259A (zh) * | 2015-05-20 | 2015-09-23 | 南京理工大学 | 一种电液位置伺服系统的增益自调节的超螺旋滑模控制方法 |
Non-Patent Citations (1)
Title |
---|
不确定时变时滞系统的自适应全局鲁棒滑模控制;杨蒲 等;《控制与决策》;20140930;第29卷(第9期);第1688-1692页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105388758A (zh) | 2016-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105388758B (zh) | 一种基于Delta算子的液位控制系统的自适应滑模控制方法 | |
CN106773691B (zh) | 基于ls‑svm的高超声速飞行器自适应时变预设性能控制方法 | |
CN102096373B (zh) | 基于增量改进bp神经网络的微波干燥pid控制方法 | |
Sofianos et al. | Robust adaptive multiple models based fuzzy control of nonlinear systems | |
CN105739310A (zh) | 一种基于多模型的伺服系统自适应控制系统 | |
CN106647283A (zh) | 一种基于改进cpso的自抗扰位置伺服系统优化设计方法 | |
CN103927451B (zh) | 一种空间机器人系统参数在轨辨识方法 | |
CN104360596B (zh) | 机电伺服系统有限时间摩擦参数辨识和自适应滑模控制方法 | |
CN107102634A (zh) | 一种基于转台伺服系统的参数估计与跟踪控制方法 | |
CN106054596B (zh) | 一种基于改进性能指标的pid控制器参数优化整定方法 | |
CN103324093B (zh) | 一种多模型自适应控制系统及其控制方法 | |
Xu et al. | Model-free adaptive command-filtered-backstepping sliding mode control for discrete-time high-order nonlinear systems | |
CN106249599A (zh) | 一种基于神经网络预测的网络化控制系统故障检测方法 | |
CN105751470A (zh) | 一种注塑机温度实时控制方法 | |
CN106842948A (zh) | 基于bp网络的hdp分子蒸馏系统的最优控制方法 | |
CN106125553A (zh) | 一种考虑状态约束的液压系统自适应控制方法 | |
CN105892296A (zh) | 一种工业加热炉系统的分数阶动态矩阵控制方法 | |
CN102509152A (zh) | 一种基于rbf神经网络的开关磁阻电机在线建模方法 | |
CN109143853B (zh) | 一种石油炼化过程中分馏塔液位自适应控制方法 | |
CN105676637B (zh) | 基于预测函数控制的分子量输出pdf控制方法 | |
Owa et al. | A wavelet neural network based non-linear model predictive controller for a multi-variable coupled tank system | |
CN108121207A (zh) | 同时满足跟踪性和鲁棒性的燃气锅炉最优分数阶pid控制方法 | |
Deshpande et al. | Design of multi scale PID controller for Temperature process | |
CN106371321A (zh) | 一种焦化炉炉膛压力系统模糊网络优化pid控制方法 | |
CN103558761A (zh) | 一种具有控制器输入饱和的非线性化学反应循环不确定时滞系统的控制方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180403 Termination date: 20211110 |
|
CF01 | Termination of patent right due to non-payment of annual fee |