Next Article in Journal
A Review of Thin-Film Magnetoelastic Materials for Magnetoelectric Applications
Next Article in Special Issue
The Hierarchic Treatment of Marine Ecological Information from Spatial Networks of Benthic Platforms
Previous Article in Journal
High Precision Low-Speed Control for Permanent Magnet Synchronous Motor
Previous Article in Special Issue
Event-Driven Data Gathering in Pure Asynchronous Multi-Hop Underwater Acoustic Sensor Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multiuser Chirp Spread Spectrum Transmission in an Underwater Acoustic Channel Applied to an AUV Fleet †

SEACom Department, L@bisen Yncréa-Ouest, 29228 Brest CEDEX 2, France
*
Author to whom correspondence should be addressed.
This paper is extension version of the conference paper: Bernard, C.; Bouvet, P.J. Multiuser Underwater Acoustic Communication for an AUV Fleet. In Proceedings of the OCEANS 2019—Marseille, Marseille, France, 17–20 June 2019.
Sensors 2020, 20(5), 1527; https://doi.org/10.3390/s20051527
Submission received: 14 February 2020 / Revised: 3 March 2020 / Accepted: 6 March 2020 / Published: 10 March 2020
(This article belongs to the Special Issue Underwater Sensor Networks)

Abstract

:
The objective of this paper is to provide a multiuser transmission technique for underwater acoustic communication in the framework of an Autonomous Underwater Vehicle (AUV) fleet. By using a variant of a Hyperbolically Frequency-Modulated (HFM) signal, we describe a new family of transmission techniques called MultiUser Chirp Spread Spectrum (MU-CSS), which allows a very simple matched-filter-based decoding. These techniques are expected to provide good resilience against multiuser interference while keeping good robustness to Underwater Acoustic (UWA) channel impairments like Doppler shift. Their implementation for the UWA scenario is described, and the performance results over a simulated shallow-water UWA channel are analyzed and compared against conventional Code-Division Multiple Access (CDMA) and Time-Division Multiple Access (TDMA) transmission. Finally, the feasibility and robustness of the proposed methods are verified over the underWater AcousTic channEl Replay benchMARK (Watermark), fed by several channel responses from sounding experiments performed in a lake.

1. Introduction

The Underwater Acoustic (UWA) channel is one of the most challenging channels for data communications. Due to the low celerity of acoustic waves (c = 1500 m· s 1 ), UWA channels are characterized by extensive multipath effects and large Doppler spreads. Moreover, frequency- dependent attenuation, temporal variations, and background noise limit the achievable data rate considerably [1,2]. On the other hand, Autonomous Underwater Vehicles (AUVs) are used for several marine applications, such as in the military field with anti-submarine warfare, in the science field with wreck exploration, or in the industrial field with offshore energy research. Nowadays, the concept of several AUVs working together within a fleet is an on-going research axis [3]. UWA communication with an AUV fleet is used to control vehicles (downlink) or to gather data from vehicles (uplink). The quality and reliability of communications are essential, mainly in shallow water areas for which the multipath effect is stronger, leading to extensive intersymbol interference.
Multiuser communication protocols in an UWA channel can be divided into two categories—random or deterministic protocols. In random protocols, the data rate cannot be predicted in advance due to the phenomenon of collisions between different users. A classical example of random protocol is ALOHA [4] and its variants [5], which use the long propagation delays to reduce the number of collisions and, consequently, to increase the data rate. Another example of random protocol is the Carrier Sense Multiple Access (CSMA) method [6], which is based on channel listening to avoid collisions. On the other side, deterministic protocols perform deterministic assignments of channel resources to the users so that their activity on the channel is predictable. The method we propose in this paper aims at building a new set of mutually orthogonal waveforms to be assigned to the users of an UWA channel so as to separate them easily at the receiver side. This consequently falls into the class of deterministic protocols. Traditional methods for deterministic multiuser UWA transmissions are inspired by radio communications and adapted to the UWA channel. As examples, we can cite the Time-Division Multiple Access (TDMA) [7], Frequency Divsion Multiple Access (FDMA) [8], Code-Division Multiple Access (CDMA) [9], and Multi-Carrier Code-Division Multiple Access (MC-CDMA) [10] transmissions. Typically, FDMA is considered inefficient, since the UWA channel has limited bandwidth and exhibits a large Doppler spread that requires guard frequency bands between users, leading to wasting of the data rate. MC-CDMA transmission schemes suffer from both time/frequency selectivity of the UWA channel and multiple-access interference, and require complex iterative equalizers. Consequently, in the following, we will focus only on the TDMA and CDMA strategies. TDMA allows several users to share the same frequency channel by dividing the signal into different time slots. Each user alternatively uses their own time slot to transmit data without interfering with other users. However, as the number of users increases, the waiting time per user increases and the user data rate decreases. In CDMA transmission, the different users transmit information data simultaneously through a different spreading sequence for each user. The disadvantage of this method lies in the multiuser interference provided by the non-orthogonality of spreading sequences, especially when the user communication channel is selective in time or in frequency. Moreover, this effect is increased when the interference power is much larger than the received signal power. This phenomenon is well known in mobile communication networks as the near–far problem. To cope with interference terms in CDMA, advanced equalization schemes can be invoked, such as multiuser detection [11] or the Multi-User Multiple-Input Multiple-Output (MU-MIMO) technique combined with Passive Phase Conjugation (PPC) [12], but at the price of a higher decoding complexity and a limited number of users. Recently, the authors of [13] proposed an alternative to CDMA and TDMA by using chirp waveforms for UWA multiuser communication. To reduce the multiuser interference, the Virtual Time Reversal Mirror (VTRM) technique is used with a Fractional Fourier Transform (FrFT) at the reception. However, this method requires an estimate of the different channels and is limited to four users because of interference.
In this paper, we describe a new transmission scheme based on Chirp Spread Spectrum (CSS), entitled MU-CSS, that we originally introduced in [14]. The basic idea consists of building a set of mutually orthogonal chirp-based waveforms that will be resistant to Doppler spread and Doppler shift. The objective is, on the one hand, to take benefit from the robustness of chirps against UWA channel impairments and, on the other hand, to use orthogonality to separate multiple users at the receiver side using a simple matched filter. With respect to [14], we derive three new methods for building MU-CSS that optimize mutual orthogonality between waveforms. By assuming an uplink scenario where a fleet of N u AUVs in motion needs to transmit data to a receiver situated at the sea surface, we provide a performance comparison of each method over simulated and experimental replay channels.
The paper is organized as follows: The system model and state of the art of multiuser transmissions are introduced in Section 2. The proposed MU-CSS multiuser schemes are presented in Section 3. Performance results of the proposed schemes against conventional multiuser transmissions are obtained in Section 4 by using a shallow-water UWA channel simulator derived from [15,16] and in Section 5 by using the Watermark replay channel [17] fed by experiments conducted in Ty-Colo lake, Saint-Renan, France. Finally, conclusions are drawn in Section 6.
In the following, | | . | | 2 denotes the euclidean norm, . the scalar product, E { . } denotes the statistical expectation, ( . ) * the complex conjugate, and u * v denotes the convolution product between u and v.

2. Multiuser Transmission

2.1. System Model

Transmitted Signal

Let d i , k be the k-th symbol transmitted by the i-th user; we assume that d i , k belong to a unit-amplitude Phase Shift Keying (PSK) alphabet, and are differentially encoded such that:
d i , k = d i , k 1 · b i , k with i [ 1 , N u ] , k [ 2 , N s ] ,
where b k is the original PSK data symbol and d i , 0 is set to 1. Beforehand, the data symbols b i , k are protected by a Forward Error Correction (FEC) code, followed by a random interleaver. In the following, the FEC code type will be a half-rate convolutive code with code generator ( 133 , 171 ) o . Moreover, N u denotes the number of users and N s the number of data symbols per frame. The choice of Differential Phase Shift Keying (DPSK) is motivated by the rapid fluctuation of the UWA channel and thus allows one to avoid the use of channel equalizers at the receiver side, which are sensitive to outdated channel estimations [18]. Thus, in an UWA communication channel with large delay spreads and rapid time variations, differential modulations are demonstrated to provide interesting performance and even outperform coherent modulation under certain conditions [19].
Let g i ( t ) be the transmit waveform associated to user i and T s the symbol duration; the baseband transmit signal for user i can be written as:
s i ( t ) = k = 1 N s d i , k g i ( t k T s ) .

2.2. Underwater Multiuser Channel

By assuming that users are mobile with relative motion v i , positive values of v i denote motion away from the receiver, while negative values denote motion toward the receiver; the received baseband signal is given by:
r ( t ) = i = 1 N u + h i ( τ , t ) s i ( 1 a i ) ( t τ ) e j 2 π f c a i ( t τ ) d τ + n ( t ) ,
where f c is the carrier frequency and a i = v i c is the Doppler scale factor. The UWA channel impulse response for the i-th user at time t is denoted by h i ( τ , t ) , and n ( t ) represents the additive noise, assumed to be Gaussian and zero-mean.

User Decoding

When the Doppler shift can be estimated at the receiver, the Doppler effect is usually removed prior to decoding by resampling the received baseband signal and compensating phase rotation as follows [1]:
z i ( t ) = r t 1 a i e j 2 π f c a i 1 a i t .
By assuming perfect time synchronization, information data of the i-th user can be estimated by matched filtering z i ( t ) with the transmitted waveform of user i, followed by integration over a symbol duration [20]:
d ^ i , k = max k T s 2 t ( k + 1 ) T s 2 + g i * ( u ) z i ( t u ) d u
= T s 2 T s 2 g i * ( t ) z i ( t + k T s ) d t
= γ i , k d i , k + η i , k + w i , k ,
where γ i , k denotes the bias of the decoder, η i , k the multiuser interference terms, and w i , k the additive noise terms; the exact expression of these three terms is provided in Appendix A.

2.3. Conventional Multiuser Transmission Schemes

2.3.1. CDMA

The objective of CDMA is to break up a finite transmission spectrum so that multiple users can access it at the same time. To accomplish time multiplexing, a code, chosen in a set of mutually orthogonal spreading codes, is assigned to each user [21]. For the i-th user, the transmitted waveform is expressed by:
g i ( t ) = c i ( t ) = l = 0 N S F 1 c i , l ϕ ( t l T c ) ,
where [ c i , 1 , c i , 2 , , c i , N S F ] is the spreading code of length N S F , T c is the chip duration, N S F is the spreading factor, and ϕ ( t ) is the pulse-shaping filter chosen as a Square Root Raised Cosine (SRRC) filter [20]. Since we are in an uplink scenario, the CDMA system is asynchronous, and spreading codes are chosen as Pseudo-Noise (PN) sequences generated pseudo-randomly such that their autocorrelation functions tend to Dirac functions as N S F grows, the mutual cross-correlation tends to zero.
At the receiver side, if T s > τ m a x , where τ m a x denotes the Root Mean Square (RMS) channel delay spread, and if the communication channel is constant over a symbol duration T s , the autocorrelation properties and quasi-orthogonality between users of PN codes lead the term η i , k in (7) to become negligible compared to γ i , k , and thus allow each user to be decoded separately [20].

2.3.2. TDMA

In a TDMA approach, the users are time-multiplexed, as depicted in Figure 1. The time slot assigned to one user is made of a frame slot of N s T s seconds, followed by a guard interval of duration T g so as to absorb multiuser interference. In order to deal with the frequency selectivity of the UWA channel, Direct Sequence Spread Spectrum (DSSS) signaling, with the same modulation parameters as CDMA, is chosen for each user such that the TDMA and CDMA approaches are equivalent in the single user scenario. The baseband received signal and the decoding process are given by particularizing (2) and (8) respectively with N u = 1 . One can note that a more spectrally efficient transmission scheme could be chosen for TDMA (see [22], for example) but at the price of higher complexity at the receiver side. Moreover, a higher spectral efficiency signaling scheme would make the comparison with CDMA difficult, especially in the single-user case.

3. MU-CSS Scheme

3.1. Generalities

Through the use of frequency-swept signals, which are resilient to the detrimental effects of the UWA channel, the CSS modulation technique offers robust performance with a very simple matched filtering-based decoder that makes such a communication scheme particularly well adapted to the UWA communication channel [23,24]. In the CSS system, a broad spectrum is occupied to modulate the information in order to achieve high processing gain and multipath resolution to the detriment of the spectral efficiency. In the following, we construct three multiuser schemes based on CSS signaling and, more precisely, on a Hyperbolically Frequency-Modulated (HFM) signal given by:
x ( t ) = cos ( 2 π ( k log ( 1 t t 0 ) + f l + f h 2 ) ) if T s 2 t T s 2 0 otherwise ,
where t 0 = T s ( f h + f l ) 2 ( f h f l ) , k = T s f l f h f h f l , f l f h , and T s is the duration of the HFM signal, whose instantaneous frequency is provided in Figure 2, with f h = B / 2 and f l = B / 2 , where B = 4 kHz and T s = 7.75 , 15.75 , and 31.75 ms.
The basic idea of MU-CSS consists of building an orthogonal basis of signals e i ( t ) thanks to the Gram–Schmidt process, where the waveform e i ( t ) is assigned to i-th user with i [ 1 , N u ] . The initial orthogonality between waveforms is brought by the combination of the HFM signals with orthogonal spreading sequences that are chosen as a Walsh–Hadamard code [21]. The set of spreading codes allows users to be differentiated at the receiver side, while the HFM waveform provides robustness against Doppler and delay spreads.

3.2. MU-CSS Gram–Schmidt Iterated

In this method, an iterative process is used to improve the mutual orthogonality between the chirp waveforms as well as the immunity against channel impairments.
Let e i ( l ) ( t ) denote the waveform corresponding to the i-th user with i { 1 , 2 , , N u } at iteration l { 1 , , N I T } . The process is based on the Gram–Schmidt method [25], as follows, for i > 0 :
e i ( l ) ( t ) = c i ( t ) + α i ( l ) e i 1 ( l ) ( t ) ,
where:
α i ( l ) = c i ( t ) , e i 1 ( l ) ( t ) | | e i 1 ( l ) ( t ) | | 2 2 = T s 2 T s 2 c i ( t ) e i 1 ( l ) * ( t ) d t | | e i 1 ( l ) ( t ) | | 2 2 ,
with c i ( t ) given by Equation (8). At the first iteration, we set e 0 ( 1 ) ( t ) = x ( t ) , where x ( t ) is defined in (9) and for l > 1 , i > 0 :
c i ( t ) = e i ( l 1 ) ( t ) .
The final waveform assigned to each user is obtained after N I T iterations of the abovementioned process by setting g i ( t ) = e i ( N I T ) ( t ) . The orthogonality between the different e i ( l ) ( t ) and the choice for the value of α i ( l ) are justified in Appendix B using the Gram–Schmidt procedure.

3.3. MU-CSS Gram–Schmidt Multiplication

In this method, the combination with the HFM is made by multiplying it with the spreading sequence while applying the Gram–Schmidt iteration process to ensure orthogonality. We start from:
e i ( t ) = c i ( t ) + α i e i 1 ( t ) with i [ 1 , N u ]
with α i defined in (11). Then, we build:
e ˜ i ( t ) = e ¯ i ( t ) + β i e ˜ i 1 ( t ) ,
where e ¯ 0 ( t ) = e ˜ 0 ( t ) = x ( t ) (this signal will be excluded from the set later) and for i > 0 :
e ¯ i ( t ) = x ( t ) e i ( t ) .
Moreover:
β i = e ¯ i ( t ) , e ˜ i 1 ( t ) | | e ˜ i 1 ( t ) | | 2 2 = T s 2 T s 2 e ¯ i ( t ) e ˜ i 1 * ( t ) d t | | e ˜ i 1 ( t ) | | 2 2 .
The final waveform assigned to each user is obtained by setting g i ( t ) = e ˜ i ( t ) .

3.4. MU-CSS Gram–Schmidt Insertion

In this last variant, we combine the previous method with the insertion of a HFM signal at regular intervals, such as:
e ¯ i ( t ) = x ( t ) if i = k p with k N * x ( t ) e i ( t ) else ,
where p is the insertion step. The idea is to try to improve the robustness of the different waveforms. To impose orthogonality between spread signals, we simply apply Equations (13) and (14), and finally get g i ( t ) = e ˜ i ( t ) .

4. Simulation Results

4.1. Underwater Acoustic Channel Simulator

For the simulation comparisons, we consider the UWA channel simulator provided by [15] based on a stochastic model. The time-varying transfer function for the i-th user is given by:
H i ( f , t ) = H ¯ i ( f ) p h i , p γ i , p ( f , t ) e j 2 π f τ i , p ( t ) ,
where H ¯ i ( f ) is the transfer function of direct path, h i , p is the relative path gain, γ i , p ( f , t ) represents the scattering coefficient modeled by a complex-valued Gaussian processes, whose statistics reflect the time coherence of the channel, and τ i , p ( t ) denotes time-varying delay of the p-th path and can be expressed as:
τ i , p ( t ) = τ ¯ i , p ( a ¯ i + a i , p ) t ,
where τ ¯ i , p is the average delay of path p, and a ¯ i represents the mean Doppler shift induced by the motion of the i-th AUV relative to the receiver. In the following, we will assume that a ¯ i is known at the receiver side and compensated. Moreover, a i , p is the residual Doppler factor that captures resulting motion-induced time scaling on the p-th path. Coefficients a i , p are assumed to be constant over a frame and to follow a zero-mean Gaussian distribution with variance σ a 2 . Time variations of γ i , p ( f , t ) and τ i , p ( t ) lead to Doppler spread effects [15].

4.2. System Parameters

The chosen model represents a short-range UWA transmission with a 10 m water depth at a center frequency of 23 kHz over a 4 kHz bandwidth. Each AUV is assumed to be at the same depth of 1 m. At the beginning of the simulations, the range between each AUV and the receiver is randomly selected in the interval [ 0.1 , 1 ] km, modeling a fleet situating in a circular area (Figure 3). The channel model parameters are summarized in Table 1, whereas the transmission system parameters are provided in Table 2. The symbol duration is set according to the channel delay spread such that T s > τ m a x and is fixed identically for all protocols. The evolution of the simulated channel impulse response | h i ( τ , t ) | over one frame is provided in Figure 4.

4.3. Orthogonality Verification

To verify the orthogonality of the proposed waveform, we compute the Signal-to-Interference-plus-Noise Ratio (SINR) obtained after matched filtering. Following (7), for user i, we have:
SINR = E | γ i , k | 2 E | η i , k | 2 + E | w k | 2 .
Simplifying (A2), (A3), and (A4) in the case of static AUV motion (i.e., a i = 0 ) and very small channel delay spread compared to the symbol duration (i.e., T s > > τ m a x ), the last equation becomes:
SINR = | T s 2 T s 2 g i * ( t ) ( + h i ( t , τ ) g i ( t τ ) d τ ) d t | 2 | j = 1 j i N u T s 2 T s 2 g i * ( t ) ( + h j ( t , τ ) g j ( t τ ) d τ ) d t | 2 + E | T s 2 T s 2 g i * ( t ) n ( t ) d t | 2 .
In Figure 5, we numerically compute the SINR by using (21) and the system parameters depicted in Table 2 over an Additive White Gaussian Noise (AWGN) channel and over the time-varying UWA channel with static users described in Section 4.1. Comparisons are performed between MU-CSS, CDMA, and TDMA transmissions. At the N u = 1 user, since there are no interference terms, all of the transmission techniques have the same SINR after matched filter decoding, which is equal to channel Signal-to-Noise Ratio (SNR) added to the spreading gain in the case of AWGN channel. Naturally, as the the number of users increases, SINR decreases due to the growing importance of the interference terms, except for the TDMA case, for which the interference terms are absent whatever the number of users, thanks to time multiplexing. In both the AWGN and UWA channels, MU-CSS transmissions outperform CDMA, demonstrating that the Gram–Schmidt-based construction method provides good orthogonality properties for MU-CSS waveforms. This SINR gap is mainly explained by the use of PN sequences in CDMA that are not perfectly orthogonal (but only quasi-orthogonal), while MU-CSS employs waveforms that are orthogonal, owing to the Gram–Schmidt process. Obviously, this SINR gap could be erased in AWGN by the use of orthogonal codes like Walsh–Hadamard sequences for CDMA; however, such codes are not suitable in the uplink scenario.

4.4. Performance Metrics

As performance metrics, we consider the average effective data rate per user, defined for each transmission technique as follows:
D e CDMA = R C log 2 M N S F . T c · ( 1 FER ) [ bps ]
D e TDMA = R C log 2 M N u N S F T c + ( N u 1 ) T g · ( 1 FER ) [ bps ]
D e MU CSS = R C log 2 M T h · ( 1 FER ) [ bps ] ,
where M is the size of the DPSK constellation, R C is channel coding rate, and FER is the Frame Error Rate. A frame is considered erroneous when at least one bit per frame after channel decoding is erroneous.

4.5. Static Channel

In a first step, we consider a static UWA channel leading to only frequency-selective fading. This yields the constant parameters γ p ( f , t ) and τ p ( t ) in time in Equation (18). The Frame Error Rate (FER) performance and effective data rate of each transmission technique over the modeled shallow-water acoustic channel are provided in Figure 6.
In the single-user scenario, the three transmission techniques have an FER of 0, and, as expected, the FER of TDMA remains unchanged when the number of users increases. With more than four users, the interfering terms of the CDMA, expressed in Equation (7) by the quantity η i , k , make the decoding of each user impossible. On the other side, the largest number of users that can be handled by the MU-CSS is eight or nine, depending on the method. The fact that MU-CSS outperforms CDMA is mainly explained by the better orthogonality properties of the MU-CSS waveforms.

4.6. Time-Varying Channel and Static Users

In a second step, we consider a time-varying channel model, where Doppler spread effect is provided in the Equation (18) by the γ i ( f , t ) and τ i , p ( t ) coefficients. In this scenario, we assume that all users are static, yielding a ¯ i = 0 in Relation (19). The performance of the time-varying channel with static users is depicted in Figure 7.
The Doppler spread effect provided by multipath time variations leads to an FER increase of both the CDMA and MU-CSS transmissions, while the decoding performance TDMA still remains error-free. In fact, the TDMA transmission is not affected by multiuser interference, but only by UWA channel time and frequency selectivity, while CDMA and MU-CSS suffer from multiuser interference in addition to the UWA channel selectivity. The MU-CSS transmissions have the best effective data rate compared to CDMA because the HFM signal makes the spreading signals resistant against channel impairments, such as Doppler spread. Among the MU-CSS transmission technique, the Gram–Schmidt iterated method appears to be slightly less robust than the other methods.

4.7. Time-Varying Channel and Mobile Users

In a last step, we consider a time-varying channel model with mobile AUVs whose speeds are randomly selected in the interval [ 2 , 2 ] m/s at each frame and for each user. The motion-induced Doppler shift is assumed to be perfectly known and compensated at the reception for each user i. According to (4), since each user has different speed, Doppler compensation of user i will increase the power of the interference terms. However, in practice, Doppler shift is unknown and must be estimated prior to decoding [26].
The performance of the time-varying UWA channel with mobile users is shown in Figure 8. In the single-user scenario, the three transmission techniques provide a FER of 0% and, as expected, the FER of TDMA remains unchanged when the number of users increases. Both the CDMA and MU-CSS transmissions are severely impacted by motion-induced Doppler shift, since Doppler shift correction for an user also applies to other users, according to Equation (4). However, the MU-CSS transmission still outperforms CDMA, which might be explained by the MU-CSS construction, which provides both an orthogonality enhancement and a better robustness against Doppler shift. Beyond six users, the TDMA approach is more efficient in terms of data rate.

5. Experimental Results

5.1. Channel Sounding

Ty-Colo Lake of Saint-Renan (France)

The sounding experiments took place in July 2019 at the lake of Ty-Colo, Saint Renan, France. The depth of the lake is around 5 m, and up to 10 transmission ranges between [ 47 , 364 ] m were sounded successively with one hydrophone at the receiver side, as depicted in Figure 9. Each channel sounding was performed over 3 min and 30 s, using a 255 Maximal Length Sequence (MLS) probe signal [27] centered on f c = 27 kHz over a 6 kHz bandwidth. Figure 10 provides an example of the delay–Doppler spread extracted from the successive estimated Channel Impulse Response (CIR). Estimated channel delay spreads and Doppler spreads are reported in Table 3.

5.2. Watermark Replay Channel

To simulate a real experiment, we consider in this section the Watermark channel [17], which is a replay channel simulator driven by measurements of the time-varying CIR. The principle of the simulator consists of distorting input waveforms by convolving them with measured channels. To simulate a multiuser communication, we sum the output of several Watermark channels fed by different CIRs and delayed by the relative range of each user. The operation of the channel replay for a static multiuser communication in the Single Input Single Output (SISO) case can be expressed in baseband as:
r ( t ) = i = 1 N u + h ^ i ( τ , t ) s i ( t τ τ ¯ i ) d τ + n ( t ) ,
where s i ( t ) is the input signal, h ^ i ( τ , t ) is the recorded CIR of the i-th user, τ ¯ i is communication delay between the i-th user and the receiver, and n ( t ) is Gaussian noise.
For mobile multiuser communication, the Doppler shift is simulated by resampling and phase-rotating the transmitted signal as follows:
r ( t ) = i = 1 N u + h ^ i ( τ , t ) s i ( 1 a i ) ( t τ τ ¯ i ) e j 2 π f c a i ( t τ ) d τ + n ( t ) .
In the following, the Doppler shift will be known by the receiver and compensated by the relation (A1). The Ty-Colo lake channel parameters are summarized in Table 3, whereas the transmission system parameters are provided in Table 2.

5.3. Performance Results

5.3.1. Static Users

Figure 11 provides the performance of multiuser transmission techniques over the Watermark channel fed by the Ty-Colo lake channel soundings. It can be noticed that the FER, and consequently the effective data rates, are worse than in the simulation. This can mainly be explained by the fact that the experimental soundings are in very shallow water ( 5 m), leading to a much more important multipath effect and, as a consequence, to higher multiple-access interference terms. Meanwhile, the FER performance of MU-CSS transmissions are still better than CDMA ones with up to six simultaneous users (except the case N u = 4 , where the CDMA is slightly ahead). Beyond this threshold, TDMA transmission is more suitable, despite its low data rate due to a large number of users.

5.3.2. Mobile Users

In Figure 12, the AUVs’ motion is emulated by adding a motion-induced Doppler scale at the output of the Watermark channel. For each frame, the speed value of each AUV is randomly selected in the interval of [ 2 , 2 ] m/s. We can see that the performances of all access schemes are degraded, except for TDMA. With one to six users, the MU-CSS transmissions remain globally more interesting in terms of the effective data rate. As seen in the simulation, the MU-CSS with the Gram–Schmidt insertion method is confirmed in experiments to provide the highest robustness of all of the MU-CSS construction methods. Beyond six users, the TDMA is demonstrated to be more advantageous.

6. Conclusions and Future Works

In this paper, we have proposed a new multiuser transmission technique based on a HFM signal, denoted MU-CSS, in the context of UWA communication within an AUV fleet. By using the Gram–Schmidt orthogonalization, we derived three construction methods for MU-CSS, allowing a very simple matched filter decoding scheme at the receiver side. Simulation comparisons against traditional CDMA, with a single user decoding over static and time-varying shallow water UWA models, demonstrate a superior effective data rate for the proposed MU-CSS scheme, even if the number of users is large and even if the users are in motion, as for an AUV fleet. The experimental results with the Watermark channel replay fed by channel soundings confirm the superiority of MU-CSS transmissions in a realistic scenario. The MU-CSS is demonstrated to be globally superior to CDMA for up to six users. Furthermore, the traditional TDMA approach is demonstrated to be more efficient. The MU-CSS approach, especially when associated with the Gram–Schmidt construction method, offers a set of waveforms that provide good orthogonal properties even in the UWA uplink channel, so that such waveforms do not require a complex multiuser decoding scheme at the receiver side. Thereby, MU-CSS transmission techniques constitute an interesting alternative to asynchronous CDMA for UWA networks.
In a future work, we will consider multi-channel decoding for MU-CSS in order to improve the number of users to be simultaneously correctly decoded, and also take into account real Doppler-shift estimation and its impact on decoding performance when AUVs have different speeds and directions.

Author Contributions

Software, C.B.; writing—original draft preparation, C.B. and P.-J.B.; writing—review and editing, A.P. and P.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partly funded by Thales DMS France in the framework of the WAVES laboratory. The APC was funded by L@bISEN Yncréa Ouest.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Calculation of γi,k, ηi,k and wi,k

The received baseband signal after Doppler pre-processing can be expressed as:
z i ( t ) = r t 1 a i e j 2 π f c ( a i 1 a i ) t    = j = 1 N u + h j τ , t 1 a i s j ( 1 a j ) t 1 a i τ e j 2 π f c a j ( t 1 a i τ ) d τ e j 2 π f c ( a i 1 a i ) t    + n t 1 a i e j 2 π f c ( a i 1 a i ) t .
Combination of (2) and (7) yields:
γ i , k = T s 2 T s 2 + h i τ , t + k T s 1 a i g i * ( t ) g i t ( 1 a i ) τ e j 2 π f c a i τ d τ d t
η i , k = n = 1 n k N s d i , n T s 2 T s 2 + h i τ , t + k T s 1 a i g i * ( t ) g i t τ ( n k ) T s e j 2 π f c a i τ d τ d t    + j = 1 j i N u n = 1 N s d j , n T s 2 T s 2 + h j τ , t + k T s 1 a i g i * ( t ) g j ( 1 a j ) t + k T s 1 a i τ n T s e j 2 π f c a i a j 1 a i k T s + t + a j τ d τ d t
and
w i , k = e j 2 π f c a i 1 a i k T s T s 2 T s 2 g i * ( t ) n t + k T s 1 a i e j 2 π f c a i 1 a i t d t

Appendix B. Justification of the MU-CSS Gram–Schmidt Construction Process

To have the orthogonality between the different e i ( t ) , we use a variant of the Gram–Schmidt process [25], which is a method for orthogonalizing a set of vectors in an inner product space. The inner product is defined by f , g L 2 ( R ) as f ( t ) , g ( t ) = T s 2 T s 2 f ( t ) g * ( t ) d t . Let { c 1 ( t ) , c 2 ( t ) } be a set of linearly independent vectors. We add the vector e 0 ( t ) to the previous family and we build an orthogonal family from vector e 0 ( t ) . By the Gram–Schmidt process, we have:
e 1 ( t ) = c 1 ( t ) + α 1 e 0 ( t ) .
Using orthogonality, the previous equation gives:
c 1 ( t ) , e 0 ( t ) + α 1 | | e 0 ( t ) | | 2 2 = 0
α 1 = c 1 ( t ) , e 0 ( t ) | | e 0 ( t ) | | 2 2 = T s 2 T s 2 c 1 ( t ) e 0 * ( t ) d t | | e 0 ( t ) | | 2 2 .
For the last vector, the Gram–Schmidt process gives:
e 2 ( t ) = c 2 ( t ) + β e 0 ( t ) + α 2 e 1 ( t ) .
We take β = 0 because that is enough to have orthogonality, and we obtain:
e 2 ( t ) = c 2 ( t ) + α 2 e 1 ( t ) .
Using orthogonality, the previous equation becomes:
c 2 ( t ) , e 1 ( t ) + α 2 | | e 1 ( t ) | | 2 2 = 0
α 2 = c 2 ( t ) , e 1 ( t ) | | e 1 ( t ) | | 2 2 = T s 2 T s 2 c 2 ( t ) e 1 * ( t ) d t | | e 1 ( t ) | | 2 2
By generalization, we deduce the Equation (10).

References

  1. Stojanovic, M. Underwater Acoustic Communications: Design Considerations on the Physical Layer. In Proceedings of the 2008 Fifth Annual Conference on Wireless on Demand Network Systems and Services, Garmisch-Partenkirchen, Germany, 23–25 January 2008; pp. 1–10. [Google Scholar] [CrossRef] [Green Version]
  2. Stojanovic, M.; Beaujean, P.P.J. Acoustic Communication. In Springer Handbook of Ocean Engineering; Dhanak, M.R., Xiros, N.I., Eds.; Springer International Publishing: Cham, Switzerland, 2016; pp. 359–386. [Google Scholar] [CrossRef]
  3. Champion, B.T.; Joordens, M.A. Underwater swarm robotics review. In Proceedings of the 2015 10th System of Systems Engineering Conference (SoSE), San Antonio, TX, USA, 17–20 May 2015; pp. 111–116. [Google Scholar] [CrossRef]
  4. Abramson, N. Development of the ALOHANET. IEEE Trans. Inf. Theory 1985, 31, 119–123. [Google Scholar] [CrossRef]
  5. Chirdchoo, N.; Soh, W.S.; Chua, K.C. Aloha-Based MAC Protocols with Collision Avoidance for Underwater Acoustic Networks. In Proceedings of the IEEE INFOCOM 2007, 26th IEEE International Conference on Computer Communications, Anchorage, AK, USA, 6–12 May 2007; pp. 2271–2275. [Google Scholar] [CrossRef]
  6. Colvin, A. CSMA with collision avoidance. Comput. Commun. 1983, 6, 227–235. [Google Scholar] [CrossRef]
  7. Otnes, R.; Asterjadhi, A.; Casari, P.; Goetz, M.; Husøy, T.; Nissen, I.; Rimstad, K.; Walree, P.v.; Zorzi, M. Underwater Acoustic Networking Techniques; Springer Briefs in Electrical and Computer Engineering; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  8. Trivedi, V.K.; Kumar, P. Carrier Interferometry Coded Single Carrier FDMA (CI/SC-FDMA) for Next Generation Underwater Acoustic Communication. Wireless Personal Commun. 2017, 95, 4747–4762. [Google Scholar] [CrossRef]
  9. Pompili, D.; Melodia, T.; Akyildiz, I.F. A CDMA-based Medium Access Control for UnderWater Acoustic Sensor Networks. IEEE Trans. Wireless Commun. 2009, 8, 1899–1909. [Google Scholar] [CrossRef]
  10. Konstantakos, D.; Tsimenidis, C.; Adams, A.; Sharif, B. Comparison of DS-CDMA and MC-CDMA techniques for dual-dispersive fading acoustic communication networks. IEE Proc. Commun. 2005, 152, 1031–1038. [Google Scholar] [CrossRef]
  11. Stojanovic, M.; Freitag, L. Multichannel Detection for Wideband Underwater Acoustic CDMA Communications. IEEE J. Oceanic Eng. 2006, 31, 685–695. [Google Scholar] [CrossRef] [Green Version]
  12. Yang, T.C. Spatially Multiplexed CDMA Multiuser Underwater Acoustic Communications. IEEE J. Oceanic Eng. 2016, 41, 217–231. [Google Scholar] [CrossRef]
  13. Yuan, F.; Wei, Q.; Cheng, E. Multiuser chirp modulation for underwater acoustic channel based on VTRM. Int. J. Naval Architect. Ocean Eng. 2016, 9. [Google Scholar] [CrossRef] [Green Version]
  14. Bernard, C.; Bouvet, P.J. Multiuser Underwater Acoustic Communication for an AUV Fleet. In Proceedings of the OCEANS 2019—Marseille, Marseille, France, 17–20 June 2019. [Google Scholar]
  15. Qarabaqi, P.; Stojanovic, M. Statistical Characterization and Computationally Efficient Modeling of a Class of Underwater Acoustic Communication Channels. IEEE J. Oceanic Eng. 2013, 38, 701–717. [Google Scholar] [CrossRef]
  16. Aval, Y.M.; Wilson, S.K.; Stojanovic, M. On the Achievable Rate of a Class of Acoustic Channels and Practical Power Allocation Strategies for OFDM Systems. IEEE J. Oceanic Eng. 2015, 40, 785–795. [Google Scholar] [CrossRef]
  17. Van Walree, P.; Socheleau, F.X.; Otnes, R.; Jenserud, T. The Watermark Benchmark for Underwater Acoustic Modulation Schemes. IEEE J. Oceanic Eng. 2017, 42, 1007–1018. [Google Scholar] [CrossRef] [Green Version]
  18. Aval, Y.M.; Stojanovic, M. Differentially Coherent Multichannel Detection of Acoustic OFDM Signals. IEEE J. Oceanic Eng. 2015, 40, 251–268. [Google Scholar] [CrossRef]
  19. Aval, Y.M.; Wilson, S.K.; Stojanovic, M. On the Average Achievable Rate of QPSK and DQPSK OFDM Over Rapidly Fading Channels. IEEE Access 2018, 6, 23659–23667. [Google Scholar] [CrossRef]
  20. Proakis, J.G.; Salehi, M. Digital Communications, 5th ed.; McGraw-Hill: Boston, MA, USA, 2008. [Google Scholar]
  21. Ziemer, R.E.; Peterson, R.L.; Borth, D.E. Introduction to Spread-Spectrum Communications; Prentice Hall: Upper Saddle River, NJ, USA, 1995. [Google Scholar]
  22. Zhou, S.; Wang, Z. OFDM for Underwater Acoustic Communications; John Wiley & Sons Ltd.: Chichester, UK, 2014. [Google Scholar] [CrossRef]
  23. Kebkal, K.G.; Bannasch, R. Sweep-spread carrier for underwater communication over acoustic channels with strong multipath propagation. J. Acoust. Soc. Am. 2002, 112, 2043–2052. [Google Scholar] [CrossRef] [PubMed]
  24. Kaminsky, E. Chirp signaling offers modulation scheme for underwater communications. SPIE Newsroom 2006. [Google Scholar] [CrossRef] [Green Version]
  25. Lang, S. Introduction to Linear Algebra; Springer Inc.: New York, NY, USA, 1985. [Google Scholar]
  26. Sharif, B.S.; Neasham, J.; Hinton, O.R.; Adams, A.E. A computationally efficient Doppler compensation system for underwater acoustic communications. IEEE J. Oceanic Eng. 2000, 25, 52–61. [Google Scholar] [CrossRef]
  27. Van Walree, P. Channel Sounding for Acoustic Communications: Techniques and Shallow-Water Examples; Tech. Rep. FFI-Rapport; Norwegian Defence Research Establishment (FFI): Kjeller, Norway, 2011; Volume 7, p. 58. [Google Scholar]
Figure 1. Scheme of Time-Division Multiple Access (TDMA).
Figure 1. Scheme of Time-Division Multiple Access (TDMA).
Sensors 20 01527 g001
Figure 2. Instantaneous frequency of the Hyperbolically Frequency-Modulated (HFM) waveform with B T s = 31 , 63, 127.
Figure 2. Instantaneous frequency of the Hyperbolically Frequency-Modulated (HFM) waveform with B T s = 31 , 63, 127.
Sensors 20 01527 g002
Figure 3. Scheme of the simulated system.
Figure 3. Scheme of the simulated system.
Sensors 20 01527 g003
Figure 4. Evolution of the simulated time-varying channel impulse response over one frame based on the parameters provided in Table 1.
Figure 4. Evolution of the simulated time-varying channel impulse response over one frame based on the parameters provided in Table 1.
Sensors 20 01527 g004
Figure 5. Average Signal-to-Interference-plus-Noise Ratio (SINR) for different waveforms over the Additive White Gaussian Noise (AWGN) and time-varying UWA channels with static users, SNR = 30 dB.
Figure 5. Average Signal-to-Interference-plus-Noise Ratio (SINR) for different waveforms over the Additive White Gaussian Noise (AWGN) and time-varying UWA channels with static users, SNR = 30 dB.
Sensors 20 01527 g005
Figure 6. Average Frame Error Rate (FER) performance (left) and effective data rate (right) versus the number of users over the static simulated UWA channel model.
Figure 6. Average Frame Error Rate (FER) performance (left) and effective data rate (right) versus the number of users over the static simulated UWA channel model.
Sensors 20 01527 g006
Figure 7. Average FER performance (left) and effective data rate (right) versus number of users over the time-varying simulated UWA channel model with static users.
Figure 7. Average FER performance (left) and effective data rate (right) versus number of users over the time-varying simulated UWA channel model with static users.
Sensors 20 01527 g007
Figure 8. Average FER performance (left) and effective date rate (right) versus number of users over the time-varying UWA channel model with mobile users.
Figure 8. Average FER performance (left) and effective date rate (right) versus number of users over the time-varying UWA channel model with mobile users.
Sensors 20 01527 g008
Figure 9. Experiment scheme on the Ty-Colo lake of St-Renan.
Figure 9. Experiment scheme on the Ty-Colo lake of St-Renan.
Sensors 20 01527 g009
Figure 10. Delay–Doppler spread function for the Ty-Colo lake.
Figure 10. Delay–Doppler spread function for the Ty-Colo lake.
Sensors 20 01527 g010
Figure 11. Average FER performance (left) and effective data rate per user (right) versus number of users for the Ty-Colo lake with the replayed channel with static users.
Figure 11. Average FER performance (left) and effective data rate per user (right) versus number of users for the Ty-Colo lake with the replayed channel with static users.
Sensors 20 01527 g011
Figure 12. Average FER performance versus number of users for the Ty-Colo lake with the replayed channel with mobile users (left), and average effective data rate per user versus number of users for the Ty-Colo lake with the replayed channel with mobile users (right).
Figure 12. Average FER performance versus number of users for the Ty-Colo lake with the replayed channel with mobile users (left), and average effective data rate per user versus number of users for the Ty-Colo lake with the replayed channel with mobile users (right).
Sensors 20 01527 g012
Table 1. Underwater Acoustic (UWA) channel model parameters.
Table 1. Underwater Acoustic (UWA) channel model parameters.
SymbolSignificationValue
f c Center frequency23 kHz
N u Number of AUVs [ 1 , 10 ]
f s Sample frequency100 kHz
BSignal bandwidth4 kHz
D i Transmission range [ 0.1 , 1 ] km
z w Water depth10 m
τ m a x RMS channel delay spread [20] [ 0.52 , 0.84 ] ms
SNRSignal to noise ratio10 dB
v i User relative speed [ 2 , 2 ] m/s
σ a Residual motion-induced Doppler spread standard deviation 10 5
Table 2. System parameters.
Table 2. System parameters.
SymbolSignificationValue
MModulation order2 (DBPSK)
N s Number of symbols per frame200
N f Number of frames5000
C FEC code typeConvolutive code
g C FEC code generator( 133 , 171 ) o
R C FEC code rate 1 2
T g Guard interval time15 ms
T h Duration of the chirp signal 7.75 ms
T c Chip duration 0.25 ms
N S F PN length code31
N I T Number of iterations1000
pInsertion step7
α Pulse shaping filter roll-off factor 0.25
T s Symbol duration 7.75 ms
Table 3. Watermark channel parameters.
Table 3. Watermark channel parameters.
SymbolSignificationValue
f c Center frequency27 kHz
f s Sampling frequency96 kHz
BSignal bandwidth4 kHz
D i Transmission range [ 47 , 364 ] m
z w Water depth5 m
SNRSignal-to-noise ratio10 dB
τ m a x RMS channel delay spread [20] [ 4.31 , 7.27 ] ms
σ m a x RMS channel Doppler spread [20] [ 0.86 , 2.51 ] Hz

Share and Cite

MDPI and ACS Style

Bernard, C.; Bouvet, P.-J.; Pottier, A.; Forjonel, P. Multiuser Chirp Spread Spectrum Transmission in an Underwater Acoustic Channel Applied to an AUV Fleet. Sensors 2020, 20, 1527. https://doi.org/10.3390/s20051527

AMA Style

Bernard C, Bouvet P-J, Pottier A, Forjonel P. Multiuser Chirp Spread Spectrum Transmission in an Underwater Acoustic Channel Applied to an AUV Fleet. Sensors. 2020; 20(5):1527. https://doi.org/10.3390/s20051527

Chicago/Turabian Style

Bernard, Christophe, Pierre-Jean Bouvet, Antony Pottier, and Philippe Forjonel. 2020. "Multiuser Chirp Spread Spectrum Transmission in an Underwater Acoustic Channel Applied to an AUV Fleet" Sensors 20, no. 5: 1527. https://doi.org/10.3390/s20051527

APA Style

Bernard, C., Bouvet, P. -J., Pottier, A., & Forjonel, P. (2020). Multiuser Chirp Spread Spectrum Transmission in an Underwater Acoustic Channel Applied to an AUV Fleet. Sensors, 20(5), 1527. https://doi.org/10.3390/s20051527

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop