A mathematical model must be established to study the motions of ships in order to control them effectively. An assessment of the model depends on the accuracy of hydrodynamic parameters. An algorithm for the parameter identification of the coupled pitch and heave motions in ships is, thus, put forward in this paper. The algorithm proposed is based on particle swarm optimization (PSO) and the opposition-based learning theory known as opposition-based particle swarm optimization (OPSO). A definition of the opposition-based learning algorithm is given first of all, with ideas on how to improve this algorithm and its process being presented next. Secondly, the design of the parameter identification algorithm is put forward, modeling the disturbing force and disturbing moment of the identification system and the output parameters of the identification system. Then, the problem involving the hydrodynamic parameters of motions is identified and the coupled pitch and heave motions of a ship described as an optimization problem with constraints. Finally, the numerical simulations of different sea conditions with unknown parameters are carried out using the PSO and OPSO algorithms. The simulation results show that the OPSO algorithm is relatively stable in terms of the hydrodynamic parameters identification of the coupled pitch and heave motions.
1. Introduction
The model concerned with the motions of ships describes the process of movement of the response characteristics of the control input, and this is the precondition to understanding the motions of ships. Widely used system identification techniques for hydrodynamic parameters and the hydrodynamic parameters themselves can be identified by observing the data relating to the motions of ships. This directly establishes the mathematical model of the ship’s hydrodynamic parameters and motions between the state models. Classic identification algorithms, such as the maximum likelihood identification method and the prediction error used for parameter estimation are too low in terms of sensitivity and are thus inadequate as methods used in order to solve true values. Abkowitz [1] extended the Kalman Filter method to estimate the hydrodynamic parameters of a ship’s motions. Clarke et al. [2] identified the nonlinear parameters of the motions of ships by using an artificial neural network. Haddara used free attenuation from the free response-signal method to identify the roll [3, 4] and extended it to the parameter identification of sway and yaw coupled equations [5] and pitch and heave coupled equations [6]. Haddara and Xu [7] additionally put forward a ship’s longitudinal motions state as a Markov chain process and, in order to simplify the longitudinal motions equation, they used a neural network for the identification of a ship’s longitudinal motions parameters. Mahfouz and Haddara [8] mixed the classic recognition algorithm and neural network methods, putting forward a hydrodynamic parameter identification method using RDLRNNT technology, a method that appears to be fairly sound.
Bhattacharyya and Haddara [9] used artificial neural networks (ANN) and spectral analysis methods to identify the hydrodynamic derivatives in the mathematical model involving ship and marine vehicle motions. ANN has also some defects, however, such as bad generalization performance, easily falling into a local minimum. Luo and Zou [10] applied support vector machines (SVM) to identify the hydrodynamic derivatives of Abkowitz’s model from the simulated free-running model test results and then used the regressive Abkowitz model to predict zigzag tests. Zhang and Zou [11] proposed a novel method of artificial intelligence technology in the shape of support vector machines in order to estimate the hydrodynamic coefficients in the mathematical models of ship maneuvering motions. A comparison between the predicted hydrodynamic forces and the test results shows that the identified hydrodynamic mathematical model has a good generalization performance.
In recent years, with the rapid development of intelligent algorithms, a number of scholars have used the swarm intelligent algorithm, applying it to hydrodynamic parameters in order to identify problems related to the motions of ships. Through the state equation of underwater bodies and observation equation, Chen et al. [12] has used intelligent recognition technology for simulation recognition of an underwater navigation body and has obtained ten hydrodynamic parameters. In 2008, Chen et al. [13] proposed a new recognition algorithm based on intelligence technology. The least-square criterion together with the Differential Swarm Intelligent (DS) algorithm is employed to identify hydrodynamic parameters. In 2011, Dai et al. [14, 15] used an improved PSO algorithm and continuous domain ant colony optimization algorithm to identify the hydrodynamic parameters of heave and pitch with satisfactory results. In order to determine a water diving device parameter test and for the theoretical calculation of difficult problems, Gao and Li [16] put forward a method based on the basic artificial bee colony algorithm and the improved artificial bee colony method used to identify a potential coefficient method. Experiments show that the use of artificial colony algorithm identification device parameters is indeed feasible.
The PSO algorithm was proposed by Kennedy and Eberhart in 1995 [17, 18]. In order to overcome the existing problems in the practical use of the PSO algorithm and improve the performance of the algorithm, an improved algorithm has been put forward by number of scholars [19–26]. Generally speaking, the improved PSO algorithm strategy currently includes two aspects, namely, the improvement of the strategy of the PSO algorithm and its fusion with other algorithms. The improved PSO algorithm mainly concerns the variation of the particles, the multipopulation cooperation, and the design of parameters. These methods can potentially prevent particle aggregation and conflict, and avoid premature convergence to local optima. PSO algorithms are integrated with other algorithms, improving the strategy of the PSO. Some scholars are currently advocating the opposition-based learning particle swarm optimization algorithm. Wang et al. [21] have introduced opposition-based learning into the PSO algorithm, then they proposed opposition-based learning and the Cauchy mutation PSO algorithm (OPSO), using opposition-based learning to initialize the group. Omran and Al-Sharhan [24] have used dynamic shrinkage factors to generate opposition-based solutions. Shahzad [25] have presented three kinds of opposition-based PSO algorithm: the first version of the OPSO algorithm using only opposition-based learning in order to initialize groups and the second version (IOPSO) in addition to opposition-based learning produced the opposition-based particles as a replacement for the worst particle in the group during each iteration. In the third version of the algorithm the initialization of opposition-based learning is removed from IOPSO. Shi and Eberhart [26] have controlled the velocity of PSO by using opposition-based learning and proposed the speed clamping PSO algorithm (OVCPSO) based on opposition-based learning, which achieved good results.
Owing to the limitations of conventional identification methods in coupled pitch and heave motions parameter identification, in this paper, we are proposing to identify ship coupled heave and pitch motions using opposition-based PSO. In order to achieve this, we have designed a model involving wave disturbance force and torque disturbance, using the design methods for an output parameter identification system, with an opposition-based PSO algorithm for parameter identification.
2. Pitch and Heave Motions Model
The coupled pitch and heave motions of a ship in a realistic sea can be described by two linear second order ordinary coupled differential equations as follows [7]:
(1)(a33+m0)z¨+b33z˙+cz33+a35θ¨+b35θ˙+c35θ=FR+F3a53z¨+b53z˙+c53z+(I5+a55)θ¨+b55θ˙+c55θ=FRXR+M5,
where z is heave; θ is pitch; I5 is the pitching moment of inertia; m0 is the ship’s mass; FR is the force of the rudder; XR is advanced from the rudder lifting center to the ship’s center of gravity; F3 is the disturbance force of heave; M5 is the disturbance moment of pitch, a33 is the added mass, a35 and a53 are mass moment; a55 is the moment of inertia; b33, b35 are the damping coefficients; c33, c35 are the resilience coefficient; b53, b55 are the damping moment coefficient; c53, c55 are the righting moment coefficients; and a33, a35, a53, a55, b33, b35, b53, b55, c33, c35, c53, and c55 are the hydrodynamic parameters.
Set x1=z, x3=θ as the state variables
(2)X=[x1x2x3x4]T=[zz˙θθ˙]T.
The system state equation is then obtained as follows:
(3)X˙=AX+Bu+CW.
Among A=E-1A*, B=E-1B*, C=E-1C*, where
(4)E=(10000a33+Δ0a3500100a530I5+a55),A*=(0100-c33-b33-c35-b350001-c53-b53-c55-b55)B=[010XR],u=FR,C=[00100001],W=[F3M5]
para =[a33,a35,a53,a55,b33,b35,b53,b55,c33,c35,c53,c55] as the parameters to be identified, W is the wave disturbing force and moment.
PSO learns from the natural phenomenon of birds looking for food made by a class of population-based stochastic global optimization techniques. With respect to birds in flight, in its initial state every bird is in a random position and flies randomly in all directions, but as time goes on, these initial random state birds form a small community, through mutual learning, mutual tracking, and self-organization, fly at the same speed in the same direction, and ultimately the entire group gathers in one place, namely, the food source.
In the PSO algorithm, each individual is called a “particle,” and each particle represents a potential solution. In continuous spatial coordinates, the PSO algorithm is described as follows.
Suppose that the size of the swarm is N and the search space is D-dimensional, then the position of the ith particle is presented as Xi=(Xi1,Xi2,…,XiD), the velocity of this particle is presented as Vi=(Vi1,Vi2,…,ViD), the T fitness value of each particle in its current position is fitnessi=fitness(Xi), its corresponding optimal value is Pi=(Pi1,Pi2,…,PiD) and the population current optimal experience value is recorded as Pg=(Pg1,Pg2,…,PgD). Each particle adjusts its speed dynamically according to the comprehensive analysis of both individual and population flying experience and flies to the best position that it and other particles have experienced. Each particle updates its speed and position according to the formula equations:
(5)Vi(t+1)=wVi(t)+c1r1(Pi(t)-Xi(t))+c2r2(Pg(t)-Xi(t)),(6)Xi(t+1)=Xi(t)+Vi(t+1),
where t is iterative times, d=1,2⋯D. r1,r2 are random numbers between 0 and 1 and C1, C2 are nonnegative constants. This is called the learning factor and each iteration step is justified accordingly.
3.2. Opposition-Based Learning
Opposition-based learning was put forward by Professor Tizhoosh [19] in 2005. He argues that intelligent algorithms are based on a random guess value concerning the initial population, with each generation coming close to the solution with the optimal solution or a close approximation of the optimal solution eventually being found. Thus, the initial guess value greatly influences the algorithm and, if the random guess value is very close to the optimal solution, the algorithm may converge quickly, but if the value is far away from the solution or is even the reverse, the algorithm will take much more time. If the current solution and its opposite are sought simultaneously, a better solution will be chosen and the efficiency of the algorithm will be greatly enhanced. According to the theory of probability, there is a 50% probability that the current solution is more remote from the optimal solution than its opposite [20].
Definition 1 (definition of the opposite number).
Let x∈R be a real number within a defined interval, where x∈[a,b]. The opposite number x0 can be defined as
(7)x0=a+b-x.
Definition 2 (definition of the opposite point).
In the high dimensional space, if p=(x1,x2,…,xD) is a set of points in the D-dimensional search space where x1,x2,…,xD∈R, xj∈[aj,bj], then the points in the opposition set po=(x1o,x2o,…,xDo) can be defined as
(8)xjo=aj+bj-xj.
Definition 3 (opposition-based optimization).
For a point in the D-dimensional space X=(x1,x2,…,xD), suppose that f(X) is the function used to measure the performance of a candidate solution, according to the opposition theorem, Xo=(x1o,x2o,…,xDo) will be the opposition set for X=(x1,x2,…,xD), If f(X0)<f(X), then the set of points X can be replaced by X0, or else X is maintained.
3.3. Opposition-Based Particle Swarm Optimization
In PSO, each particle adjusts its search direction on the basis of the optimum location of all particles. In the initial stage, the algorithm converges quickly but slows down later on or even stops. These particles lose the ability to evolve when the speed of all the particles approaches zero and the algorithm is thought to represent convergence. Sometimes the algorithm does not converge to global extreme values, however, not even local extreme values. This is because the high aggregate and deficiency diversity of the particles takes a long time or an infinite time to skip from the focusing point.
In order to solve this problem, the opposition-based learning mechanism is introduced into the basic PSO, and a new random optimization algorithm is constructed, the opposition-based particle swarm optimization (OPSO) algorithm. In the OPSO algorithm, a variable is set and referred to as conNum. If the global best fitness is not updated during a single iteration, then the conNum = conNum + 1, and when conNum reaches a constant set number setNum, it shows that there is a high concentration of particles, and the algorithm cannot find a better solution in the current position and speed. At this point, opposition-based learning is brought into play. According to (8) the position-based position Xo of the particle’s current position is produced, the smaller values N are selected from X and Xo which have 2N locations and a new swarm is then formed. If the conNum is less than setNum, then according to (5) and (6) the particle velocity and position are updated. The process of OPSO is shown in Algorithm 1.
Algorithm 1: The process of OPSO.
%% Initialization
(1) Initialize swarm size N and constant number C1 and C2; space dimension D; maximum
(2) Initialize the iteration number for opposition calculation conNum=0, setnum=100;
(3) for i=1: particleNum
(4) initialize Xi randomly with the search range (Xmin,Xmax)
(5) initialize Vi randomly with the velocity range (Vmin,Vmax)
(6) End for
(7) Evaluate each particle’s fitness fi=fitness(Xi), and the best fitness PbestValuei=fi
(8) Identify the best particle’s position Pg and its fitness GbestValue=min(PbestValuei); %%Loop
(9) While (t<nmax)
(10) If conNum>setNum then
(11) for i=1 to N
(12) Calculate the opposite particle X0 using (8);
(13) Evaluate fitness in opposite vectorfio=fitness(Xio);
(14) end for
(15) Select N fittest particles Newfitness from F(f1,f1,…,fN) and F0(f1o,f1o,…,fNo) to
(16) create a population of size N;
(17) Else
(18) for i=1 to N
(19) Calculate particle velocity Vi(t+1) using (5)
(20) update particle position Xi(t+1) using (6)
(21) Evaluate newfitnessi=fitness(Xi(t+1))
(22) End for
(23) End if
(24) update Pbesti
(25) for i=1 to N
(26) if (newFitness<GbestValue)
(27) update Gbest
(28) conNum=0;
(29) end if
(30) end for
(31) conNum=conNum+1;
(32) end while
4. Design of Pitch and Heave Parameter Identification
The value of pitch and heave can be measured, so select state x1 and x3 in (2) as the measurement state, with the observation equation as
(9)Y=HX+V,H=(10000010),
where Y is an observation vector and V is the two-dimensional measurement of noise, and this is generally considered to be white noise.
Before computer simulation, it was necessary to disperse the state equation and observation equation. This paper does not consider the steering input but only the wave disturbance on the ship, with the state equation being formulated as follows:
(10)Xk+1=𝚽Xk+ΓWk,Yk+1=HXk+1+Vk+1,
where, 𝚽=eATs, Γ=∫0TseAtdtC.
In this paper, the single-parameter ITTC spectrum is chosen for wave disturbance simulation, with the expression being formulated as follows:
(11)Sζ(ω)=8.1×10-3g2ω5exp(-3.11h1/32ω4),
where g is the acceleration of gravity, w is the natural angular frequency, and h1/3 is a third significant wave height, for a level 4 and 5 sea condition, and the ITTC recommended value is 2.5 m and 3.75 m, respectively.
By (10), the observed quantity is the function of X, para, namely:
(12)yi=g(xi,parai,Wi)+Vi,i=1,2,…,dataNum,
where dataNum is the times of observation, g(·) is the model output vector, and yi is the observation vector. The minimum mean square errors between the observed data and identical value are setting as follows:
(13)Fj=sqrt(1dataNum∑i=1dataNum(yi-g(xi,parai,Wi))2),111111111111111111111111111.11111111111j=1,2.
The hydrodynamic parameters of our requirements should mean that the two components in (13) obtain their minimum value simultaneously under constrained conditions (10). We must set different weight coefficients to ensure that each parameter has the same effect on the objective function, taking into consideration the different order of magnitude of pitch and heave. The fitness function is shown as follows:
(14)Fit=F1+αF2,
where α is the weight coefficient.
5. Experiment Testing and Comparisons
Here the experimental ship’s parameters are as follows: the hull quality m is 442000 kg, the beam B is 7.2 m, and the draft is 2.25 m. The waterline Lpp is 60 m. I5=(0.25*LPP), m0=99450000 (kg·m2), XR=25.8 m, and ZR=1.32 m. The number of rudder n is 2 and the hydrostatic resilience factor C44 is 3370000. For the ship, C33, C35, C55, and C53 are constants and can be obtained by calculating the structural parameters of the ship. V is Gaussian white noise V is Gaussian white noise and its covariance matrix can be taken as so that the accuracy of the sensor can be taken as QVV=diag[20.3e-42.25e-6].
A comparative algorithm is the OPSO algorithm proposed in this paper and the PSO-w (PSO with inertia weight [26]), with the parameters being set as shown in Table 1.
Parameters and their range of values used in our proposed algorithm.
Parameters
Value/range
Population size of swarm
100
Acceleration constants c1,c2
c1=c2=1.4962
Inertia weight w
[0.4, 0.9]
Random number r1,r2
[0, 1]
Maximum iteration times nmax
1000
Range of velocity Vmax
(xMax–xMin)/2
setNum for OPSO
100
Sampling dot number dataNum
200
The identification results by PSO and OPSO are shown in Tables 2, 3, and 4, for level 4 sea condition, speed 18 kn and course angle 90°, 135°, and 180°, respectively.
Parameter identification of pitch and heave motions for level 4 sea condition, ship speed 18 kn and course angle 90°.
Parameters
Theoretical values
PSO-w
OPSO
Identification value
Relative error
Identification value
Relative error
a33
8.045E+05
9.446E+05
17.42%
8.141E+05
1.19%
b33
9.369E+05
7.901E+05
15.66%
9.447E+05
0.83%
a35
7.753E+06
8.902E+06
14.82%
8.037E+06
3.66%
b35
1.159E+07
1.480E+07
27.72%
1.222E+07
5.39%
a55
2.595E+08
2.830E+08
9.05%
2.706E+08
4.27%
b55
2.885E+08
4.077E+08
41.31%
2.892E+08
0.26%
a53
1.479E+07
1.777E+07
20.17%
1.500E+07
1.39%
b53
3.324E+06
2.296E+06
30.91%
3.068E+06
7.70%
Parameter identification of pitch and heave motions for level 4 sea condition, ship speed 18 kn and course angle 135°.
Parameters
Theoretical values
PSO-w
OPSO
Identification value
Relative error
Identification value
Relative error
a33
6.311E+05
5.569E+05
11.77%
6.177E+05
2.12%
b33
8.033E+05
7.635E+05
4.96%
8.245E+05
2.64%
a35
3.281E+06
5.248E+06
59.95%
3.057E+06
6.84%
b35
9.429E+06
7.879E+06
16.44%
9.598E+06
1.80%
a55
1.525E+08
1.541E+08
1.01%
1.533E+08
0.52%
b55
2.091E+08
2.716E+08
29.88%
2.115E+08
1.14%
a53
7.865E+06
9.509E+06
20.90%
7.772E+06
1.18%
b53
2.267E+06
1.616E+06
28.72%
2.423E+06
6.87%
Parameter identification of pitch and heave motions for level 4 sea condition, ship speed 18 kn and course angle 180°.
Parameters
Theoretical values
PSO-w
OPSO
Identification value
Relative error
Identification value
Relative error
a33
6.101E+05
5.815E+05
4.69%
6.148E+05
0.77%
b33
7.400E+05
7.782E+05
−5.16%
7.366E+05
0.46%
a35
2.663E+06
2.396E+06
10.03%
2.758E+06
3.56%
b35
8.951E+06
8.835E+06
1.30%
8.873E+06
0.87%
a55
1.377E+08
1.348E+08
2.14%
1.367E+08
0.74%
b55
1.873E+08
1.944E+08
−3.79%
1.838E+08
1.87%
a53
6.611E+06
6.594E+06
0.25%
6.568E+06
0.65%
b53
2.356E+06
3.123E+06
−32.57%
2.113E+06
10.32%
As we can see from Table 2 to Table 4, the parameters of pitch and heave are correctly identified by using the PSO algorithm, with the OPSO algorithm clearly obtaining better results than the PSO, especially for level 4 sea condition, speed 18 kn, and course angle 135°, and the relative error of the OPSO is smaller than for the PSO algorithm.
From the identification of the parameters of pitch and heave motions carried out by using the OPSO algorithm, we can obtain the mathematical model for pitch and heave which is constructed according to the identification of the hydrodynamic parameters. The model in which the pitch and heave for sea condition 4, speed 18 kn, and course angle 90°, 135°, and 180° is shown in Figures 1, 2, and 3, respectively. There are three curves in each diagram; the first curve “—□—” represents the observed values of the coupled pitch and heave motions, the second curve “—▼—” represents the model constructed by PSO-w, and the last curve “—●—” represents the model output values constructed by OPSO.
Pitch and heave observed value. Model identified by PSO-w and model identified by OPSO for level 4 sea condition, ship speed 18 kn and course angle 90°.
Pitch
Heave
Pitch and heave observed value. Model identified by PSO-w and model identified by OPSO for level 4 sea condition, ship speed 18 kn and course angle 135°.
Pitch
Heave
Pitch and heave observed value. Model identified by PSO-w and model identified by OPSO for level 4 sea condition, ship speed 18 kn and course angle 180°.
Pitch
Heave
It is obvious in these figures that there was an agreement between the pitch and heave model estimated by the PSO-w model and observed values for a level 4 sea condition, ship speed 18 kn, and course angle 90° and 180°, but when the course angle was 90°, the results were not accurate. However, the model estimated by the OPSO algorithm tallies completely with the observed values.
Tables 5, 6, and 7 and Figures 4, 5, and 6 show similar results for a level 5 sea condition. It is clear in Tables 5–7 that the parameters identified by the OPSO algorithm are close to the real ones while those identified by PSO-w are not, especially the results for level 5 sea condition, ship speed 18 kn, and course angle 90 degrees, as the parameters have greater relative errors. Because of this, the corresponding pitch and heave models do not tally with the observed values. However, for other course angles, both the PSO-w and OPSO algorithms agree with observed values.
Parameter identification for pitch and heave motions for level 5 sea condition, ship speed 18 kn and course angle 90°.
Parameters
Theoretical values
PSO-w
OPSO
Identification value
Relative error
Identification value
Relative error
a33
9.375E+05
1.178E+06
25.64%
9.111E+05
2.82%
b33
9.089E+05
8.644E+05
4.89%
9.635E+05
6.01%
a35
1.137E+07
2.102E+07
−84.86%
1.030E+07
9.39%
b35
1.274E+07
1.831E+07
43.71%
1.295E+07
1.65%
a55
3.613E+08
4.017E+08
11.19%
3.532E+08
2.25%
b55
3.202E+08
4.593E+08
43.43%
3.428E+08
7.07%
a53
1.942E+07
2.430E+07
25.14%
1.964E+07
1.11%
b53
4.641E+06
3.829E+06
17.50%
5.348E+06
−15.24%
Parameter identification of pitch and heave motions for level 5 sea condition, ship speed 18 kn and course angle 135°.
Parameters
Theoretical values
PSO-w
OPSO
Identification value
Relative error
Identification value
Relative error
a33
7.202E+05
7.515E+05
−4.34%
7.106E+05
1.34%
b33
8.553E+05
8.283E+05
3.15%
8.921E+05
4.31%
a35
5.539E+06
5.923E+06
−6.93%
5.426E+06
2.05%
b35
1.046E+07
9.616E+06
8.10%
1.043E+07
0.35%
a55
2.051E+08
2.086E+08
−1.70%
2.079E+08
1.34%
b55
2.451E+08
2.473E+08
−0.90%
2.467E+08
0.64%
a53
1.124E+07
1.115E+07
0.82%
1.167E+07
3.85%
b53
2.876E+06
2.297E+06
20.13%
2.895E+06
0.64%
Parameter identification of pitch and heave motions for level 5 sea condition, ship speed 18 kn and course angle 180°.
Parameters
Theoretical values
PSO-w
OPSO
Identification value
Relative error
Identification value
Relative error
a33
6.787E+05
6.812E+05
−0.38%
6.762E+05
0.37%
b33
8.182E+05
7.952E+05
2.82%
7.939E+05
2.97%
a35
4.428E+06
3.406E+06
23.09%
4.638E+06
4.74%
b35
9.914E+06
9.902E+06
0.12%
9.849E+06
0.66%
a55
1.783E+08
1.798E+08
−0.89%
1.742E+08
2.30%
b55
2.246E+08
2.317E+08
−3.17%
2.284E+08
1.70%
a53
9.511E+06
9.738E+06
−2.38%
9.351E+06
1.69%
b53
2.664E+06
2.228E+06
16.37%
2.974E+06
11.62%
Pitch and heave observed value. Model identified by PSO-w and model identified by OPSO for level 5 sea condition, ship speed 18 kn and course angle 90°.
Pitch
Heave
Pitch and heave observed value. Model identified by PSO-w and model identified by OPSO for level 5 sea condition, ship speed 18 kn and course angle 135°.
Pitch
Heave
Pitch and heave observed value. Model identified by PSO-w and model identified by OPSO for level 5 sea condition, ship speed 18 kn and course angle 180°.
Pitch
Heave
6. Conclusions
The identification of the hydrodynamic parameters of ships is an important way of obtaining these parameters. In this paper, we have used OPSO to design the method of identifying the hydrodynamic parameters of the coupled pitch and heave motions of ships. This paper introduces in detail the opposition-basedlearning algorithm and puts forward an improved idea and process for the opposition-based algorithm. In addition, this paper introduces the process involving the hydrodynamic parameter identification algorithm. Here, we have established wave disturbance as the model input, with the algorithm’s fitness function being the output model. The hydrodynamic parameter identification problem was then converted into a constrained optimization problem and the OPSO algorithm was used to find the optimal solution. Finally, we made use of computer simulation, with the simulation results showing that the OPSO algorithm is relatively stable in terms of identifying the hydrodynamic parameters connected with the problem of the coupled pitch and heave motions of ships. In addition, the identified coupled pitch and heave model values and the observed values are consistent. This method may provide a new solution for the identification of coupled pitch and heave motions. This paper has not taken the disturbance of the rudder angle into consideration, something which the author intends to research in due course.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work has been financially supported by the National Natural Science Foundation of China under the Grant 51109041 and the fundamental research funds for the central universities.
AbkowitzM. A.Measurement of hydrodynamic characteristics from ship manoeuvring trials by system identification198088283318ClarkeD. P.GedlingP.HineG.The application of manoeuvring criteria in hull design using linear theory198212545682-s2.0-0020721324HaddaraM. R.On the use of neural network techniques for the identification of ship stability parameters at sea2Proceedings of the 14th International Conference on Offshore Mechanics and Arctic Engineering19951271352-s2.0-0029225227HaddaraM. R.On the random decrement for nonlinear rolling motion2Proceedings of the 12th International Conference on Offshore Mechanics and Arctic Engineering19922832882-s2.0-0027091585HaddaraM. R.WangY.Parametric identification of coupled sway and yaw motions1Proceedings of the 15th International Conferece on Offshore Mechanics and Arctic Engineering1996267273HaddaraM. R.XuJ. S.On the use of random decrement in the identification of two degrees of freedom systems4Proceedings of the Canadian Society for Mechanical Engineering Forum (CSME '98)1998Ryerson Polytechnic University499507HaddaraM. R.XuJ.On the identification of ship coupled heave-pitch motions using neural networks19982653814002-s2.0-034460980310.1016/S0029-8018(98)00005-5MahfouzA. B.HaddaraM. R.Effect of the damping and excitation on the identification of the hydrodynamic parameters for an underwater robotic vehicle2003308100510252-s2.0-003740899310.1016/S0029-8018(02)00080-XBhattacharyyaS. K.HaddaraM. R.Parametric identification for nonlinear ship maneuvering20065031972072-s2.0-33749494099LuoW. L.ZouZ. J.Parametric identification of ship maneuvering models by using support vector machines200953119302-s2.0-68949108705ZhangX. G.ZouZ. J.Estimation of the hydrodynamic coefficients from captive model test results by using support vector machines201373253110.1016/j.oceaneng.2013.07.007ChenW. Q.YanK.ShiG. J.WangS. T.LiuZ. Y.Research of hydrodynamic parameter identification for underwater vehicle using swarm intelligence algorithm200711140462-s2.0-34047226368ChenW. Q.YanK.WangB. S.Intelligent identification of hydrodynamic parameters of navigating body2008122204210DaiY. T.ZhaoX. R.LiuL. Q.Parameter identification of ship longitudinal motions using continuous domains ant colony algorithm201026162324DaiY. T.ZhaoX. R.LiuL. Q.Parameter identification of ship vertical motions using improved particle swarm optimization2010141-244502-s2.0-77949888216GaoM. S.LiS. M.Identification on hydrodynamic coefficients of underwater vehicle with the ABC algorithm201334815KennedyJ.EberhartR. C.Particle swarm optimizationProceedings of the IEEE International Conference on Neural Networks1995194219482-s2.0-0029535737EberhartR.KennedyJ.A new optimizer using particle swarm theoryProceedings of the 6th International Symposium on Micro Machine and Human Science199539432-s2.0-0029517385TizhooshH. R.Opposition-based learning: a new scheme for machine intelligence1Proceedings of the International Conference on Computational Intelligence for Modelling, Control and Automation (CIMCA '05)2005Vienna, Austria6957012-s2.0-3384722147910.1109/CIMCA.2005.1631345RahnamayanS.TizhooshH.R.SalamaM. M.Opposition-based differential evolution200812164792-s2.0-4024908470610.1109/TEVC.2007.894200WangH.LiH.LiuY.LiC.ZengS.Opposition-based particle swarm algorithm with cauchy mutationProceedings of the IEEE Congress on Evolutionary Computation (CEC '07)2007Singapore475047562-s2.0-7795609408010.1109/CEC.2007.4425095KennedyJ.Particle swarm: social adaptation of knowledgeProceedings of the IEEE International Conference on Evolutionary Computation19973033082-s2.0-0030645460LinH.XingshiH.A novel opposition-based particle swarm optimization for noisy problemsProceedings of the 3rd International Conference on Natural Computation (ICNC '07)2007Haikou, China6246292-s2.0-3804902261410.1109/ICNC.2007.119OmranM. G. H.Al-SharhanS.Using opposition-based learning to improve the performance of particle swarm optimizationProceedings of the IEEE Swarm Intelligence Symposium (SIS '08)2008Saint Louis, Mo, USA162-s2.0-5764917820110.1109/SIS.2008.4668288ShahzadF.BaigA. R.MasoodS.KarmanM.NaveedN.Opposition-based particle swarm optimization with velocity clamping (OVCPSO)2009116339348Advances in Intelligent and Soft Computing10.1007/978-3-642-03156-4_34ShiY.EberhartR. C.A modified particle swarm optimizerProceedings of the IEEE International Conference on Evolutionary Computation, IEEE World Congress on Computational Intelligence1998Anchorage, Alaska, USA697310.1109/ICEC.1998.699146ZBL1191.68833