M-ary nonlinear sine chirp spread spectrum for underwater acoustic communication based on virtual time-reversal mirror method

Linear chirp spread spectrum technique is widely used in underwater acoustic communication because of their resilience to high multipath and Doppler shift. Linear frequency modulated signal requires a high spreading factor to nearly reach orthogonality between two pairs of signals. On the other hand, nonlinear chirp spread spectrum signals can provide orthogonality at a low spreading factor. As a result, it improves spectral efficiency and is more insensitive to Doppler spread than the linear counterpart. To achieve a higher data rate, we propose two variants (half cycle sine and full cycle sine) of the M-ary nonlinear sine chirp spread spectrum technique based on virtual time-reversal mirror (VTRM). The proposed scheme uses different frequency bands to transmit chirp, and VTRM is used to improve the bit error rate due to high multipath. Its superior Doppler sensitivity makes it suitable for underwater acoustic communication. Furthermore, the proposed method uses a simple, low-power bank of matched filters; thus, it reduces the overall system complexity. Simulations are performed in different underwater acoustic channels to verify the robustness of the proposed scheme.

used matched filter to detect chirp signals. However, the linear chirp signal's main problem is that it requires high bandwidth-time (BT) to achieve orthogonality between two chirp signals. As a result, it requires a large bandwidth, which is one of the valuable resources for any communication system, especially in UWA communication, because bandwidth is limited. To increases the data rate of the underwater communication system and to solve the problem of high SF, spectral efficiency, and inter-symbol interference caused by underwater acoustic channels, M-ary Sine Chirp Spread Spectrum (MS-CSS) based on VTRM is proposed in this paper. MS-CSS is a combination of two techniques used in UWA communication. The nonlinear sine upchirp and downchirp is used to represent data bits, and M-ary technique is utilized to achieve a higher bit rate by transmitting at different frequency bands, whereas correlators are used to detect the symbols. In contrast, VTRM implementation helps improve the bit error rate caused by the multipath environment.
The rest of the paper is organized as follows: Sect. 2 introduces the analysis of linear and nonlinear chirp signals. Section 3 presents the chirp spread spectrum concept using the M-ary technique, whereas the system model of the M-ary sine chirp spread spectrum is given in Sect. 4. The channel estimation using matching pursuit and virtual time-reversal mirror methods is provided in Sect. 5 and Sect. 6, respectively. Section 7 validates the proposed scheme's simulation results and Sect. 8 finally concludes the research article.

Linear chirp signal
In traditional chirp, signal frequency increases linearly (upchirp) and decreases (downchirp) over a given time period. The LFM signal can be expressed as where A is the amplitude of the chirp, f o is the starting frequency, T LC is the chirp duration, and µ is the frequency sweep rate and can be defined as µ = B/T LC where B is the bandwidth and can be defined as B = f final − f initial . When the bit "0" is to be sent, the function c 0 (t) that is upchirp is transmitted, and when bit "1" is to be sent, function c 1 (t) downchirp is transmitted. The matched filter-based receiver is used for detection. The expression for the autocorrelation and cross-correlation in [28] and it can be given as, where C(x) and S(x) are defined as the Fresnel function. Figure 1 shows that the output function of a matched filter has most of its energy existing from −1/B ≤ τ ≤ 1/B of period, and the peak of the matched filter is amplified up to √ T LC B . This pulse compression property of LFM makes it resilient to multipath and Doppler shifts. The cross-correlation coefficient of LFM is defined by Figure 2 shows that the cross-correlation coefficient ρ l depends upon the product of time period of linear chirp T LC , and the bandwidth B of the chirp signal. As T LC B increases, ρ l decreases, but it never reaches 0, even if the T LC B is equal to 100. Therefore, linear chirps are not entirely orthogonal, even at a high T LC B value and can be considered quasi-orthogonal [29]. Fig. 1 Linear frequency-modulated correlator output

Nonlinear chirp signal analysis
There are many candidates of nonlinear chirps as a carrier, but one of the preferable contenders for UWA communication can be sine chirp because of its orthogonality between the pair of chirp under certain conditions. The sine chirp can be mathematically written as where T S is the duration of the sine chirp, and η represents the number of cycles. For example, if η is set to "1" in Eq. (5), half-cycle sine chirp is obtained, whereas if η is set to "2" in Eq. (5), full-cycle sine chirp is obtained. The autocorrelation property of the sine chirp is similar to the linear chirp signal and can be expressed using the equation.
where θ = πη 2T s , k m s 1 = 2BT s η cos πη 2T s τ , k m s 2 = 2BT s η sin πη 2T s τ and J 2n is the first kind Bessel function. In Eq. (6), the first and third terms have a negligible effect on the autocorrelation coefficient. Hence Eq. (6) is written as The following equation gives the cross-correlation coefficient of pair of sine chirp where θ = πη 2T s , k m s 1 = 2BT s η cos πη 2T s τ , k m s 2 = 2BT s η sin πη 2T s τ and J n is the first kind Bessel function. The cross-correlation coefficient at time shift τ = 0 can be simplified as in [30] and given by the equation below. Figure 4 shows that the cross-correlation coefficient and bandwidth time product of nonlinear sine chirps (half and full cycle) are periodically zero; hence, they are entirely orthogonal at the specific condition. For example, the half-sine period has a cross-correlation coefficient zero at various periodic bandwidth-time product points, including 4.4 sHz, 9 sHz, and 20 sHz, while the full sine period has zero cross-correlation coefficient at bandwidth-time 5.5 sHz, 8.65 sHz, and 18 sHz. In contrast, the LFM signal never reaches zero. Figure 3a shows that the autocorrelation between the pair of one complete sine cycle and the half-sine cycle is similar to the linear chirp. Still, side lobes of both the nonlinear chirps are greater than linear counterparts, which can be reduced using any side lobe reduction scheme. Figure 3b verifies the relationship between the cross-correlation coefficient and bandwidth time product of nonlinear chirps. The crosscorrelation between the pair to full and half-period sine chirp at 18sHz and 20sHz in Fig. 3b shows that at time shift ( τ = 0) the value of cross-correlation is zero for both cycles. However, the full period sine chirp has a large value at non zero time shift, resulting in false detection in the UWA dispersive channel. As a result, the overall performance of the system can be affected. Similar results can be obtained for BT's different values; as the value of BT is decreased, the side lobe of auto-correlation increases for both linear and nonlinear chirp signals (Fig. 4).

M-ary chirp spread spectrum
M-ary chirp spread spectrum can be implemented by transmitting upchirp and downchirp in two different frequency bands, as shown in Fig. 5. The main issue with linear chirp is that it requires high BT products to achieve acceptable performance. As the number of chirps are increased to increase the data rate, the system's overall performance is worsening. The sine chirp is used in this paper, for which the pair of chirp is orthogonal at low SF. In Fig. 6, the time-frequency graph of 8-ary sine chirp at four different frequency bands is shown. When the symbol to be transmitted is "000", the sine chirp in the first frequency band with a positive gradient is transmitted similarly. If "001" is transmitted, the frequency band is the same but having a negative gradient. As a result, there is a 180 phase shift in sine chirp, as shown in Fig. 6. For the next two symbols, the frequency band ranges from f 2 to, f 3 and upchirp and downchirp is transmitted

System model of MS-CSS
The block diagram of the M-ary sine chirp spread spectrum is shown in Fig. 7. First, the stream of serial data bits x(n) is converted from serial to parallel, then chirp is generated by the chirp generator according to the given symbol; for example, if the symbol to be transmitted is "000", then upchirp c s0 (t) is generated in the frequency band 1 as shown in Fig. 7, then it is pass through the hydrophone for transmitting it through an underwater acoustic channel in the presence of Additive White Gaussian Noise (AWGN). The transmit signal is received by the transducer and passed through the bank of matched filters. The matched filter works on the principle of correlation. The received signal y(t) is multiplied with a time-shifted version of the initially generated chirp c s0 (t − τ )(in this case) and integrated over the chirp period T s . The matched signal gives maximum output while all other unmatched output has low values. As a result, the maximum value detector determines the given chirp signal, and the corresponding symbol is mapped against it. Finally, the output data is received by converting the given symbols from parallel to serial.

Channel estimation
The channel estimation is one of the requisite components of VTRM-based systems. Its performance is highly dependent on how accurately the channel is estimated. The sparse nature of the UWA channel results in a considerable delay spread with minimal nonzero values. Compressed sensing algorithms are widely used to solve the problem of sparsity [31]. In many underwater acoustic communication systems, compressed sensing algorithms are used. In [32], a compressed sensing algorithm and equalization are used to evaluate orthogonal signal division multiplexing performance. In contrast, matching pursuit helps resolve the problem of sparse channel estimation. It is widely used due to its robustness and computational efficiency [33]. Matching pursuit along with VTRM used in [26,34], 35 proves to be very useful. The channel is estimated by sending a training signal through a sparse UWA channel. The training signal received can be given by the following equation Let y(n) is the received training signal and l its length, where the number of channel taps. The above equation can be rewritten as (10) It can be concluded from the above equation that most of the channel taps are either zero or have minimal values.
The iteration process estimates the channel until the given threshold is achieved. Initially, the received signal Y is considered as a residual ϒ 0 and ϒ 0 is appropriately matched with the column of X , which is denoted by x 1 . If the threshold is not achieved, the irritation continues. At m th iteration, the residual ϒ m−1 for which the X has the maximal rank-one projection can be expressed as x m and it is given by the following equation X k contains the column of X and the superscript H is Hermitian Transpose. The value of ∧ h m which can be defined as the components h associated with X k can be computed as Finally, the residual vector can be calculated as

Virtual time-reversal mirror method
The virtual time-reversal mirror method is widely used for equalization in many types of communication channels, including underwater acoustic communication channel [36]. VTRM-based equalization had proven to be effective in underwater acoustic single . . .  [37] proposed orthogonal chirp carrier, and to mitigate the effect of inters symbol interference VTRM is used. Moreover, VTRM proves to be useful in multiuser chirp spread spectrum communication proposed in [38]. Furthermore, VTRM is also used in time-varying underwater acoustic communication channel [39,40]. The idea behind the working principle of the VTRM technique is to estimate the channel impulse response by processing the received training signal. The authors in [41] showed that the performance of VTRM is dependent on the use of appropriate training signals and the duration of the time window. The VTRM system design is shown in Fig. 8; first, the estimated channel impulse response h ′ (t) is time-reversed h ′ (−t) , and then it is convolved with the received signal [42]. Let x r (t) be the transmit signal passed through the UWA channel can be given as Then after virtual time reversal, the received signal can be expressed as The above equation h(t) ⊗ h ′ (−t) is VTRM based on equalized channel. It is a concurrence between the real channel impulse response and estimated channel impulse response. If the channel is estimated precisely, the direct path's energy concentration is much greater than the other reflected paths in the matched channel. As a result, it not only reduces inter-symbol interference but also improves the SNR.

Simulation using BELL HOP ray tracing algorithm
The parameters used for the simulation are shown in Table 1. The simulation setup of the UWA channel generated using the BELLHOP ray tracing algorithm is shown in Fig. 9. The actual impulse response of a channel and estimated channel response when the distance between transmitter and receiver is 1 km is shown in Fig. 10. The simulation results in Fig. 11 depict that when the chirp period is 10 ms and Spread factor is 9sHz for half cycle sine (HCS) chirp and LFM, whereas 8.65sHz for full cycle sine (FCS) chirp without channel estimation and VTRM, the BER of HCS chirp is zero after − 4 dBs. For the case of FCS chirp and LFM, the BER is zero after − 2 dB. The UWA channel exhibits high multipath. Consequently, the system's overall performance is affected, but using VTRM improves the bit error rate. As shown in Fig. 11, using the VTRM method, all three types of Chirps' performance is improved while HCS Chirp outperforms the other two types. To investigate further, the BT was decreased from 9sHz to 4.4sHz for HCS chirp and LFM, while for FCS chirp, it was decreased from 8.65sHz to 5.5sHz as at the mentioned values, nonlinear chirps are orthogonal. The simulation   Underwater acoustic actual and estimated channel impulse response at a distance of 1 km using BELLHOP a actual impulse response, b estimated impulse response using matching pursuit results shown in Fig. 12 indicate that as the SF is decreased, linear chirp performance deteriorates compared to nonlinear chirp. To increase the bit rate, the chirp period was decreased from 10 to 5 ms. The result shown in Fig. 13 indicates that the performance of the nonlinear sine full cycle is better than the other two types. In contrast, at low SF the full cycle and half-cycle have similar performance but outperform the linear chirp, as illustrated in Fig. 14. The simulation results validate that the MS-CSS solves the problem of increasing the information rate at low spreading factor and improving the system's overall performance influenced by the high multipath due to underwater acoustic channel by using the VTRM technique.

Simulation using watermark reply channel
The watermark replay channel is a time-varying UWA channel; it is widely used to evaluate the performance of different types of the modulation scheme of UWA communication. The watermark channel uses the equation below to evaluate the performance of the system.
where y(t) is the output of the time-varying UWA channel recorded impulse response ∧ h (t, τ ) of the real channel h (t, τ ) , while x(t) is the input signal and n(t) is the noise. The watermark offers many different types of channels. The channel used in this simulation is NOF1 (Norway Oslofjord). The details of the simulation parameters are shown in (19) Table 2. It is to be noted that both the transmitter and receiver are static and time-varying Doppler is compensated using resampling. The impulse response of the channel is shown in Fig. 15. The performance of the HCS and LFM chirp was evaluated at BT 9sHz and 4.4sHz, while the performance of FCS at 8.65 sHz and 5.5 sHz. The results in Fig. 16 shows that the performance of both the nonlinear chirps are better than the linear chirp, but NOF1 channel requires higher SNR to achieve satisfactory performance compared to the BELLHOP channel. In addition, VTRM is used to improve the performance of the system. Moreover, the performance of the linear and nonlinear chirp is also evaluated at 5 ms and SF 9 sHz for HCS and LFM and 8.65 sHz for FCS in NOF1 channel. The results show that linear chirp performance is acceptable after 18 dBs at BT = 9sHz and improves after VTRM processing while the performance of HCS and FCS at 5 ms remains superior then LFM and further improves after using VTRM as shown in Fig. 17a. Furthermore, when the BT of the LFM and HCS was decreased to 4.4sHz, and for FCS 5.5sHz, the LFM showed disappointing results, on the other hand performance of both the nonlinear chirp also degraded, however after VTRM processing the nonlinear chirp showed satisfactory results, while HCS was better than FCS as shown in Fig. 17b.

Conclusion and future work
In this paper, MS-CSS-VTRM is proposed. The mentioned technique uses nonlinear Sine chirp, which can perform better at low spreading factor, which plays a significant role in hydro-acoustic communication's overall performance. However, to increase the information rate, M-ary technique is presented by transmitting the chirp at different frequency bands. At the same time, VTRM equalizes the UWA channel to minimize the effect of multipath interference. As a result, the overall reliability of the system is improved. Firstly, the paper presented the sine chirp and proved that at low SF it is orthogonal. Secondly, the system model of the M-ary sine chirp spread spectrum is presented. Then, the concept of VTRM in underwater acoustic communication to mitigate inter-symbol interference is presented. Finally, the MS-CSS scheme proposed in this a b