0 引　言

在海洋环境中，船舶受到波浪扰动，会产生横摇、纵摇、首摇、升沉、横荡、纵荡6个自由度上的运动^[1]，利用波浪补偿技术对这些运动进行补偿，抵消波浪的影响，可以提高人员换乘过程的效率和安全性。

目前对于在单自由度的升沉补偿已经有深入研究和较广泛的应用。文献[2-5]对升沉补偿装置进行了结构设计和仿真验证；文献[6]建立了电动升沉补偿模型并采用了PID控制方法；文献[7]设计了一种并联补偿机构并采用了模糊PID控制方法；文献[8]提出了一种自适应鲁棒控制器，使系统在参数不确定和存在外部干扰情况下有良好的控制精度和鲁棒性；文献[9]提出了一种基于运动学分析和能量耗散的控制方法，实现了对海上吊机的升沉补偿。

现代船舶一般带有动力定位系统，可以补偿船舶的横摇、纵摇、首摇三自由度运动，因此人员换乘作业时，需要对横荡、纵荡、升沉3个自由度的运动进行补偿。在多自由度波浪补偿技术领域，文献[10]对波浪补偿舷梯进行了运动学建模与仿真；文献[11]对波浪补偿舷梯进行了机械结构设计，采取了PID控制方法；文献[12]对三自由度稳定平台进行了建模并提出了超螺旋滑模控制方法；文献[13]对波浪补偿舷梯进行了运动学建模并提出了一种模型预测控制方案；文献[14]对波浪补偿舷梯进行了运动学和动力学建模，并提出了一种逆动力学控制方案。

本文针对三自由度波浪补偿舷梯建立舷梯的运动学模型、动力学模型以及液压执行器的模型，考虑海浪作用下船舶横荡、纵荡、升沉的三自由度运动对舷梯的未知时变影响和舷梯本身的动态不确定性，构造扩张状态观测器（extended state observer，ESO），并进一步基于舷梯的逆动力学分析，提出含模型辅助的自抗扰控制（active disturbance rejection control，ADRC）作为位置外环控制器，对舷梯的横荡、纵荡、升沉3个自由度进行控制。选择线性自抗扰控制方案进行液压执行器的内环力矩控制，并通过仿真模型验证控制策略的有效性。

1 问题描述与预备知识

波浪补偿舷梯由底座、立柱、梯架、回转机构、俯仰机构、伸缩机构组成，其中底座固定在船舶甲板上，立柱轴线垂直于甲板平面，梯架一端固定在底座上，另一端（舷梯末端）将与目标平台或船舶连接，回转机构、俯仰机构、伸缩机构为3个可动机构。回转机构由液压马达驱动使舷梯绕立柱轴线旋转，俯仰机构由液压缸驱动使梯架绕梯架固定端旋转，伸缩机构由液压马达驱动使梯架第2节平移。波浪补偿舷梯结构如图1所示。

定义坐标系如下：

北东坐标系 $\left\{NED\right\}$ ：其中坐标系原点 ${O}_{NED}$ 一般选为地球表面一固定点， ${x}_{n}$ 轴指向正北方向， ${y}_{n}$ 轴指向正东方向， ${z}_{n}$ 轴方向由右手法则确定，方向向下。北东地坐标系跟随地球自转，进行船舶控制时可以看作一个惯性坐标系；附体坐标系 $\left\{B\right\}$ ：其中坐标原点 ${O}_{b}$ 位于舷梯底座与船舶甲板连接面的中心， ${x}_{b}$ 轴指向船首方向， ${y}_{b}$ 轴指向船舶右舷方向， ${z}_{b}$ 轴方向由右手法则确定，方向垂直船舶甲板向下。附体坐标系中的运动可以表示船舶在海浪影响下的运动。

根据D-H参数法^[15]，将舷梯看作串联机械臂，定义其连杆结构如图2所示。

定义连杆0与附体坐标系重合，定义连杆1的原点为舷梯底座中心， $z$ 轴为回转轴，定义连杆2的原点为梯架固定端， $z$ 轴为俯仰轴，定义连杆3的原点为梯架末端， $z$ 轴为舷梯伸缩方向。

记船舶甲板到立柱顶端距离为 ${L}_{1}$ ，梯架首端点与立柱中轴线距离为 ${L}_{2}$ ，梯架第1段长度为 ${L}_{3}$ ，梯架第2段伸出长度为 ${L}_{4}$ ；记舷梯回转角度为 ${\theta }_{1}$ ，舷梯俯仰角度为 ${\theta }_{2}$ ，定义舷梯系统连杆坐标系如图3所示。

根据Lagrange方程，建立波浪补偿舷梯动力学方程如下式：

$I\left(\lambda \right)\ddot{\lambda }+D\left(\lambda ,\dot{\lambda }\right)\dot{\lambda }+G\left(\lambda \right)=\tau 。$

(1)

式中： $\lambda ={\left[\begin{array}{ccc}{\lambda }_{1}& {\lambda }_{2}& {\lambda }_{3}\end{array}\right]}^{{\rm{T}}}={\left[\begin{array}{ccc}{\theta }_{1}& {\theta }_{2}& {L}_{4}\end{array}\right]}^{{\rm{T}}}$ 表示广义坐标，由回转角 ${\theta }_{1}$ 、俯仰角 ${\theta }_{2}$ 和伸缩长度 ${L}_{4}$ 组成； $\tau ={\left[\begin{array}{ccc}{\tau }_{1}& {\tau }_{2}& {\tau }_{3}\end{array}\right]}^{{\rm{T}}}$ 表示作用在关节上的广义力； $I\left(\lambda \right)$ 为系统惯性矩阵，其表达式如下：

$I\left(\lambda \right)=\left[\begin{array}{ccc}{i}_{11}& 0& 0\\ 0& {i}_{22}& 0\\ 0& 0& {i}_{33}\end{array}\right] 。$

(2)

式中： ${i}_{11}={I}_{z1}+\left({I}_{y2}+{I}_{z3}\right){{c}_{2}}^{2}+({I}_{x2}+{I}_{x3}+{m}_{3}{\left({Z}_{l3}+{\lambda }_{3}\right)}^{2}+ {m}_{2}{{Y}_{l2}}^{2}){{s}_{2}}^{2}$ ， ${i}_{22}={I}_{z2}+{I}_{y3}+{m}_{2}{{Y}_{l2}}^{2}+{m}_{3}{\left({Z}_{l3}+{\lambda }_{3}\right)}^{2}$ ， ${i}_{33}= {m}_{3}$ 。其中 ${I}_{xi}\mathrm{，}{I}_{yi}\mathrm{，}{I}_{zi}$ 为杆 $i$ 对于 $x\mathrm{，}y\mathrm{，}z$ 轴的转动惯量。

$D\left(\lambda ,\dot{\lambda }\right)$ 的形式如下：

$D\left(\lambda ,\dot{\lambda }\right)=\left[\begin{array}{ccc}0& 0& 0\\ 0& {\dot{\lambda }}^{T}{d}_{1}\left(\lambda \right)& 0\\ 0& 0& {\dot{\lambda }}^{T}{d}_{2}\left(\lambda \right)\end{array}\right]。$

(3)

式中： ${d}_{1}\left(\lambda \right) = \mathrm{diag}(2(({I}_{x2} + {I}_{x3}-{I}_{y2}-{I}_{z3}+{m}_{2}{{Y}_{l2}}^{2}+{m}_{3}({Z}_{l3}+ {\lambda }_{3})^{2}){c}_{2}{s}_{2},0,0)$ ， ${d}_{2}\left(\lambda \right)=\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left(2{m}_{3}\left({Z}_{l3}+{\lambda }_{3}\right),\right.2{m}_{3}({Z}_{l3}+ {\lambda }_{3}) {{s}_{2}}^{2},0)$ 。

$G\left(\lambda \right)$ 为系统势能，表达式如下：

$G\left(\lambda \right)=\left[\begin{array}{c}0\\ -\left({m}_{2}{Y}_{l2}+{m}_{3}\left({Z}_{l3}+{\lambda }_{3}\right)\right)g{s}_{2}\\ {m}_{3}g{c}_{2}\end{array}\right]，$

(4)

式中， $g$ 为重力加速度。

根据D-H参数法^[15]对舷梯进行运动学建模，可以得到广义坐标 $\lambda ={\left[\begin{array}{ccc}{\theta }_{1}& {\theta }_{2}& {L}_{4}\end{array}\right]}^{{\rm{T}}}$ 与舷梯端点位置坐标 ${P}_{tip}={\left[\begin{array}{ccc}X& Y& Z\end{array}\right]}^{{\rm{T}}}$ 变换关系如下式：

$\left\{\begin{array}{c}X={L}_{2}{c}_{1}+\left({L}_{3}+{L}_{4}\right){c}_{1}{s}_{2}，\\ Y={L}_{2}{s}_{1}+\left({L}_{3}+{L}_{4}\right){s}_{1}{s}_{2}，\\ Z={L}_{1}-\left({L}_{3}+{L}_{4}\right){c}_{2}。\end{array}\right.$

(5)

$\left\{\begin{array}{l}{\theta }_{1}={\rm{arctan}}\dfrac{Y}{X}，\\ {\theta }_{2}={\rm{arccot}}\dfrac{{L}_{1}-Z}{\dfrac{X}{\mathrm{c}\mathrm{o}\mathrm{s}{\theta }_{1}}-{L}_{2}}，\\ {L}_{4}=\dfrac{X-{L}_{2}\cdot \mathrm{c}\mathrm{o}\mathrm{s}{\theta }_{1}}{\mathrm{c}\mathrm{o}\mathrm{s}{\theta }_{1}\cdot \mathrm{s}\mathrm{i}\mathrm{n}{\theta }_{2}}-{L}_{3}。\end{array}\right.$

(6)

其中， ${c}_{\mathrm{*}}$ 和 ${s}_{\mathrm{*}}$ 指 $\mathrm{c}\mathrm{o}\mathrm{s}\left({\lambda }_{\mathrm{*}}\right)$ 和 $\mathrm{s}\mathrm{i}\mathrm{n}\left({\lambda }_{\mathrm{*}}\right)$ 。

根据舷梯液压传动系统结构，可以得到液压执行器输出力/力矩与关节力/力矩的换算关系：

$\tau =\left[\begin{array}{c}{\tau }_{1}\\ {\tau }_{2}\\ {\tau }_{3}\end{array}\right]=\left[\begin{array}{ccc}26& 0& 0\\ 0& 4\mathrm{c}\mathrm{o}\mathrm{s}\dfrac{{\lambda }_{2}}{2}& 0\\ 0& 0& 15.625\end{array}\right]\left[\begin{array}{c}{q}_{1}\\ {q}_{2}\\ {q}_{3}\end{array}\right] 。$

(7)

其中： ${q}_{1}\mathrm{，}{q}_{2}\mathrm{，}{q}_{3}$ 为液压执行器输出力/力矩。

根据液压系统的流量方程、流量连续性方程和力平衡方程^[16]，建立液压系统传递函数模型：

$\begin{split}{G}_{i}\left(s\right)=&\frac{{q}_{i}}{{x}_{vi}}=\left[\frac{{k}_{qi}{k}_{svi}}{{A}_{pi}}\left({M}_{ti}{{s}_{i}}^{2}+{B}_{pi}{s}_{i}+{K}_{si}\right)\right]\Biggr/\\ &\left[\frac{{V}_{ti}{M}_{ti}}{4{\beta }_{ei}{{A}_{pi}}^{2}}{{s}_{i}}^{3}+\left(\frac{{k}_{cei}{M}_{ti}}{{{A}_{pi}}^{2}}+\frac{{V}_{ti}{B}_{pi}}{4{\beta }_{ei}{{A}_{pi}}^{2}}\right){{s}_{i}}^{2}+\right.\\ &\left.\left(1+\frac{{k}_{cei}{B}_{pi}}{{{A}_{pi}}^{2}}+\frac{{V}_{ti}{K}_{si}}{4{\beta }_{ei}{{A}_{pi}}^{2}}\right){s}_{i}+\frac{{k}_{cei}{K}_{si}}{{{A}_{pi}}^{2}}\right] 。\end{split}$

(8)

其中： $i$ 取1，2，3时，方程8分别为 ${q}_{1}\mathrm{，}{q}_{2}\mathrm{，}{q}_{3}$ 对应的液压系统模型。

假设1　 $I\left(\lambda \right)$ ， $D\left(\lambda ,\dot{\lambda }\right)$ ， $G\left(\lambda \right)$ 均不确定。

假设2　液压系统传递函数 ${G}_{i}\left(s\right)$ 中的参数不确定。

假设3　海浪对舷梯姿态的扰动影响是未知时变的，且该扰动的变化率有界。

其中，波浪补偿舷梯的模型参数难以准确获得，且会随舷梯姿态变化而变化；液压系统的组成器件参数会由于生产情况而在一定范围内随机变化，难以获得准确值；实际情况中的海浪是随机且不可知的，且海浪能量是有限的。因此假设1、假设2和假设3合理。

本文的控制目标为：在假设1、假设2和假设3条件下，设计控制器，补偿舷梯在横荡、纵荡、升沉3个自由度的运动，使舷梯末端位置满足惯性坐标系的期望位置。

2 控制律设计

将波浪补偿舷梯的动态不确定性和海浪扰动视为“总扰动”^[17-18]，同时以舷梯动力学模型进行辅助，设计舷梯位置的自抗扰控制方案；然后将液压系统的参数不确定性视为“总扰动”，设计液压执行器的力和力矩自抗扰控制方案。

2.1 含模型辅助的位置线性自抗扰控制

根据自抗扰控制设计方法，将控制对象变换为积分串联型系统，即 ${y}^{\left(n\right)}=bu$ ^[19]。

其中系统阶次 $n$ 应根据实际物理意义进行选取^[20]。对于波浪补偿舷梯动力学系统，根据方程1可知，系统输入为广义力，系统输出为广义坐标，故取 $n=2$ ，系统阶次为二阶。在舷梯动力学系统方程中，取 ${y}_{1}= {x}_{11}= \lambda$ ， ${x}_{12}=\dot{\lambda }$ ， ${x}_{13}={f}_{1}=\left(I{\left(\lambda \right)}^{-1}-{b}_{1}\right)\tau -\left(I{\left(\lambda \right)}^{-1}D\left(\lambda ,\dot{\lambda }\right)\dot{\lambda }\right.+$ $\left.I{\left(\lambda \right)}^{-1}G\left(\lambda \right)\right)$ ， ${h}_{1}={\dot{f}}_{1}$ ， ${u}_{1}=\tau$ ，则式(1)可以写成：

$\left\{\begin{array}{l}{\dot{x}}_{1}={A}_{1}{x}_{1}+{B}_{1}{u}_{1}+{E}_{1}{h}_{1}，\\ {y}_{1}={C}_{1}{x}_{1}。\end{array}\right.$

(9)

其中： ${A}_{1}=\left[\begin{array}{ccc}O& I& O\\ O& O& I\\ O& O& O\end{array}\right]$ ， ${B}_{1}=\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left(O,{b}_{1},O\right)$ ， ${C}_{1}= \mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left(I,O,O\right)$ ， ${E}_{1}=\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left(O,O,I\right)$ ， $I$ 为三阶单位矩阵， $O$ 为三阶零矩阵， ${b}_{1}$ 是一个三阶对角矩阵， $\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left(\mathrm{*}\right)$ 表示以 $\mathrm{*}$ 为对角线元素的对角矩阵。

对波浪补偿舷梯动力学系统设计二阶线性自抗扰控制器，结构如图4所示。

在经典的自抗扰控制器设计方法中， $b$ 为估计的系统增益，配置为一固定参数。对于舷梯动力学系统， ${b}_{1}$ 会随着舷梯运动而发生变化，记为 ${b}_{1}=b+\Delta b$ ，由式（9）可知， $\Delta bu$ 会被归为“总扰动”中，可能导致扩张状态观测器负担过大以致失效。对经典自抗扰控制器设计方法进行改进，在设计时引入模型辅助，令

${b}_{1}=\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left(\frac{1}{{i}_{11}{{'}}},\frac{1}{{i}_{12}{{'}}},\frac{1}{{i}_{13}{{'}}}\right)，$

(10)

其中： ${i}_{11}\mathrm{{'}}$ ， ${i}_{12}\mathrm{{'}}$ ， ${i}_{13}\mathrm{{'}}$ 为 ${i}_{11}$ ， ${i}_{12}$ ， ${i}_{13}$ 的估计值。

设计扩张状态观测器（ESO）如下：

${\dot{z}}_{1}={A}_{1}{z}_{1}+{B}_{1}{u}_{1}-{L}_{1}\left({z}_{11}-{y}_{1}\right)，$

(11)

其中， ${L}_{1}=\mathrm{d}\mathrm{i}\mathrm{a}\mathrm{g}\left({\beta }_{11}I,{\beta }_{12}I,{\beta }_{13}I\right)$ 。

定义ESO观测误差为 ${e}_{1i}={z}_{1i}-{x}_{1i}$ ，则观测误差状态方程可以表示为：

${\dot{e}}_{1}=\left[\begin{array}{ccc}{\beta }_{11}I& I& O\\ {\beta }_{12}I& O& I\\ {\beta }_{13}I& O& O\end{array}\right]{e}_{1}。$

(12)

通过“带宽法”设计观测器带宽为 ${\omega }_{1o}= [{\omega }_{1o1} \quad {\omega }_{1o2}\quad {\omega }_{1o3}]$ ，可以得到： ${\ \beta }_{11}=3{\omega }_{1o}$ ， ${\ \beta }_{12}=3[{{\omega }_{1o1}}^{2} \quad {{\omega }_{1o2}}^{2} \quad {{\omega }_{1o3}}^{2}]$ ， ${\ \beta }_{13}=[\begin{array}{ccc}{{\omega }_{1o1}}^{3}& {{\omega }_{1o2}}^{3}& {{\omega }_{1o3}}^{3}\end{array}]$ 。

设计线性状态误差反馈控制律（LSEF）如下：

${u}_{1}={{b}_{1}}^{-1}\left[{k}_{1p}\left({v}_{11}-{z}_{11}\right)+{k}_{1d}\left({v}_{12}-{z}_{12}\right)-{z}_{13}\right] 。$

(13)

其中， ${k}_{1p}\mathrm{和}{k}_{1d}$ 为控制器参数。

2.2 液压力矩自抗扰控制

液压系统由3个独立的液压驱动装置组成，以 ${q}_{1}$ 对应的液压马达为例，进行自抗扰控制方案设计， ${q}_{2}$ 和 ${q}_{3}$ 对应部分的设计方法可以类推。

取阀控液压马达系统的阶次为二阶，将简化后的系统方程改写如下：

${k}_{1}{\ddot{F}}_{g}+{k}_{2}{\dot{F}}_{g}+{k}_{3}{F}_{g}=g{\dot{x}}_{v}。$

(14)

其中： ${k}_{1}=\dfrac{{V}_{t}{M}_{t}}{4{\beta }_{e}{{A}_{p}}^{2}}$ ， ${k}_{2}=\dfrac{{k}_{ce}{M}_{t}}{{{A}_{p}}^{2}}$ ， ${k}_{3}=1$ ， $g=\dfrac{{k}_{q}}{{A}_{p}}{M}_{t}$ 。

取 $u = {\dot{x}}_{v}$ ， ${y}_{2} = {x}_{21} = {F}_{g}$ ， ${x}_{22} = {\dot{F}}_{g}$ ， ${x}_{23} = {f}_{2} = -\dfrac{{k}_{3}}{{k}_{1}}{F}_{g}- \dfrac{{k}_{2}}{{k}_{1}}{\dot{F}}_{g}+\left(\dfrac{g}{{k}_{1}}-{b}_{2}\right){\dot{x}}_{v}$ ，则系统方程可以写为：

$\left\{\begin{array}{l}{\dot{x}}_{2}={A}_{2}{x}_{2}+{B}_{2}{u}_{2}+{E}_{2}{h}_{2}，\\ {y}_{2}={C}_{2}{x}_{2}。\end{array}\right.$

(15)

其中： ${A}_{2}=\left[\begin{array}{ccc}0& 1& 0\\ 0& 0& 1\\ 0& 0& 0\end{array}\right]$ ， ${B}_{2}=\left[\begin{array}{c}0\\ {b}_{2}\\ 0\end{array}\right]$ ， ${E}_{2}=\left[\begin{array}{c}0\\ 0\\ 1\end{array}\right]$ ， ${C}_{2}= \left[\begin{array}{ccc}1& 0& 0\end{array}\right]$ ， ${b}_{2}$ 为对系统增益的估计值。

设计扩张状态观测器如下：

${\dot{z}}_{2}={A}_{2}{z}_{2}+{B}_{2}{u}_{2}-{L}_{2}\left({z}_{21}-{y}_{2}\right)，$

(16)

其中， ${L}_{2}={\left[{\beta }_{21},{\beta }_{22},{\beta }_{23}\right]}^{{\rm{T}}}$ 。

定义ESO观测误差为 ${e}_{2i}={x}_{2i}-{z}_{2i}$ ，则观测误差状态方程可以表示为：

${\dot{e}}_{2}=\left[\begin{array}{ccc}{\beta }_{21}& 1& 0\\ {\beta }_{22}& 0& 1\\ {\beta }_{23}& 0& 0\end{array}\right]{e}_{2}。$

(17)

通过“带宽法”设计观测器带宽为 ${\omega }_{2o}=[{\omega }_{2o1} \quad {\omega }_{2o2}\quad{\omega }_{2o3}]$ ，可以得到： ${\ \beta }_{21}=3{\omega }_{2o}$ ， ${\ \beta }_{22}=3[{{\omega }_{2o1}}^{2} \quad {{\omega }_{2o2}}^{2} \quad {{\omega }_{2o3}}^{2}]$ ， ${\ \beta }_{23}=[{{\omega }_{2o1}}^{3}\quad {{\omega }_{2o2}}^{3}\quad {{\omega }_{2o3}}^{3}]$ 。

设计误差反馈控制律如下：

${u}_{2}={{b}_{2}}^{-2}\left[{k}_{2p}\left({v}_{21}-{z}_{21}\right)+{k}_{2d}\left({v}_{22}-{z}_{22}\right)-{z}_{23}\right]。$

(18)

其中， ${k}_{2p}\mathrm{和}{k}_{2d}$ 为控制器参数。

3 仿真研究 3.1 仿真模型与参数

仿真模型中利用Fossen等^[21]开发的Marine System Simulator（MSS）工具箱获得海浪引起的船舶升沉、横荡、纵荡3个自由度上的位移对舷梯模型造成的影响运动，将此运动反向并与期望位置相加作为总参考输入。

仿真采用ITTC海浪谱，选取海浪参数如下：有义波高4 m，平均波向60°，波浪扩散因子3，对应5级海况。

舷梯参数选取如表1所示，设置舷梯初始位置为 $\left[\begin{array}{ccc}{x}_{0}& {y}_{0}& {z}_{0}\end{array}\right]=\left[\begin{array}{ccc}10.5& 0& 6\end{array}\right]$ ，设置舷梯期望位置为 $\left[\begin{array}{ccc}{x}_{d}& {y}_{d}& {z}_{d}\end{array}\right]=\left[\begin{array}{ccc}17& 3& 10\end{array}\right]$ ，控制器参数如表2所示。

得到海浪对舷梯三自由度运动的影响如图5所示。

3.2 仿真比较

分别采用带有重力补偿的PID控制方案和串级线性自抗扰控制方案对舷梯进行控制，得到仿真结果如图6所示。图中曲线分别为舷梯末端在横荡、纵荡、升沉3个方向的位置变化曲线。2种控制方案的动态性能指标对比见表3。

可以看出，串级线性自抗扰控制方案对波浪补偿舷梯有良好的控制效果，可以使舷梯快速地达到期望位置并在海浪环境中稳定保持在目期望位置，并对海浪影响下的船舶运动干扰可以有效抑制。同时串级线性自抗扰控制方案在控制效果上优于PID方案，串级自抗扰控制方案的调节速度、超调量都更小，受到干扰时产生的位置误差也更小，证明所设计的串级自抗扰控制方案对外界未知扰动和系统参数不确定性有着更好的鲁棒性。

4 结　语

1）本文建立主动式波浪补偿舷梯的模型，对其运动学建模、动力学建模和液压系统建模；

2）针对主动式波浪补偿舷梯模型的非线性、模型不确定性等特点，提出了以模型辅助的位置线性自抗扰控制为外环，力矩线性自抗扰控制为内环的串级自抗扰控制方案；

3）根据“带宽法”配置线性自抗扰控制方案的参数，仿真结果表明设计的控制方案对海洋干扰和参数不确定性的有效抑制作用，并验证了其相较于传统PID控制方法更好的控制效果和更强的抗干扰能力。