張心光，鄒早建，王巖松

（1.上海工程技術(shù)大學(xué) 汽車工程學(xué)院，上海201620；2.上海交通大學(xué) 船舶海洋與建筑工程學(xué)院，上海 200240；3.上海交通大學(xué) 海洋工程國家重點實驗室，上海 200240）

基于支持向量回歸機和粒子群算法的船舶操縱運動模型辨識

張心光1，鄒早建2，3，王巖松1

（1.上海工程技術(shù)大學(xué) 汽車工程學(xué)院，上海201620；2.上海交通大學(xué) 船舶海洋與建筑工程學(xué)院，上海 200240；3.上海交通大學(xué) 海洋工程國家重點實驗室，上海 200240）

基于仿真的Z形試驗數(shù)據(jù)，應(yīng)用支持向量回歸機對船舶操縱運動響應(yīng)模型進行了機理建模，從核函數(shù)結(jié)構(gòu)中得到了模型中的操縱性指數(shù)，并利用建立的響應(yīng)模型進行了Z形試驗的運動預(yù)報，同時引入粒子群算法對懲罰因子C值進行尋優(yōu)，以減少懲罰因子C值選擇的任意性對船舶操縱運動模型辨識精度產(chǎn)生的不利影響。通過將運動預(yù)報結(jié)果同仿真試驗數(shù)據(jù)進行比較，驗證了文中方法的有效性。

船舶操縱；響應(yīng)模型；參數(shù)辨識；支持向量回歸機；粒子群算法

0 引 言

應(yīng)用船舶操縱運動數(shù)學(xué)模型加計算機模擬的方法進行船舶操縱性預(yù)報的關(guān)鍵是精確確定數(shù)學(xué)模型中的水動力導(dǎo)數(shù)和操縱性指數(shù)。目前采用的船舶操縱運動數(shù)學(xué)模型主要包括兩類：一類是水動力模型[1-2]；另一類是響應(yīng)模型[3]。

目前，基于自航模試驗或?qū)嵈囼灥南到y(tǒng)辨識方法是求取船舶操縱運動水動力導(dǎo)數(shù)和操縱性指數(shù)的有效方法之一。極大似然（Maximal Likelihood，ML）估計法、模型參考自適應(yīng)（Model Reference Adaptive，MRA）法等傳統(tǒng)的系統(tǒng)辨識方法曾被廣泛應(yīng)用于船舶操縱運動建模研究中[4]。但是，這些系統(tǒng)辨識方法有一些固有缺陷，如，對數(shù)學(xué)模型的依賴性過大，對初值的依賴性過大等。基于人工智能技術(shù)的人工神經(jīng)網(wǎng)絡(luò)能克服傳統(tǒng)系統(tǒng)辨識方法的缺陷，已被廣泛應(yīng)用于船舶操縱運動建模研究中，例如，Ebada和Abdel-Maksoud應(yīng)用BP神經(jīng)網(wǎng)絡(luò)進行了回轉(zhuǎn)操縱性能預(yù)報研究[5]；Guedes Soares應(yīng)用遞歸神經(jīng)網(wǎng)絡(luò)進行了船舶操縱運動建模研究[6]。然而，人工神經(jīng)網(wǎng)絡(luò)也有一些缺陷，如，泛化性能不足，存在局部極值問題和維數(shù)災(zāi)難問題等。近二十幾年來，在系統(tǒng)辨識領(lǐng)域出現(xiàn)了一種新的人工智能方法—支持向量機（Support Vector Machines，SVM）[7-8]，該方法解決了在人工神經(jīng)網(wǎng)絡(luò)中無法避免的局部極值問題和維數(shù)災(zāi)難問題并具有很好的泛化能力。

當SVM用于參數(shù)辨識時，稱為支持向量回歸機（Support Vector Regression，SVR）。SVR包括最小二乘支持向量回歸機（Least Square-SVR，LS-SVR）、ε—支持向量回歸機（ε-SVR）、ν—支持向量回歸機（ν-SVR）等。Luo等[9]曾應(yīng)用LS-SVR對操縱運動響應(yīng)模型進行了參數(shù)辨識，并利用建立的響應(yīng)模型對10°/10°和20°/20°Z形試驗進行了運動預(yù)報，通過將預(yù)報結(jié)果和自航模試驗結(jié)果進行比較，對所采用方法的有效性進行了驗證。SVR的結(jié)構(gòu)參數(shù)包括懲罰因子C和不敏感因子ε。文獻[10]討論了不敏感因子ε對船舶操縱運動建模的影響，并給出了不敏感因子ε取值與SVR回歸精度的關(guān)系。但是，懲罰因子C值的選取目前尚缺乏理論依據(jù)，且對懲罰因子C值選取的研究，鮮有文獻報道，這嚴重影響了基于SVR的船舶操縱運動模型辨識精度。

本文將粒子群算法應(yīng)用于懲罰因子C值的尋優(yōu)，同時基于仿真得到的10°/10°、15°/15°Z形試驗數(shù)據(jù)，應(yīng)用ε-SVR辨識了船舶操縱運動響應(yīng)模型參數(shù)，并采用所建立的響應(yīng)模型對10°/10°、15°/15°Z形試驗進行了運動預(yù)報，通過將首向角和轉(zhuǎn)艏角速度的預(yù)報結(jié)果同仿真試驗數(shù)據(jù)進行對比，驗證了本文方法的有效性。

1 方法介紹

本文應(yīng)用ε-SVR辨識的響應(yīng)模型包括一階線性響應(yīng)模型和二階線性響應(yīng)模型：

將方程（1）和（2）分別用前向差分離散化，得到：

式中：Δt為采樣間隔，a1、a2、b1和b2是待辨識的參數(shù)。

令T1T2=h，T1+T2=g，則a1、a2、b1、b2和h、g及操縱性指數(shù)K、T的關(guān)系為：

由方程（3）和（4）可分別得到訓(xùn)練樣本對：

一階線性響應(yīng)模型：

二階線性響應(yīng)模型：

其中l(wèi)為樣本數(shù)。

應(yīng)用ε-SVR對上述響應(yīng)模型進行參數(shù)辨識，其算法如下：

（2）選擇適當?shù)膽土P因子C和不敏感因子ε，并選擇適當?shù)暮?/p>

（3）構(gòu)造并求解下面的最優(yōu)化問題：

其中：αi為拉格朗日乘子；由以上最優(yōu)化問題得到最優(yōu)解

（4）構(gòu)造決策函數(shù)：

在粒子群算法（Particle Swarm Optimization，PSO）中，每一個優(yōu)化問題的解都是搜索空間中的一只鳥，稱之為“粒子”。PSO采用速度—位置搜索模型，且所有粒子都有一個根據(jù)優(yōu)化目標定義的適應(yīng)度值，及一個用于決定它們飛翔方向和距離的速度。

剛開始時，PSO隨機初始化一群粒子，其中第s個粒子在d維解空間的位置和速度分別表示為Xs和通過每一次迭代，粒子動態(tài)跟蹤個體極值和全局極值并根據(jù)公式（12）和公式（13）實現(xiàn)速度和位置的更新。

文中利用粒子群算法對懲罰因子C值進行優(yōu)化，其實現(xiàn)步驟包括：①初始化粒子群，包括群體規(guī)模N，每個粒子的位置Xs和速度Vs；②計算Z形試驗的運動預(yù)報誤差并將其作為每個粒子的適應(yīng)度；③ 根據(jù)粒子的適應(yīng)度更新Ps和g；④ 根據(jù)公式（12）和公式（13）更新粒子群的位置和速度；⑤ 如果達到最大迭代次數(shù)或者滿足最小錯誤標準，終止程序；否則，跳轉(zhuǎn)到步驟②。

2 方法驗證

表1 一階線性響應(yīng)模型操縱性指數(shù)結(jié)果對比Tab.1 Comparison of manoeuvring indices of the first-order linear response model

圖1 一階線性響應(yīng)模型10°/10°Z形試驗預(yù)報結(jié)果和仿真試驗值的比較Fig.1 Prediction results of 10°/10°zig-zag manoeuvre compared with simulated data by using the first-order linear response model

本文首先利用由自航模試驗得到的操縱性指數(shù)建立的響應(yīng)模型[4]進行10°/10°和15°/15°Z形試驗仿真，由仿真試驗得到用于辨識的訓(xùn)練樣本，其次是設(shè)置用于懲罰因子C值尋優(yōu)的粒子群算法參數(shù)，包括：（1）群體規(guī)模N。通常，種群太小則不能提供足夠的采樣點，以致算法性能很差；種群太大盡管可以增加優(yōu)化信息，阻止早熟收斂的發(fā)生，但無疑會增加計算量，造成收斂時間太長，表現(xiàn)為收斂速度緩慢。因此，種群規(guī)模一般設(shè)為100～1 000，文中選擇群體規(guī)模N為100；（2）最大迭代次數(shù)。迭代次數(shù)越多，越能保證解的收斂性，但是影響運算速度，文中選取的最大迭代次數(shù)為1 000次；（3）學(xué)習(xí)因子c1和c2。學(xué)習(xí)因子分別用于控制粒子指向自身或鄰域最佳位置的運動。通常，c1+c2≤4，文中取c1=c2=2；（4）粒子維數(shù)d。粒子維數(shù)d取決于待優(yōu)化的參數(shù)個數(shù)，而文中待優(yōu)化的參數(shù)只有懲罰因子C值，故文中的粒子維數(shù)d為1。

依據(jù)文獻[10]和懲罰因子C值的尋優(yōu)結(jié)果，選取C=104和ε=0，在每次辨識時，從每個Z形試驗的數(shù)據(jù)中選取前100個訓(xùn)練樣本，運用ε-SVR并采用線性核分別對船舶操縱運動一階線性響應(yīng)模型和二階線性響應(yīng)模型進行參數(shù)辨識，得到K、T等操縱性指數(shù)。然后，用龍格—庫塔法數(shù)值求解這些響應(yīng)模型，進行Z形試驗的運動預(yù)報，得到首向角和轉(zhuǎn)艏角速度預(yù)報值。

對于一階線性響應(yīng)模型，利用10°/10°Z形試驗仿真數(shù)據(jù)，辨識結(jié)果同由自航模試驗得到的操縱性指數(shù)“試驗值”的對比見表1。將辨識得到的K、T值代入方程（1），進行Z形試驗的運動預(yù)報，得到在一定舵角規(guī)律下的首向角和轉(zhuǎn)艏角速度的預(yù)報，其結(jié)果和仿真試驗數(shù)據(jù)的比較如圖1所示。對于二階線性響應(yīng)模型，利用15°/15°Z形試驗數(shù)據(jù)，辨識結(jié)果同由自航模試驗得到的操縱性指數(shù)“試驗值”的對比見表2。將辨識得到的K、T值代入方程（2），進行Z形試驗的運動預(yù)報，得到在一定舵角規(guī)律下的首向角和轉(zhuǎn)艏角速度的預(yù)報，其結(jié)果和仿真試驗數(shù)據(jù)的比較如圖2所示。

表2 二階線性響應(yīng)模型操縱性指數(shù)結(jié)果對比Tab.2 Comparison of manoeuvring indices of the second-order linear response model

圖2 二階線性響應(yīng)模型15°/15°Z形試驗預(yù)報結(jié)果和仿真試驗值的比較Fig.2 Prediction results of 15°/15°zig-zag manoeuvre compared with simulated data by using the second-order linear response model

由上述各圖所顯示的首向角和轉(zhuǎn)艏角速度預(yù)報值同仿真試驗值的比較結(jié)果，可以得出以下結(jié)論：（1）ε-SVR是一種進行船舶操縱運動響應(yīng)模型參數(shù)辨識的有效方法；（2）粒子群算法是一種有效的懲罰因子C值優(yōu)化方法，為提高船舶操縱運動響應(yīng)模型參數(shù)辨識精度提供了一個有效途徑。

3 結(jié) 論

本文將粒子群算法用于懲罰因子C值的尋優(yōu)，同時基于Z形試驗仿真數(shù)據(jù)，應(yīng)用ε-SVR辨識了船舶操縱運動響應(yīng)模型參數(shù)，并進行了Z形試驗的運動預(yù)報。由預(yù)報結(jié)果同仿真試驗數(shù)據(jù)的對比結(jié)果，可以得出以下結(jié)論：（1）ε-SVR是一種進行船舶操縱運動響應(yīng)模型參數(shù)辨識的有效方法；（2）粒子群算法可有效進行懲罰因子C值的優(yōu)化，為提高船舶操縱運動響應(yīng)模型參數(shù)辨識精度提供了一個有效方法。

[1]Abkowitz M A.Lectures on ship hydrodynamics-steering and manoeuvrability[R].Hydro-and Aerodynamics Laboratory, Report No.Hy-5,Lyngby,Denmark,1964.

[2]Ogawa A,Kasai H.On the mathematical model of manoeuvring motion of ship[J].International Shipbuilding Progress, 1978,25(292):306-319.

[3]Nomoto K,Taguchi T,Honda K,Hirano S.On the steering qualities of ships[J].International Shipbuilding Progress,1957, 4(35):354-370.

[4]賈欣樂,楊鹽生.船舶運動數(shù)學(xué)模型—機理建模與辨識建模[M].大連:大連海事大學(xué)出版社,1999. Jia Xinle,Yang Yansheng.Mathematical model of ship motion-Mechanism modeling and identification modeling[M]. Dalian:Dalian Maritime University Press,1999.

[5]Ebada A,Abdel-Maksoud M.Applying artificial intelligence(A.I.)to predict the limits of ship turning manoeuvres[J]. Jahrbuch der Schiffbautechnischen Gesellschaft,99 Band,2005:132-139.

[6]Moreira L,Guedes S C.Dynamic model of manoeuvrability using recursive neural networks[J].Ocean Engineering,2003, 30(13):1669-1697.

[7]Zhang Xinguang,Zou Zaojian.Estimation of the hydrodynamic coefficients from captive model test results by using support vector machines[J].Ocean Engineering,2013,73:25-31.

[8]Zhang Huaxiang,Cao Linlin,Gao Shuang.A locality correlation preserving support vector machine[J].Pattern Recognition,2014,47(9):3168-3178.

[9]Luo Weilin,Zou Zaojian,Li Tieshan.Application of support vector machine to ship steering[J].Journal of Shanghai Jiaotong University(Science),2009,14(2):462-466.

[10]張心光,鄒早建.支持向量回歸機不敏感因子ε對運動建模的影響[J].中國造船,2011,52(3):46-51. Zhang Xinguang,Zou Zaojian.The effect of SVR insensitive factor on modeling of motion[J].Shipbuilding of China,2011, 52(3):46-51.

Identification of models of ship manoeuvring motion using Support Vector Regression and Particle Swarm Optimization

ZHANG Xin-guang1,ZOU Zao-jian2,3,WANG Yan-song1
(1.Automotive Engineering College,Shanghai University of Engineering Science,Shanghai 201620,China;2.School of Naval Architecture,Ocean and Civil Engineering,Shanghai Jiao Tong University,Shanghai 200240,China;3.State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University,Shanghai 200240,China)

By analyzing the simulated zig-zag test data,Support Vector Regression is applied to conduct mechanism modeling of response models of ship manoeuvring motion.The parameters of response models of ship manoeuvring motion are obtained from the structure of kernel function.The established response models are used to predict zig-zag manoeuvres.To reduce the negative effects produced by choosing randomly the penalty factor C,Particle Swarm Optimization is used to optimize the penalty factor C.The prediction results are compared with the simulated test data to verify the effectiveness of the proposed method.

ship manoeuvring;response model;parameter identification;Support Vector Regression; Particle Swarm Optimization

U661.33

A

10.3969/j.issn.1007-7294.2016.07.009

1007-7294（2016）11-1427-06

2016-08-19

國家自然科學(xué)基金項目（51609132）；上海高校青年教師培養(yǎng)資助計劃（ZZGCD15044）；校啟動基金（2015-66）

張心光（1982-），男，講師，通訊作者，E-mail：zxg416@126.com；鄒早建（1956-），男，教授，博士生導(dǎo)師。