-
在海上航行的船舶为尽快到达目的港,降低航行途中的燃料消耗,最为理想的航行方式就是沿着计划航线航行。但因风、浪、流等的影响,一般船舶很难沿计划好的航线航行。因此,设计具有强鲁棒性的控制器以克服外部扰动,对于保证船舶沿设定的航线航行十分必要。
针对这一问题,众多学者提出了很多先进的控制算法。其中,反步递推算法是解决非线性控制问题的一种典型方法,也是近几年船舶运动非线性控制领域的研究热点[1],在该研究中,出现了许多采用混合方法增加其适应性或鲁棒性的研究成果。例如,张显库等[2]将反步递推算法与鲁棒控制算法相结合,使得算法的鲁棒性得以提高,在将上述算法应用到船舶航向保持控制中后,取得了较好的仿真结果。在使用反步递推算法设计非线性控制器时,一般分为两步或是多步。如果被控对象复杂,最后的非线性控制器未定参数多,则设计过程复杂,用常规的反步递推方法设计出的非线性控制率一般都对消了系统的非线性项[3-4]。董文瀚等[5]提出了一种直接模型,其通过参考反步递推自适应控制系统结构,克服了传统自适应律引入规范化信号后使系统过渡过程品质下降的缺点[6]。张显库[7]在使用反步递推方法设计船舶航向保持控制器的问题上做出了贡献。随后,张显库[8]通过构造简单的Lyapunov函数,将非线性控制器的设计过程由两步简化为了一步,并且设计出的非线性控制律保留了系统的非线性项,但其仿真实验仅用简单的白噪声模拟了风浪的影响。本文运用文献[8]设计的控制器在复杂风浪模型下进行仿真实验,发现其控制效果不是特别理想。为了改进控制效果,减少系统能耗,本文拟在文献[8]所设计控制器的基础上,引入非线性反馈,即使用指数函数处理后的航向误差代替航向误差本身,设计基于Lyapunov稳定性原理的具有非线性反馈的控制器,并将所设计的控制器用于“育鹏”轮船舶模型,实施航向保持控制。
-
以大连海事大学的校船“育鹏”轮为例,参考文献[1, 9],建立“育鹏”轮Nomoto船舶运动响应模型,如图 1所示。图中:δ为舵角;δr为舵角输入;δe为舵角误差;$ \dot \delta $为转舵速率;$ {\dot \delta _{\max }} $为最大转舵速率;δD为干扰;ψ为航向;$ \dot \psi $为转向速率;s为拉普拉斯算子;k为旋回性指数;T为追随性指数。
其非线性二阶Nomoto船舶运动响应模型如式(1)所示。
$$ \ddot \psi + \frac{K}{T}H\left( {\dot \psi } \right) = \frac{K}{T}\delta $$ (1) 式中,$ H\left( {\dot \psi } \right) = \alpha \dot \psi + \beta {\dot \psi ^3}$,其中α和β为通过最小二乘法系统辨识得到的参数[9]。
采用表 1所示“育鹏”轮的船舶压载状态数据,利用Matlab软件的Simulink工具箱,对建立的“育鹏”轮非线性Nomoto船舶运动响应模型进行压载右旋回仿真实验,并与实船压载右旋回测试进行对比,结果如图 2所示。
表 1 “育鹏”轮船舶参数
Table 1. Parameters of ship Yu Peng
参数 压载状态 满载状态 两柱间长L/m 189 189 船宽B/m 27.8 27.8 重心距船舯距离Xc/m -4.04 -1.80 舵叶面积AR/m2 31.67 38.00 航速v/kn 17.3 17.3 排水体积∇/m3 22 036.7 29 268.3 船舯吃水dM/m 6.3 11.0 方形系数Cb 0.661 0.720 旋回性指数K/s-1 0.21 0.08 追随性指数T/s-1 107.78 39.54 最大舵角δmax/(°) 35 35 最大转舵速率$ \dot \delta $max/(°·s-1) ±5 ±5 α 13.14 18.80 β 16 212.5 21 459.9 仿真得到的横向战术直径为3.10L,实船测试所得横向战术直径为3.34L;二者所得的纵向战术直径基本相等。经对比可得,仿真结果的横向符合度为92.8%,整体符合度为96.0%,实验结果验证了模型的有效性[10]。
-
设计船舶非线性控制律,使实际船艏方向ψ跟随期望航向ψR。令$ {x_1} = \psi $,$ {x_2} = {\dot x_1} = r = \dot \psi $,则
$$ {{\dot x}_2} = f\left( {{x_2}} \right) + bu,\;\;\;\;y = {x_1} $$ (2) $$ f\left( {{x_2}} \right) = - KH\left( {\dot \psi } \right)/T,H\left( {\dot \psi } \right) = \alpha \dot \psi + \beta {{\dot \psi }^3} $$ (3) $$ b = K/T $$ (4) 式中:y∈R,为系统输出;u为要设计的控制律,u=δ。
令
$$ {z_1} = {x_1} - {\psi _{\rm{R}}} $$ (5) $$ {z_2} = {x_2} $$ (6) 若最后的控制器镇定状态变量为z1和z2,则原系统在平衡点x1=ψR,x2=r=0处达到一致、渐近稳定。
构造一个Lyapunov函数:
$$ {V_1} = \frac{1}{2}z_2^2 $$ (7) 由于z1和z2之间存在微分关系,若z2被镇定到平衡点0处,则z1同时也被镇定[8]。这是因为若z2被镇定到平衡点0处,即$ {x_2} = {\dot x_1} = 0 $,则x1为定值,同时ψR也为定值,由此可知z1同时被镇定。
z1和z2之间的这种微分关系在实际系统中有一定的普遍性。为简化设计过程,构造V1时没有包含z1的信息,但要保证适当选择控制律,使$ {\dot V_1} $中含有z1。最后,同时镇定z1和z2。
$$ {{\dot V}_1} = {z_2}{{\dot z}_2} $$ (8) $$ {{\dot z}_2} = {{\dot x}_2} = f\left( {{x_2}} \right) + bu $$ (9) 为使$ {\dot V_1} \le 0 $,设计控制器:
$$ u = \frac{1}{b}\left[ {f\left( {{x_2}} \right) - {k_1}\left( {{\omega ^{{z_1}}} - 1} \right)} \right] $$ (10) 式中,k1 > 0,ω > 1,两者均为控制器设计参数。
将ωz1展开成z1的幂级数[11],并忽略高阶小项,则有
$$ {\omega ^{{z_1}}} - 1 \approx {z_1}\ln \omega + \frac{1}{2}z_1^2{\left( {\ln \omega } \right)^2} + \frac{1}{6}z_1^3{\left( {\ln \omega } \right)^3} $$ (11) 将式(11)代入式(10),由式(2)~式(9),可得
$$ \begin{array}{*{20}{c}} {{{\dot V}_1} = {z_2}\left[ {f\left( {{x_2}} \right) + bu} \right] \approx }\\ {{z_2}\left\{ {f\left( {{x_2}} \right) + b \cdot \frac{1}{b}\left[ {f\left( {{x_2}} \right) - {k_1}{z_1}\ln \omega - \frac{{{k_1}}}{2}z_1^2{{\left( {\ln \omega } \right)}^2} - } \right.} \right.}\\ {\left. {\left. {\frac{{{k_1}}}{6}z_1^3{{\left( {\ln \omega } \right)}^3}} \right]} \right\} = 2{x_2}f\left( {{x_2}} \right) - {k_1}{z_1}{z_2}\ln \omega - }\\ {\frac{{{k_1}}}{2}z_1^2{z_2}{{\left( {\ln \omega } \right)}^2} - \frac{{{k_1}}}{6}z_1^3{z_2}{{\left( {\ln \omega } \right)}^3}} \end{array} $$ (12) 根据Young式不等式[12],下列不等式成立:
$$ - z_1^2{z_2} \le \frac{1}{2}z_1^4 + \frac{1}{2}z_2^2 $$ (13) 将式(13)代入式(12),有
$$ \begin{array}{*{20}{c}} {{{\dot V}_1} \le 2{x_2}f\left( {{x_2}} \right) - {k_1}{z_1}{z_2}\ln \omega + }\\ {\frac{{{k_1}}}{2}\left( {\frac{1}{2}z_1^4 + \frac{1}{2}z_2^2} \right){{\left( {\ln \omega } \right)}^2} - \frac{{{k_1}}}{6}z_1^3{z_2}{{\left( {\ln \omega } \right)}^3} = }\\ { - 2b\left( {\alpha x_2^2 + \beta x_2^4} \right) - {k_1}\frac{{{x_1} - {\psi _{\rm{R}}}}}{h}h{x_2}\ln \omega + }\\ {\frac{{{k_1}}}{2}{{\left( {\ln \omega } \right)}^2}\left[ {\frac{1}{2}{{\left( {\frac{{{x_1} - {x_{\rm{R}}}}}{h}h} \right)}^4} + \frac{1}{2}x_2^2} \right] - }\\ {\frac{{{k_1}}}{6}{{\left( {\frac{{{x_1} - {x_{\rm{R}}}}}{h}h} \right)}^3}{x_2}{{\left( {\ln \omega } \right)}^3} \approx - 2b\left( {\alpha x_2^2 + \beta x_2^4} \right) - }\\ {{k_1}hx_2^2\ln \omega + \frac{{{k_1}}}{4}{{\left( {\ln \omega } \right)}^2}\left[ {{h^4}x_2^4 + x_2^2} \right] - \frac{{{k_1}}}{6}{h^3}{{\left( {\ln \omega } \right)}^3}x_2^4 = }\\ { - \left( {2b\alpha + {k_1}h\ln \omega - \frac{{{k_1}}}{4}{{\left( {\ln \omega } \right)}^2}} \right)x_2^2 - }\\ {\left( {2b\beta - \frac{{{k_1}}}{4}{{\left( {\ln \omega } \right)}^2}{h^4} + \frac{{{k_1}}}{6}{h^3}{{\left( {\ln \omega } \right)}^3}} \right)x_2^4} \end{array} $$ (14) 式中,h为系统采样周期。式(14)中,用x2来近似$ {\frac{{{x_1} - {\psi _{\rm{R}}}}}{h}} $并非严格遵循数学上导数的定义,但其在控制工程中是可以接受的。
因为b, α, β均大于0,所以只要下列不等式成立,
$$ {k_1}h\ln \omega - \frac{{{k_1}}}{4}{\left( {\ln \omega } \right)^2} > 0 $$ (15) $$ - \frac{{{k_1}}}{4}{\left( {\ln \omega } \right)^2}{h^4} + \frac{{{k_1}}}{6}{h^3}{\left( {\ln \omega } \right)^3} > 0 $$ (16) 则$ {\dot V_1} < 0 $恒成立。
通过求解不等式(15)~式(16),有
$$ \frac{1}{4}\ln \omega < h < \frac{2}{3}\ln \omega $$ (17) $$ {{\rm{e}}^{1.5h}} < \omega < {{\rm{e}}^{4h}} $$ (18) 仿真实验的采样周期一般满足h≤1 s,本文选取的采样时间h=0.2 s,满足要求。此时,系统设计参数ω在区间(1.35,2.23)内必能保证$ {\dot V_1} < 0 $恒成立。
由Lyapunov稳定性定理可知,本文所设计的控制律(式10)能够使整个系统被镇定,从而系统在平衡点x1=ψR,x2=0处达到一致渐近稳定。式(10)的非线性控制律求解过程只构造了一个Lyapunov函数,且控制律中只有2个待定参数k1,ω,控制器设计过程和参数选定均较简单。
-
以大连海事大学的“育鹏”轮为例进行仿真实验,图 3给出了未加非线性反馈的系统仿真框图。在图 3所示的系统中,操舵伺服系统是由舵机延迟、最大舵角饱和限制、最大舵速限制和积分环节组成。考虑到船舶实际航行时航向受风、浪的影响,仿真时加入了6级风浪干扰,得到了有海浪干扰的仿真结果。航行中风的情况比较复杂,为了准确描述风对船舶航行的影响,本文将风的干扰等效描述为一个白噪声和代表相应风级的等效舵角的合成[13]。根据文献[1, 13],计算得到“育鹏”轮在6级风且风舷角为30°的情况下,风的等效舵角δwind=0.8°。6级风引起的海浪干扰用白噪声驱动的二阶振荡环节[0.419 8s/(s2+0.363 8s+0.367 5)]描述。
应用图 3所示的系统进行仿真实验,输入期望航向60°,为了得到更好的控制效果,调整控制器参数k1。仿真实验中,当控制器参数k1=0.003 5时,控制效果最优,如图 4和图 5中蓝线所示,分别表示反步递推非线性控制器控制下的船舶舵角效果和航向保持效果。由此可知,船舶转向至60°,调节时间约为250 s,操舵最大幅值35°,超调量8°。
舵的能耗体现在操舵平稳性、动舵(0.5°以上)次数、作用时间和舵叶转动幅度等上[14]。本文应用平均舵角${\overline \delta _1} = \frac{1}{{{t_{{n_1}}} - {t_0}}}\int {\left| {\delta \left( t \right)} \right|} {\rm{d}}t $(其中t0=0,tn1为调节时间)来考虑系统能耗[15],得到系统的平均舵角δ1=9.82°。
为优化控制,解决系统能耗较多的问题,受文献[16]的启发,保持控制器不变,在系统中引入了非线性函数反馈,即使用非线性指数函数处理后的航向误差代替航向误差本身,来作为控制器的输入信号。系统设计如图 6所示(图中,ωu-1为指数函数的反馈项,d为干扰)。进行仿真,并调节指数函数的参数ω以得到最优的控制效果。当取参数ω=1.7时,利用Simulink工具箱进行仿真控制的效果最优,如图 4和图 5中红线所示,分别表示具有非线性反馈的控制器控制下的船舶舵角效果和航向保持效果。由此可知,船舶转向至60°,调节时间约为200 s,船舶操舵最大幅值35°,超调量5°,计算得到平均舵角δ2= 8.24°。为量化两者间的能耗差异,引入了节能比P的概念。
$$ P = \frac{{{{\bar \delta }_1} - {{\bar \delta }_2}}}{{{{\bar \delta }_1}}} \times 100\% $$ (19) 由两次仿真结果的对比可知,引入指数函数非线性反馈后,控制器的控制效果更好,其节能比达16.1%。
为了验证非线性反馈控制器的控制效果及节能情况,又分别对20°,40°,60°,80°和100°下的航向保持进行了仿真实验,得到的数据如表 2所示。由表 2中数据可知,引入指数函数非线性反馈控制后,在整个控制过程中,调节时间缩短,平均舵角更小。说明改进后的控制器其控制效果更加优异,系统更加节能。
表 2 2种控制方式对比
Table 2. The comparison of control effects between two methods
稳定航向
ψ/(°)控制类型 调节时间tn2/s 最大舵角
δmax/(°)平均舵角
δ/(°)节能比
P/%20 常规反步递推 200 28 6.66 35.1 非线性反步递推 150 17 4.32 40 常规反步递推 200 35 10.08 19.4 非线性反步递推 150 35 8.12 60 常规反步递推 250 35 9.82 16.1 非线性反步递推 200 35 8.24 80 常规反步递推 250 35 11.47 10.1 非线性反步递推 200 35 10.31 100 常规反步递推 250 35 16.02 21.8 非线性反步递推 200 35 12.53 在船舶实际航行过程中,存在着频繁转向的情况。为研究改进后控制器的控制效果及节能情况,本文进行了船舶航向跟踪仿真实验。仿真环境与图 3所示系统相同,输入航向为(0±20)°,频率0.009 rad/s。仿真结果如图 7和图 8所示,其中蓝线表示反步递推控制器控制下的船舶舵角效果和航向保持效果,红线表示具有非线性反馈控制器控制下的船舶舵角效果和航向保持效果。由图 7可知,引入指数函数非线性反馈后,系统的超调量减小,稳定时间缩短,说明改进后的控制器的控制效果更加优异。由图 8可计算得到,引入指数函数非线性反馈前,整个控制系统的平均舵角δ3=5.59°,引入指数函数非线性反馈后,平均舵角δ4=3.25°,节能比高达41%。由此进一步证明,引入指数函数非线性反馈改进后的控制器的控制效果更优,且系统更加节能。
为了验证改进后的非线性反馈控制器的鲁棒性,经查阅船舶资料,通过计算得到“育鹏”轮满载时的船舶数据及其非线性Nomoto模型的参数,详见表 1。利用表 1所示数据进行船舶转向至60°的仿真实验,结果如图 9和图 10所示。其中蓝线为反步递推控制器的控制效果:整个控制过程中最大舵角为35°,调节时间约130 s,超调量3°;红线为改进的非线性反馈控制器的控制效果:整个控制过程中最大舵角为35°,调节时间约130 s,无超调现象。由此可知,在船舶装载状态发生改变的情况下,本文所设计的改进的控制器的控制效果依然良好,且系统具有鲁棒性。
-
为了设计控制效果更加优异的控制器,本文将指数函数非线性反馈引入到了反步递推控制算法中,从理论上证明了引入指数函数设计的控制器满足Lyapunov稳定性要求,同时,整个控制器的设计过程比较简单。文章以“育鹏”轮为例进行了仿真研究,结果表明,具有非线性反馈的控制器在风浪模型下仍具有较好的控制效果。在整个仿真过程中,控制器能很好地响应操舵指令,操舵平均舵角较小,很好地保护了舵机,节约了能量。同时,“育鹏”轮的满载仿真实验结果表明,具有非线性反馈的控制器具有很好的鲁棒性能。在后续的实验中,作者应用本文设计的控制器对其他种类多艘船舶的航向保持控制进行仿真实验,皆取得了很好的控制效果。
Improved nonlinear control for ship course-keeping based on Lyapunov stability
-
摘要:
目的 船舶航向保持是船舶自动导航的核心,其实现离不开非线性控制算法的设计。为了简化控制器的设计过程,提高系统鲁棒性,节省控制器的能量输出,提出一种基于Lyapunov函数改进的具有非线性反馈的控制器。 方法 首先,通过构建Lyapunov函数,简化控制器设计过程,避免非线性控制律对消系统的非线性项;然后,用指数函数处理后的航向误差代替航向误差本身的非线性反馈项,优化控制系统设计,减少能耗;最后,以实际船舶“育鹏”轮为例进行系统仿真。 结果 结果显示,采用上述方法设计的控制器控制效果良好,更加节能,且系统具有鲁棒性。 结论 研究表明,非线性反馈项的引入可以优化控制系统,其能节能的优点在实际工程应用中具有重要意义。 -
关键词:
- 船舶航向保持 /
- Lyapunov稳定性 /
- 非线性反馈
Abstract:Objectives Course-keeping for ships is the core of automatic navigation in marine transportation, and it cannot be achieved without the design of the nonlinear controller. In order to simplify the design process, improve the robustness of the system and save the energy output of the controller, this paper proposes a nonlinear feedback controller based on the Lyapunov function. Methods By constructing a simplified Lyapunov function, the design process of the nonlinear controller is shortened and the cancellation of the nonlinear term of the control law avoided. The design of the control system can be optimized by replacing the course error with a nonlinear function of itself in the feedback loop. This optimization reduces energy consumption and improves the robustness of the system. Results Taking the ship 'Yu Peng' as an example, the simulation results show that the improved nonlinear control algorithm possesses good performance and robustness, and reduces energy consumption. Conclusions The results of this paper show that the introduction of a nonlinear feedback term can optimize the control system, and the energy saving advantages of such a measure are of great significance for practical engineering application. -
Key words:
- course-keeping for ships /
- Lyapunov stability /
- nonlinear feedback
-
表 1 “育鹏”轮船舶参数
Table 1. Parameters of ship Yu Peng
参数 压载状态 满载状态 两柱间长L/m 189 189 船宽B/m 27.8 27.8 重心距船舯距离Xc/m -4.04 -1.80 舵叶面积AR/m2 31.67 38.00 航速v/kn 17.3 17.3 排水体积∇/m3 22 036.7 29 268.3 船舯吃水dM/m 6.3 11.0 方形系数Cb 0.661 0.720 旋回性指数K/s-1 0.21 0.08 追随性指数T/s-1 107.78 39.54 最大舵角δmax/(°) 35 35 最大转舵速率$ \dot \delta $max/(°·s-1) ±5 ±5 α 13.14 18.80 β 16 212.5 21 459.9 表 2 2种控制方式对比
Table 2. The comparison of control effects between two methods
稳定航向
ψ/(°)控制类型 调节时间tn2/s 最大舵角
δmax/(°)平均舵角
δ/(°)节能比
P/%20 常规反步递推 200 28 6.66 35.1 非线性反步递推 150 17 4.32 40 常规反步递推 200 35 10.08 19.4 非线性反步递推 150 35 8.12 60 常规反步递推 250 35 9.82 16.1 非线性反步递推 200 35 8.24 80 常规反步递推 250 35 11.47 10.1 非线性反步递推 200 35 10.31 100 常规反步递推 250 35 16.02 21.8 非线性反步递推 200 35 12.53 -
[1] 张显库.船舶运动简捷鲁棒控制[M].北京:科学出版社, 2012:1-10. Zhang X K. Simple and robust control of ship motion[M]. Beijing:Science Press, 2012:1-10(in Chinese). [2] 张显库, 郭晨, 杜佳璐.船舶航向不对称信息理论与非线性逆推鲁棒控制算法[J].交通运输工程学报, 2006, 6(2):47-50. doi: 10.3321/j.issn:1671-1637.2006.02.011 Zhang X K, Guo C, Du J L. Asymmetric information theory and nonlinear backstepping robust control algorithm of ship navigation[J]. Journal of Traffic and Transportation Engineering, 2006, 6(2):47-50(in Chinese). doi: 10.3321/j.issn:1671-1637.2006.02.011 [3] Krstic M, Kanellakopoulos I, Kokotovic P V. Nonlinear and adaptive control design[M]. New York:Wiley, 1995:89-108. [4] Do K D, Pan J. Global robust adaptive path following of underactuated ships[J]. Automatica, 2006, 42(10):1713-1722. doi: 10.1016/j.automatica.2006.04.026 [5] 董文瀚, 孙秀霞, 林岩, 等.一类直接模型参考Backstepping自适应控制[J].控制与决策, 2008, 23(9):981-986, 993. doi: 10.3321/j.issn:1001-0920.2008.09.004 Dong W H, Sun X X, Lin Y, et al. Direct model reference Backstepping adaptive control[J]. Control and Decision, 2008, 23(9):981-986, 993(in Chinese). doi: 10.3321/j.issn:1001-0920.2008.09.004 [6] 郭源博, 张晓华, 陈宏钧.基于双线性系统理论的PWM整流器非线性控制[J].西南交通大学学报, 2009, 44(2):232-237. doi: 10.3969/j.issn.0258-2724.2009.02.016 Guo Y B, Zhang X H, Chen H J. Nonlinear control of PWM rectifier based on bilinear system theory[J]. Journal of Southwest Jiaotong University, 2009, 44(2):232-237(in Chinese). doi: 10.3969/j.issn.0258-2724.2009.02.016 [7] 张显库.船舶航向保持的非线性逆推鲁棒控制算法[J].大连海事大学学报, 2007, 33(2):80-83. doi: 10.3969/j.issn.1006-7736.2007.02.019 Zhang X K. Nonlinear backstepping robust control algorithm for ship's course-keeping[J]. Journal of Dalian Maritime University, 2007, 33(2):80-83(in Chinese). doi: 10.3969/j.issn.1006-7736.2007.02.019 [8] 张显库.基于Lyapunov稳定性的船舶航向保持非线性控制[J].西南交通大学学报, 2010, 45(1):140-143. doi: 10.3969/j.issn.0258-2724.2010.01.024 Zhang X K. Nonlinear control for ship course-keeping based on Lyapunov stability[J]. Journal of Southwest Jiao Tong University, 2010, 45(1):140-143(in Chinese). doi: 10.3969/j.issn.0258-2724.2010.01.024 [9] 贾欣乐, 张显库.船舶运动智能控制与H∞鲁棒控制[M].大连:大连海事大学出版社, 2002:273-275. Jia X L, Zhang X K. Ship motion intelligent control and H∞ robust control[M]. Dalian:Dalian Maritime University Press, 2002:273-275(in Chinese). [10] 马超, 张显库.基于指数函数非线性反馈的船舶航向保持控制[J].船舶工程, 2017, 39(增刊1):162-165. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QKC20172017081800007898 Ma C, Zhang X K. Ship course keeping control based on nonlinear feedback of exponential function[J]. Ship Engineering, 2017, 39(Supp 1):162-165(in Chinese). http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QKC20172017081800007898 [11] 同济大学应用数学系.高等数学[M].北京:高等教学出版社, 2005:215-220. Department of Applied Mathematics, Tongji University. Advanced mathematics[M]. Beijing:Higher Education Press, 2005:215-220(in Chinese). [12] Deng H, Krsti Ć M. Stochastic nonlinear stabilization-Ⅰ:a backstepping design[J]. Systems & Control Letters, 1997, 32(3):143-150. http://d.old.wanfangdata.com.cn/Periodical/zdhxb200105003 [13] 贾欣乐, 张显库.风力等效舵角及应用[J].中国造船, 2015, 56(3):116-123. doi: 10.3969/j.issn.1000-4882.2015.03.013 Jia X L, Zhang X K. Equivalent rudder angle to wind force and its application[J]. Shipbuilding of China, 2015, 56(3):116-123(in Chinese). doi: 10.3969/j.issn.1000-4882.2015.03.013 [14] 涂武.船舶航向最优控制与仿真研究[D].哈尔滨: 哈尔滨工程大学, 2008. http://cdmd.cnki.com.cn/Article/CDMD-10217-2009058868.htm Tu W. Research of optimal control of ship heading and using Simulink for simulation[D]. Harbin: Harbin Engineering University, 2008(in Chinese). http://cdmd.cnki.com.cn/Article/CDMD-10217-2009058868.htm [15] 张显库, 杨光平, 张强.一种双极性S函数修饰的非线性船舶航向保持算法[J].大连海事大学学报, 2016, 42(3):15-19. http://d.old.wanfangdata.com.cn/Periodical/dlhsdxxb201603003 Zhang X K, Yang G P, Zhang Q. A kind of bipolar Sigmoid function decorated nonlinear ship course keeping algorithm[J]. Journal of Dalian Maritime University, 2016, 42(3):15-19(in Chinese). http://d.old.wanfangdata.com.cn/Periodical/dlhsdxxb201603003 [16] Zhang X K, Zhang G Q. Design of ship coursekeeping autopilot using a sine function-based nonlinear feedback technique[J]. The Journal of Navigation, 2016, 69(2):246-256. doi: 10.1017/S0373463315000612 -