sensors
Article Biomimicking Covert Communication by Time-Frequency Shift Modulation for Increasing Mimicking and BER Performances
Jongmin Ahn 1 , Hojun Lee 1 , Yongcheol Kim 1, Wanjin Kim 2 and Jaehak Chung 1,*
1 Department of Electronics Engineering, INHA University, Inhceon 22201, Korea; [email protected] (J.A.); [email protected] (H.L.); [email protected] (Y.K.) 2 Agency of Defense Development, Changwon-si 51682, Korea; [email protected] * Correspondence: [email protected]; Tel.: +82-032-860-7421
Abstract: Underwater acoustic (UWA) biomimicking communications have been developed for covert communications. For the UWA covert communications, it is difficult to achieve the bit error rate (BER) and the degree of mimic (DoM) performances at the same time. This paper proposes a biomimicking covert communication method to increase both BER and DoM (degree of mimic) perfor- mances based on the Time Frequency Shift Keying (TFSK). To increase DoM and BER performances, the orthogonality requirements of the time- and frequency-shifting units of the TFSK are theoretically derived, and the whistles are multiplied by the sequence with a large correlation. Two-step DoM assessments are also developed for the long-term whistle signals. Computer simulations and practical lake and ocean experiments demonstrate that the proposed method increases the DoM by 35% and attains a zero BER at −6 dB of Signal to Noise Ratio (SNR). Keywords: convert communication; bio-mimetic; signal processing; underwater communication; Citation: Ahn, J.; Lee, H.; Kim, Y.; MOS test Kim, W.; Chung, J. Biomimicking Covert Communication by Time-Frequency Shift Modulation for Increasing Mimicking and BER 1. Introduction Performances. Sensors 2021, 21, 2184. Underwater acoustic (UWA) covert communication requires covertness and confiden- https://doi.org/10.3390/s21062184 tiality. The conventional UWA covert communication schemes achieve the covertness using the spread spectrum that spreads out the communication signal over a wide frequency Academic Editor: Theodore band, which is considered as background noise [1–7]. However, the narrow bandwidth of E. Matikas the UWA communication cannot utilize a large spreading factor, and a low transmit power for the covertness causes a short available communication range [1–7]. As an alternate, Received: 5 February 2021 biomimetic communications have been developed, which mimic the dolphin whistles to Accepted: 17 March 2021 Published: 20 March 2021 increase the covertness [8–13]. The idea of biomimetic communications is to make the enemy confuse the communication signals with dolphin sounds [8–14].
Publisher’s Note: MDPI stays neutral The CV-CFM (continuously varying-carrier frequency modulation) method was de- with regard to jurisdictional claims in veloped for biomimetic communication [12]. The CV-CFM selects one whistle among published maps and institutional affil- many whistles, and divides the whistle into many short-time periods, and modulates iations. binary bits to the periods using the conventional digital communication modulations (i.e., Chirp Spread Spectrum (CSS), Frequency Shift Keying (FSK), and Phase Shift Keying (PSK)) [9,10,12]. The short division of the whistle increased the frequency bandwidth of the whistle pattern, which decreased the degree of mimic (DoM) and Bit Error Ratio (BER) performances. To increase the DoM and BER performances, the time-frequency shift Copyright: © 2021 by the authors. keying (TFSK) method was researched [13]. The TFSK method did not allocate the binary Licensee MDPI, Basel, Switzerland. This article is an open access article information bits to the divided whistles, but to the shifted locations in the time-frequency of distributed under the terms and the whistles [13]. However, the BER performance of the TFSK varied with the transmitted conditions of the Creative Commons whistle patterns, and the time- and frequency-shifting units that satisfy the orthogonality Attribution (CC BY) license (https:// in the time- and frequency-domains were not derived. If a few whistle patterns with a creativecommons.org/licenses/by/ low BER are used for a long data sequence, the DoM of Ref. [13] decreases by the repeated 4.0/). similar whistles. For the practical biomimetic covert communication, the large DoM and
Sensors 2021, 21, 2184. https://doi.org/10.3390/s21062184 https://www.mdpi.com/journal/sensors Sensors 2021, 21, 2184 2 of 18
the low BER need to be simultaneously achieved and to be evaluated according to vari- ous parameters of the real dolphin sounds, e.g., long-term signal duration, bandwidth, etc. [8–14]. In general, the time duration and the frequency bandwidth of the dolphin whistles vary from several hundred milliseconds to two seconds and from several hundred Hz to tens of kHz, respectively [15–19]. Thus, the DoM evaluation needs to be performed for a long-term period. This paper proposes a TFSK-based biomimetic communication method increasing the DoM and BER performances. We theoretically derive the orthogonality requirements in the time- and frequency-domains and utilizes various whistle patterns with sequences with the large autocorrelation to increases the DoM and BER performances. To evaluate the long-term DoM of the proposed method, a two-step DoM assessment is developed. The computer simulation and the ocean experiments were executed to show that the proposed method demonstrated better DoM and BER performances than the conventional convert communications. The main contributions of the paper are summarized as follows: 1. The time- and frequency-shifting unit requirements of the TFSK are theoretically derived for the orthogonality in the time- and frequency-domains. The requirements guarantee the low BER. 2. For the large DoM and the low BER, the sequence with a large correlation is multi- plied to the whistles, which enables to use of various dolphin whistles without any restrictions. In addition, the orthogonality requirements for the proposed method are also derived. 3. Since the sequence makes the whistle spread in the frequency domain, the DoM assessments are conducted to find the unrecognizable spreading parameter. Thus, a two-step DoM assessment is proposed: The 1st step is conducted to find the best length of the sequence. The 2nd step is executed to confirm whether the long-term whistles signal with the selected sequence length is acceptable for the covert commu- nication. 4. The computer simulations and the practical lake and ocean experiments were con- ducted and demonstrated the proposed method had the large DoM and the lower BER compared with the conventional covert communication methods. This paper is organized as follows. In Section2, the orthogonal requirement of the TFSK according to a whistle pattern is derived, and the proposed biomimetic communica- tion method is described. Section3 explains the DoM assessment method for the proposed biomimetic communication signals. In Section4, the BER and DoM performances of the proposed method are analyzed. In Section5, the proposed biomimetic communication method demonstrated the large DoM and the low BER through computer simulations and practical lake and ocean experiments. Section6 concludes the paper.
2. Proposed Method The conventional TFSK scheme modulates binary bits to the shifted time-frequency position of the whistles [13]. Ref. [13] does not provide the requirements of the time- and frequency-shifting units which guarantee the orthogonality of whistles, and if zero-chirp rate whistles are transmitted, the detection performance degrades by the time ambiguity. Thus, few selected whistles may be utilized to attain the large BER. Figure1 shows the spectrograms of the real whistles. In Figure1, the zero-chirp rate whistles are frequently found. If these zero-chirp rate whistles are not transmitted to avoid the detection ambiguity of the TFSK, the DoM of the TFSK decreases for the long-term observation. Sensors 2021, 21, x FOR PEER REVIEW 3 of 19
in the time- and frequency-domains and utilizes various whistle patterns with sequences with the large autocorrelation to increases the DoM and BER performances. To evaluate the long-term DoM of the proposed method, a two-step DoM assessment is developed. The computer simulation and the ocean experiments were executed to show that the pro- posed method demonstrated better DoM and BER performances than the conventional convert communications. The main contributions of the paper are summarized as follows: 1. The time- and frequency-shifting unit requirements of the TFSK are theoretically de- rived for the orthogonality in the time- and frequency-domains. The requirements guarantee the low BER. 2. For the large DoM and the low BER, the sequence with a large correlation is multi- plied to the whistles, which enables to use of various dolphin whistles without any restrictions. In addition, the orthogonality requirements for the proposed method are also derived. 3. Since the sequence makes the whistle spread in the frequency domain, the DoM as- sessments are conducted to find the unrecognizable spreading parameter. Thus, a two-step DoM assessment is proposed: The 1st step is conducted to find the best length of the sequence. The 2nd step is executed to confirm whether the long-term whistles signal with the selected sequence length is acceptable for the covert commu- nication. 4. The computer simulations and the practical lake and ocean experiments were con- ducted and demonstrated the proposed method had the large DoM and the lower BER compared with the conventional covert communication methods. This paper is organized as follows. In Section 2, the orthogonal requirement of the TFSK according to a whistle pattern is derived, and the proposed biomimetic communi- cation method is described. Section 3 explains the DoM assessment method for the pro- posed biomimetic communication signals. In Section 4, the BER and DoM performances of the proposed method are analyzed. In Section 5, the proposed biomimetic communica- tion method demonstrated the large DoM and the low BER through computer simulations and practical lake and ocean experiments. Section 6 concludes the paper.
2. Proposed Method The conventional TFSK scheme modulates binary bits to the shifted time-frequency position of the whistles [13]. Ref. [13] does not provide the requirements of the time- and frequency-shifting units which guarantee the orthogonality of whistles, and if zero-chirp rate whistles are transmitted, the detection performance degrades by the time ambiguity. Thus, few selected whistles may be utilized to attain the large BER. Figure 1 shows the spectrograms of the real whistles. In Figure 1, the zero-chirp rate whistles are frequently Sensors 2021, 21, 2184 3 of 18 found. If these zero-chirp rate whistles are not transmitted to avoid the detection ambigu- ity of the TFSK, the DoM of the TFSK decreases for the long-term observation.
Sensors 2021, 21, x FOR PEER REVIEW 4 of 19
Figure 1. Spectrograms of the various realreal Delphinus delphisdelphis whistles.whistles.
We theoreticallyWe derive theoretically the orthogonality derive the orthogonalityrequirements requirementsin time- and frequency-shift- in time- and frequency-shifting ing units and proposeunits and a novel propose scheme a novel that scheme multiplies that the multiplies sequence the with sequence the large with autocor- the large autocorrela- relation to the tionwhistles to the to whistles increase to the increase BER performance the BER performance in time. inThus, time. all Thus, whistles all whistles are are utilized utilized for the fortransmission the transmission whistles, whistles, which whichincreases increases the DoM. the DoM.
2.1. TFSK Performance Analysis According to Whistle Pattern 2.1. TFSK Performance Analysis According to Whistle Pattern For the derivation of the time- and frequency-shifting unit requirements, the whistles For the derivationare mathematically of the time- and modeled, frequenc andy-shifting the orthogonality unit requirements, requirements the whistles are derived. Assume are mathematically modeled, and the orthogonality requirements are derived. Assume that the whistle pattern is modeled as a function ( fw(t)), which is presented as [12]: that the whistle pattern is modeled as a function ( ( )), which is presented as [12]: Z w(t) = cos f (t)dt . (1) ( ) = ( ) . w (1)
Since the conventionalSince the TFSK conventional modulation TFSK sc modulationheme shifts scheme the whistle shifts thepatterns whistle by patterns the by the time time and the frequencyand the frequency according according to transmit to transmitthe binary the bits, binary let the bits, time- let the and time- frequency- and frequency-shifting shifting units beunits ∆ beand∆t ∆ and, respectively,∆ f , respectively, and and the the total total number number of of time time and and frequencyfrequency grids be M and grids be andN , respectively., respectively. If theIf the whistle whistle is modulatedis modulated with with input input bits, bits, a reference a reference whistle is shifted M− M− N− N− whistle is shiftedfrom from−∆t −∆ 1 to ∆ tto ∆ 1 in the in time the domain time domain and from and−∆ fromf 1 −∆ to ∆ f 1toin the frequency 2 2 2 2 domain, respectively. Gray coding to the time-frequency mapping may be utilized to ∆ in the frequency domain, respectively. Gray coding to the time-frequency map- increase the BER. For M and N grids, the TFSK modulation conveys ( log2(N) + log2(M)) ping may be utilized to increase the BER. For and grids, the TFSK modulation con- bits per one whistle. ( ( ) + ( )) veys If the whistle bits per is modulatedone whistle. with m and n grid, which is an integer less than M and M, If the whistle is modulated with and grid, which is an integer less than and respectively, the TFSK modulated whistle (STF(t)) is expressed as: , respectively, the TFSK modulated whistle ( ( )) is expressed as: j2πn∆ f t STF(t) = [δ(t − m∆t ∆ ) ⊗ w(t)] × e , (2) ( ) = [ ( − ∆ ) ⊗ ( )] × , (2) where ⊗ denoteswhere a convolutional⊗ denotes a operation. convolutional Figure operation. 2 shows the Figure two different2 shows TFSK the two mod- different TFSK ulated whistle examples.modulated whistle examples.
01 bit modulated 11 bit modulated 01 bit modulated 11 bit modulated whistle whistle Ambiguity whistle whistle 10 10 ∆ ∆ 11 Frequency Frequency ∆ ∆ 11 01 0 0 01 00 00 00 01 11 10 Bits 00 01 11 10 Bits Time Time
(a) (b)
Figure 2. ExamplesFigure 2. ofExamples Time Frequency of Time Frequency Shift Keying Shift (TFSK) Keying modulated (TFSK) modulated signals: (a )signals: zero-chirp (a) zero-chirp rate whistle rate pattern (chirp whistle pattern (chirp rate: 0 Hz/s), (b) Up-chirp whistle pattern (chirp rate is larger than 0 Hz/s). rate: 0 Hz/s), (b) Up-chirp whistle pattern (chirp rate is larger than 0 Hz/s).
In Figure 2, Inand Figure are2, fourM and eachN andare foura total each of four and abits total are of allocated. four bits The are allocated.gray- The gray- colored whistlecolored denotes whistle a reference denotes whistle, a reference and blue whistle, and red and whistles blue and denote red whistles the TFSK denote the TFSK modulated ones.modulated In Figure ones. 2a, a Inzero-chirp Figure2a, ra ate zero-chirp whistle is ratemodulated. whistle isThe modulated. red and blue The red and blue whistles are modulated by 01 and 11 in time, respectively, and by 10 in frequency. In Fig- ure 2b, a whistle with a large-chirp rate is modulated. The red and blue ones are modu- lated by 01 and 11 in time, respectively, and by 10 in frequency. In Figure 2a, the blue and the red whistles are overlapped in the time-frequency domain because of the zero-chirp rate whistle and the small ∆ , and the receiver is unable to correctly determine the trans- mitted bits by the overlapped region. If ∆ is larger than that of the whistle length, two whistles are separated in the time domain and the detection of the whistle has no ambi- guity. Then, the receiver correctly decodes the transmitted bits. In Figure 2b, however, even though the same ∆ is utilized, two modulated whistles are orthogonal in the time- frequency domain, and the receiver has a low BER. Thus, the BER performance of the TFSK method depends on the whistle pattern, and to increase the BER, finding the good
Sensors 2021, 21, 2184 4 of 18
whistles are modulated by 01 and 11 in time, respectively, and by 10 in frequency. In Figure2b , a whistle with a large-chirp rate is modulated. The red and blue ones are modulated by 01 and 11 in time, respectively, and by 10 in frequency. In Figure2a, the blue and the red whistles are overlapped in the time-frequency domain because of the zero-chirp rate whistle and the small ∆t, and the receiver is unable to correctly determine the transmitted bits by the overlapped region. If ∆t is larger than that of the whistle length, two whistles are separated in the time domain and the detection of the whistle has no ambiguity. Then, the receiver correctly decodes the transmitted bits. In Figure2b , Sensors 2021, 21, x FOR PEER REVIEWhowever, even though the same ‘∆t is utilized, two modulated whistles are orthogonal5 of 19 in
the time-frequency domain, and the receiver has a low BER. Thus, the BER performance of the TFSK method depends on the whistle pattern, and to increase the BER, finding the ∆ good and∆ ∆ t and satisfying∆ f satisfying the orthogonal the orthogonal requirements requirements of the TFSK of the needs TFSK for needs the forgiven the whis- given tlewhistle pattern. pattern. For ∆ f , the calculation of the orthogonality requirement for an arbitrary whistle For ∆ , the calculation of the orthogonality requirement for an arbitrary whistle pat- pattern is difficult. However, if the receiver demodulates the received signal using the tern is difficult. However, if the receiver demodulates the received signal using the mul- multiplication of the complex conjugate to the transmitted whistle, the demodulated signal tiplication of the complex conjugate to the transmitted whistle, the demodulated signal has a zero-chirp rate pattern. Then, the orthogonality requirement of ∆ f becomes similar has a zero-chirp rate pattern. Then, the orthogonality requirement of ∆ becomes similar to that of the conventional FSK modulation. When the time length of the whistle is L , the to that of the conventional FSK modulation. When the time length of the whistle isw , ∆ f satisfying orthogonal requirements of the TFSK is given as: the ∆ satisfying orthogonal requirements of the TFSK is given as: ∆ f > 1/L . (3) ∆ >1/ w. (3)
ForFor ∆ ∆,t ,more more manipulation manipulation is is need needed.ed. Firstly, Firstly, assume assume that that ∆ ∆t isis lessless thanthan L w .. IfIf∆ ∆ t is islarger larger than than L w ,, no no overlapoverlap andand nono misdetection occur, occur, but but the the data data rate rate decreases. decreases. Thus, Thus, thisthis case case is is not not considered inin thisthis paper. paper. Since Since the the whistle whistle patterns patterns are are varied varied and and non-linear non- linearfunction, function, the derivation the derivation of the of orthogonality the orthogonality requirement requirement of ∆t for of the∆ forwhistles the whistles is difficult. is difficult.If the non-linear If the non-linear whistle patternwhistle is pattern divided is into divided short-intervals, into short-intervals, a piece of thea piece whistle of the can whistlesimply can be modeledsimply be as modeled a linear functionas a linear ( fw (functiont)), e.g., 1( f w(( t))),= e.g.,at + b(with ) = + a chirp with rate (a ), chirpwhich rate is defined( ), which as a =is ddefinedf /∆t where as =d /∆ ∆t and d f wheredenote ∆ the and time- d and denote frequency-differences, the time- and frequency-differences,respectively. Thus, if ∆respectively.t is derived fromThus, the if ∆ smallest is derived chirp-rate from ofthe the smallest whistle, chirp-rate the derived of the∆t whistle,satisfies the the derived orthogonality ∆ satisfies of the the whistle orthogonality with the length of the (whistleLw). with the length ( ). AssumeAssume that that θ is is thethe declinedecline of of the the whistle whistle in in the the time-frequency time-frequency domain domain and and defined de- − finedas θ as= tan = 1(d f /(d /∆ )∆t). To. calculateTo calculate the the orthogonality orthogonality requirement requirement of of the the whistle whistle with witha , we, we rotate rotate two two whistles whistles to toθ clockwise clockwise in thein the time-frequency time-frequency domain domain in Figurein Figure3. The 3. The two twomodulated modulated whistles whistles are are separated separated in the in time-frequencythe time-frequency domain domain keeping keeping the orthogonality. the orthog- onality.In Figure In 3Figurea, if one 3a, whistle if one whistle is shifted is shifted by ∆t inby the∆ time in the domain, time domain, the overlapped the overlapped time is ( − ) timegiven is given as Lw as ( ∆t and−∆ ) the and frequency the frequency difference difference is given is asgivend f Hz. as In Figure Hz. In3 b,Figure after 3b, the rotation by θ, the time duration of (L − ∆t)/cos θ is overlapped in the time domain and after the rotation by , the time durationw of ( −∆ )/ is overlapped in the time domainthe frequency and the gap frequency of cosθ ×gapd fofis obtained × in is the obtained frequency in the domain. frequency domain.
overlapped region