CN115309058A - 一种动力定位船的有限时间复合学习控制方法 - Google Patents

一种动力定位船的有限时间复合学习控制方法 Download PDF

Info

Publication number
CN115309058A
CN115309058A CN202211138901.8A CN202211138901A CN115309058A CN 115309058 A CN115309058 A CN 115309058A CN 202211138901 A CN202211138901 A CN 202211138901A CN 115309058 A CN115309058 A CN 115309058A
Authority
CN
China
Prior art keywords
matrix
ship
dynamic positioning
finite time
coordinate system
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
Application number
CN202211138901.8A
Other languages
English (en)
Other versions
CN115309058B (zh
Inventor
李鸣阳
徐骄
谢文博
王海滨
陈海城
刘洋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Publication of CN115309058A publication Critical patent/CN115309058A/zh
Application granted granted Critical
Publication of CN115309058B publication Critical patent/CN115309058B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive 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/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)
  • Navigation (AREA)

Abstract

一种动力定位船的有限时间复合学习控制方法,涉及船舶运动控制技术领域。本发明是为了解决现有动力定位船运动时系统参数不确定以及船舶进行特殊作业时运动输出轨迹受限的问题。本发明采用复合学习控制方法在有限激励条件下能在线辨识系统未知参数,保证各参数向量分量的收敛速度不受激励水平的影响且可灵活调整。采用非对称积分障碍李雅普诺夫函数直接将期望行为指标强加于输出轨迹上,避免了系统保守性增加和计算量增大等问题。本发明基于有限时间理论设计复合学习控制器,进一步提高了控制效果。

Description

一种动力定位船的有限时间复合学习控制方法
技术领域
本发明属于船舶运动控制技术领域。
背景技术
随着不断向海洋领域的拓展和延伸,海洋事业的发展逐渐得到重视,合理开发、利用并切实保护和管理海洋已成为目前人类生存、发展和强盛的重大战略问题。在此背景下,各海洋工程装备领域都得到了飞速发展,具有定位精度高、灵活性好和定位不受海洋深度限制等优点的动力定位船在海洋工程作业中得到了广泛的应用。
动力定位船本身体积大、重量大,同时在其作业过程中受负载变化、燃料消耗和恶劣海洋环境等因素影响,就会导致其系统惯量矩阵和水动力参数矩阵往往是不确定的。但现有方法大多需假设系统参数矩阵在船舶整个作业过程中精确已知,因此增强了所设计控制器的保守性。此外,现有大多数动力定位船运动控制方法只关注于保证系统的输出误差能够收敛到一个残差集内,即保证系统的稳态性能,却忽略了瞬态性能对系统的影响。然而,当动力定位船通过一些狭窄航道、完成一些绕平台作业或其他某些特殊作业任务时,为确保作业任务的顺利完成并保证船舶的安全性,通常会对船舶整个作业过程中的输出轨迹施加一定的约束限制。
发明内容
本发明是为了解决现有动力定位船运动时系统参数不确定以及船舶进行特殊作业时运动输出轨迹受限的问题,现提供一种动力定位船的有限时间复合学习控制方法。
一种动力定位船的有限时间复合学习控制方法,包括以下步骤:
步骤一:参照动力定位船的期望位置ηd对其实际位置η施加约束条件:
Figure BDA0003852579720000011
并建立动力定位船控制系统的误差方程:
Figure BDA0003852579720000012
其中,z1为η和ηd之间的误差,
Figure BDA0003852579720000013
为η的一阶导数,z2
Figure BDA0003852579720000014
与虚拟控制率α之间的误差;-δ
Figure BDA0003852579720000015
分别为实际位置η的下限和上限,η=[η123]T,ηd=[ηd1d2d3]T,α=[α123]T
Figure BDA0003852579720000021
Figure BDA0003852579720000022
ai为误差增益系数,Λi为滑模项增益系数,z1i=ηidi,η1和η2分别为北东坐标系下船舶实际位置的横纵坐标,ηd1和ηd2分别为北东坐标系下船舶期望位置的横纵坐标,η3和ηd3分别为船舶实际与期望的姿态角;
步骤二:基于实际位置η的约束条件构造非对称积分障碍李雅普诺夫函数V1和李雅普诺夫函数V2
Figure BDA0003852579720000023
其中,
Figure BDA0003852579720000024
M*=R(η3)MRT3),R(η3)为北东坐标系和船体坐标系之间的转换矩阵,RT3)=R-13),M为动力定位船控制系统的惯性矩阵,
Figure BDA0003852579720000025
τ为积分变量;
步骤三:对V2求一阶导数:
Figure BDA0003852579720000026
其中,f为动力定位船控制系统的不确定部分,
Figure BDA0003852579720000027
Figure BDA0003852579720000028
D*=R(η3)DRT3),D为动力定位船控制系统的阻尼矩阵,
Figure BDA0003852579720000029
Figure BDA00038525797200000210
Figure BDA00038525797200000211
分别为α和R(η3)的一阶导数,U为动力定位船的控制输入向量;
步骤四:对f进行线性化处理得:f=Φθ,
其中,Φ为初始回归矩阵,θ为未知参数向量,
Figure BDA0003852579720000031
θ=[M11 M22 M23 M33 D11 D22 D23 D32 D33]T
步骤五:在初始回归矩阵Φ的基础上构建增广回归矩阵W:
Figure BDA0003852579720000032
其中,z2f为z2的滤波值,
Figure BDA0003852579720000033
z2f通过如下滤波器获得:
Figure BDA0003852579720000034
c为滤波器时间常数;
步骤六:分别建立增广回归矩阵W和动力定位船的控制输入向量U的稳定滤波器
Figure BDA0003852579720000035
Figure BDA0003852579720000036
Figure BDA0003852579720000037
Figure BDA0003852579720000038
其中,Uf和Wf分别为U和W的滤波值;
步骤七:将步骤六中的两个滤波器代入
Figure BDA0003852579720000039
的一阶导数中求解获得:
Figure BDA00038525797200000310
对Uf进行参数线性化获得:
Uf=Waθ,
其中,Wa为线性化后的回归矩阵;
步骤八:对Uf=Waθ等号左右两侧分别乘以
Figure BDA00038525797200000311
获得:
Figure BDA00038525797200000312
利用Laplace算子
Figure BDA00038525797200000313
建立信息矩阵
Figure BDA00038525797200000314
和辅助矩阵
Figure BDA00038525797200000315
步骤九:基于并行学习思想和动态回归扩张与混合过程设计有限时间复合学习律:
Figure BDA0003852579720000041
其中,γ为正定自适应增益矩阵,Γ为正定调谐参数矩阵,te为信息矩阵的更新截止时刻,Δe为te时刻信息矩阵的行列式,Ψe为te时刻的辅助矩阵,
Figure BDA0003852579720000042
为θ的估计值,
Figure BDA0003852579720000043
Figure BDA0003852579720000044
的一阶导;
步骤十:利用有限时间复合学习律并基于非对称积分障碍李雅普诺夫函数设计有限时间复合学习控制器:
Figure BDA0003852579720000045
其中,K为误差增益矩阵,H为滑模项增益矩阵,
利用有限时间复合学习控制器输出的控制输入向量U控制动力定位船。
进一步的,上述步骤七中求解Uf的具体方法为:
对z2求一阶导数并在其一阶导数的等号两边同时乘以M*获得:
Figure BDA0003852579720000046
建立北东坐标系下动力定位船的数学模型:
Figure BDA0003852579720000047
利用回归滤波方法建立滤波系统动力学方程并执行仿射参数线性化,获得:
Figure BDA0003852579720000048
Figure BDA0003852579720000049
公式二和公式三均代入公式一中获得:
Figure BDA00038525797200000410
将公式四以及稳定滤波器
Figure BDA00038525797200000411
Figure BDA00038525797200000412
代入
Figure BDA00038525797200000413
的一阶导数中获得:
Figure BDA00038525797200000414
Figure BDA00038525797200000415
并对
Figure BDA00038525797200000416
求一阶导,将公式五代入
Figure BDA00038525797200000417
中,对其求解获得:
Figure BDA00038525797200000418
进一步的,上述步骤八中,设计有限时间复合学习律的具体方法为:
分别对信息矩阵Ω和辅助矩阵N积分,获得:
Figure BDA0003852579720000051
Figure BDA0003852579720000052
其中,l为Laplace算子的滤波时间常数,t为积分时间上限,χ为积分变量,
结合公式六和七获得:
N(t)=Ω(t)θ 公式八,
在公式八等号的两边均乘以adj(Ω)获得:
Ψ(t)=Δ(t)θ,
其中,Ψ(t)=N(t)adj(Ω),Δ(t)=Ω(t)adj(Ω),adj(Ω)为信息矩阵Ω的伴随矩阵,
Figure BDA0003852579720000053
其中,
Figure BDA0003852579720000054
为信息矩阵Ω的积分上限,m为待确定向量中未知参数的个数,
将te代入Ψ(t)和Δ(t)分别获得:
Δe=Δ(te),Ψe=Ψ(te),
利用te、Δe和Ψe建立有限时间复合学习律。
本发明所述的一种动力定位船的有限时间复合学习控制方法,考虑了存在系统输出约束和参数不确定性的动力定位船有限时间复合学习控制方法与现有技术相比的优点有:
(1)本发明考虑了船舶作业过程中普遍存在的系统参数不确定这一问题,基于回归滤波、并行学习思想和动态回归扩张与混合过程设计复合学习律,在有限激励条件下确保参数估计收敛于其真值,放宽了参数估计对持续激励条件的依赖。此外,得益于动态回归扩张与混合过程的设计步骤,所有参数向量分量的收敛速度可以灵活调整,且不受激励水平的影响。
(2)本发明在实现有限激励条件下参数估计的同时,基于非对称积分障碍李雅普诺夫函数直接将期望行为指标强加于输出轨迹上,实现了对系统输出的直接约束。一方面避免了传统障碍李雅普诺夫函数在施加约束时需要将输出约束问题转化为误差约束问题从而导致的系统保守性增加,另一方面将积分障碍李雅普诺夫函数的适用范围扩展到了非对称边界条件下。
(3)本发明在实现输出约束和参数估计的同时,将并行学习推广到有限时间控制领域,提出了一种有限时间复合学习控制器,保证了参数辨识和轨迹跟踪的有限时间收敛,切实提高动力定位船作业的安全性和有效性,具有很强的工程意义。
本发明主要应用于作业过程中输出受限且存在刚体惯量和水动力参数未知的船舶动力定位系统。
附图说明
图1为本发明所述一种动力定位船的有限时间复合学习控制方法的流程框图;
图2为北东坐标系和船体坐标系示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其它实施例,都属于本发明保护的范围。需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。
针对动力定位船在完成某些特殊作业任务时存在系统输出受限这一问题,以及船舶作业过程中受负载变化和燃料消耗等原因引起的刚体惯量不确定和水动力参数不确定这一问题,提出一种基于非对称积分障碍李雅普诺夫函数的有限时间复合学习控制方法,它是一种能够预先设定系统输出轨迹的性能边界,并在不使用加速度信息的情况下在线辨识系统未知参数的控制方法。基于回归滤波、并行学习思想和动态回归扩张与混合过程设计复合学习律,所有参数矢量分量的收敛速度可以灵活调整,且不受激励水平的影响。此外,利用非对称积分障碍李雅普诺夫函数直接将期望行为指标强加于输出轨迹上,确保船舶的实际输出轨迹满足作业任务的要求。此方法,可保证在有限激励条件下,参数估计误差和轨迹跟踪误差在有限时间内收敛到零,提高了动力定位船作业的安全性和有效性。具体如下:
具体实施方式一:参照图1和图2具体说明本实施方式,本实施方式所述的一种动力定位船的有限时间复合学习控制方法,包括以下步骤:
步骤一:建立北东坐标系XEYEZE和船体坐标系XBYBZB
北东坐标系XEYEZE为原点OE位于水面、XE轴和YE轴分别朝向北方和东方、ZE轴垂直于水面的空间直角坐标系。船体坐标系XBYBZB为原点OB位于船体重心,XB轴为船体的纵向、YB轴为船体的横向且ZB轴垂直于船面的空间直角坐标系。
以两个坐标系为基础建立动力定位船的运动学和动力学模型:
运动学模型:
Figure BDA0003852579720000071
动力学模型:
Figure BDA0003852579720000072
其中,η为动力定位船的实际位置,η=[η123]T,η1和η2分别为北东坐标系下船舶实际位置的XE轴和YE轴坐标,η3为船舶实际的姿态角,
Figure BDA0003852579720000073
为η的一阶导数;
υ=[u,v,r]T,u和v分别为动力定位船XB轴和YB轴的线速度,r为动力定位船的艏向回转率;
R(η3)为北东坐标系和船体坐标系之间的转换矩阵、且满足RT3)=R-13),
Figure BDA0003852579720000074
U为动力定位船的控制输入向量,M为动力定位船控制系统的惯性矩阵,D为动力定位船控制系统的阻尼矩阵,
Figure BDA0003852579720000075
Figure BDA0003852579720000076
参照动力定位船的期望位置ηd=[600cos(t/200)+200;600sin(t/200)+200;0]以及初始状态η(0)=[200;200;0]和υ(0)=[0;0;0],对其实际位置η施加约束条件:
Figure BDA0003852579720000077
并建立动力定位船控制系统的误差方程:
Figure BDA0003852579720000078
其中,z1为η和ηd之间的误差,z2
Figure BDA0003852579720000079
与虚拟控制率α之间的误差;-δ
Figure BDA00038525797200000710
分别为实际位置η的下限和上限,且δ=[425;430;5×10-3];
Figure BDA00038525797200000711
ηd=[ηd1d2d3]T,α=[α123]T
Figure BDA0003852579720000081
Figure BDA0003852579720000082
ai为误差增益系数,Λi为滑模项增益系数,本实施方式中a1=0.9,a2=0.85,a3=1,Λ1=5,Λ2=Λ3=1。z1i=ηi-ηdi,ηd1和ηd2分别为北东坐标系下船舶期望位置的XE轴和YE轴坐标,ηd3为船舶期望的姿态角。
步骤二:基于实际位置η的约束条件构造非对称积分障碍李雅普诺夫函数V1和李雅普诺夫函数V2
Figure BDA0003852579720000083
Figure BDA0003852579720000084
其中,
Figure BDA0003852579720000085
M*=R(η3)MRT3),τ为积分变量。
步骤三:对V2求一阶导数获得:
Figure BDA0003852579720000086
其中,f为动力定位船控制系统的不确定部分,其具体表达式为:
Figure BDA0003852579720000087
Figure BDA0003852579720000088
D*=R(η3)DRT3),
Figure BDA0003852579720000089
Figure BDA00038525797200000810
分别为α和R(η3)的一阶导数。
步骤四:对动力定位船控制系统的不确定部分f进行线性化处理获得:
f=Φθ (8),
其中,Φ为初始回归矩阵,θ为未知参数向量,
Figure BDA00038525797200000811
θ=[M11 M22 M23 M33 D11 D22 D23 D32 D33]T
步骤五:在初始回归矩阵Φ的基础上构建增广回归矩阵W:
Figure BDA0003852579720000091
其中,z2f为z2的滤波值,
Figure BDA0003852579720000092
z2f通过如下滤波器获得:
Figure BDA0003852579720000093
Figure BDA0003852579720000094
滤波器时间常数c=1。
步骤六:分别建立增广回归矩阵W和动力定位船的控制输入向量U的稳定滤波器
Figure BDA0003852579720000095
Figure BDA0003852579720000096
Figure BDA0003852579720000097
Figure BDA0003852579720000098
其中,Uf和Wf分别为U和W的滤波值,Wf(0)=0,Uf(0)=0。
步骤七:对z2求一阶导数并在其一阶导数的等号两边同时乘以M*获得:
Figure BDA0003852579720000099
建立北东坐标系下动力定位船的数学模型:
Figure BDA00038525797200000910
对上式进行变换获得:
Figure BDA00038525797200000911
利用回归滤波方法建立滤波系统动力学方程并执行仿射参数线性化获得:
Figure BDA00038525797200000912
将式(7)代入式(14)获得:
Figure BDA00038525797200000913
对上式进行变换获得:
Figure BDA00038525797200000914
将式(13)和(15)代入式(12)中整理获得:
Figure BDA0003852579720000101
由于
Figure BDA0003852579720000102
所以
Figure BDA0003852579720000103
则有:
Figure BDA0003852579720000104
将式(16)、(10)和(11)代入
Figure BDA0003852579720000105
的一阶导数中获得:
Figure BDA0003852579720000106
Figure BDA0003852579720000107
并对
Figure BDA0003852579720000108
求一阶导,将公式五代入
Figure BDA0003852579720000109
中,对其求解获得:
Figure BDA00038525797200001010
对Uf进行参数线性化获得:
Uf=Waθ (19),
其中,Wa为线性化后的回归矩阵,其具体表达式为:
Figure BDA00038525797200001011
其中,
Figure BDA00038525797200001012
步骤八:已知传统的复合学习律普遍存在秩亏问题,即仅当线性化后的回归矩阵Wa满足持续激励条件时才能保证参数收敛到其真值。因此,为了放宽这一必要条件,基于并行学习思想和动态回归扩张与混合过程来设计复合学习律,具体为:
首先,执行回归扩张,对Uf=Waθ等号左右两侧分别乘以
Figure BDA00038525797200001013
获得:
Figure BDA00038525797200001014
利用Laplace算子
Figure BDA00038525797200001015
建立信息矩阵
Figure BDA00038525797200001016
和辅助矩阵
Figure BDA00038525797200001017
且两矩阵的一阶导数分别为:
Figure BDA0003852579720000111
Figure BDA0003852579720000112
Laplace算子
Figure BDA0003852579720000113
的具体表达式为:
Figure BDA0003852579720000114
其中,s为算子变量,Laplace算子的滤波时间常数l=0.01。
步骤九:分别对信息矩阵Ω和辅助矩阵N进行积分,获得:
Figure BDA0003852579720000115
Figure BDA0003852579720000116
其中,t为积分时间上限,χ为积分变量。
则结合式(24)和(25)获得:
N(t)=Ω(t)θ (26)。
在式(26)等号的两边均乘以信息矩阵Ω的伴随矩阵adj(Ω)获得:
N(t)adj(Ω)=Ω(t)adj(Ω)θ (27),
使Ψ(t)=N(t)adj(Ω),Δ(t)=Ω(t)adj(Ω),则有:
Ψ(t)=Δ(t)θ (28)。
此时,只要Wa满足更易于实现的有限激励条件,即可使信息矩阵Ω达到满秩,进而保证未知参数的收敛,从而放宽了参数收敛的条件。值得注意的是,由于(22)中的信息矩阵Ω采用了指数遗忘设计,所以在激励信号消失后信息矩阵Ω将趋近于0,这将导致Δ(t)趋近于0从而大大降低参数收敛的性能。因此,为了在有限激励下获得良好的参数估计性能,不应使用整个时间间隔[0,t]内的所有信息来生成信息矩阵Ω。为此,设信息矩阵的更新截止时刻:
Figure BDA0003852579720000117
其中,
Figure BDA0003852579720000118
为信息矩阵Ω的积分上限,m=9为待确定向量中未知参数的个数。
将te代入Ψ(t)和Δ(t)分别获得:
te时刻信息矩阵的行列式:Δe=Δ(te),
te时刻更新的辅助矩阵:Ψe=Ψ(te)。
综上,为进一步提高收敛速度,利用te、Δe和Ψe基于并行学习思想和动态回归扩张与混合过程设计有限时间复合学习律:
Figure BDA0003852579720000121
其中,γ=0.5I9×9为正定自适应增益矩阵,Γ=5I9×9为正定调谐参数矩阵,
Figure BDA0003852579720000122
为θ的估计值,
Figure BDA0003852579720000123
Figure BDA0003852579720000124
的一阶导,
Figure BDA0003852579720000125
步骤十:为确保在存在输出约束和参数不确定的情况下,系统输出能快速准确的跟踪期望轨迹,进一步提高动力定位船作业的可靠性和有效性,基于非对称积分障碍李雅普诺夫函数和有限时间复合学习律,设计有限时间复合学习控制器:
Figure BDA0003852579720000126
其中,K=diag{1×104,1×104,1×104}为误差增益矩阵,H=diag{2000,2000,2000}为滑模项增益矩阵,
利用有限时间复合学习控制器输出的控制输入向量U控制动力定位船。
通过Matlab仿真,本实施方式的基于非对称积分障碍李雅普诺夫函数的动力定位船有限时间复合学习控制方法,可以在存在参数不确定的情况下实现对系统输出的约束,并保证参数估计误差和系统误差在有限时间内收敛到零。具有较强的灵活性和鲁棒性且能够有效提高动力定位船作业的可靠性和有效性。
本实施方式中所采用的复合学习控制方法相比于其他控制方法,具有在有限激励条件下在线辨识系统未知参数的能力,保证各参数向量分量的收敛速度不受激励水平的影响且可灵活调整。同时,因采用非对称积分障碍李雅普诺夫函数直接将期望行为指标强加于输出轨迹上,避免了采用传统障碍李雅普诺夫函数完成类似功能时所带来的系统保守性增加和计算量增大等问题。此外,本实施方式基于有限时间理论设计复合学习控制器,进一步提高了控制效果。
虽然在本文中参照了特定的实施方式来描述本发明,但是应该理解的是,这些实施例仅仅是本发明的原理和应用的示例。因此应该理解的是,可以对示例性的实施例进行许多修改,并且可以设计出其他的布置,只要不偏离所附权利要求所限定的本发明的精神和范围。应该理解的是,可以通过不同于原始权利要求所描述的方式来结合不同的从属权利要求和本文中所述的特征。还可以理解的是,结合单独实施例所描述的特征可以使用在其它所述实施例中。

Claims (8)

1.一种动力定位船的有限时间复合学习控制方法,其特征在于,包括以下步骤:
步骤一:参照动力定位船的期望位置ηd对其实际位置η施加约束条件:
Figure FDA0003852579710000011
并建立动力定位船控制系统的误差方程:
Figure FDA0003852579710000012
其中,z1为η和ηd之间的误差,
Figure FDA0003852579710000013
为η的一阶导数,z2
Figure FDA0003852579710000014
与虚拟控制率α之间的误差;-δ和
Figure FDA0003852579710000015
分别为实际位置η的下限和上限,η=[η123]T,ηd=[ηd1d2d3]T,α=[α123]T
Figure FDA0003852579710000016
Figure FDA0003852579710000017
ai为误差增益系数,Λi为滑模项增益系数,z1i=ηidi,η1和η2分别为北东坐标系下船舶实际位置的横纵坐标,ηd1和ηd2分别为北东坐标系下船舶期望位置的横纵坐标,η3和ηd3分别为船舶实际与期望的姿态角;
步骤二:基于实际位置η的约束条件构造非对称积分障碍李雅普诺夫函数V1和李雅普诺夫函数V2
Figure FDA0003852579710000018
其中,
Figure FDA0003852579710000019
M*=R(η3)MRT3),R(η3)为北东坐标系和船体坐标系之间的转换矩阵,RT3)=R-13),M为动力定位船控制系统的惯性矩阵,
Figure FDA00038525797100000110
τ为积分变量;
步骤三:对V2求一阶导数:
Figure FDA0003852579710000021
其中,f为动力定位船控制系统的不确定部分,
Figure FDA0003852579710000022
Figure FDA0003852579710000023
D*=R(η3)DRT3),D为动力定位船控制系统的阻尼矩阵,
Figure FDA0003852579710000024
Figure FDA0003852579710000025
Figure FDA0003852579710000026
分别为α和R(η3)的一阶导数,U为动力定位船的控制输入向量;
步骤四:对f进行线性化处理得:f=Φθ,
其中,Φ为初始回归矩阵,θ为未知参数向量,
Figure FDA0003852579710000027
θ=[M11 M22 M23 M33 D11 D22 D23 D32 D33]T
步骤五:在初始回归矩阵Φ的基础上构建增广回归矩阵W:
Figure FDA0003852579710000028
其中,z2f为z2的滤波值,
Figure FDA0003852579710000029
z2f通过如下滤波器获得:
Figure FDA00038525797100000210
c为滤波器时间常数;
步骤六:分别建立增广回归矩阵W和动力定位船的控制输入向量U的稳定滤波器
Figure FDA00038525797100000211
Figure FDA00038525797100000212
Figure FDA00038525797100000213
Figure FDA00038525797100000214
其中,Uf和Wf分别为U和W的滤波值;
步骤七:将步骤六中的两个滤波器代入
Figure FDA00038525797100000215
的一阶导数中求解获得:
Figure FDA00038525797100000216
对Uf进行参数线性化获得:
Uf=Waθ,
其中,Wa为线性化后的回归矩阵;
步骤八:对Uf=Waθ等号左右两侧分别乘以
Figure FDA0003852579710000031
获得:
Figure FDA0003852579710000032
利用Laplace算子
Figure FDA0003852579710000033
建立信息矩阵
Figure FDA0003852579710000034
和辅助矩阵
Figure FDA0003852579710000035
步骤九:基于并行学习思想和动态回归扩张与混合过程设计有限时间复合学习律:
Figure FDA0003852579710000036
其中,γ为正定自适应增益矩阵,Γ为正定调谐参数矩阵,te为信息矩阵的更新截止时刻,Δe为te时刻信息矩阵的行列式,Ψe为te时刻的辅助矩阵,
Figure FDA0003852579710000037
为θ的估计值,
Figure FDA0003852579710000038
Figure FDA0003852579710000039
的一阶导;
步骤十:利用有限时间复合学习律并基于非对称积分障碍李雅普诺夫函数设计有限时间复合学习控制器:
Figure FDA00038525797100000310
其中,K为误差增益矩阵,H为滑模项增益矩阵,
利用有限时间复合学习控制器输出的控制输入向量U控制动力定位船。
2.根据权利要求1所述的一种动力定位船的有限时间复合学习控制方法,其特征在于,北东坐标系XEYEZE为原点OE位于水面、XE轴和YE轴分别朝向北方和东方且ZE轴垂直于水面的空间直角坐标系。
3.根据权利要求2所述的一种动力定位船的有限时间复合学习控制方法,其特征在于,船体坐标系XBYBZB为原点OB位于船体重心,XB轴为船体的纵向、YB轴为船体的横向且ZB轴垂直于船面的空间直角坐标系。
4.根据权利要求3所述的一种动力定位船的有限时间复合学习控制方法,其特征在于,北东坐标系和船体坐标系之间的转换矩阵具体为:
Figure FDA0003852579710000041
5.根据权利要求1所述的一种动力定位船的有限时间复合学习控制方法,其特征在于,步骤七中求解Uf的具体方法为:
对z2求一阶导数并在其一阶导数的等号两边同时乘以M*获得:
Figure FDA0003852579710000042
建立北东坐标系下动力定位船的数学模型:
Figure FDA0003852579710000043
利用回归滤波方法建立滤波系统动力学方程并执行仿射参数线性化,获得:
Figure FDA0003852579710000044
Figure FDA0003852579710000045
公式二和公式三均代入公式一中获得:
Figure FDA0003852579710000046
将公式四以及稳定滤波器
Figure FDA0003852579710000047
Figure FDA0003852579710000048
代入
Figure FDA0003852579710000049
的一阶导数中获得:
Figure FDA00038525797100000410
Figure FDA00038525797100000411
并对
Figure FDA00038525797100000412
求一阶导,将公式五代入
Figure FDA00038525797100000413
中,对其求解获得:
Figure FDA00038525797100000414
6.根据权利要求1或5所述的一种动力定位船的有限时间复合学习控制方法,其特征在于,步骤七中线性化后的回归矩阵Wa的具体表达式为:
Figure FDA00038525797100000415
其中,
Figure FDA00038525797100000416
7.根据权利要求1所述的一种动力定位船的有限时间复合学习控制方法,其特征在于,步骤八中,设计有限时间复合学习律的具体方法为:
分别对信息矩阵Ω和辅助矩阵N积分,获得:
Figure FDA0003852579710000051
Figure FDA0003852579710000052
其中,l为Laplace算子的滤波时间常数,t为积分时间上限,χ为积分变量,
结合公式六和七获得:
N(t)=Ω(t)θ 公式八,
在公式八等号的两边均乘以adj(Ω)获得:
Ψ(t)=Δ(t)θ,
其中,Ψ(t)=N(t)adj(Ω),Δ(t)=Ω(t)adj(Ω),adj(Ω)为信息矩阵Ω的伴随矩阵,
Figure FDA0003852579710000053
其中,
Figure FDA0003852579710000054
为信息矩阵Ω的积分上限,m为待确定向量中未知参数的个数,
将te代入Ψ(t)和Δ(t)分别获得:
Δe=Δ(te),Ψe=Ψ(te),
利用te、Δe和Ψe建立有限时间复合学习律。
8.根据权利要求1或7所述的一种动力定位船的有限时间复合学习控制方法,其特征在于,
Laplace算子
Figure FDA0003852579710000055
的具体表达式为:
Figure FDA0003852579710000056
其中,s为算子变量。
CN202211138901.8A 2022-06-10 2022-09-19 一种动力定位船的有限时间复合学习控制方法 Active CN115309058B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2022106571913 2022-06-10
CN202210657191 2022-06-10

Publications (2)

Publication Number Publication Date
CN115309058A true CN115309058A (zh) 2022-11-08
CN115309058B CN115309058B (zh) 2023-06-30

Family

ID=83866690

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211138901.8A Active CN115309058B (zh) 2022-06-10 2022-09-19 一种动力定位船的有限时间复合学习控制方法

Country Status (1)

Country Link
CN (1) CN115309058B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116088309A (zh) * 2023-01-03 2023-05-09 哈尔滨理工大学 一种基于故障辨识的水面船复合学习容错控制方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004118658A (ja) * 2002-09-27 2004-04-15 Advanced Telecommunication Research Institute International 物理系の制御方法および装置ならびに物理系の制御のためのコンピュータプログラム
CN107861384A (zh) * 2017-11-05 2018-03-30 西北工业大学 基于复合学习的mems陀螺仪快速启动方法
CN108828957A (zh) * 2018-08-20 2018-11-16 西北工业大学 基于切换机制的飞行器全局有限时间神经网络控制方法
CN112394740A (zh) * 2019-08-12 2021-02-23 天津大学 一种无人驾驶碾压机的复合抗扰轨迹跟踪控制算法
CN112612209A (zh) * 2020-12-14 2021-04-06 哈尔滨工业大学 基于指令滤波神经网络控制器的全驱动船舶轨迹跟踪控制方法及系统
US20210133555A1 (en) * 2019-10-30 2021-05-06 Microsoft Technology Licensing, Llc Distributed learning of composite machine learning models
CN113157000A (zh) * 2021-05-06 2021-07-23 西北工业大学 基于虚拟结构和人工势场的飞行编队协同避障自适应控制方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004118658A (ja) * 2002-09-27 2004-04-15 Advanced Telecommunication Research Institute International 物理系の制御方法および装置ならびに物理系の制御のためのコンピュータプログラム
CN107861384A (zh) * 2017-11-05 2018-03-30 西北工业大学 基于复合学习的mems陀螺仪快速启动方法
CN108828957A (zh) * 2018-08-20 2018-11-16 西北工业大学 基于切换机制的飞行器全局有限时间神经网络控制方法
CN112394740A (zh) * 2019-08-12 2021-02-23 天津大学 一种无人驾驶碾压机的复合抗扰轨迹跟踪控制算法
US20210133555A1 (en) * 2019-10-30 2021-05-06 Microsoft Technology Licensing, Llc Distributed learning of composite machine learning models
CN112612209A (zh) * 2020-12-14 2021-04-06 哈尔滨工业大学 基于指令滤波神经网络控制器的全驱动船舶轨迹跟踪控制方法及系统
CN113157000A (zh) * 2021-05-06 2021-07-23 西北工业大学 基于虚拟结构和人工势场的飞行编队协同避障自适应控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
戴永强,傅佩琛,马秀娟,强文义: "基于神经网络的复合学习方案及在机器人控制中的应用", 哈尔滨工业大学学报, no. 04 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116088309A (zh) * 2023-01-03 2023-05-09 哈尔滨理工大学 一种基于故障辨识的水面船复合学习容错控制方法
CN116088309B (zh) * 2023-01-03 2023-07-21 哈尔滨理工大学 一种基于故障辨识的水面船复合学习容错控制方法

Also Published As

Publication number Publication date
CN115309058B (zh) 2023-06-30

Similar Documents

Publication Publication Date Title
CN114047773B (zh) 一种基于扩张状态观测器的水下集矿机器人反步滑模自适应姿态控制方法
Zhang et al. Linear reduction of backstepping algorithm based on nonlinear decoration for ship course-keeping control system
CN112612209B (zh) 基于指令滤波神经网络控制器的全驱动船舶轨迹跟踪控制方法及系统
CN110377036A (zh) 一种基于指令约束的无人水面艇航迹跟踪固定时间控制方法
CN105242683B (zh) 一种飞艇神经网络终端滑模航迹控制方法
CN111198502B (zh) 基于干扰观测器和模糊系统的无人艇航迹跟踪控制方法
CN113110511B (zh) 一种基于广义模糊双曲模型的智能船舶航向控制方法
CN105974930A (zh) 一种基于非线性模型预测控制的uuv对运动母船的跟踪方法
CN110703605B (zh) 一种面向智能船舶自动舵系统的自适应模糊最优控制方法及系统
CN115309058A (zh) 一种动力定位船的有限时间复合学习控制方法
CN111948937A (zh) 多智能体系统的多梯度递归强化学习模糊控制方法及系统
Morishita et al. Modified observer backstepping controller for a dynamic positioning system
Das et al. GESO based robust output tracking controller for marine vessels
Miao et al. DOPH∞-based path-following control for underactuated marine vehicles with multiple disturbances and constraints
Wang et al. Fuzzy iterative sliding mode control applied for path following of an autonomous underwater vehicle with large inertia
CN110968961A (zh) 一种连续回转电液伺服马达参数辨识方法
Fu et al. A unified barrier function-based trajectory tracking control for hovercraft with asymmetric dynamic safety constraints
Do et al. Global waypoint tracking control of underactuated ships under relaxed assumptions
Hwang et al. Applications of polytopic approximations of reachable sets to linear dynamic games and a class of nonlinear systems
CN108614922B (zh) 一种基于应力杂交化后处理的板壳边缘效应的计算方法
Timchenko et al. Robust stabilization of marine mobile objects on the basis of systems with variable structure of feedbacks
CN112904719B (zh) 一种适用于水下机器人位置环形区域跟踪控制方法
CN111856941B (zh) 一种基于主动抗扰的自适应终端动态滑模控制方法
CN110263400B (zh) 一种无人艇非线性动力学模型的积分辨识方法
Udin et al. Using Pontryagin maximum principle for parametrical identification of ship maneuvering mathematical model

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