Next Article in Journal
Components of Artificial Neural Networks Realized in CMOS Technology to be Used in Intelligent Sensors in Wireless Sensor Networks
Previous Article in Journal
A Robust Approach to TDMA Synchronization in Aerial Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Chirp-Based FHSS Receiver with Recursive Symbol Synchronization for Underwater Acoustic Communication

1
Department of the Radio Communication Engineering, Korea Maritime and Ocean University, Busan 49112, Korea
2
Agency for Defense Development, Changwon 51678, Korea
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(12), 4498; https://doi.org/10.3390/s18124498
Submission received: 18 October 2018 / Revised: 10 December 2018 / Accepted: 11 December 2018 / Published: 19 December 2018
(This article belongs to the Section Physical Sensors)

Abstract

:
In this paper, we propose a covert underwater acoustic communication method that is robust to fading using a chirp signal combined with a frequency-hopping spread spectrum scheme. A fractional Fourier transform, which estimates the slope of the signal frequency variation, is applied to the receiver to enable a robust and reliable symbol estimation with respect to the frequency and irregular phase variations. In addition, since the recursive symbol synchronization can be implemented using a chirp signal, compression and expansion effects due to the Doppler shift can be mitigated. Simulation and lake trials were performed to verify the performance of the proposed method. The simulation was performed by two different methods.

1. Introduction

The upcoming 5G communication standard has been an active area of research in recent years, focusing on terrestrial communications. However, progress in underwater communications has been comparatively much slower. The underwater acoustic environment is much more challenging, with far greater dispersion and multipath effects, limiting the potential use of the classical 5G scheme in the underwater domain [1]. Some research in developing underwater technology has taken place with a focus on orthogonal frequency division multiplexing (OFDM) and single-carrier frequency division multiple access (SC-FDMA) schemes [2,3]. In general, the underwater acoustic communication community has focused on increasing transmission distance, improving bandwidth efficiency, and reducing the rate of error in the complex time–space-varying underwater channel [4,5]. Covert underwater acoustic communication has been studied for military purposes to hide positioning information or data information. Although it is not an acoustic signal, underwater optical communication using a blue–green laser has also been aimed at enhancing covertness through the use of the strong directivity of lasers [6]. In addition, biomimetic signaling has been studied as a scheme for underwater acoustic communications, but the subject is in the early stages [7,8]. On the other hand, the spread spectrum schemes are widely used in the field to achieve covert communication. There are two typical spread spectrum schemes. One is a direct-sequence spread spectrum (DSSS) method in which the data are multiplied by a spread factor [9,10,11,12]. The other is a frequency-hopping spread spectrum (FHSS) method in which the frequency band is shifted using a spread code [13]. The FHSS method is randomly hopped according to the hopping pattern. By this method, an interceptor who does not know the hopping pattern cannot recover the signal. Furthermore, if any jamming is executed with malicious intent, the signal performance will be maintained because the center frequency will be continuously changed. This capability is called anti-jamming [14,15]. A simple multiple-access underwater acoustic communication protocol called JANUS was designed and tested as a NATO standard. It was a binary frequency shift keying (BFSK)-based FH method [16]. A covert underwater acoustic communication scheme between an unmanned underwater vehicle (UUV) and a distant mother platform using an OFDM had been studied [17]. However, the previous FHSS method was of the form in which a single center frequency was hopped. A chirp-based method utilizes a data-modulated signal, which has its energy spread over a bandwidth that is greater than the rate of information being sent, like all spread spectrum systems. The chirp-based method is one in which the binary information modulates the slope of a linear chirp, with chirps representing zeros or ones. This method alone can be considered a covert underwater acoustic communication technology. Due to the wide bandwidth in the chirp-based method, immunity against various frequency-selective fading and a frequency diversity effect may be expected [18,19].
Fractional Fourier transform (FrFT) is a generalized analysis method of the existing Fourier transform concept [20,21]. In particular, when a chirp signal has a linear characteristic in the time–frequency domain, robust detection is possible in a noisy or reverberant environment, as the FrFT results can obtain a concentrated energy spectrum. The FrFT has been applied in various fields, such as radar, image processing and communication [22,23,24,25]. Since underwater acoustic communication signals must typically pass through poor underwater channels, the coherent methods that are sensitive to phase shifts complicate the structure of the receiver. On the other hand, the non-coherent methods, such as the FrFT, can mitigate the fading effects with a relatively simple receiver structure.
The chirp signal combined with FHSS and demodulated by the FrFT method was studied in recent times [26]. This method simultaneously transmitted two chirp signals in one symbol duration and demodulated the received frequency-hopped signals in the FrFD (fractional Fourier domain) of the FrFT spectrum. As the up-slope chirp signal and the down-slope chirp signal were applied in one symbol duration, a parallel receiver should have retrieved the different transform orders.
In this paper, we propose the single chirp-based FHSS method demodulated by the non-coherent FrFT method. In the proposed method, either the up-slope chirp signal or the down-slope chirp signal is transmitted within one symbol duration, and does not require the parallel demodulator. There is also a difference from the conventional method in symbol demapping using the FrFT spectrum. The conventional method was independently demodulated using different orders for the up-slope chirp and the down-slope chirp without considering the negative frequency domain in the fractional frequency domain. On the other hand, since the proposed method uses the symmetrical characteristic of the frequency (having the positive frequency and negative frequency), there is a difference in the form represented by the FrFT due to the slope inversion of the positive and the negative frequency domains. In addition, as only a single transform order is required, a simpler receiver structure can be used.
Due to the movement, a Doppler shift can occur to the received signal, and the Doppler frequency will be relatively large due to the slow acoustic velocity in water compared to the propagation velocity on the ground [27]. Since the received signal is expanded and compressed by the Doppler shift, a method of implementing symbol synchronization recursively using an autocorrelation function of the estimated chirp signal should be employed. To consider the performance of the proposed method, we conducted simulations based on an underwater channel propagation model and conducted a lake trial to verify its validity.
The remainder of this paper is organized as follows: In the next section, the principle of the single chirp-based FHSS scheme is presented. Section 3 describes the receiver based on the FrFT. Section 4 and Section 5 then show the performance of the proposed method as evaluated via simulations and lake trials. Section 6 provides the conclusion.

2. Chirp-Based FHSS Scheme

The basic chirp signal used in the proposed method is expressed in Equation (1) as
s ( t ) = exp [ j 2 π ( f 0 t + k 2 t 2 ) ] ,
where f 0 is the initial frequency and k is the chirp rate. Equation (1) can be rewritten to express the chirp signal as follows,
s u ( t ) = exp [ j 2 π t ( f m i n + f m a x f m i n T t ) ] ,
s d ( t ) = exp [ j 2 π t ( f m a x + f m i n f m a x T t ) ] ,
where f m i n is the minimum frequency and f m a x is the maximum frequency. T is the duration of the symbol. Equation (2) represents a signal that has an up-slope frequency with an initial frequency f m i n . Equation (3) represents a signal that has a down-slope frequency with an initial frequency f m a x . To express the chirp signal in real numbers, the above equations can be expressed as
S u ( t ) = R e { s u ( t ) } ,
S d ( t ) = R e { s d ( t ) } .
The chirp-based scheme selects one of the chirp symbols of Equation (4) using the binary sequence b i . Therefore, the chirp-based signal s c ( t ) with the symbol period T is expressed as [28]
s c ( t ) = { S u ( t ) ,   i f   b i = 0 S d ( t ) ,   i f   b i = 1 ,   n T t ( n + 1 ) T .
This equation implies assigning a binary value “1” to a signal with an up-slope frequency, and assigning a binary value “0” to a signal with a down-slope frequency. During the duration of the symbol, the chirp signal is mapped by the binary bit stream. The received signal is usually used by a correlator receiver.
Figure 1 shows the basic structure of the hopping scheme in a chirp-based FHSS. In this case, since the hopping duration is longer than the symbol duration, there are several symbols in the single hopping interval.

3. Receiver Design Based on Fractional Fourier Transform

Figure 2 shows a block diagram of the proposed method. Firstly, a hopping pattern is generated and multiplied by the modulated transmission signal to perform frequency hopping, and the receiver utilizes a bandpass filter with a known hopping pattern. Then, the compression or expansion effect of the signal due to the Doppler shift can be corrected by introducing the matched filter results between the symbol estimation through the FrFT spectrum and the received signal.
In the example shown in Figure 1, there are four symbols in a single hopping frequency interval. Each symbol is divided into a low-frequency band and a high-frequency band. The symbol with an up-slope in the low-frequency band is called “Symbol Mapping #1”, and the symbol with a down-slope is called “Symbol Mapping #2”. The symbol with up-slope in the high-frequency band is called “Symbol Mapping #3”, and the symbol with down-slope is called “Symbol Mapping #4”.
Figure 3a shows the matched filter result, which is from the auto-correlation function of the “Symbol Mapping #1” in any hopping pattern. It can be seen that the output of the matched filter has a maximum value at the ideal point. Figure 3b shows the matched filter result, which is the cross-correlation function between the “Symbol Mapping #1” and “Symbol Mapping #2” in any hopping pattern. Since “Symbol Mapping #1” and “Symbol Mapping #2” have different slopes, the maximum point of the matched filter result cannot be correctly represented, and the magnitude is less than the ideal value [29].
Figure 4a shows the matched filter result from the cross-correlation function between “Symbol Mapping #1” and “Symbol Mapping #3” in any hopping pattern. The mapping types have the same slopes but show a lower magnitude because the frequency bands are different. Figure 4b shows the matched filter result as a cross-correlation function between “Symbol Mapping #1” and “Symbol Mapping #4” in any hopping pattern. The mapping types have the same frequency, but show lower values due to their different slopes.
In this paper, we propose a method to estimate the symbol in the received signal by applying the FrFT scheme. The FrFT is a generalization of the conventional Fourier transform concept, which makes it possible to convert the time domain signal into a time–frequency intermediate domain. That is, when the time domain is represented by the x-axis and the frequency domain by the y-axis, the FrFT is the result obtained by rotating an arbitrary angle ∅ in the counterclockwise direction. The FrFT can obtain the spectra with various characteristics depending on the transform orders defined as α = 2∅/π. If the different signals have different frequency slopes, the signals will have peaks at different positions with different transform orders on the FrFT spectrum. In this way, the correct transform order with peaks in each FrFT spectrum is called the optimal transform order, which is denoted by α o p t [27,30]. The FrFT equation is
F a f ( u ) = K α ( u , x ) f ( x ) d x ,
where the term K α ( u , x ) is expressed as
K α ( u , x ) = { A α exp { i π [ cot ( α ) ( x 2 + u 2 ) 2 csc ( α ) u x ] } ,   ( α n π ) δ ( u x )   ( α = 2 n π ) δ ( u + x )   ( α = ( 2 n ± 1 ) π ) ,
where A α is the amplitude factor defined as 1 i c o t ( α ) . If = π / 2 , the FrFT is equivalent to the Fourier transform. Next, the optimal transform order in a discrete region can be expressed as
α o p t = 2 π = 2 π c o t 1 ( M k f s 2 ) ,
where f s is the sampling frequency, k is the chirp rate, and M represents the sampling point. If the sampling points and sampling frequency are fixed, the optimal transform order is determined from the value k .
Figure 5a shows a conceptual diagram of the FrFT for a chirp signal. In this case, we can assume that the chirp signal x u ( t ) with the up-slope characteristic exists for the analysis period T , and the chirp signal x u ( t ) has the initial frequency f 0 at the start time ( k τ T / 2 ) and the frequency f 1 at the end time ( k τ + T / 2 ) (generally, the frequency f 1 is higher than the frequency f 0 ). Therefore, the point o u which is the middle point between p u = [ ( k τ T / 2 ) / s ,   f 0 s ] and point q u = [ ( k τ + T / 2 ) / s ,   f 1 s ] can be expressed as [ k τ / s ,   f m s ], where f m is the center frequency and s is the rescaling factor, and is defined as s = T / f s . The k is the symbol number. The vector u 0 , which is the point o u that corresponds to the u -axis, is
u 0 = | o u | sin ( α o p t π 2 ) = [ f m s k τ s ] sin ( α o p t π 2 ) .
When the starting point for the signal x u ( t ) is τ = 0 , the vector u 0 can be written as:
u 0 = ( f m s ) sin ( α o p t π 2 ) .
When considered in the frequency domain, the signal also has a negative frequency that is inversely proportional to the positive frequency. The negative frequency exists symmetrically with the positive frequency and forms a relatively wide and low spectrum on the FrFT as [ p u   t o   q u ] in Figure 5a. On the other hand, in Figure 5b, the chirp signal x d ( t ) with the down-slope characteristic will appear symmetrically with the signal x u ( t ) . Therefore, the vector u 1 can be written as
u 1 = ( f m s ) sin ( α o p t π 2 ) ,
where f m is the center frequency of the signal x d ( t ) . If the signal x u ( t ) and x d ( t ) are in the same frequency band and have the same slope, the vector u 1 can be rewritten as
u 1 = ( f m s ) sin ( α o p t π 2 ) = u 0 .
The opposite characteristics of FrFT can be seen from Equation (13). If the center frequency and the frequency variation are same, the positions of the vectors u 0 and u 1 are the opposite and the same value.
Figure 6a shows the symmetrical characteristics in the fractional domain to distinguish between the up-slope and the down-slope. This is an example of the FrFT spectrum for Figure 5a. In the proposed method, with respect to the middle line of the fractional spectrum, the left side is the positive frequency domain, and the right side is the negative frequency domain (theoretically existing). Therefore, in the up-slope state (like Figure 5a), the energies are concentrated on the left side, and in the down-slope state (like Figure 5b), the energies are concentrated on the right side. On the other hand, Figure 6b shows the FrFT spectrum of the conventional method that coexists with the up-slope and the down-slope in the symbol duration. Unlike the proposed method, the conventional method is processed only in the positive frequency domain and distinguishes between the up-slope and the down-slope by using the two transform orders which are dependent on each symbol. This means that the conventional receiver must have the parallel structures. In this paper, we propose the FrFT receiver using only one transform order.
In the proposed method, when the frequency-hopping interval is f w d , the vector u h n regarding the hopped frequency f h n is written as
f h n = f m + n f w d .
Then,
        u h n = ( f h n s ) sin ( α o p t π 2 ) = [ ( f m + n f w d ) s ] sin ( α o p t π 2 ) = u 0 + ( n f w d s ) sin ( α o p t π 2 ) ,
where n represents the number of hopped intervals. The signals hop due to f w d , which means that they are spaced at regular intervals in the FrFT spectrum. Using the above property, regardless of the frequency band of the hopping, a ‘decision section’ for estimating a symbol can be formed because the FrFT bin formed for each symbol is always regular. That is, the FrFT bin is obtained for the received signal through Equation (15), and the symbol is estimated by the demapping method which distinguishes the FrFT bin in the decision section. Since the estimated symbol can be implemented as an ideal form of the signal, it is possible to perform a matched filter between the received signal and the ideal signal. The synchronization error can be obtained between the results of the above procedure and the results of the matched filter in the ideal case. When compensating for the next symbol along with the weight, it is possible to implement symbol synchronization to correct the errors associated with the Doppler shift.
Figure 7 is a plot of the FrFT spectrum for all frequency-hopping bands in a single space. From the relation given in Equation (15), it can be seen that the spectrum is formed regularly in proportion to the frequency-hopping interval. Based on this, even if an arbitrary frequency hopping is performed, it will be completed in the above spectrum section. That is, in the case that any hopped symbol is input, it can be estimated through the section of the spectrum.

4. Simulation and Results

We performed the simulation to analyze the performance of the proposed method. Figure 8 shows a packet structure of a transmission signal designed for the simulation. We used 1080 bits per frame, a center frequency of 16 kHz, a sampling frequency of 192 kHz, and a data rate of 100 bps. One chirp symbol band is 1.2 kHz and the total band is approximately 10 kHz. The LFM (linear frequency modulation) signals are used to distinguish the packets. The signal length is 0.1 s and the bandwidth is 4 kHz. A guard interval of 1 s is arranged to prevent an interference due to the signal. The pseudo-noise code uses 512 bits for 0.512 s, and no channel coding technique is applied.
There are several software packages for simulating channel transfer characteristics in underwater environments. Norwegian FFI provides the Watermark software that simulates transmission information in several underwater channels [31]. However, to verify the performance of the proposed method, a channel is simulated by applying the sound velocity distribution measured in the South Sea, South Korea, to the VirTEX (Virtual Time Series Experiment), which is a Bellhop-based underwater channel modeling program developed and released by the Scripps Institution of Oceanography (SIO) [32,33]. The simulation environment is shown in Figure 9. In the simulation, it is assumed that the sea floor is inclined at a depth of 30 m from the transmitter to 40 m from the receiver. The position of the receiver is placed at a depth of 30 m and the transmitter is assumed to be moving horizontally with a fixed depth of 5 m. It is assumed that the distance between the transmitter and receiver moves from 300 m to 200 m linearly at approximately 12 knots towards the receiver and then moves in the opposite direction at the same speed from 200 m to 300 m.
The channel response characteristics for the simulation are shown in Figure 10. The channel impulse response shows that not only the direct path, but also the sub-path, are reflected from the surface and the floor.
Figure 11a shows the result of the matched filter for the termination by overlapping the received signal and the ideal signal when approaching between the transmitter and the receiver. As shown in Figure 11a, the signal is compressed since it is positioned backward, even though the packet synchronization is initially adjusted due to the Doppler frequency, which is caused by the mobility of the transmitter. The synchronization error due to the compression is approximately 0.056 s, and assuming a transmission rate of 100 bps, the error is approximately 5–6 symbols.
On the other hand, Figure 11b shows an overlapping of the received signal and the ideal signal in the case of retreating, which expands the signal due to the Doppler frequency. The synchronization error due to the expansion is approximately 0.055 s and has an error of approximately 5 to 6 symbols, as in the previous case.
Table 1 shows the bit error rate according to the signal-to-noise ratio (SNR) in the presence or absence of the symbol synchronization. The symbol length changes continuously due to the compression and expansion of the signal. Accordingly, if the synchronization is not performed, a large error occurs. On the other hand, when the symbol synchronization is performed using the proposed method, the uncoded bit error rate (BER) is considerably low. The simulation results show that the proposed method is robust to multipath delays and variability due to the Doppler shift.
The second simulation was performed using the underwater channel software released by FFI. Watermark is a benchmark for physical-layer schemes of underwater acoustic communications. It contains data with five different channel characteristics, and we used a channel (named as NCS1) measured in Norway. This is the channel in which broad Doppler spreading occurs and repeats for 60 cycles. We set the Doppler shift higher to make more visible the performance difference of the proposed method. Other detailed information is provided in the reference [31]. Figure 12 shows the bit error rate. The horizontal axis of the figure refers to the number of cycles. The bit error rate of the proposed method is approximately 0.07, which shows better performance than the conventional method.

5. Lake Trials

To investigate the performance of the proposed method in a real environment, lake trials were conducted at Kyungcheon Lake, Mungyeong City. The first and second experiments were conducted at the beginning and end of May 2018, respectively. The water depth of the lake was approximately 30–40 m, and the transmitters were located approximately 5 m from the boat. The receiver was located at a depth of 30 m, and the boat was steered randomly from 100 m to approximately 400 m from the receiver. The boat moved at speeds of 2 knots and 4 knots, and the configuration of the signals used was the same as those of the previous simulations. However, in the second experiment, which was conducted at the end of May, the hopping rate designed a half-spacing signal for each hopping frequency. Because of the overlap of the hopping frequencies, the ISI can degrade the communication performance due to multipath delay. However, this can be overcome by using the FrFT and bandpass filtering. Additionally, by designing with the above method, the frequency band can be used more efficiently. A Neptune D/17/BB model projector was used as the transmitter and a B&K 8106 hydrophone was used as the receiver. The Neptune D/17/BB model that was used as a transmitter is a broadband projector with a bandwidth of approximately 10 kHz [34,35].
Figure 13 shows the channel impulse response and Figure 14 shows the scattering function used for the Doppler estimation. Figure 13a shows the measured channel response in the first experiment performed at the beginning of May, and Figure 13b is the result of the second experiment performed at the end of May. As shown in Figure 13, the difference in the reception time occurs as the transmitter moves, which generates a slope. Figure 13b shows there are more complex delay profiles. Figure 14a shows that the Doppler frequency is approximately −16 Hz to −22 Hz, and Figure 14b shows that the Doppler frequency is approximately −10 Hz to −5 Hz, where the (−) sign means that the distance between the transceivers gradually increased. However, there seems to be an issue with the scattering function in Figure 14b, as a main Doppler shift cannot be detected and the multipath arrival is entirely smeared. This was a result of the increased level of noise due to minor equipment problems and because the channel conditions in the second experiment were relatively poor due to stronger winds.
Figure 15a and Figure 16a show spectrograms of the transmitted signal, and Figure 15b and Figure 16b show spectrograms of the received signal through the lake experiments. From the figures, we can see that there is considerable fading in the received signals due to multipath propagation delay.
Since the received signal experiences substantial fading as it passes through the underwater channel, the delay due to multipath propagation or the slope change due to the Doppler frequency should be considered. Therefore, it is necessary to calculate the range adjacent to the optimal transform order for the application of the FrFT scheme at the receiver. Figure 17 shows the accumulation of the spectra of the received signal that is obtained from the optimal transform order to the adjacent range. When the optimal transform order is α o p t = 0.996 , the spectrum is spread by the delay caused by multipath propagation in the range adjacent to the order. However, due to the property of the chirp signal, the multipath has little effect on the estimation. Figure 17a shows the correctly estimated symbol, which is on the right side in the FrFT bin and fixed on the number 240 of the FrFT bin point. We can determine from this symbol that the down-slope chirp signal has a center frequency of 1.2 kHz. On the other hand, Figure 17b shows the incorrectly estimated symbol. The receiver estimates to a different point. The transmitted signal has the transform order α o p t = 0.996 , but the received signal has the approximate value “1”. This means that the receiver cannot distinguish between left and right. Therefore, we determined this symbol as the correct center frequency and the incorrect slope. In this case, the received symbol was affected by the serious Doppler effect beyond the Doppler tolerance of the chirp signal.
Figure 18 shows the cumulative error rate per elapsed symbol to analyze the synchronizing effect of the proposed method. If the symbol synchronization is not performed, the cumulative errors increase because the symbol is compressed or expanded over time. Conversely, the performance of the synchronized symbol is good. As a result, the FrFT receiver can be a system used to efficiently demodulate chirp signals that are robust and resistant to fading and Doppler effects.
Table 2 and Table 3 show the results of the uncoded bit error rate obtained by demodulating the received signal from the lake experiment, where the transmitted signal is not applied to the channel coding. Table 2 and Table 3 show the cases where the moving speed of the transmitter is approximately 2 knots and 4 knots, respectively. Some experiments were performed in an environment where the weather condition was relatively poor, and the corresponding performance may be lower than those of other cases. Nevertheless, according to the proposed method, the bit error rate is reduced to 36.4% after iterative symbol synchronization. Therefore, we verified that the proposed method improves the performance. The chirp signal has properties that are robust to fading and Doppler effects, and the proposed method showed that it was possible to maintain these advantages through efficient demodulation. To compare the computations, we compared the computation times for each of the arbitrary symbols designed with the proposed method and the conventional method. The parameters of symbols were set to be similar to those shown in Figure 6. In this case, by using the MATLAB function “Tic − Toc”, the running program time was measured and compared. The proposed method used a single receiver and scanned the value around the transform order to obtain the spectrum (see Figure 17). On the other hand, the conventional method used at least two receivers, which means that two transform orders must be calculated. As a result of a relative computation time comparison with a computer equipped with an Intel Core i7-6700/4 GHz processor, the proposed method took 0.2277 s while the conventional method took 0.4717 s.

6. Conclusions

In this paper, we proposed a covert acoustic communication scheme that is robust to fading and Doppler effects. It consists of a single chirp-based FHSS that is resistant to delay propagation and the Doppler shift. Additionally, the receiver demodulates to the FrFT spectrum. To verify the performance of the proposed method, we estimated the signal in an environment with a multipath delay, and the Doppler frequency through simulations and lake experiments, to determine the uncoded BER. Since the chirp signal has a relatively good orthogonality, it is robust against the instability of the phase. It was also possible to estimate the slope using the FrFT spectrum, even in the case of frequency variation due to the Doppler shift. This means that the proposed method can be demodulated while maintaining the properties of the chirp signal. As a result, we could obtain high-reliability communication performance. In the conventional method, the up-sloped chirp signal and the down-sloped chirp signal existed within one symbol simultaneously, and independent demodulators with different orders were used. However, in the proposed method, either the up-sloped or down-sloped chirp signal was transmitted within one symbol, and a parallel structure, as used with the conventional method, was not required in the demodulator. To verify the performance of the proposed method, we performed various simulations. Specifically, we used the Watermark software released by FFI. In the future, it will be necessary to conduct research for efficient FrFT spectrum formation and optimization of the receiver structure. Sea trials will be performed.

Author Contributions

All authors contributed significantly to the work presented in this manuscript. G.L. and K.K. proposed the method described in this paper; G.L., T.K. and K.K. conceived and designed the experiments; G.L., T.K., W.-J.K. and W.P. performed the experiments and analyzed the data; W.-J.K. provided some valuable advice.

Funding

This work was supported by the Agency for Defense Development, South Korea under Grant UD170022DD and Low Frequency Underwater Research Laboratory.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Catipovic, J.A. Performance limitations in underwater acoustic telemetry. IEEE J. Ocean. Eng. 1990, 15, 205–216. [Google Scholar] [CrossRef]
  2. Yaacoub, E.; Dawy, Z. A survey on uplink resource allocation in OFDMA wireless networks. IEEE Commun. Surv. Tutor. 2012, 14, 322–337. [Google Scholar] [CrossRef]
  3. Tsiropoulou, E.E.; Kapoukakis, A.; Papavassiliou, S. Uplink resource allocation in SC-FDMA wireless networks: A survey and taxonomy. Comput. Netw. 2016, 96, 1–28. [Google Scholar] [CrossRef]
  4. Stojanovic, M. Underwater acoustic communications. In Proceedings of the IEEE Electro/International, Boston, MA, USA, 21–23 June 1995; pp. 435–440. [Google Scholar]
  5. Stojanovic, M.; Catipovic, J.; Proakis, J. Phase coherent digital communications for underwater acoustic channels. IEEE J. Ocean. Eng. 1994, OE-16, 100–111. [Google Scholar] [CrossRef]
  6. Kaushal, H.; Kaddoum, G. Underwater optical wireless communication. IEEE J. Mag. 2016, 4, 1518–1547. [Google Scholar] [CrossRef]
  7. ElMoslimany, A.; Zhou, M.; Duman, T.M.; Suppapplola, A.P. A new signaling scheme for underwater acoustic. In Proceedings of the IEEE Computer Society Conference on OCEANS 2013 MTS/IEEE, San Diego, CA, USA, 23–27 September 2013; pp. 1–5. [Google Scholar]
  8. Liu, S.; Qiao, G.; Ismail, A.; Liu, B.; Zhang, L. Covert underwater acoustic communication using whale noise masking on DSSS signal. In Proceedings of the OCEANS 2013 MTS/IEEE Bergen: The Challenges of the Northern Dimension, Bergen, Norway, 10–14 June 2013; pp. 1–6. [Google Scholar]
  9. Stojanovic, M.; Proakis, J.G.; Rice, J.A.; Green, M.D. Spread spectrum underwater acoustic telemetry. Proc. MTS/OCEANS 1998, 2, 650–654. [Google Scholar]
  10. Yang, T.C.; Yang, W.B. Performance analysis of direct-sequence spread-spectrum underwater acoustic communications with low signal-to-noise-ratio input signals. J. Acoust. Soc. Am. 2008, 2, 842–855. [Google Scholar] [CrossRef] [PubMed]
  11. Yang, T.C.; Yang, W.B. Low probability of detection underwater acoustic communications using direct-sequence spread spectrum. J. Acoust. Soc. Am. 2008, 6, 3632–3647. [Google Scholar] [CrossRef]
  12. Ling, J.; He, H.; Li, J.; Roberts, W.; Stoica, P. Covert underwater acoustic communications. J. Acoust. Soc. Am. 2010, 128, 2898–2909. [Google Scholar] [CrossRef] [PubMed]
  13. Freitag, L.; Stojanovic, M.; Singh, S.; Johnson, M. Analysis of channel effects on direct-sequence and frequency-hopped spread-spectrum acoustic communication. IEEE J. Ocean. Eng. 2001, 26, 586–593. [Google Scholar] [CrossRef] [Green Version]
  14. Kaleh, G.K. Performance comparison of frequency-diversity and frequency-hopping spread-spectrum systems. IEEE Trans. Commun. 1997, 45, 910–912. [Google Scholar] [CrossRef]
  15. Yang, W.B.; Ynag, T.C. High-frequency FH-FSK underwater acoustic communications: The environmental effect and signal processing. Proc. High Freq. Ocean Acoust. Conf. 2004, 728, 106–113. [Google Scholar]
  16. Potter, J.; Alves, J.; Green, D.; Zappa, G.; Nissen, I.; McCoy, K. The JANUS underwater communications standard. In Proceedings of the 2014 Underwater Communications and Networking (UComms), Sestri Levante, Italy, 3–5 September 2014. [Google Scholar]
  17. Walree, P.; Ludwig, T.; Solberg, C.; Sangfelt, E.; Laine, A.; Bertolotto, G.; Ishoy, A. UUV Covert Acoustic Communications. Available online: https://pdfs.semanticscholar.org/ee86/f2f5e30a241831f7cdc32e697a541f5f1c83.pdf (accessed on 28 September 2018).
  18. Marchetti, L.; Reggiannini, R. An efficient receiver structure for sweep-spread-carrier underwater acoustic links. IEEE J. Ocean. Eng. 2016, 41, 440–449. [Google Scholar]
  19. Cook, C.E. Linear FM signal formats for beacon and communication systems. IEEE Trans. Aerosp. Electr. Syst. 1974, AES-10, 471–478. [Google Scholar] [CrossRef]
  20. Ozaktas, H.M.; Arikan, O.; Kutay, M.A.; Bozdagi, G. Digital computation of the fractional Fourier transform. IEEE Trans. Signal Process. 1996, 44, 2141–2150. [Google Scholar] [CrossRef] [Green Version]
  21. Almeida, L.B. The fractional Fourier transform and time-frequency representations. IEEE Trans. Signal Process. 1994, 42, 3084–9091. [Google Scholar] [CrossRef]
  22. Gaglione, D.; Clemente, C.; Llioudis, C.V.; Persico, A.R.; Proudler, I.K.; Soraghan, J.J. Fractional Fourier based waveform for a joint radar-communication system. In Proceedings of the IEEE Radar Conference, Philadelphia, PA, USA, 2–6 May 2016. [Google Scholar]
  23. Qiu, X.; Sha, X.J.; Shen, Q.H.; Wu, X.L. Fractional Fourier transform based secure communication system. In Proceedings of the 2010 6th International Conference on Wireless Communications Networking and Mobile Computing, Chengdu, China, 23–25 September 2010. [Google Scholar]
  24. Yuan, F.; Wei, Q.; Cheng, E. Multiuser chirp modulation for underwater acoustic channel based on VTRM. Int. J. Navel Archit. Ocean Eng. 2017, 9, 256–265. [Google Scholar] [CrossRef]
  25. Chen, Y.; Clemente, C.; Soraghan, J.J.; Weiss, S. Partial fractional Fourier transform (PFRFT)-OFDM for underwater acoustic communication. In Proceedings of the 23rd European Signal Processing Conference, Nice, France, 31 August–4 September 2015. [Google Scholar]
  26. Tu, X.; Xu, X.; Zou, Z.; Yang, L.; Wu, J. Fractional Fourier domain hopped communication method based on chirp modulation for underwater acoustic channels. J. Syst. Eng. Electr. 2017, 28, 449–456. [Google Scholar]
  27. Sharif, B.S.; Neasham, J.; Hinton, O.R.; Adams, A.E. A computationally efficient Doppler compensation system for underwater acoustic communications. IEEE J. Ocean. Eng. 2000, 25, 52–61. [Google Scholar] [CrossRef]
  28. Kaminsky, E.J.; Simanjuntak, L. Chirp slope keying for underwater communications. Proc. SPIE Prog. Biomed. Opt. Imaging 2005, 5778, 894–905. [Google Scholar] [Green Version]
  29. Plumb, R.G. Matched filter response of a linear array with time-varying phase weights. IEEE Trans. Aerosp. Electr. Syst. 1993, 29, 1046–1050. [Google Scholar] [CrossRef]
  30. Kaihan, Z.; Keyu, C.; Fei, Y.; Jinwang, Y. Chirp FSK based on FRFT for underwater acoustic communication. In Proceedings of the 7th IEEE International Conference on Electronics Information and Emergency Communication (ICEIEC), Macau, China, 21–23 July 2017; pp. 49–52. [Google Scholar]
  31. FFI. Available online: https://www.ffi.no/en/research-projects/Sider/Watermark.aspx (accessed on 28 September 2018).
  32. Kim, J.S.; Song, H.C.; Hodgkiss, W.S. Virtual time series experiment (VIRTEX) simulation tool for underwater acoustic communications. J. Acoust. Soc. Am. 2009, 10, 2174–3647. [Google Scholar] [CrossRef]
  33. Siderius, M.; Oirter, M.B. Modeling broadband ocean acoustic transmissions with time-varying sea surface. J. Acoust. Soc. Am. 2008, 124, 137–150. [Google Scholar] [CrossRef] [PubMed]
  34. NEPTUNE SONAR. Available online: http://www.neptune-sonar.co.uk/wp-content/uploads/2016/03/Model-D17.pdf (accessed on 1 August 2018).
  35. Bruel & Kjaer. Available online: https://www.bksv.com/en/products/transducers/acoustic/microphones/hydrophones/8106 (accessed on 1 August 2018).
Figure 1. Basic scheme of the chirp-based FHSS.
Figure 1. Basic scheme of the chirp-based FHSS.
Sensors 18 04498 g001
Figure 2. Block diagram of the proposed method.
Figure 2. Block diagram of the proposed method.
Sensors 18 04498 g002
Figure 3. Matched filter output in the same frequency band. (a) The same slope symbol; (b) different slope symbols.
Figure 3. Matched filter output in the same frequency band. (a) The same slope symbol; (b) different slope symbols.
Sensors 18 04498 g003aSensors 18 04498 g003b
Figure 4. Matched filter output in different frequency bands. (a) The same slope symbol; (b) different slope symbols.
Figure 4. Matched filter output in different frequency bands. (a) The same slope symbol; (b) different slope symbols.
Sensors 18 04498 g004
Figure 5. Relation between the chirp signal and FrFT coefficient α . (a) The up-slope chirp signal; (b) the down-slope chirp signal.
Figure 5. Relation between the chirp signal and FrFT coefficient α . (a) The up-slope chirp signal; (b) the down-slope chirp signal.
Sensors 18 04498 g005
Figure 6. FrFT spectrum comparison between (a) the proposed method, and (b) the conventional method.
Figure 6. FrFT spectrum comparison between (a) the proposed method, and (b) the conventional method.
Sensors 18 04498 g006
Figure 7. Aggregate FrFT spectra for all frequency bins.
Figure 7. Aggregate FrFT spectra for all frequency bins.
Sensors 18 04498 g007
Figure 8. Packet design.
Figure 8. Packet design.
Sensors 18 04498 g008
Figure 9. Sound velocity profile for the simulation.
Figure 9. Sound velocity profile for the simulation.
Sensors 18 04498 g009
Figure 10. Channel impulse response for simulation.
Figure 10. Channel impulse response for simulation.
Sensors 18 04498 g010
Figure 11. Matched filtering output for the received data and the ideal data. (a) The approach; (b) the retreat.
Figure 11. Matched filtering output for the received data and the ideal data. (a) The approach; (b) the retreat.
Sensors 18 04498 g011
Figure 12. Simulation result using Watermark benchmark model.
Figure 12. Simulation result using Watermark benchmark model.
Sensors 18 04498 g012
Figure 13. Measured channel impulse responses for (a) the first lake trial and (b) the second lake trial.
Figure 13. Measured channel impulse responses for (a) the first lake trial and (b) the second lake trial.
Sensors 18 04498 g013
Figure 14. Scattering function for the Doppler frequency estimation for (a) the first lake trial and (b) the second lake trial.
Figure 14. Scattering function for the Doppler frequency estimation for (a) the first lake trial and (b) the second lake trial.
Sensors 18 04498 g014
Figure 15. Spectrogram of (a) the transmitted signal, and (b) the received signal in the first trial.
Figure 15. Spectrogram of (a) the transmitted signal, and (b) the received signal in the first trial.
Sensors 18 04498 g015
Figure 16. Spectrogram of (a) the transmitted signal, and (b) the received signal in the second trial.
Figure 16. Spectrogram of (a) the transmitted signal, and (b) the received signal in the second trial.
Sensors 18 04498 g016
Figure 17. Color intensity of the FrFT results in (a) the correctly estimated symbol and (b) the incorrectly estimated symbol.
Figure 17. Color intensity of the FrFT results in (a) the correctly estimated symbol and (b) the incorrectly estimated symbol.
Sensors 18 04498 g017
Figure 18. Synchronization performance.
Figure 18. Synchronization performance.
Sensors 18 04498 g018
Table 1. Uncoded bit error rate in the simulation.
Table 1. Uncoded bit error rate in the simulation.
SNR (dB)Non-Sync MethodSync Method
150.490.000
100.480.002
50.490.004
00.480.050
Table 2. Uncoded bit error rate at 2 knots.
Table 2. Uncoded bit error rate at 2 knots.
TrialFirst TrialSecond Trial
Non-Sync MethodSync MethodNon-Sync MethodSync Method
10.230.090.190.08
20.240.090.350.11
30.280.0010.100.09
40.120.090.080.07
50.220.110.100.01
60.090.010.070.06
Table 3. Uncoded bit error rate at 4 knots.
Table 3. Uncoded bit error rate at 4 knots.
TrialNon-Sync MethodSync Method
10.170.04
20.260.05

Share and Cite

MDPI and ACS Style

Lee, G.; Park, W.; Kang, T.; Kim, K.; Kim, W. Chirp-Based FHSS Receiver with Recursive Symbol Synchronization for Underwater Acoustic Communication. Sensors 2018, 18, 4498. https://doi.org/10.3390/s18124498

AMA Style

Lee G, Park W, Kang T, Kim K, Kim W. Chirp-Based FHSS Receiver with Recursive Symbol Synchronization for Underwater Acoustic Communication. Sensors. 2018; 18(12):4498. https://doi.org/10.3390/s18124498

Chicago/Turabian Style

Lee, Geunhyeok, Woongjin Park, Taewoong Kang, Kiman Kim, and Wanjin Kim. 2018. "Chirp-Based FHSS Receiver with Recursive Symbol Synchronization for Underwater Acoustic Communication" Sensors 18, no. 12: 4498. https://doi.org/10.3390/s18124498

APA Style

Lee, G., Park, W., Kang, T., Kim, K., & Kim, W. (2018). Chirp-Based FHSS Receiver with Recursive Symbol Synchronization for Underwater Acoustic Communication. Sensors, 18(12), 4498. https://doi.org/10.3390/s18124498

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