# Heave Motion Measurement by Adaptive Filter Based on Longuet-Higgins Wave Model

Authors: Jiazhen Lu; Qiuwei Luo; Yanqiang Yang
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

A method is proposed to obtain heave motion information based on the Longuet-Higgins wave model. The Longuet-Higgins wave model which is closer to the sea wave is introduced. Based on it, random process of the noise is analyzed and the highpass filter is designed to reduce errors. Then it is the key point in this article that an adaptive algorithm is put forward because of the complexity of the waves. The algorithm adjusts the cutoff frequency to reduce the amplitude attenuation of the filter by analyzing the wave. For the same reason the comprehensive parameter of the phase compensation can be also obtained by the algorithm. Simulation measurement results show that under the rough sea situation the maximum value of absolute error is 0.4942 m according to the normal method, the method is 0.1170 m, and the average error ratio of the rough sea test reduces to 3.89% from 12.54%, which demonstrates that the adaptive filter is more effective in measuring heave movement. A variety of simulation cases show that the adaptive filter can also improve the precision of the heave motion under different sea situations.


---

## Body

## 1. Introduction

Heave motion is generally used in the research of floating offshore engineering. For example, the introduction of heave information is beneficial for ship replenishment, landing of shipboard aircraft, and the design of heave compensation of offshore oil drilling platform. Moreover, the measurement of heave motion is used to provide a reference for the safety of vessels sailing. For example, “sheng lu” ship capsized due to the factor of the storm in 1991.How to get heave information has been mentioned in many studies. Basically there are two main types of filtering methods for heave motion. One is the Kalman filter. For example, El-Hawary estimated heave component with Kalman filtering which established a handled value [1]. Triantafyllou et al. used Kalman filtering techniques to estimate the ship motion [2]. Aminzadeh and El-Hawary used the Generalized Kalman to solve the problem of heave compensation in underwater applications [3]. Fossen and Perez introduced the ship motion-control models and used Kalman filter to eliminate the oscillatory [4]. But the method has two main drawbacks. One drawback is that the model of the heave dynamics cannot be equivalent to the actual model. The other drawback is that Kalman filtering requires GPS or other external reference information. However, the result of GPS height is not accurate enough to compute heave motion value. In this paper according to the requirment of navigation system design, a long time pure inertial navigation is desired, so Kalman filter is not applicable in this case.The other type is the highpass filter. For example, in the 1990s, an adaptive filter algorithm for heave motion was developed by Godhavn and was implanted in the seatex MRU [5]. Yang et al. did active heave compensation experiments to validate the validity of active heave compensation system in 2008 [6]. Simulations and the tests with MTi unit were done by Yang et al. and the results have proven that the adaptive highpass filter was effective in 2009 [7]. Wei and Feng used highpass filter to eliminate Schuler errors to get heave information in 2012 [8]. Hu and Tao put forward the weighted-frequency Fourier linear combiner algorithm for phase compensation in 2013 [9]. However, little attention has been focused on the sea wave model which decides the cutoff frequency of the filter. It is often assumed that the wave model is a simple sine function. In fact, ship motion in waves is a complex random process as a result of sea tides, currents, and wind effects.In order to get accurate heave motion which is based on strapdown inertial navigation system (SINS) composed of inertial sensors with high accuracy, it is necessary to select a sea model which can better reflect the actual situation to analysis. Thus, this paper chooses the Longuet-Higgins wave model which is closer to the actual situation and analyzes the random process of the narrowband noise. On the basis of this analysis, an algorithm is proposed to adjust the cutoff frequency to reduce the amplitude error caused by the filter. In addition, the same algorithm is used to change all pass system parameters to compensate phase error.The structure of this paper is given as follows. In Section2 the Longuet-Higgins wave model and the accelerometer error model are introduced to analyze the noise. In Section 3, the method is presented and analyzed. Moreover, the steps of heave measurement are shown. Simulation measurement results are given in Section 4. Finally, the article ends with concluding remarks in Section 5.


## 2. Review and Analysis

Referring to previous literatures, the sea model and the inertial unit error model [10–12] are set up. In addition, the correlated noise is analyzed based on the model.


### 2.1. Mathematical Model

#### 2.1.1. Longuet-Higgins Wave Model

In order to describe the characteristics of the actual wave, Longuet-Higgins wave model is established on the basis of the measurement data [13]. According to the assumption of the model, the wave motion is a random process. The fixed point fluctuation of the wave is composed of some different phases and different amplitude sine wave.The wave surface displacement is expressed as(1)ηt=∑n=1∞ancos⁡wnt+εn.In the equation above, an and wn are the amplitude and frequency of the nth wave component, respectively, t represents time, and εn is the random phase uniformly distributed in the (0,2π). Furthermore, the sea wave only has a dominant component under most condition.


#### 2.1.2. Inertial Unit Error Model

The accelerometer error model is given as follows [14]:(2)Ea=n+k0+kiGi+kiiGi2+kiiiGi3+kioGiGo+kipGiGp,where Ea is the accelerometer error, n is the white Gaussian noise, k0 is a constant bias, ki is linear error coefficient, kii, kiii are, respectively, second-order and third-order nonlinear error coefficient, and kio, kip are cross-axis coupling coefficient, respectively. Gi, G0, and Gp are, respectively, the specific force of input axis, output axis, and pendulum axis.The relationship of two models is given as follows:(3)η¨=T·A-Ea,where A is the measured value of accelerometers; T is coordinate transformation matrix of SINS.On the other hand, the gyroscope error contains the constant error and random drift. Affected by inertial unit errors, the coordinate transformation matrixT is not accurate. When SINS works for a long time, the 6-degree-of-freedom navigation results are in oscillation because of incorrect T. According to the error analysis of SINS, under the undamped case the results include three kinds of different cycles of persistent oscillation [15]. They are, respectively, Schuler oscillation period Ts, earth oscillation period Te, and Foucault oscillation period TF.(4)Ts=2πRg≈84.8min,Te=2πΩ=24h,TF=2πΩsin⁡φ,where R is the earth radius; Ω is self-rotation angle velocity of the earth; φ is latitude; g is gravitational acceleration. Foucault oscillation period decreases with increased latitude. In addition, the period of earth oscillation and the period of Foucault oscillation are greater than the period of the heave motion, so this paper mainly considers the effects of Schuler oscillation for heave motion measurement.


#### 2.1.3. Highpass Filter

Schuler oscillation whose period is 84.8 min in the process of SINS is caused by the constant bias of inertial sensors with high accuracy. The existence of Schuler oscillation makes heave information divergent. The frequency of Schuler oscillation is less than the frequency of heave motion. The highpass filter is designed to eliminate Schuler oscillation and reduce noise [16].In this paper a fourth-order highpass filter is given by the following equation:(5)Hs=η∬Adtdt=s4s2+2ξwcs+wc22,where ξ=1/2 is the damping factor, s is the Laplace variable, and wc is the cutoff frequency of the highpass filter which can be tuned.


### 2.2. Analysis of Narrowband Noise

Due to the influence of random phase of the wave, the high order nonlinear error part sensed by the accelerometer can be equivalent to a narrowband noise after two integrations [17]. The white Gaussian noise is also in the form of narrowband noise in the heave information because of the output with a range of sampling frequency. The random process of narrowband noise in the wave model is discussed as follows.Without considering the impact of other errors which are the constant bias and linear errors of the accelerometer, narrowband noise in the random phase wave model is analyzed. The model can be expressed as follows:(6)Yt=∫∫Adtdt≈ηt+Nt,where ηt is the sine wave with the random phase(7)ηt=acos⁡w0t+β=acos⁡βcos⁡w0t-asin⁡βsin⁡w0t,where a is the known amplitude, β is the random phase uniformly distributed in the (0,2π), and w0 is the known angular frequency.In (5), narrowband noise Nt is a narrowband stationary Gaussian random process. The mean is zero and the variance is σ2. It is given as follows:(8)Nt=Qtcos⁡w0t+φnt=Qctcos⁡w0t-Qstsin⁡w0t.According to (6) and (7), (5) is rewritten as(9)Yt=Bctcos⁡w0t-Bstsin⁡w0t=Btcos⁡w0t+Φt.In the equation above,(10)Bct=acos⁡β+Qct,Bst=asin⁡β+Qst,Bt=Bc2t+Bs2t1/2,Φt=arctan⁡QstQct.Under the condition of a given signal phase θ and the time of t, with the value range of B and θ being B≥0, θ≤2π, the probability density function of the envelope is given by(11)pB∣θ=Bσ2exp⁡-a2+B22σ2I0aBσ2B≥0,where I0aB/σ2 is the zero-order Bessel function and aB/σ2 is the number of Bessel functions.Assuming thatb=a/σ, v=B/σ, (10) can be expressed as(12)pv=vσexp⁡-12v2+b2I0bvv≥0.In the actual situation, the useful signal of the accelerometer is much larger than the noise; therefore, bv≥1(13)pv≈vσexp⁡-12b2+v21vb2πexp⁡bv≈1σ2πexp⁡-12v-b2.With B≈a and v=b, the mean and variance of Bt are given by(14)EBt=a,DBt=σ2.However, some of frequency of the noise is the same as frequency of heave motion. So the important thing is to select the appropriate cutoff frequency of the filter to minimize errors.


## 2.1. Mathematical Model

### 2.1.1. Longuet-Higgins Wave Model

In order to describe the characteristics of the actual wave, Longuet-Higgins wave model is established on the basis of the measurement data [13]. According to the assumption of the model, the wave motion is a random process. The fixed point fluctuation of the wave is composed of some different phases and different amplitude sine wave.The wave surface displacement is expressed as(1)ηt=∑n=1∞ancos⁡wnt+εn.In the equation above, an and wn are the amplitude and frequency of the nth wave component, respectively, t represents time, and εn is the random phase uniformly distributed in the (0,2π). Furthermore, the sea wave only has a dominant component under most condition.


### 2.1.2. Inertial Unit Error Model

The accelerometer error model is given as follows [14]:(2)Ea=n+k0+kiGi+kiiGi2+kiiiGi3+kioGiGo+kipGiGp,where Ea is the accelerometer error, n is the white Gaussian noise, k0 is a constant bias, ki is linear error coefficient, kii, kiii are, respectively, second-order and third-order nonlinear error coefficient, and kio, kip are cross-axis coupling coefficient, respectively. Gi, G0, and Gp are, respectively, the specific force of input axis, output axis, and pendulum axis.The relationship of two models is given as follows:(3)η¨=T·A-Ea,where A is the measured value of accelerometers; T is coordinate transformation matrix of SINS.On the other hand, the gyroscope error contains the constant error and random drift. Affected by inertial unit errors, the coordinate transformation matrixT is not accurate. When SINS works for a long time, the 6-degree-of-freedom navigation results are in oscillation because of incorrect T. According to the error analysis of SINS, under the undamped case the results include three kinds of different cycles of persistent oscillation [15]. They are, respectively, Schuler oscillation period Ts, earth oscillation period Te, and Foucault oscillation period TF.(4)Ts=2πRg≈84.8min,Te=2πΩ=24h,TF=2πΩsin⁡φ,where R is the earth radius; Ω is self-rotation angle velocity of the earth; φ is latitude; g is gravitational acceleration. Foucault oscillation period decreases with increased latitude. In addition, the period of earth oscillation and the period of Foucault oscillation are greater than the period of the heave motion, so this paper mainly considers the effects of Schuler oscillation for heave motion measurement.


### 2.1.3. Highpass Filter

Schuler oscillation whose period is 84.8 min in the process of SINS is caused by the constant bias of inertial sensors with high accuracy. The existence of Schuler oscillation makes heave information divergent. The frequency of Schuler oscillation is less than the frequency of heave motion. The highpass filter is designed to eliminate Schuler oscillation and reduce noise [16].In this paper a fourth-order highpass filter is given by the following equation:(5)Hs=η∬Adtdt=s4s2+2ξwcs+wc22,where ξ=1/2 is the damping factor, s is the Laplace variable, and wc is the cutoff frequency of the highpass filter which can be tuned.


## 2.1.1. Longuet-Higgins Wave Model

In order to describe the characteristics of the actual wave, Longuet-Higgins wave model is established on the basis of the measurement data [13]. According to the assumption of the model, the wave motion is a random process. The fixed point fluctuation of the wave is composed of some different phases and different amplitude sine wave.The wave surface displacement is expressed as(1)ηt=∑n=1∞ancos⁡wnt+εn.In the equation above, an and wn are the amplitude and frequency of the nth wave component, respectively, t represents time, and εn is the random phase uniformly distributed in the (0,2π). Furthermore, the sea wave only has a dominant component under most condition.


## 2.1.2. Inertial Unit Error Model

The accelerometer error model is given as follows [14]:(2)Ea=n+k0+kiGi+kiiGi2+kiiiGi3+kioGiGo+kipGiGp,where Ea is the accelerometer error, n is the white Gaussian noise, k0 is a constant bias, ki is linear error coefficient, kii, kiii are, respectively, second-order and third-order nonlinear error coefficient, and kio, kip are cross-axis coupling coefficient, respectively. Gi, G0, and Gp are, respectively, the specific force of input axis, output axis, and pendulum axis.The relationship of two models is given as follows:(3)η¨=T·A-Ea,where A is the measured value of accelerometers; T is coordinate transformation matrix of SINS.On the other hand, the gyroscope error contains the constant error and random drift. Affected by inertial unit errors, the coordinate transformation matrixT is not accurate. When SINS works for a long time, the 6-degree-of-freedom navigation results are in oscillation because of incorrect T. According to the error analysis of SINS, under the undamped case the results include three kinds of different cycles of persistent oscillation [15]. They are, respectively, Schuler oscillation period Ts, earth oscillation period Te, and Foucault oscillation period TF.(4)Ts=2πRg≈84.8min,Te=2πΩ=24h,TF=2πΩsin⁡φ,where R is the earth radius; Ω is self-rotation angle velocity of the earth; φ is latitude; g is gravitational acceleration. Foucault oscillation period decreases with increased latitude. In addition, the period of earth oscillation and the period of Foucault oscillation are greater than the period of the heave motion, so this paper mainly considers the effects of Schuler oscillation for heave motion measurement.


## 2.1.3. Highpass Filter

Schuler oscillation whose period is 84.8 min in the process of SINS is caused by the constant bias of inertial sensors with high accuracy. The existence of Schuler oscillation makes heave information divergent. The frequency of Schuler oscillation is less than the frequency of heave motion. The highpass filter is designed to eliminate Schuler oscillation and reduce noise [16].In this paper a fourth-order highpass filter is given by the following equation:(5)Hs=η∬Adtdt=s4s2+2ξwcs+wc22,where ξ=1/2 is the damping factor, s is the Laplace variable, and wc is the cutoff frequency of the highpass filter which can be tuned.


## 2.2. Analysis of Narrowband Noise

Due to the influence of random phase of the wave, the high order nonlinear error part sensed by the accelerometer can be equivalent to a narrowband noise after two integrations [17]. The white Gaussian noise is also in the form of narrowband noise in the heave information because of the output with a range of sampling frequency. The random process of narrowband noise in the wave model is discussed as follows.Without considering the impact of other errors which are the constant bias and linear errors of the accelerometer, narrowband noise in the random phase wave model is analyzed. The model can be expressed as follows:(6)Yt=∫∫Adtdt≈ηt+Nt,where ηt is the sine wave with the random phase(7)ηt=acos⁡w0t+β=acos⁡βcos⁡w0t-asin⁡βsin⁡w0t,where a is the known amplitude, β is the random phase uniformly distributed in the (0,2π), and w0 is the known angular frequency.In (5), narrowband noise Nt is a narrowband stationary Gaussian random process. The mean is zero and the variance is σ2. It is given as follows:(8)Nt=Qtcos⁡w0t+φnt=Qctcos⁡w0t-Qstsin⁡w0t.According to (6) and (7), (5) is rewritten as(9)Yt=Bctcos⁡w0t-Bstsin⁡w0t=Btcos⁡w0t+Φt.In the equation above,(10)Bct=acos⁡β+Qct,Bst=asin⁡β+Qst,Bt=Bc2t+Bs2t1/2,Φt=arctan⁡QstQct.Under the condition of a given signal phase θ and the time of t, with the value range of B and θ being B≥0, θ≤2π, the probability density function of the envelope is given by(11)pB∣θ=Bσ2exp⁡-a2+B22σ2I0aBσ2B≥0,where I0aB/σ2 is the zero-order Bessel function and aB/σ2 is the number of Bessel functions.Assuming thatb=a/σ, v=B/σ, (10) can be expressed as(12)pv=vσexp⁡-12v2+b2I0bvv≥0.In the actual situation, the useful signal of the accelerometer is much larger than the noise; therefore, bv≥1(13)pv≈vσexp⁡-12b2+v21vb2πexp⁡bv≈1σ2πexp⁡-12v-b2.With B≈a and v=b, the mean and variance of Bt are given by(14)EBt=a,DBt=σ2.However, some of frequency of the noise is the same as frequency of heave motion. So the important thing is to select the appropriate cutoff frequency of the filter to minimize errors.


## 3. Design and Error Analysis of Adaptive Algorithm

According to Section2, it is important to design the filter to get the heave information. While there are the phase errors and the amplitude attenuation in the filter, it is necessary to compensate them.


### 3.1. Adaptive Phase Compensation of Filter

The peak to peak error which is induced by the filter is presented as follows:(15)1-Hs=1-w4w4+wc4=wc4w4+wc4⟶wcw4forw≫wc.From (14), the peak to peak error cannot be ignored. Phase compensation is realized to reduce the peak to peak error by using first-order all pass system in this paper [18, 19]. It is given by(16)Hpz=z-1-k1-kz-10<k<1.The amplitude response of the system is not affected by the frequency so that all frequency components of the input signal can be passed without amplitude attenuation. The phase error caused by the filter can be compensated by choosing the appropriatek value.Then thez plane is mapped to the s plane through the expression z=esT, so(17)z=2/T+s2/T-s.With (15) and (16), we can get(18)Hps=s-2-2k/T1+k-s-2-2k/T1+k.In (17), T is the sampling period.By multiplying the highpass filterHfjw and the all pass system Hpjw, the expression is given as follows:(19)Hjw=Hfjw·Hpjw.w is the frequency of the heave motion. By making the phase zero, we can determine the parameter k value. But the heave motion is composed of several different phase waves. So with comprehensive consideration, an adaptive algorithm is proposed to select the parameter to minimize the phase error.(20)K=1∑n=1Nan∑n=1Nankn,where K is the comprehensive parameter; an is the amplitude of the nth wave component of the sea, and kn is the phase compensation parameter of the nth wave.


### 3.2. Adaptive Selection of The Cutoff Frequency

The amplitude attenuation of the filter is given by(21)1--Hs=wc8+8w6wc2w4+wc4⟶22wcwforw≫wc.The magnitude of errors is closely related to the cutoff frequency and the input frequency of the filter. Furthermore, according to analysis of narrowband noise, the variance ofBt is σ2, so the narrowband noise error is given as follows:(22)Eη¯n2t=12π∫-∞∞Hsjw2σ2dw=2-7/2σ2wc3.The magnitude of the errors and the noise are given as follows:(23)Eη¯2t=fwc,a,w,σ≤2-7/2σ2wc3+8a2wcw2.According to the change of the frequency w, selecting the optimal cutoff frequency wc minimizes the magnitude of errors. The method is given by (24)Ewc=Eη¯2t,∂Ewc∂wc=0,(25)wc=0.0166σ2w2a21/5.Only with one frequency of the signal, a single optimal cutoff frequency can be obtained by (24). However, the heave motion is composed of some sine waves. Therefore, on this basis, an adaptive algorithm is proposed to find the comprehensive optimal cutoff frequency. It is given by(26)Wc=∑n=1Nanwcn∑n=1Nan,where Wc is the comprehensive optimal cutoff frequency and wcn represents the single optimal cutoff frequency of the nth wave component of the sea wave.


### 3.3. Error Analysis

The method can estimate heave motion, but there still exist errors at heave information. It is necessary to analyze range of errors. So(27)EP2t≤∑n=1N2-7/2σ2Wc3+8a2Wcwn2.In the above equation, wn is the frequency of the nth wave component. They are determined by the sea wave. Even though the sea wave motion is changing, the range of error can be obtained by (26).


### 3.4. The Steps of Heave Measure

The steps of heave measure are shown in Figure1. Firstly, the displacement of altitude channel is calculated by SINS. Secondly, the frequency and the amplitude of each wave component are obtained by FFT spectrum analysis. Next, the comprehensive parameter K and the comprehensive optimal cutoff frequency Wc are, respectively, calculated by the adaptive algorithm. Finally, the heave motion information can be obtained by the adaptive filter and phase compensation.Figure 1
The steps of heave measurement.


## 3.1. Adaptive Phase Compensation of Filter

The peak to peak error which is induced by the filter is presented as follows:(15)1-Hs=1-w4w4+wc4=wc4w4+wc4⟶wcw4forw≫wc.From (14), the peak to peak error cannot be ignored. Phase compensation is realized to reduce the peak to peak error by using first-order all pass system in this paper [18, 19]. It is given by(16)Hpz=z-1-k1-kz-10<k<1.The amplitude response of the system is not affected by the frequency so that all frequency components of the input signal can be passed without amplitude attenuation. The phase error caused by the filter can be compensated by choosing the appropriatek value.Then thez plane is mapped to the s plane through the expression z=esT, so(17)z=2/T+s2/T-s.With (15) and (16), we can get(18)Hps=s-2-2k/T1+k-s-2-2k/T1+k.In (17), T is the sampling period.By multiplying the highpass filterHfjw and the all pass system Hpjw, the expression is given as follows:(19)Hjw=Hfjw·Hpjw.w is the frequency of the heave motion. By making the phase zero, we can determine the parameter k value. But the heave motion is composed of several different phase waves. So with comprehensive consideration, an adaptive algorithm is proposed to select the parameter to minimize the phase error.(20)K=1∑n=1Nan∑n=1Nankn,where K is the comprehensive parameter; an is the amplitude of the nth wave component of the sea, and kn is the phase compensation parameter of the nth wave.


## 3.2. Adaptive Selection of The Cutoff Frequency

The amplitude attenuation of the filter is given by(21)1--Hs=wc8+8w6wc2w4+wc4⟶22wcwforw≫wc.The magnitude of errors is closely related to the cutoff frequency and the input frequency of the filter. Furthermore, according to analysis of narrowband noise, the variance ofBt is σ2, so the narrowband noise error is given as follows:(22)Eη¯n2t=12π∫-∞∞Hsjw2σ2dw=2-7/2σ2wc3.The magnitude of the errors and the noise are given as follows:(23)Eη¯2t=fwc,a,w,σ≤2-7/2σ2wc3+8a2wcw2.According to the change of the frequency w, selecting the optimal cutoff frequency wc minimizes the magnitude of errors. The method is given by (24)Ewc=Eη¯2t,∂Ewc∂wc=0,(25)wc=0.0166σ2w2a21/5.Only with one frequency of the signal, a single optimal cutoff frequency can be obtained by (24). However, the heave motion is composed of some sine waves. Therefore, on this basis, an adaptive algorithm is proposed to find the comprehensive optimal cutoff frequency. It is given by(26)Wc=∑n=1Nanwcn∑n=1Nan,where Wc is the comprehensive optimal cutoff frequency and wcn represents the single optimal cutoff frequency of the nth wave component of the sea wave.


## 3.3. Error Analysis

The method can estimate heave motion, but there still exist errors at heave information. It is necessary to analyze range of errors. So(27)EP2t≤∑n=1N2-7/2σ2Wc3+8a2Wcwn2.In the above equation, wn is the frequency of the nth wave component. They are determined by the sea wave. Even though the sea wave motion is changing, the range of error can be obtained by (26).


## 3.4. The Steps of Heave Measure

The steps of heave measure are shown in Figure1. Firstly, the displacement of altitude channel is calculated by SINS. Secondly, the frequency and the amplitude of each wave component are obtained by FFT spectrum analysis. Next, the comprehensive parameter K and the comprehensive optimal cutoff frequency Wc are, respectively, calculated by the adaptive algorithm. Finally, the heave motion information can be obtained by the adaptive filter and phase compensation.Figure 1
The steps of heave measurement.


## 4. Simulation and Measurement Results

Experiments based on the Longuet-Higgins wave model are performed to validate the effectiveness of the adaptive algorithm. To simulate different sea conditions, the model is built by the computer simulation. It is assumed that the latitude is 30° and the longitude is 120°. The constant error and random drift of the gyroscope are, respectively, 0.01°/h and 0.005°/(h·√Hz). The constant error of the accelerometer is 100 ug and the random drift is 80 ug/(h·√Hz).


### 4.1. Simulation Test under Certain Sea Condition

To get closer to the actual sea vessel heave motion, the simulation of true heave motion based on the Longuet-Higgins wave model is implemented. What is more, in order to confirm the adaptive algorithm’s effectiveness, the single optimal cutoff frequencywc and the comprehensive optimal cutoff frequency Wc are, respectively, selected as the cutoff frequency of the filter to estimate heave motion information.The Sea Model Setup. The main frequency of true heave motion is 0.1 Hz and the amplitude is 3 m; others are 0.14 Hz, 0.1 m; 0.08 Hz, 0.2 m; 0.15 Hz, 0.02 m; 0.07 Hz, 0.01 m; 0.15 Hz, 0.008 m; 0.06 Hz, 0.009 m; 0.011 Hz, 0.005 m; 0.09 Hz, 0.002 m, and so on. According to Longuet-Higgins wave model, the wave is composed of much trigonometric function, but the amplitude or energy of much wave component is low, so they can be ignored. The simulation of the sea wave is shown in Figure 2.Figure 2
Simulation of true heave motion.In Figure3, the cutoff frequency of the filter is the single optimal cutoff frequency wc. In Figure 4, the cutoff frequency of the filter is the comprehensive optimal cutoff frequency Wc.Figure 3
Measurement with the single optimal cutoff frequency.Figure 4
Measurement with the comprehensive optimal cutoff frequency.According to (24), the single optimal cutoff frequency of the dominate wave component of the heave movement is obtained wc=0.0168 rad/s. As shown in Figure 3, the black line represents filtered but uncompensated heave, the blue line represents filtered and compensated heave, and the red line represents the true value of the heave simulation. After filtering, the maximum error is 0.4942 m, the average value is −0.0372 m, and the standard deviation is 0.1874 m. The error ratio that the maximum error is divided by the amplitude of the heave is 16.43%. There is a big error by comparing the black with the red and only a little amplitude error by comparing the blue with the red.As shown in Figure4, the comprehensive optimal cutoff frequency wc=0.0198 rad/s is obtained from (25). The maximum error is 0.1170 m, the average value is 8.929e-05 m, and the standard deviation is 0.0453 m. The error ratio is 3.89%. We can see that the blue line almost coincides with the red line.By comparing two simulation experiment results, the maximum error has been reduced by 0.3772 m. The average value has been almost 0 m. The standard deviation has been reduced to 0.0453 m. The error ratio has been reduced by 12.54%. It is clear that the results of the comprehensive optimal cutoff frequency are superior to the results of the single optimal cutoff frequency in this case.According to wave model, the phase of the wave component is a random process. Four-group simulation experiments with different phases are conducted to validate the proposed algorithm. The maximum error, the average value, the standard deviation, and the error ratio are shown in Table1 according to the single and the comprehensive optimal cutoff frequency. The experiments show that the adaptive algorithm is also suitable for the different phases in the same case.Table 1
Four-group simulation experiments.
Simulation number Cutoff frequency Maximum value of absolute error (m) Average value of errors (m) Standard deviation of errors (m) Error ratio (%) 1 w c 0.4925 - 0.0370 0.1869 16.37 W c 0.1226 - 2.877 e - 05 0.0453 4.076 2 w c 0.5006 - 0.0367 0.1868 16.64 W c 0.1206 6.709 e - 05 0.0453 4.009 3 w c 0.4929 - 0.0367 0.1872 16.39 W c 0.1181 3.548 e - 05 0.0453 3.926 4 w c 0.5000 - 0.0369 0.1871 16.62 W c 0.1263 6.872 e - 05 0.0453 4.199


### 4.2. Test Results under Different Sea Cases

Considering different sea states, the waves are divided into three types which are, respectively, rough sea, moderate sea, and slight sea, respectively. The amplitude range of rough sea is 2 m to 3 m. The amplitude range of moderate sea is 1 m to 2 m. The amplitude range of slight sea is 0 m to 1 m. Therefore, three tests are conducted to validate the validity of the adaptive algorithm under different sea height circumstances.Test 1. The model of rough sea is built that the main frequency of true heave motion is 0.1 Hz and the amplitude is 2.5 m; others are 0.08 Hz, 0.15 m, and 0.14 Hz, 0.1 m.Test 2. The model of moderate sea is built that the main frequency of true heave motion is 0.1 Hz and the amplitude is 1.6 m; others are 0.08 Hz, 0.05 m, and 0.14 Hz, 0.05 m.Test 3. The model of slight sea is built that the main frequency of true heave motion is 0.1 Hz and the amplitude is 0.6 m; others are 0.08 Hz, 0.02 m, and 0.14 Hz, 0.01 m.Tables2, 3, and 4 are, respectively, the measurement results of Test 1, Test 2, and Test 3.Table 2
Rough sea test.
Simulation number Cutoff frequency Maximum value of absolute error (m) Average value of errors (m) Standard deviation of errors (m) Error ratio (%) 1 w c 0.2945 0.0182 0.0971 11.74 W c 0.1181 1.326 e - 05 0.0422 4.712 2 w c 0.2969 0.0184 0.0975 11.85 W c 0.1257 3.084 e - 05 0.0422 5.015 3 w c 0.2865 0.0182 0.0974 11.43 W c 0.1193 2.117 e - 05 0.0422 4.760 4 w c 0.2948 0.0185 0.0974 11.76 W c 0.1236 2.117 e - 05 0.0422 4.931Table 3
Moderate sea test.
Simulation number Cutoff frequency Maximum value of absolute error (m) Average value of errors (m) Standard deviation of errors (m) Error ratio (%) 1 w c 0.2601 0.0183 0.0924 16.24 W c 0.0852 - 3.007 e - 05 0.0235 5.320 2 w c 0.2582 0.0186 0.0928 16.12 W c 0.0817 - 1.140 e - 5 0.0235 5.101 3 w c 0.2570 0.0186 0.0926 16.05 W c 0.0818 3.896 e - 05 0.0235 5.107 4 w c 0.2574 0.0184 0.0921 16.07 W c 0.0881 - 3.474 e - 05 0.0235 5.500Table 4
Slight sea test.
Simulation number Cutoff frequency Maximum value of absolute error (m) Average value of errors (m) Standard deviation of errors (m) Error ratio (%) 1 w c 0.2387 0.0185 0.0915 39.76 W c 0.0600 3.552 e - 04 0.0134 9.993 2 w c 0.2354 0.0182 0.0909 39.21 W c 0.0557 3.479 e - 04 0.0134 9.277 3 w c 0.2355 0.0185 0.0910 39.22 W c 0.0548 3.418 e - 04 0.0134 9.127 4 w c 0.2370 0.0184 0.0914 39.47 W c 0.0583 3.452 e - 04 0.0134 9.710As Table2 shows, the average error ratio of the rough sea test reduces to 4.855% from 11.70%. As Table 3 shows, the average error ratio of the moderate sea test reduces to 5.255% from 16.12%. As Table 4 shows, the average error ratio of the slight sea test reduces to 9.527% from 39.42%. Furthermore, the maximum error, the average value, and the standard deviation obviously have decreased by the adaptive algorithm in all tests.The measure results show that the adaptive algorithm is effective in getting better accuracy estimation of heave motion under different sea cases. In addition, the error radio of the rough sea is the smallest, then followed by the error radio of the moderate sea; the worst is the error radio of the slight sea. Thus, the adaptive algorithm is more suitable for the rough sea situation.


## 4.1. Simulation Test under Certain Sea Condition

To get closer to the actual sea vessel heave motion, the simulation of true heave motion based on the Longuet-Higgins wave model is implemented. What is more, in order to confirm the adaptive algorithm’s effectiveness, the single optimal cutoff frequencywc and the comprehensive optimal cutoff frequency Wc are, respectively, selected as the cutoff frequency of the filter to estimate heave motion information.The Sea Model Setup. The main frequency of true heave motion is 0.1 Hz and the amplitude is 3 m; others are 0.14 Hz, 0.1 m; 0.08 Hz, 0.2 m; 0.15 Hz, 0.02 m; 0.07 Hz, 0.01 m; 0.15 Hz, 0.008 m; 0.06 Hz, 0.009 m; 0.011 Hz, 0.005 m; 0.09 Hz, 0.002 m, and so on. According to Longuet-Higgins wave model, the wave is composed of much trigonometric function, but the amplitude or energy of much wave component is low, so they can be ignored. The simulation of the sea wave is shown in Figure 2.Figure 2
Simulation of true heave motion.In Figure3, the cutoff frequency of the filter is the single optimal cutoff frequency wc. In Figure 4, the cutoff frequency of the filter is the comprehensive optimal cutoff frequency Wc.Figure 3
Measurement with the single optimal cutoff frequency.Figure 4
Measurement with the comprehensive optimal cutoff frequency.According to (24), the single optimal cutoff frequency of the dominate wave component of the heave movement is obtained wc=0.0168 rad/s. As shown in Figure 3, the black line represents filtered but uncompensated heave, the blue line represents filtered and compensated heave, and the red line represents the true value of the heave simulation. After filtering, the maximum error is 0.4942 m, the average value is −0.0372 m, and the standard deviation is 0.1874 m. The error ratio that the maximum error is divided by the amplitude of the heave is 16.43%. There is a big error by comparing the black with the red and only a little amplitude error by comparing the blue with the red.As shown in Figure4, the comprehensive optimal cutoff frequency wc=0.0198 rad/s is obtained from (25). The maximum error is 0.1170 m, the average value is 8.929e-05 m, and the standard deviation is 0.0453 m. The error ratio is 3.89%. We can see that the blue line almost coincides with the red line.By comparing two simulation experiment results, the maximum error has been reduced by 0.3772 m. The average value has been almost 0 m. The standard deviation has been reduced to 0.0453 m. The error ratio has been reduced by 12.54%. It is clear that the results of the comprehensive optimal cutoff frequency are superior to the results of the single optimal cutoff frequency in this case.According to wave model, the phase of the wave component is a random process. Four-group simulation experiments with different phases are conducted to validate the proposed algorithm. The maximum error, the average value, the standard deviation, and the error ratio are shown in Table1 according to the single and the comprehensive optimal cutoff frequency. The experiments show that the adaptive algorithm is also suitable for the different phases in the same case.Table 1
Four-group simulation experiments.
Simulation number Cutoff frequency Maximum value of absolute error (m) Average value of errors (m) Standard deviation of errors (m) Error ratio (%) 1 w c 0.4925 - 0.0370 0.1869 16.37 W c 0.1226 - 2.877 e - 05 0.0453 4.076 2 w c 0.5006 - 0.0367 0.1868 16.64 W c 0.1206 6.709 e - 05 0.0453 4.009 3 w c 0.4929 - 0.0367 0.1872 16.39 W c 0.1181 3.548 e - 05 0.0453 3.926 4 w c 0.5000 - 0.0369 0.1871 16.62 W c 0.1263 6.872 e - 05 0.0453 4.199


## 4.2. Test Results under Different Sea Cases

Considering different sea states, the waves are divided into three types which are, respectively, rough sea, moderate sea, and slight sea, respectively. The amplitude range of rough sea is 2 m to 3 m. The amplitude range of moderate sea is 1 m to 2 m. The amplitude range of slight sea is 0 m to 1 m. Therefore, three tests are conducted to validate the validity of the adaptive algorithm under different sea height circumstances.Test 1. The model of rough sea is built that the main frequency of true heave motion is 0.1 Hz and the amplitude is 2.5 m; others are 0.08 Hz, 0.15 m, and 0.14 Hz, 0.1 m.Test 2. The model of moderate sea is built that the main frequency of true heave motion is 0.1 Hz and the amplitude is 1.6 m; others are 0.08 Hz, 0.05 m, and 0.14 Hz, 0.05 m.Test 3. The model of slight sea is built that the main frequency of true heave motion is 0.1 Hz and the amplitude is 0.6 m; others are 0.08 Hz, 0.02 m, and 0.14 Hz, 0.01 m.Tables2, 3, and 4 are, respectively, the measurement results of Test 1, Test 2, and Test 3.Table 2
Rough sea test.
Simulation number Cutoff frequency Maximum value of absolute error (m) Average value of errors (m) Standard deviation of errors (m) Error ratio (%) 1 w c 0.2945 0.0182 0.0971 11.74 W c 0.1181 1.326 e - 05 0.0422 4.712 2 w c 0.2969 0.0184 0.0975 11.85 W c 0.1257 3.084 e - 05 0.0422 5.015 3 w c 0.2865 0.0182 0.0974 11.43 W c 0.1193 2.117 e - 05 0.0422 4.760 4 w c 0.2948 0.0185 0.0974 11.76 W c 0.1236 2.117 e - 05 0.0422 4.931Table 3
Moderate sea test.
Simulation number Cutoff frequency Maximum value of absolute error (m) Average value of errors (m) Standard deviation of errors (m) Error ratio (%) 1 w c 0.2601 0.0183 0.0924 16.24 W c 0.0852 - 3.007 e - 05 0.0235 5.320 2 w c 0.2582 0.0186 0.0928 16.12 W c 0.0817 - 1.140 e - 5 0.0235 5.101 3 w c 0.2570 0.0186 0.0926 16.05 W c 0.0818 3.896 e - 05 0.0235 5.107 4 w c 0.2574 0.0184 0.0921 16.07 W c 0.0881 - 3.474 e - 05 0.0235 5.500Table 4
Slight sea test.
Simulation number Cutoff frequency Maximum value of absolute error (m) Average value of errors (m) Standard deviation of errors (m) Error ratio (%) 1 w c 0.2387 0.0185 0.0915 39.76 W c 0.0600 3.552 e - 04 0.0134 9.993 2 w c 0.2354 0.0182 0.0909 39.21 W c 0.0557 3.479 e - 04 0.0134 9.277 3 w c 0.2355 0.0185 0.0910 39.22 W c 0.0548 3.418 e - 04 0.0134 9.127 4 w c 0.2370 0.0184 0.0914 39.47 W c 0.0583 3.452 e - 04 0.0134 9.710As Table2 shows, the average error ratio of the rough sea test reduces to 4.855% from 11.70%. As Table 3 shows, the average error ratio of the moderate sea test reduces to 5.255% from 16.12%. As Table 4 shows, the average error ratio of the slight sea test reduces to 9.527% from 39.42%. Furthermore, the maximum error, the average value, and the standard deviation obviously have decreased by the adaptive algorithm in all tests.The measure results show that the adaptive algorithm is effective in getting better accuracy estimation of heave motion under different sea cases. In addition, the error radio of the rough sea is the smallest, then followed by the error radio of the moderate sea; the worst is the error radio of the slight sea. Thus, the adaptive algorithm is more suitable for the rough sea situation.


## 5. Conclusions

This paper discusses the obtainment of the heave information by the adaptive filter. An adaptive algorithm based on the Longuet-Higgins wave model and random noise analysis is proposed to reduce the wave information errors. It adjusts the optimal cutoff frequency by analyzing the wave and is more suitable for the real ocean wave. The adaptive filter is useful to improve accuracy of heave measure under three sea states. The method in this paper is based on normal wave condition with only one primary wave. If the wave condition is very complicated especially with more than two primary waves, we need to design special filters.



---

*Source: 6547312-2017-02-28.xml*