Prediction Method of Underwater Acoustic Transmission Loss Based on Deep Belief Net Neural Network

: The prediction of underwater acoustic transmission loss in the sea plays a key role in generating situational awareness in complex naval battles and assisting underwater operations. However, the traditional classical underwater acoustic transmission loss models do not consider the regional hydrological elements, and the performance of underwater acoustic transmission loss prediction under complex environmental conditions in a wide range of sea areas is limited. In order to solve this problem, we propose a deep learning-based underwater acoustic transmission loss prediction method. First, we studied the application domains of typical underwater acoustic transmission loss models (ray model, normal model, fast ﬁeld program model, parabolic equation model), analyzed the constraint rules of its characteristic parameters, and constructed a dataset according to the rules. Then, according to the characteristics of the dataset, we built a DBN (deep belief net) neural network model and used DBN to train and learn the dataset. Through the DBN method, the adaptation and calculation of the underwater acoustic transmission loss model under different regional hydrological elements were carried out in a simulation environment. Finally, the new method was veriﬁed with the measured transmission loss data of acoustic sea trials in a certain sea area. The results show that the RMSE error between the underwater acoustic transmission loss calculated by the new method and the measured data was less than 6.5 dB, the accuracy was higher than that of the traditional method, and the prediction speed was faster, the result was more accurate, and had a wide range of adaptability in complex seas.


Introduction
The ocean contains abundant natural resources, and all countries are intensifying their development of the ocean and enhancing the maintenance of maritime rights and interests in those territorial waters. Future naval warfare will have a high degree of informatization. Obtaining information superiority on the battlefield is the basis for winning modern wars and defending maritime rights and interests. Since the Germans published their first paper on underwater acoustics in 1919 [1], people have gradually established a theoretical system of underwater acoustic modeling and constructed a classic underwater acoustic transmission loss model, mainly including the ray model (RM), normal model (NM), fast field program (FFP) model, parabolic equation (PE) model, and some hybrid algorithms [1][2][3]. Each hydroacoustic model algorithm has limitations in the applicable situations. For example, the NM and FFP models are suitable for distance-independent sound field calculations; the RM and PE model are suitable for distance-dependent sound field calculation. In addition, the RM is suitable for high-frequency calculations, and it is more difficult to handle low-frequency and caustic calculations; the PE model is suitable for low-frequency and narrow-angle calculations, and the calculation efficiency is not very high for high-frequency and deep-sea environments. In different marine environments, 2 of 19 regional physical parameters and signal frequency ranges may change. It is necessary to select an adapted hydroacoustic model to calculate the most accurate hydroacoustic propagation loss prediction value. This is an urgent problem in the research of marine underwater acoustics. Jensen (1982Jensen ( , 1984 proposed an underwater acoustic transmission loss model matching scheme that classified five acoustic models according to the depth (deep and shallow), frequency (high and low), and the relationship between ocean environment and distance (related and unrelated) adaptation [3]. This method has no restrictions on frequency and water depth, has few judged features, and relies on a large amount of prior experience and expert knowledge for feature extraction. Based on the ray model and the parabolic equation model, Chen Zhiguo and others [4] studied the energy propagation mode of ocean fronts and vortices and made more accurate predictions for four different ocean sound fields. Li Jiaxun and others [5] designed a system with the Kraken normal model as the core including the Bellhop ray model, the Scooter finite element model, and the Sparc time-domain stepping model to calculate several typical sound field transmission losses. The above forecasting method still relies on expert experience to realize it. However, the large dynamic changes of the marine environment make the hydrological conditions have complex time-space-frequency characteristics, and the underwater acoustic field also presents complicated time-space random fluctuations, the environment is uncertain, and the parameters are uncertain. In a dynamically changing environment, it is difficult for traditional methods that rely on expert experience to accurately extract features so that the acoustic model can effectively and adaptively match the changing marine environment in real time. Therefore, when predicting the underwater acoustic transmission loss, the adaptive matching of the underwater acoustic transmission loss model needs to find a new method. In recent years, in the context of big data and deep learning, deep learning has been successfully applied to multiple research fields such as speech recognition, image processing, wireless communication, medical treatment, and social decision-making [6][7][8], bringing new ideas in solving traditional underwater acoustic transmission loss prediction problems. The deep belief net (DBN) proposed by Hinton is a neural network containing multiple hidden layers [9]. Compared with shallow models, it can effectively avoid problems such as overfitting and insufficient classification capabilities [10][11][12]. DBN is composed of several layers of restricted Boltzmann machine (RBM) and a layer of back propagation (BP) neural network [13][14][15]. It not only has a good classification effect but also has an efficient feature extraction ability. We propose a prediction method of underwater acoustic transmission loss based on deep learning. Based on the structure of the DBN network, according to the marine environment and signal characteristics, we designed a deep belief neural network model (DBN) and generated a training dataset to train the DBN neural network model and evaluate the accuracy of the model. The new method could accurately adapt the underwater acoustic transmission loss model and accurately predict the underwater acoustic transmission loss in the absence of prior knowledge, multiple frequency bands, and complex marine environmental conditions. In order to verify the practicability of the underwater acoustic transmission loss prediction method proposed in this paper, we used sea test data in a certain sea area in May 2020 to calculate the measured underwater acoustic transmission loss. We compared and analyzed the measured underwater acoustic transmission loss with the underwater acoustic transmission loss predicted by the new method. The results show that the RMSE error of the predicted transmission loss and the sea trial transmission loss was stable below 6.5 dB, which was lower than the transmission loss prediction error predicted by the classical method, which verifies the reliability and accuracy of the new method.

Principle of DBN Algorithm
DBN is a generative model of deep learning, which can be used for unsupervised learning and as an autoencoder, and can also be used for supervised learning and used as a classifier.
From the perspective of unsupervised learning, its purpose is to retain the characteristics of the original features as much as possible while reducing the dimension of the features; from the perspective of supervised learning, its purpose is to make the classification error rate as small as possible. As the basic building block of DBN, RBM is a probabilistic modeling method based on energy function. It has strong unsupervised learning ability and is a key model in the process of constructing and training DBN. The model performs continuous state transitions on the input features of the sample, so that the model can finally fit the distribution characteristics of the sample set, can learn arbitrarily complex rules in the data, and is suitable for the training and learning of large-scale dynamic data features in complex sea areas. The core of the algorithm is to optimize the connection weights between layers with a layer-by-layer greedy learning algorithm [16][17][18]. As shown in Figure 1, each RBM contains a pair of visible layers (v 1 , v 2 , . . . , v m ) and hidden layer (h 1 , h 2 , . . . , h n ), where each layer is connected by a fully connected layer [19] and w is the weight value of the connection between the two layers. It is a random neural network model without feedback.

Principle of DBN Algorithm
DBN is a generative model of deep learning, which can be used for unsupervised learning and as an autoencoder, and can also be used for supervised learning and used as a classifier. From the perspective of unsupervised learning, its purpose is to retain the characteristics of the original features as much as possible while reducing the dimension of the features; from the perspective of supervised learning, its purpose is to make the classification error rate as small as possible. As the basic building block of DBN, RBM is a probabilistic modeling method based on energy function. It has strong unsupervised learning ability and is a key model in the process of constructing and training DBN. The model performs continuous state transitions on the input features of the sample, so that the model can finally fit the distribution characteristics of the sample set, can learn arbitrarily complex rules in the data, and is suitable for the training and learning of large-scale dynamic data features in complex sea areas. The core of the algorithm is to optimize the connection weights between layers with a layer-by-layer greedy learning algorithm [16][17][18]. As shown in Figure 1, each RBM contains a pair of visible layers ( 1, 2, … , ) and hidden layer (ℎ1, ℎ2, … , ℎ ), where each layer is connected by a fully connected layer [19] and w is the weight value of the connection between the two layers. It is a random neural network model without feedback. Suppose that in a restricted Boltzmann machine network, the number of visible units is m, the state of visible units set in a certain layer is represented by a vector v, and the state of hidden units in a certain layer is represented by a vector h. and ℎ represent the states of the i visible unit and the j hidden unit, respectively. When the state of the nodes in the visible layer and the hidden layer is determined, the energy function of the RBM [20] can be expressed as: in which = ( , , ) is the parameter of RBM: is the connection weight between the visible unit i and the hidden unit j; is the bias of the visible unit i; and is the bias of the hidden unit j. When is known, since the state of RBM obeys the regular distribution, the joint probability distribution of any set of visible and hidden unit node states ( , ℎ) is: Suppose that in a restricted Boltzmann machine network, the number of visible units is m, the state of visible units set in a certain layer is represented by a vector v, and the state of hidden units in a certain layer is represented by a vector h. V i and h j represent the states of the i visible unit and the j hidden unit, respectively. When the state of the nodes in the visible layer and the hidden layer is determined, the energy function of the RBM [20] can be expressed as: in which θ = (w, a, b) is the parameter of RBM: W ij is the connection weight between the visible unit i and the hidden unit j; a i is the bias of the visible unit i; and b j is the bias of the hidden unit j. When θ is known, since the state of RBM obeys the regular distribution, the joint probability distribution of any set of visible and hidden unit node states (v, h) is: in which, Z(θ) = ∑ v,h e −E(v,h|θ) is the normalization factor. For the distribution P(v|θ) of the input data v, it can be determined by the likelihood function of the joint probability distribution P(v, h|θ) : Because the connection between RBM layers is restricted, regardless of whether it is a visible layer or a hidden layer, there is no connection between nodes in the layer, and its activation state is also independent. Therefore, when the visible unit is given, the activation probability of the j hidden unit is: 1+e −x is the sigmoid activation function. Similarly, when the hidden unit is given, the activation probability of the i visible element is: When the state v of the visible layer node is known, the above formula can be used to calculate P(h|v) to determine the state of the hidden layer h; then the state of the visible layer V 1 can be obtained by calculating P(v|h).
Take S = {v 1 , v 2 , . . . , v s } (capacity s) as the training sample, and obtain the model parameter θ by maximizing the log likelihood function L(θ) of the RBM on the sample and fit the training sample, The hidden layer can be regarded as a feature of the input data of the visual layer.
The log-likelihood function is differentiated by the contrast divergence (CD) [14] algorithm. For example, the current weight is w ij : Gibbs sampling is used to determine the expectation E v i h j ,data of the training samples and the expectation E v i h j ,model defined by the RBM model. Up to now, w i can be updated as follows: In the formula, ϑ is the momentum and η is the learning rate. Similarly, a j and b j can also be updated in this way.

DBN Model Structure
The structure of the DBN model constructed in this paper is shown in Figure 2. The first visible layer and the second hidden layer constitute the first RBM, and the second hidden layer and the third hidden layer form the second RBM, by analogy, a total of five RBM are formed. Input the data of the fifth RBM into a supervised classifier, and use BP to fine-tune the relevant parameters of the DBN to form a complete DBN to classify the input data. X = {x 1 , x 2 , . . . , x n } is the input of the first visual layer, representing the current marine environment and signal characteristic parameters (as shown in Table 1). h x,y are the states of y unit in the x hidden layer; {x 1 , x 2 , x 3 , x 4 , x 5 } are the connection weight vectors between adjacent layers. The top-level classifier network structure uses the BP backpropagation neural network and w bp is the weight vector of the BP neural network. If the output result cannot get the expected tag data, the BP backpropagation adjustment is used to continuously approach the expected output of the representative label. agation neural network and is the weight vector of the BP neural network. If the output result cannot get the expected tag data, the BP backpropagation adjustment is used to continuously approach the expected output of the representative label.   The thickness of the deposited layer As shown in Figure 3, the training DBN was divided into three stages. The first stage is the model initialization stage. Regarding any two adjacent layers in the DBN as an RBM, input the output vector of the previous RBM to the next RBM, and initialize the parameters in the DBN layer by layer. The second stage is the unsupervised forward learning stage. The preprocessed marine environment features and signal features are input to the first visual layer. The greedy unsupervised learning method is used to train layer by layer,  The thickness of the deposited layer As shown in Figure 3, the training DBN was divided into three stages. The first stage is the model initialization stage. Regarding any two adjacent layers in the DBN as an RBM, input the output vector of the previous RBM to the next RBM, and initialize the parameters in the DBN layer by layer. The second stage is the unsupervised forward learning stage. The preprocessed marine environment features and signal features are input to the first visual layer. The greedy unsupervised learning method is used to train layer by layer, and the activation function of each unit is sigmoid. In the second stage, the shallow original features are expressed at a high level. Since the unsupervised algorithm cannot optimize the entire network parameters, it is necessary to enter the third stage, the reverse supervised learning stage, which uses the BP supervised algorithm to reversely fine-tune the DBN to optimize the underwater acoustic transmission loss prediction under the large-scale marine environment and the changeable signal frequency band. This training mode directly maps data from input to output through unsupervised training, shortening the time of supervised training, and is also the key to a powerful feature rule extraction capability. mize the entire network parameters, it is necessary to enter the third stage, the reverse supervised learning stage, which uses the BP supervised algorithm to reversely fine-tune the DBN to optimize the underwater acoustic transmission loss prediction under the large-scale marine environment and the changeable signal frequency band. This training mode directly maps data from input to output through unsupervised training, shortening the time of supervised training, and is also the key to a powerful feature rule extraction capability.

Model initialization
Unsupervised learning Supervised learning The first stage The second stage The third phase

Feature Parameter Rules and Sample Set Construction
This article mainly studied the prediction method of underwater acoustic transmission loss, so it focused on the influence of different underwater acoustic characteristic parameters on the prediction results of the acoustic propagation model. The problem of sound transmission under water is mainly studied by two methods: one method is the wave theory, which studies the changes in the amplitude and phase of the acoustic signal in the sound field; and the other method is the ray theory, which can be used in a highfrequency range, sound wave is regarded as a beam of rays, and the change of the sound intensity in the sound field with the beam of rays is usually studied. From these two theories, several classic underwater acoustic transmission models have been developed to characterize the underwater acoustic propagation process including NM, RM, PE model, etc. [3].
In underwater acoustics, the standard measure of sound field signal strength versus distance is the transmission loss (TL), which is the ratio of the sound intensity ( , ) at a certain point in the sound field to the sound intensity I0 at a distance of 1 m from the sound source. It can be expressed as Equation (9): where ( , ) is the sound pressure at a certain point in the sound field; r is the distance from the sound source position to a certain point in the sound field; z is the receiving depth of a certain point in the sound field; 0 P is the sound pressure value at a distance of 1 m from the sound source; and TL is the value of the transmission loss of the sound field in decibels (dB). This paper selected the Bellhop model in ray mode for research. This model was written by Porter and Bueker [21] in 1987 using the Gaussian beam tracking method to calculate the sound field in a horizontally non-uniform environment (conditions related to distance). The experimental data of this model in the frequency range of 600 Hz-30 KHz were consistent with the theoretical model [22]. It is suitable for near-field conditions. The important characteristic length of the water depth or sound velocity profile should be greater than 20 wavelengths [23].

Feature Parameter Rules and Sample Set Construction
This article mainly studied the prediction method of underwater acoustic transmission loss, so it focused on the influence of different underwater acoustic characteristic parameters on the prediction results of the acoustic propagation model. The problem of sound transmission under water is mainly studied by two methods: one method is the wave theory, which studies the changes in the amplitude and phase of the acoustic signal in the sound field; and the other method is the ray theory, which can be used in a high-frequency range, sound wave is regarded as a beam of rays, and the change of the sound intensity in the sound field with the beam of rays is usually studied. From these two theories, several classic underwater acoustic transmission models have been developed to characterize the underwater acoustic propagation process including NM, RM, PE model, etc. [3].
In underwater acoustics, the standard measure of sound field signal strength versus distance is the transmission loss (TL), which is the ratio of the sound intensity I(r, z) at a certain point in the sound field to the sound intensity I 0 at a distance of 1 m from the sound source. It can be expressed as Equation (9): where p(r, z) is the sound pressure at a certain point in the sound field; r is the distance from the sound source position to a certain point in the sound field; z is the receiving depth of a certain point in the sound field; P 0 is the sound pressure value at a distance of 1 m from the sound source; and TL is the value of the transmission loss of the sound field in decibels (dB). This paper selected the Bellhop model in ray mode for research. This model was written by Porter and Bueker [21] in 1987 using the Gaussian beam tracking method to calculate the sound field in a horizontally non-uniform environment (conditions related to distance). The experimental data of this model in the frequency range of 600 Hz-30 KHz were consistent with the theoretical model [22]. It is suitable for near-field conditions. The important characteristic length of the water depth or sound velocity profile should be greater than 20 wavelengths [23].
For the normal model, this paper selected the Kraken model for research. This model was jointly developed by the U.S. Naval Ocean System Center (NOSC) and the U.S. Naval Research Laboratory (NRL) using the finite difference method to solve the normal wave equation. It is suitable for low-frequency (lower than 500 Hz) shallow sea and farfield calculation under uniform horizontal (distance-independent) marine environmental conditions, and the seabed structure must be given in the calculation, that is, the parameters of Cs, Kd, p, and Mt in Table 1.
For the fast field program model (FFP model), we choose the Scooter model for research. This model applies fast Fourier transform (FFT) technology to solve the wave equation. It is limited to the condition of a horizontally layered homogeneous medium that is not related to distance, near-field environment, and is suitable for calculations in the low-frequency range (below 500 Hz).
For the parabolic equation (PE) model, the RAMGeo model was selected. This model is an approximate solution of the wave equation with a narrow angle (<±20 • ) considering the coupling effect of acoustic diffraction and normal waves of various numbers. It is suitable for calculating narrow angles. The calculation speed for low frequency (lower than 500 Hz) and shallow sea is very fast. When the frequency increases and the depth increases, the calculation time increases rapidly, therefore, this model is not applicable.
The above far-field and near-field rules are constrained by the characteristic parameter horizontal distance (R). Suppose that the total spacing of the uniform linear array element D = m × d; m is the number of array element spacing; the highest frequency wavelength of the sound source (the minimum wavelength of the sound source) is λ min , if the sound source is at the level of the array Distance R > 2D 2 λ min , it is the far-field model, otherwise it is the near-field model. Under the near-field model, sound waves are regarded as spherical waves, and under the far-field model, sound waves are regarded as plane waves. The calculation of underwater acoustic transmission loss mainly includes characteristic parameters such as signal frequency, sea depth, and seabed environment. The specific parameters and parameter ranges are shown in Table 1.
The optimal characteristic parameter ranges of the above four types of underwater acoustic transmission loss models are shown in Table 2. Sample data generation is very important for neural network training, and the distribution and statistical characteristics of the sample will also directly affect the training efficiency and prediction effect of the neural network. Based on the above constraints, we built a dataset with a normal distribution to train the neural network model. The training set generation method is as follows: (1) Select an underwater acoustic propagation model, and generate the data on the principle of a normal distribution within the constraint range according to the abovementioned constraint conditions of the parameters of the selected model. (2) Set labels for the samples in the established training set. The labels of each sample are different. Usually, a combination of 0-K (K is the number of models) is used to set the corresponding label. The labels are formed into a new set corresponding to the input data. (3) Combine different marine environment characteristics and signal characteristics according to the above steps to form a large training set. Then, after making a one-to-one correspondence between the input data and the label, randomly shuffle the order of the input data, and recombine the training set. The purpose of shuffling the data is to prevent too much noise data of the same type from clustering together, causing over-fitting.
Different feature parameters have different weights in the adaptation model. In the DBN neural network, the weight of each feature parameter is automatically optimized according to the training data, and finally, the optimal weight is achieved through iterative training. If we want the neural network to learn these rules in-depth, we will build more data for the network's learning of the feature parameters with larger weight, and the rest of the feature parameters can be iterated by the neural network to obtain the optimal weight. As shown in Table 2, we extracted a total of nine feature parameters. Compared with other models, Fs is a significant feature parameter in the ray mode and has a larger weight, so we increased the number of Fs parameters in the ray mode training. For the ray model, we increased the number of samples for the feature parameter Fs with a larger weight, and constructed 50 different values within the constraint range. The remaining eight different characteristic parameters were used as a group, and a total of 20 groups of different values were generated within the constraint range. Through the combination of the characteristic parameter Fs and other characteristic parameters, the total number of data samples of the ray model was 1000. For normal mode, it can be seen that H and R had a larger weight. We first generated 50 data with different values for the characteristic parameter H, and the remaining eight characteristic parameters generated 20 sets of sample data within the constraint range, and combined them to generate 1000 sets of data. Then, we generated 50 different values for the characteristic parameter R, and the remaining eight characteristic parameters generated 20 sets of sample data within the constraint range, and combined them to generate 1000 sets of data. A total of 2000 sets of sample data were generated for training the normal wave model data. For the fast field model, like the normal model, the H parameter and the R parameter had greater weight, but the constraint range was different, so 2000 sets of sample data were also generated. For the parabolic model (PE), it can be seen that the characteristic parameter SA had a greater weight, the same as the ray model, and 1000 sets of data were generated. Finally, 6000 sets of sample data were generated for DBN training.

Model Initialization and Training
First, we initialized the model to determine the momentum factor, learning rate, maximum number of iterations, and the number of hidden layer nodes related to the RBM training process. In addition, because DBN does not clearly define the most suitable layer number, the determination of the layer number needs to determine the depth of the DBN model through experimental evaluation results. Then, we set the relevant parameters of RBM, set the initial value of learning rate to 0.1 according to experience, and set the initial momentum to 0.5, where the number of input nodes is the number of characteristic parameters in Table 1, and the number of output nodes is four nodes equal to the number of acoustic models. The number of pre-training iterations (epochs) was selected as 200, and the number of global reverse fine-tuning was 100. In order to prevent over-fitting, we randomly selected 3/4 of each type of data as training samples, and the rest as test samples. The number of the random sample training batch size was 20. Table 3 shows the performance of the DBN neural network at different depths. It can be seen that as the depth of the DBN model increases, the time cost also doubles. When the depth is greater than 5, the classification accuracy does not increase significantly. The training time has been increasing by a large margin, so for comprehensive training time and classification accuracy, five layers is a more reasonable depth choice. The training of DBN mainly consists of two parts. The first part is unsupervised bottom-up layer-by-layer pre-training, which transfers the feature vector from the input to the output. The second part is the supervised top-down reverse fine-tuning. The output value is compared with the given data label to get the error, and the error is back propagated from the output end to the input end to achieve the purpose of adjusting the DBN parameters. The overall training process is shown in Figure 4: The training of DBN mainly consists of two parts. The first part is unsupervised bottom-up layer-by-layer pre-training, which transfers the feature vector from the input to the output. The second part is the supervised top-down reverse fine-tuning. The output value is compared with the given data label to get the error, and the error is back propagated from the output end to the input end to achieve the purpose of adjusting the DBN parameters. The overall training process is shown in Figure 4: The specific training steps are as follows: Step 1: Preprocess the dataset constructed in Section 3.1 and input data into the model in batches for training.
Step 2: Set the bias of the visible layer and the hidden layer and the initial value of the weight between each neuron = { , , }. represents the bias of the input data layer or visual layer, represents the bias of each hidden layer or data output layer, and represents the connection weight of each layer. According to the research of Hinton et al., the connection weight of each layer is initialized to obey the normal distribution N (0,1), and the initial value of the bias of the visible layer and the bias of the hidden layer is set to 0.
Step 3: Next, we perform iterative training on DBN. In the DBN model, each layer of RBM learns the characteristics of the input data forward and unsupervised, and uses the greedy unsupervised learning algorithm to update the bias value of each unit of each layer and the weight between the hidden layer.
Step 4: After completing the network pre-training, we add the Softmax output layer, and use fine-tuning to transform the entire network from a generative model to a discriminative model. This network loss function can be defined as: The specific training steps are as follows: Step 1: Preprocess the dataset constructed in Section 3.1 and input data into the model in batches for training.
Step 2: Set the bias of the visible layer and the hidden layer and the initial value of the weight between each neuron θ = {b, c, w}. b represents the bias of the input data layer or visual layer, c represents the bias of each hidden layer or data output layer, and w represents the connection weight of each layer. According to the research of Hinton et al., the connection weight w of each layer is initialized to obey the normal distribution N (0,1), and the initial value of the bias b of the visible layer and the bias c of the hidden layer is set to 0.
Step 3: Next, we perform iterative training on DBN. In the DBN model, each layer of RBM learns the characteristics of the input data forward and unsupervised, and uses the greedy unsupervised learning algorithm to update the bias value of each unit of each layer and the weight between the hidden layer.
Step 4: After completing the network pre-training, we add the Softmax output layer, and use fine-tuning to transform the entire network from a generative model to a discriminative model. This network loss function can be defined as: In the formula, 1{label(n) = 1} means that the condition in brackets is set to 1 when the condition is satisfied, and 0 is not satisfied. W 0 (n) and b (n) are the connection weight vector of the nth node in the output layer and its offset; x is the response of the upper node. The fine-tuning process generally uses the gradient descent method, the goal is to minimize the loss function and use the error back propagation (BP) to fine-tune the network. The fine-tune process is to compare the correct classification labels corresponding to different feature parameters (Table 1) in turn. If the expected output value is not obtained, the difference between the actual output and the expected output (error) is calculated, and the error is successively modified by modifying the weight of each neuron propagate to the input layer.
Step 5: When iterations end, the training of the entire model is completed.
After the training of the DBN model is completed, all training samples can be accurately classified under the correct label assigned, that is, the error between the predicted label and the real label converges to 0. Use loss L(λ) to evaluate the degree of convergence of the model: Loss tending to 0 means that the model converges well. As shown in the loss curve of Figure 5, as the number of iterations increases, the loss tends to 0.
to different feature parameters (Table 1) in turn. If the expected output value is not ob-tained, the difference between the actual output and the expected output (error) is calculated, and the error is successively modified by modifying the weight of each neuron propagate to the input layer.
Step 5: When iterations end, the training of the entire model is completed.
After the training of the DBN model is completed, all training samples can be accurately classified under the correct label assigned, that is, the error between the predicted label and the real label converges to 0. Use loss ( ) to evaluate the degree of convergence of the model: Loss tending to 0 means that the model converges well. As shown in the loss curve of Figure 5, as the number of iterations increases, the loss tends to 0. This training uses 75% of the data to train the model, and its loss function and accuracy are shown in Figure 5. It can be seen from the figure that as the number of training increases, the model prediction accuracy gets higher and higher, and the error loss curve also declines, and finally remains at a stable value. At this time, an increase in the number of training may lead to overfitting. Finally, the 25% data in the sample set is used for testing, and the accuracy of model adaptation is obtained, which is 94.86% in combination. It shows that the recognition accuracy is already very high, close to the distribution of the sample data, and then the underwater acoustic transmission loss is calculated through the model matching and adaptive underwater acoustic propagation model.

Sources of Data
To verify the accuracy of the new forecasting method, we conducted a sea trial experiment in May 2020 to analyze and process the obtained sea trial data to obtain the actual underwater acoustic transmission loss data in this sea area. Then, we compared and This training uses 75% of the data to train the model, and its loss function and accuracy are shown in Figure 5. It can be seen from the figure that as the number of training increases, the model prediction accuracy gets higher and higher, and the error loss curve also declines, and finally remains at a stable value. At this time, an increase in the number of training may lead to overfitting. Finally, the 25% data in the sample set is used for testing, and the accuracy of model adaptation is obtained, which is 94.86% in combination. It shows that the recognition accuracy is already very high, close to the distribution of the sample data, and then the underwater acoustic transmission loss is calculated through the model matching and adaptive underwater acoustic propagation model.

Sources of Data
To verify the accuracy of the new forecasting method, we conducted a sea trial experiment in May 2020 to analyze and process the obtained sea trial data to obtain the actual underwater acoustic transmission loss data in this sea area. Then, we compared and analyzed the measured transmission loss data with the data predicted by this method. The sea trial site is shown in Figure 6.
The experimental area is shown in Figure 7. The experimental ship traveled from Q1 ( Figure 7) to Q2 (Figure 7) for the experiments. The ship carried the towed array sonar and continuously emitted a set of signals, and the 20-element vertical receiving array at the ZS (Figure 7) station received the signals emitted by the towed array sonar. We placed a CTD at the receiving array for sound velocity profile measurement. The measured sound velocity profile is shown in Figure 7. The planned total voyage of the experimental ship was 160 km, the starting Q1 was 274 km away from the receiving array, and the ending Q2 point was 420 km away from the receiving array. Among them, we selected one of the sea trial data with a voyage to establish transmission loss data for verification. analyzed the measured transmission loss data with the data predicted by this method. The sea trial site is shown in Figure 6. The experimental area is shown in Figure 7. The experimental ship traveled from Q1 ( Figure 7) to Q2 (Figure 7) for the experiments. The ship carried the towed array sonar and continuously emitted a set of signals, and the 20-element vertical receiving array at the ZS (Figure 7) station received the signals emitted by the towed array sonar. We placed a CTD at the receiving array for sound velocity profile measurement. The measured sound velocity profile is shown in Figure 7. The planned total voyage of the experimental ship was 160 km, the starting Q1 was 274 km away from the receiving array, and the ending Q2 point was 420 km away from the receiving array. Among them, we selected one of the sea trial data with a voyage to establish transmission loss data for verification. The signals continuously transmitted by the experimental ship during the experiment are shown in Table 4, where Continuous Wave (CW) represents a single-frequency signal, and Linear Frequency Modulated (LFM) represents a frequency-modulation signal. In each group of signals, single-frequency or FM signals are transmitted according to the law shown in Table 4, and a group of signals is transmitted with a period of 160 S. As  The experimental area is shown in Figure 7. The experimental ship traveled from Q ( Figure 7) to Q2 (Figure 7) for the experiments. The ship carried the towed array sona and continuously emitted a set of signals, and the 20-element vertical receiving array a the ZS (Figure 7) station received the signals emitted by the towed array sonar. We place a CTD at the receiving array for sound velocity profile measurement. The measured soun velocity profile is shown in Figure 7. The planned total voyage of the experimental ship was 160 km, the starting Q1 was 274 km away from the receiving array, and the endin Q2 point was 420 km away from the receiving array. Among them, we selected one of th sea trial data with a voyage to establish transmission loss data for verification. The signals continuously transmitted by the experimental ship during the experi ment are shown in Table 4, where Continuous Wave (CW) represents a single-frequenc signal, and Linear Frequency Modulated (LFM) represents a frequency-modulation signa In each group of signals, single-frequency or FM signals are transmitted according to th law shown in Table 4, and a group of signals is transmitted with a period of 160 S. A  Table 4, where Continuous Wave (CW) represents a single-frequency signal, and Linear Frequency Modulated (LFM) represents a frequency-modulation signal. In each group of signals, single-frequency or FM signals are transmitted according to the law shown in Table 4, and a group of signals is transmitted with a period of 160 S. As shown in Figure 8, it is a time-frequency domain diagram of the transmitted signal received at a certain moment for the 311 m depth array element in the receiving array. As shown in Figure 8a, although it is impossible to analyze useful information on the time domain diagram, it can be seen in the frequency domain diagram that there are significantly higher peaks at the frequency of 360 Hz than the other frequency points. This phenomenon indicates that the 360 Hz signal transmitted by the experimental ship is received on the receiving array element, and the energy intensity of the received signal can be obtained in the frequency domain. Similarly, as shown in Figure 8b, the 320-400 Hz LFM signal can be found in the frequency domain, and its center frequency can be selected to obtain the energy intensity of the received signal. The remaining CW and LFM signals use the same method to obtain the received signal energy intensity. Next, according to the known transmitted signal energy intensity and received signal energy intensity, we can calculate the actual underwater acoustic transmission loss for this distance and generate the underwater acoustic transmission loss dataset.   Based on the transmitted signals in Table 4, this paper constructed a sea triainl transmission loss dataset. The sound source level (SL) of the CW signal was 192 dB, the SL of the LFM signal was 174 dB, and the emission depth was 140 m. As shown in Figure 9, it can be clearly seen in the LOFAR spectrum that the high-brightness area represents the transmitted signal received at the current moment. The LOFAR spectrum of the 360 Hz CW signal is shown in Figure 9a. It can be seen that the duration was 10 S, which was the same as the description in Table 4, indicating that the transmitted CW signal was received at this time. The LOFAR spectrum of the 320-400 Hz LFM signal is shown in Figure 9b. It can be seen that its duration was also 10 S, which was the same as the description in Table 4, indicating that the transmitted LFM signal was received at this time. The rest of the signal can be found in the same way, and the received energy can be analyzed. For this experiment, the calculation of the center frequency of the CW signal and the LFM signal was used. According to the definition of the sound source level and signal level, the transmission loss calculated by Equation (9) can be simplified as: where TL i represents the transmission loss of the i transmitted pulse, and represents the signal level of the transmitted pulse at the receiving element end. At each signal frequency, each array element collects 30 data points, so each signal frequency forms a sea trial transmission loss data set of 250 data points. The five transmitted signals form a sea trial transmission loss dataset of 1250 data points in total.

Sea Area Characteristic Parameters
The currently known environmental parameters are as follows: the depth of the sound source was 141 m, the distance between the sound source and the receiving array was 274~340 km, and 200 sampling points were taken at equal intervals in this range. The signal frequency was 360 Hz, the sea depth was 2000 m, the sedimentary layer thickness was 50 m, the sea bottom was sand-mud-clay, that is, the sedimentary sound velocity was 1579 m/s, the sedimentary density was 1.596 g/cm 3 , and the sedimentary attenuation was 0.113 dB/m-kHz. We placed a CTD at the receiving array to measure the sound velocity profile. The measured sound velocity profile is shown in Figure 10.

Sea Area Characteristic Parameters
The currently known environmental parameters are as follows: the depth of the sound source was 141 m, the distance between the sound source and the receiving array was 274~340 km, and 200 sampling points were taken at equal intervals in this range. The signal frequency was 360 Hz, the sea depth was 2000 m, the sedimentary layer thickness was 50 m, the sea bottom was sand-mud-clay, that is, the sedimentary sound velocity was 1579 m/s, the sedimentary density was 1.596 g/cm 3 , and the sedimentary attenuation was 0.113 dB/m-kHz. We placed a CTD at the receiving array to measure the sound velocity profile. The measured sound velocity profile is shown in Figure 10.
sound source was 141 m, the distance between the sound source and the receiving array was 274~340 km, and 200 sampling points were taken at equal intervals in this range. The signal frequency was 360 Hz, the sea depth was 2000 m, the sedimentary layer thickness was 50 m, the sea bottom was sand-mud-clay, that is, the sedimentary sound velocity was 1579 m/s, the sedimentary density was 1.596 g/cm 3 , and the sedimentary attenuation was 0.113 dB/m-kHz. We placed a CTD at the receiving array to measure the sound velocity profile. The measured sound velocity profile is shown in Figure 10.  Figure 10 shows that it is a typical deep-sea sound speed profile curve. It can be seen that below 1000 m, the temperature distribution was mainly affected by sunlight, which presents a relatively stable negative gradient. Above 1000 m, the seawater temperature is constant, which is called the deep-sea isothermal layer. The sound velocity distribution of  Figure 10 shows that it is a typical deep-sea sound speed profile curve. It can be seen that below 1000 m, the temperature distribution was mainly affected by sunlight, which presents a relatively stable negative gradient. Above 1000 m, the seawater temperature is constant, which is called the deep-sea isothermal layer. The sound velocity distribution of this layer mainly depends on the increase in static pressure, presenting an almost linear positive gradient (0.16 m/s/m) distribution.

Forecast Result Analysis
To make a more scientific assessment of the overall performance of the underwater acoustic transmission loss prediction by the new method, we applied the root mean square error (RMSE) to evaluate the transmission loss and the transmission loss of the underwater acoustic model mismatch prediction by this method. The actual transmission loss of sea trials was compared with errors to evaluate the forecast accuracy of this method. The specific definition is as Equation (12): whereŷ i represents the predicted transmission loss; y i represents the actual transmission loss; and N represents the predicted transmission loss point. The smaller the value of RMSE, the smaller the difference between the predicted transmission loss data and the sea trial transmission loss data, that is, the more accurate the forecasting ability. According to the above-mentioned marine environmental characteristic parameters, we selected the characteristic parameters mentioned in Table 1, and input them into the DBN neural network constructed in this article for underwater acoustic model adaptation and underwater acoustic transmission loss prediction. As shown in Figure 10, Comparing the prediction results of the model selected by the new method with the prediction results of other underwater acoustic models, it can be seen that the prediction results of the model selected by the new method were the most consistent. In the figure, the x-axis is the transmission distance in m, and the y-axis is the energy intensity in dB. Figure 11a shows that the ocean transmission loss predicted by the new method using this system is basically in line with the actual transmission loss. Although there are error fluctuations between the predicted curve and the actual curve, it can be seen in Figure 11 that the system has the smallest prediction error. The calculated RMSE error is shown in Comparing Figure 11c,d, it is obvious that these two models are not suitable for the current marine environmental conditions. If this model is selected incorrectly for calculation, it will cause the calculation result to deviate too much from the actual result.
The results verified by actual data showed that the new method had higher matching for the model selection of transmission loss prediction, and the prediction of transmission loss had higher accuracy.     Figure 11b shows that the transmission loss curve predicted by Bellhop had a certain error deviation from the actual transmission loss curve. The calculated RMSE error is shown in Table 5. The average RMSE error of the transmission loss data under each array element was 11.56 dB.
Comparing Figure 11a with Figure 11b, it can be seen that the predicted transmission loss curve was actually very close, and in Figure 12, it can be seen that the error curve was also very close. In this case, in actual application, if there is no rich expert experience, it is difficult to choose the best matching underwater acoustic propagation model. The RMSE of the transmission loss result predicted by this method was compared with the transmission loss RMSE predicted by Bellhop, and the error was reduced by 5.25 dB. The results show that this method can easily match a suitable underwater acoustic propagation model and calculate accurate underwater acoustic transmission loss prediction in a complex ocean environment.
(c) RAMGeo prediction and actual transmission loss (d) Scooter prediction and actual transmission loss Figure 11. Comparison of system prediction TL and sea trial TL. Figure 12. Comparison of absolute error between system prediction and classic model prediction.
As shown in Figure 11c, it can be seen that there were obvious errors between the actual sea test data and the data predicted by the RAMGeo model, and there was a model mismatch. Figure 12 shows that the error curve was obviously too large. The calculated RMSE error is shown in Table 5, and the average RMSE error of the transmission loss data under each array element was 17.594 dB.
As shown in Figure 11d, it can be seen that there was a large error between the actual sea test data and the data predicted by the Scooter model. It is obvious that this model does not fit under the current marine environmental conditions, and the calculated average RMSE error was 55.334 dB.
Comparing Figure 11c,d, it is obvious that these two models are not suitable for the current marine environmental conditions. If this model is selected incorrectly for calculation, it will cause the calculation result to deviate too much from the actual result.
The results verified by actual data showed that the new method had higher matching for the model selection of transmission loss prediction, and the prediction of transmission loss had higher accuracy.

Discussion
Aiming at the prediction of underwater acoustic transmission loss, we proposed a new method for the prediction of underwater acoustic transmission loss. Under ideal marine environment and constraint conditions, we constructed sample data to train and test the accuracy of the new method and finally verified by sea trial transmission loss data. Under the ideal marine environment and constraint conditions, the parameter changes are controllable, the model applicable parameter boundaries are clear, and simulation testing is easy. In the simulation test, the adaptation accuracy of the underwater acoustic transmission loss model in the new method reached 94.86%, indicating that the recognition rate is already very high. Under constrained conditions, the sample dataset generated based on the normal distribution conforms to the actual marine environment parameter distribution, which is conducive to training stable and reliable results.
Compared with the ideal environment condition where the characteristic parameter boundary during simulation training is clear, the characteristic parameter of the sea test environment has a certain degree of complexity, and the boundary overlaps each other, which increases the difficulty in adapting the acoustic model. Without in-depth research on marine acoustics, it is difficult to quickly adapt a suitable underwater acoustic transmission loss model. Using the new method proposed in this paper, only the currently measurable characteristic parameters of the marine environment (shown in Table 1) are needed, then a suitable underwater acoustic transmission loss model can be quickly matched to calculate the transmission loss prediction. We used the measured data of a certain sea area in 2020 to verify, calculate the measured transmission loss under different depth array elements, and compare it with the transmission loss predicted by the new method. Finally, all underwater acoustic transmission loss models were used for calculation, and RMSE was used to evaluate the accuracy. The RMSE results are shown in Tables 5-9. It can be seen that the transmission loss error of the acoustic model selected by the new method under different signals and elements was the smallest. In practical applications, there are still many complicated marine environments, and the new method needs to be tested for reliability in more areas.

Conclusions
In this paper, we constructed a new prediction underwater acoustic transmission loss method based on the DBN neural network. The new method solves the problem that the traditional underwater acoustic transmission loss prediction cannot accurately and adaptively select the underwater acoustic transmission loss model in the complex and changeable marine environment.
First of all, based on the characteristics of the DBN network, we designed a DBN model according to the marine environment and signal characteristics. We studied the applicable boundary conditions of various underwater acoustic transmission loss models ray model, normal model, fast field program model, parabolic equation model). Based on these conditions, we constructed rules to generate datasets for training and testing DBN neural networks. Then, we evaluated the accuracy of the DBN neural network model and obtained the accuracy of the model adaptation, which was 94.86%. It shows that the recognition rate was already very high. Finally, the new method was verified by a certain sea trial data in 2020. The results showed that the new method achieved accurate prediction of the underwater acoustic propagation model in a complex ocean environment, and the average RMSE error of the prediction results under different array elements was less than 6.5 dB. The new method reduces the difficulty in matching underwater acoustic transmission loss models in complex marine environments and improves the speed and accuracy of underwater acoustic transmission loss calculations. The new method can accurately calculate the transmission loss of multi-band signals under different marine environmental conditions and has the generalization ability for the calculation of transmission loss in a wide range of sea areas.