# Parameter Identifiability of Ship Manoeuvring Modeling Using System Identification

Authors: Weilin Luo
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

To improve the feasibility of system identification in the prediction of ship manoeuvrability, several measures are presented to deal with the parameter identifiability in the parametric modeling of ship manoeuvring motion based on system identification. Drift of nonlinear hydrodynamic coefficients is explained from the point of view of regression analysis. To diminish the multicollinearity in a complicated manoeuvring model, difference method and additional signal method are employed to reconstruct the samples. Moreover, the structure of manoeuvring model is simplified based on correlation analysis. Manoeuvring simulation is performed to demonstrate the validity of the measures proposed.


---

## Body

## 1. Introduction

Prediction of ship manoeuvrability at the ship design stage is an important method for evaluating the manoeuvring performances of ships [1]. In general, two ways are available to obtain the manoeuvrability parameters (e.g., advance, transfer, tactical diameter, and overshoots), which include simulation-free and simulation-based methods. For the simulation-free method, one can obtain the manoeuvrability parameters directly from database, full-scale trials, or free-running model tests. For the simulation-based method, either system based manoeuvring simulation or computational fluid dynamics (CFD) based manoeuvring simulation can be employed [2]. No matter which kind of simulation is used, usually the ship manoeuvring model is indispensable, which also determines the accuracy of the prediction of manoeuvrability. Generally, there are two ways to obtain a ship manoeuvring model. The most used way is determination of hydrodynamic derivatives (coefficients) in the mathematical model of ship manoeuvring. The second way is where artificial intelligence techniques are employed to obtain an implicit model (or so-called a black-box model) that describes the input-output mapping characteristics of ship dynamics. Examples are the uses of artificial neural networks (ANN) [3, 4] and support vector machines (SVM) [5, 6].To determine the hydrodynamic coefficients in the mathematical model of ship manoeuvring, several methods can be used, including database, empirical formula, captive model test, CFD calculation, and system identification (SI) combined with free-running model test or full-scale trials. Among them, the SI based method provides an effective and practical way to identify the hydrodynamic parameters in the ship manoeuvring models. When combined with free-running model tests, this method requires relatively low cost because it is easy to generate more manoeuvres after a first set of free-running model tests [2]. When combined with full-scale trials, the so-called scale effects due to different Reynolds numbers between that of the ship model and full-scale can be avoided. During last decades, the SI based method has made good progresses and many techniques have been developed in the application to ship manoeuvring modeling. Examples are least squares (LS) regression [7, 8], model reference method (MRM) [9, 10], extend Kalman filter (EKF) [11, 12], maximum likelihood (ML) estimation [13, 14], recursive prediction error (RPE) method [15, 16], frequency spectrum analysis (FSA) method [17, 18], particle swarm optimization (PSO) [19, 20], genetic algorithm (GA) [21, 22], and support vector regression (SVR) [23, 24]. Generally, there have been many SI applications to the parameter identification of ship manoeuvring models.Like the other means to come to manoeuvring predictions, SI has its own advantages and disadvantages when applied to ship manoeuvring modeling. The main disadvantage is the problem of parameter identifiability. It can be explained from two aspects. One kind of parameter identifiability is whether the parameters can be identified or not. When solving a state space equation with respect to the ship manoeuvring motion, some parameters cannot be identified for instance the added masses or acceleration derivatives, no matter what input-output samples are provided and no matter which SI technique is adopted. The other kind of parameter identifiability is the effect of parameter drift. The identified coefficients might be not physically correct but be mathematically correct [2], which implies that even if the identified parameters deviate from their true values, the manoeuvring model can predict the manoeuvring motion well, compared with the test results. However, this prediction ability is usually restricted to the prediction of the manoeuvre which provides the data samples for system identification. For example, a manoeuvring model for a ship based on the identification of 10°/10° zigzag manoeuvre can predict the same manoeuvre well (i.e., 10°/10° zigzag manoeuvre) but fails when predicting other manoeuvres, for instance, 20°/20° zigzag manoeuvre or turning circle manoeuvres. The effect of parameter drift happens more frequently and severely for a complicated ship manoeuvring model in which a number of hydrodynamic coefficients are contained, for instance, the Abkowitz model.SI application to the modeling of ship manoeuvring had ever been widely studied in the 1980s with the development of SI techniques for instance the applications of EKF, ML, and RPE. Correspondingly, the issues on parameter identifiability were addressed and some measures had been proposed [11, 25, 26]. In the 1990s, studies on the application of SI to the ship manoeuvring had been paid relatively less attention to due to a new charming method (well known as the CFD method) and the limitations of test technology and SI techniques available, including the system observability, the sensitivity to the initial estimation of parameters and noises, and the parameter identifiability especially the effect of parameter drift. With the development of modern measurement equipment and the introduction of advanced SI techniques, SI based modeling of ship manoeuvring received renewed attention during the past decade. On one hand, some classical SI techniques have been continuously applied and developed, for instance, the LS method [8], the EKF method [12], the RPE method [16], and some hybrid methods [27, 28]. On the other hand, several novel and interesting SI methods have been proposed. Representative examples are the artificial intelligence techniques including PSO [19, 20], GA [21, 22], SVM [23, 24], and ANN [29–31]. However, it should be noted that much work (including the work on classical SI or the work on new SI) concerned little about the problem of parameter identifiability especially for the effect of parameter drift or avoided dealing with the problem by selecting a simple manoeuvring model (e.g., the Nomoto model or the response model) in which the parameter drift is weak because few parameters are involved.In this paper, some issues on the parameter identifiability are summarized. Furthermore, several measures to deal with the identifiability problem are presented especially for the diminishment of the parameter drift in the identification of a complicated ship manoeuvring model, the Abkowitz model. Combined with free-running model test results, SVM based system identification is employed to obtain the hydrodynamic coefficients. Based on the identified manoeuvring model, manoeuvring simulation is conducted and prediction results are compared with the test results, which demonstrates the validity of the proposed methods.


## 2. Determination of Added Masses and Added Moment of Inertia

Three kinds of models are available for the prediction of ship manoeuvring, including the Abkowitz model [32], the MMG model [33], and the response model [34]. Comparatively, the first two models are frequently used because they give a more comprehensive understanding of the ship manoeuvring motion. On the other hand, a number of hydrodynamic coefficients are contained in the models, which makes it challenging to determine these coefficients accurately. From the point of view of system identification, one of the difficulties is that not all coefficients can be identified.Taking a linear manoeuvring model as an example, the state space equation can be written as(1)m-Xu˙000m-Yv˙mxG-Yr˙0mxG-Nv˙Izz-Nr˙u˙v˙r˙=Xu000YvYr-mu00NvNr-mxGu0uvr+0YδNδδ,where u is the surge speed, v is the sway speed, r is the yaw rate, δ is the rudder angle, m is the mass of the ship, Izz is the moment of inertia of the ship, xG longitudinal coordinate of the ship’s center of gravity, u0 is the service speed of the ship, Xu˙ et al. are the hydrodynamic derivatives among which the acceleration derivatives -Xu˙, -Yv˙, -Yr˙, -Nv˙ are named as the added masses, and -Nr˙ is the added moment of inertia. Denote(2)X=uvr,P=m-Xu˙000m-Yv˙mxG-Yr˙0mxG-Nv˙Izz-Nr˙,A=Xu000YvYr-mu00NvNr-mxGu0,B=0YδNδ.Equation (1) can be written as(3)X˙=P-1AX+P-1Bδ,if P is invertible.Using SI, two matrix combinations in (3), that is, P-1A and P-1B, can be determined; however, the three coefficient matrices P, A, and B cannot be obtained separately. One solution is to determine one coefficient matrix first in another way, for example, the slender-body theory, empirical formula, captive model test, or CFD calculation. Usually but not always, the mass matrix P is selected to be predetermined in that the abovementioned methods have been successfully applied to determine the added masses and added moment of inertia for example [35–38]. Among these methods, the slender-body theory provides a simple but effective way. Table 1 lists the comparison results of nondimensional added masses and added moment of inertia, between the calculation based on slender-body theory [39] and captive model tests for a Mariner ship [40]. It is assumed that the ship is symmetrical about yz plane in calculation, as such Yr˙′=Nv˙′=0.Table 1
Added masses and added moment of inertia (×10-5).
Coefficient Slender-body theory Captive model test X u ˙ ′ −39.9 −42 Y v ˙ ′ −728 −748 Y r ˙ ′ 0 −9.354 N v ˙ ′ 0 4.646 N r ˙ ′ −60.7 −43.8As can be seen from the comparison results, the slender-body theory gives a satisfactory estimation of three important hydrodynamic derivativesXu˙′, Yv˙′, and Nr˙′. Combined with free-running model tests, SVM based SI was used to identify the hydrodynamic derivatives in an Abkowitz model, except for the added masses and added moment of inertia that were obtained by slender-body method [24].


## 3. Parameter Drift and Diminishment Measures

Another kind of parameter identifiability, referred to as the parameter drift or simultaneous drift of parameters, is more difficult to deal with, especially for the identification of a complicated manoeuvring model such as an Abkowitz model or a nonlinear MMG model. As aforementioned, when parameter drift happens, the identified parameters might be mathematically correct but physically incorrect. Even worse, the “mathematically correct” parameters cannot be obtained.The reason for the drift of hydrodynamic coefficients has not been made fully clear physically. Hwang applied the slender-body theory to explain the mechanism of the drift of linear hydrodynamic coefficients which is due to the dynamic cancellation [25]. However, for the drift of nonlinear hydrodynamic coefficients, it can only be explained from the point of view of regression analysis, not physically but statistically. The drift of nonlinear hydrodynamic coefficients is due to the so-called multicollinearity [23]. In fact, many nonlinear hydrodynamic coefficients in the mathematical model of ship manoeuvring motion have no clear physical meaning.


### 3.1. Simultaneous Drift of Linear Hydrodynamic Coefficients

Brief explanation of the simultaneous drift of linear hydrodynamic coefficients is given as follows. More details can be found in [11, 25].For a slender ship moving on the horizontal plane, the viscous sway force and yaw moment can be calculated as [39](4)Y=-12ρπd2u0v+12Lr,N=-14ρπd2Lu0v+12Lr,where ρ is the fluid density, d is the draft of the ship, u0 is the service speed of the ship assuming u≈u0, and L the length of ship. The above two expressions can be rewritten as(5)Y=Yvv+Yrr,N=Nvv+Nrr,in terms of hydrodynamic derivatives and motion variables. It can be inferred that(6)YvYr=NvNr=2L.During a ship manoeuvre, the pivot of the ship is near to the bow; therefore, it follows that(7)v+12Lr≈0,which indicates that(8)Yvv+Yrr≈0,(9)Nvv+Nrr≈0.Combined with (6), it can be inferred that two pairs of parameters, that is, Yv and Yr, Nv, and Nr, are allowed to become larger or smaller simultaneously. Figure 1 presents the simulation results of nondimensional linear components of sway force and yaw moment for the 10∘/10∘ zigzag manoeuvre of aMariner vessel [40]. As can be seen, the sway force and yaw moment are small at the helming stage when rudder angle switches. At the steady stage when rudder angle keeps unchanged, the following approximations hold:(10)Yv′v′+Yr′r′≈k1Yδ′δ′,Nv′v′+Nr′r′≈k2Nδ′δ′,v′≈k3δ′,r′≈k4δ′,where k1~4 are constants. Based on the above approximations, two linear expressions can be obtained as(11)k3Yv′+k4Yr′≈k1Yδ′k3Nv′+k4Nr′≈k2Nδ′which means that the elements within the parameter set {Yv′,Yr′,Yδ′} or {Nv′,Nr′,Nδ′} possibly drift simultaneously.Figure 1
Nondimensional linear components of sway force and yaw moment for 10°/10° zigzag manoeuvre of aMariner vessel.


### 3.2. Drift of Nonlinear Hydrodynamic Coefficients

From the point of view of regression analysis, the drift of hydrodynamic coefficients, linear and/or nonlinear, results from the so-called multicollinearity. It can be explained that if the input variables of a mathematical model are strongly linearly dependent on each other, their corresponding coefficients in the regression model might be incorrect, even if the predicted results based on the regression model do agree well with the desired. The case is more possible and serious especially when the system contains a large number of parameters to be identified. Equations (8) and (9) present simple examples of linear dependence. It can be inferred that v and r are linearly dependent on each other. This implies that there exists multicollinearity in the model, which results in the parameter drift of Yv and Yr for (8) and Nv and Nr for (9).It is noted that multicollinearity is a commonplace in the regression analysis [41], not a special effect that happens in the system identification of ship manoeuvring. Moreover, this phenomenon cannot be eliminated but only moderated or diminished, because the input variables are always linearly dependent on each other, more or less.It is important to assess the degree of multicollinearity in a system before diminishment measures are taken. Common used assessment tools includeF- and t-tests, Variance Inflation Factor (VIF) or tolerance (reciprocal of VIF), eigenanalysis, determination coefficient, and correlation coefficient. The last method provides a simple and effective way. Using this method, a correlation coefficient between two variables x1 and x2 can be defined as(12)ρx1x2=cov⁡x1,x2Dx1Dx2,where D(x1) and D(x2) denote the variance and cov⁡(x1,x2) denotes the covariance. Usually, the nearer to one the absolute value of ρx1x2 approximates, the more serious the linear dependence between x1 and x2 is. An example is given with respect to a nonlinear ship manoeuvring model, in which the nondimensional yaw-equation can be described as(13)m′xG′-Nv˙′Δv˙′+Izz′-Nr˙′Δr˙′=Δf3′Δu′,Δv′,Δr′,Δδ′,where Δu′, Δv′, Δr′, and Δδ′ are small perturbations from nominal values of surge speed, sway speed, yaw rate, and rudder angle, respectively. The nonlinear function Δf3′ on the right-hand side of (13) can be expressed as [42] (14)Δf3′=Nv′Δv′+Nr′Δr′+Nδ′Δδ′+Nvv′Δv′Δv′+Nvr′Δv′Δr′+Nrr′Δr′Δr′+Nvr′Δv′Δr′+Nδδδ′Δδ′3+Nvvδ′Δv′2Δδ′+Nvvδ′Δv′Δδ′2+Nrδδ′Δr′Δδ′2+Nrrδ′Δr′2Δδ′+Nrvδ′Δr′Δv′Δδ′+N0′.Combined with the test data of 25∘/5∘ zigzag manoeuvre of a KVLCC2 ship which was conducted in the Hamburg Ship Model Basin (HSVA), the correlation coefficients of variables in (14) can be calculated, as listed in Table 2. As can be seen, the degree of linear dependence is heavy for some variables for instance Δv′ and Δv′Δδ′2, with the correlation coefficient 0.99. A heavy multicollinearity might exist with respect to the system (14) which results in the parameter drift, not only for the linear parameters but also the nonlinear parameters.Table 2
Matrix of correlation coefficients.
N v ′ N r ′ N v v ′ N v r ′ N δ ′ N δ δ δ ′ N v δ δ ′ N v v δ ′ N r r ′ N v r ′ N r δ δ ′ N r r δ ′ N r v δ ′ N v ′ 1.00 0.45 0.98 0.67 0.91 0.92 0.99 0.97 0.39 0.84 0.54 0.54 0.68 N r ′  1.00 0.42 0.84 0.53 0.54 0.47 0.43 0.97 0.66 0.91 0.75 0.81 N v v ′   1.00 0.66 0.83 0.85 0.98 1.00 0.32 0.77 0.50 0.42 0.67 N v r ′    1.00 0.62 0.64 0.68 0.66 0.75 0.87 0.92 0.73 1.00 N δ ′     1.00 0.99 0.91 0.83 0.50 0.83 0.61 0.71 0.62 N δ δ δ ′      1.00 0.92 0.85 0.50 0.83 0.62 0.69 0.64 N v δ δ ′       1.00 0.98 0.40 0.83 0.55 0.53 0.69 N v v δ ′        1.00 0.33 0.77 0.50 0.42 0.67 N r r ′         1.00 0.62 0.86 0.81 0.72 N v r ′          1.00 0.78 0.82 0.89 N r δ δ ′           1.00 0.83 0.91 N r r δ ′            1.00 0.73 N r v δ ′             1.00


### 3.3. Measures to Diminish the Drift of Hydrodynamic Coefficients

From the above discussion, it can be concluded that multicollinearity is the main reason for the drift of hydrodynamic coefficients especially for the nonlinear ones. Therefore, the alleviation of multicollinearity is vital for the diminishment of parameter drift.Several remedies are available to moderate the degree of multicollinearity for instance difference method, pruning-variable method [43], principal component regression (PCR) analysis [44], partial least squares (PLS) [45], and a group method of data handling (GMDH) [46]. Some of them were developed to diminish the drift of hydrodynamic coefficients. For example, Hwang proposed parallel processing, exaggerated over- and underestimation, parameter transformation to moderate the drift of linear hydrodynamic coefficients [25]. For nonlinear hydrodynamic coefficients, ridge regression and modification of input scenario were proposed [47, 48]. However, for given samples from real manoeuvring tests, such a modification of excitation might be not applicable; in other words, this approach might only work for simulation data.Since the multicollinearity is an inherent characteristic of a regression model, besides the regression tool (e.g., LS and EKF), modification of the model is a direct and effective way to alleviate the multicollinearity, which can be achieved by altering the structure of the model and reconstructing the input variables in the model, with the aim to decrease the linear dependence of variables. Based on these ideas, three approaches are presented to diminish the drift of hydrodynamic coefficients, including simplification of the manoeuvring model, difference of the samples, and the method of additional signal. Studies are performed mainly combined with the Abkowitz model.


#### 3.3.1. Simplification of the Manoeuvring Model

Although the Abkowitz model is often preferred in the prediction of ship manoeuvring, such a “large” parametric system makes it inconvenient to use in practice. Actually some nonlinear derivatives are difficult to be determined although they could be unimportant. In practical application of the Abkowitz model, simplification had been performed to improve the feasibility of this model. Examples are the work by Abkowitz [11] and Norrbin Nils [37]. Physical meanings were the main concern in the simplification. The simplification can also be performed from the point of view of regression analysis [24]. Taking the yaw-equation (13) as an example, the original expression of (14) is(15)Δf3′=N0u′Δu′+N0uu′Δu′2+Nv′Δv′+Nr′Δr′+Nδ′Δδ′+Nvvv′Δv′3+Nδδδ′Δδ′3+Nvvr′Δv′2Δr′+Nvvδ′Δv′2Δδ′+Nvδδ′Δv′Δδ′2+Nδu′Δδ′Δu′+Nvu′Δv′Δu′+Nru′Δr′Δu′+Nδuu′Δδ′Δu′2+Nrrr′Δr′3+Nvrr′Δv′Δr′2+Nvuu′Δv′Δu′2+Nruu′Δr′Δu′2+Nrδδ′Δr′Δδ′2+Nrrδ′Δr′2Δδ′+Nrvδ′Δr′Δv′Δδ′+N0′,in which 22 hydrodynamic coefficients are contained [40]. To simplify it, all the nonlinear terms related to the surge perturbation Δu′ are eliminated for two reasons. First, since the ship total velocity U(t) is taken as the nondimensional factor, the effect of velocity loss has been taken into account in every velocity term. Second, the linear dependence between the linear and nonlinear terms relating to Δu′ is significant, which will result in the parameter drift. For example, the third and twelfth terms on the right hand side of (15) at the kth sampling time can be written as(16)Nv′ΔvkUk+Nvu′ΔvkΔukby considering the prime nondimensional definition(17)Δu′=ΔuU,Δu˙′=Δu˙U2/L.For moderate manoeuvres, it holds that(18)Uk≈u0+Δuk.The expression (16) can be written as(19)Nv′Δvku0+Nv′ΔvkΔuk+Nvu′ΔvkΔuk.Obviously, the second term is linearly dependent on the third term. In consideration of the significance of the linear derivative Nv′, the nonlinear term Nvu′Δv(k)Δu(k) is removed. In the same way, the other six terms related to the derivatives N0uu′, Nδu′, Nru′, Nδuu′, Nvuu′, and Nruu′ can be removed from (15). Ignoring the effect of surge speed u on the change of yaw moment, N0u′Δu′ can be removed furthermore. Finally, by introducing the cross flow models to replace the third-order nonlinear terms, model (15) is simplified to model (14), with the number of hydrodynamic coefficients decreasing from 22 to 14.Sensitivity analysis provides another way to simplify the mathematical model of ship manoeuvring [49–53]. Using this method, the effect of the hydrodynamic coefficient on the prediction of ship manoeuvring is studied by observing the change of manoeuvring (e.g., velocity, angular velocity, advance, transfer, overshoots, and tactical diameter) to the change of hydrodynamic coefficients. Based on the results of sensitivity analysis, some trivial hydrodynamic coefficients are removed from the mathematical model. The main difficulties with this method are the determination of two factors: (i) the criteria of importance for a hydrodynamic coefficient which decides this coefficient will be retained or removed from the model; (ii) the reference value of a hydrodynamic coefficient when the real experimental data are used for identification, instead of the simulation data.


#### 3.3.2. Difference Method

Besides simplification of the mathematical model of ship manoeuvring, reconstruction of the input variables or samples provides another effective way to alleviate the parameter drift. It is noted that such a reconstruction does not require a new manoeuvring to obtain new input-output samples for system identification but modify the obtained samples. Difference method provides a simple way to reconstruct the samples [24].Taking (13) as an example, the following discrete equation can be obtained:(20)Δrk+1=Δrk+θ1ΔvkUk+θ2ΔrkUk+θ3ΔδkU2k+θ4ΔvkΔvk+θ5ΔvkΔrk+θ6ΔrkΔrk+θ7ΔvkΔrk+θ8Δδ3kU2k+θ9Δv2kΔδk+θ10ΔvkΔδ2kUk+θ11ΔrkΔδ2kUk+θ12Δr2kΔδk+θ13ΔrkΔvkΔδk+θ14U2k,where θ1~14 are the linear combinations of the hydrodynamic coefficients with respect to the sway and yaw motion [23]. Define vectors(21)Θ=θ1,θ2,…,θ1414×1,Xk=ΔvkUk,ΔrkUk,ΔδkU2k,ΔvkΔvk,ΔvkΔrk,ΔrkΔrk,ΔvkΔrk,Δδ3kU2k,Δv2kΔδk,ΔvkΔδ2kUk,ΔrkΔδ2kUk,Δr2kΔδk,ΔrkΔvkΔδk,U2k14×1T.Equation (20) can be rewritten as(22)Δrk+1=Δrk+ΘXk.At a neighbour sampling time, the above equation becomes(23)Δrk=Δrk-1+ΘXk-1.Combined with (22), the following equality holds: (24)zk=zk-1+ΘYk-1,by definition of(25)zk=Δrk+1-Δrk,Yk-1=Xk-Xk-1.Comparing (23) with (24), it can be seen that although the input-output pair are altered, the coefficient matrix remains the same. Usually, if the degree of linear dependence of input variables in the original system is serious, it could be moderated after difference of samples. As a comparison, Table 3 lists the correlation coefficients after the difference method is employed, as shown below the principal diagonal (elements above the principal diagonal are the correlation coefficients without difference method shown in Table 2). As can be seen, the linear dependence of the input variables has been obviously decreased.Table 3
Comparison of correlation coefficients using difference method.
N v ′ N r ′ N v v ′ N v r ′ N δ ′ N δ δ δ ′ N v δ δ ′ N v v δ ′ N r r ′ N v r ′ N r δ δ ′ N r r δ ′ N r v δ ′ N v ′ 1.00 0.45 0.98 0.67 0.91 0.92 0.99 0.97 0.39 0.84 0.54 0.54 0.68 N r ′ 0.44 1.00 0.42 0.84 0.53 0.54 0.47 0.43 0.97 0.66 0.91 0.75 0.81 N v v ′ 0.83 0.31 1.00 0.66 0.83 0.85 0.98 1.00 0.32 0.77 0.50 0.42 0.67 N v r ′ 0.34 0.21 0.38 1.00 0.62 0.64 0.68 0.66 0.75 0.87 0.92 0.73 1.00 N δ ′ 0.64 0.44 0.23 0.01 1.00 0.99 0.91 0.83 0.50 0.83 0.61 0.71 0.62 N δ δ δ ′ 0.50 0.33 0.25 0.13 0.74 1.00 0.92 0.85 0.50 0.83 0.62 0.69 0.64 N v δ δ ′ 0.84 0.41 0.77 0.36 0.48 0.66 1.00 0.98 0.40 0.83 0.55 0.53 0.69 N v v δ ′ 0.81 0.31 0.99 0.39 0.22 0.29 0.82 1.00 0.33 0.77 0.50 0.42 0.67 N r r ′ 0.34 0.88 0.03 0.07 0.54 0.35 0.28 0.04 1.00 0.62 0.86 0.81 0.72 N v r ′ 0.57 0.27 0.21 0.46 0.58 0.41 0.40 0.19 0.30 1.00 0.78 0.82 0.89 N r δ δ ′ 0.19 0.17 0.08 0.58 0.21 0.29 0.09 0.06 0.05 0.30 1.00 0.83 0.91 N r r δ ′ 0.29 0.32 0.15 0.16 0.72 0.44 0.09 0.18 0.49 0.62 0.50 1.00 0.73 N r v δ ′ 0.32 0.11 0.33 0.90 0.01 0.09 0.26 0.32 0.16 0.49 0.73 0.24 1.00


#### 3.3.3. Additional Signal

Besides the difference method, another approach to reconstruct the input-output samples can be used, which refers to the method of additional signal. Such a method had been successfully applied to the identification of ship manoeuvring by analyzing simulation data [23] and captive model test results [54]. Taking the system (22) as an example, an additional ramp signal can be designed as(26)yak+1=yak+λ,where λ is a constant. Adding this signal to (22), a new parameter system can be obtained as(27)Δrk+1+yak+1=Δrk+yak+ΘXk+λ,where the output variable and the first input variable change; however, the parameter matrix Θ related to hydrodynamic coefficients does not. The additional signal can be viewed as an additional excitation to the original system. As a result, the manoeuvring characteristics of the original system, which might cause the parameter drift as described in the Section 3.1, are changed and the parameter drift can be alleviated resultantly. From the point of view of regression analysis, the introduction of the additional signal alters the input variables, which might alleviate the degree of linear dependence between the changed variables and others. Taking model (27) as an example, assuming λ=1, Table 4 lists the correlation coefficients before and after the additional signal is attached to the samples, in which N1=N1′=1 indicates that coefficient of the first input variable Δr(k)+ya(k) in (27) is one. As can be seen from the comparison results, the degree of linear dependence remarkably decreases for most of variables after an additional signal is added.Table 4
Comparison of correlation coefficients using additional signal.
N 1 N 1 ′ N v ′ 0.465 0.081 N r ′ 0.996 0.011 N v v ′ 0.450 0.026 N v r ′ 0.849 0.002 N δ ′ 0.552 0.028 N δ δ δ ′ 0.556 0.022 N v δ δ ′ 0.491 0.005 N v v δ ′ 0.455 0.029 N r r ′ 0.023 0.834 N v r ′ 0.970 0.039 N r δ δ ′ 0.665 0.015 N r r δ ′ 0.912 0.013 N r v δ ′ 0.771 0.024It is noted that the design of additional signal depends on the characteristic of the data samples or the manoeuvring motion. For periodically nonlinear motion such as zigzag manoeuvre or turning circle manoeuvring, the linear ramp signal like (26) works well [23]; for a linear motion such as oblique towing test, a nonlinear additional signal is more appropriate [54].


## 3.1. Simultaneous Drift of Linear Hydrodynamic Coefficients

Brief explanation of the simultaneous drift of linear hydrodynamic coefficients is given as follows. More details can be found in [11, 25].For a slender ship moving on the horizontal plane, the viscous sway force and yaw moment can be calculated as [39](4)Y=-12ρπd2u0v+12Lr,N=-14ρπd2Lu0v+12Lr,where ρ is the fluid density, d is the draft of the ship, u0 is the service speed of the ship assuming u≈u0, and L the length of ship. The above two expressions can be rewritten as(5)Y=Yvv+Yrr,N=Nvv+Nrr,in terms of hydrodynamic derivatives and motion variables. It can be inferred that(6)YvYr=NvNr=2L.During a ship manoeuvre, the pivot of the ship is near to the bow; therefore, it follows that(7)v+12Lr≈0,which indicates that(8)Yvv+Yrr≈0,(9)Nvv+Nrr≈0.Combined with (6), it can be inferred that two pairs of parameters, that is, Yv and Yr, Nv, and Nr, are allowed to become larger or smaller simultaneously. Figure 1 presents the simulation results of nondimensional linear components of sway force and yaw moment for the 10∘/10∘ zigzag manoeuvre of aMariner vessel [40]. As can be seen, the sway force and yaw moment are small at the helming stage when rudder angle switches. At the steady stage when rudder angle keeps unchanged, the following approximations hold:(10)Yv′v′+Yr′r′≈k1Yδ′δ′,Nv′v′+Nr′r′≈k2Nδ′δ′,v′≈k3δ′,r′≈k4δ′,where k1~4 are constants. Based on the above approximations, two linear expressions can be obtained as(11)k3Yv′+k4Yr′≈k1Yδ′k3Nv′+k4Nr′≈k2Nδ′which means that the elements within the parameter set {Yv′,Yr′,Yδ′} or {Nv′,Nr′,Nδ′} possibly drift simultaneously.Figure 1
Nondimensional linear components of sway force and yaw moment for 10°/10° zigzag manoeuvre of aMariner vessel.


## 3.2. Drift of Nonlinear Hydrodynamic Coefficients

From the point of view of regression analysis, the drift of hydrodynamic coefficients, linear and/or nonlinear, results from the so-called multicollinearity. It can be explained that if the input variables of a mathematical model are strongly linearly dependent on each other, their corresponding coefficients in the regression model might be incorrect, even if the predicted results based on the regression model do agree well with the desired. The case is more possible and serious especially when the system contains a large number of parameters to be identified. Equations (8) and (9) present simple examples of linear dependence. It can be inferred that v and r are linearly dependent on each other. This implies that there exists multicollinearity in the model, which results in the parameter drift of Yv and Yr for (8) and Nv and Nr for (9).It is noted that multicollinearity is a commonplace in the regression analysis [41], not a special effect that happens in the system identification of ship manoeuvring. Moreover, this phenomenon cannot be eliminated but only moderated or diminished, because the input variables are always linearly dependent on each other, more or less.It is important to assess the degree of multicollinearity in a system before diminishment measures are taken. Common used assessment tools includeF- and t-tests, Variance Inflation Factor (VIF) or tolerance (reciprocal of VIF), eigenanalysis, determination coefficient, and correlation coefficient. The last method provides a simple and effective way. Using this method, a correlation coefficient between two variables x1 and x2 can be defined as(12)ρx1x2=cov⁡x1,x2Dx1Dx2,where D(x1) and D(x2) denote the variance and cov⁡(x1,x2) denotes the covariance. Usually, the nearer to one the absolute value of ρx1x2 approximates, the more serious the linear dependence between x1 and x2 is. An example is given with respect to a nonlinear ship manoeuvring model, in which the nondimensional yaw-equation can be described as(13)m′xG′-Nv˙′Δv˙′+Izz′-Nr˙′Δr˙′=Δf3′Δu′,Δv′,Δr′,Δδ′,where Δu′, Δv′, Δr′, and Δδ′ are small perturbations from nominal values of surge speed, sway speed, yaw rate, and rudder angle, respectively. The nonlinear function Δf3′ on the right-hand side of (13) can be expressed as [42] (14)Δf3′=Nv′Δv′+Nr′Δr′+Nδ′Δδ′+Nvv′Δv′Δv′+Nvr′Δv′Δr′+Nrr′Δr′Δr′+Nvr′Δv′Δr′+Nδδδ′Δδ′3+Nvvδ′Δv′2Δδ′+Nvvδ′Δv′Δδ′2+Nrδδ′Δr′Δδ′2+Nrrδ′Δr′2Δδ′+Nrvδ′Δr′Δv′Δδ′+N0′.Combined with the test data of 25∘/5∘ zigzag manoeuvre of a KVLCC2 ship which was conducted in the Hamburg Ship Model Basin (HSVA), the correlation coefficients of variables in (14) can be calculated, as listed in Table 2. As can be seen, the degree of linear dependence is heavy for some variables for instance Δv′ and Δv′Δδ′2, with the correlation coefficient 0.99. A heavy multicollinearity might exist with respect to the system (14) which results in the parameter drift, not only for the linear parameters but also the nonlinear parameters.Table 2
Matrix of correlation coefficients.
N v ′ N r ′ N v v ′ N v r ′ N δ ′ N δ δ δ ′ N v δ δ ′ N v v δ ′ N r r ′ N v r ′ N r δ δ ′ N r r δ ′ N r v δ ′ N v ′ 1.00 0.45 0.98 0.67 0.91 0.92 0.99 0.97 0.39 0.84 0.54 0.54 0.68 N r ′  1.00 0.42 0.84 0.53 0.54 0.47 0.43 0.97 0.66 0.91 0.75 0.81 N v v ′   1.00 0.66 0.83 0.85 0.98 1.00 0.32 0.77 0.50 0.42 0.67 N v r ′    1.00 0.62 0.64 0.68 0.66 0.75 0.87 0.92 0.73 1.00 N δ ′     1.00 0.99 0.91 0.83 0.50 0.83 0.61 0.71 0.62 N δ δ δ ′      1.00 0.92 0.85 0.50 0.83 0.62 0.69 0.64 N v δ δ ′       1.00 0.98 0.40 0.83 0.55 0.53 0.69 N v v δ ′        1.00 0.33 0.77 0.50 0.42 0.67 N r r ′         1.00 0.62 0.86 0.81 0.72 N v r ′          1.00 0.78 0.82 0.89 N r δ δ ′           1.00 0.83 0.91 N r r δ ′            1.00 0.73 N r v δ ′             1.00


## 3.3. Measures to Diminish the Drift of Hydrodynamic Coefficients

From the above discussion, it can be concluded that multicollinearity is the main reason for the drift of hydrodynamic coefficients especially for the nonlinear ones. Therefore, the alleviation of multicollinearity is vital for the diminishment of parameter drift.Several remedies are available to moderate the degree of multicollinearity for instance difference method, pruning-variable method [43], principal component regression (PCR) analysis [44], partial least squares (PLS) [45], and a group method of data handling (GMDH) [46]. Some of them were developed to diminish the drift of hydrodynamic coefficients. For example, Hwang proposed parallel processing, exaggerated over- and underestimation, parameter transformation to moderate the drift of linear hydrodynamic coefficients [25]. For nonlinear hydrodynamic coefficients, ridge regression and modification of input scenario were proposed [47, 48]. However, for given samples from real manoeuvring tests, such a modification of excitation might be not applicable; in other words, this approach might only work for simulation data.Since the multicollinearity is an inherent characteristic of a regression model, besides the regression tool (e.g., LS and EKF), modification of the model is a direct and effective way to alleviate the multicollinearity, which can be achieved by altering the structure of the model and reconstructing the input variables in the model, with the aim to decrease the linear dependence of variables. Based on these ideas, three approaches are presented to diminish the drift of hydrodynamic coefficients, including simplification of the manoeuvring model, difference of the samples, and the method of additional signal. Studies are performed mainly combined with the Abkowitz model.


### 3.3.1. Simplification of the Manoeuvring Model

Although the Abkowitz model is often preferred in the prediction of ship manoeuvring, such a “large” parametric system makes it inconvenient to use in practice. Actually some nonlinear derivatives are difficult to be determined although they could be unimportant. In practical application of the Abkowitz model, simplification had been performed to improve the feasibility of this model. Examples are the work by Abkowitz [11] and Norrbin Nils [37]. Physical meanings were the main concern in the simplification. The simplification can also be performed from the point of view of regression analysis [24]. Taking the yaw-equation (13) as an example, the original expression of (14) is(15)Δf3′=N0u′Δu′+N0uu′Δu′2+Nv′Δv′+Nr′Δr′+Nδ′Δδ′+Nvvv′Δv′3+Nδδδ′Δδ′3+Nvvr′Δv′2Δr′+Nvvδ′Δv′2Δδ′+Nvδδ′Δv′Δδ′2+Nδu′Δδ′Δu′+Nvu′Δv′Δu′+Nru′Δr′Δu′+Nδuu′Δδ′Δu′2+Nrrr′Δr′3+Nvrr′Δv′Δr′2+Nvuu′Δv′Δu′2+Nruu′Δr′Δu′2+Nrδδ′Δr′Δδ′2+Nrrδ′Δr′2Δδ′+Nrvδ′Δr′Δv′Δδ′+N0′,in which 22 hydrodynamic coefficients are contained [40]. To simplify it, all the nonlinear terms related to the surge perturbation Δu′ are eliminated for two reasons. First, since the ship total velocity U(t) is taken as the nondimensional factor, the effect of velocity loss has been taken into account in every velocity term. Second, the linear dependence between the linear and nonlinear terms relating to Δu′ is significant, which will result in the parameter drift. For example, the third and twelfth terms on the right hand side of (15) at the kth sampling time can be written as(16)Nv′ΔvkUk+Nvu′ΔvkΔukby considering the prime nondimensional definition(17)Δu′=ΔuU,Δu˙′=Δu˙U2/L.For moderate manoeuvres, it holds that(18)Uk≈u0+Δuk.The expression (16) can be written as(19)Nv′Δvku0+Nv′ΔvkΔuk+Nvu′ΔvkΔuk.Obviously, the second term is linearly dependent on the third term. In consideration of the significance of the linear derivative Nv′, the nonlinear term Nvu′Δv(k)Δu(k) is removed. In the same way, the other six terms related to the derivatives N0uu′, Nδu′, Nru′, Nδuu′, Nvuu′, and Nruu′ can be removed from (15). Ignoring the effect of surge speed u on the change of yaw moment, N0u′Δu′ can be removed furthermore. Finally, by introducing the cross flow models to replace the third-order nonlinear terms, model (15) is simplified to model (14), with the number of hydrodynamic coefficients decreasing from 22 to 14.Sensitivity analysis provides another way to simplify the mathematical model of ship manoeuvring [49–53]. Using this method, the effect of the hydrodynamic coefficient on the prediction of ship manoeuvring is studied by observing the change of manoeuvring (e.g., velocity, angular velocity, advance, transfer, overshoots, and tactical diameter) to the change of hydrodynamic coefficients. Based on the results of sensitivity analysis, some trivial hydrodynamic coefficients are removed from the mathematical model. The main difficulties with this method are the determination of two factors: (i) the criteria of importance for a hydrodynamic coefficient which decides this coefficient will be retained or removed from the model; (ii) the reference value of a hydrodynamic coefficient when the real experimental data are used for identification, instead of the simulation data.


### 3.3.2. Difference Method

Besides simplification of the mathematical model of ship manoeuvring, reconstruction of the input variables or samples provides another effective way to alleviate the parameter drift. It is noted that such a reconstruction does not require a new manoeuvring to obtain new input-output samples for system identification but modify the obtained samples. Difference method provides a simple way to reconstruct the samples [24].Taking (13) as an example, the following discrete equation can be obtained:(20)Δrk+1=Δrk+θ1ΔvkUk+θ2ΔrkUk+θ3ΔδkU2k+θ4ΔvkΔvk+θ5ΔvkΔrk+θ6ΔrkΔrk+θ7ΔvkΔrk+θ8Δδ3kU2k+θ9Δv2kΔδk+θ10ΔvkΔδ2kUk+θ11ΔrkΔδ2kUk+θ12Δr2kΔδk+θ13ΔrkΔvkΔδk+θ14U2k,where θ1~14 are the linear combinations of the hydrodynamic coefficients with respect to the sway and yaw motion [23]. Define vectors(21)Θ=θ1,θ2,…,θ1414×1,Xk=ΔvkUk,ΔrkUk,ΔδkU2k,ΔvkΔvk,ΔvkΔrk,ΔrkΔrk,ΔvkΔrk,Δδ3kU2k,Δv2kΔδk,ΔvkΔδ2kUk,ΔrkΔδ2kUk,Δr2kΔδk,ΔrkΔvkΔδk,U2k14×1T.Equation (20) can be rewritten as(22)Δrk+1=Δrk+ΘXk.At a neighbour sampling time, the above equation becomes(23)Δrk=Δrk-1+ΘXk-1.Combined with (22), the following equality holds: (24)zk=zk-1+ΘYk-1,by definition of(25)zk=Δrk+1-Δrk,Yk-1=Xk-Xk-1.Comparing (23) with (24), it can be seen that although the input-output pair are altered, the coefficient matrix remains the same. Usually, if the degree of linear dependence of input variables in the original system is serious, it could be moderated after difference of samples. As a comparison, Table 3 lists the correlation coefficients after the difference method is employed, as shown below the principal diagonal (elements above the principal diagonal are the correlation coefficients without difference method shown in Table 2). As can be seen, the linear dependence of the input variables has been obviously decreased.Table 3
Comparison of correlation coefficients using difference method.
N v ′ N r ′ N v v ′ N v r ′ N δ ′ N δ δ δ ′ N v δ δ ′ N v v δ ′ N r r ′ N v r ′ N r δ δ ′ N r r δ ′ N r v δ ′ N v ′ 1.00 0.45 0.98 0.67 0.91 0.92 0.99 0.97 0.39 0.84 0.54 0.54 0.68 N r ′ 0.44 1.00 0.42 0.84 0.53 0.54 0.47 0.43 0.97 0.66 0.91 0.75 0.81 N v v ′ 0.83 0.31 1.00 0.66 0.83 0.85 0.98 1.00 0.32 0.77 0.50 0.42 0.67 N v r ′ 0.34 0.21 0.38 1.00 0.62 0.64 0.68 0.66 0.75 0.87 0.92 0.73 1.00 N δ ′ 0.64 0.44 0.23 0.01 1.00 0.99 0.91 0.83 0.50 0.83 0.61 0.71 0.62 N δ δ δ ′ 0.50 0.33 0.25 0.13 0.74 1.00 0.92 0.85 0.50 0.83 0.62 0.69 0.64 N v δ δ ′ 0.84 0.41 0.77 0.36 0.48 0.66 1.00 0.98 0.40 0.83 0.55 0.53 0.69 N v v δ ′ 0.81 0.31 0.99 0.39 0.22 0.29 0.82 1.00 0.33 0.77 0.50 0.42 0.67 N r r ′ 0.34 0.88 0.03 0.07 0.54 0.35 0.28 0.04 1.00 0.62 0.86 0.81 0.72 N v r ′ 0.57 0.27 0.21 0.46 0.58 0.41 0.40 0.19 0.30 1.00 0.78 0.82 0.89 N r δ δ ′ 0.19 0.17 0.08 0.58 0.21 0.29 0.09 0.06 0.05 0.30 1.00 0.83 0.91 N r r δ ′ 0.29 0.32 0.15 0.16 0.72 0.44 0.09 0.18 0.49 0.62 0.50 1.00 0.73 N r v δ ′ 0.32 0.11 0.33 0.90 0.01 0.09 0.26 0.32 0.16 0.49 0.73 0.24 1.00


### 3.3.3. Additional Signal

Besides the difference method, another approach to reconstruct the input-output samples can be used, which refers to the method of additional signal. Such a method had been successfully applied to the identification of ship manoeuvring by analyzing simulation data [23] and captive model test results [54]. Taking the system (22) as an example, an additional ramp signal can be designed as(26)yak+1=yak+λ,where λ is a constant. Adding this signal to (22), a new parameter system can be obtained as(27)Δrk+1+yak+1=Δrk+yak+ΘXk+λ,where the output variable and the first input variable change; however, the parameter matrix Θ related to hydrodynamic coefficients does not. The additional signal can be viewed as an additional excitation to the original system. As a result, the manoeuvring characteristics of the original system, which might cause the parameter drift as described in the Section 3.1, are changed and the parameter drift can be alleviated resultantly. From the point of view of regression analysis, the introduction of the additional signal alters the input variables, which might alleviate the degree of linear dependence between the changed variables and others. Taking model (27) as an example, assuming λ=1, Table 4 lists the correlation coefficients before and after the additional signal is attached to the samples, in which N1=N1′=1 indicates that coefficient of the first input variable Δr(k)+ya(k) in (27) is one. As can be seen from the comparison results, the degree of linear dependence remarkably decreases for most of variables after an additional signal is added.Table 4
Comparison of correlation coefficients using additional signal.
N 1 N 1 ′ N v ′ 0.465 0.081 N r ′ 0.996 0.011 N v v ′ 0.450 0.026 N v r ′ 0.849 0.002 N δ ′ 0.552 0.028 N δ δ δ ′ 0.556 0.022 N v δ δ ′ 0.491 0.005 N v v δ ′ 0.455 0.029 N r r ′ 0.023 0.834 N v r ′ 0.970 0.039 N r δ δ ′ 0.665 0.015 N r r δ ′ 0.912 0.013 N r v δ ′ 0.771 0.024It is noted that the design of additional signal depends on the characteristic of the data samples or the manoeuvring motion. For periodically nonlinear motion such as zigzag manoeuvre or turning circle manoeuvring, the linear ramp signal like (26) works well [23]; for a linear motion such as oblique towing test, a nonlinear additional signal is more appropriate [54].


## 3.3.1. Simplification of the Manoeuvring Model

Although the Abkowitz model is often preferred in the prediction of ship manoeuvring, such a “large” parametric system makes it inconvenient to use in practice. Actually some nonlinear derivatives are difficult to be determined although they could be unimportant. In practical application of the Abkowitz model, simplification had been performed to improve the feasibility of this model. Examples are the work by Abkowitz [11] and Norrbin Nils [37]. Physical meanings were the main concern in the simplification. The simplification can also be performed from the point of view of regression analysis [24]. Taking the yaw-equation (13) as an example, the original expression of (14) is(15)Δf3′=N0u′Δu′+N0uu′Δu′2+Nv′Δv′+Nr′Δr′+Nδ′Δδ′+Nvvv′Δv′3+Nδδδ′Δδ′3+Nvvr′Δv′2Δr′+Nvvδ′Δv′2Δδ′+Nvδδ′Δv′Δδ′2+Nδu′Δδ′Δu′+Nvu′Δv′Δu′+Nru′Δr′Δu′+Nδuu′Δδ′Δu′2+Nrrr′Δr′3+Nvrr′Δv′Δr′2+Nvuu′Δv′Δu′2+Nruu′Δr′Δu′2+Nrδδ′Δr′Δδ′2+Nrrδ′Δr′2Δδ′+Nrvδ′Δr′Δv′Δδ′+N0′,in which 22 hydrodynamic coefficients are contained [40]. To simplify it, all the nonlinear terms related to the surge perturbation Δu′ are eliminated for two reasons. First, since the ship total velocity U(t) is taken as the nondimensional factor, the effect of velocity loss has been taken into account in every velocity term. Second, the linear dependence between the linear and nonlinear terms relating to Δu′ is significant, which will result in the parameter drift. For example, the third and twelfth terms on the right hand side of (15) at the kth sampling time can be written as(16)Nv′ΔvkUk+Nvu′ΔvkΔukby considering the prime nondimensional definition(17)Δu′=ΔuU,Δu˙′=Δu˙U2/L.For moderate manoeuvres, it holds that(18)Uk≈u0+Δuk.The expression (16) can be written as(19)Nv′Δvku0+Nv′ΔvkΔuk+Nvu′ΔvkΔuk.Obviously, the second term is linearly dependent on the third term. In consideration of the significance of the linear derivative Nv′, the nonlinear term Nvu′Δv(k)Δu(k) is removed. In the same way, the other six terms related to the derivatives N0uu′, Nδu′, Nru′, Nδuu′, Nvuu′, and Nruu′ can be removed from (15). Ignoring the effect of surge speed u on the change of yaw moment, N0u′Δu′ can be removed furthermore. Finally, by introducing the cross flow models to replace the third-order nonlinear terms, model (15) is simplified to model (14), with the number of hydrodynamic coefficients decreasing from 22 to 14.Sensitivity analysis provides another way to simplify the mathematical model of ship manoeuvring [49–53]. Using this method, the effect of the hydrodynamic coefficient on the prediction of ship manoeuvring is studied by observing the change of manoeuvring (e.g., velocity, angular velocity, advance, transfer, overshoots, and tactical diameter) to the change of hydrodynamic coefficients. Based on the results of sensitivity analysis, some trivial hydrodynamic coefficients are removed from the mathematical model. The main difficulties with this method are the determination of two factors: (i) the criteria of importance for a hydrodynamic coefficient which decides this coefficient will be retained or removed from the model; (ii) the reference value of a hydrodynamic coefficient when the real experimental data are used for identification, instead of the simulation data.


## 3.3.2. Difference Method

Besides simplification of the mathematical model of ship manoeuvring, reconstruction of the input variables or samples provides another effective way to alleviate the parameter drift. It is noted that such a reconstruction does not require a new manoeuvring to obtain new input-output samples for system identification but modify the obtained samples. Difference method provides a simple way to reconstruct the samples [24].Taking (13) as an example, the following discrete equation can be obtained:(20)Δrk+1=Δrk+θ1ΔvkUk+θ2ΔrkUk+θ3ΔδkU2k+θ4ΔvkΔvk+θ5ΔvkΔrk+θ6ΔrkΔrk+θ7ΔvkΔrk+θ8Δδ3kU2k+θ9Δv2kΔδk+θ10ΔvkΔδ2kUk+θ11ΔrkΔδ2kUk+θ12Δr2kΔδk+θ13ΔrkΔvkΔδk+θ14U2k,where θ1~14 are the linear combinations of the hydrodynamic coefficients with respect to the sway and yaw motion [23]. Define vectors(21)Θ=θ1,θ2,…,θ1414×1,Xk=ΔvkUk,ΔrkUk,ΔδkU2k,ΔvkΔvk,ΔvkΔrk,ΔrkΔrk,ΔvkΔrk,Δδ3kU2k,Δv2kΔδk,ΔvkΔδ2kUk,ΔrkΔδ2kUk,Δr2kΔδk,ΔrkΔvkΔδk,U2k14×1T.Equation (20) can be rewritten as(22)Δrk+1=Δrk+ΘXk.At a neighbour sampling time, the above equation becomes(23)Δrk=Δrk-1+ΘXk-1.Combined with (22), the following equality holds: (24)zk=zk-1+ΘYk-1,by definition of(25)zk=Δrk+1-Δrk,Yk-1=Xk-Xk-1.Comparing (23) with (24), it can be seen that although the input-output pair are altered, the coefficient matrix remains the same. Usually, if the degree of linear dependence of input variables in the original system is serious, it could be moderated after difference of samples. As a comparison, Table 3 lists the correlation coefficients after the difference method is employed, as shown below the principal diagonal (elements above the principal diagonal are the correlation coefficients without difference method shown in Table 2). As can be seen, the linear dependence of the input variables has been obviously decreased.Table 3
Comparison of correlation coefficients using difference method.
N v ′ N r ′ N v v ′ N v r ′ N δ ′ N δ δ δ ′ N v δ δ ′ N v v δ ′ N r r ′ N v r ′ N r δ δ ′ N r r δ ′ N r v δ ′ N v ′ 1.00 0.45 0.98 0.67 0.91 0.92 0.99 0.97 0.39 0.84 0.54 0.54 0.68 N r ′ 0.44 1.00 0.42 0.84 0.53 0.54 0.47 0.43 0.97 0.66 0.91 0.75 0.81 N v v ′ 0.83 0.31 1.00 0.66 0.83 0.85 0.98 1.00 0.32 0.77 0.50 0.42 0.67 N v r ′ 0.34 0.21 0.38 1.00 0.62 0.64 0.68 0.66 0.75 0.87 0.92 0.73 1.00 N δ ′ 0.64 0.44 0.23 0.01 1.00 0.99 0.91 0.83 0.50 0.83 0.61 0.71 0.62 N δ δ δ ′ 0.50 0.33 0.25 0.13 0.74 1.00 0.92 0.85 0.50 0.83 0.62 0.69 0.64 N v δ δ ′ 0.84 0.41 0.77 0.36 0.48 0.66 1.00 0.98 0.40 0.83 0.55 0.53 0.69 N v v δ ′ 0.81 0.31 0.99 0.39 0.22 0.29 0.82 1.00 0.33 0.77 0.50 0.42 0.67 N r r ′ 0.34 0.88 0.03 0.07 0.54 0.35 0.28 0.04 1.00 0.62 0.86 0.81 0.72 N v r ′ 0.57 0.27 0.21 0.46 0.58 0.41 0.40 0.19 0.30 1.00 0.78 0.82 0.89 N r δ δ ′ 0.19 0.17 0.08 0.58 0.21 0.29 0.09 0.06 0.05 0.30 1.00 0.83 0.91 N r r δ ′ 0.29 0.32 0.15 0.16 0.72 0.44 0.09 0.18 0.49 0.62 0.50 1.00 0.73 N r v δ ′ 0.32 0.11 0.33 0.90 0.01 0.09 0.26 0.32 0.16 0.49 0.73 0.24 1.00


## 3.3.3. Additional Signal

Besides the difference method, another approach to reconstruct the input-output samples can be used, which refers to the method of additional signal. Such a method had been successfully applied to the identification of ship manoeuvring by analyzing simulation data [23] and captive model test results [54]. Taking the system (22) as an example, an additional ramp signal can be designed as(26)yak+1=yak+λ,where λ is a constant. Adding this signal to (22), a new parameter system can be obtained as(27)Δrk+1+yak+1=Δrk+yak+ΘXk+λ,where the output variable and the first input variable change; however, the parameter matrix Θ related to hydrodynamic coefficients does not. The additional signal can be viewed as an additional excitation to the original system. As a result, the manoeuvring characteristics of the original system, which might cause the parameter drift as described in the Section 3.1, are changed and the parameter drift can be alleviated resultantly. From the point of view of regression analysis, the introduction of the additional signal alters the input variables, which might alleviate the degree of linear dependence between the changed variables and others. Taking model (27) as an example, assuming λ=1, Table 4 lists the correlation coefficients before and after the additional signal is attached to the samples, in which N1=N1′=1 indicates that coefficient of the first input variable Δr(k)+ya(k) in (27) is one. As can be seen from the comparison results, the degree of linear dependence remarkably decreases for most of variables after an additional signal is added.Table 4
Comparison of correlation coefficients using additional signal.
N 1 N 1 ′ N v ′ 0.465 0.081 N r ′ 0.996 0.011 N v v ′ 0.450 0.026 N v r ′ 0.849 0.002 N δ ′ 0.552 0.028 N δ δ δ ′ 0.556 0.022 N v δ δ ′ 0.491 0.005 N v v δ ′ 0.455 0.029 N r r ′ 0.023 0.834 N v r ′ 0.970 0.039 N r δ δ ′ 0.665 0.015 N r r δ ′ 0.912 0.013 N r v δ ′ 0.771 0.024It is noted that the design of additional signal depends on the characteristic of the data samples or the manoeuvring motion. For periodically nonlinear motion such as zigzag manoeuvre or turning circle manoeuvring, the linear ramp signal like (26) works well [23]; for a linear motion such as oblique towing test, a nonlinear additional signal is more appropriate [54].


## 4. Example Study

The abovementioned measures on how to diminish the drift of hydrodynamic coefficients, including difference method and additional signal, are employed in identifying a simplified Abkowitz model combined with free-running model tests of a KVLCC2. SVM based SI is used to obtain the hydrodynamic coefficients. The details of the mathematical model of ship manoeuvring, the SVM methodology, the ship, and the model tests can be found in [24] and relevant literature. Table 5 lists the identification results of the hydrodynamic coefficients, in which the added masses and added moment of inertia are calculated by slender-body theory, as stated in Section 2. Based on the regression model, 25°/5° zigzag manoeuvring motion is simulated and the predicted results are compared with the experimental results, as shown in Figure 2. Comparisons with the experimental results of other two manoeuvres are presented in Figure 3 (20°/5° zigzag manoeuvre) and Figure 4 (20°/10° zigzag manoeuvre), respectively. As can be seen, good agreement between prediction and experiment is achieved. It is owing not only to the use of SVM method, a powerful modeling tool, but also to the measures taken to diminish the parameter drift.Table 5
Hydrodynamic coefficients of KVLCC2 (×10-5).
Value X-coef. Value Y-coef. Value N-coef. Value X u ˙ ′ −95.4 X u ′ −93 Y v ′ −358 N v ′ −6.3 Y v ˙ ′ −1283 X v v ′ 149 Y r ′ 1980 N r ′ −74 Y r ˙ ′ 0 X r r ′ −108 Y δ ′ 448 N δ ′ −165 N v ˙ ′ 0 X δ δ ′ −105 Y v v ′ −3.5 N v v ′ 20 N r ˙ ′ −107 X v r ′ −237 Y v r ′ 30 N v r ′ −9.9 x G ′ 3486 X v δ ′ 172 Y r r ′ 679 N r r ′ −17 m ′ 1908 X r δ ′ −433 Y v r ′ 158 N v r ′ −43 I z ′ 119 X 0 ′ −48 Y δ δ δ ′ −94 N δ δ δ ′ 22 Y v δ δ ′ −3.6 N v δ δ ′ 17 Y v v δ ′ −7.3 N v v δ ′ 12 Y r δ δ ′ −23 N r δ δ ′ 22 Y r r δ ′ 134 N r r δ ′ −40 Y r v δ ′ −28 N r v δ ′ 3.9 Y 0 ′ −6.6 N 0 ′ 0.2Figure 2
Prediction of 25°/5° zigzag manoeuvre of a KCLCC2 ship.Figure 3
Prediction of 20°/5° zigzag manoeuvre of a KCLCC2 ship.Figure 4
Prediction of 20°/10° zigzag manoeuvre of a KCLCC2 ship.


## 5. Conclusions

Parameter identifiability is an important issue for the SI application to the modeling of ship manoeuvring. Some parameters in the manoeuvring model cannot be obtained using SI, while some parameters deviate from true values even if they can be derived. Several measures are proposed to deal with the problem of parameter identifiability and the validities are demonstrated by examples. To determine the added masses and added moment of inertia, slender-body theory provides a simple estimation method. To diminish the parameter drift effect, the mathematical model of ship manoeuvring motion is simplified on the basis of regression analysis; difference of the samples and additional signal also take effect. These measures can be viewed as developed approaches that are commonly used in regression analysis to diminish the multicollinearity. For instance, the simplification of a model can be viewed as a kind of pruning-variable method or PCR analysis since the main components in the model remained while some trivial variables are removed. The additional signal method can be viewed as the opposition to pruning-variable method.Future work will be devoted to the study on a more accurate, practical and universal estimation method for added masses and added moment of inertia since slender-body method gives satisfactory results only for slender ships and “mild” manoeuvres. Moreover, how to design an appropriate additional signal according to the characteristics of samples will be studied.



---

*Source: 8909170-2016-11-17.xml*