CN111025909A - 船舶运动控制系统的Kalman三自由度解耦滤波方法 - Google Patents
船舶运动控制系统的Kalman三自由度解耦滤波方法 Download PDFInfo
- Publication number
- CN111025909A CN111025909A CN201911337502.2A CN201911337502A CN111025909A CN 111025909 A CN111025909 A CN 111025909A CN 201911337502 A CN201911337502 A CN 201911337502A CN 111025909 A CN111025909 A CN 111025909A
- Authority
- CN
- China
- Prior art keywords
- ship
- matrix
- kalman
- freedom
- motion
- 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
Images
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/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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
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)
Abstract
本发明属于船舶运动控制技术领域,具体涉及一种船舶运动控制系统的Kalman三自由度解耦滤波方法。本发明利用DP控制系统获取的传感器测量信息和连续型Kalman滤波算法,根据DP船舶的船舶模型和海洋环境参数,设计了Kalman三自由度解耦滤波方法,对船舶的运动信息进行重组和计算,实现位置及航向信息的滤波,并估计出DP船舶的实际位置、航向及一阶波浪力干扰力。通过仿真实现,证明本发明滤波效果明显,有效地消除了一阶波浪力的干扰。本发明对于DP船舶的运动控制滤波来说具有重要的研究意义和工程价值,有利于降低船舶运动控制器的设计难度,减少推进系统的无效推力输出及旋转。
Description
技术领域
本发明属于船舶运动控制技术领域,具体涉及一种船舶运动控制系统的Kalman三自由度解耦滤波方法。
背景技术
当今全球各国的经济正在高速的发展中,各国对能源的消耗也在与日俱增,众多科学家早已将目光投向具有丰富自然矿藏资源的海洋领域。面对海洋环境的复杂多变,必须使得船舶拥有先进的控制系统,因此一个新的系泊方式—动力定位系统诞生了。动力定位系统可以使船舶在深水海域以及复杂环境海域下保持自身位置的恒定以及运动航迹的稳定,同时可以迅速进入以及离开工程地点。
DP船舶在海洋中运动时,周围海洋环境的扰动会明显地影响船舶的运动与控制,而且扰动具有较强的波动性和随机性,其对于研究船舶的运动以及响应特性来说,会产生一定的干扰。船舶运动过程中,主要影响船舶运动的环境因素有风、浪、流,其中海浪导致引起产生的一阶波浪力会导致船舶发生高频运动,并且这种运动是均值为零的往复运动,船舶不应抵抗的运动,如果船舶抵抗响应这种高频运动就会引起很多不必要的控制,进而使得船舶的推进系统会不断的响应这种高频运动,加速推进系统的磨损老化,消耗过多的能量。
如果在DP控制系统中加入了合适的滤波器,滤波器能将位置及航向信号中的的高频运动滤除并估计,同时保留低频信号,并估计船舶的实际位置,从而计算出测量系统中的高斯白噪声,使船舶仅仅响应低频运动的扰动,必然会提升系统的控制以响应能力,同时减少不必要的推力输出和推进器回转,提升动力定位系统的能力。
发明内容
本发明的目的在于提供实现船舶的高精度运动滤波控制的一种船舶运动控制系统的Kalman三自由度解耦滤波方法。
本发明的目的通过如下技术方案来实现:包括以下步骤:
步骤1:依据DP船舶的船舶模型获取无因次质量矩阵M和阻尼矩阵D;
步骤2:依据DP船舶在港口不运动时获取的数据记录的样本协方差,获取测量噪声协方差矩阵R;
步骤3:设定具有正可调参数的对角线的状态噪声协方差Q;
步骤4:设置误差协方差矩阵的初始值为P0=I15×15,通过DP船舶实际位置确定初始位置x0;
步骤5:依据海浪强度的参数、阻尼系数、主导海浪频率及船舶输出推力矩计算确定偏差矩阵b、矩阵E和矩阵A;
其中:
步骤6:依据DP船舶需要滤除一阶波浪力干扰的要求,确定矩阵H;
H=[Cω I3×3 03×3 03×3]
Cω=[03×3 I3×3]
步骤7:依据连续型Kalman滤波方法求解总的Kalman增益,并进行解耦获得单自由度Kalman增益;
步骤9:依据矩阵A、DP船舶动力学模型和过程噪声向量,重组船舶的运动结构,获取船舶的估计运动位置和一阶波浪力干扰;
步骤11:依据每一时刻DP船舶的输出推力(矩)和传感器测量的船舶位置信息,使用连续型Kalman滤波算法迭代计算相关矩阵,确定Kalman增益,重组船舶运动结构,获得船舶的估计运动位置和一阶波浪力;在连续型Kalman滤波器中,状态估计迭代式和误差协方差迭代式如下:
初始条件:P(0)=P0;
其中,y=Hx+υ;y为传感器测量出的纵荡、横荡位置及航向角;x为状态变量,ηP=RT(ψ)η;η为在北东坐标系下DP船舶三自由度向量,η=[N E ψ]T,分别表示纵荡、横荡和艏摇;R(ψ)为转换矩阵,具体为:
v为船体坐标系下的速度向量,v=[u v r]T;bP=RT(ψ)b。
本发明的有益效果在于:
本发明设计了一种船舶运动控制系统的Kalman三自由度解耦滤波方法,利用DP控制系统获取的传感器测量信息和连续型Kalman滤波算法,根据DP船舶的船舶模型和海洋环境参数,设计了Kalman三自由度解耦滤波方法,对船舶的运动信息进行重组和计算,实现位置及航向信息的滤波,并估计出DP船舶的实际位置、航向及一阶波浪力干扰力。通过仿真实现,证明本发明滤波效果明显,有效地消除了一阶波浪力的干扰。本发明对于DP船舶的运动控制滤波来说具有重要的研究意义和工程价值,有利于降低船舶运动控制器的设计难度,减少推进系统的无效推力输出及旋转。
附图说明
图1为三自由度Kalman增益计算及其解耦仿真图。
图2为每个单自由度的滤波方法原理图。
图3为单自由度的滤波方法仿真图。
图4为DP船舶运动控制滤波的仿真图。
图5为利用Matlab仿真得到的北向位置滤波估计曲线图。
图6为利用Matlab仿真得到的东向位置滤波估计曲线图。
图7为利用Matlab仿真得到的航向角度滤波估计曲线图。
图8为利用Matlab仿真得到的100s至150s航向角度滤波估计放大曲线图。
图9为利用Matlab仿真得到的一阶波浪力估计曲线图。
具体实施方式
下面结合附图对本发明做进一步描述。
本发明涉及一种船舶运动控制系统的Kalman三自由度解耦滤波方法,特别是在DP船舶在海洋中进行定位时,利用DP控制系统获取的测量信息,基于Kalman三自由度解耦的滤波器对测量的位置及航向信息进行滤波,估计出DP船舶的实际位置及航向。本发明充分利用DP船舶获取的船舶位置及航向的测量信息,使用DP船舶的线性模型、连续Kalman滤波算法和相关滤波器参数的设定,计算出DP船舶三自由度的Kalman增益,进行解耦得出纵荡、横荡及航向上的Kalman增益,实现单自由度上的滤波及控制,重组DP船舶的运动,估计并滤除一阶高频波浪力,使DP船舶抵抗船舶的低频运动,实现Kalman三自由度解耦滤波,有利于DP船舶精准的定位控制,减少船舶推进系统不必要的控制响应。
一种船舶运动控制系统的Kalman三自由度解耦滤波方法,包括以下步骤:
步骤1:依据DP船舶的船舶模型获取无因次质量矩阵M和阻尼矩阵D;
步骤2:依据DP船舶在港口不运动时获取的数据记录的样本协方差,获取测量噪声协方差矩阵R;
步骤3:设定具有正可调参数的对角线的状态噪声协方差Q;
步骤4:设置误差协方差矩阵的初始值为P0=I15×15,通过DP船舶实际位置确定初始位置x0;
步骤5:依据海浪强度的参数、阻尼系数、主导海浪频率及船舶输出推力矩计算确定偏差矩阵b、矩阵E和矩阵A;
其中:
步骤6:依据DP船舶需要滤除一阶波浪力干扰的要求,确定矩阵H;
H=[Cω I3×3 03×3 03×3]
Cω=[03×3 I3×3]
步骤7:依据连续型Kalman滤波方法求解总的Kalman增益,并进行解耦获得单自由度Kalman增益;
步骤9:依据矩阵A、DP船舶动力学模型和过程噪声向量,重组船舶的运动结构,获取船舶的估计运动位置和一阶波浪力干扰;
步骤11:依据每一时刻DP船舶的输出推力(矩)和传感器测量的船舶位置信息,使用连续型Kalman滤波算法迭代计算相关矩阵,确定Kalman增益,重组船舶运动结构,获得船舶的估计运动位置和一阶波浪力;在连续型Kalman滤波器中,状态估计迭代式和误差协方差迭代式如下:
初始条件:P(0)=P0;
其中,y=Hx+υ;y为传感器测量出的纵荡、横荡位置及航向角;x为状态变量,ηP=RT(ψ)η;η为在北东坐标系下DP船舶三自由度向量,η=[N E ψ]T,分别表示纵荡、横荡和艏摇;R(ψ)为转换矩阵,具体为:
v为船体坐标系下的速度向量,v=[u v r]T;bP=RT(ψ)b。
本发明设计了一种船舶运动控制系统的Kalman三自由度解耦滤波方法,利用DP控制系统获取的传感器测量信息和连续型Kalman滤波算法,根据DP船舶的船舶模型和海洋环境参数,设计了Kalman三自由度解耦滤波方法,对船舶的运动信息进行重组和计算,实现位置及航向信息的滤波,并估计出DP船舶的实际位置、航向及一阶波浪力干扰力。在Matlab中进行了仿真验证,其滤波效果明显,有效地消除了一阶波浪力的干扰,证明了一种船舶运动控制系统的Kalman三自由度解耦滤波方法的有效性和真实性。本发明对于DP船舶的运动控制滤波来说具有重要的研究意义和工程价值,有利于降低船舶运动控制器的设计难度,减少推进系统的无效推力输出及旋转。
本发明的推导步骤如下:
第一步,对于DP船舶动力定位系统来说,一般来说只考虑水面上的三个自由度(纵荡、横荡、艏摇)。选择两个DP系统常用的参考系,北东坐标系和船体坐标系,定义在北东坐标系下船舶三自由度向量η和船体坐标系下的速度向量v分别为:
η=[N E ψ]T (1)
v=[u v r]T (2)
两坐标系的转换关系如下:
式中:R(ψ)为转换矩阵。
因此,使用船舶平行坐标来表示船舶DP模型,有下式成立:
ηP=RT(ψ)η (5)
bP=RT(ψ)b (6)
所以运动学模型可以近似为线性模型:
选用DP的线性运动模型为:
式中:M为包括了附加质量的运动系统的惯性矩阵;D为代表了阻尼系数的矩阵;τ为船舶的推进系统输出的力及力矩矢量矩阵;bp为力(矩)的偏差向量;Tb为含有正偏差时间常数的对角矩阵;ωi(i=1,2,3)为过程噪声变量。
第三步,建立线性连续时间系统的状态空间方程和测量方程(传感器系统)为:
y=Hx+υ (14)
如果上述系统方程是可观测的,状态向量x∈Rn可以是通过测量向量y∈Rn和控制输入向量u∈Rn重新递推求解计算的。
连续时间的Kalman滤波算法如下:
设计矩阵:Q(t)=QT(t)>0(通常为常数) (15)
R(t)=RT(t)>0(通常为常数) (16)
卡尔曼增益矩阵迭代式:K(t)=P(t)HT(t)R-1(t) (19)
H=[Cω I3×3 03×3 03×3] (25)
系统矩阵中Aω、Eω和Eω的设定如下:
Cω=[03×3 I3×3] (28)
式中,ω0为波浪PM谱中的主导海洋频率;λ为相对阻尼系数;σ为与波浪的强度相关的常数值。
第五步,在连续型Kalman滤波器中,状态估计迭代式和误差协方差迭代式如下:
式中,协方差矩阵Q=QT∈R9×9和测量方差矩阵R=RT∈R3×3必须由使用者指定;其中测量协方差矩阵R的的方法可以使用下述矩阵:
一般来说,状态噪声协方差Q可以选择为具有正可调参数的对角线。这些通常是通过试验和错误发现的。将这个协方差矩阵选为对角矩阵,即:
矩阵Q1是噪声ω1的协方差,用来驱动代表线性波浪运动的噪声滤波器,其可与船舶WF运动模型的参数一起根据船舶运行前和运行中测量的数据进行估计;Q2是噪声ω2的协方差,代表运动方程中的不确定性,同时Q2也被视为位置测量噪声方差的一部分;Q3是噪声ω3的协方差,代表了模拟其他环境力量的偏差项的不确定性,Q3可以给滤波器提供模式各部分不确定性的适当平衡。
第六步,根据以上方法计算出的Kalman增益K是15×1的矩阵,如状态估计迭代矩阵公式(29)所示,这里对其进行解耦,将三自由度分开进行运算,纵荡、横荡和艏摇分别进行单独的滤波和状态估计。即:将计算出的Kalman增益K分成15组变量,第1、第4、第7、第10、第13,5个变量重新组成一组作为纵荡上的Kalman增益KN;第2、第5、第8、第11、第14,5个变量重新组成一组作为横荡上的Kalman增益KE;第3、第6、第9、第12、第15,5个变量重新组成一组作为艏摇上的Kalman增益KP,再将三个增益KN、KE、KP分别与计算出的纵荡、横荡、艏摇上的噪声相乘,从而将三自由度解耦成3个单自由度来进行滤波,其三自由度Kalman增益计算及其解耦仿真图如图1所示,此方法极大的简化了计算过程以及简便了程序设计的难度。
第七步,每个单自由度的滤波方法如下如图2所示,由此可以获得船舶的实际位置信息、三自由度上的速度(角速度)以及一阶波浪力估计结果。
由图2可知,在单自由度的滤波过程中,由于船舶模型的原因,运动重组过程中需要进行三自由度的合成,再进行矩阵运算,然后解耦,进行矩阵运算的相关矩阵为M-1、D-1、R-1(ψ)和R(ψ),从而进行船舶运动状态的重组以及一阶高频运动的估计,滤除单自由度上的一阶高频运动,输出单自由度上的低频运动,得出测量系统中的估计高斯白噪声,实现单自由度的滤波过程,其仿真图如图3所示。
将三个单自由度上的位置、航向和速度(角速度)信息进行合成,输入到控制器中,作为DP船舶实际运动的信息进行控制计算。
本发明方法给出了船舶运动控制系统的Kalman三自由度解耦滤波方法,根据发明步骤并对其进行仿真验证,对于DP船舶的高精度运动控制具有重要的理论意义和工程价值。
对本发明有益的说明:
1.仿真条件设置
仿真条件下,设海浪强度的参数为0.5、阻尼系数为0.1和主导海浪频率为0.8,为已知量,DP船舶的船舶模型的无因次质量矩阵M和阻尼矩阵D为:
误差协方差矩阵的初始值为P0=I15×15,测量噪声协方差矩阵R、偏差矩阵b以及状态噪声协方差矩阵Q1、Q2和Q3如下所示:
设定初始位置为北东坐标系下,纵荡位置为0m,横荡位置为0m,航向角为0°,仿真时间为200s。
2.仿真实验过程
根据DP船舶运动控制滤波的仿真模型在Matlab中进行仿真分析,具体过程如下:
将上述参数输入到连续型Kalman滤波器中,再Simulink中进行DP船舶控制系统的滤波和状态估计仿真分析,DP船舶运动控制滤波的仿真图如图4所示。
设定北东坐标系下,期望纵荡位置为0m,期望横荡位置为0m,期望航向角为20°。
将期望位置、航向及估计速度输入到DP船舶控制器中,通过PID控制器计算出DP船舶所需三自由度推力(矩)。
将推力(矩)输入到DP船舶的动力学模型中,计算出船舶的三自由度运动位置信息,加上均值为零的一阶波浪力干扰,作为传感器获得的测量信息。
将传感器获得的测量信息输入到连续型Kalman滤波器中,经过上述方法计算,可得估计的DP船舶实际位置、航向及估计的一阶波浪力干扰。
将估计的DP船舶实际位置和运动速度反馈到DP船舶控制器中,进行下一次计算。
将滤波前后DP船舶的运动状态及估计的一阶波浪力进行保存,输出仿真曲线。
将根据滤波前DP船舶的位置、航向及一阶波浪力的仿真曲线进行对比及结果分析,得出结论。
3.仿真实验结果
上述仿真结果表明:
(1)本发明一种船舶运动控制系统的Kalman三自由度解耦滤波方法可以在纵荡、横荡及航向上有效的滤除一阶波浪力的干扰,估计出DP船舶的实际运动情况。
(2)本发明一种船舶运动控制系统的Kalman三自由度解耦滤波方法可以滤除超90%的一阶波浪力,估计的一阶波浪力与实际的一阶波浪力基本吻合。
(3)本发明一种船舶运动控制系统的Kalman三自由度解耦滤波方法可以获得实际的DP船舶运动情况,有利于船舶精准的定位控制,减少推进器的无效响应。
本发明的具体实施方式为:
第一步,依据DP船舶的船舶模型获取无因次质量矩阵M和阻尼矩阵D。
第二步,依据DP船舶在港口不运动时获取的数据记录的样本协方差来估计矩阵R;
第三步,通过试验和错误发现具有正可调参数的对角线的状态噪声协方差Q,一般取值如下:
第四步,设置误差协方差矩阵的初始值为P0=I15×15,通过DP船舶实际位置确定初始位置x0;
第五步,依据海浪强度的参数、阻尼系数、主导海浪频率及船舶输出推力(矩)计算确定偏差矩阵b、矩阵EQET和矩阵A;
第六步,依据DP船舶需要滤除一阶波浪力干扰的要求,确定矩阵H;
第七步,依据连续型Kalman滤波方法求解总的Kalman增益,并进行解耦获得单自由度Kalman增益;
第九步,依据矩阵A、DP船舶动力学模型公式(8)至(11)和过程噪声向量,重组船舶的运动结构,获取船舶的估计运动位置和一阶波浪力干扰;
第十一步,依据每一时刻DP船舶的输出推力(矩)和传感器测量的船舶位置信息,使用连续型Kalman滤波算法迭代计算相关矩阵,确定Kalman增益,重组船舶运动结构,获得船舶的估计运动位置和一阶波浪力。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种船舶运动控制系统的Kalman三自由度解耦滤波方法,其特征在于,包括以下步骤:
步骤1:依据DP船舶的船舶模型获取无因次质量矩阵M和阻尼矩阵D;
步骤2:依据DP船舶在港口不运动时获取的数据记录的样本协方差,获取测量噪声协方差矩阵R;
步骤3:设定具有正可调参数的对角线的状态噪声协方差Q;
步骤4:设置误差协方差矩阵的初始值为P0=I15×15,通过DP船舶实际位置确定初始位置x0;
步骤5:依据海浪强度的参数、阻尼系数、主导海浪频率及船舶输出推力矩计算确定偏差矩阵b、矩阵E和矩阵A;
其中:
步骤6:依据DP船舶需要滤除一阶波浪力干扰的要求,确定矩阵H;
H=[Cω I3×3 03×3 03×3]
Cω=[03×3 I3×3]
步骤7:依据连续型Kalman滤波方法求解总的Kalman增益,并进行解耦获得单自由度Kalman增益;
步骤9:依据矩阵A、DP船舶动力学模型和过程噪声向量,重组船舶的运动结构,获取船舶的估计运动位置和一阶波浪力干扰;
步骤11:依据每一时刻DP船舶的输出推力(矩)和传感器测量的船舶位置信息,使用连续型Kalman滤波算法迭代计算相关矩阵,确定Kalman增益,重组船舶运动结构,获得船舶的估计运动位置和一阶波浪力;在连续型Kalman滤波器中,状态估计迭代式和误差协方差迭代式如下:
初始条件:P(0)=P0;
其中,y=Hx+υ;y为传感器测量出的纵荡、横荡位置及航向角;x为状态变量,ηP=RT(ψ)η;η为在北东坐标系下DP船舶三自由度向量,η=[N E ψ]T,分别表示纵荡、横荡和艏摇;R(ψ)为转换矩阵,具体为:
v为船体坐标系下的速度向量,v=[u v r]T;bP=RT(ψ)b。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911337502.2A CN111025909B (zh) | 2019-12-23 | 2019-12-23 | 船舶运动控制系统的Kalman三自由度解耦滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911337502.2A CN111025909B (zh) | 2019-12-23 | 2019-12-23 | 船舶运动控制系统的Kalman三自由度解耦滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111025909A true CN111025909A (zh) | 2020-04-17 |
CN111025909B CN111025909B (zh) | 2023-02-14 |
Family
ID=70211684
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911337502.2A Active CN111025909B (zh) | 2019-12-23 | 2019-12-23 | 船舶运动控制系统的Kalman三自由度解耦滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111025909B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112083663A (zh) * | 2020-09-14 | 2020-12-15 | 大连海事大学 | 船舶推进电机在四种海况环境下控制系统仿真建模方法 |
CN114111673A (zh) * | 2021-12-08 | 2022-03-01 | 中国船舶科学研究中心 | 一种六自由度相对运动测量系统 |
CN117572408A (zh) * | 2024-01-17 | 2024-02-20 | 中国人民解放军海军航空大学 | 一种大测量误差下雷达高速高机动目标跟踪方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090043436A1 (en) * | 2005-04-15 | 2009-02-12 | Kazuyuki Igarashi | Automatic Vessel Position Holding Control Method and Controller |
CN101920762A (zh) * | 2009-06-09 | 2010-12-22 | 同济大学 | 一种基于噪声矩阵实时修正的船舶动力定位方法 |
CN103217160A (zh) * | 2013-04-22 | 2013-07-24 | 哈尔滨工程大学 | 一种船舶动力定位卡尔曼滤波方法 |
CN105867165A (zh) * | 2016-04-01 | 2016-08-17 | 哈尔滨工程大学 | 基于扩展卡尔曼滤波器的动力定位船舶波频模型参数估计系统 |
CN107272415A (zh) * | 2017-07-25 | 2017-10-20 | 青岛科技大学 | 一种船舶动力定位控制系统非线性滤波器设计方法 |
CN107742026A (zh) * | 2017-10-16 | 2018-02-27 | 江苏科技大学 | 一种船舶动力定位系统非线性状态估计方法 |
CN108388738A (zh) * | 2018-03-01 | 2018-08-10 | 青岛科技大学 | 一种船舶动力定位系统噪声和状态实时估计自适应滤波方法 |
CN108445762A (zh) * | 2018-03-21 | 2018-08-24 | 青岛科技大学 | 基于ukf滤波和约束模型预测控制的动力定位控制方法 |
CN108460210A (zh) * | 2018-03-01 | 2018-08-28 | 青岛科技大学 | 一种基于噪残差和协方差匹配的动力定位系统噪声特性实时估计方法 |
US20180290719A1 (en) * | 2017-04-10 | 2018-10-11 | Mitsubishi Electric Corporation | Motion control device and motion control method for ship |
-
2019
- 2019-12-23 CN CN201911337502.2A patent/CN111025909B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090043436A1 (en) * | 2005-04-15 | 2009-02-12 | Kazuyuki Igarashi | Automatic Vessel Position Holding Control Method and Controller |
CN101920762A (zh) * | 2009-06-09 | 2010-12-22 | 同济大学 | 一种基于噪声矩阵实时修正的船舶动力定位方法 |
CN103217160A (zh) * | 2013-04-22 | 2013-07-24 | 哈尔滨工程大学 | 一种船舶动力定位卡尔曼滤波方法 |
CN105867165A (zh) * | 2016-04-01 | 2016-08-17 | 哈尔滨工程大学 | 基于扩展卡尔曼滤波器的动力定位船舶波频模型参数估计系统 |
US20180290719A1 (en) * | 2017-04-10 | 2018-10-11 | Mitsubishi Electric Corporation | Motion control device and motion control method for ship |
CN107272415A (zh) * | 2017-07-25 | 2017-10-20 | 青岛科技大学 | 一种船舶动力定位控制系统非线性滤波器设计方法 |
CN107742026A (zh) * | 2017-10-16 | 2018-02-27 | 江苏科技大学 | 一种船舶动力定位系统非线性状态估计方法 |
CN108388738A (zh) * | 2018-03-01 | 2018-08-10 | 青岛科技大学 | 一种船舶动力定位系统噪声和状态实时估计自适应滤波方法 |
CN108460210A (zh) * | 2018-03-01 | 2018-08-28 | 青岛科技大学 | 一种基于噪残差和协方差匹配的动力定位系统噪声特性实时估计方法 |
CN108445762A (zh) * | 2018-03-21 | 2018-08-24 | 青岛科技大学 | 基于ukf滤波和约束模型预测控制的动力定位控制方法 |
Non-Patent Citations (4)
Title |
---|
IVAN POPOV ET AL.: ""Adaptive Kalman Filtering for Dynamic Positioning of Marine Vessels"", 《IFAC-PAPERSONLINE》 * |
杜佳璐 等: ""船舶动力定位系统的非线性输出反馈控制"", 《华南理工大学学报(自然科学版)》 * |
钱佳 等: ""船舶动力定位非线性自适应滤波算法"", 《舰船科学技术》 * |
陈勇嘉 等: ""船舶动力定位控制系统的非线性观测器设计"", 《舰船电子工程》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112083663A (zh) * | 2020-09-14 | 2020-12-15 | 大连海事大学 | 船舶推进电机在四种海况环境下控制系统仿真建模方法 |
CN114111673A (zh) * | 2021-12-08 | 2022-03-01 | 中国船舶科学研究中心 | 一种六自由度相对运动测量系统 |
CN114111673B (zh) * | 2021-12-08 | 2023-05-12 | 中国船舶科学研究中心 | 一种六自由度相对运动测量系统 |
CN117572408A (zh) * | 2024-01-17 | 2024-02-20 | 中国人民解放军海军航空大学 | 一种大测量误差下雷达高速高机动目标跟踪方法 |
CN117572408B (zh) * | 2024-01-17 | 2024-04-16 | 中国人民解放军海军航空大学 | 一种大测量误差下雷达高速高机动目标跟踪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111025909B (zh) | 2023-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107168312B (zh) | 一种补偿uuv运动学和动力学干扰的空间轨迹跟踪控制方法 | |
CN111025909B (zh) | 船舶运动控制系统的Kalman三自由度解耦滤波方法 | |
CN111650948B (zh) | 一种可底栖式auv的水平面轨迹快速跟踪控制方法 | |
CN104898688B (zh) | Uuv四自由度动力定位自适应抗扰滑模控制系统及控制方法 | |
Hassani et al. | Robust dynamic positioning of offshore vessels using mixed-μ synthesis modeling, design, and practice | |
Banazadeh et al. | Frequency domain identification of the Nomoto model to facilitate Kalman filter estimation and PID heading control of a patrol vessel | |
CN111158383B (zh) | 基于干扰观测器和rbfnn的无人艇航迹跟踪控制方法 | |
CN109917657B (zh) | 动力定位船舶的抗干扰控制方法、装置与电子设备 | |
CN109739248B (zh) | 基于adrc的船载三自由度并联稳定平台稳定控制方法 | |
CN105867417B (zh) | 一种dvl测速失效时的uuv反步滑模动力定位控制方法 | |
CN114880874B (zh) | 一种水面无人船参数自适应鲁棒估计方法与系统 | |
CN113608534A (zh) | 一种无人艇跟踪控制方法及系统 | |
CN109521798A (zh) | 基于有限时间扩张状态观测器的auv运动控制方法 | |
Garcia-Garcia et al. | Wave filtering for heading control of an AUV based on passive observer | |
CN113608541B (zh) | 一种水下航行器姿态控制方法和系统 | |
Zhao et al. | Nonlinear controller based ADRC for dynamic positioned vessels | |
Maurya et al. | Complementary terrain/single beacon-based AUV navigation | |
Belanger et al. | Submarine autopilot performance optimization with system identification | |
Rodiana et al. | Software and hardware in the loop simulation of navigation system design based on state observer using Kalman filter for autonomous underwater glider | |
Zhang et al. | A Nonlinear Observer for Ship Dynamic Positioning System | |
Lin et al. | Nonlinear passive robust observer design for a DP system based on ACA | |
Xie | Kalman Filter for Dynamic Positioning Ships using position and acceleration feedback | |
CN112904719B (zh) | 一种适用于水下机器人位置环形区域跟踪控制方法 | |
Liu et al. | Dynamic positioning of AUVs in shallow water environment: Observer and controller design | |
An et al. | Dynamic Positioning Observer Design Using Exogenous Kalman Filter |
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 |