0 引　言

潜艇运动为六自由度空间运动，表现出大惯性、强非线性、各运动平面强耦合的复杂特点，通常难以获得精确的数学模型，且用于潜艇控制设计的方程系数会由于环境、水深及潜艇状态等的变化而变化，使系统工作点偏离设计状态，影响控制效果，使控制器的设计难度大大增加。

国内外众多学者基于经典控制理论和现代控制理论设计了潜艇控制器。基于伯德图和根轨迹法设计的控制器设计从频域角度入手，结构简单，容易实现，但对参数摄动和扰动鲁棒性较差，尤其是难以处理复杂非线性问题^[1]。针对非线性系统的控制问题，多名学者进行了深入研究^{[2 – 3]}，但控制器结构复杂，难以实现。而通过模糊方法设计非线性鲁棒补偿器，一定程度上可以简化控制器的复杂度^[4]。近年来，滑模控制（SMC，sliding mode control）对满足匹配条件的任意摄动和外界干扰具有较强鲁棒性，且其算法简单、响应快速、易于工程实现，受到广泛关注^{[5 – 10]}。传统的滑模控制方法采用极点配置法、最优控制等方法设计线性切换函数，利用基于边界层方法的准滑模控制降低抖振，设计潜艇垂直面准滑模控制器。潜艇垂直面操纵运动为双输入双输出过程，深度控制通道与纵倾控制通道相互耦合，且耦合程度由潜艇的艇型与升降舵的配置密切相关，因此针对此类双入双出的多变量系统，一般可相对增益分析法分析其各回路的耦合程度^[11]，再进行解耦处理^[12]，在弱耦合和确保鲁棒性的情况下可忽略耦合^{[13 – 14]}。不同机动情况下，潜艇水动力差异巨大，对控制模型的设计和控制器的鲁棒性提出较高要求。文献[15]采用线性化模型研究了潜艇大攻角下的滑模控制，起到了较好的控制效果。同时在不增加模型复杂度的情况下也可结合模糊控制等方法设计滑模控制器以降低抖振^{[16 – 17]}。

基于以上考虑，为了使控制器模型简便且确保控制系统的鲁棒性，本文采用潜艇垂直面线性设计模型，基于双幂次趋近律设计了一个滑模控制器。为了估计系统不可测量的状态，引入了Luenberger状态观测器。将该控制器与经典PID控制器进行仿真对比，并应用于有部分参数摄动及有外界干扰情况下的垂直面控制。

1 潜艇垂直面运动模型

坐标系、名词术语、符号规则均采用ITTC和SNAME术语公报的体系^{[18 – 19]}。固定坐标系（定系）E-ξηζ和运动坐标系（动系）O-xyz如图1所示。

潜艇垂直面非线性操纵运动模型由潜艇标准运动方程简化得到^[20]：

1）纵向方程

$ \begin{split} &m\left( {\dot u - vr + wq} \right) = \\ &\frac{1}{2}\rho {L^4}\left( {X_{qq}'{q^2} + X_{rr}'{r^2} + X_{rp}'rp} \right) + \\ & \frac{1}{2}\rho {L^3}\left( {X_{\dot u}'\dot u + X_{vr}'vr + X_{wq}'wq} \right) + \\ &\frac{1}{2}\rho {L^2}\left( {X_{uu}'{u^2} + X_{vv}'{v^2} + X_{ww}'{w^2}} \right) + \\ &\frac{1}{2}\rho {L^2}{u^2}\left( {X_{{\delta _r}{\delta _r}}'\delta _r^2 + X_{{\delta _b}{\delta _b}}'\delta _b^2 + X_{{\delta _s}{\delta _s}}'\delta _s^2} \right) + \\ &\frac{1}{2}\rho {L^2}\left( {{a_T}{u^2} + {b_T}u{u_c} + {c_T}u_c^2} \right){\text{。}} \end{split} $

(1)

2）垂向方程

$ \begin{split} &m\left( {\dot w - uq + vp} \right) = \\ &\frac{1}{2}\rho {L^4}\left( {Z_{\dot q}'\dot q + Z_{rr}'{r^2} + Z_{rp}'rp} \right) + \\ &\frac{1}{2}\rho {L^3}\left[ {Z_{\dot w}'\dot w + Z_{vr}'vr + Z_{vp}'vp + Z_{\left| q \right|{\delta _s}}'u\left| q \right|{\delta _s}} \right] + \\ &\frac{1}{2}\rho {L^3}\left[ {Z_q'uq + Z_{w\left| q \right|}'\frac{w}{{\left| w \right|}}\left| {{{\left( {{v^2} + {w^2}} \right)}^{{1 / 2}}}} \right|\left| q \right|} \right] + \\ &\frac{1}{2}\rho {L^2}\left[ {Z_0'{u^2} + Z_w'uw + Z_{vv}'{v^2} + Z_{w\left| w \right|}'w\left| {{{\left( {{v^2} + {w^2}} \right)}^{{1 / 2}}}} \right|} \right] + \\ & \frac{1}{2}\rho {L^2}\left( {Z_{{\delta _b}}'{u^2}{\delta _b} + Z_{{\delta _s}}'{u^2}{\delta _s}} \right){\text{。}} \end{split} $

(2)

3）纵倾方程

$ \begin{split} & {I_y}\dot q + \left( {{I_x} - {I_z}} \right)rp = \\ & \frac{1}{2}\rho {L^5}\left( {M_{\dot q}'\dot q + M_{rr}'{r^2} + M_{rp}'rp} \right) + \\ & \frac{1}{2}\rho {L^4}\left[ {M_q'uq + M_{\left| w \right|q}'\left| {{{\left( {{v^2} + {w^2}} \right)}^{{1 / 2}}}} \right|q} \right] + \\ & \frac{1}{2}\rho {L^4}\left[ {M_{\dot w}'\dot w + M_{vr}'vr + M_{vp}'vp + M_{\left| q \right|{\delta _s}}'u\left| q \right|{\delta _s}} \right] + \\ & \frac{1}{2}\rho {L^3}\left[ {M_0'{u^2} + M_w'uw + M_{w\left| w \right|}'w\left| {{{\left( {{v^2} + {w^2}} \right)}^{{1 / 2}}}} \right| + M_{vv}'{v^2}} \right] + \\ & \frac{1}{2}\rho {L^3}\left( {M_{{\delta _b}}'{u^2}{\delta _b} + M_{{\delta _s}}'{u^2}{\delta _s}} \right) - mgh\sin \theta {\text{，}}\quad \qquad \qquad \;\;\;(3) \end{split} $

式中，h为潜艇稳心高。

4）运动方程

$\left\{ \begin{array}{l} \dot \zeta = w\cos \theta - u\sin \theta {\text{，}}\\ \dot \theta = q{\text{。}} \end{array} \right.$

(4)

由于非线性方程不便于控制系统分析，本文作如下假设：1）纵向速度恒定；2）忽略所有非线性项；3）纵倾、横摇很小， $\sin \theta \approx \theta $ ， $\sin \phi \approx \phi $ 。4）忽略舵力与角速度项的耦合，略去 $Y_{\left| r \right|{\delta _r}}'\left| r \right|u{\delta _r}$ 、 $Z_{\left| q \right|{\delta _s}}'u\left| q \right|{\delta _s}$ 、 $M_{\left| q \right|{\delta _s}}'u\left| q \right|{\delta _s}$ 、 ${N_{\left| r \right|{\delta _r}}}u\left| r \right|{\delta _r}$ ，由此得到垂直面线性设计模型：

$\left\{ \begin{array}{l} \dot w = {f_w} + {b_{wb}}{\delta _b} + {b_{ws}}{\delta _s}{\text{，}}\\ \dot q = {f_q} + {b_{qb}}{\delta _b} + {b_{qs}}{\delta _s}{\text{，}}\\ \dot \zeta = - u\theta + w{\text{，}}\\ \dot \theta = q{\text{。}} \end{array} \right.$

(5)

式中： ${f_w} = {{\left( {{a_4}{f_{4w}} - {a_2}{f_{5q}}} \right)} / {\left( {{a_1}{a_4} - {a_2}{a_3}} \right)}}$ ， ${f_q} = {{\left( { - {a_3}{f_{4w}} + {a_1}{f_{5q}}} \right)} / {\left( {{a_1}{a_4} - {a_2}{a_3}} \right)}}$ ， ${f_{4w}} = muq + \displaystyle\frac{1}{2}\rho {L^3}Z_q'uq + \displaystyle\frac{1}{2}\rho {L^2}Z_w'uw$ ， ${f_{5q}} = \displaystyle\frac{1}{2}\rho {L^4}M_q'uq + \displaystyle\frac{1}{2}\rho {L^3}M_w'uw - mgh\theta $ ， ${b_{2b}} = \displaystyle\frac{1}{2}\rho {L^2}Z_{{\delta _b}}'{u^2}$ ， ${b_{2s}} = \displaystyle\frac{1}{2}\rho {L^2}Z_{{\delta _s}}'{u^2}$ ， ${b_{3b}} = \displaystyle\frac{1}{2}\rho {L^3}M_{{\delta _b}}'{u^2}$ ， ${b_{3s}} = \displaystyle\frac{1}{2}\rho {L^3}M_{{\delta _s}}'{u^2}$ ， ${a_1} = m - 0.5\rho {L^3}Z_{\dot w}'$ ， ${a_2} = - 0.5\rho {L^4}Z_{\dot q}'$ ${a_3} = - 0.5\rho {L^4}M_{\dot w}'$ ， ${a_4} \!=\! {I_y} - 0.5\rho {L^5}M_{\dot q}'$ ， ${a_5} \!=\! m - 0.5\rho {L^3}Y_{\dot v}'$ ， ${a_6} = - 0.5\rho {L^4}Y_{\dot p}'$ ${a_7} = - 0.5\rho {L^4}Y_{\dot r}'$ ， ${a_8} = - 0.5\rho {L^4}K_{\dot v}'$ ， ${a_9} = {I_x} - 0.5\rho {L^5}K_{\dot p}'$ ， ${a_{10}} = - 0.5\rho {L^5}K_{\dot r}'$ ${a_{11}} = - 0.5\rho {L^4}N_{\dot v}'$ ， ${a_{12}} = - 0.5\rho {L^5}N_{\dot p}'$ ， ${a_{13}} = {I_z} - 0.5\rho {L^5}N_{\dot r}'$ ， ${k_u} = m - 0.5\rho {L^3}X_{\dot u}'{\text{。}}$

式（9）可进一步写为：

$\left\{ \begin{gathered} \ddot \zeta = {f_\zeta } + {b_{wb}}{\delta _b} + {b_{ws}}{\delta _s} \hfill {\text{，}}\\ \ddot \theta = {f_\theta } + {b_{qb}}{\delta _b} + {b_{qs}}{\delta _s} \hfill{\text{。}} \\ \end{gathered} \right.$

2）舵机模型

首尾舵舵机模型均采用如下模型^[13]：

${T_E}\dot \delta = {K_E}\left( {{\delta _E} - \delta } \right){\text{。}}$

(6)

式中： ${T_E}$ 为舵机时间常数； ${K_E}$ 为比例系数。舵角及舵速限制条件：

首舵 $\left| {{\delta _b}} \right| \leqslant {25^\circ }$ ，尾舵 $\left| {{\delta _s}} \right| \leqslant {30^\circ }$ ；舵速 $\left| {{{\dot \delta }_b}} \right| \leqslant 3{^\circ / {s }}$ ， $\left| {{{\dot \delta }_s}} \right| \leqslant 3{^\circ / {s }}$ 。

2 基于新型趋近律的滑模控制器设计 2.1 潜艇垂直面控制策略

一般潜艇深度自动控制有如下几种模式：1）尾舵自动深度模式，即尾舵单独控制深度和纵倾，指令舵角的解算与深度的偏差和纵倾的偏差均有关，在给定深度变化的同时，纵倾角偏差反馈给尾舵，产生相应舵角以控制纵倾角不超过规定范围；2）首、尾舵联合自动深度模式，即深度和纵倾分别由首舵和尾舵控制。本文综合考虑首舵和尾舵和控制效果，选取第2种模式设计控制器。针对首舵和尾舵控制的耦合，采取策略为：首舵主要控制深度，视尾舵为干扰；尾舵主要控制纵倾，视首舵为干扰。

潜艇垂直面运动分变深机动和定深机动。其中变深步骤为：1）给定纵倾指令 ${\theta _d}$ 和深度指令 ${\zeta _d}$ ，自动控制潜艇下潜或上浮；2）在接近指令深度的某一转换深度 ${\zeta _z}$ 上，将指令纵倾转换为 ${\theta _d} = {0^\circ}$ ，控制潜艇进入指令深度上的无纵倾航行状态。转换深度可根据变深的范围合理取值。分深度和纵倾控制2个子系统，如图2所示。

2.2 传统滑模控制器设计

深度控制系统，定义深度误差 ${e_\zeta } = \zeta - {\zeta _d}$ ，深度误差变化率 ${\dot e_\zeta } = \dot \zeta - {\dot \zeta _d}$ ；定义切换函数： ${\sigma _\zeta } = {\dot e_\zeta } + {k_\zeta }{e_\zeta }$ 。纵倾控制系统，定义纵倾误差 ${e_\theta } = \theta - {\theta _d}$ ，纵倾角误差变化率 ${\dot e_\theta } = \dot \theta - {\dot \theta _d}$ ；定义切换函数： ${\sigma _\theta } = {\dot e_\theta } + {k_\theta }{e_\theta }$ 。 ${k_\zeta }$ ， ${k_\theta }$ ， ${\lambda _\zeta }$ ， ${\lambda _\theta }$ 为正常数。

滑模控制控制由等效控制部分 ${u_{eq}}$ 和切换控制部分 ${u_{vss}}$ 组成。等效控制保证系统状态位于滑模面上；切换控制实现对不确定性和外加干扰的鲁棒控制，其保证到达条件成立并迫使系统状态在滑模上滑动。系统到达滑动面上时，在理想状态下，存在

$\dot \sigma = \ddot e + k\dot e = 0{\text{，}}$

(7)

由式（5）和式（7）可得等效控制：

$\left\{ \begin{array}{l} {\delta _b} = \displaystyle\frac{1}{{{b_{wb}}}}\left( {{{\ddot \zeta }_d} - {k_\zeta }{{\dot e}_\zeta } - {f_\zeta } - {b_{ws}}{\delta _s}} \right){\text{，}}\\ {\delta _s} = \displaystyle\frac{1}{{{b_{qs}}}}\left( {{{\ddot \theta }_d} - {k_\theta }{{\dot e}_\theta } - {f_\theta } - {b_{qb}}{\delta _b}} \right){\text{。}} \end{array} \right.$

(8)

切换控制部分采用准滑模控制^[10]：

${u_{vss}} = K\operatorname{sgn} \left( \sigma \right){\text{，}}$

(9)

式中， $K$ 为切换增益。将符号函数替换为连续光滑的双曲正切函数，可有效降低抖振^[21]：

$\tanh \left( {\mu \sigma } \right) = {{\left( {{e^{\mu \sigma }} - {e^{ - \mu \sigma }}} \right)} / {\left( {{e^{\mu \sigma }} + {e^{ - \mu \sigma }}} \right)}}{\text{，}}$

其中μ>0为待定调节系数，其值大小决定了双曲正切函数拐点的变化快慢。总的滑模控制律为：

$\begin{aligned} & {\delta _b} = \frac{1}{{{b_{wb}}}}\left[ {{{\ddot \zeta }_d} - {k_\zeta }{{\dot e}_\zeta } - {f_\zeta } - {b_{ws}}{\delta _s}} \right] + {K_\zeta }\tanh \left( {{\mu _\zeta }{\sigma _\zeta }} \right){\text{，}}\\ & {\delta _s} = \frac{1}{{{b_{qs}}}}\left[ {{{\ddot \theta }_d} - {k_\theta }{{\dot e}_\theta } - {f_\theta } - {b_{qb}}{\delta _b}} \right] + {K_\theta }\tanh \left( {{\mu _\theta }{\sigma _\theta }} \right){\text{。}} \end{aligned}$

(10)

考虑首舵、尾舵动态响应并忽略相互间耦合影响，首、尾舵控制律 ${\delta _{bc}}$ ， ${\delta _{sc}}$ 分别为：

$\left\{ \begin{array}{l} {\delta _{bc}} = {{{T_{E1}}} / {{K_{E1}}{{\dot \delta }_b}}} + {\delta _b}{\text{，}}\\ {\delta _{sc}} = {{{T_{E2}}} / {{K_{E2}}{{\dot \delta }_s}}} + {\delta _s}{\text{。}} \end{array} \right.$

滑模运动的稳定性完全取决于滑模参数的选取，滑模的存在性问题转化为式（9）系统的原点镇定问题，选取Lyapunov函数为：

$ \begin{align} & {V_\zeta } = \frac{1}{2}\sigma _\zeta ^2,\;\;{\sigma _\zeta } \ne 0;\;\;{V_\theta } = \frac{1}{2}\sigma _\theta ^2,\;\;{\sigma _\theta } \ne 0{\text{；}}\\ & {\sigma _\zeta }{{\dot \sigma }_\zeta } = {\sigma _\zeta }\left( {{f_\zeta } + {b_{wb}}{\delta _b} + {b_{ws}}{\delta _s} - {{\ddot \zeta }_d} + {k_\zeta }{{\dot e}_\zeta }} \right){\rm{ = }}\\ & {K_\zeta }{\sigma _\zeta }\tanh \left( {{\mu _\zeta }{\sigma _\zeta }} \right){\text{；}}\\ &{\sigma _\theta }{{\dot \sigma }_\theta } = {\sigma _\theta }\left( {{f_\theta } + {b_{qb}}{\delta _b} + {b_{qs}}{\delta _s} - {{\ddot \theta }_d} + {k_\theta }{{\dot e}_\theta }} \right){\rm{ = }}\\ & {K_\theta }{\sigma _\theta }\tanh \left( {{\mu _\theta }{\sigma _\theta }} \right){\text{。}} \end{align} $

只要使滑模参数 ${K_\zeta } > 0$ 、 ${K_\theta } > 0$ ，即可满足全局渐近稳定性条件 $\left( {\sigma \to 0} \right)$ ，保证滑动模态的存在性、可达性、稳定性。

2.3 双幂次趋近律设计

为进一步降低抖振，采用双幂次趋近律^[22]：

$\dot \sigma = - \varepsilon {\left| \sigma \right|^\gamma }\tanh \left( {\mu \sigma } \right) - \lambda {\left| \sigma \right|^\chi }\tanh \left( {\mu \sigma } \right){\text{，}}$

(11)

式中， $\gamma > 1$ ， $0 < \chi < 1$ ， $\varepsilon > 0$ ， $\lambda > 0$ 。当系统状态远离滑动模态 $\left( {\left| \sigma \right| > 1} \right)$ 时，第1项起主导作用；当系统状态接近滑动模态 $\left( {\left| \sigma \right| < 1} \right)$ 时，第2项起主导作用，结合后保证了滑动模态趋近品质。

根据滑模可达性，由式（11）可得：

$\begin{array}{l} \sigma \dot \sigma = \sigma \left[ { - \varepsilon {{\left| \sigma \right|}^\gamma }\tanh \left( {\mu \sigma } \right) - \lambda {{\left| \sigma \right|}^\chi }\tanh \left( {\mu \sigma } \right)} \right]{\rm{ = }}\\ - \varepsilon {\left| \sigma \right|^{1 + \gamma }} - \lambda {\left| \sigma \right|^{1 + \chi }} < 0{\text{。}} \end{array}$

滑动态可在有限时间内趋于平衡零点。适当增大参数 $\varepsilon $ 和 $\gamma $ 可加快远离滑动态的趋近速度；适当增大参数 $\lambda $ 和 $\chi $ 可以提高接近滑动态的趋近速度。

当系统存在不确定性有界干扰时，状态 $\sigma $ 和 $\dot \sigma $ 不能保证在有限时间内收敛到 $\left( {0,0} \right)$ ，而仅能收敛到 $\left( {0,0} \right)$ 的一个邻域内。证明参考文献[22]。

结合式（5）、式（7）和式（9）可得控制律：

$ \!\!\!\!\!\begin{aligned} & {\delta _b} \!=\! \displaystyle\frac{1}{{{b_{wb}}}}\left[\!\!\!\!\! \begin{array}{l} {{\ddot \zeta }_d} - {k_\zeta }{{\dot e}_\zeta } \!-\! {f_\zeta } \!-\! {b_{ws}}{\delta _s} \!-\! \\ {\varepsilon _\zeta }{\left| {{\sigma _\zeta }} \right|^{{\gamma _{_\zeta }}}}\!\tanh \left( {{\mu _\zeta }{\sigma _\zeta }} \right) \!-\! {\lambda _\zeta }{\left| {{\sigma _\zeta }} \right|^{{\chi _{_\zeta }}}}\!\tanh \left( {{\mu _\zeta }{\sigma _\zeta }} \right) \end{array} \!\!\!\!\!\right]{\text{，}}\\ & {\delta _s} \!=\! \displaystyle\frac{1}{{{b_{qs}}}}\left[\!\!\!\!\! \begin{array}{l} {{\ddot \theta }_d} - {k_\theta }{{\dot e}_\theta } \!-\! {f_\theta } - {b_{qb}}{\delta _b} \!-\! \\ {\varepsilon _\theta }{\left| {{\sigma _\theta }} \right|^{{\gamma _{_\theta }}}}\!\tanh \left( {{\mu _\theta }{\sigma _\theta }} \right) \!-\! {\lambda _\theta }{\left| {{\sigma _\theta }} \right|^{{\chi _{_\theta }}}}\tanh \left( {{\mu _\theta }{\sigma _\theta }} \right) \end{array}\!\!\!\!\!\right]{\text{。}} \end{aligned} $

(12)

选取Lyapunov函数为：

$ \begin{aligned} & {V_\zeta } = \frac{1}{2}\sigma _\zeta ^2,\;\;{\sigma _\zeta } \ne 0;\;\;{V_\theta } = \frac{1}{2}\sigma _\theta ^2,\;\;{\sigma _\theta } \ne 0{\text{；}}\\ & {\sigma _\zeta }{{\dot \sigma }_\zeta } = {\sigma _\zeta }\left( {{f_\zeta } + {b_{wb}}{\delta _b} + {b_{ws}}{\delta _s} - {{\ddot \zeta }_d} + {k_\zeta }{{\dot e}_\zeta }} \right){\rm{ = }}\\ & {\sigma _\zeta }\left[ { - {\varepsilon _\zeta }{{\left| {{\sigma _\zeta }} \right|}^{{\gamma _{_\zeta }}}}\tanh \left( {{\mu _\zeta }{\sigma _\zeta }} \right) - {\lambda _\zeta }{{\left| {{\sigma _\zeta }} \right|}^{{\chi _{_\zeta }}}}\tanh \left( {{\mu _\zeta }{\sigma _\zeta }} \right)} \right] < 0{\text{；}}\\ & {\sigma _\theta }{{\dot \sigma }_\theta } = \left( {{f_\theta } + {b_{qb}}{\delta _b} + {b_{qs}}{\delta _s} - {{\ddot \theta }_d} + {k_\theta }{{\dot e}_\theta }} \right){\rm{ = }}\\ & {\sigma _\theta }\left[ { - {\varepsilon _\theta }{{\left| {{\sigma _\theta }} \right|}^{{\gamma _{_\theta }}}}\tanh \left( {{\mu _\theta }{\sigma _\theta }} \right) - {\lambda _\theta }{{\left| {{\sigma _\theta }} \right|}^{{\chi _{_\theta }}}}\tanh \left( {{\mu _\zeta }{\sigma _\zeta }} \right)} \right] < 0{\text{。}} \end{aligned} $

该控制律满足全局渐近稳定性条件。

3 Luenberger状态观测器设计

$w$ ， $q$ 等状态量无法直接测得，需设计观测器进行估计。系统方程和量测方程：

$\left\{ \begin{array}{l} \dot X = AX + BU + DF{\text{，}}\\ Y = CX{\text{。}} \end{array} \right.$

(13)

式中： $X = \left[ w \ q \ \theta \ \zeta \right]^{\rm{T}}$ 为潜艇状态变量； $U = \left[ {{\delta _b}}\ {{\delta _s}} \right]^{\rm{T}}$ ； $F = {\left[ {{Z_F}}\ {{M_F}} \right]^{\rm{T}}}$ 为总的水动干扰力和干扰力矩。 $\left( {A,C} \right)$ 为可观测对，设计一个任置极点的全维观测器：

$ \hat X = \left( {A - LC} \right)\hat X + BU + LY{\text{，}}$

(14)

式中， $L$ 为状态观测器的反馈阵。

引入状态观测器后，滑模控制为：

$ \begin{align} & {\delta _b} = \displaystyle\frac{1}{{{b_{wb}}}}\left[ \begin{array}{l} {{\ddot \zeta }_d} - {k_\zeta }{{\dot {\hat e}}_\zeta } - {{\hat f}_\zeta } - {b_{ws}}{\delta _s} - \\ {\varepsilon _\zeta }{\left| {{{\hat \sigma }_\zeta }} \right|^{{\gamma _{_\zeta }}}}\tanh \left( {{\mu _\zeta }{{\hat \sigma }_\zeta }} \right) - {\lambda _\zeta }{\left| {{{\hat \sigma }_\zeta }} \right|^{{\chi _{_\zeta }}}}\tanh \left( {{\mu _\zeta }{{\hat \sigma }_\zeta }} \right) \end{array} \right]{\text{，}}\\ & {\delta _s} = \displaystyle\frac{1}{{{b_{qs}}}}\left[ \begin{array}{l} {{\ddot \theta }_d} - {k_\theta }{{\dot {\hat e}}_\theta } - {{\hat f}_\theta } - {b_{qb}}{\delta _b} - \\ {\varepsilon _\theta }{\left| {{{\hat \sigma }_\theta }} \right|^{{\gamma _{_\theta }}}}\tanh \left( {{\mu _\theta }{{\hat \sigma }_\theta }} \right) - {\lambda _\theta }{\left| {{{\hat \sigma }_\theta }} \right|^{{\chi _{_\theta }}}}\tanh \left( {{\mu _\theta }{{\hat \sigma }_\theta }} \right) \end{array} \right]{\text{。}} \end{align} $

(15)

选取Lyapunov函数为：

$ \begin{array}{l} {V_\zeta } = \frac{1}{2}\hat \sigma _\zeta ^2,\;\;{{\hat \sigma }_\zeta } \ne 0;\;\;{V_\theta } = \frac{1}{2}\hat \sigma _\theta ^2,\;\;{{\hat \sigma }_\theta } \ne 0 {\text{；}}\\ {{\hat \sigma }_\zeta }{{\dot {\hat \sigma} }_\zeta } = {{\hat \sigma }_\zeta }\left[ - {\varepsilon _\zeta }{{\left| {{{\hat \sigma }_\zeta }} \right|}^{{\gamma _{_\zeta }}}}\tanh \left( {{\mu _\zeta }{{\hat \sigma }_\zeta }} \right)\right.- \\ \quad \left. {\lambda _\zeta }{{\left| {{{\hat \sigma }_\zeta }} \right|}^{{\chi _{_\zeta }}}}\tanh \left( {{\mu _\zeta }{{\hat \sigma }_\zeta }} \right) \right] < 0 {\text{；}}\\ {{\hat \sigma }_\theta }{{\dot {\hat \sigma} }_\theta } = {{\hat \sigma }_\theta }\left[ - {\varepsilon _\theta }{{\left| {{{\hat \sigma }_\theta }} \right|}^{{\gamma _{_\theta }}}}\tanh \left( {{\mu _\theta }{{\hat \sigma }_\theta }} \right) \right. -\\ \quad \left. {\lambda _\theta }{{\left| {{{\hat \sigma }_\theta }} \right|}^{{\chi _{_\theta }}}}\tanh \left( {{\mu _\theta }{{\hat \sigma }_\theta }} \right) \right] < 0{\text{。}} \end{array} $

该控制律满足全局渐近稳定性条件。

4 仿真

为验证控制器的正确性，针对潜艇定深和变深运动控制进行仿真和分析。潜艇仿真方程为完整的非线性运动方程。采用文献[1]的潜艇参数。控制器设计参数为： ${k_\zeta } = 0.1$ ， ${k_\theta } = 0.08$ ， ${\varepsilon _\zeta } = 5$ ， ${\varepsilon _\theta } = 15$ ， ${\lambda _\zeta } = 100$ ， ${\lambda _\theta } = 200$ ， ${\mu _\zeta } = 0.000\;51$ ， ${\mu _\theta } = 0.000\;81$ ， ${\gamma _\zeta } = 3$ ， ${\gamma _\theta } = 4$ ， ${\chi _\zeta } = 0.1$ ， ${\chi _\theta } = 0.1$ 。观测器极点配置： $\left[ { { - 7}\ { - 4.7}\ { - 6}\ { - 3} } \right]$ 。

潜艇初始航速 ${u_c} = 10$ kn，初始深度 ${\zeta _0} = 100$ m，指令深度 ${\zeta _d} = 10$ m，转换深度 ${\zeta _z} = 20$ m，纵倾角指令 ${\theta _d} = {5^\circ}$ ，初始潜艇已均衡。

分3种情况将SMC的控制效果和经典PID控制器的控制效果对比仿真。

1）当潜艇参数无摄动和无外界扰动时，图3和图4分别给出了深度、纵倾响应曲线与舵角曲线。由图可见，SMC比PID控制使深度能更快速无超调趋近指令深度，达到指令深度后沿期望深度稳定航行，纵倾符合限制条件。滑模控制首、尾舵打舵比PID控制更为平滑。

2）系统部分参数摄动但无外界扰动。由于控制性能主要受垂向方程和纵倾方程有关参数影响，由式（9）知， ${f_w}$ ， ${f_q}$ 包含了主要系统参数，未验证滑模控制器的鲁棒性，设系统参数进行大范围摄动，令 $f_w^*$ ， $f_q^*$ 为系统实际参数，选取范围： $0.5{f_w} \leqslant f_w^* \leqslant 2{f_w}$ ， $0.5{f_q} \leqslant $ $ f_q^* \leqslant 2{f_q}$ 。控制器参数不变，当系统参数分别在标称模型的50%和200%范围摄动时，SMC比PID控制使深度能更快速无超调趋近指令深度。滑模控制首、尾舵打舵比PID控制更为平滑，如图5和图6所示。

3）当无摄动但有外界扰动时，为便于分析，任意选取正弦扰动信号： ${d_\zeta } = 50\;000\sin \left( {0.1{\text{π}} t} \right)$ ， ${d_\theta } = 50\,000 $ $\cos \left( {0.2{\text{π}} t} \right)$ 。

控制参数及初始条件不变。图7和图8分别给出了深度、纵倾响应曲线与舵角曲线。由图可见，SMC比PID控制使深度能更快速无超调趋近指令深度。滑模控制打舵次数与幅度明显小于PID控制，且收敛更快速。经测试，当外界扰动力量级小于 ${10^6}$ 时，系统仍能保持较好的控制性能，深度仍可快速无超调到达稳态，纵倾符合限制条件，舵角出现较小抖振。而当扰动量级大于 ${10^7}$ ，两类控制器控制下的系统均会出现大幅振荡，舵角出现大幅高频打舵现象，无法满足控制条件。

表1分别给出了3种情况下，在SMC和PID控制下的潜艇深度和纵倾的超调量、稳态误差，SMC的控制效果均优于传统的PID控制效果。

5 结　语

本文基于潜艇垂直面线性模型，忽略首尾舵之间的耦合，采用首舵控制深度，尾舵控制纵倾的策略，采用新型趋近律设计滑模控制器。仿真结果表明，在参数无摄动和有摄动及有无外界干扰情况下滑模控制器控制效果均优于PID控制器，且无明显抖振，为了确保系统的稳定性，对控制器参数的整定优化应做进一步研究，因此可采用垂直面线性模型设计垂直面滑模控制器，从而既降低了控制器复杂度，又提高了控制器的鲁棒性，对工程应用及研究具有重要意义。