0 引　言

水下无人航行器（UUV）自主、智能、灵活且能长时间工作需要具备良好的控制性能，。但是因为水下环境复杂，水流多变，干扰不可控制，航向稳定性及其精度无法得到保证仍是一个需要解决的问题。

目前实现UUV航向控制的主要方法包括PID控制、反步法控制、模糊控制及滑模控制等^[1]。PID控制算法较为简便且控制性能较强，在航向控制方面受到广泛应用。但对于模型复杂的控制对象而言，PID控制无法满足其设计需求，通常需要与其他的控制方法相结合来实现较复杂系统的控制，这会使得其设计变得更为复杂^[2]。反步法控制的本质是先拆分，再集成，其主要思想是将高复杂度、高阶数的非线性系统逐一拆分成若干个低阶数子系统，再根据李雅普诺夫稳定性原理对每个子系统分别设计模拟控制量，再将各部分子系统集成，最终实现总体系统控制^[3]。该方法针对模型较为简单的系统设计精确，而且当系统存在不确定项时，可以放宽匹配条件。但是设计过程较为繁琐，在航向控制问题中，系统控制量求解时随着系统阶数的增加，微分计算越来越复杂，较难在实际控制中应用。模糊控制方法是一种借助模糊规则和论域分割来实现将复杂系统简单化的控制方法，其不完全依赖于系统数学模型，从而具有较强鲁棒性^[4]，但由于模糊控制自身局限性，难以实现复杂系统的稳定性分析及整体化设计，使得无法精确控制系统，使得其动态品质变差。

滑模控制相较于其他控制方法而言，其系统结构可以根据当前的状态有针对性地不断优化，且其对外界干扰及模型参数变化不敏感，从而表现出较强的自适应能力和鲁棒性^[5]。本文采用滑模控制方法，在模糊运动学模型的基础上，采用干扰观测器模拟外界干扰，对系统进行干扰补偿，结合连续滑模控制器实现水下无人航行器的航向跟踪。

1 水下无人航行器建模 1.1 坐标系构建

将UUV的空间运动分解成两部分，第1部分忽略其形变将其看成在空间中运动的质点，第2部分为UUA自身绕重心定点转动。此次建模主要应用了地面坐标系^[6]（固定坐标系）和潜航器坐标系（运动坐标系）对其空间运动进行分析。

1.1.1 地面坐标系（固定坐标系） $\boldsymbol{E}-\mathit{\xi }\mathit{\eta }\mathit{\zeta }$

为了分析UUV在整个空间中的运动，首先建立一个地面坐标系。在地面适合处选择固定点E将其设置为坐标原点，ζ轴垂直于水平面向下，称为铅垂轴；ξ轴平行于水平面，且UUV前进方向为正方向；η轴垂直于另外两轴，整个坐标轴满足右手系，如图1所示。

1.1.2 UUV坐标系（运动坐标系）O-xyz

为了分析UUV在水下的俯仰、首摇、转动问题，建立一个UUV坐标系^[7]（运动坐标系）。坐标原点设置在UUV的重心处，x轴水平于UUV艇身，前进方向为正方向；z轴垂直x轴，下为正方向；y轴与另两轴垂直且整个坐标系满足右手系。

1.2 运动参数 1.2.1 姿态角

偏航角 $\varPsi$ ：坐标系的x轴在水平面上的投影与地面坐标系ξ轴的夹角。

俯仰角 $ \theta $ ：x轴与水平面的夹角。

横滚角 $ \phi $ ：z轴与地面坐标系ζ轴的夹角。

1.2.2 坐标系之间的转化

一般UUV运动问题需要用到2个坐标系之间相互转换，适用于航向追踪问题的坐标系转换满足方程^[8]

$ \left[\begin{array}{c}\xi \\ \eta \\ \zeta \end{array}\right] ={{T}} \left[\begin{array}{c}x\\ y\\ z\end{array}\right], $

(1)

$\begin{split} {\rm{T}}=&\left[\begin{array}{ccc}{\rm{cos}}\varPsi {\rm{cos}}\theta & {\rm{cos}}\varPsi {\rm{sin}}\theta {\rm{sin}}\phi -{\rm{sin}}\varPsi {\rm{cos}}\varPsi \\ {\rm{sin}}\varPsi {\rm{cos}}\theta & {\rm{sin}}\varPsi {\rm{sin}}\theta {\rm{sin}}\phi +{\rm{cos}}\varPsi {\rm{cos}}\phi\\ -{\rm{sin}}\theta & {\rm{cos}}\theta {\rm{sin}}\phi \end{array}\right. \\ & \left.\begin{array}{ccc}{\rm{cos}}\varPsi sin\theta {\rm{cos}}\varPsi +{\rm{sin}}\varPsi {\rm{sin}}\phi \\ {\rm{sin}}\varPsi {\rm{sin}}\theta {\rm{cos}}\phi -{\rm{cos}}\varPsi {\rm{sin}}\phi \\ {\rm{cos}}\theta{\rm{cos}}\phi \end{array}\right] 。\end{split}$

(2)

式中：T=[n s a]= $ \left[\begin{array}{ccc}{r}_{11}&{r}_{12}&{r}_{13}\\ {r}_{21}&{r}_{22}&{r}_{23}\\ {r}_{31}&{r}_{32}&{r}_{33}\end{array}\right] $ ，满足正交矩阵且行列式值为1。

1.3 数学模型

在UUV系统中，依然满足牛顿第二定律。忽略UUV的形变将其看做一个刚体，并对其整体进行分析，满足基础动力学方程：

$ \sum _{i}{m}_{i}\frac{{{\rm{d}}v}_{i}}{{\rm{d}}t} = \sum _{t}{F}_{i}。$

(3)

式中： ${{m}}_{{i}}$ 为刚体质量； ${{v}}_{{i}}$ 为表示刚体质心的速度矢量； ${{F}}_{{i}}$ 为刚体所受的外力。

1.3.1 运动学模型（固定坐标系下）

UUV行器通常有6个自由度，并且它的运动可以看做是刚体运动^[9]。在UUV运动过程中，令 $ v=\left[{v}_{x}0 0 \right]T $ 表示潜航器只有沿着x轴的线速度，其他方向线速度为0， $ p=[x,y,z]T $ 。令 $ \omega =\left[{\omega }_{x},{\omega }_{y},{\omega }_{z}\right]T $ 为沿x，y，z轴方向的角速度。推导可得 $ \dot{p}={J}_{1}\left(\eta \right){v}_{x},\dot{\eta }= $ ${J}_{2}\left(\eta \right)\omega 。$

将上式转化为广义向量的矩阵形式为：

$\begin{split} \left[\begin{array}{c}\dot{x}\left(t\right)\\ \dot{y}\left(t\right)\\ \dot{z}\left(t\right)\\ \dot{\varphi }\left(t\right)\\ \dot{\theta }\left(t\right)\\ \dot{\varPsi }\left(t\right)\end{array}\right]= &\left[\begin{array}{c}{\rm{cos}}\varPsi {\rm{cos}}\theta \\ {\rm{sin}}\varPsi {\rm{cos}}\theta \\ -{\rm{sin}}\theta \\ 0\\ 0\\ 0\end{array}\right]{v}_{x}+ \left[\begin{array}{c}0\\ 0\\ 0\\ 1\\ 0\\ 0\end{array}\right]{\mathrm{\omega }}_{x}+ \\ &\left[\begin{array}{c}0\\ 0\\ 0\\ {\rm{sin}}\phi {\rm{tan}}\theta \\ {\rm{cos}}\phi \\ {\rm{sin}}\phi {\rm{sec}}\theta \end{array}\right]{\mathrm{\omega }}_{y}+ \left[\begin{array}{c}0\\ 0\\ 0\\ {\rm{cos}}\phi {\rm{t}}an\theta \\ -{\rm{sin}}\phi \\ {\rm{cos}}\phi {\rm{sec}}\theta \end{array}\right]{\mathrm{\omega }}_{z} ，\end{split} $

(4)

由于UUV运动学模型是在线速度存在非完整约束的条件下得到的^[10]，减少了速度的维数。统有2个非完整约束，分别为绕y轴的线速度和绕z轴的线速度。因为速度维度的减少，所以速度在这2个方向为0。这2个约束可表示为：

$ {s}^{{\rm{T}}}\dot{p}=0 ，$

(5)

$ {a}^{{\rm{T}}}\dot{p} =0，$

(6)

将s，a，p展开并且用A(q) $ \dot{q} $ =0表示为：

$ \begin{split}{{A(q)}} = &\left[ \begin{array}{*{20}{c}} {\rm{cos}}\varPsi {\rm{sin}}\theta {\rm{sin}}\phi - {\rm{sin}}\varPsi {\rm{cos}}\phi\\ {\rm{cos}}\varPsi {\rm{sin}}\theta {\rm{cos}}\phi + {\rm{sin}}\varPsi {\rm{sin}}\phi\end{array}\right.\\ & \left.\begin{array}{*{20}{c}} {\rm{sin}}\varPsi {\rm{sin}}\theta {\rm{sin}}\phi +{\rm{cos}}\varPsi {\rm{cos}}\phi & {\rm{cos}}\theta sin\phi & 0 & 0 & 0\\ {\rm{sin}}\varPsi {\rm{sin}}\theta {\rm{cos}}\phi -{\rm{cos}}\varPsi {\rm{sin}}\phi & {\rm{cos}}\theta {\rm{cos}}\phi & 0 & 0 & 0\end{array} \right]。\end{split} $

(7)

用向量场f₁(g)，f₂(g)，f₃(g)和f₄(g)线性组合称行的矩阵A(q)来描述速度，则建立如下运动模型：

$ \dot{q}\left(t\right)= {f}_{1}\left(q\right){v}_{1} + {f}_{2}\left(q\right){v}_{2}+ {f}_{3}\left(q\right){v}_{3}+ {f}_{4}\left(q\right){v}_{4} 。$

(8)

式中： $ {v}_{1}={v}_{x} $ ； $ {v}_{2}={\mathrm{\omega }}_{x} $ ； $ {v}_{3}={\mathrm{\omega }}_{y} $ ； $ {v}_{4}={\mathrm{\omega }}_{z} $ 。

$ \left\{\begin{array}{*{20}{l}} {f}_{1}\left(q\right)=[{\rm{cos}}\varPsi {\rm{cos}}\theta \;\; {\rm{cos}}\theta {\rm{sin}}\varPsi \;\; -{\rm{sin}}\theta \;\; 0 \;\; 0 \;\; 0]，\\ {f}_{2}\left(q\right)=[0 \;\; 0 \;\; 0 \;\; 1 \;\; 0 \;\; 0]，\\ {f}_{3}\left(q\right)=[0 \;\; 0 \;\; 0 \;\; {\rm{sin}}\phi \mathrm{tan}\theta \;\; {\rm{cos}}\phi \;\; {\rm{sin}}\phi sec\theta ]，\\ {f}_{4}\left(q\right)=[0 \;\; 0 \;\; 0 \;\; {\rm{cos}}\phi \mathrm{tan}\theta \;\; -{\rm{sin}}\phi \;\; {\rm{cos}}\phi \mathrm{sec}\theta ]。\end{array}\right. $

(9)

1.4 模型简化 1.4.1 纵向运动模型

对于UUV，纵向运动可独立存在，将UUV纵向运动单独拆分出来，为空间的垂直平面运动。对于纵向运动，依然遵循牛顿定律和基础动力学方程^[11]，则有：

$ {m}\dot{v}={-A}_{x}{v}^{2}+{T}\mathrm{cos}\alpha -{P}\mathrm{sin}\varTheta，$

(10)

$\begin{split} {m}{v}\dot{\varTheta }=&{A}_{y}^{2}{v}^{2}\mathrm{\alpha }+{A}_{y}^{\delta }{v}^{2}{\mathrm{\delta }}_{\mathrm{e}}{+A}_{y}^{\omega }v{\omega }_{{z}_{1}}+\\ & {\lambda }_{22}\left(\dot{v}\alpha +v\dot{\alpha }\right)-{\lambda }_{26}{\dot{\omega }}_{{z}_{1}}+{T}\mathrm{sin}\mathrm{\alpha }-{P}\mathrm{cos}\varTheta，\end{split}$

(11)

$\begin{split} {J}_{{z}_{1}}{\dot{\omega }}_{{z}_{1}}=&{A}_{{m}_{z}}^{a}{v}^{2}a-{A}_{y}^{\delta }{x}_{e}{v}^{2}{\delta }_{e}-{A}_{{m}_{z1}}^{\omega }v{\omega }_{{z}_{1}}+\\ &{\lambda }_{26}\left(\dot{v}a+v\dot{a}\right)-{\lambda }_{66}{\dot{\omega }}_{{z}_{1}}- \\&B\left(h\mathrm{sin}\theta -{x}_{b}\mathrm{cos}\theta \right)-Th。\end{split}$

(12)

式中： $ \dot{\theta }={\omega }_{{z}_{1}} $ ； $ {\dot{y}}_{e} $ = ${v}\mathrm{sin}\varTheta$ ； ${\dot{x}}_{e}=v\mathrm{cos}\varTheta$ ； $\varTheta =\theta -a$ ； $ {\delta }_{e} $ 为横舵偏角。

1.4.2 侧向运动模型

由于纵向参数在任何时候都不会完全为0，所以侧向运动并不会单独存在。假设UUV是理想的直鳍舵上下对称，那么此时的侧向运动不含横滚，此时的侧向运动方程为：

$ m{v}_{0}\dot{\kappa } = \left({A}_{z}^{\beta }{v}_{0}^{2} + T\right)\beta + {A}_{z}^{\delta }{v}_{0}^{2}{\delta }_{r} + {A}_{z}^{\omega }{v}_{0}{\omega }_{{y}_{1}} + {\lambda }_{33}{v}_{0}\dot{\beta } + {\lambda }_{35}{\dot{\omega }}_{{y}_{1}}，$

(13)

$\begin{split} {J}_{{y}_{1}}{\dot{\omega }}_{{y}_{1}}=&{A}_{{m}_{y}}^{\beta }{v}_{0}^{2}\beta -{A}_{z}^{\delta }{x}_{r}{v}_{0}^{2}{\delta }_{r}-{A}_{{m}_{y1}}^{\omega }{v}_{0}{\omega }_{{y}_{1}}-\\ &{\lambda }_{55}{\dot{\omega }}_{{y}_{1}}-{\lambda }_{35}{v}_{0}\dot{\beta }+T{z}_{b}。\end{split}$

(14)

式中： $ \dot{\kappa }={\omega }_{{y}_{1}} $ ； ${\varLambda }=\kappa -\beta$ ； $ {\delta }_{r} $ 为直舵偏向角； $ {v}_{0} $ 为潜航器的额定速度，可由方程解出。

1.4.3 横滚运动模型

UUV绕x轴转动被称之为横滚，用滚角 $ \mathrm{\phi } $ 来表示其程度大小。分析可得包括横滚的侧向运动模型^[12]关于速度 ${v}$ 的运动模型为：

$\begin{split} & {J}_{{y}_{1}}{\dot{\omega }}_{{y}_{1}}+\left({J}_{{z}_{1}}{-J}_{{y}_{1}}\right){\omega }_{{y}_{1}}{\omega }_{{z}_{1}}={A}_{{m}_{x}}^{\beta }{v}^{2}\beta -{A}_{{m}_{x}}^{\delta }{v}^{2}{\delta }_{d}-\\ & {A}_{{m}_{x1}}^{\beta }v{\omega }_{{x}_{1}}+{A}_{{m}_{x}p}{v}^{2}-{\lambda }_{44}{\dot{\omega }}_{{x}_{1}}-\\& B\left({Z}_{b}\mathrm{cos}\phi +h\mathrm{sin}\phi \right)\mathrm{cos}\theta，\end{split} $

(15)

$ \dot{\phi }={\omega }_{{x}_{1}}-\left({\omega }_{{y}_{1}}{{\rm{cos}}}\phi -{\omega }_{{z}_{1}}{{\rm{si}}n}\phi \right){{\rm{tan}}}\theta ，$

(16)

$\begin{split} {{\rm{cos}}}\varTheta {{\rm{sin}}}\mathrm{\varPhi }=& {{\rm{cos}}}\beta {{\rm{cos}}}\theta {{\rm{sin}}}\phi -{{\rm{sin}}}\alpha {{\rm{sin}}}\beta {{\rm{cos}}}\theta {{\rm{cos}}}\phi +\\ & {{\rm{cos}}}\alpha {{\rm{sin}}}\beta {{\rm{sin}}}\theta。\end{split}$

(17)

式中， $ {\delta }_{d} $ 为差动舵偏角。

2 基于数学模型的航向滑模控制

航向控制可以不考虑横滚运动，但侧向运动不能单独存在^[13]，所以将UUV看做一个以重心为质点的刚体，考虑其纵向及包含横滚的侧向运动。

将简化模型中的纵向运动模型和包含横滚的侧向运动模型合并，得

$\begin{aligned} & u={J}_{{z}_{1}}{\dot{\omega }}_{{z}_{1}}+{J}_{{y}_{1}}{\dot{\omega }}_{{y}_{1}}+\left({J}_{{z}_{1}}{-J}_{{y}_{1}}\right){\omega }_{{y}_{1}}{\omega }_{{z}_{1}} ={A}_{{m}_{z}}^{a}{v}^{2}a-\\ & {A}_{y}^{\delta }{x}_{e}{v}^{2}{\delta }_{e}-{A}_{{m}_{z1}}^{\omega }v{\omega }_{{z}_{1}}+{\lambda }_{26}\left(\dot{v}a+v\dot{a}\right)-{\lambda }_{66}{\dot{\omega }}_{{z}_{1}} -\\ & B\left(h{{\rm{sin}}}\theta -{x}_{b}\mathit{{\rm{cos}}}\theta \right)-Th+{A}_{{m}_{x}}^{\beta }{v}^{2}\beta -{A}_{{m}_{x}}^{\delta }{v}^{2}{\delta }_{d}-{A}_{{m}_{x1}}^{\beta }v{\omega }_{{x}_{1}}+\\ & {A}_{{m}_{x}p}{v}^{2}-{\lambda }_{44}{\dot{\omega }}_{{x}_{1}}-B\left({Z}_{b}\mathrm{cos}\phi +h\mathrm{sin}\phi \right)\mathrm{cos}\theta ，\\[-10pt]\end{aligned}$

(18)

航向控制问题^[14]可以表现为：

$ {\ddot{x}}_{1}=-b{\dot{x}}_{1}+au-F。$

(19)

式中： $ {x}_{1} $ 为航向角；a，b为时间常数；u为简单系统模型；F为不确定性总干扰。

2.1 干扰观测器设计

建立一个标称模型^[15]，把建模误差，噪声，还有实际的外部干扰都给增广为一个新的状态变量，将新的状态变量等效输入到控制端，进行干扰补偿，实现对干扰的完全控制。

$ \left\{\begin{array}{l}\dot{\widehat{F}}={k}_{1}\left(\widehat{\omega }-{\dot{x}}_{1}\right)，\\ \dot{\widehat{\omega }}=-\widehat{F}+au-{k}_{2}\left(\widehat{\omega }-{\dot{x}}_{1}\right)-b{\dot{x}}_{1}。\end{array}\right. $

(20)

式中： $ \dot{\widehat{F}} $ 为对F的估计， $ \widehat{\omega } $ 为对 $ {\dot{x}}_{1} $ 的估计， $ {k}_{1} > $ 0， $ {k}_{2} > $ 0。

$ F=-{\ddot{x}}_{1}-b{\dot{x}}_{1}+au ，$

(21)

$ \widehat{F}=-\dot{\widehat{\omega }}-b{\dot{x}}_{1}+au-{k}_{2}\left(\widehat{\omega }-{\dot{x}}_{1}\right)，$

(22)

则 $\stackrel{~}{F}=\left(\dot{\widehat{\omega }}+{\ddot{x}}_{1}\right)+{k}_{2}\left(\widehat{\omega }-{\dot{x}}_{1}\right)=-\dot{\stackrel{~}{\omega }}-{k}_{2}\stackrel{~}{\omega } 。$

定义Lyapunov函数：

$ {V}_{1}=\frac{1}{2{k}_{1}}{\stackrel{~}{F}}^{2}+\frac{1}{2}{\stackrel{~}{\omega }}^{2}，$

(23)

$ \stackrel{~}{F}=F-\widehat{F}\text{，}\stackrel{~}{\omega }={\dot{x}}_{1}-\widehat{\omega }。$

(24)

合适的参数选取能够更加精准实现外界干扰补偿^[16]，增强系统鲁棒性。取a=5，b=0.15，d=150 sin（0.1 t）， ${k}_{1}=500，{k}_{2}=200 。$

图2为设计的干扰观测器干扰及观测仿真结果，虚线为实际干扰值，点线为观测器估计值。仿真结果显示虚线与点线很好的贴合，说明设计的干扰观测器可以很好实时估计干扰值。

2.2 基于干扰观测器的连续滑模控制

将期望航向角设为 $ {x}_{d} $ ，误差取 $ {Z}_{1} ={x}_{d}-{x}_{1}$ ，则滑模函数为：

$ S={\dot{Z}}_{1}+c{Z}_{1}\text{，} c > 0，$

(25)

$ \dot{{S}}={\ddot{Z}}_{1}+c{\dot{Z}}_{1} ={\ddot{x}}_{d}-{\ddot{x}}_{1}+c{\dot{Z}}_{1} = {\ddot{x}}_{d}+b{\dot{x}}_{1}-au+F+c{\dot{Z}}_{1}，$

(26)

基于干扰补偿的滑模控制器设计为：

$ {\rm{U}}=\frac{1}{a}\left[{\ddot{x}}_{d}+b{\dot{x}}_{1}+c{\dot{Z}}_{1}+\widehat{F}+\eta {\rm{sgn}}\left(s\right)\right] ，$

(27)

其中 $ \stackrel{~}{F}=F-\widehat{F} $ ， $\left|\stackrel{~}{F}\right|\leqslant \eta。$

选择Lyapunov函数为：

$ {V}_{2}=\frac{1}{2}{S}^{2}，$

(28)

$ \begin{aligned} & {\dot{V}}_{2}=S\dot{S}=\\ & S\left({\ddot{x}}_{d}+b{\dot{x}}_{1}-au+F+c{\dot{Z}}_{1}\right)=\\ &{\rm{S}} \left({\ddot{x}}_{d}+b{\dot{x}}_{1}-\left({\ddot{x}}_{d}+b{\dot{x}}_{1}+c{\dot{Z}}_{1}+\widehat{F}+\eta {\rm{sgn}}\left(s\right)\right)+F+c{\dot{Z}}_{1}\right) =\\ &{\rm{S}} \left(F-\widehat{F}-\eta {\rm{sgn}}\left(s\right)\right)= \widehat{F}S-\eta \left|S\right|\leqslant 0 ，\end{aligned}$

(29)

$ {{V}}={{V}}_{1}+{{V}}_{2}=\frac{1}{2{K}_{1}}{\stackrel{~}{F}}^{2} + \frac{1}{2}{\stackrel{~}{\omega }}^{2}+\frac{1}{2}{S}^{2}。$

(30)

仿真中，用 ${\rm{ sat}}\left(t\right)$ 饱和函数代替符号函数 ${\rm{sgn}}\left(s\right)$ 来降低抖动影响^[17]。取a=5，b=0.1，d=150 ${\rm{ sin}}\left(0.1t\right)，{x}_{d}=\mathrm{sin}t，{k}_{1}=500，{k}_{2}=200，\eta =5.0，$ c=15，∆ = 0.10。

图3为干扰值及其误差对比图。可以看出，实际干扰值越大，那么观测器估计值就越精确，实际干扰值越小，观测器估计值的偏差越大，但误差最大不超过5%，符合设计要求。

当预设航向为一个固定值即为一条直线 $ {x}_{d} $ =1时，图4为航向跟踪仿真结果。实线为预设航向，虚线为航向跟踪。由仿真结果可以看出，随着时间增加，航向跟踪曲线快速和预设航向曲线贴合，说明当预设航向为固定值时，采用基于干扰补偿的航向滑模控制器可以实现对设定简单航向的稳定跟踪，并且对干扰具有强撸棒性和自适应性。

当预设航向角设为一个正弦波 $ {x}_{d}=3\mathrm{s}\mathrm{i}\mathrm{n}(3.15t+1.5) $ 时，更方便看出航向角追踪是否出现偏差。图5为预设航向及航向追踪仿真结果，实线为预设航向，虚线为航向跟踪。由仿真结果可以看出，随着时间增加，预设航向曲线和航向跟踪曲线逐渐贴合，说明采用基于干扰补偿的航向滑模控制器可以实现对设定较复杂航向的稳定跟踪，并且对干扰具有强撸棒性和自适应性。

当预设航向为一个自定函数x_d= acot(cot(θ)/sin(Ψ)+tan(Ψ)/sin(θ))，Ψ=atan(xact(2,: ))， ${\theta}=-a{\rm{tan}}\left(x{\rm{act}}\right(4,:)./ {\rm{sec}}(\varPsi\left)\right)$ 。由图6的仿真结果可知，当预设航向角出现航向变化速度过快，变化角度过大时，会使得反馈控制信号具有不连续性，具体体现为高频率切换函数。这对整个控制系统带来了不利影响，导致跟踪结果出现较大偏差。但随时间最终预设航向与实际航向会重合。由以上3种不同函数类型的预设航向仿真结果表明，该连续滑模控制算法在跟踪方面效果较好。

3 结　语

本文对水下无人航行器进行模糊运动学建模，并对包含不确定因素的系统模型设计一种变结构控制和干扰补偿的连续滑模控制器进行航向跟踪仿真研究。当预设航线为一条直线时，系统可以快速调整航向，并与预设航向重合。接着以正弦波预设航向，改变固定频率时，系统能够快速调整航向与预定航向重合，虽然其间会产生一定误差但在系统可控范围。最后设置自定义无规律预设航向，若预设航向变化速度过快或变化角度过大时，系统由于反馈控制信号的不连续性导致无法及时调整航向，但随着时间增加最终会与预设航向重合。以上3个仿真实验证明基于干扰补偿的滑模控制器可以实现对于航向的稳定跟踪，并且具有较强的鲁棒性。