一种基于非线性模型预测控制的UUV对运动母船的跟踪方法技术

技术编号:13776471 阅读:297 留言:0更新日期:2016-09-30 23:50
一种基于非线性模型预测控制的UUV对运动母船的跟踪方法,涉及欠驱动UUV对运动母船的跟踪控制方法。为了解决现有的针对线性模型进行处理的方法对UUV进行控制时存在准确性不高的问题,本发明专利技术首先建立欠驱动UUV水平面的预测模型和相应的约束条件;然后选取综合性能指标,将欠驱动UUV对运动母船的跟踪控制问题转化为约束条件下的优化问题;并利用泰勒级数展开和李导数对步骤三中约束条件下的优化问题进行处理,并将求得的解析解作为UUV跟踪系统的控制输入;在每一预测时域,用新的状态值刷新步骤四中的解析解,以此不断运行,直至UUV完成对运动母船跟踪作业。本发明专利技术适用于欠驱动UUV对运动母船的跟踪控制。

【技术实现步骤摘要】

本专利技术涉及欠驱动UUV对运动母船的跟踪控制方法。
技术介绍
无人水下航行器(Unmanned Underwater Vehicle,UUV)具有活动范围大、潜水深、机动性好、智能化、运行和维护费用低等优点,作为人类在海洋活动中,特别是深海活动中的重要替代者和执行者,已被广泛应用于科学考察,深海作业等领域,在作业完成时,需要对其进行回收,尤其在采取运动母船背驼搭载回收过程中,需要UUV对运动母船进行准确跟踪。目前,国内外学者针对UUV跟踪控制问题做了大量的研究,相应提出了很多控制方法,如滑模控制,反步控制等等。滑模控制器不仅具有的对系统模型参数不确定及外界扰动的鲁棒性,而且极大地提高了系统状态在平衡点附近的收敛速度。反步法计算简单,实时性好,响应快,对传感器要求低,在工程应用中得到了广泛的应用,但该方法中存在虚拟控制量的高阶导数问题。大多都未能够充分考虑UUV在约束条件下的跟踪控制问题,且在快速性上有一定的局限,而模型预测控制(Model predictive control,MPC)可以很好地解决此类问题,因此,提出将非线性模型预测控制应用到UUV对运动母船跟踪中具有重要的意义。2014年5月Proceedings of the Institution of Mechanical Engineers Part M:Journal of Engineering for the Maritime Environment第228卷中的“Model predictive control of a hybrid autonomous underwater vehicle with experimental verification”提出将模型预测控制用于AUV悬停状态时深度控制,但该方法只是针对线性模型进行处理,若将该方法迁移到UUV对运动母船跟踪控制中,在对模型进行线性化时将影响控制的准确性,同时该方法为在线计算,将影响控制的快速性。
技术实现思路
本专利技术是为了解决现有的针对线性模型进行处理的方法对UUV进行控制时存在准确性不高的问题。一种基于非线性模型预测控制的UUV对运动母船的跟踪方法,包括如下步骤:步骤一:建立欠驱动UUV水平面的预测模型;步骤二:根据实际,建立UUV水平面的预测模型相应的约束条件;步骤三:选取综合性能指标,将欠驱动UUV对运动母船的跟踪控制问题转化为约束 条件下的优化问题;步骤四:利用泰勒级数展开和李导数对步骤三中约束条件下的优化问题进行处理,并将求得的解析解作为UUV跟踪系统的控制输入;步骤五:在每一预测时域,用新的状态值刷新步骤四中的解析解,以此不断运行,直至UUV完成对运动母船跟踪作业。本专利技术具有以下效果:本专利技术提供了一种基于非线性模型预测控制的UUV对运动母船的跟踪方法,该方法将UUV跟踪运动母船控制问题转换为一定约束条件下的优化问题,并利用泰勒级数展开和李导数对优化问题进行处理,最终求得每一预测时域中控制输入的解析解,可以很好的解决针对UUV非线性系统对运动母船的快速跟踪控制。该方法获得的解析解具备离线优化的特点,不仅可以减小海流不确定性因素的累积作用,同时可以增强跟踪控制系统的快速性。附图说明图1为UUV的水平面建模图;图2为基于非线性模型预测控制的UUV对运动母船跟踪控制框图;图3为基于非线性模型预测控制的UUV对运动母船跟踪控制流程图;图4为UUV跟踪正弦运动母船时轨迹;图5为UUV跟踪正弦运动母船时纵向速度和横向速度;其中,图5a为UUV跟踪正弦运动母船时纵向速度曲线图,图5b为UUV跟踪正弦运动母船时横向速度曲线图;图6为UUV跟踪正弦运动母船时艏向角及其加速度;其中,图6a为UUV艏向角的加速度曲线图,图6b为UUV的艏向曲线图;图7为UUV跟踪正弦运动母船时纵向推力和转艏力矩;其中,图7a为UUV执行机构的推力曲线图,图7b为UUV执行机构的转艏力矩曲线图;图8为UUV跟踪正弦运动母船时轨迹误差;其中,图8a为UUV在跟踪母船时x轴的坐标误差曲线图,图8b为UUV在跟踪母船时y轴坐标误差曲线图;图9为UUV跟踪圆运动母船时轨迹;图10为UUV跟踪圆运动母船时纵向速度和横向速度;其中,图10a为UUV跟踪圆运动母船时纵向速度曲线图,图10b为UUV跟踪圆运动母船时横向速度曲线图;图11为UUV跟踪圆运动母船时艏向角及其加速度;其中,图11a为UUV跟踪圆运动母船时加速度曲线图,图11b为UUV跟踪圆运动母船时艏向角曲线图;图12为UUV跟踪圆运动母船时纵向推力和转艏力矩;其中,图12a为UUV跟踪圆运动母船时纵向推力曲线图,图12b为UUV跟踪圆运动母船时转艏力矩曲线图;图13为UUV跟踪圆运动的母船轨迹误差;其中,图13a为UUV跟踪圆运动的母船x轴的坐标误差曲线图,图13b为UUV跟踪圆运动的母船y轴坐标误差曲线图。具体实施方式具体实施方式一:结合图1至图3说明本实施方式,一种基于非线性模型预测控制的UUV对运动母船的跟踪方法,包括如下步骤:步骤一:建立欠驱动UUV水平面的预测模型;步骤二:根据实际,建立UUV水平面的预测模型相应的约束条件;步骤三:选取综合性能指标,将欠驱动UUV对运动母船的跟踪控制问题转化为约束条件下的优化问题;步骤四:利用泰勒级数展开和李导数对步骤三中约束条件下的优化问题进行处理,并将求得的解析解作为UUV跟踪系统的控制输入;步骤五:在每一预测时域,用新的状态值刷新步骤四中的解析解,以此不断运行,直至UUV完成对运动母船跟踪作业。具体实施方式二:本实施方式步骤一所述的建立欠驱动UUV水平面的预测模型的具体过程如下:结合图1,建立了UUV的坐标系统,E-ξηζ是固定坐标系,O-XYZ是运动坐标系;固定坐标系E-ξηζ的原点任取大地或海洋上一点,Eξ保持水平并指向地理北向,被选为UUV的主航向;Eη指向地理东向,代表UUV的横向;Eζ指向地心,代表UUV的垂向;运动坐标系O-XYZ置于UUV上,并随着其一起运动,其原点O置于UUV重心的位置;Ox指向UUV前行的方向为正,绕Ox旋转产生横摇角速度p,形成横摇角φ;Oy指向UUV的右舷方向为正,绕Oy旋转产生纵倾角速度q,形成纵倾角θ;Oz指向UUV底部为正,绕Oz旋转产生艏向角角速度r,形成艏向角ψ,根据该坐标系统建立UUV水平面系统;假设UUV轨迹(x,y,ψ)和母船轨迹(xd,yd,ψd)都是连续的且满足可以作多次微分运算;UUV轨迹中的x、y为UUV在固定坐标系中的位置,ψ为UUV的艏向角,母船轨 迹中的xd、yd为母船在固定坐标系中的位置,ψd为母船在固定坐标系中的艏向角;定义x=[u v r x y ψ]T为状态变量,定义y为UUV轨迹;其中,UUV的状态包括:x、y、ψ、纵向速度u、横向速度v、艏向角角速度r;粗体x为状态变量,细体x为UUV在固定坐标系中的位置;粗体y为UUV轨迹,细体y为UUV在固定坐标系中的位置;进一步令x=[u v r x y ψ]T=[x1 x2 x3 x4 x5 x6]T,得到UUV水平面预测模型:其中,x∈R6,u∈R2,y∈R3,u为控制输入,由u1、u2共同构成,且u1=Xprop,u2=τr;粗体u为控制输入,细体u为UUV纵向速度;本文档来自技高网
...

【技术保护点】
一种基于非线性模型预测控制的UUV对运动母船的跟踪方法,其特征在于包括如下步骤:步骤一:建立欠驱动UUV水平面的预测模型;步骤二:建立UUV水平面的预测模型相应的约束条件;步骤三:选取综合性能指标,将欠驱动UUV对运动母船的跟踪控制问题转化为约束条件下的优化问题;步骤四:利用泰勒级数展开和李导数对步骤三中约束条件下的优化问题进行处理,并将求得的解析解作为UUV跟踪系统的控制输入;步骤五:在每一预测时域,用新的状态值刷新步骤四中的解析解,以此不断运行,直至UUV完成对运动母船跟踪作业。

【技术特征摘要】
1.一种基于非线性模型预测控制的UUV对运动母船的跟踪方法,其特征在于包括如下步骤:步骤一:建立欠驱动UUV水平面的预测模型;步骤二:建立UUV水平面的预测模型相应的约束条件;步骤三:选取综合性能指标,将欠驱动UUV对运动母船的跟踪控制问题转化为约束条件下的优化问题;步骤四:利用泰勒级数展开和李导数对步骤三中约束条件下的优化问题进行处理,并将求得的解析解作为UUV跟踪系统的控制输入;步骤五:在每一预测时域,用新的状态值刷新步骤四中的解析解,以此不断运行,直至UUV完成对运动母船跟踪作业。2.根据权利要求1所述的一种基于非线性模型预测控制的UUV对运动母船的跟踪方法,其特征在于步骤一所述的建立欠驱动UUV水平面的预测模型的具体过程如下:假设UUV轨迹(x,y,ψ)和母船轨迹(xd,yd,ψd)都是连续的且满足作多次微分运算;UUV轨迹中的x、y为UUV在固定坐标系中的位置,ψ为UUV的艏向角,母船轨迹中的xd、yd为母船在固定坐标系中的位置,ψd为母船在固定坐标系中的艏向角;定义x=[u v r x y ψ]T为状态变量,定义y为UUV轨迹;其中,UUV的状态包括:x、y、ψ、纵向速度u、横向速度v、艏向角角速度r;粗体x为状态变量,细体x为UUV在固定坐标系中的位置;粗体y为UUV轨迹,细体y为UUV在固定坐标系中的位置;进一步令x=[u v r x y ψ]T=[x1 x2 x3 x4 x5 x6]T,得到UUV水平面预测模型: x · = f ( x ) + g ( x ) u y = h ( x ) - - - ( 6 ) ]]>其中,x∈R6,u∈R2,y∈R3,u为控制输入,由u1、u2共同构成,且u1=Xprop,u2=τr;粗体u为控制输入,细体u为UUV纵向速度;同时给出式(6)中每项的具体形式: f ( x ) = a 1 x 2 x 3 + a 2 x 1 b 1 x 1 x 3 + b 2 x 2 c 1 x 1 x 2 + c 2 x 3 x 1 cos x 6 - x 2 sin x 6 x 1 sin x 6 + x 2 cos x 6 x 3 , g ( x ) = g 1 ( x ) g 2 ( x ) = a 3 0 0 0 0 0 0 0 c 3 0 0 0 T ]]>u=[u1 u2]T=[Xprop τr]T,y=[h1 h2 h3]T=[x4 x5 x6]T a 1 = m 2 m 1 , a 2 = - d 1 m 1 , a 3 = 1 m 1 , b 1 = - m 1 m 2 , b 2 = - d 2 m 2 , c 1 = m 1 - m 2 m 3 , c 2 = - d 3 m 3 , c 3 = 1 m 3 ]]>其中,Xprop,τr分别为主推进器推力和转艏力矩;d1为前进方向的水动力粘性项,d2横移方向的水动力粘性项,d3转艏方向的水动力粘性系数,m1为前进方向水动力附加质量,m2横移方向的水动力附加质量,m3为转艏方向的水动力附加质量。3.根据权利要求2所述的一种基于非线性模型预测控制的UUV对运动母船的跟踪方法,其特征在于步骤二所述的UUV水平面的预测模型相应的约束条件如下:对于UUV水平面预测模型,在某个预测时域[t,t+T]中,将满足约束: x ^ · ( t + τ ) = f ( x ^ ( t + τ ) ) + g ( x ^ ( t + τ ) ) u ^ ( t + τ ) y ^ ( t + τ ) = h ( x ^ ( t + τ ) ) , τ ∈ [ 0 , T ] - - - ( 7 ) ]]>其中τ为时间变量,且τ∈[0,T];为UUV在时域[t,t+T]中状态变量的预测值;为UUV在时域[t,t+T]中轨迹的预测值;为UUV在时域[t,t+T]中控制输入;同时,在[t,t+T]预测时域中,状态变量预测值的初始值满足: x ^ ( t + τ ) = x ( t ) , τ = 0 - - - ( 8 ) ]]>其中,t表示时刻,T表示预测时域的时间。4.根据权利要求3所述的一种基于非线性模型预测控制的UUV对运动母船的跟踪方法,其特征在于步骤三具体过程如下:在每个预测时域中,选取综合性能指标为: J = 1 2 μ 1 [ y ^ ( t + T ) - y ^ d ( t + T ) ] T [ y ^ ( t + T ) - y ^ d ( t + T ) ] + 1 2 ∫ 0 T [ μ 2 ( y ^ ( t + τ ) - y ^ d ( t + τ ) ) T ( y ^ ( t + τ ) - y ^ d ( t + τ ) ) + μ 3 u ^ T ( t + τ ) u ^ ( t + τ ) ] d τ - - - ( 9 ) ]]>其中,为UUV在终端时刻t+T轨迹预测值,为母船在终端时刻t+T的轨迹,为母船在时域[t,t+T]中的轨迹;通过最小化式(9)所给出的性能指标J,获得跟踪控制系统的最优控制输入并根据在步骤二中UUV水平面预测模型需要满足的约束,将UUV对运动母船的跟踪问题转化为优化问题: m i n u ^ ( t ) J - - - ( 10 ) . ]]>s.t.式(7)和(8)约束5.根据权利要求4所述的一种基于非线性模型预测控制的UUV对运动母船的跟踪方法,其特征在于步骤四具体过程如下:假设控制输入在预测时域[t,t+T]内为一个常值: u ^ ( t + τ ) = u ( t ) = c o n s t , τ ∈ [ 0 , T ] - - - ( 11 ) ]]>其中,const表示常数;在此条件下,得到预测时域[t,t+T]中跟踪控制系统控制输入的各阶导数皆为零;利用泰勒级数展开和李导数对优化问题进行处理,UUV轨迹的预测输出以及母船轨迹展开为: y ^ ( t + τ ) = τ T ( τ ) y ^ ( t ) y ^ d ( t + τ ) = τ T ( τ ) y ^ d ( t ) - - - ( 12 ) ]]>粗体τ为以时间变量τ为自变量的函数;式中,所有带角标d的参数均表示母船对应参数; y ^ ( t + τ ) = h ^ 1 ( t + τ ) h ^ 2 ( t + τ ) h ^ 3 ( t + τ ) = x ^ ( t + τ ) y ^ ( t + τ ) ψ ^ ( t + τ ) , y ^ d ( t + τ ) = h ^ 1 d ( t + τ ) h ^ 2 d ( t + τ ) h ^ 3 d ( t + τ ) = x ^ d ( t + τ ) y ^ d ( t + τ ) ψ ^ d ( t + τ ) ]]> τ ( τ ) = τ 1 τ 2 τ 3 , τ 1 = τ 2 = τ 3 = 1 τ ... τ N N ! T ]]> y ^ ( t ) = h ^ 1 ( t ) h ^ · 1 ( t ) ... h ^ 1 [ N ] ( t ) ...

【专利技术属性】
技术研发人员:张伟郭毅滕延斌梁志成孟德涛周佳加严浙平
申请(专利权)人:哈尔滨工程大学
类型:发明
国别省市:黑龙江;23

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1