Performance Analysis of Target Depth Classification Algorithm Based on Sea Experiment Data

,


Introduction
The movement of underwater platforms in channel encounters the security threats caused by various targets.To achieve the passive localization of the target without the target prior information is very essential.The target depth classification results directly related to the target category in the passive localization process are particularly important.Effective classification of targets (surface targets or underwater targets) has a profound impact on the security and concealment of underwater platforms.
The main research focus of this paper is to study a method for target depth classification, which is actually the method for target category classification.In recent years, this area of research has developed huge interests and a large number of research work are published in the field of target depth classification.
The initial purpose of the target category is to classify the target depth, rather than accurately estimating the target depth.There are numerous algorithms for depth classification using array signals.Golden [1] performed maximumlikelihood estimation for localization.Liang et al. [2] and Premus et al. [3] proposed mode filtering approaches to acoustic source depth discrimination.Based on either modal decomposition technique [4] or data-based method [5], source depth estimation was realized by Nicolas and Yang.Target depth classification was achieved using a waveguide invariant adaptive matched filter by Goldhahn [6].There are various researches about target depth estimation on the basis of matched field processing [7][8][9][10][11][12].Premus and Backman [13] proposed a matched subspace [14] approach for depth discrimination.The concept of modified modal scintillation index was introduced for source depth classification by An et al. [15].By using cepstrum techniques, the determination of source depth was carried out by Wu, Zhang et al., and Mitchell and Bedford [16][17][18].Source range and depth are obtained using a ray-based backpropagation algorithm by Felisberto et al. [19].These depth classification algorithms have high accuracy, but most of these methods are not feasible to be applied in real application due to their high complexity and time consumption.In addition, considering the space limitations and security requirements of underwater platforms, it may be impossible to place or tow an array with a large aperture.However, these research hold high theoretical value but need to be further researched for implying in real application scenario for target depth classification by placing one or two sensors.
The cross-spectrum calculation is performed using signals acquired by a single 3D vector sensor [20,21] or two 2D vector sensors [22,23], and the target depth classification can be achieved based on the characteristics of the crossspectrum signals.However, most of the above algorithms are theoretical studies [20][21][22][23], and the actual performance of the algorithm is rarely analyzed through the analysis of experiment data processing results [24,25].Based on the sea experiment data processing results, this paper analyzes the influence of the ideal receiving depth concept on the performance of the algorithm while verifying the feasibility of the target depth classification algorithm.It shows the real application value of the ideal receiving depth concept.Moreover, the paper makes an in-depth research about the affection of the sound velocity profile parameters (negative gradient, thermocline intensity, thermocline thickness, upboundary depth, etc.) on the algorithm performance.

Theoretical Model
The acoustic propagation characteristics of shallow water field on the basis of the two-layer fluid-fluid model is shown in Figure 1.H is the sea depth; seawater density ρ 1 is constant, and the sound velocity of seawater c 1 is constant or variable with depth z only.The sediment layer is homogeneous fluid half-space whose sound velocity c 2 and density ρ 2 are constant.A point source of unit amplitude and angular frequency ω is located at (0, z 0 ) in cylindrical coordinates.The receiver is located at (r, z).The sea surface is an absolutely soft interface whose depth is z = 0 [26].
If c 1 is constant, this model is known as the Pekeris model.A Pekeris channel of depth H in the shallow water is essentially the same as a waveguide with two pressurerelease boundaries whose depth is H e [27].H e is called the effective depth.
is the critical grazing angle of the bottom interface.Thus, the pressure field can be expressed as follows: where n is the serial number of the normal mode; is the nth order eigenvalue; ξ n and β n are also called horizontal and vertical wave number, respectively.Ψ n ðzÞ = sin ðβ n zÞ is the mode depth function.In the effective depth model, Pðz 0 , z 1 , rÞ and Pðz 0 , z 2 , rÞ are signals collected by receiving pressure sensors at different depths z 1 and z 2 .The active component of their cross-spectrum signal is defined as I pA : In the equation, A in ðz 0 , 2, the superscript * represents the complex conjugate operation; Re ½⋅ represents the operation of getting real part from the complex number.
The cutoff frequency of the nth order normal mode can be calculated as S rz represents the component of I pA whose sign changes in the horizontal and vertical direction.S z is denoted for the component of I pA whose sign changes in the vertical direction only.According to equation ( 5), the peak value of S rz is ∑ n,n≠m ∑ m A 1n A 2m .S rz resembles cosine signal but its peak value varies with source and receiver depth.
The frequency can only excite the first two normal modes, and target depth classification can be realized using sign distribution of I pA .Based on effective depth model, when source frequency can only excite the first two modes, there are ideal receiving depths z 1 + z 2 = H e resulting in ðA 11 A 22 + A 12 A 21 Þ cos ðΔξ 12 rÞ ≈ 0, so that S rz is tending to zero.The impact of S rz on the sign distribution of I pA can be ignored.If z 1 + z 2 = H e , it can guarantee that the sign distribution of I pA will not change in the horizontal direction but change in the vertical direction only.
The critical depth h s is defined as the position of boundary of I pA sign when z 1 + z 2 = H e .All the signs of I pA above h s are identical, but different from the value lesser than h s .The approximate expression of h s is The flow chart of target depth classification is presented in Figure 2. According to the simulation parameters used in this paper, when the source depth is lesser than h s , the sign of I pA is negative, and sources located at these positions can be identified as surface targets.When the source depth is greater than h s , the sign of I pA is positive, and the source at these positions can be identified as underwater targets.Corresponding to different simulation parameters, the sign of surface targets' I pA can be positive and the sign of underwater targets' I pA can be negative.
Based on the above theory which have already been published in [22,23], we can conclude that the value of critical depth changes with different receiving depth, source frequency, sound velocity profile, and seabed.The reason why these factors affect algorithm performance has never been researched.We have to develop a model of these parameters affecting the algorithm effectiveness.The selection method of receiving depths has also been described briefly.
Only if ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ffi ξ 2 /16ξ 1 cos 2 ðγz 1 Þ p ≤ 1, the existence of critical depth h s can be guaranteed.Thus, the maximum allowable value of the first sensor depth is b⋅c represents that fractions are rounded down.
If z 1 + z 2 ≠ H e , the sign of S rz changes in the horizontal direction whose variation law depends on cos ðΔξ 12 rÞ component.The sign of cos ðΔξ 12 rÞ will hop every other π/Δξ 12 in the horizontal direction.The sign of S rz changes in the vertical direction whose variation law depends on ðA 11 A 22 + A 12 A 21 Þ component.The sign of ðA 11 A 22 + A 12 A 21 Þ will hop at about H e /2.It can be found that whether in the horizontal or vertical direction, the sign hopping of S rz is both not connected with the value of z 1 , z 2 , but the amplitude of S rz varies with z 1 , z 2 .So, the sign of I pA will vary with z 1 , z 2 .At most depths, S rz ≠ 0, so that S rz will exert significant influence on the sign distribution of I pA .On the premise of the existence of the critical depth, the existence of S rz will make the boundary of I pA sign change from straight line to curve that resembles cosine distribution whose peak value changes with z 1 , z 2 .
Summarizing our presented algorithm, the acoustic field interference structure I pA changes with the variation of receiving depths, source frequency, sound velocity profile, and seabed, which alters the value of h s and the amplitude of S rz accordingly.The change of h s and S rz leads to fluctuate the algorithm performance.The ideal receiving depths corresponding to various conditions are different.So the performance can be improved through proper adjustment of the receiving depths according to the change of these key factors.

Simulation Results
To verify the feasibility of our presented technique, we conducted a simulation experiment.Sea depth is selected as H = 22 m; densities of seawater and sediment are ρ 1 = 1:0g/cm 3 , ρ 2 = 1:75g/cm 3 ; sound velocity in the sediment is c 2 = 1750m/s; source frequency is f = 120 Hz, and the range of source depth is 1~22 m.The range of horizontal distance is 2~8 km.f 2 = 100Hz, f 3 = 165Hz; thus, the source frequency can only excite the first two modes.The impact of sound velocity profile on the algorithm performance has been analyzed considering three common types of sound velocity profile (isovelocity, negative gradient, and thermocline) in shallow water.Corresponding to isovelocity profile, sound velocity in the seawater is c 1 = 1500m/s.

Journal of Sensors
With respect to different negative gradient profiles, sound velocity in the seawater is c 1 = 1540m/s at z = 0m and c 1 = ½1520, 1500, 1480m/s at z = 22m; the gradient is negative and constant, namely, the gradient of different conditions is G c = ½0:9, 1:8, 2:7s −1 .
Thermocline profiles are effected mainly due to thermocline intensity, thermocline thickness, and up-boundary depth.The sound velocity in the seawater is c 1 = 1540m/s at z = 0m, and c 1 = ½1530, 1520, 1510, 1500m/s at z = 22m corresponding to different conditions.Thus, the thermocline intensity can be G c = ½2:5, 5s −1 , the thermocline thickness may be considered as 4m, 8m, 12m, and the up-boundary depth may be 5m, 10m.

Impact of Receiving Depth on Algorithm Performance.
Under the condition of isovelocity profile, according to equation (1), H e ≈ 29m.According to equation ( 8), only if z 1 ≤ 12m, target depth classification can be realized.The relationship between S rz , h s and z 1 , z 2 is shown in Figure 3.It can be witnessed from Figure 3 that with the increase of z 1 , the peak value of S rz gradually increases.The peak value of S rz will also vary with the distance 4 Journal of Sensors between z 1 + z 2 and H e , so as to affect the sign distribution of I pA ; h s will reduce with the increase of z 1 .Consequently, in order to get a better classification of target depth, the value of z 1 , z 2 should be considered comprehensively.Firstly, the ideal range of z 1 , z 2 should be obtained.Then, conferring to the requirement of h s , the ideal range of z 1 should be confirmed.Finally, in conditions permitting, the value of z 1 should be reduced appropriately to improve the performance of the depth classification algorithm.Figure 4 illustrates the impact of z 1 on the sign distribution of S z and I pA .In this article, all sign distributions are indicated by white and grey as positive sign and black indicates negative sign.It can be observed from Figure 4 that the critical depth h s is inversely proportional to z 1 .When z 1 ≤ 12m, the target depth classification can be achieved by using the sign distribution of I pA .If z 1 > 12m, the sign of I pA does not change in the horizontal direction, but the signs are all positive.The sign boundary does not exist and the target depth classification cannot be performed.
Figure 5 expresses the impact of z 1 + z 2 on the sign distribution of S z and I pA when z 1 = 11m.As it can be witnessed from Figure 5, when z 1 + z 2 = ½26, 27m, the sign of I pA changes with horizontal distance r, which coincides with the change rule shown in Figure 3.With the change of the distance between z 1 + z 2 and H e , S rz does not tend to 0. S rz affects the sign distribution of I pA far more than S z resulting in the sign distributions as shown in Figure 5 when z 1 + z 2 = ½26, 27m.If z 1 + z 2 = ½28, 29m, S rz tends to 0, the sign of I pA hardly changes with r. Figure 6 is the sign distribution of S rz under the condition of isovelocity profile.The position of the sign hopping in the vertical direction is about 14 m, roughly equal to H e /2.

Impact of Source Frequency on Algorithm Performance.
When the source frequency can only excite the first two modes, the change of the source frequency will affect the target depth classification effect.It can be known from equation (1) that, in case of isovelocity profile, the effective depth H e is inversely proportional to the frequency.The maximum depth of the first receiving sensor (z 1 max ) is related to H e , ξ 1 , ξ 2 , and the change of the frequency causes the change of the above three variables.Therefore, there is no obvious changing law of z 1 max along with frequency.Similarly, the critical depth h s is also related to H e , ξ 1 , ξ 2 , so the relationship between h s and frequency does not follow any obvious regular pattern.When f ∈ ½100, 165, H e ∈ ð26:9, 30:2Þ, the change scope of H e is small.h s and z 1 max are both directly related to ξ 2 /ξ 1 .When f ∈ ½100, 165, ξ 2 /ξ 1 ∈ ð0:8875, 0:9520Þ, the change scope of ξ 2 /ξ 1 is also not large.Therefore, as the frequency changes, h s and z 1 max will vary, but the impact of frequency on the It can be witnessed from Figure 7 that the change of the frequency affects the value of h s and z 1 max and the sign distribution of I pA changes too.However, by adjusting the receiving depth appropriately, based on the sign distribution of I pA , the depth classification algorithm can still obtain ideal classification results under different frequency conditions.Therefore, the change in frequency has minor impact on the performance of the algorithm.

Impact of Signal-Noise Ratio on Algorithm Performance.
In this section, we will discuss the effect of SNR on the performance of our presented algorithm.The research object changes from point source to radiation noise of surface or underwater targets.The parameters considered are as follows: heading angle is ψ = 50 ∘ , target velocity is v t = 11m/s, closest distance is r min = 1200m, initial distance is r 0 = 2000 m, platform velocity is v s = 2m/s, and range of SNR (signalnoise ratio) is SNR = −30~0dB.Assuming that the target radiated noise is a single-frequency line spectrum superimposed on the continuum, two pressure sensors are used to receive signals at depths of 45 m and 90 m, respectively, and the total sailing time is 400 s.The target is kept in motion away from the receiving platform.With respect to the abovementioned simulation parameters, the influence of SNR on the performance of the algorithm is analyzed.Figure 8 expresses the relationship between the performance of the target depth classification algorithm and SNR under isovelocity condition.In the process of target depth classification, the index that implies the performance of the algorithm is the probability of identified as underwater target in percent-age.The probability of identified as underwater target is defined as the proportion of the cases in which the sign of the pressure cross-spectrum active component is positive.It is equal to is the total number of times of Monte Carlo simulation; M 1 is the total sampling numbers in each Monte Carlo simulation; L i (i = 1, 2, ⋯, N 1 ) is the number of times of the cases in which the sign is positive when discriminating the sign of the pressure cross-spectrum active component in each simulation, 0 ≤ L i ≤ M 1 .It can be found from Figure 8 that when SNR≥−20dB, the target depth classification effect is good, and there is an obvious sign boundary at about 7 m.When SNR = −30dB, the algorithm can achieve the target depth classification, but the classification effect is poor.It can be witnessed from Figure 9 that the sound velocity gradient is inversely proportional to the depth classification effect and is directly proportional to the critical depth.Even in the case of strong negative gradient, the robustness of the algorithm is high.

Thermocline Existing in Sound Velocity
Profile.When the thermocline thickness is fixed to 8m and the upboundary depth is fixed to 5m, the sound velocity profiles and the sign distributions of I pA under different conditions of G c = ½2:5, 5s −1 are displayed in Figure 10.When the thermocline intensity is fixed to G c = 2:5s −1 and the thermocline thickness is fixed to 8m, the sound velocity profiles and the sign distributions of I pA when the up-boundary depths are 5m, 10m are shown in Figure 10.
It can be witnessed from Figure 10 that the algorithm in this paper is also applicable when thermocline exists in the sound velocity profile.The depth classification effect of the algorithm deteriorates as thermocline intensity increases, and the critical depth value is nearly constant.While the depth classification effect of the algorithm deteriorates as the depth of the up-boundary increases, the critical depth value almost keeps the same.In the presence of thermocline, both the thermocline intensity and up-boundary depth have a large impact on the depth classification performance of the algorithm, but the critical depth value remains almost the same.Even with the presence of a strong thermocline, the depth classification performance of the algorithm is good.
When the thermocline intensity is fixed to G c = 2:5s −1 and the up-boundary depth is fixed to 5m, the sound velocity profiles and the sign distributions of I pA when the thermocline thicknesses are 4m, 8m, 12m are displayed in Figure 11.
It can be seen from Figure 11 that the depth classification effect of the algorithm deteriorates with the increase of the thermocline thickness and the critical depth value increases with the increase of the thermocline thickness.Therefore, the thermocline thickness has a great impact on the depth classification performance of the algorithm.

Experiment Results
In order to examine the influence of the ideal receiving depth concept on the performance of target depth classification algorithm, an experiment was carried out in Yellow Sea in Qingdao, China.The layout of the sea trial is displayed in Figure 12.The source is suspended into the water by the sending ship at about 8 m.The sound velocity profile near the source position is calculated by CTD (Conductivity-Temperature-Depth System).The CTD device is shown in Figure 13(a), and the collected sound velocity profile is graphed in Figure 14.It can be witnessed from Figure 14 that Four receiving pressure sensors are equally spaced into the water by the receiving ship, and one TD (depth sensor) is fixed at each sensor.The signal collection device is fixed between the second and third sensors, and all sensors are connected to the signal collection device.The effect of the ocean current on the experiment results is reduced by fixing the weight on the rope.The receiving device is shown in Figure 13(b).
According to the CTD, the sea depth is approximately 22 m.It can be obtained from equation ( 1) that H e ∈ ð29, 30 Þm.According to equation ( 8), z 1 max < 13m.The length of the pressure sensor cable is limited; in order to ensure the stability of the sensor power supply, the overall length of pressure array length cannot exceed 6 m.The four array elements are laid at equal intervals to keep the interval between two neighboring elements of 2 m.The designed depths of the four array elements are 11, 13, 15, and 17 m.According to the TD record, the actual depths of the four elements are ½z 1 , z 2 , z 3 , z 4 ≈ ½10:8, 12:5, 14:6, 16:5m.Due to the existence of the sea current in the sea, the four elements are deviated slightly from their designed depth.When the sending ship is about 3, 5, and 7 km away from the receiving ship, the single-frequency continuous sine transmitted by the source is repeated for two times where its frequency is 120 Hz and its time length is 1 min.The serial number of the corresponding signals is signal 1 to signal 6.As there are a total of four array elements, during the processing of the transmitted signal, signals collected by two elements are selected to perform the target depth classification using the algorithm described in this paper.There are six different combinations of elements.Different combinations of elements can result in different target depth classification results.A sort is defined for the six element combinations, and the serial number of element combinations is 1 to 6, corresponding to If the parameters of the negative gradient sound velocity profile in the experiment are approximated as c 1 = 1540m/s at z = 0m and c 1 = 1520m/s at z = 22m, the simulation analyzes the effect of receiving depths on depth classification performance under negative gradient conditions, as shown in Figure 15.
The algorithm presented in this paper is applicable to the isovelocity condition as well as to the negative gradient condition.Based on the signals collected by the third and fifth element combinations, the sign distributions of the calculated I pA can be used for depth classification of both surface targets and underwater targets.The other four combinations can only accurately classify the depth of the underwater targets but cannot accurately classify the depth of the surface targets.In addition, although the critical depth of the fifth combination is smaller than that of the third combination, the sign boundary of the third combination has a smaller  13 Journal of Sensors fluctuation than the fifth combination.If the above two combinations can both meet the actual depth classification requirements, it can be considered that the depth classification effect of the third combination is better than that of the fifth combination.The change of the sound velocity profile type will have a certain impact on the depth classification performance, but the impact of the receiving depths on the depth classification performance is similar to that of the isovelocity condition.An improved depth classification effect can be obtained by appropriately adjusting the depths of the two receivers.In addition, according to the simulation results, the source used in this experiment should theoretically be identified as an underwater target.
The PSD has been achieved through smoothed FFT spectrum calculation.The Hamming window is the corresponding weighted window used for smoothing.Its window length and step are 4 s and 2 s, respectively.
The spectra of the signal 1 to signal 6 received by the four elements in the actual sea experiment are graphed in Figure 16.It can be witnessed from Figure 16 that when the horizontal distance between the source and the receiver is r ≈ 3km, each spectrum of signal 1 received by four elements has a clear line spectrum at f = 120Hz, but only the spectra of signal 2 received by the first three elements have a clear line spectrum at f = 120Hz.There is not a clear line spectrum at f = 120Hz in the signal 2 bandwidth received by the fourth element.When r ≈ 5km, only the spectra of the signals 3 and 4 received by the first element have a clear line spectrum at f = 120Hz.The spectra of the signals 3 and 4 received by the other three elements do not have a clear line spectrum at f = 120Hz.When r ≈ 7km, the spectra of signals 5 and 6 received by all the four elements do not have a clear line spectrum at f = 120Hz.
Line spectrum intensity and received SNR have a large impact on the depth classification effect.The target depth classification is performed using the signals 1 to 6, and the probability that the source is accurately discriminated as an underwater target is shown in Table 1.The serial number in the table is the serial number of element combinations.
According to the data in Table 1, when r ≈ 3km, the source can be accurately identified as an underwater target using signals collected by each of the six array element combinations.This outcome is identical to the simulation results.When r ≈ 5, 7km, signals collected by only some of the element combinations are useful for target depth classification, and the corresponding classification effects are poor.This is mainly because the source level of the low-frequency source is limited, and the amplification factor of the power amplifier is also limited.So, when the transmitted signal propagates to a distant place the signal becomes weak.The low SNR affects the performance of the target depth classification algorithm.In addition, when r ≈ 3km, the obtained depth classification results cannot prove the feasibility of the method.This method accurately classifies underwater targets and perfectly classifies surface targets at the same time.In this experiment, only one source depth was tested.
In the actual experiment data processing, there are surface targets such as fishing boats and patrol boats in the sea area.On the contrary, the source carried by the sending ship is the only underwater target which can radiate strong line spectrum signal.Therefore, by performing spectrum analysis on the received signal, if there are clear continuous line spectrum outside the time range of the underwater target radiated signals, they are the spectrum of the signals radiated by surface targets.After time-frequency analysis of the signal in the whole observation time range, it is found that there are evident continuous line spectrum corresponding to the four different segments of received signal as it can be witnessed in Figure 17.
The first segment of signal contains two line spectra whose frequencies are ½ f 1 , f 2 = ½106, 121:5Hz.The line spectrum frequency of the second signal is f 3 = 142:5Hz.There are two line spectra in the third signal, whose frequencies are ½ f 4 , f 5 = ½119:8, 159:5Hz.½ f 6 , f 7 , f 8 , f 9 = ½123, 138:5, 153:5, 169Hz are line spectrum frequencies of the fourth signal.Each segment of signal containing multiple clear continuous line spectrum may be harmonic signals radiated by the same surface target, but it may also be single or harmonic signals radiated by different surface targets.Therefore, the order of these line spectra is represented by the line spectrum number f 1 ~f9 .The corresponding depth classification results are shown in Table 2.The data in the table is the probability that the target is accurately classified as the surface target.The serial number in Table 2 is the serial number of the element combinations.
According to the data in Table 2, except for the third line spectrum, only using signals collected by the third and fifth elements can accurately identify the target as a surface target.It is only based on signals collected by these two combinations, whether it is a surface target or an underwater target, its depth can be accurately classified.Among the six combinations, the depth summation of these two combinations is closest to 29 m exactly, so that the correctness and practical application value of the ideal receiving depth concept have 14 Journal of Sensors been verified.In addition, the depth classification effect based on the signals collected by the third element combination is better than the fifth array combination.It is mainly because the depth classification effect is related to the depth summation of the element combination and related to the depth of the first element.Extremely large depth of the first element results in poor classification effect and even it is not able to classify the target depth.Considering the influence of the depth summation of the element combination and the depth of its first element on the performance of the algorithm, the depth classification effect using the third element combination is the best.The above results are consistent with the outcome obtained from the simulation results.However, using the third line spectrum signal for target depth classification, a good classification effect cannot be obtained because this segment of signal is less continuous than the other received signals, and the received SNR is lower resulting in poor performance of the source depth classification algorithm.So, accurate depth classification results cannot be obtained.
Through the analysis of the data in Table 2 and Figure 17, it can be found that the target depth classification effect is directly related to the line spectrum continuity and the received SNR.The better line spectrum continuity or the larger received signal-to-noise ratio results in a better target depth classification effect.

Conclusions
This paper analyzes the impact of receiving depth on the performance of the target depth classification algorithm through simulation and sea experiment.After a sufficient theoretical study on the effective depth model, it is verified that when the sensors are at the ideal receiving depth, the performance of the target depth classification algorithm based on the first two modes is significantly better than that of the sensors placed at the nonideal receiving depth.The depth summation of the sensors and the depth value of the first sensor will affect the depth classification effect and critical depth value.The impact of the source frequency on the performance of

Figure 2 :
Figure 2: Flow chart of target depth classification.

z 1 z 1 Figure 5 :Figure 6 :
Figure 5: Impact of z 1 + z 2 on the sign distribution of S z and I pA when z 1 = 11m.(a) Relationship between S z and z 1 + z 2 .(b) Relationship between I pA and z 1 + z 2 .

Figure 7 :Figure 8 : 3 . 4 .
Figure 7: The impact of frequency on the sign distribution of I pA under condition of isovelocity profile.(a) f = 110Hz.(b) f = 150Hz.

Figure 9 :Figure 10 :Figure 11 :
Figure 9: The impact of gradient on the algorithm performance under condition of negative gradient sound velocity profile.(a) Sound velocity profile.(b) The sign distribution of I pA .

Figure 15 :Figure 16 Figure 16 Figure 16 :
Figure 15: Relationship between the effect of depth classification and receiving depths under the condition of negative gradient profile.(a) First element combination.(b) Second element combination.(c) Third element combination.(d) Fourth element combination.(e) Fifth element combination.(f) Sixth element combination.

Figure 17 :
Figure 17: Time-frequency distributions of surface target radiated signals.(a) First segment of signal.(b) Second segment of signal.(c) Third segment of signal.(d) Fourth segment of signal.

Table 1 :
Depth classification accuracy of underwater target.

Table 2 :
Depth classification accuracy of surface target.