Modelling and Simulation of Pseudo-Noise Sequence-Based Underwater Acoustic OSDM Communication System

underwater acoustic communication applications. Abstract: Orthogonal signal division multiplex (OSDM) is an emerging signal modulation technology which has a lower peak-to-average power ratio (PAPR) and a ﬂexible subcarrier system architecture. Particularly, it can be seen as a bridge between the single-carrier modulation and the orthogonal frequency division multiplexing (OFDM) modulation in the frequency domain. Aiming at the development trend and demand of underwater acoustic hybrid and adaptive modulation communication technology, a pseudo-noise (PN) sequence-based underwater acoustic OSDM communication system is proposed in this paper. A data frame structure with PN sequence is designed to solve the multipath and Doppler e ﬀ ect of underwater acoustic channel. On the basis of the PN sequence, a compressive sensing method based on the orthogonal matching pursuit (OMP) algorithm and the minimum mean square error (MMSE) algorithm is designed for channel estimation and equalization. On the basis of the system construction, the relationship among the OSDM vector length M, the OSDM subcarrier number N, and the underwater acoustic channel length is further studied for adaptive modulation of underwater acoustic communication. Finally, the proposed system is veriﬁed by simulation. The OSDM system has lower and controllable PAPR. When the OSDM vector length M is bigger than the channel length, and the system subcarrier ﬂexibility is guaranteed, the bit error rate (BER) of the OSDM system is lower than that of the OFDM system and the single-carrier system. The PN sequence-based compressive sensing channel estimation and equalization with the OMP and MMSE algorithms has a good performance to resist the multipath e ﬀ ect of underwater acoustic channel.

single-carrier communication system. In general, the emerging OSDM technology is an interesting modulation scheme which can realize SC-FDE modulation, OFDM modulation and different vector length OSDM modulation by changing the vector length M.
The OSDM modulation technology has potential application prospects in underwater acoustic hybrid single-and multi-carrier modulations and adaptive modulation. For this reason, a pseudo-noise (PN) sequence-based underwater acoustic OSDM communication system (UWAOSDM/PN) is proposed in this paper. A data frame structure with PN sequence is designed to solve the multipath and Doppler effect of underwater acoustic channel. On the basis of the PN sequence, a compressive sensing method based on the orthogonal matching pursuit (OMP) algorithm and the minimum mean square error (MMSE) algorithm is designed for channel estimation and equalization. On the basis of the system construction, the relationship among the OSDM vector length M, the OSDM subcarrier number N, and the underwater acoustic channel length is further studied for adaptive modulation of underwater acoustic communication.
The paper is organized as follows. In Section 2, the system design of UWAOSDM/PN is described. In Section 3, the key modules are introduced. In Section 4, the simulation results are introduced and discussed. In Section 5, the system parameters and performance are compared and in Section 6 the conclusions are summarized.

System Design
The system structure of the proposed UWAOSDM/PN is shown in Figure 1. Transmitter: First of all, the transmitted data are encoded by convolutional code. Then, QPSK mapping is used to modulate the coded data. After that, the data is modulated by OSDM method. Then, PN sequences are added to get data frame. Finally, the baseband OSDM signal is modulated into the frequency band by carrier modulation.
Receiver: When the maximum correlation value between the received signal and the local PN sequence reaches the threshold value, the effective system synchronization is obtained. Then, the Doppler estimate and compensation is done. After that, the baseband signal is obtained by carrier demodulation. Then, the channel estimation and equalization are done to remove the multipath effect. After that, OSDM demodulation is used to get the transmitted data. Finally, through the QPSK demodulation and convolutional code decoding, the received data is obtained. The OSDM modulation technology has potential application prospects in underwater acoustic hybrid single-and multi-carrier modulations and adaptive modulation. For this reason, a pseudonoise (PN) sequence-based underwater acoustic OSDM communication system (UWAOSDM/PN) is proposed in this paper. A data frame structure with PN sequence is designed to solve the multipath and Doppler effect of underwater acoustic channel. On the basis of the PN sequence, a compressive sensing method based on the orthogonal matching pursuit (OMP) algorithm and the minimum mean square error (MMSE) algorithm is designed for channel estimation and equalization. On the basis of the system construction, the relationship among the OSDM vector length M, the OSDM subcarrier number N, and the underwater acoustic channel length is further studied for adaptive modulation of underwater acoustic communication.
The paper is organized as follows. In Section 2, the system design of UWAOSDM/PN is described. In Section 3, the key modules are introduced. In Section 4, the simulation results are introduced and discussed. In Section 5, the system parameters and performance are compared and in Section 6 the conclusions are summarized.

System Design
The system structure of the proposed UWAOSDM/PN is shown in Figure 1. Transmitter: First of all, the transmitted data are encoded by convolutional code. Then, QPSK mapping is used to modulate the coded data. After that, the data is modulated by OSDM method. Then, PN sequences are added to get data frame. Finally, the baseband OSDM signal is modulated into the frequency band by carrier modulation.
Receiver: When the maximum correlation value between the received signal and the local PN sequence reaches the threshold value, the effective system synchronization is obtained. Then, the Doppler estimate and compensation is done. After that, the baseband signal is obtained by carrier demodulation. Then, the channel estimation and equalization are done to remove the multipath effect. After that, OSDM demodulation is used to get the transmitted data. Finally, through the QPSK demodulation and convolutional code decoding, the received data is obtained.

OSDM Theory
First, the data X, in the 1 × MN dimension vector, is grouped into N groups of 1 × M dimension vector as follows:

OSDM Theory
First, the data X, in the 1 × MN dimension vector, is grouped into N groups of 1 × M dimension vector as follows: As we all know, inverse discrete Fourier transform (IDFT) formula is: IDFT matrix is: where, Using the Kronecker product between formulae (1) and (4), we can get: Formula (5) can be expanded to The orthogonality between any two S k can be demonstrated [11]. By summing S k , we can get the modulated signal s which is the baseband OSDM signal. This is a good mathematical expression for the principle of OSDM modulation.
Signal modulated by Formulae (1-6) also can be modulated by grouping and inverse fast Fourier transform (IFFT). As shown in Figure 2, first, all data are divided into M groups in a matrix denoted as d.
Do IFFT by column We define permutation matrix: To do OSDM modulation, the above modulation process can be expressed as: is row-wise write operations, N-point IDFT and row-wise read operations, respectively, and d is the left matrix in Figure 2. As shown in Figure 3, the data of the same group of IFFT are dispersed on the time scale. It is obvious that the OSDM modulation separates the original IFFT-transformed symbols into the time domain by (M-1) bits information in the same vector. In this way, not only the advantage of orthogonal carrier in the frequency domain is utilized, but also the time-scale diversity in the time domain is used. It can resist the influence of multipath of a certain delay on the data in the same IFFT group. Thus, OSDM can be considered as a hybrid modulation method in the time-frequency domain. At the receiver, in the OSDM demodulation, only the inverse process of the modulation process needs to be performed, and then the modulated signal can be demodulated. The demodulation can be expressed as: where, r means the data after channel equalization.

Frame Structure with PN Sequence
Due to the particularity of underwater acoustic channel, the data frame structure needs to be designed carefully in order to resist the adverse effect of UWA channel multipath and Doppler effect. The design of data frames must have the functions of synchronization, channel estimation, and We define permutation matrix: To do OSDM modulation, the above modulation process can be expressed as: where, matrix P H N,M , I M ⊗ F H N and P N,M is row-wise write operations, N-point IDFT and row-wise read operations, respectively, and d is the left matrix in Figure 2.
As shown in Figure 3, the data of the same group of IFFT are dispersed on the time scale. It is obvious that the OSDM modulation separates the original IFFT-transformed symbols into the time domain by (M-1) bits information in the same vector. In this way, not only the advantage of orthogonal carrier in the frequency domain is utilized, but also the time-scale diversity in the time domain is used. It can resist the influence of multipath of a certain delay on the data in the same IFFT group. Thus, OSDM can be considered as a hybrid modulation method in the time-frequency domain.
Do IFFT by column We define permutation matrix: To do OSDM modulation, the above modulation process can be expressed as: , , is row-wise write operations, N-point IDFT and row-wise read operations, respectively, and d is the left matrix in Figure 2. As shown in Figure 3, the data of the same group of IFFT are dispersed on the time scale. It is obvious that the OSDM modulation separates the original IFFT-transformed symbols into the time domain by (M-1) bits information in the same vector. In this way, not only the advantage of orthogonal carrier in the frequency domain is utilized, but also the time-scale diversity in the time domain is used. It can resist the influence of multipath of a certain delay on the data in the same IFFT group. Thus, OSDM can be considered as a hybrid modulation method in the time-frequency domain.  At the receiver, in the OSDM demodulation, only the inverse process of the modulation process needs to be performed, and then the modulated signal can be demodulated. The demodulation can be expressed as: where, r means the data after channel equalization.

Frame Structure with PN Sequence
Due to the particularity of underwater acoustic channel, the data frame structure needs to be designed carefully in order to resist the adverse effect of UWA channel multipath and Doppler effect. The design of data frames must have the functions of synchronization, channel estimation, and At the receiver, in the OSDM demodulation, only the inverse process of the modulation process needs to be performed, and then the modulated signal can be demodulated. The demodulation can be expressed as: where, r means the data after channel equalization.

Frame Structure with PN Sequence
Due to the particularity of underwater acoustic channel, the data frame structure needs to be designed carefully in order to resist the adverse effect of UWA channel multipath and Doppler effect. The design of data frames must have the functions of synchronization, channel estimation, and Doppler estimation. To solve the above problems, the frame structure with three pseudo-noise (PN) sequences is designed.
PN sequences are set at the head, middle, and tail of the data frame, respectively. Eight OSDM symbols with cyclic prefix are inserted among three PN sequences to form the data frame. The data frame format is shown in Figure 4.
The PN sequence consists of two identical short PN sequences to overcome multipath. PN sequences have three roles in this system. The header PN sequence can be used for signal detection and synchronization. At the same time, the header PN sequence and the middle PN sequence can be regarded as a block pilot, respectively, which is used as the pilot to do channel estimation for the subsequent OSDM symbols. Last but not least, the three PN sequences are used to do a Doppler estimation by the relative position of the three correlation peaks of sliding cross-correlation operation of the received signal.
The PN sequence not only has good randomness and cross-correlation which make it easy to be separated from other signals, but also the deterministic and repeatable nature of the PN sequence makes it have good anti-interference performance. More significantly, for the pseudo randomness, when the Toeplitz method is used to construct an observation matrix, the PN sequence can be used as initial vectors to create observation matrix φ. In this way, the matrix φ satisfies the important RIP metric. Doppler estimation. To solve the above problems, the frame structure with three pseudo-noise (PN) sequences is designed. PN sequences are set at the head, middle, and tail of the data frame, respectively. Eight OSDM symbols with cyclic prefix are inserted among three PN sequences to form the data frame. The data frame format is shown in Figure 4.
The PN sequence consists of two identical short PN sequences to overcome multipath. PN sequences have three roles in this system. The header PN sequence can be used for signal detection and synchronization. At the same time, the header PN sequence and the middle PN sequence can be regarded as a block pilot, respectively, which is used as the pilot to do channel estimation for the subsequent OSDM symbols. Last but not least, the three PN sequences are used to do a Doppler estimation by the relative position of the three correlation peaks of sliding cross-correlation operation of the received signal.
The PN sequence not only has good randomness and cross-correlation which make it easy to be separated from other signals, but also the deterministic and repeatable nature of the PN sequence makes it have good anti-interference performance. More significantly, for the pseudo randomness, when the Toeplitz method is used to construct an observation matrix, the PN sequence can be used as initial vectors to create observation matrix ϕ. In this way, the matrix ϕ satisfies the important RIP metric. Under this frame structure design, the system data rate can be calculated as: where, c r represents the coding rate. R means the coding rate. N is the number of symbol and D is the amount of each OSDM symbol, s T , cp T and PN T indicates the symbol duration, the cyclic prefix duration, and PN sequence duration, respectively.

Doppler Estimation and Compensation
Due to the water flows and the relative motion between the transmitters and receivers, the system needs to correct the Doppler effect on the received signal. The replicated PN sequences are placed at the beginning, middle, and end of the data frame, respectively. When the signal is affected by Doppler, the receiver signal length will be changed. Through the correlation calculation of the PN sequence at the receiver, the relative positions of the correlation peaks are detected to estimate the Doppler [20].
From the data frame structure, three different maximum values can be obtained in the search process, and the signal length influenced by Doppler can be estimated by the position of the first and the last maximum value. Since the PN sequence is affected by Doppler, the good autocorrelation property of the PN sequence is severely damaged. Therefore, the Doppler estimation should not use the cross-correlation calculation of the local PN sequence. On the contrary, using sliding autocorrelation between the received signal and the first PN sequence that intercepts from the received signal, three distinct correlation peaks can be obtained. Both the length of the signal at the receiver L  and the length of the signal at the transmitter L can be calculated by the relative position of correlation peaks, respectively.
Assume that the Doppler factor is α, by comparing L with L  , the receiver infers how the received signal has been compressed or extended by the UWA channel. Under this frame structure design, the system data rate can be calculated as: where, r c represents the coding rate. R means the coding rate. N is the number of symbol and D is the amount of each OSDM symbol, T s , T cp and T PN indicates the symbol duration, the cyclic prefix duration, and PN sequence duration, respectively.

Doppler Estimation and Compensation
Due to the water flows and the relative motion between the transmitters and receivers, the system needs to correct the Doppler effect on the received signal. The replicated PN sequences are placed at the beginning, middle, and end of the data frame, respectively. When the signal is affected by Doppler, the receiver signal length will be changed. Through the correlation calculation of the PN sequence at the receiver, the relative positions of the correlation peaks are detected to estimate the Doppler [20].
From the data frame structure, three different maximum values can be obtained in the search process, and the signal length influenced by Doppler can be estimated by the position of the first and the last maximum value. Since the PN sequence is affected by Doppler, the good autocorrelation property of the PN sequence is severely damaged. Therefore, the Doppler estimation should not use the cross-correlation calculation of the local PN sequence. On the contrary, using sliding autocorrelation between the received signal and the first PN sequence that intercepts from the received signal, three distinct correlation peaks can be obtained. Both the length of the signal at the receiver L and the length of the signal at the transmitter L can be calculated by the relative position of correlation peaks, respectively.
Assume that the Doppler factor is α, by comparing L with L, the receiver infers how the received signal has been compressed or extended by the UWA channel. The receiver uses the estimated Doppler factor α to resample the received signal, so that the compression or extension caused by the Doppler can be compensated.

Channel Estimation and Compensation
In recent years, many studies have shown that UWA channels have sparse characteristics in many situations, that is to say, most channel taps are zero or very small and the number of channel taps is far less than the maximum channel delay. Due to the sparse characteristics, using the compressive sensing (CS) method to do channel estimation is very suitable for an underwater acoustic channel. On the one hand, the minimum mean square error (MMSE) method which has superior equalization performance should be adopted. The MMSE takes noise effects into account. As compared with the zero forcing (ZF) method, the MMSE can improve system performance. Also, the MMSE method can get higher diversity gain than the ZF method.
In the UWAOSDM/PN system, channel estimation is conducted in the time domain and equalization is calculated in the frequency domain. We assume that the UWA channel is sufficiently stable to allow prediction. After synchronizing the received signal using Doppler estimation and compensation, channel estimation is conducted according to the above principle. Figure 5 shows the channel estimation and equalization process. First, the PN sequence is used as the initial vector and the Toeplitz format is employed to construct the observation matrix. Because of the pseudo randomness of the PN sequence, the constructed matrix meets the RIP principle [21]. Then, the orthogonal matching pursuit (OMP) method is used as the recovery method to estimate the sparse underwater acoustic channel in the time domain, and then, transforming the time domain h to the frequency domain by FFT. Last, calculating the channel equalization coefficients C n according to the MMSE criterion.
The receiver uses the estimated Doppler factor α to resample the received signal, so that the compression or extension caused by the Doppler can be compensated.

Channel Estimation and Compensation
In recent years, many studies have shown that UWA channels have sparse characteristics in many situations, that is to say, most channel taps are zero or very small and the number of channel taps is far less than the maximum channel delay. Due to the sparse characteristics, using the compressive sensing (CS) method to do channel estimation is very suitable for an underwater acoustic channel. On the one hand, the minimum mean square error (MMSE) method which has superior equalization performance should be adopted. The MMSE takes noise effects into account. As compared with the zero forcing (ZF) method, the MMSE can improve system performance. Also, the MMSE method can get higher diversity gain than the ZF method.
In the UWAOSDM/PN system, channel estimation is conducted in the time domain and equalization is calculated in the frequency domain. We assume that the UWA channel is sufficiently stable to allow prediction. After synchronizing the received signal using Doppler estimation and compensation, channel estimation is conducted according to the above principle. Figure 5 shows the channel estimation and equalization process. First, the PN sequence is used as the initial vector and the Toeplitz format is employed to construct the observation matrix. Because of the pseudo randomness of the PN sequence, the constructed matrix meets the RIP principle [21]. Then, the orthogonal matching pursuit (OMP) method is used as the recovery method to estimate the sparse underwater acoustic channel in the time domain, and then, transforming the time domain h to the frequency domain by FFT. Last, calculating the channel equalization coefficients n C according to the MMSE criterion.
Serial-to-parallel FFT IFFT Parallel-to-serial

Orthogonal Matching Pursuit Channel Estimation
The UWA channel is a time domain sparse channel whose length is L and sparsity is K, so we can use the OMP algorithm of compressive sensing to do channel estimation.
After adding the cyclic prefix, the linear convolution will be changed to circular convolution, and the circular convolution can be expressed as:

Orthogonal Matching Pursuit Channel Estimation
The UWA channel is a time domain sparse channel whose length is L and sparsity is K, so we can use the OMP algorithm of compressive sensing to do channel estimation.
After adding the cyclic prefix, the linear convolution will be changed to circular convolution, and the circular convolution can be expressed as: . . .

Minimum Mean Square Error Channel Equalization
The minimum mean square error function can be expressed as: Derivation of Equation (14) and let the partial derivative be 0 and the equalization coefficient is shown as: where, H n is estimated by OMP that have been convert into frequency domain, δ means the channel noise.

Simulation Conditions
In this section, simulation results are provided to illustrate the performance of the proposed UWAOSDM/PN system over UWA channels. In the simulation, the duration T s of each OSDM symbol is 213 ms, the data bits D of each OSDM symbol is a constant 1024, that is to say, the product of M and N is a fixed value D. Table 1 shows the simulation parameters of the UWAOSDM/PN system. As shown in Equation (16), the time-domain model of the UWA channel is composed of the attenuating component and the time-delay component. Assuming that the channel response is wide-sense stationary within one data frame.
where, p represents the number of multipaths, A p denotes the attenuation coefficient of each path, δ means impulse response function and τ p stands the delay time of the pth path. Figure 6 shows two UWA channel models adopted in this paper. Both models have been normalized. As shown in the figure, the horizontal axis is delay time, the vertical axis is normalized amplitude, and Figure 6a shows the channel model 1 of BELLHOP [22] which contains five paths. The first path is direct path. The maximum delay spread is about 6.9 ms. Figure 6b shows the channel model 2 which contains seven paths [23]. The first path is a direct path. Its maximum delay spread is about 12.5 ms.  (16) where, p represents the number of multipaths, p A denotes the attenuation coefficient of each path, δ means impulse response function and p τ stands the delay time of the pth path. Figure 6 shows two UWA channel models adopted in this paper. Both models have been normalized. As shown in the figure, the horizontal axis is delay time, the vertical axis is normalized amplitude, and Figure 6a shows the channel model 1 of BELLHOP [22] which contains five paths. The first path is direct path. The maximum delay spread is about 6.9 ms. Figure 6b shows the channel model 2 which contains seven paths [23]. The first path is a direct path. Its maximum delay spread is about 12.5 ms.  Figure 7 shows the performance comparison of the UWAOSDM/PN system with different vector lengths under different UWA channels. According to the channel model adopted in this paper, the maximum delay time of the channel model 1 is 6.9 ms and that of the model 2 is 12.5 ms. The system sampling rate is 96 kHz and the up-sample coefficient is 16. With these parameters, we can calculate the vector duration of different length vectors, such that the vector duration is 5.34 ms when the vector length is 32, the vector duration is 10.67 ms when the vector length is 64, and the vector duration is 21.33 ms when the length is 128. The simulation results in Figures 7a and b show that the system bit error rate (BER) of OSDM is lower than that of OFDM when the vector length M of OSDM system is larger than the maximum channel delay spread. With the increase of the vector length, the system can obtain more diversity gain. Simultaneously, the simulation results in Figures 7a and b show that the performance of a single-carrier system is also lower than that of the OSDM system. Therefore, the increase of vector length is limited in a scope. The vector length can be expressed as:  Figure 7 shows the performance comparison of the UWAOSDM/PN system with different vector lengths under different UWA channels. According to the channel model adopted in this paper, the maximum delay time of the channel model 1 is 6.9 ms and that of the model 2 is 12.5 ms. The system sampling rate is 96 kHz and the up-sample coefficient is 16. With these parameters, we can calculate the vector duration of different length vectors, such that the vector duration is 5.34 ms when the vector length is 32, the vector duration is 10.67 ms when the vector length is 64, and the vector duration is 21.33 ms when the length is 128. The simulation results in Figure 7a,b show that the system bit error rate (BER) of OSDM is lower than that of OFDM when the vector length M of OSDM system is larger than the maximum channel delay spread. With the increase of the vector length, the system can obtain more diversity gain. Simultaneously, the simulation results in Figure 7a,b show that the performance of a single-carrier system is also lower than that of the OSDM system. Therefore, the increase of vector length is limited in a scope. The vector length can be expressed as:

System Parameters Analysis of UWAOSDM/PN
where, D represents channel maximum delay spread, S is system sampling frequency, λ is up-sample coefficient and means ceil the number.
According to Equation (17), the vector length of the OSDM system should be larger than the maximum channel delay and keep as many subcarriers as possible to ensure the system carrier flexibility.
where, D represents channel maximum delay spread, S is system sampling frequency, λ is up-sample coefficient and   means ceil the number.
According to Equation (17), the vector length of the OSDM system should be larger than the maximum channel delay and keep as many subcarriers as possible to ensure the system carrier flexibility.

System Peak-to-Average Power Ratio
The peak-to-average power ratio is defined as:  In addition, the PAPR reflects the advantage of the OSDM from another perspective. High PAPR has low signal energy utilization and short propagation range, while low PAPR has high signal energy utilization and long propagation range. Therefore, the relationship between the vector length

System Peak-to-Average Power Ratio
The peak-to-average power ratio is defined as: Figure 8 shows the complementary cumulative distribution function (CCDF) curves of the PAPR of the UWAOSDM/PN system with different vector length M. As shown in the figure, the peak-to-average power ratio (PAPR) decreases with the increase of system vector length.
where, D represents channel maximum delay spread, S is system sampling frequency, λ is up-sample coefficient and   means ceil the number.
According to Equation (17), the vector length of the OSDM system should be larger than the maximum channel delay and keep as many subcarriers as possible to ensure the system carrier flexibility.

System Peak-to-Average Power Ratio
The peak-to-average power ratio is defined as:  In addition, the PAPR reflects the advantage of the OSDM from another perspective. High PAPR has low signal energy utilization and short propagation range, while low PAPR has high signal energy utilization and long propagation range. Therefore, the relationship between the vector length In addition, the PAPR reflects the advantage of the OSDM from another perspective. High PAPR has low signal energy utilization and short propagation range, while low PAPR has high signal energy utilization and long propagation range. Therefore, the relationship between the vector length and propagation range can be shown in Figure 9. The appropriate vector length can be selected according to propagation range, that is to say, OSDM is a PAPR tunable modulation method. and propagation range can be shown in Figure 9. The appropriate vector length can be selected according to propagation range, that is to say, OSDM is a PAPR tunable modulation method.

System Power of ICI
The UWA multi-carrier system is more susceptible to Doppler frequency offset. The OSDM system signal is similar to a traditional OFDM [24]. Therefore, the power of ICI caused by the Doppler spread can be bounded by: where, 1 The formula (19) can be expressed as: The relationship between ICI P and relative velocity of transmitter and receiver is shown in Figure 10, low ICI P , is achieved by reducing the subcarrier number N. However, the spectrum flexibility decreases with the subcarrier number decreasing.

System Power of ICI
The UWA multi-carrier system is more susceptible to Doppler frequency offset. The OSDM system signal is similar to a traditional OFDM [24]. Therefore, the power of ICI caused by the Doppler spread can be bounded by: where, (19) can be expressed as: The relationship between P ICI and relative velocity of transmitter and receiver is shown in Figure 10, low P ICI , is achieved by reducing the subcarrier number N. However, the spectrum flexibility decreases with the subcarrier number decreasing. and propagation range can be shown in Figure 9. The appropriate vector length can be selected according to propagation range, that is to say, OSDM is a PAPR tunable modulation method.

System Power of ICI
The UWA multi-carrier system is more susceptible to Doppler frequency offset. The OSDM system signal is similar to a traditional OFDM [24]. Therefore, the power of ICI caused by the Doppler spread can be bounded by: where, 1 The formula (19) can be expressed as: The relationship between ICI P and relative velocity of transmitter and receiver is shown in Figure 10, low ICI P , is achieved by reducing the subcarrier number N. However, the spectrum flexibility decreases with the subcarrier number decreasing.  Figure 10. The relationship between P ICI and relative velocity of transmitter and receiver. Figure 11a shows the system performance with different equalization methods and different vector lengths. The diversity gain obtained by the ZF method is smaller than the MMSE method. With the increase of SNR, the diversity gain of the ZF method will decrease. On the contrary, the MMSE method can achieve good diversity gains with different vector lengths. Figure 11b shows the system performance with the MMSE equalization method and different channel estimation methods. Compared with least squares (LS) and DFT-based least squares (LS-DFT), the system using the OMP channel estimation method can obtain the lowest BER. It is worth noting that OMP and LS-DFT have similar BER performance. In the OMP method, the signal recovery and signal denoising is combined. LS-DFT needs to know the noise threshold to achieve better performance. vector lengths. The diversity gain obtained by the ZF method is smaller than the MMSE method. With the increase of SNR, the diversity gain of the ZF method will decrease. On the contrary, the MMSE method can achieve good diversity gains with different vector lengths. Figure 11b shows the system performance with the MMSE equalization method and different channel estimation methods. Compared with least squares (LS) and DFT-based least squares (LS-DFT), the system using the OMP channel estimation method can obtain the lowest BER. It is worth noting that OMP and LS-DFT have similar BER performance. In the OMP method, the signal recovery and signal denoising is combined. LS-DFT needs to know the noise threshold to achieve better performance.

System Complexity
For the SC-FDE, OFDM, and OSDM systems, system complexity is mainly determined by IFFT and FFT points. The modulation complexity of SC-FDE, OFDM, and OSDM are compared in Table 2. At the transmitter, D-point IFFT which complexity is O(DlogD) is used to modulate the OFDM signal, and M times N-point IFFT which complexity is O (MNlogN) is used for the OSDM modulation. Given that the product of M and N is a fixed value D, the OSDM and OFDM modulation complexity is merely influenced by logarithmic term. With the increase of vector length M (the decrease of subcarriers N), the modulation complexity of the OSDM system decreases with fewer subcarriers. Generally speaking, the transmitter complexity is sorted from low to high as SC-FDE < OSDM < OFDM (21) At the receiver, block equalization is considered in this paper. OFDM equalization and demodulation can be implement at the same time, which complexity is O(DlogD + DLK). But both OSDM and SC-FDE are needed to be transformed by IFFT and FFT for equalization which complexity is O(2DlogD + DLK). Besides that, M times N-point FFT which complexity is O(MNlogN) is needed for OSDM demodulation. The receiver complexity is sorted from low to high as OFDM < SC-FDE < OSDM (22)

System Complexity
For the SC-FDE, OFDM, and OSDM systems, system complexity is mainly determined by IFFT and FFT points. The modulation complexity of SC-FDE, OFDM, and OSDM are compared in Table 2. At the transmitter, D-point IFFT which complexity is O(DlogD) is used to modulate the OFDM signal, and M times N-point IFFT which complexity is O (MNlogN) is used for the OSDM modulation. Given that the product of M and N is a fixed value D, the OSDM and OFDM modulation complexity is merely influenced by logarithmic term. With the increase of vector length M (the decrease of subcarriers N), the modulation complexity of the OSDM system decreases with fewer subcarriers. Generally speaking, the transmitter complexity is sorted from low to high as SC-FDE < OSDM < OFDM (21)

Modulation Scheme Transmitter Complexity Receiver Complexity
At the receiver, block equalization is considered in this paper. OFDM equalization and demodulation can be implement at the same time, which complexity is O(DlogD + DLK). But both OSDM and SC-FDE are needed to be transformed by IFFT and FFT for equalization which complexity is O(2DlogD + DLK). Besides that, M times N-point FFT which complexity is O(MNlogN) is needed for OSDM demodulation. The receiver complexity is sorted from low to high as OFDM < SC-FDE < OSDM Table 3 shows a performance comparison of the proposed UWAOSDM/PN system with the other UWA OSDM systems. The comparison is based on simulation results. As shown in the comparison, the UWAOSDM/PN has a relatively high communication rate, a lower system BER and, a low receiver complexity. The effectiveness of the proposed PN sequence-based compressive sensing channel estimation and equalization with the OMP and MMSE algorithm is demonstrated. The main contributions of this paper are as follows: (1) Provide a solution of combining PN sequence, OMP, and MMSE for underwater acoustic OSDM communication. Compared with the other UWA OSDM systems, the proposed UWAOSDM/PN has a relatively high communication rate, a lower system BER, and a low receiver complexity. (2) PN sequences are specially selected and designed for the system. PN sequences have three roles in this system. The header PN sequence can be used for signal detection and synchronization. The header PN sequence and the middle PN sequence can be regarded as a block pilot, respectively, which is used as the pilot to do channel estimation for the subsequent OSDM symbols. In particular, PN sequence is suitable for the OMP algorithm. The Toeplitz method is used to construct an observation matrix, and the PN sequence can be used as initial vectors to create observation matrix φ for the pseudo randomness. In this way, the matrix φ satisfies the important RIP metric. The three PN sequences are used for the Doppler estimation by the relative position of the three correlation peaks of sliding cross-correlation operation of the received signal.

Comparison with Related Systems
(3) For the potential application prospects in the underwater acoustic hybrid single-and multi-carrier modulations and adaptive modulation, the relationship among the OSDM vector length M, the OSDM subcarrier number N, and the underwater acoustic channel length is further studied.

Conclusions
In this paper, an UWA communication system based on OSDM technology is proposed and designed, and a special data frame format based on PN sequence is designed. The OSDM signals with different vector lengths are analyzed. For the problems of long delay, strong multipath, and large-scale Doppler frequency shift in underwater channels, OMP-based compressive sensing channel estimation, MMSE equalization, sliding correlation Doppler estimation and resampling compensation methods are adopted. As shown in the simulation results, compared with a traditional OFDM system, the proposed UWAOSDM/PN system has low PAPR and ICI When the OSDM vector length M is bigger than the channel length, the system subcarrier flexibility is guaranteed, and the bit error rate of the OSDM system is lower than that of the OFDM system and the single-carrier system. The PN sequence-based compressive sensing channel estimation and equalization using the OMP and MMSE algorithms has a good performance to resist the multipath effect of underwater acoustic channel.
The OSDM modulation technology has potential application prospects in underwater acoustic hybrid single-and multi-carrier modulations and adaptive modulation. In the future, on the basis of the UWA OSDM system studied in this paper, an adaptive multimode UWA communication technology study based on the OSDM modulation will be carried out.