CN108008628A - 一种不确定欠驱动无人艇***的预设性能控制方法 - Google Patents
一种不确定欠驱动无人艇***的预设性能控制方法 Download PDFInfo
- Publication number
- CN108008628A CN108008628A CN201711143413.5A CN201711143413A CN108008628A CN 108008628 A CN108008628 A CN 108008628A CN 201711143413 A CN201711143413 A CN 201711143413A CN 108008628 A CN108008628 A CN 108008628A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mover
- mfrac
- unmanned boat
- 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/0265—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion
- G05B13/027—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion using neural networks only
-
- 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
-
- 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/048—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 using a predictor
Landscapes
- Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Feedback Control In General (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
本发明公开了一种不确定欠驱动无人艇***的预设性能控制方法,具体包括以下步骤:建立无人艇动态模型,包括运动学和动力学模型;提出预设性能控制法处理跟踪误差约束条件;应用动态面控制技术处理虚拟控制器;采用神经网络技术在线学习***的不确定动态;应用李雅普诺夫稳定性理论设计符合要求的跟踪控制器。本发明的预设性能控制方法,能够避免控制器设计奇异的问题,并且能够满足预先指定的跟踪误差暂态及稳态性能,应用自适应神经网络控制技术在线学习***的不确定动态,并设计控制器完成轨迹跟踪控制目标,提高了控制方法的实用性。
Description
技术领域
本发明属于无人艇的轨迹跟踪控制领域,特别涉及一种不确定欠驱动无人艇***的预设性能控制方法。
背景技术
无人艇是一种用于海上搜救,海洋资源勘探,科学考察,深海作业等领域的水面舰艇,它可以按照预先设定的轨迹及程序完成环境探测等一系列任务。尤其对于一些恶劣的环境或危险的海域,使用无人艇可以有效保障人员安全,符合实际工程应用需求。目前,使用数量最为广泛的无人艇都是欠驱动的,这意味着他们仅拥有为纵向与转向运动而装备的螺旋桨与船舵,而没有直接控制横荡运动的执行器,这对无人艇的跟踪控制设计带来了极大的挑战。由于无人艇控制具有极大的理论挑战性和广泛的实际应用性,无人艇的轨迹跟踪控制逐渐成为一个国内外研究的重要前沿问题。
欠驱动无人艇的轨迹跟踪控制的难点主要在于它是一个三个自由度的机械***,而它的独立控制输入却只有两个(该***的独立控制变量个数小于***自由度个数)。因此,相对于全驱动无人艇的控制设计,欠驱动无人艇的跟踪误差通常需要满足一定的约束条件以避免跟踪控制器设计中可能产生的奇异问题。为了使跟踪误差的约束条件能够在无人艇整个运动过程中得到满足,本专利提出一种预设性能控制的方法,该方法不仅能够满足预先指定的***跟踪误差的暂态及稳态性能,这对于高性能、高精度的无人艇的控制***至关重要,而且有效地避免了跟踪控制器设计中可能产生的奇异问题。此外,由于海洋环境的复杂性与多变性,对无人艇的动力学模型通常难以获得精确的数学模型,尤其是与外界环境息息相关的水动力阻尼非线性动态的准确建模是极其困难的,此外,外界的风、浪以及洋流对无人艇控制***的影响通常也是不可忽视的。本专利在跟踪控制器设计中,充分考虑了***的不确定动态和外部时变扰动的情况,应用自适应神经网络控制技术和扰动观测器方法分别对它们进行在线估计,并在反馈的前向通道中对***的不确定动态及外部时变扰动进行补偿设计。
发明内容
本发明的目的是针对上述现有技术的不足,提出一种不确定欠驱动无人艇***的预设性能控制方法,能够有效避免跟踪控制器设计中可能产生的奇异问题,提供预先指定的跟踪误差的暂态及稳态性能,并结合自适应神经网络控制技术和扰动观测器方法对***的不确定动态和外部时变扰动进行估计和补偿,最终完成跟踪控制器设计。
本发明目的可以通过如下技术方案实现:
本发明公开了一种不确定欠驱动无人艇***的预设性能控制方法,包含如下步骤:
步骤1)、建立无人艇动态模型:包括无人艇的运动学和动力学模型;
步骤2)、提出预设性能控制法处理跟踪误差约束条件:定义无人艇实际轨迹与理想轨迹之间的跟踪误差,根据避免跟踪控制器奇异的条件设定跟踪误差需满足的约束条件,然后提出预设性能控制技术把受到约束条件的跟踪误差转化为不受约束条件的误差;
步骤3)、应用动态面控制技术处理虚拟控制器:采用动态面技术避免虚拟控制器的求导;
步骤4)、采用神经网络在线学习***的不确定动态:利用RBF神经网络对任意光滑函数以任意精度逼近的特性,分别构造第一RBF神经网络和第二RBF神经网络对***不确定动态fr(r)和fu(u)进行在线学习,其中,神经网络的权值更新率设计如下:
其中,和分别为不确定动态fr(r)和fu(u)的神经网络权值向量;和分别为不确定动态fr(r)和fu(u)的神经网络权值向量的转置向量;Γ1和Γ2分别为和所设计的对称正定矩阵;re和ue分别为r,u与各自滤波虚拟控制器之间的误差,所述滤波虚拟控制器为步骤3)中虚拟控制器通过动态控制技术处理后得到的变量;S1(Z1)和S2(Z2)均为高斯函数;σ1和σ2分别为和设计的正常数;r和u分别是无人艇的航向角速度和纵向角速度;Z1是第一RBF神经网络的输入量,数值等于r,Z2是第二RBF神经网络的输入量,数值等于u;Ω1和Ω2分别表示r和u的紧集;
步骤5)、设计跟踪控制器:应用李雅普诺夫稳定性理论设计符合要求的跟踪控制器:
其中,m11,m22和m33是无人艇在u,v,r各自速度方向上的质量,包括附加的质量,u为无人艇的纵向角速度,v为无人艇的横荡速度,r为无人艇的航向角速度,k1和k2分别为跟踪控制器τr和τu的控制增益,z1和z2为轨迹跟踪误差经过步骤2)处理后得到的转化误差,χ1和χ2分别为跟踪控制器τr和τu包含的中间变量,ψe为无人艇实际航向角与参考轨迹航向角的误差,和分别为跟踪控制器τr和τu中用于抵消时变扰动的扰动观测器,yr和yu分别为虚拟控制器αr和αu采用动态面技术而引入的层面误差,c1和c2分别为虚拟控制器αr和αu应用动态面技术时引入的设计参数。
作为优选的技术方案,步骤1)中,所述的无人艇动态模型为:
其中,(x,y)为无人艇的位置,ψ为航向角,u为无人艇的纵向角速度,v为无人艇的横荡速度,r为无人艇的航向角速度,表示无人艇的实际轨迹关于时间t的导数,表示无人艇的各个速度关于时间t的导数,(τu,τr)为无人艇的控制输入,(τwu,τwv,τwr)为外界时变扰动,(fu(u),fv(v),fr(r))为无人艇动力学模型中的不确定动态,具有如下形式:
du表示u方向的第一个水动力阻尼系数,dui表示u方向的第i个水动力阻尼系数,dv表示v方向的第一个水动力阻尼系数,dvi表示v方向的第i个水动力阻尼系数,dr表示r方向的第一个水动力阻尼系数,dri表示r方向的第i个水动力阻尼系数。
作为优选的技术方案,步骤2)中,
所述理想轨迹即参考轨迹满足如下条件:
其中,所述参考轨迹是一个圆,(xd,yd)为参考轨迹的位置,是参考轨迹关于时间t的导数,rd表示参考轨迹的航向角速度,ud表示参考轨迹的纵向角速度;
为了避免跟踪控制器奇异,跟踪误差ψe(t)和ze(t)需满足的约束条件如下:
ψe(t)具体为无人艇实际航向角与参考轨迹航向角的误差,即为ψe,ze(t)具体为无人艇实际位置与参考轨迹位置的误差,即为ze;为了满足上述约束条件,同时能够进一步刻画***的跟踪误差暂态性能及稳态性能要求,将上述的约束条件常数上下界进一步描述成预先指定的时变函数上下界;上述的约束条件刻画成如下更加严格的约束条件:
其中,e 1(t), e 2(t),分别为跟踪误差ψe(t)和ze(t)约束条件的上下界函数,它们是一种按指数衰减的单调递减函数,且满足:
其中,上下界函数具体为:
引入一个严格单调递增的光滑转化函数R(z),将受到约束的跟踪误差转化为不受约束的误差:
设计了如下的转化函数:
其中,γ1为第一个转换函数R1(z1)的中间变量,γ2为第二个转换函数R2(z2)的中间变量,得到如下的转化误差:
作为优选的技术方案,步骤3)中,所述应用动态面控制技术处理虚拟控制器,具体如下:
其中,αfr,αfu分别为r,u的滤波虚拟控制,αr,αu分别为r,u的虚拟控制器,c1,c2分别是αfr,αfu选定的滤波器时间常数;构造虚拟控制器αr和αu如下:
其中,ψd为参考轨迹的航向角,ψd=atan2(ye,xe);为参考轨迹关于时间t的导数;和分别为虚拟控制器αr和αu的设计参数,且 φ1和φ2分别为虚拟控制器αr和αu的中间变量:
其中,分别表示跟踪误差ψe(t)和ze(t)约束条件的上下界函数关于时间t的导数;分别表示γ1、γ2关于时间t的导数。
作为优选的技术方案,步骤4)中,所述Γ1和Γ2分别为和所设计的对称正定矩阵,具体为:Γ1=Γ2=1;所述σ1和σ2分别为和设计的正常数,具体为:σ1=σ2=1;所述re和ue分别为r,u与各自滤波虚拟控制器之间的误差,具体为:re=r-αfr,ue=u-αfu,其中αfr,αfu分别为r,u的滤波虚拟控制。
作为优选的技术方案,步骤5)中,扰动观测器和设计如下:
其中,ξ1和ξ2分别为扰动观测器和的状态;kdr和kdu分别为扰动观测器和的设计参数,且kdr>0,kdu>0。
作为优选的技术方案,步骤5)中,所述χ1和χ2分别为跟踪控制器τr和τu包含的中间变量,具体为:
其中,即为表示跟踪误差ψe(t)约束条件的上界函数;即为表示跟踪误差ze(t)约束条件的上界函数;所述yr和yu分别为虚拟控制器αr和αu采用动态面技术而引入的层面误差,具体为:yr=αfr-αr,yu=αfu-αu,其中,αfr,αfu分别为r,u的滤波虚拟控制,αr,αu分别为r,u的虚拟控制器。
本发明与现有技术相比,具有如下优点和有益效果:
1、本发明提出采用预设性能控制方法,不仅避免了跟踪控制器设计奇异的问题,而且能够满足预先指定的***跟踪误差的暂态及稳态性能,提高了控制***的***的跟踪速度及精度。
2、本发明通过引入动态面控制设计的方法,避免了后推设计第二步对速度的求导,使得所提出跟踪控制器更加易于实现。
3、本发明应用自适应神经网络控制技术在线学习***的不确定动态,并结合扰动观测器来在线估计外部时变扰动,然后在反馈的前向通道中对***的不确定动态及外部时变扰动进行补偿设计,增强了控制***的鲁棒性。
附图说明
图1为发明的无人艇结构示意图。
图2为发明的无人艇的轨迹跟踪控制***结构框图。
图3为发明的无人艇的跟踪误差ψe(t)的变化示意图。
图4为发明的无人艇的跟踪误差ze(t)的变化示意图。
图5为发明的无人艇的神经网络权值估计的二范数示意图。
图6为发明的无人艇轨迹跟踪控制的输出轨迹示意图。
图7为发明的无人艇控制***的跟踪控制器输入τr的示意图。
图8为发明的无人艇控制***的跟踪控制器输入τu的示意图。
图9为发明的无人艇扰动观测器与实际扰动比对的示意图。
图10为发明的无人艇扰动观测器与实际扰动比对的示意图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例
在欠驱动无人艇轨迹跟踪控制中,考虑***模型具有不确定动态的情况,无人艇结构示意图如图1所示,无人艇的轨迹跟踪控制***结构框图如图2所示,一种不确定欠驱动无人艇***的预设性能控制方法,其详细实施过程包括:
步骤1)、建立无人艇动态模型,包括无人艇的运动学和动力学模型:
无人艇动态模型为:
其中,(x,y)为无人艇的位置,ψ为航向角,u为纵向速度,v为无人艇的横荡速度,r为航向角速度,(τu,τr)为无人艇的控制输入,(τwu,τwv,τwr)为外界时变扰动,(fu(u),fv(v),fr(r))为无人艇动力学模型中的不确定动态,且(fu(u),fv(v),fr(r))具有如下形式:
在本实例中,选取的无人艇动态模型的***参数分别为:
m11=120×103kg,m22=177.9×103kg,m33=636×105kg,
du=215×102kg,du2=0.2×du,du3=0.1×du,
dv=147×103kg,dv2=0.2×dv,dv3=0.1×dv,
dr=802×104kg,dr2=0.2×dr,dr3=0.1×dr,
τwu=m11×(1.5+sin(0.1t)+cos(0.01t)),
τwv=m22×(1.5+cos(0.1t)+sin(0.01t)),
τwr=m33×(1.5-sin(0.1t)+cos(0.01t)),
艇体长度为L=38m。
步骤2)、提出预设性能控制法处理跟踪误差约束条件:定义无人艇实际轨迹与理想轨迹之间的跟踪误差,根据避免跟踪控制器奇异的条件给出跟踪误差需满足的约束条件,然后提出预设性能控制技术把受到约束条件的跟踪误差转化为不受约束条件的误差;
理想轨迹即参考轨迹满足如下条件:
其中,xd(0)=0,yd(0)=0,ud=5,rd=0.05,该参考轨迹是一个半径100的圆。无人艇的初始位置及航向角分别选择为[x(0),y(0),ψ(0)]=[0,20,0.6],初始速度选择为[u(0),v(0),r(0)]=[0,0,0]。
为了避免跟踪控制器设计的奇异问题,所述的跟踪误差的满足如下的约束条件:
为了满足上述约束条件,同时能够进一步刻画***的跟踪误差暂态性能和稳态性能要求,将上述的约束条件常数上下界进一步描述成预先指定的时变函数上下界。上述的约束条件刻画成如下更加严格的约束条件:
其中,e 1(t),e2(t),分别为跟踪误差ψe(t)和ze(t)约束条件的上下界,它们是一种按指数衰减的单调递减函数,且满足:
选择如下:
图3,图4分别为无人艇的跟踪误差ψe(t)和ze(t)的变化示意图,可以看到误差的暂态性能和稳态性能都在预先指定的范围内。
然后引入一个严格单调递增的光滑转化函数R(z),将受到约束的跟踪误差转化为不受约束的误差:
设计了如下的转化函数:
其中,之后得到如下的转化误差:
步骤3)、应用动态面控制技术处理虚拟控制器:虚拟控制器求导过程中会引入速度的动态,增加跟踪控制器实现的难度,采用动态面技术可以避免虚拟控制器的求导。引入动态面控制技术处理虚拟控制器:
其中,αfr,αfu分别为r,u的滤波虚拟控制输入,αr,αu分别为r,u的虚拟控制器,c1,c2分别是αfr,αfu选定的滤波器时间常数,c1=c2=0.001,构造虚拟控制器αr和αu如下:
其中,ψd为参考轨迹的航向角,ψd=atan2(ye,xe),为参考轨迹的导数,和分别为虚拟控制器αr和αu的设计参数,φ1和φ2分别为虚拟控制器αr和αu的中间变量:
步骤4)、采用神经网络在线学习***的不确定动态:利用RBF神经网络可以对任意光滑函数以任意精度逼近的特性,分别构造第一RBF神经网络和第二RBF神经网络去逼近***不确定动态fr(r)和fu(u),第一RBF神经网络使用88个节点,中心点分布在[-0.1,0.1]范围内,宽度选为0.03,第二RBF神经网络使用63个节点,中心点分布在[-5,5]范围内,宽度选为0.55,和分别为逼近fr(r)和fu(u)的神经网络权值,权值更新率设计如下:
其中,r和u分别是无人艇的航向角速度和纵向速度,S1(Z1)和S2(Z2)都是高斯函数,Γ1=Γ2=1,σ1=σ2=1,re和ue分别为r,u与各自滤波虚拟控制器之间的误差,re=r-αfr,ue=u-αfu,权值更新率初值选为 图5无人艇的神经网络权值估计的二范数。
步骤5)、设计跟踪控制器:应用李雅普诺夫稳定性理论设计符合要求的跟踪控制器:
其中,m11,m22和m33是无人艇在u,v,r方向上的质量(包括附加的质量),k1=2和k2=1分别为跟踪控制器τr和τu的控制增益,z1和z2为轨迹跟踪误差经过步骤2)处理后得到的转化误差,χ1和χ2分别为跟踪控制器τr和τu包含的中间变量,ψe为实际航向角与参考航向角的误差,ψe=ψ-ψd,和分别为跟踪控制器τr和τu中用于抵消时变扰动的扰动观测器,yr和yu分别为虚拟控制器αr和αu应用动态面技术而引入的层面误差,yr=αfr-αr,yu=αfu-αu。
扰动观测器和设计如下:
其中,ξ1和ξ2分别为扰动观测器和的状态,ξ1(0)=3和ξ2(0)=-4,kdr=4和kdu=1分别为扰动观测器和的设计参数。
图6为无人艇轨迹跟踪控制的输出轨迹示意图。图7、图8分别为无人艇控制***的跟踪控制器输入τr和τu的示意图。图9、图10分别为无人艇扰动观测器与实际扰动比对的示意图。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。
Claims (7)
1.一种不确定欠驱动无人艇***的预设性能控制方法,其特征在于,包含如下步骤:
步骤1)、建立无人艇动态模型:包括无人艇的运动学和动力学模型;
步骤2)、提出预设性能控制法处理跟踪误差约束条件:定义无人艇实际轨迹与理想轨迹之间的跟踪误差,根据避免跟踪控制器奇异的条件设定跟踪误差需满足的约束条件,然后提出预设性能控制技术把受到约束条件的跟踪误差转化为不受约束条件的误差;
步骤3)、应用动态面控制技术处理虚拟控制器:采用动态面技术避免虚拟控制器的求导;
步骤4)、采用神经网络在线学习***的不确定动态:利用RBF神经网络对任意光滑函数以任意精度逼近的特性,分别构造第一RBF神经网络和第二RBF神经网络对***不确定动态fr(r)和fu(u)进行在线学习,其中,神经网络的权值更新率设计如下:
<mrow>
<msub>
<mover>
<mover>
<mi>W</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&Gamma;</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>r</mi>
<mi>e</mi>
</msub>
<msub>
<mi>S</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<msub>
<mi>Z</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mn>1</mn>
</msub>
<mo>|</mo>
<msub>
<mi>r</mi>
<mi>e</mi>
</msub>
<mo>|</mo>
<msub>
<mover>
<mi>W</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>&ForAll;</mo>
<msub>
<mi>Z</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>r</mi>
<mo>&Element;</mo>
<msub>
<mi>&Omega;</mi>
<mn>1</mn>
</msub>
</mrow>
<mrow>
<msub>
<mover>
<mover>
<mi>W</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&Gamma;</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>u</mi>
<mi>e</mi>
</msub>
<msub>
<mi>S</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<msub>
<mi>Z</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>&sigma;</mi>
<mn>2</mn>
</msub>
<mo>|</mo>
<msub>
<mi>u</mi>
<mi>e</mi>
</msub>
<mo>|</mo>
<msub>
<mover>
<mi>W</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>&ForAll;</mo>
<msub>
<mi>Z</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mi>u</mi>
<mo>&Element;</mo>
<msub>
<mi>&Omega;</mi>
<mn>2</mn>
</msub>
</mrow>
其中,和分别为不确定动态fr(r)和fu(u)的神经网络权值向量;和分别为不确定动态fr(r)和fu(u)的神经网络权值向量的转置向量;Γ1和Γ2分别为和所设计的对称正定矩阵;re和ue分别为r,u与各自滤波虚拟控制器之间的误差,所述滤波虚拟控制器为步骤3)中虚拟控制器通过动态控制技术处理后得到的变量;S1(Z1)和S2(Z2)均为高斯函数;σ1和σ2分别为和设计的正常数;r和u分别是无人艇的航向角速度和纵向角速度;Z1是第一RBF神经网络的输入量,数值等于r,Z2是第二RBF神经网络的输入量,数值等于u;Ω1和Ω2分别表示r和u的紧集;
步骤5)、设计跟踪控制器:应用李雅普诺夫稳定性理论设计符合要求的跟踪控制器:
<mrow>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>=</mo>
<msub>
<mi>m</mi>
<mn>33</mn>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>1</mn>
</msub>
<msub>
<mi>r</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&chi;</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<mn>11</mn>
</msub>
<mo>-</mo>
<msub>
<mi>m</mi>
<mn>22</mn>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mn>33</mn>
</msub>
</mfrac>
<mi>u</mi>
<mi>v</mi>
<mo>-</mo>
<msubsup>
<mover>
<mi>W</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>S</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<msub>
<mi>Z</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>-</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>w</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>-</mo>
<mfrac>
<msub>
<mi>y</mi>
<mi>r</mi>
</msub>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&tau;</mi>
<mi>u</mi>
</msub>
<mo>=</mo>
<msub>
<mi>m</mi>
<mn>11</mn>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mn>2</mn>
</msub>
<msub>
<mi>u</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&chi;</mi>
<mn>2</mn>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>m</mi>
<mn>22</mn>
</msub>
<msub>
<mi>m</mi>
<mn>11</mn>
</msub>
</mfrac>
<mi>v</mi>
<mi>r</mi>
<mo>-</mo>
<msubsup>
<mover>
<mi>W</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>S</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<msub>
<mi>Z</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>-</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>w</mi>
<mi>u</mi>
</mrow>
</msub>
<mo>-</mo>
<mfrac>
<msub>
<mi>y</mi>
<mi>u</mi>
</msub>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中,m11,m22和m33是无人艇在u,v,r各自速度方向上的质量,包括附加的质量,u为无人艇的纵向角速度,v为无人艇的横荡速度,r为无人艇的航向角速度,k1和k2分别为跟踪控制器τr和τu的控制增益,z1和z2为轨迹跟踪误差经过步骤2)处理后得到的转化误差,χ1和χ2分别为跟踪控制器τr和τu包含的中间变量,ψe为无人艇实际航向角与参考轨迹航向角的误差,和分别为跟踪控制器τr和τu中用于抵消时变扰动的扰动观测器,yr和yu分别为虚拟控制器αr和αu采用动态面技术而引入的层面误差,c1和c2分别为虚拟控制器αr和αu应用动态面技术时引入的设计参数。
2.根据权利要求1所述的一种不确定欠驱动无人艇***的预设性能控制方法,其特征在于:步骤1)中,所述的无人艇动态模型为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>u</mi>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mi>&psi;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mi>&psi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>u</mi>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mi>&psi;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>v</mi>
<mi>cos</mi>
<mrow>
<mo>(</mo>
<mi>&psi;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>&psi;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>r</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>u</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfrac>
<msub>
<mi>m</mi>
<mn>22</mn>
</msub>
<msub>
<mi>m</mi>
<mn>11</mn>
</msub>
</mfrac>
<mi>v</mi>
<mi>r</mi>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>u</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>m</mi>
<mn>11</mn>
</msub>
</mfrac>
<msub>
<mi>&tau;</mi>
<mi>u</mi>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>m</mi>
<mn>11</mn>
</msub>
</mfrac>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>w</mi>
<mi>u</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>v</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>m</mi>
<mn>11</mn>
</msub>
<msub>
<mi>m</mi>
<mn>22</mn>
</msub>
</mfrac>
<mi>u</mi>
<mi>r</mi>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>v</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>m</mi>
<mn>22</mn>
</msub>
</mfrac>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>w</mi>
<mi>v</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>r</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<mn>11</mn>
</msub>
<mo>-</mo>
<msub>
<mi>m</mi>
<mn>22</mn>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mn>33</mn>
</msub>
</mfrac>
<mi>u</mi>
<mi>v</mi>
<mo>+</mo>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>m</mi>
<mn>33</mn>
</msub>
</mfrac>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>m</mi>
<mn>33</mn>
</msub>
</mfrac>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>w</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,(x,y)为无人艇的位置,ψ为航向角,u为无人艇的纵向角速度,v为无人艇的横荡速度,r为无人艇的航向角速度,表示无人艇的实际轨迹关于时间t的导数,表示无人艇的各个速度关于时间t的导数,(τu,τr)为无人艇的控制输入,(τwu,τwv,τwr)为外界时变扰动,(fu(u),fv(v),fr(r))为无人艇动力学模型中的不确定动态,具有如下形式:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mi>u</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>d</mi>
<mi>u</mi>
</msub>
<msub>
<mi>m</mi>
<mn>11</mn>
</msub>
</mfrac>
<mi>u</mi>
<mo>-</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
<mn>3</mn>
</munderover>
<mfrac>
<msub>
<mi>d</mi>
<mrow>
<mi>u</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>m</mi>
<mn>11</mn>
</msub>
</mfrac>
<mo>|</mo>
<mi>u</mi>
<msup>
<mo>|</mo>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>u</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mi>v</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>v</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>d</mi>
<mi>v</mi>
</msub>
<msub>
<mi>m</mi>
<mn>22</mn>
</msub>
</mfrac>
<mi>v</mi>
<mo>-</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
<mn>3</mn>
</munderover>
<mfrac>
<msub>
<mi>d</mi>
<mrow>
<mi>v</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>m</mi>
<mn>22</mn>
</msub>
</mfrac>
<mo>|</mo>
<mi>v</mi>
<msup>
<mo>|</mo>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>v</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>d</mi>
<mi>r</mi>
</msub>
<msub>
<mi>m</mi>
<mn>33</mn>
</msub>
</mfrac>
<mi>r</mi>
<mo>-</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
<mn>3</mn>
</munderover>
<mfrac>
<msub>
<mi>d</mi>
<mrow>
<mi>r</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>m</mi>
<mn>33</mn>
</msub>
</mfrac>
<mo>|</mo>
<mi>r</mi>
<msup>
<mo>|</mo>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>r</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
du表示u方向的第一个水动力阻尼系数,dui表示u方向的第i个水动力阻尼系数,dv表示v方向的第一个水动力阻尼系数,dvi表示v方向的第i个水动力阻尼系数,dr表示r方向的第一个水动力阻尼系数,dri表示r方向的第i个水动力阻尼系数。
3.根据权利要求1所述的一种不确定欠驱动无人艇***的预设性能控制方法,其特征在于:步骤2)中,
所述理想轨迹即参考轨迹满足如下条件:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mo>=</mo>
<msub>
<mi>u</mi>
<mi>d</mi>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>d</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mo>=</mo>
<msub>
<mi>u</mi>
<mi>d</mi>
</msub>
<mi>sin</mi>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>d</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,所述参考轨迹是一个圆,(xd,yd)为参考轨迹的位置,是参考轨迹关于时间t的导数,rd表示参考轨迹的航向角速度,ud表示参考轨迹的纵向角速度;
为了避免跟踪控制器奇异,跟踪误差ψe(t)和ze(t)需满足的约束条件如下:
<mrow>
<msub>
<mi>&psi;</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&NotEqual;</mo>
<mo>&PlusMinus;</mo>
<mi>&pi;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>t</mi>
<mo>></mo>
<mn>0</mn>
</mrow>
<mrow>
<msub>
<mi>z</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&NotEqual;</mo>
<mn>0</mn>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>t</mi>
<mo>></mo>
<mn>0</mn>
</mrow>
ψe(t)具体为无人艇实际航向角与参考轨迹航向角的误差,即为ψe,ze(t)具体为无人艇实际位置与参考轨迹位置的误差,即为ze;为了满足上述约束条件,同时能够进一步刻画***的跟踪误差暂态性能及稳态性能要求,将上述的约束条件常数上下界进一步描述成预先指定的时变函数上下界;上述的约束条件刻画成如下更加严格的约束条件:
<mrow>
<mo>-</mo>
<msub>
<munder>
<mi>e</mi>
<mo>&OverBar;</mo>
</munder>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<msub>
<mi>&psi;</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>t</mi>
<mo>></mo>
<mn>0</mn>
</mrow>
<mrow>
<msub>
<munder>
<mi>e</mi>
<mo>&OverBar;</mo>
</munder>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<msub>
<mi>z</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mo>&ForAll;</mo>
<mi>t</mi>
<mo>></mo>
<mn>0</mn>
</mrow>
其中,分别为跟踪误差ψe(t)和ze(t)约束条件的上下界函数,它们是一种按指数衰减的单调递减函数,且满足:
<mrow>
<mn>0</mn>
<mo><</mo>
<msub>
<munder>
<mi>e</mi>
<mo>&OverBar;</mo>
</munder>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<mi>&pi;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
<mo><</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<mi>&pi;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>,</mo>
<mn>0</mn>
<mo><</mo>
<msub>
<munder>
<mi>e</mi>
<mo>&OverBar;</mo>
</munder>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mn>0</mn>
<mo><</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,上下界函数具体为:
<mrow>
<msub>
<munder>
<mi>e</mi>
<mo>&OverBar;</mo>
</munder>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>&pi;</mi>
<mo>/</mo>
<mn>2</mn>
<mo>-</mo>
<mn>0.02</mn>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mn>0.1</mn>
<mi>t</mi>
</mrow>
</msup>
<mo>+</mo>
<mn>0.02</mn>
</mrow>
<mrow>
<msub>
<munder>
<mi>e</mi>
<mo>&OverBar;</mo>
</munder>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>40</mn>
<mo>-</mo>
<mn>0.5</mn>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mn>0.1</mn>
<mi>t</mi>
</mrow>
</msup>
<mo>+</mo>
<mn>0.5</mn>
<mo>;</mo>
</mrow>
引入一个严格单调递增的光滑转化函数R(z),将受到约束的跟踪误差转化为不受约束的误差:
<mrow>
<msub>
<mi>&psi;</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>z</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
设计了如下的转化函数:
<mrow>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>e</mi>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
</msup>
<mo>-</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
</mrow>
</msup>
</mrow>
<mrow>
<msup>
<mi>e</mi>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
</msup>
<mo>+</mo>
<msubsup>
<mi>&gamma;</mi>
<mn>1</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
</mrow>
</msup>
</mrow>
</mfrac>
<mo>,</mo>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msup>
<mi>e</mi>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
</msup>
<mo>+</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
</mrow>
</msup>
</mrow>
<mrow>
<msup>
<mi>e</mi>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
</msup>
<mo>+</mo>
<msubsup>
<mi>&gamma;</mi>
<mn>2</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
</mrow>
</msup>
</mrow>
</mfrac>
</mrow>
其中,γ1为第一个转换函数R1(z1)的中间变量,γ2为第二个转换函数R2(z2)的中间变量,得到如下的转化误差:
<mrow>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mi>e</mi>
</msub>
<mo>/</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<munder>
<mi>e</mi>
<mo>&OverBar;</mo>
</munder>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&psi;</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>&psi;</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msubsup>
<mi>R</mi>
<mn>2</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mi>e</mi>
</msub>
<mo>/</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<munder>
<mi>e</mi>
<mo>&OverBar;</mo>
</munder>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>z</mi>
<mi>e</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
4.根据权利要求1所述一种不确定欠驱动无人艇***的预设性能控制方法,其特征在于:步骤3)中,所述应用动态面控制技术处理虚拟控制器,具体如下:
<mrow>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
<msub>
<mover>
<mi>&alpha;</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>f</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>f</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&alpha;</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>f</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&alpha;</mi>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
<msub>
<mover>
<mi>&alpha;</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>f</mi>
<mi>u</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>f</mi>
<mi>u</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&alpha;</mi>
<mi>u</mi>
</msub>
<mo>,</mo>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>f</mi>
<mi>u</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>&alpha;</mi>
<mi>u</mi>
</msub>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,αfr,αfu分别为r,u的滤波虚拟控制,αr,αu分别为r,u的虚拟控制器,c1,c2分别是αfr,αfu选定的滤波器时间常数;构造虚拟控制器αr和αu如下:
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>r</mi>
</msub>
<mo>=</mo>
<msub>
<mover>
<mi>&psi;</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>&chi;</mi>
<mn>1</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
</msub>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&phi;</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>u</mi>
</msub>
<mo>=</mo>
<msup>
<mi>cos</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&lsqb;</mo>
<mo>-</mo>
<msubsup>
<mi>&chi;</mi>
<mn>2</mn>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>k</mi>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
</msub>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&phi;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>v</mi>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mi>e</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mi>d</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&psi;</mi>
<mi>d</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
其中,ψd为参考轨迹的航向角,ψd=atan2(ye,xe);为参考轨迹关于时间t的导数;和分别为虚拟控制器αr和αu的设计参数,且 φ1和φ2分别为虚拟控制器αr和αu的中间变量:
<mrow>
<msub>
<mi>&phi;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<munder>
<mover>
<mi>e</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>&OverBar;</mo>
</munder>
<mn>1</mn>
</msub>
<mrow>
<msub>
<munder>
<mi>e</mi>
<mo>&OverBar;</mo>
</munder>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&psi;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mover>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mrow>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&psi;</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mover>
<mi>&gamma;</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&phi;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<munder>
<mover>
<mi>e</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>&OverBar;</mo>
</munder>
<mn>2</mn>
</msub>
<mrow>
<msub>
<munder>
<mi>e</mi>
<mo>&OverBar;</mo>
</munder>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>z</mi>
<mi>e</mi>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mover>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mrow>
<msub>
<mi>z</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mover>
<mi>&gamma;</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中,分别表示跟踪误差ψe(t)和ze(t)约束条件的上下界函数关于时间t的导数;分别表示γ1、γ2关于时间t的导数。
5.根据权利要求1所述一种不确定欠驱动无人艇***的预设性能控制方法,其特征在于:步骤4)中,所述Γ1和Γ2分别为和所设计的对称正定矩阵,具体为:Г1=Γ2=1;所述σ1和σ2分别为和设计的正常数,具体为:σ1=σ2=1;所述re和ue分别为r,u与各自滤波虚拟控制器之间的误差,具体为:re=r-αfr,ue=u-αfu,其中αfr,αfu分别为r,u的滤波虚拟控制。
6.根据权利要求1所述一种不确定欠驱动无人艇***的预设性能控制方法,其特征在于:步骤5)中,扰动观测器和设计如下:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>&xi;</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>r</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>d</mi>
<mi>r</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<mn>11</mn>
</msub>
<mo>-</mo>
<msub>
<mi>m</mi>
<mn>22</mn>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mn>33</mn>
</msub>
</mfrac>
<mi>u</mi>
<mi>v</mi>
<mo>+</mo>
<msubsup>
<mover>
<mi>W</mi>
<mo>^</mo>
</mover>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>S</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<msub>
<mi>Z</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>m</mi>
<mn>33</mn>
</msub>
</mfrac>
<msub>
<mi>&tau;</mi>
<mi>r</mi>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>w</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<msub>
<mi>y</mi>
<mi>r</mi>
</msub>
<msub>
<mi>c</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>w</mi>
<mi>r</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&xi;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>d</mi>
<mi>r</mi>
</mrow>
</msub>
<msub>
<mi>r</mi>
<mi>e</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>&xi;</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>u</mi>
<mi>e</mi>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>d</mi>
<mi>u</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>m</mi>
<mn>22</mn>
</msub>
<msub>
<mi>m</mi>
<mn>11</mn>
</msub>
</mfrac>
<mi>v</mi>
<mi>r</mi>
<mo>+</mo>
<msubsup>
<mover>
<mi>W</mi>
<mo>^</mo>
</mover>
<mn>2</mn>
<mi>T</mi>
</msubsup>
<msub>
<mi>S</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<msub>
<mi>Z</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>m</mi>
<mn>11</mn>
</msub>
</mfrac>
<msub>
<mi>&tau;</mi>
<mi>u</mi>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>w</mi>
<mi>u</mi>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<msub>
<mi>y</mi>
<mi>u</mi>
</msub>
<msub>
<mi>c</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>d</mi>
<mo>^</mo>
</mover>
<mrow>
<mi>w</mi>
<mi>u</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>&xi;</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>d</mi>
<mi>u</mi>
</mrow>
</msub>
<msub>
<mi>u</mi>
<mi>e</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,ξ1和ξ2分别为扰动观测器和的状态;kdr和kdu分别为扰动观测器和的设计参数,且kdr>0,kdu>0。
7.根据权利要求1所述一种不确定欠驱动无人艇***的预设性能控制方法,其特征在于:步骤5)中,所述χ1和χ2分别为跟踪控制器τr和τu包含的中间变量,具体为:
<mrow>
<msub>
<mi>&chi;</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>1</mn>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>&gamma;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&chi;</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>&gamma;</mi>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>R</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
其中,即为表示跟踪误差ψe(t)约束条件的上界函数;即为表示跟踪误差ze(t)约束条件的上界函数;所述yr和yu分别为虚拟控制器αr和αu采用动态面技术而引入的层面误差,具体为:yr=αfr-αr,yu=αfu-αu,其中,αfr,αfu分别为r,u的滤波虚拟控制,αr,αu分别为r,u的虚拟控制器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711143413.5A CN108008628B (zh) | 2017-11-17 | 2017-11-17 | 一种不确定欠驱动无人艇***的预设性能控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711143413.5A CN108008628B (zh) | 2017-11-17 | 2017-11-17 | 一种不确定欠驱动无人艇***的预设性能控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108008628A true CN108008628A (zh) | 2018-05-08 |
CN108008628B CN108008628B (zh) | 2020-02-18 |
Family
ID=62052732
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711143413.5A Active CN108008628B (zh) | 2017-11-17 | 2017-11-17 | 一种不确定欠驱动无人艇***的预设性能控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108008628B (zh) |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109143867A (zh) * | 2018-09-26 | 2019-01-04 | 上海海事大学 | 一种基于神经网络控制的混合动力船舶的能量管理方法 |
CN109240091A (zh) * | 2018-11-13 | 2019-01-18 | 燕山大学 | 一种基于强化学习的水下机器人控制方法及其进行跟踪的控制方法 |
CN109240316A (zh) * | 2018-11-15 | 2019-01-18 | 哈尔滨工程大学 | 考虑推进器输出饱和的海底飞行节点预设性能轨迹跟踪控制方法 |
CN109283941A (zh) * | 2018-11-15 | 2019-01-29 | 哈尔滨工程大学 | 基于扰动观测器的预设性能海底飞行节点轨迹跟踪控制方法 |
CN109407671A (zh) * | 2018-12-11 | 2019-03-01 | 大连海事大学 | 一种欠驱动无人船目标包围控制器结构及设计方法 |
CN109633605A (zh) * | 2018-12-29 | 2019-04-16 | 安徽优思天成智能科技有限公司 | 一种船舶废气监测激光雷达的自适应鲁棒随动控制方法 |
CN109739088A (zh) * | 2019-01-07 | 2019-05-10 | 大连海事大学 | 一种无人船有限时间收敛状态观测器及其设计方法 |
CN110262494A (zh) * | 2019-06-26 | 2019-09-20 | 华南理工大学 | 一种同构多无人艇***的协同学习与编队控制方法 |
CN110362075A (zh) * | 2019-06-26 | 2019-10-22 | 华南理工大学 | 一种具有预设性能的无人艇输出反馈编队控制设计方法 |
CN110618611A (zh) * | 2019-10-23 | 2019-12-27 | 哈尔滨工程大学 | 一种基于回转率约束的无人艇轨迹跟踪安全控制方法 |
CN111158383A (zh) * | 2020-01-19 | 2020-05-15 | 浙江大学 | 基于干扰观测器和rbfnn的无人艇航迹跟踪控制方法 |
CN111221335A (zh) * | 2020-01-17 | 2020-06-02 | 大连海事大学 | 一种智能船舶自动舵***的模糊自适应输出反馈有限时间控制方法及*** |
CN111290387A (zh) * | 2020-02-21 | 2020-06-16 | 大连海事大学 | 一种智能船舶自动舵***的模糊自适应输出反馈指定性能控制方法及*** |
CN111880413A (zh) * | 2020-08-12 | 2020-11-03 | 东南大学 | 一种针对船舶航向保持的自适应动态面算法 |
CN112083654A (zh) * | 2020-09-16 | 2020-12-15 | 交通运输部东海航海保障中心连云港航标处 | 一种航标船抗干扰轨迹跟踪控制方法 |
CN113093735A (zh) * | 2021-03-18 | 2021-07-09 | 华南理工大学 | 随机扰动下输出受限欠驱动无人艇的横截函数控制方法、***及介质 |
CN113220000A (zh) * | 2021-05-11 | 2021-08-06 | 华中科技大学 | 一种面向水下探测作业的无人艇路径跟踪预设性能控制方法及*** |
CN113835338A (zh) * | 2021-09-08 | 2021-12-24 | 哈尔滨工程大学 | 基于积分滑模的欠驱动无人艇的有限时间跟踪控制方法及装置 |
CN114326756A (zh) * | 2021-09-08 | 2022-04-12 | 哈尔滨工程大学 | 一种考虑输入量化的自适应预设性能轨迹跟踪控制方法 |
CN114706298A (zh) * | 2021-09-08 | 2022-07-05 | 哈尔滨工程大学 | 一种基于预设性能的usv鲁棒无模型轨迹跟踪控制器设计方法 |
CN116400691A (zh) * | 2023-03-29 | 2023-07-07 | 大连海事大学 | 一种新型离散时间指定性能强化学习无人艇航向跟踪控制方法及*** |
CN117452827A (zh) * | 2023-12-20 | 2024-01-26 | 广东海洋大学 | 一种欠驱动无人艇轨迹跟踪控制方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020771A (zh) * | 2014-06-13 | 2014-09-03 | 大连海事大学 | 一种基于动态虚拟小船制导算法的欠驱动船舶路径跟踪规划方法 |
CN104076821A (zh) * | 2014-06-19 | 2014-10-01 | 哈尔滨工程大学 | 基于模糊自适应观测器的欠驱动水面艇轨迹跟踪控制*** |
CN104698840A (zh) * | 2015-02-02 | 2015-06-10 | 中国民航大学 | 一种抗干扰和参数摄动的四旋翼uav渐近稳定控制方法 |
CN105511480A (zh) * | 2016-02-25 | 2016-04-20 | 哈尔滨工程大学 | 一种基于人工蜂群优化的模型预测的欠驱动uuv深度控制方法 |
CN106292287A (zh) * | 2016-09-20 | 2017-01-04 | 哈尔滨工程大学 | 一种基于自适应滑模控制的uuv路径跟踪方法 |
CN107015562A (zh) * | 2017-04-26 | 2017-08-04 | 华南理工大学 | 一种满足预设跟踪性能的欠驱动水面船舶控制方法 |
CN107065847A (zh) * | 2016-03-31 | 2017-08-18 | 北京航空航天大学 | 一种执行器非对称饱和的水面舰船轨迹跟踪控制方法 |
CN107085427A (zh) * | 2017-05-11 | 2017-08-22 | 华南理工大学 | 一种基于领导跟随结构的无人水面艇编队控制方法 |
CN107121928A (zh) * | 2017-05-22 | 2017-09-01 | 华南理工大学 | 一种无人水面艇的扰动补偿控制方法 |
CN107168312A (zh) * | 2017-05-17 | 2017-09-15 | 哈尔滨工程大学 | 一种补偿uuv运动学和动力学干扰的空间轨迹跟踪控制方法 |
CN107255923A (zh) * | 2017-06-14 | 2017-10-17 | 哈尔滨工程大学 | 基于rbf辨识的ica‑cmac神经网络的欠驱动无人艇航迹跟踪控制方法 |
US20170302364A1 (en) * | 2016-02-04 | 2017-10-19 | Proxy Technologies, Inc. | Unmanned vehicle, system and method for transmitting signals |
CN107315348A (zh) * | 2017-07-04 | 2017-11-03 | 哈尔滨工程大学 | 一种基于惩罚式小波网络的无人潜航器位姿控制装置及方法 |
-
2017
- 2017-11-17 CN CN201711143413.5A patent/CN108008628B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104020771A (zh) * | 2014-06-13 | 2014-09-03 | 大连海事大学 | 一种基于动态虚拟小船制导算法的欠驱动船舶路径跟踪规划方法 |
CN104076821A (zh) * | 2014-06-19 | 2014-10-01 | 哈尔滨工程大学 | 基于模糊自适应观测器的欠驱动水面艇轨迹跟踪控制*** |
CN104698840A (zh) * | 2015-02-02 | 2015-06-10 | 中国民航大学 | 一种抗干扰和参数摄动的四旋翼uav渐近稳定控制方法 |
US20170302364A1 (en) * | 2016-02-04 | 2017-10-19 | Proxy Technologies, Inc. | Unmanned vehicle, system and method for transmitting signals |
CN105511480A (zh) * | 2016-02-25 | 2016-04-20 | 哈尔滨工程大学 | 一种基于人工蜂群优化的模型预测的欠驱动uuv深度控制方法 |
CN107065847A (zh) * | 2016-03-31 | 2017-08-18 | 北京航空航天大学 | 一种执行器非对称饱和的水面舰船轨迹跟踪控制方法 |
CN106292287A (zh) * | 2016-09-20 | 2017-01-04 | 哈尔滨工程大学 | 一种基于自适应滑模控制的uuv路径跟踪方法 |
CN107015562A (zh) * | 2017-04-26 | 2017-08-04 | 华南理工大学 | 一种满足预设跟踪性能的欠驱动水面船舶控制方法 |
CN107085427A (zh) * | 2017-05-11 | 2017-08-22 | 华南理工大学 | 一种基于领导跟随结构的无人水面艇编队控制方法 |
CN107168312A (zh) * | 2017-05-17 | 2017-09-15 | 哈尔滨工程大学 | 一种补偿uuv运动学和动力学干扰的空间轨迹跟踪控制方法 |
CN107121928A (zh) * | 2017-05-22 | 2017-09-01 | 华南理工大学 | 一种无人水面艇的扰动补偿控制方法 |
CN107255923A (zh) * | 2017-06-14 | 2017-10-17 | 哈尔滨工程大学 | 基于rbf辨识的ica‑cmac神经网络的欠驱动无人艇航迹跟踪控制方法 |
CN107315348A (zh) * | 2017-07-04 | 2017-11-03 | 哈尔滨工程大学 | 一种基于惩罚式小波网络的无人潜航器位姿控制装置及方法 |
Non-Patent Citations (4)
Title |
---|
SHI-LU DAI.ET: "《Dynamic Learning From Adaptive Neural Network Control of a Class of Nonaffine Nonlinear Systems》", 《IEEE TRANSACTIONS ON NEURAL NETWORKS AND LEARNING SYSTEMS》 * |
SHI-LU DAI.ET: "《Identification and Learning Control of Ocean Surface Ship Using Neural Networks》", 《IEEE TRANSACTIONS ON INDUSTRIAL INFORMATICS》 * |
SHI-LU DAI.ET: "《Neural Learning Control of Marine Surface Vessels With Guaranteed Transient Tracking Performance》", 《IEEE TRANSACTIONS ON INDUSTRIAL ELECTRONICS》 * |
戴诗陆等: "《一类不确定***的鲁棒可靠跟踪控制及其在飞行控制中的应用》", 《自动化学报》 * |
Cited By (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109143867A (zh) * | 2018-09-26 | 2019-01-04 | 上海海事大学 | 一种基于神经网络控制的混合动力船舶的能量管理方法 |
CN109240091A (zh) * | 2018-11-13 | 2019-01-18 | 燕山大学 | 一种基于强化学习的水下机器人控制方法及其进行跟踪的控制方法 |
CN109240091B (zh) * | 2018-11-13 | 2020-08-11 | 燕山大学 | 一种基于强化学习的水下机器人控制方法及其进行跟踪的控制方法 |
CN109283941B (zh) * | 2018-11-15 | 2021-08-13 | 哈尔滨工程大学 | 基于扰动观测器的预设性能海底飞行节点轨迹跟踪控制方法 |
CN109240316A (zh) * | 2018-11-15 | 2019-01-18 | 哈尔滨工程大学 | 考虑推进器输出饱和的海底飞行节点预设性能轨迹跟踪控制方法 |
CN109283941A (zh) * | 2018-11-15 | 2019-01-29 | 哈尔滨工程大学 | 基于扰动观测器的预设性能海底飞行节点轨迹跟踪控制方法 |
CN109407671A (zh) * | 2018-12-11 | 2019-03-01 | 大连海事大学 | 一种欠驱动无人船目标包围控制器结构及设计方法 |
CN109407671B (zh) * | 2018-12-11 | 2021-06-29 | 大连海事大学 | 一种欠驱动无人船目标包围控制器结构及设计方法 |
CN109633605A (zh) * | 2018-12-29 | 2019-04-16 | 安徽优思天成智能科技有限公司 | 一种船舶废气监测激光雷达的自适应鲁棒随动控制方法 |
CN109739088B (zh) * | 2019-01-07 | 2021-08-27 | 大连海事大学 | 一种无人船有限时间收敛状态观测器及其设计方法 |
CN109739088A (zh) * | 2019-01-07 | 2019-05-10 | 大连海事大学 | 一种无人船有限时间收敛状态观测器及其设计方法 |
CN110362075A (zh) * | 2019-06-26 | 2019-10-22 | 华南理工大学 | 一种具有预设性能的无人艇输出反馈编队控制设计方法 |
CN110262494A (zh) * | 2019-06-26 | 2019-09-20 | 华南理工大学 | 一种同构多无人艇***的协同学习与编队控制方法 |
CN110362075B (zh) * | 2019-06-26 | 2020-09-22 | 华南理工大学 | 一种具有预设性能的无人艇输出反馈编队控制设计方法 |
CN110618611A (zh) * | 2019-10-23 | 2019-12-27 | 哈尔滨工程大学 | 一种基于回转率约束的无人艇轨迹跟踪安全控制方法 |
CN111221335A (zh) * | 2020-01-17 | 2020-06-02 | 大连海事大学 | 一种智能船舶自动舵***的模糊自适应输出反馈有限时间控制方法及*** |
CN111158383A (zh) * | 2020-01-19 | 2020-05-15 | 浙江大学 | 基于干扰观测器和rbfnn的无人艇航迹跟踪控制方法 |
CN111158383B (zh) * | 2020-01-19 | 2021-04-20 | 浙江大学 | 基于干扰观测器和rbfnn的无人艇航迹跟踪控制方法 |
CN111290387A (zh) * | 2020-02-21 | 2020-06-16 | 大连海事大学 | 一种智能船舶自动舵***的模糊自适应输出反馈指定性能控制方法及*** |
CN111880413A (zh) * | 2020-08-12 | 2020-11-03 | 东南大学 | 一种针对船舶航向保持的自适应动态面算法 |
CN112083654A (zh) * | 2020-09-16 | 2020-12-15 | 交通运输部东海航海保障中心连云港航标处 | 一种航标船抗干扰轨迹跟踪控制方法 |
CN112083654B (zh) * | 2020-09-16 | 2022-09-16 | 交通运输部东海航海保障中心连云港航标处 | 一种航标船抗干扰轨迹跟踪控制方法 |
CN113093735B (zh) * | 2021-03-18 | 2022-07-26 | 华南理工大学 | 随机扰动下输出受限欠驱动无人艇的横截函数控制方法、***及介质 |
CN113093735A (zh) * | 2021-03-18 | 2021-07-09 | 华南理工大学 | 随机扰动下输出受限欠驱动无人艇的横截函数控制方法、***及介质 |
CN113220000B (zh) * | 2021-05-11 | 2022-12-02 | 华中科技大学 | 一种面向水下探测作业的无人艇路径跟踪预设性能控制方法及*** |
CN113220000A (zh) * | 2021-05-11 | 2021-08-06 | 华中科技大学 | 一种面向水下探测作业的无人艇路径跟踪预设性能控制方法及*** |
CN114326756A (zh) * | 2021-09-08 | 2022-04-12 | 哈尔滨工程大学 | 一种考虑输入量化的自适应预设性能轨迹跟踪控制方法 |
CN114706298A (zh) * | 2021-09-08 | 2022-07-05 | 哈尔滨工程大学 | 一种基于预设性能的usv鲁棒无模型轨迹跟踪控制器设计方法 |
CN113835338A (zh) * | 2021-09-08 | 2021-12-24 | 哈尔滨工程大学 | 基于积分滑模的欠驱动无人艇的有限时间跟踪控制方法及装置 |
CN114706298B (zh) * | 2021-09-08 | 2023-08-04 | 哈尔滨工程大学 | 一种基于预设性能的usv鲁棒无模型轨迹跟踪控制器设计方法 |
CN113835338B (zh) * | 2021-09-08 | 2024-03-01 | 哈尔滨工程大学 | 积分滑模的欠驱动无人艇的有限时间跟踪控制方法及装置 |
CN116400691A (zh) * | 2023-03-29 | 2023-07-07 | 大连海事大学 | 一种新型离散时间指定性能强化学习无人艇航向跟踪控制方法及*** |
CN116400691B (zh) * | 2023-03-29 | 2023-11-21 | 大连海事大学 | 一种新型离散时间指定性能强化学习无人艇航向跟踪控制方法及*** |
CN117452827A (zh) * | 2023-12-20 | 2024-01-26 | 广东海洋大学 | 一种欠驱动无人艇轨迹跟踪控制方法 |
CN117452827B (zh) * | 2023-12-20 | 2024-04-05 | 广东海洋大学 | 一种欠驱动无人艇轨迹跟踪控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108008628B (zh) | 2020-02-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108008628A (zh) | 一种不确定欠驱动无人艇***的预设性能控制方法 | |
Peng et al. | Adaptive dynamic surface control for formations of autonomous surface vehicles with uncertain dynamics | |
Xiang et al. | Robust fuzzy 3D path following for autonomous underwater vehicle subject to uncertainties | |
CN105138001B (zh) | 一种四旋翼飞行器姿态控制方法 | |
Zhang et al. | Constrained coordinated path-following control for underactuated surface vessels with the disturbance rejection mechanism | |
CN107085427A (zh) | 一种基于领导跟随结构的无人水面艇编队控制方法 | |
Chen et al. | A reconfiguration scheme for quadrotor helicopter via simple adaptive control and quantum logic | |
Yuan et al. | Sliding mode observer-based heading control for a gliding robotic dolphin | |
CN104850122B (zh) | 基于可变船长比的抵抗侧风无人水面艇直线路径跟踪方法 | |
Liu et al. | Coordinated path following of multiple underacutated marine surface vehicles along one curve | |
CN106292287A (zh) | 一种基于自适应滑模控制的uuv路径跟踪方法 | |
CN110262494B (zh) | 一种同构多无人艇***的协同学习与编队控制方法 | |
CN105676641A (zh) | 针对四旋翼无人机非线性模型的基于反步和滑模控制技术的非线性鲁棒控制器的设计方法 | |
CN105785762A (zh) | 一种基于自适应反步滑模的两轴惯性稳定平台高精度控制方法 | |
CN104238357A (zh) | 一种近空间飞行器的容错滑模控制方法 | |
CN104898688A (zh) | Uuv四自由度动力定位自适应抗扰滑模控制***及控制方法 | |
CN110377036A (zh) | 一种基于指令约束的无人水面艇航迹跟踪固定时间控制方法 | |
Wu et al. | Global tracking controller for underactuated ship via switching design | |
Zhang et al. | Disturbance observer-based composite neural learning path following control of underactuated ships subject to input saturation | |
CN108170151A (zh) | 一种水下机器人的自适应运动控制装置及其方法 | |
CN103901775B (zh) | 一种基于t-s模型带有输入约束的舵减横摇模糊控制器及其控制方法 | |
Selfridge et al. | Multivariable output feedback MRAC for a quadrotor UAV | |
Zemalache et al. | Controlling a drone: Comparison between a based model method and a fuzzy inference system | |
Gao et al. | Adaptive integral backstepping control for a 3-DOF helicopter | |
Lin et al. | Finite-time formation control of unmanned vehicles using nonlinear sliding mode control with disturbances |
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 |