CN109870914A - 一种机动浮标在随机海浪干扰下的自守卫控制方法 - Google Patents

Abstract

本发明属于海洋浮标领域，具体涉及一种机动浮标在随机海浪干扰下的自守卫控制方法。包括机动浮标自守卫航行时的海浪方向估计、考虑波浪增阻问题的能量最优航路规划。海浪方向估计是在线实时预测机动浮标在自守卫区域内的浪向角，进而得到机动浮标自守卫航行时与海浪的遭遇角；在机动浮标自守卫区域内设置虚拟节点，采用A*算法对机动浮标进行能量最优航路规划；根据机动浮标所受到的波浪增阻与遭遇角之间的函数关系，将机动浮标在规划的航路上航行时的能量损耗作为目标函数，使得机动浮标在规划的航路上航行时的能量损耗最低。本发明可以有效地降低机动浮标在波浪增阻干扰下自守卫航行时的能量损耗，提高机动浮标的航行效率，应用前景广泛。

Description

一种机动浮标在随机海浪干扰下的自守卫控制方法

技术领域

本发明属于海洋浮标领域，具体涉及一种机动浮标在随机海浪干扰下的自守卫控制方法。

背景技术

海洋浮标是以锚定在海上的观测浮标为主体组成的海洋水文水质气象自动观测站。它能按规定要求长期、连续地收集所需海洋水文水质气象资料，特别是能收集到调查船难以恶劣天气及海况的资料。机动浮标是一种特殊的海洋浮标，在受到外界环境影响而漂离到给定的守卫区域外时，它可以依靠自身的动力，自主航行到给定的守卫点，继续执行自守卫测量任务。但是，当机动浮标在近水面自守卫航行时，在随机海浪的干扰下，由于海浪与机动浮标运动之间的相互作用产生的波浪增阻会降低机动浮标的主机推进功率，使得机动浮标在自守卫航行过程中的能量严重损耗，进而降低机动浮标的航行效率。

发明内容

本发明的目的在于提供一种机动浮标在随机海浪干扰下的自守卫控制方法，其能有效减少机动浮标在自守卫航行时的波浪增阻，进而减少能量损耗，提高机动浮标的航行效率。

一种机动浮标在随机海浪干扰下的自守卫控制方法，具体包括以下步骤：

步骤1、根据机动浮标在随机海浪干扰下的纵摇、横摇、垂荡运动的时间序列，采用MAR(MultiAuto Regressive model)模型法计算机动浮标在上述三自由度上运动的交叉谱；

步骤2、根据机动浮标自身的模型参数通过MATLAB中的MSS工具包计算出浮标在上述三自由度上运动的RAO(ResponseAmplitude Operator)数据；

步骤3、根据步骤1中获得的交叉谱和步骤2中获得的RAO数据计算得到机动浮标的海浪方向谱，三者之间的关系为：

式中，φ_ij(ω)为通过机动浮标的时间序列获取的交叉谱，RAO_i(ω,θ)表示在第i个自由度上，波浪入射角为θ，频率为ω时的响应幅值算子，S(ω,θ)为海浪方向谱；

步骤4、根据机动浮标的实际航行情况，在给定的自守卫区域内，设置虚拟节点，每个虚拟节点中都包含此时刻的浪向角信息；

步骤5、根据机动浮标在随机海浪干扰下受到的波浪增阻与浪向角之间的函数关系，将机动浮标在自守卫航行时的能量损耗作为目标函数；

步骤6、采用A*算法对机动浮标在随机海浪干扰下自守卫航行时，从守卫区域边缘点到达给定守卫点进行路径规划，使得机动浮标在规划的航路上航行时能量损耗最低。

所述一种机动浮标在随机海浪干扰下的自守卫控制方法，步骤1中MAR模型建立方法为：

步骤1、假设{Z_k}是k维随机平稳的零均值时间序列，表示如下：

Z_k＝[z₁ z₂ ... z_k]^T

其中，z_i(s·Δt),s＝1,2,...,N，i＝1,2,...,k,其中Δt为采样时间，令采样时间Δt＝1，则z_i(s·Δt)＝z_i(s)；

步骤2、由于所需要分析的是机动浮标三自由度的运动时间序列，在每个时间序列中，其前p个时间序列值对第s个时间序列进行建模，则有：

Z(s)-A₁Z(s-1)-...-A_p(s-p)＝ε(s)

其中，A_i为k×k维的权重系数矩阵；ε(s)＝[ε₁ ε₂ ... ε_k]^T表示均值为0，方差为σ_i的高斯白噪声序列，由于ε_i与z(s)的取值无关，因此ε(s)与ε(s-p)的协方差为零；由此可得白噪声向量的协方差矩阵如下：

式中，σ_ij表示ε_i与ε_j之间的协方差；将上式改写如下：

令：V(s)＝[z(s-1) z(s-2) ... z(s-p)]^T

Λ_p＝[Λ₁ Λ₂ ... Λ_p](k×(k·p))

则MAR模型的最终形式为：

Z(s)＝Λ_pV(s)+ε(s)

所述一种机动浮标在随机海浪干扰下的自守卫控制方法，对上述权重系数矩阵A_i求解时采用Whittle递推法；由于MAR模型的权值系数A_j(j＝1,2,...,p)随着模型阶数p的变化而改变；因此，令：Λ_pj＝(A_p1 A_p1 ... A_pj)；j＝1,2,...,p，R_i(i＝1,2,...,p)为相关矩阵，对于不同迟滞l相关矩阵中的第i,j个元素的计算公式表示为：

因此，递推公式为：

针对上述递推法求出的P项权重系数矩阵采用最优阶数判断准则来选取最有效的一组权重系数矩阵；该准则的描述为：对于长度为N的k维平稳随机序列{Z_k}，其某一步预报误差方差矩阵的估计量用D_p来表示，最终的预报误差为其行列式的值：

所述一种机动浮标在随机海浪干扰下的自守卫控制方法，步骤1中计算机动浮标运动交叉谱的方法包括以下步骤：

根据上述建立的MAR模型，选择适当的阶数，就可以准确地使用前P个历史值准确预测所求摇荡时序的当前值。

步骤1、由交叉谱密度的基本定义得到白噪声的交叉谱密度

式中f为变换后的频率，由于当且仅当s＝t时，即等式右边只是将σ_st进行傅里叶变换，而σ_st是一个常数项，对于常数项的傅里叶变换仍然是其本身；可得

步骤2、考虑残差的影响，令A_ij(0)＝-δ_ij，其中：

根据上式可得：

其中，ε_i(s)是σ_j(s)的线性变换，模型当维数退化到一维时，A_ij(p)为脉冲响应函数：

此时，原始信号经过脉冲响应函数的线性变化后，在其g频点的功率谱密度p_HH(g)为：

p_HH(g)＝|B(g)|²p_yy(g),-∞＜g＜∞

式中，p_yy(g)是y(t)的功率谱密度，为脉冲响应函数h(τ)的频率响应函数；

分析可知，多维与一维的情况相同，但是多维误差的交叉谱为复数，可知：

其中，α,β＝1,2,...k；p(f)即为所求的交叉谱；考虑到其矩阵形式表示为：

Π＝A(f)p(f)A(f)^H

式中，Π为由σ_st组成的矩阵；假定矩阵A(f)为非奇异矩阵，上式两端同时乘以A(f)的逆矩阵和A(f)^H(即A(f)的共轭转置矩阵)，可得：

p(f)＝A(f)^-1ΠA(f)^-H

综上，便可获得机动浮标三自由度摇荡运动的交叉谱，最终可得海浪方向谱。

所述一种机动浮标在随机海浪干扰下的自守卫控制方法，步骤4中虚拟节点的设置原则为：根据机动浮标在湖上和海上试验的实际情况，自守卫区域半径一般取50m左右，守卫内径在5m左右，而机动浮标自身的长度为1.35m。因此，本发明设置的虚拟节点是以自守卫内圆为中心点，在守卫区域内以5m的间距均匀设置虚拟节点。每个虚拟节点都包含此时的浪向角信息。

所述一种机动浮标在随机海浪干扰下的自守卫控制方法，机动浮标在随机海浪干扰下产生的波浪增阻f(β)与浪向即遭遇角β之间的函数关系为：

f(β)＝p₀+p₁β+...+p_kβ_i ^k

式中，(p₀,...,p_k)表示所需要求出的多项式系数，在一级海况、二级海况、三级海况下，采用最小二乘法对试验数据进行拟合仿真，求得的系数矩阵分别为：p＝[2.73×10^-10-1.22×10^-7 1.85×10^-4 -7.44×10^-4]^T，p＝[-1.24×10^-9 8.85×10^-7 -2.22×10^-40.027]^T，p＝[-1.13×10^-9 1.53×10^-6 -4.92×10^-4 0.07]^T。

所述一种机动浮标在随机海浪干扰下的自守卫控制方法，步骤6具体包括以下步骤：

步骤6.1、设计代价函数，A*算法中的代价函数f(x)表示为：

f(x)＝g(x)+h(x)

式中，f(x)表示从起始节点N₀开始到终止节点的最佳路径的代价总值；g(x)表示从起始节点N₀开始到某一个节点x时的最佳路径代价值，h(x)表示从某一节点x到终止节点的实际最佳路径代价值，h(x)又被称为启发函数，当h(x)的值小于某一节点x到终止节点的实际距离时，可以得到最优解；

步骤6.2、根据波浪增阻与遭遇角之间的函数关系表达式，可得机动浮标航行至某一节点处的能量消耗值，设机动浮标的航速为v，起始节点的位置为(x₀,y₀)，在某一节点N_i(x_i,y_i)节点处，航行至下一节点N_i+1(x_i+1,y_i+1)的时间为t，则在这两个相邻节点的路径长度为l，则：

机动浮标在这段路径上航行的能量损失ΔP_loss(N_i,N_i+1)表示为：

把在整个规划的路径上损耗的总能量P_loss作为代价函数，表示为：

步骤6.3、设计启发函数，选择欧几里得距离法作为机动浮标路径规划的启发方法，其对应的函数表达式表示为：

所述一种机动浮标在随机海浪干扰下的自守卫控制方法，步骤6中A*算法具体包括以下步骤：

步骤1、建立三个表，分别为open表、closed表、parent表，其中，open表用来存放当前根据海流和海浪方向信息规划的节点信息、closed表用来存放已经选取的节点信息、parent表用来存放父节点信息；

步骤2、对算法进行初始化设置，即将open表、closed表和parent表清空，然后设置起始节点N₀、终止节点N_n的位置信息；

步骤3、将起始节点N₀放入open表中，令g(N₀)＝0，则f＝h；

步骤4、判断open表是否为空，若open表为空，则所需规划的路径不存在，算法结束，若open表不为空，则执行步骤5；

步骤5、把open表中的代价函数值最小的节点定义为N_best，判断N_best是否为终止节点N_n，若N_best＝N_n，则表示当前路径规划完成，算法结束，如果N_best≠N_n，则执行步骤6；

步骤6、计算open表中所存入的节点的代价函数值P_loss，并比较相对应的代价函数值的大小，把代价函数值最小的节点即N_best从open表中赋给closed表中；

步骤7、判断是否满足能量损耗最低的约束条件，若满足在该节点的能量损耗最低，则计算该节点的代价值，并与实际代价值进行比较，将代价值最小的节点添加到open表和parent表中，返回步骤6，若不满足能量损耗最低的约束条件，则返回步骤4。

本发明的有益效果在于：

本发明首先根据机动浮标在随机海浪干扰下的纵摇、横摇、垂荡运动的时间序列采用MAR模型法获取机动浮标相对应运动的交叉谱φ_ij(ω)，使用MATLAB获取机动浮标相应运动的RAO数据，结合交叉谱和RAO数据得到海浪方向谱S(ω,θ)，即可得到机动浮标自守卫航行时的实时浪向信息。然后根据机动浮标自守卫航行时的实际情况，合理设置虚拟节点，针对波浪增阻问题，通过设计A*算法对机动浮标自守卫航行进行能量最优航路规划。本发明可以有效减少机动浮标自守卫航行时的波浪增阻，进而减少能量损耗，提高航行效率。

附图说明

图1为本发明的机动浮标自守卫区域示意图；

图2为本发明的机动浮标的结构示意图；

图3为本发明的机动浮标在随机海浪干扰下的自守卫航行控制方法流程图；

图4为本发明的虚拟节点设置的示意图；

图5为本发明的A*算法的设计流程图。

具体实施方式

下面结合附图对本发明做进一步描述：

如附图1所示，为本发明的机动浮标自守卫区域示意图，如附图2所示，为本发明的机动浮标的结构示意图；机动浮标停机在守卫内圆中进行测量任务时，在外界未知环境干扰下，会逐渐漂流到守卫外圆之外。此时，机动浮标开启自守卫任务，即测量并计算从当前位置到达守卫点的目标航向与机动浮标艏向的差值，启动推进装置，包括主推进器1与辅助推进器2，通过操纵方向舵3改变航向进行自守卫航行，直至到达守卫内圆中，然后停机继续进行测量任务。

如附图3所示，为本发明的机动浮标在随机海浪干扰下的自守卫航行控制方法流程图，主要包括以下内容：

步骤1、根据机动浮标在随机海浪干扰下的纵摇、横摇、垂荡运动的时间序列，采用MAR(MultiAuto Regressive model)模型法计算机动浮标在上述三自由度上运动的交叉谱；

步骤2、根据机动浮标自身的模型参数通过MATLAB中的MSS工具包计算出浮标在上述三自由度上运动的RAO(ResponseAmplitude Operator)数据；

步骤3、根据步骤1中获得的交叉谱和步骤2中获得的RAO数据计算得到机动浮标的海浪方向谱，三者之间的关系为：

式中，φ_ij(ω)为通过机动浮标的时间序列获取的交叉谱；RAO_i(ω,θ)表示在第i个自由度上，波浪入射角为θ，频率为ω时的响应幅值算子，S(ω,θ)为海浪方向谱；

步骤4、根据机动浮标的实际航行情况，在给定的自守卫区域内，设置虚拟节点，每个虚拟节点中都包含此时刻的浪向角信息；

步骤5、根据机动浮标在随机海浪干扰下受到的波浪增阻与浪向角之间的函数关系，将机动浮标在自守卫航行时的能量损耗作为目标函数；

步骤6、采用A*算法对机动浮标在随机海浪干扰下自守卫航行时，从守卫区域边缘点到达给定守卫点进行路径规划，使得机动浮标在规划的航路上航行时能量损耗最低。

所述一种机动浮标在随机海浪干扰下的自守卫控制方法，步骤1中MAR模型建立方法为：

步骤1、假设{Z_k}是k维随机平稳的零均值时间序列，表示如下：

Z_k＝[z₁ z₂ ... z_k]^T

其中，z_i(s·Δt),s＝1,2,...,N，i＝1,2,...,k,其中Δt为采样时间，令采样时间Δt＝1，则z_i(s·Δt)＝z_i(s)；

步骤2、由于所需要分析的是机动浮标三自由度的运动时间序列，在每个时间序列中，其前p个时间序列值对第s个时间序列进行建模，则有：

Z(s)-A₁Z(s-1)-...-A_p(s-p)＝ε(s)

其中，A_i为k×k维的权重系数矩阵；ε(s)＝[ε₁ ε₂ ... ε_k]^T表示均值为0，方差为σ_i的高斯白噪声序列，由于ε_i与z(s)的取值无关，因此ε(s)与ε(s-p)的协方差为零；由此可得白噪声向量的协方差矩阵如下：

式中，σ_ij表示ε_i与ε_j之间的协方差；将上式改写如下：

令：V(s)＝[z(s-1) z(s-2) ... z(s-p)]^T

Λ_p＝[Λ₁ Λ₂ ... Λ_p](k×(k·p))

则MAR模型的最终形式为：

Z(s)＝Λ_pV(s)+ε(s)

所述一种机动浮标在随机海浪干扰下的自守卫控制方法，对上述权重系数矩阵求解时采用Whittle递推法；由于MAR模型的权值系数A_j(j＝1,2,...,p)随着模型阶数p的变化而改变；因此，令：Λ_pj＝(A_p1 A_p1 ... A_pj)；j＝1,2,...,p，R_i(i＝1,2,...,p)为相关矩阵，对于不同迟滞l相关矩阵中的第i,j个元素的计算公式表示为：

因此，递推公式为：

针对上述递推法求出的P项权重系数矩阵采用最优阶数判断准则来选取最有效的一组权重系数矩阵；该准则的描述为：对于长度为N的k维平稳随机序列{Z_k}，其某一步预报误差方差矩阵的估计量用D_p来表示，最终的预报误差为其行列式的值：

所述一种机动浮标在随机海浪干扰下的自守卫控制方法，步骤1中计算机动浮标运动交叉谱的方法包括以下步骤：

根据上述建立的MAR模型，选择适当的阶数，就可以准确地使用前P个历史值准确预测所求摇荡时序的当前值。

步骤1、由交叉谱密度的基本定义得到白噪声的交叉谱密度

式中f为变换后的频率，由于当且仅当s＝t时，即等式右边只是将σ_st进行傅里叶变换，而σ_st是一个常数项，对于常数项的傅里叶变换仍然是其本身；可得

步骤2、考虑残差的影响，令A_ij(0)＝-δ_ij，其中：

根据上式可得：

其中，ε_i(s)是σ_j(s)的线性变换，模型当维数退化到一维时，A_ij(p)为脉冲响应函数：

此时，原始信号经过脉冲响应函数的线性变化后，在其g频点的功率谱密度p_HH(g)为：

p_HH(g)＝|B(g)|²p_yy(g),-∞＜g＜∞

式中，p_yy(g)是y(t)的功率谱密度，为脉冲响应函数h(τ)的频率响应函数；

分析可知，多维与一维的情况相同，但是多维误差的交叉谱为复数，可知：

其中，α,β＝1,2,...k；p(f)即为所求的交叉谱；考虑到其矩阵形式表示为：

Π＝A(f)p(f)A(f)^H

式中，Π为由σ_st组成的矩阵；假定矩阵A(f)为非奇异矩阵，上式两端同时乘以A(f)的逆矩阵和A(f)^H(即A(f)的共轭转置矩阵)，可得：

p(f)＝A(f)^-1ΠA(f)^-H

综上，便可获得机动浮标三自由度摇荡运动的交叉谱，最终可得海浪方向谱。

所述一种机动浮标在随机海浪干扰下的自守卫控制方法，机动浮标在随机海浪干扰下产生的波浪增阻f(β)与浪向即遭遇角β之间的函数关系为：

f(β)＝p₀+p₁β+...+p_kβ_i ^k

式中，(p₀,...,p_k)表示所需要求出的多项式系数，在一级海况、二级海况、三级海况下，采用最小二乘法对试验数据进行拟合仿真，求得的系数矩阵分别为：p＝[2.73×10^-10-1.22×10^-7 1.85×10^-4 -7.44×10^-4]^T，p＝[-1.24×10^-9 8.85×10^-7 -2.22×10^-40.027]^T，p＝[-1.13×10^-9 1.53×10^-6 -4.92×10^-4 0.07]^T。

所述一种机动浮标在随机海浪干扰下的自守卫控制方法，步骤6具体包括以下步骤：

步骤6.1、设计代价函数，A*算法中的代价函数f(x)表示为：

f(x)＝g(x)+h(x)

式中，f(x)表示从起始节点N₀开始到终止节点的最佳路径的代价总值；g(x)表示从起始节点N₀开始到某一个节点x时的最佳路径代价值，h(x)表示从某一节点x到终止节点的实际最佳路径代价值，h(x)又被称为启发函数，当h(x)的值小于某一节点x到终止节点的实际距离时，可以得到最优解；

步骤6.2、根据波浪增阻与遭遇角之间的函数关系表达式，可得机动浮标航行至某一节点处的能量消耗值，设机动浮标的航速为v，起始节点的位置为(x₀,y₀)，在某一节点N_i(x_i,y_i)节点处，航行至下一节点N_i+1(x_i+1,y_i+1)的时间为t，则在这两个相邻节点的路径长度为l，则：

机动浮标在这段路径上航行的能量损失ΔP_loss(N_i,N_i+1)表示为：

把在整个规划的路径上损耗的总能量P_loss作为代价函数，表示为：

步骤6.3、设计启发函数，选择欧几里得距离法作为机动浮标路径规划的启发方法，其对应的函数表达式表示为：

如附图4所示，为本发明用于机动浮标在随机海浪干扰下的自守卫航行控制方法的虚拟节点设置示意图，根据机动浮标在湖上和海上试验的实际情况，自守卫区域半径一般取50m左右，守卫内径在5m左右，而机动浮标自身的长度为1.35m。因此，本发明设置的虚拟节点是以自守卫内圆为中心点，在守卫区域内以5m的间距均匀设置虚拟节点。每个虚拟节点都包含此时的浪向角信息。其中，规定机动浮标开始自守卫航行时的点为起始节点N₀，到达设置的守卫目标点为终止节点N_n。

如附图5所示，为本发明用于机动浮标在随机海浪干扰下的自守卫航行控制方法的A*算法设计流程图，对其设计步骤详细描述如下：

本发明对于虚拟节点的设置与选取是依据机动浮标在随机海浪干扰下，产生的波浪增阻最小为准则；当机动浮标航行至某一节点时，其艏向与海浪方向之间的夹角即遭遇角的变化，影响波浪增阻的变化；通过航向控制器调整艏向来改变遭遇角的大小，进而使波浪增阻减少。此时选取多个虚拟节点，通过计算这些到达节点处的波浪增阻大小，通过比较来判断选取最优节点，进而完成能量最优路径规划。

步骤1、建立三个表，分别为open表、closed表、parent表，其中，open表用来存放当前根据海流和海浪方向信息规划的节点信息、closed表用来存放已经选取的节点信息、parent表用来存放父节点信息；

步骤2、对算法进行初始化设置，即将open表、closed表和parent表清空，然后设置起始节点N₀、终止节点N_n的位置信息；

步骤3、将起始节点N₀放入open表中，令g(N₀)＝0，则f＝h；

步骤4、判断open表是否为空，若open表为空，则所需规划的路径不存在，算法结束，若open表不为空，则执行步骤5；

步骤5、把open表中的代价函数值最小的节点定义为N_best，判断N_best是否为终止节点N_n，若N_best＝N_n，则表示当前路径规划完成，算法结束，如果N_best≠N_n，则执行步骤6；

步骤6、计算open表中所存入的节点的代价函数值P_loss，并比较相对应的代价函数值的大小，把代价函数值最小的节点即N_best从open表中赋给closed表中；

步骤7、判断是否满足能量损耗最低的约束条件，若满足在该节点的能量损耗最低，则计算该节点的代价值，并与实际代价值进行比较，将代价值最小的节点添加到open表和parent表中，返回步骤6，若不满足能量损耗最低的约束条件，则返回步骤4。

Claims (8)

1.一种机动浮标在随机海浪干扰下的自守卫控制方法，其特征在于，具体包括以下步骤：

步骤1、根据机动浮标在随机海浪干扰下的纵摇、横摇、垂荡运动的时间序列，采用MAR模型法计算机动浮标在纵摇、横摇、垂荡三自由度上运动的交叉谱；

步骤2、根据机动浮标自身的模型参数通过MATLAB中的MSS工具包计算出浮标在上述三自由度上运动的RAO数据；

步骤3、根据步骤1中获得的交叉谱和步骤2中获得的RAO数据计算得到机动浮标的海浪方向谱，三者之间的关系为：

式中，φ_ij(ω)为通过机动浮标的时间序列获取的交叉谱；RAO_i(ω,θ)表示在第i个自由度上，波浪入射角为θ，频率为ω时的响应幅值算子，S(ω,θ)为海浪方向谱；

步骤4、根据机动浮标的实际航行情况，在给定的自守卫区域内，设置虚拟节点，每个虚拟节点中都包含此时刻的浪向角信息；

步骤5、根据机动浮标在随机海浪干扰下受到的波浪增阻与浪向角之间的函数关系，将机动浮标在自守卫航行时的能量损耗作为目标函数；

步骤6、采用A*算法对机动浮标在随机海浪干扰下自守卫航行时，从守卫区域边缘点到达给定守卫点进行路径规划，使得机动浮标在规划的航路上航行时能量损耗最低。

2.根据权利要求1所述一种机动浮标在随机海浪干扰下的自守卫控制方法，其特征在于，步骤1中所述MAR模型建立方法为：

步骤1、设{Z_k}是k维随机平稳的零均值时间序列，表示为：

Z_k＝[z₁ z₂...z_k]^T

其中，z_i(s·Δt),s＝1,2,...,N，i＝1,2,...,k,其中Δt为采样时间，令采样时间Δt＝1，则z_i(s·Δt)＝z_i(s)；

步骤2、由于所需要分析的是机动浮标三自由度的运动时间序列，在每个时间序列中，其前p个时间序列值对第s个时间序列进行建模，则有：

Z(s)-A₁Z(s-1)-...-A_p(s-p)＝ε(s)

其中，A_i为k×k维的权重系数矩阵，ε(s)＝[ε₁ ε₂...ε_k]^T表示均值为0，方差为σ_i的高斯白噪声序列，由于ε_i与z(s)的取值无关，因此ε(s)与ε(s-p)的协方差为零；由此可得白噪声向量的协方差矩阵如下：

式中，σ_ij表示ε_i与ε_j之间的协方差；将上式改写为：

令：V(s)＝[z(s-1) z(s-2)...z(s-p)]^T

Λ_p＝[Λ₁ Λ₂...Λ_p](k×(k·p))

则MAR模型的最终形式为：

Z(s)＝Λ_pV(s)+ε(s)。

3.根据权利要求1或权利要求2所述一种机动浮标在随机海浪干扰下的自守卫控制方法，其特征在于：所述权重系数矩阵A_i求解时采用Whittle递推法；由于MAR模型的权值系数A_j(j＝1,2,...,p)随着模型阶数p的变化而改变；因此，令：Λ_pj＝(A_p1 A_p1...A_pj)；j＝1,2,...,p，R_i(i＝1,2,...,p)为相关矩阵，对于不同迟滞l相关矩阵中的第i,j个元素的计算公式表示为：

因此，递推公式为：

针对上述递推法求出的P项权重系数矩阵采用最优阶数判断准则来选取最有效的一组权重系数矩阵；该准则的描述为：对于长度为N的k维平稳随机序列{Z_k}，其某一步预报误差方差矩阵的估计量用D_p来表示，最终的预报误差为其行列式的值：

4.根据权利要求1所述一种机动浮标在随机海浪干扰下的自守卫控制方法，其特征在于，步骤1中所述计算机动浮标运动交叉谱的方法包括以下步骤：

步骤1、由交叉谱密度的定义得到白噪声的交叉谱密度

式中f为变换后的频率，由于当且仅当s＝t时，即等式右边只是将σ_st进行傅里叶变换，而σ_st是一个常数项，对于常数项的傅里叶变换仍然是其本身；可得

步骤2、考虑残差的影响，令A_ij(0)＝-δ_ij，其中：

根据上式可得：

其中，ε_i(s)是σ_j(s)的线性变换，模型当维数退化到一维时，A_ij(p)为脉冲响应函数：

此时，原始信号经过脉冲响应函数的线性变化后，在其g频点的功率谱密度p_HH(g)为：

pHH(_g)＝|B(g)|²p_yy(g),-∞＜g＜∞

式中，p_yy(g)是y(t)的功率谱密度，为脉冲响应函数h(τ)的频率响应函数；可知多维与一维的情况相同，但是多维误差的交叉谱为复数，可知：

其中，p(f)即为所求的交叉谱；考虑到其矩阵形式表示为：

Π＝A(f)p(f)A(f)^H

式中，Π为由σ_st组成的矩阵；假定矩阵A(f)为非奇异矩阵，上式两端同时乘以A(f)的逆矩阵和A(f)^H，可得：

p(f)＝A(f)^-1ΠA(f)^-H。

5.根据权利要求1所述一种机动浮标在随机海浪干扰下的自守卫控制方法，其特征在于，步骤4中所述虚拟节点的设置原则为：以自守卫内圆为中心点，在守卫区域内以5m的间距均匀设置虚拟节点，每个虚拟节点都包含此时的浪向角信息。

6.根据权利要求1所述一种机动浮标在随机海浪干扰下的自守卫控制方法，其特征在于，

步骤5所述机动浮标在随机海浪干扰下产生的波浪增阻f(β)与浪向即遭遇角β之间的函数关系为：

f(β)＝p₀+p₁β+...+p_kβ_i ^k

式中，(p₀,...,p_k)表示所需要求出的多项式系数，在一级海况、二级海况、三级海况下，采用最小二乘法对试验数据进行拟合仿真，求得的系数矩阵分别为：p＝[2.73×10^-10 -1.22×10^-7 1.85×10^-4 -7.44×10^-4]^T，p＝[-1.24×10^-9 8.85×10^-7 -2.22×10^-4 0.027]^T，p＝[-1.13×10^-9 1.53×10^-6 -4.92×10^-4 0.07]^T。

7.根据权利要求1所述一种机动浮标在随机海浪干扰下的自守卫控制方法，其特征在于，所述步骤6具体包括以下步骤：

步骤6.1、设计代价函数，A*算法中的代价函数f(x)表示为：

f(x)＝g(x)+h(x)

式中，f(x)表示从起始节点N₀开始到终止节点的最佳路径的代价总值；g(x)表示从起始节点N₀开始到某一个节点x时的最佳路径代价值，h(x)表示从某一节点x到终止节点的实际最佳路径代价值，h(x)又被称为启发函数，当h(x)的值小于某一节点x到终止节点的实际距离时，得到最优解；

步骤6.2、根据波浪增阻与遭遇角之间的函数关系表达式，可得机动浮标航行至某一节点处的能量消耗值，设机动浮标的航速为v，起始节点的位置为(x₀,y₀)，在某一节点N_i(x_i,y_i)节点处，航行至下一节点N_i+1(x_i+1,y_i+1)的时间为t，则在这两个相邻节点的路径长度为l，则：

机动浮标在这段路径上航行的能量损失ΔP_loss(N_i,N_i+1)表示为：

把在整个规划的路径上损耗的总能量P_loss作为代价函数，表示为：

步骤6.3、设计启发函数，选择欧几里得距离法作为机动浮标路径规划的启发方法，其对应的函数表达式表示为：

8.根据权利要求1所述一种机动浮标在随机海浪干扰下的自守卫控制方法，其特征在于，步骤6中所述A*算法具体包括以下步骤：

步骤1、建立三个表，分别为open表、closed表、parent表，其中，open表用来存放当前根据海流和海浪方向信息规划的节点信息、closed表用来存放已经选取的节点信息、parent表用来存放父节点信息；

步骤2、对算法进行初始化设置，即将open表、closed表和parent表清空，然后设置起始节点N₀、终止节点N_n的位置信息；

步骤3、将起始节点N₀放入open表中，令g(N₀)＝0，则f＝h；

步骤4、判断open表是否为空，若open表为空，则所需规划的路径不存在，算法结束，若open表不为空，则执行步骤5；

步骤5、把open表中的代价函数值最小的节点定义为N_best，判断N_best是否为终止节点N_n，若N_best＝N_n，则表示当前路径规划完成，算法结束，如果N_best≠N_n，则执行步骤6；

步骤6、计算open表中所存入的节点的代价函数值P_loss，并比较相对应的代价函数值的大小，把代价函数值最小的节点即N_best从open表中赋给closed表中；

步骤7、判断是否满足能量损耗最低的约束条件，若满足在该节点的能量损耗最低，则计算该节点的代价值，并与实际代价值进行比较，将代价值最小的节点添加到open表和parent表中，返回步骤6，若不满足能量损耗最低的约束条件，则返回步骤4。