CN105227246B - It is a kind of to utilize the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals - Google Patents

It is a kind of to utilize the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals Download PDF

Info

Publication number
CN105227246B
CN105227246B CN201510676931.8A CN201510676931A CN105227246B CN 105227246 B CN105227246 B CN 105227246B CN 201510676931 A CN201510676931 A CN 201510676931A CN 105227246 B CN105227246 B CN 105227246B
Authority
CN
China
Prior art keywords
signal
mrow
lfm
msub
section
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510676931.8A
Other languages
Chinese (zh)
Other versions
CN105227246A (en
Inventor
刘凇佐
乔钢
马天龙
刘冰洁
聂东虎
马璐
尹艳玲
周锋
孙宗鑫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201510676931.8A priority Critical patent/CN105227246B/en
Publication of CN105227246A publication Critical patent/CN105227246A/en
Application granted granted Critical
Publication of CN105227246B publication Critical patent/CN105227246B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

The underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals is utilized the invention discloses a kind of.The characteristics of being changed over time according to dolphin whistle signal frequency, utilize multiple LFM signal subsections combine analog dolphin whistle signals, modulates information to LFM signals is realized using BOK methods, and then realizes the bionical underwater sound communication modulation based on dolphin whistle signal and LFM signals.The centre frequency and frequency modulation rate of every section of LFM signal in the reception data after equilibrium are estimated in receiving terminal, and the dolphin whistle signal sample used during based on modulation, it is demodulated, obtains modulation intelligence.The underwater acoustic communication method of dolphin whistle signal is imitated using segmentation LFM, it is possible to achieve the bionic camouflage underwater sound communication of high fidelity is carried out using simple linear FM signal, strengthen the disguise and practicality of communication.

Description

It is a kind of to utilize the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals
Technical field
The water that LFM signals imitate dolphin whistle signal is segmented the invention belongs to field of underwater acoustic communication, more particularly to a kind of utilize Acoustic communication method.
Background technology
In recent years, with the development of modern Detection Techniques, concealed requirement is proposed to underwater sound communication.
Traditional concealed communication method is signal of communication is hidden in ambient noise, i.e., hidden under Low SNR Communication.But it is due to time extension, frequency expansion and the serious attenuation characteristic of underwater acoustic channel, the signal of identical sound source level is not The covert communications in different distance being realized, especially when there is audiomonitor between transmitting-receiving node, will be unable to realize hidden lead to Letter.It is different from the hidden underwater sound communication of low signal-to-noise ratio, utilize the camouflage underwater sound communication side for being segmented the imitative dolphin whistle signal of LFM signals Method, is used as information carrier, it is allowed to which signal of communication is detected using dolphin whistle signal intrinsic in ocean, but in identification process In be taken as marine biological noise to exclude, reach the effect of covert communications, trial solves the problems, such as covert communications under water
A kind of bionical underwater sound based on difference Pattern pattern time-delay codings is disclosed in Chinese patent CN103368660A to lead to Letter method.Dolphin cry signal is applied in difference Pattern pattern time-delay coding communication systems by this method, bionical hidden to reach Cover the purpose of communication.Disclose a kind of bionical based on Virtual time reversal mirror M members in Chinese patent specification CN103401619A The underwater acoustic communication method of Signal coding, this method is equally modulated using dolphin cry signal, realizes covert communications.But, Both communication systems are encoded using time delay, rather than BOK modulation;Meanwhile, both, which communicate, uses the ticktack of dolphin Acoustical signal, not whistle signal.Real dolphin cry signal is communicated, sample quality of the communication efficiency dependent on use, Using inconvenience.Chinese patent specification CN102355440A discloses a kind of underwater sound based on very small linear frequency modulation keying modulation Communication means, this method realizes high-rate information transmission using the ultra-narrow band signal of very small linear frequency modulation keying modulation, but simultaneously Do not emphasize the concealment of communication.
The content of the invention
It is an object of the invention to provide a kind of simple and easy to apply, amount of calculation is less, dolphin whistle is imitated using LFM signals are segmented The underwater acoustic communication method of signal.
It is a kind of to be comprised the following steps using the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals,
Step one:Dolphin whistle signal is chosen, Short Time Fourier Transform is carried out to dolphin whistle signal, its time-frequency spectrum is obtained Contour curve figure;
Step 2:The time shaft of dolphin whistle time-frequency spectral curve equidistantly take a little according to symbol width, Mei Geshi Between put respectively correspond to a dolphin whistle signal instantaneous frequency values fi, obtain one group of frequency values f corresponding with time point;
Step 3:Coding to be sent is obtained by modulating, if coding number is not the class frequency value f intermediate frequencies obtained The integral multiple of the number of rate, then supplement code element behind coding to be sent so that encode number after supplement to obtain one group The integral multiple of the number of frequency in frequency values f;
Step 4:BOK modulation is carried out to sent coding, the frequency modulation rate value corresponding to each code element is determined, obtains Tune group of frequencies k after one group of modulation;
Step 5:Each frequency values in one class frequency value f of acquisition are circulated in order and are used as every section of LFM signal Reference frequency, using symbol width T as the duration of every section of LFM signal, by each the tune frequency values adjusted in group of frequencies k in order As the frequency modulation rate of every section of LFM signal, every section of LFM signal is produced, by obtained LFM signals by frequency values f order in the time On be combined, form one section of segmentation LFM and imitate whistle data-signal;
Step 6:The multistage segmentation LFM of acquisition is imitated into whistle data-signal to be combined in order, every section of imitative whistle data Blank protection interval is added between signal, total imitative whistle data-signal is obtained;
Step 7:Synchronizing signal is the dolphin whistle signal that obtains in step one, by synchronizing signal and total imitative whistle number It is believed that number between one section of blank protection interval of addition be combined, obtain bionical communication transmitting frame signal;
Step 8:Bionical communication transmitting frame signal is launched, after underwater acoustic channel, signal is received in receiving terminal;
Step 9:Signal will be received and carry out copy related operation with dolphin whistle signal, find and corresponded at peak-peak Time point be synchronizing signal starting point, add the duration of synchronizing signal and the duration of blank protection interval in step 7, really Determine the start time point of bionic data signal;
Step 10:To having determined that the bionic data signal of start time point carries out equilibrium, the bionical number after being equalized It is believed that number;
Step 11:To the bionic data signal after equilibrium, intercepted by length of symbol width, and exclude sky White protection interval;
Step 12:Every section of LFM signal after interception is demodulated, modulation intelligence is obtained.
The present invention is a kind of using the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals, can also include
1st, every section of LFM signal is:
Wherein, T is symbol width, fiFor i-th of frequency values, kiCorresponding frequency modulation rate.
2nd, the bionic data signal after equilibrium is:
R (t)=sr(t)*h'(-t)
=[s (t) * h (t)] * h'(-t)+n (t) * h'(-t)
=s (t) * [h (t)] * h'(-t)+n (t) * h'(-t)
Wherein, sr(t) it is the signal received, s (t) is that frame signal is launched in bionical communication, and n (t) is noise signal, h (t) For the impulse response function of channel, h'(t) it is the channel impulse response function estimated, pass through the synchronizing signal to receiving and original The whistle signal of beginning makees the related acquisition of copy.
3rd, the method being demodulated to every section of LFM signal after interception is Fourier Transform of Fractional Order parameter Estimation demodulation method Either match correlation demodulation method.
4th, Fourier Transform of Fractional Order parameter Estimation demodulation method is:
Bionic data signal is segmented according to symbol width, every section of LFM signals y is extracted successivelyi(t) y, is soughti(t) Fourier Transform of Fractional Order Y under different ordersp(u), parameter Estimation is carried out to obtain
It is rightCarry out polarity judgement:IfThen i-th section of LFM signal demodulation resultIfThen i-th section of LFM Signal demodulation result
Finally all demodulation results are carried out and closed, final modulation intelligence is obtained
5th, matching correlation demodulation method is:
Bionic data signal is segmented according to symbol width, every section of LFM signals y is extracted successivelyi(t), produce respectively It is raw to correspond to " 1 " every section of LFM signals s1i(t) and corresponding to every section of LFM signals s of " 0 "0i(t) y, is then obtained respectivelyi(t) and s1i(t) correlation function R1i(τ), and yiAnd s (t)0i(t) correlation function R0i(τ):
Then R is contrasted1i(τ) and R0iThe size of (τ) relevant peaks, to determine the corresponding code element of i-th section of LFM signal:If R1i (τ)>R0i(τ), then i-th section of LFM signal demodulation resultIf R1i(τ)<R0i(τ), then i-th section of LFM signal demodulation result
Finally the demodulation result of every section of LFM signal is combined, final demodulation code is obtained
Beneficial effect:
By emulation experiment and tank experiments, the feasibility of this method is verified.- 10dB is realized in emulation experiment Zero bit error rate communicates under the conditions of signal to noise ratio, and tank experiments realize zero mistake under the conditions of 19 rice grain pattern communication distances, 20 milliseconds of symbol widths Bit rate communication.
Using the communication means of the present invention communicated using the imitative dolphin whistle signals of LFM, it is possible to use segmentation LFM Frequency signal is modulated demodulation as carrier wave, and simple and easy to apply, amount of calculation is smaller, and traffic rate is moderate.Meanwhile, the imitative sea after synthesis Globefish whistle signal is disguised higher, it is possible to achieve camouflage covert communications.
Brief description of the drawings
Fig. 1 is using being segmented the flow chart that the imitative dolphin whistle signals of LFM are communicated;
Fig. 2 is true dolphin whistle signal time-frequency spectrum curve map;
Fig. 3 is the true dolphin whistle signal fundamental wave time-frequency spectrum curve map of one extracted;
Fig. 4 is modulation principle figure;
Fig. 5 is bionical segmentation LFM spectrum diagrams;
Fig. 6 is transmission signal frame structure schematic diagram;
Fig. 7 is transmission signal waveform time-frequency spectrum curve map;
Fig. 8 is subjunctive mood time reversal mirror schematic diagram;
Fig. 9 is Fourier Transform of Fractional Order demodulation principle figure;
Figure 10 is matching relative theory figure.
Embodiment
The present invention is described in further details below in conjunction with accompanying drawing.
The technical problem to be solved in the present invention is to propose that a kind of utilize is segmented the imitative dolphin whistle signal of LFM signals, is realized imitative The method of raw hidden underwater sound communication.
The object of the invention technical scheme is realized, as shown in Figure 1:
It is a kind of to utilize the camouflage underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals, it is characterised in that:
In transmitting terminal:
Step 1:Suitable dolphin whistle signal is chosen, Short Time Fourier Transform is carried out to dolphin whistle signal, it is obtained Time-frequency spectrum contour curve figure;
Step 2:The time shaft of the dolphin whistle time-frequency spectral curve described in step 1 is equidistantly taken according to symbol width Point, each time point corresponds to the instantaneous frequency values f of a dolphin whistle signali, so that obtain one group it is corresponding with time point Frequency values f;
Step 3:Obtain modulation intelligence;
Step 4:The modulation intelligence that step 3 is obtained carries out BOK modulation, determines that the frequency modulation rate corresponding to each code element takes Value k0Or k1, and its putting in order according to symbol order, the tune group of frequencies k after one group of modulation is obtained, for example:0 corresponds to k0 =-40000Hz/s, 1 corresponds to k1=+40000Hz/s, is encoded to 011010, then the tune group of frequencies k obtained is: k0k1k1k0k1k0
Step 5:It regard each frequency values in the class frequency value f that step 2 is obtained as every section of LFM signal in order Reference frequency, each using symbol width T as the duration of every section of LFM signal, in the tune group of frequencies k obtained using step 4 Adjust frequency values in order as the frequency modulation rate of every section of LFM signal, produce every section of LFM signal, these signals are obtained by step 2 Frequency values f order is combined in time, is formed one section of segmentation LFM and is imitated whistle data-signal;
Step 6:Total bionic data signal is obtained, content includes
(1) product that coding total duration is equal to he number and symbol width in the modulation intelligence described in step 3 is defined, it is fixed Adopted code length is total he number of the coding described in step 3;
(2) if coding number is less than the number of frequency in the class frequency value f that step 2 is obtained, that is, total duration is encoded small In the length of whistle signal, then to ensure the integrality of imitative whistle signal, supplement one that should be appropriate behind original code element The code element of fixed number amount so that code length is identical with the number of frequency in the class frequency value f that step 2 is obtained after supplement, then Step 5 is carried out again;
(3) if code length is more than the number of frequency in the class frequency value f that step 2 is obtained, that is, total duration is encoded big In the length of whistle signal, then the frequency values f that can be obtained step 2 is recycled, and repeats step 5, until surplus Remaining code length is less than after f number, (2) in repeat step 6;
(4) when supplementing code element, trend changed over time according to dolphin whistle signal time-frequency spectrum contour curve in step 1 and entered Row complement code, for example:If time-frequency spectrum contour curve is on a declining curve with the time, 0 yard should be supplemented;If time-frequency spectrum contour curve is at any time Between it is in rising trend, then should mend 1 yard;
(5) the multistage segmentation LFM of acquisition is imitated into whistle data-signal to be combined in order, every section of imitative whistle data-signal Between add blank protection interval, obtain total imitative whistle data-signal;
Step 7:Synchronizing signal is the real dolphin whistle signal described in step 1, and synchronizing signal is obtained with step 6 Imitative whistle data-signal between one section of blank protection interval of addition be combined, obtain bionical communication transmitting frame signal;
Step 8:The transmitting frame signal that step 7 is obtained is launched, after underwater acoustic channel, is received and is believed in receiving terminal Number;
In receiving terminal:
Step 9:The docking collection of letters number is synchronized, and method is will to receive signal with the dolphin whistle signal described in step 1 to enter Row copy related operation, finds peak-peak place corresponding time point, the as starting point of synchronizing signal, and then plus step The duration of synchronizing signal described in 7 and the duration of blank protection interval, determine the start time point of bionic data signal;
Step 10:Equilibrium is carried out to the bionic data signal for having determined time starting point described in step 9, is equalized Bionic data signal afterwards;
Step 11:Bionic data signal after the equilibrium obtained to step 10, is being intercepted by length of symbol width, And exclude blank protection interval;
Step 12:Every section of LFM signal after interception described in step 10 is demodulated, modulation intelligence is obtained.
Above-mentioned is imitated in the communication means that dolphin whistle signal is communicated using LFM, the dolphin whistle used in step 1 Signal is the whistle signal for all kinds of dolphins that actual acquisition is arrived, and will be entered on a timeline with symbol width frequency spectrum profiles curve at that time Row equidistantly divide after choose correspondence time point frequency values as every section of LFM signal reference frequency, using symbol width as The duration of every section of LFM signal, is modulated using the slope of every section of LFM signal, can obtain listening domain from time-frequency domain and human ear The upper imitative whistle data-signal similar to dolphin whistle.Original whistle signal is added before bionic data signal again as synchronous letter Number, you can form bionical communication frame.In transmitting procedure, even if the signal of communication is intercepted and captured by enemy, can also identification when quilt As biological cry intrinsic in ocean so as to be excluded, the purpose of bionical covert communications has been reached.It is first when receiving terminal is demodulated Equilibrium is carried out first with subjunctive mood time reversal mirror, is then demodulated using Fourier Transform of Fractional Order or matching correlation, Obtain modulation intelligence.
Some explanation is done to the present invention below:Dolphin whistle signal in step 1 refers to the whistle signal of all kinds of dolphins. Used time-frequency spectrum contour curve is dolphin whistle signal fundamental wave time-frequency spectrum contour curve in step 2.Every section described in step 5 The reference frequency of LFM signals refers to the centre frequency of every section of LFM signal.When the equalization methods used in step 10 is subjunctive moods Between time reversal equalization methods.Blank protection interval length described in step 6 is 50ms.Blank protection interval described in step 7 is long Spend for 1s.Synchronizing signal described in step 7 is the whistle signal described in step 1.The demodulation method used in step 12 can be with For Fourier Transform of Fractional Order parameter Estimation demodulation method.The demodulation method used in step 12 can be matching correlation demodulation side Method.
The camouflage underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals is utilized the present invention relates to a kind of.
Dolphin main signal of communication during the whistle signal of dolphin, its video envelope is similar with chirp signals, therefore can be with Dolphin whistle signal is simulated using LFM signals are segmented, the communication information is carried using LFM signals, is communicated, is realized imitative Dolphin whistle signal covert communications.
In transmitting terminal:
Step 1:Suitable dolphin whistle signal is chosen, Short Time Fourier Transform is carried out to dolphin whistle signal, it is obtained Time-frequency spectrum contour curve figure, as shown in Figures 2 and 3.
Step 2:The time shaft of the dolphin whistle time-frequency spectral curve described in step 1 is equidistantly taken according to symbol width Point, each time point corresponds to the instantaneous frequency values f of a dolphin whistle signali, so that obtain one group it is corresponding with time point Frequency values f;
Step 3:Modulation intelligence is obtained, as shown in Figure 4;
Step 4:The modulation intelligence that step 3 is obtained carries out BOK modulation, determines that the frequency modulation rate corresponding to each code element takes Value k0Or k1, and its putting in order according to symbol order, the tune group of frequencies k after one group of modulation is obtained, for example:0 corresponds to k0 =-40000Hz/s, 1 corresponds to k1=+40000Hz/s, is encoded to 011010, then the tune group of frequencies k obtained is: k0k1k1k0k1k0
If unipolar binary to be modulated is encoded to a={ 1,0,1,1,0,0,1 ... }, BOK modulation is carried out to it, correspondence Tune group of frequencies k={ k1,k0,k1,k1,k0,k0,k1…}。
Step 5:It regard each frequency values in the class frequency value f that step 2 is obtained as every section of LFM signal in order Reference frequency, each using symbol width T as the duration of every section of LFM signal, in the tune group of frequencies k obtained using step 4 Adjust frequency values in order as the frequency modulation rate of every section of LFM signal, produce every section of LFM signal, these signals are obtained by step 2 Frequency values f order is combined in time, is formed one section of segmentation LFM and is imitated whistle data-signal, as shown in Figure 5;
If i-th of frequency values is fi, corresponding frequency modulation rate is ki, then i-th section of linear FM signal s can be obtainedi(t) it is:
Wherein, T is symbol width.The instantaneous frequency of the signal is:
Step 6:Total bionic data signal is obtained, content includes
(1) product that coding total duration is equal to he number and symbol width in the modulation intelligence described in step 3 is defined, it is fixed Adopted code length is total he number of the coding described in step 3;
(2) if coding number is less than the number of frequency in the class frequency value f that step 2 is obtained, that is, total duration is encoded small In the length of whistle signal, then to ensure the integrality of imitative whistle signal, supplement one that should be appropriate behind original code element The code element of fixed number amount so that code length is identical with the number of frequency in the class frequency value f that step 2 is obtained after supplement, then Step 5 is carried out again;
(3) if code length is more than the number of frequency in the class frequency value f that step 2 is obtained, that is, total duration is encoded big In the length of whistle signal, then the frequency values f that can be obtained step 2 is recycled, and repeats step 5, until surplus Remaining code length is less than after f number, (2) in repeat step 6;
(4) when supplementing code element, trend changed over time according to dolphin whistle signal time-frequency spectrum contour curve in step 1 and entered Row complement code, for example:If time-frequency spectrum contour curve is on a declining curve with the time, 0 yard should be supplemented;If time-frequency spectrum contour curve is at any time Between it is in rising trend, then should mend 1 yard;
(5) the multistage segmentation LFM of acquisition is imitated into whistle data-signal to be combined in order, every section of imitative whistle data-signal Between add blank protection interval, obtain total imitative whistle data-signal;
Step 7:Synchronizing signal is the real dolphin whistle signal described in step 1, and synchronizing signal is obtained with step 6 Imitative whistle data-signal between one section of blank protection interval of addition be combined, obtain bionical communication transmitting frame signal, such as Fig. 6 With shown in Fig. 7;
Step 8:The transmitting frame signal that step 7 is obtained is launched, after underwater acoustic channel, is received and is believed in receiving terminal Number;
In receiving terminal:
Step 9:The docking collection of letters number is synchronized, and method is will to receive signal with the dolphin whistle signal described in step 1 to enter Row copy related operation, finds peak-peak place corresponding time point, the as starting point of synchronizing signal, and then plus step The duration of synchronizing signal described in 7 and the duration of blank protection interval, determine the start time point of bionic data signal;
Step 10:Equilibrium is carried out to the bionic data signal for having determined time starting point described in step 9, is equalized Bionic data signal afterwards;
Equalization methods use subjunctive mood time reversal mirror balancing technique.As shown in figure 8, in this patent, detectable signal and Synchronizing signal is same signal.In receiving terminal, to the detectable signal p receivedr(t) copied with original detectable signal p (t) Correlation, so as to estimate the impulse response function h'(t of channel), then being carried out time reversal obtains h'(-t).It will receive The data letter s arrivedrNumber (t) and h'(-t) make convolution, obtain the signal r (t) finally virtually received, its waveform is similar to former letter Waveform is ceased, portfolio effect is reached.Virtual time reversal mirror processing procedure expression formula is as follows:
sr(t)=s (t) * h (t)+n (t) (3)
R (t)=sr(t)*h'(-t)
=[s (t) * h (t)] * h'(-t)+n (t) * h'(-t) (4)
=s (t) * [h (t)] * h'(-t)+n (t) * h'(-t)
Note:Referred to as virtual time reversal channel, can be considered as that signal eventually passes through is " total Channel ".After signal is by virtual time reversal channel, direct signal is more much higher than multi-path signalses energy, it is suppressed that underwater acoustic channel The intersymbol interference that multi-path effect is brought, simultaneously as noise is non-coherent addition, the technology can also improve signal to noise ratio.
Step 11:Bionic data signal after the equilibrium obtained to step 10, is being intercepted by length of symbol width, And exclude blank protection interval;
Step 12:Every section of LFM signal after interception described in step 10 is demodulated, modulation intelligence is obtained.
Demodulation method can either match correlation demodulation method for Fourier Transform of Fractional Order parameter Estimation demodulation method.
(1) Fourier Transform of Fractional Order parameter Estimation demodulation method:
In demodulation, it is first according to symbol width and bionic data signal is segmented, every segment signal is extracted successively.Such as Shown in Fig. 9.Assuming that extracting i-th section of LFM signals y in bionic data signali(t) y, is soughti(t) fraction under different orders Rank Fourier transformation Yp(u), then utilize
Formula (5) carries out parameter Estimation to it, obtainsIt is rightCarry out polarity judgement:IfThen i-th section of LFM signal solution Adjust resultIfThen i-th section of LFM signal demodulation resultFinally all demodulation results are carried out and closed, are obtained To final modulation intelligence
(2) correlation demodulation method is matched
Because LFM signals have good correlated performance, the correlation peak of its auto-correlation function is very precipitous, similar to δ Function, therefore can be demodulated using correlation technique is matched, as shown in Figure 10.
In demodulation, it is first according to symbol width and bionic data signal is segmented, every segment signal is extracted successively. Receiving terminal, the frequency waypoint obtained according to step 2 produces correspond to " 1 " and every section of LFM signal corresponding to " 0 " respectively.Assuming that carrying Take out i-th section of LFM signals y in bionic data signali(t), i-th section of LFM signal for corresponding to " 1 " is s1i(t), i-th section it is right Ying Yu " 0 " LFM signals are s0i(t) y, is then obtained respectivelyiAnd s (t)1i(t) correlation function R1i(τ), and yi(t) and s0i(t) correlation function R0i(τ), its calculation formula is as follows:
Then R is contrasted1i(τ) and R0iThe size of (τ) relevant peaks, to determine the corresponding code element of i-th section of LFM signal:If R1i (τ)>R0i(τ), then i-th section of LFM signal demodulation resultIf R1i(τ)<R0i(τ), then i-th section of LFM signal demodulation resultFinally the demodulation result of every section of LFM signal is combined, final demodulation code is obtained
The camouflage underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals, its feature are utilized the present invention relates to a kind of It is:The characteristics of being changed over time according to dolphin whistle signal frequency, utilizes multiple LFM signal subsections combine analog dolphin whistles Signal, is realized to the modulates informations of LFM signals using BOK methods, and then is realized imitative based on dolphin whistle signal and LFM signals Unboiled water sound communication is modulated.In receiving terminal by Fourier Transform of Fractional Order to every section of LFM signal in the reception data after equilibrium Centre frequency and frequency modulation rate are estimated, and the dolphin whistle signal sample used during based on modulation, are demodulated, are modulated Information.The underwater acoustic communication method of dolphin whistle signal is imitated using segmentation LFM, it is possible to achieve utilize simple linear FM signal The bionic camouflage underwater sound communication of high fidelity is carried out, strengthens the disguise and practicality of communication.

Claims (6)

1. a kind of utilize the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals, it is characterised in that:Including following step Suddenly,
Step one:Dolphin whistle signal is chosen, Short Time Fourier Transform is carried out to dolphin whistle signal, frequency spectrum profiles at that time are obtained Curve map;
Step 2:The time shaft of dolphin whistle time-frequency spectral curve equidistantly take a little according to symbol width, each time point The instantaneous frequency values f of a dolphin whistle signal is corresponded to respectivelyi, obtain one group of frequency values f corresponding with time point;
Step 3:Coding to be sent is obtained by modulating, if coding number is not frequency in the class frequency value f obtained The integral multiple of number, then supplement code element behind coding to be sent so that the class frequency for encoding number after supplement to obtain The integral multiple of the number of frequency in value f;
Step 4:BOK modulation is carried out to sent coding, the frequency modulation rate value corresponding to each code element is determined, obtains one group Tune group of frequencies k after modulation;
Step 5:Each frequency values in one class frequency value f of acquisition are circulated to the reference as every section of LFM signal in order Frequency, using symbol width T as the duration of every section of LFM signal, using adjust group of frequencies k in each adjust frequency values in order as The frequency modulation rate of every section of LFM signal, produces every section of LFM signal, obtained LFM signals is entered in time by frequency values f order Row combination, forms one section of segmentation LFM and imitates whistle data-signal;
Step 6:The multistage segmentation LFM of acquisition is imitated into whistle data-signal to be combined in order, every section of imitative whistle data-signal Between add blank protection interval, obtain total imitative whistle data-signal;
Step 7:Synchronizing signal is the dolphin whistle signal obtained in step one, and synchronizing signal is believed with total imitative whistle data One section of blank protection interval is added between number to be combined, and obtains bionical communication transmitting frame signal;
Step 8:Bionical communication transmitting frame signal is launched, after underwater acoustic channel, signal is received in receiving terminal;
Step 9:Signal and dolphin whistle signal will be received and carry out copy related operation, find peak-peak place it is corresponding when Between point be the starting point of synchronizing signal, the duration of synchronizing signal and the duration of blank protection interval in step 7 are added, it is determined that imitative The start time point of raw data-signal;
Step 10:To having determined that it is balanced that the bionic data signal of start time point is carried out, the bionic data letter after being equalized Number;
Step 11:To the bionic data signal after equilibrium, intercepted by length of symbol width, and excluding blank guarantor Shield interval;
Step 12:Every section of LFM signal after interception is demodulated, modulation intelligence is obtained.
2. according to claim 1 a kind of using the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals, it is special Levy and be:Every section of described LFM signal is:
<mrow> <msub> <mi>s</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mn>2</mn> <mi>&amp;pi;</mi> <mi>t</mi> <mo>(</mo> <mrow> <msub> <mi>f</mi> <mi>i</mi> </msub> <mo>-</mo> <mfrac> <mrow> <msub> <mi>k</mi> <mi>i</mi> </msub> <mi>T</mi> </mrow> <mn>2</mn> </mfrac> </mrow> <mo>)</mo> <mo>+</mo> <msub> <mi>&amp;pi;k</mi> <mi>i</mi> </msub> <msup> <mi>t</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow>
Wherein, T is symbol width, fiFor i-th of frequency values, kiCorresponding frequency modulation rate.
3. according to claim 1 a kind of using the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals, it is special Levy and be:It is described it is balanced after bionic data signal be:
R (t)=sr(t)*h'(-t)
=[s (t) * h (t)] * h'(-t)+n (t) * h'(-t)
=s (t) * [h (t)] * h'(-t)+n (t) * h'(-t)
Wherein, sr(t) it is the signal received, s (t) is that frame signal is launched in bionical communication, and n (t) is noise signal, and h (t) is letter The impulse response function in road, h'(t) for estimation channel impulse response function, by the synchronizing signal to receiving with it is original Whistle signal makees the related acquisition of copy.
4. according to claim 1 a kind of using the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals, it is special Levy and be:The described method being demodulated to every section of LFM signal after interception is Fourier Transform of Fractional Order parameter Estimation solution Tune method either matches correlation demodulation method.
5. according to claim 4 a kind of using the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals, it is special Levy and be:Described Fourier Transform of Fractional Order parameter Estimation demodulation method is:
Bionic data signal is segmented according to symbol width, every section of LFM signals y is extracted successivelyi(t) y, is soughti(t) not With the Fourier Transform of Fractional Order Y under orderp(u), parameter Estimation is carried out to obtain
<mrow> <mo>{</mo> <mi>p</mi> <mo>,</mo> <msub> <mi>u</mi> <mn>0</mn> </msub> <mo>}</mo> <mo>=</mo> <mi>arg</mi> <munder> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mrow> <mi>p</mi> <mo>,</mo> <mi>u</mi> </mrow> </munder> <mo>|</mo> <msub> <mi>Y</mi> <mi>p</mi> </msub> <mrow> <mo>(</mo> <mi>u</mi> <mo>)</mo> </mrow> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow>
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mover> <mi>k</mi> <mo>^</mo> </mover> <mo>=</mo> <mo>-</mo> <mi>cot</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>p</mi> <mi>&amp;pi;</mi> </mrow> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mi>f</mi> <mo>^</mo> </mover> <mn>0</mn> </msub> <mo>=</mo> <msub> <mi>u</mi> <mn>0</mn> </msub> <mi>csc</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>p</mi> <mi>&amp;pi;</mi> </mrow> <mn>2</mn> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
It is rightCarry out polarity judgement:IfThen i-th section of LFM signal demodulation resultIfThen i-th section of LFM signal Demodulation result
Finally all demodulation results are combined, final modulation intelligence is obtained
6. according to claim 5 a kind of using the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals, it is special Levy and be:Described matching correlation demodulation method is:
Bionic data signal is segmented according to symbol width, every section of LFM signals y is extracted successivelyi(t) correspondence, is produced respectively In " 1 " every section of LFM signals s1i(t) and corresponding to every section of LFM signals s of " 0 "0i(t) y, is then obtained respectivelyiAnd s (t)1i(t) Correlation function R1i(τ), and yiAnd s (t)0i(t) correlation function R0i(τ):
<mrow> <msub> <mi>R</mi> <mrow> <mn>1</mn> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>|</mo> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <mi>&amp;infin;</mi> </mrow> <mrow> <mo>+</mo> <mi>&amp;infin;</mi> </mrow> </msubsup> <msub> <mi>y</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msub> <mi>s</mi> <mrow> <mn>1</mn> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> <mo>|</mo> </mrow>
<mrow> <msub> <mi>R</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>|</mo> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <mi>&amp;infin;</mi> </mrow> <mrow> <mo>+</mo> <mi>&amp;infin;</mi> </mrow> </msubsup> <msub> <mi>y</mi> <mi>i</mi> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msub> <mi>s</mi> <mrow> <mn>0</mn> <mi>i</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mo>-</mo> <mi>&amp;tau;</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> <mo>|</mo> </mrow>
Then R is contrasted1i(τ) and R0iThe size of (τ) relevant peaks, to determine the corresponding code element of i-th section of LFM signal:If R1i(τ) > R0i(τ), then i-th section of LFM signal demodulation resultIf R1i(τ) < R0i(τ), then i-th section of LFM signal demodulation result
Finally the demodulation result of every section of LFM signal is combined, final demodulation code is obtained
CN201510676931.8A 2015-10-13 2015-10-13 It is a kind of to utilize the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals Active CN105227246B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510676931.8A CN105227246B (en) 2015-10-13 2015-10-13 It is a kind of to utilize the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510676931.8A CN105227246B (en) 2015-10-13 2015-10-13 It is a kind of to utilize the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals

Publications (2)

Publication Number Publication Date
CN105227246A CN105227246A (en) 2016-01-06
CN105227246B true CN105227246B (en) 2017-10-31

Family

ID=54995982

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510676931.8A Active CN105227246B (en) 2015-10-13 2015-10-13 It is a kind of to utilize the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals

Country Status (1)

Country Link
CN (1) CN105227246B (en)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105833410B (en) * 2016-05-27 2022-03-29 无锡丹辰智能科技有限公司 Children autism auxiliary treatment equipment for realizing dolphin-like sound ultrasonic device
CN106503336B (en) * 2016-10-21 2019-09-27 哈尔滨工程大学 A kind of modeling of dolphin ticktack acoustical signal and synthetic method
CN106789794B (en) * 2016-11-30 2019-05-28 厦门大学 A kind of bionical communication means modulated using dolphin signal of communication frequency modulation
CN106899357B (en) * 2017-03-13 2020-04-03 哈尔滨工程大学 Disguised and concealed underwater communication device simulating dolphin whistle
CN108055087B (en) * 2017-12-30 2024-04-02 天津大学 Communication method and device for coding by using number of long-limb piloting whale sound harmonics
CN108667532B (en) * 2018-04-11 2021-09-28 哈尔滨工程大学 Burst type underwater acoustic communication method
CN108957412A (en) * 2018-07-12 2018-12-07 西安电子科技大学 Radar detection and method for communication transmission based on piece-wise linear FM signal
CN109085543B (en) * 2018-08-24 2020-08-18 斯凯瑞利(北京)科技有限公司 Linear frequency modulation signal generator
CN109347568B (en) * 2018-09-05 2021-04-20 哈尔滨工程大学 Dolphin whistle-imitating continuous phase multi-element frequency modulation underwater acoustic communication method
CN109347569B (en) * 2018-09-05 2021-05-14 哈尔滨工程大学 Disguised and concealed underwater acoustic communication method based on discrete cosine transform
CN111431625B (en) * 2020-02-25 2021-07-27 天津大学 Synthetic and modification method for vocals of cetacea animals
CN112953652B (en) * 2021-01-27 2022-04-15 浙江大学 Whale whistle-like communication method based on segmented time-frequency contour time delay modulation
CN113259016B (en) * 2021-04-23 2022-09-13 西北工业大学 Underwater sound concealed communication method based on dolphin whistle parameter model
CN115037386A (en) * 2022-05-16 2022-09-09 广西民族大学 Bionic communication signal simulation test method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3987455B2 (en) * 2003-05-06 2007-10-10 日本電信電話株式会社 Digital modulation transmitter
CN101813771A (en) * 2009-12-08 2010-08-25 中国科学院声学研究所 Dolphin biomimetic sonar signal processing method
CN104967489A (en) * 2015-06-08 2015-10-07 哈尔滨工程大学 Biomimetic underwater acoustic communication method of MSK signal embedded dolphin whistle signal

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030072456A1 (en) * 2001-10-17 2003-04-17 David Graumann Acoustic source localization by phase signature

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3987455B2 (en) * 2003-05-06 2007-10-10 日本電信電話株式会社 Digital modulation transmitter
CN101813771A (en) * 2009-12-08 2010-08-25 中国科学院声学研究所 Dolphin biomimetic sonar signal processing method
CN104967489A (en) * 2015-06-08 2015-10-07 哈尔滨工程大学 Biomimetic underwater acoustic communication method of MSK signal embedded dolphin whistle signal

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
M元仿海豚叫声隐蔽水声通信;刘凇佐等;《哈尔滨工程大学学报》;20131115;第35卷(第1期);全文 *
一种利用海豚叫声的仿生水声通信方法;刘凇佐等;《物理学报》;20130528;第62卷(第14期);全文 *

Also Published As

Publication number Publication date
CN105227246A (en) 2016-01-06

Similar Documents

Publication Publication Date Title
CN105227246B (en) It is a kind of to utilize the underwater acoustic communication method for being segmented the imitative dolphin whistle signal of LFM signals
Leus et al. Multiband OFDM for covert acoustic communications
Liu et al. Biologically inspired covert underwater acoustic communication by mimicking dolphin whistles
CN109347568B (en) Dolphin whistle-imitating continuous phase multi-element frequency modulation underwater acoustic communication method
Jia-jia et al. Bio-inspired steganography for secure underwater acoustic communications
CN103501201B (en) A kind of frequency hopping arteries and veins position based on linear FM signal coding underwater acoustic communication method
CN105391500B (en) A kind of imitative dolphin ticktack sound underwater acoustic communication method based on ultra-broadband signal
CN109412993A (en) A kind of having based on deep learning inhibits the full duplex cognition water sound communication signal detection method of self-interference
Liu et al. Covert underwater acoustic communication using whale noise masking on DSSS signal
CN112953652B (en) Whale whistle-like communication method based on segmented time-frequency contour time delay modulation
CN107359899B (en) Orthogonal frequency division multiplexing spread-spectrum underwater sound communication is without pilot tone judgment feedback channel estimation method under the conditions of condition of sparse channel
CN105515683A (en) Differential chaos shift keying communication method based on hybrid system
CN106899357B (en) Disguised and concealed underwater communication device simulating dolphin whistle
Lou et al. Underwater Communications and Networks
Jiang et al. Covert underwater communication based on combined encoding of diverse time-frequency characteristics of sperm whale clicks
Liu et al. Biologically inspired covert underwater acoustic communication using high frequency dolphin clicks
He et al. Passive time reversal communication with cyclic shift keying over underwater acoustic channels
CN114362835A (en) Dolphin whistle-imitated underwater acoustic communication method
CN111355677A (en) Multi-carrier underwater high-speed communication system based on filter bank
CN103441980A (en) Sideband-information-free shallow sea underwater acoustic communication pattern selection peak-to-average ratio restraining algorithm based on frequency reversal mirror technology
CN106910508B (en) Hidden underwater acoustic communication method for imitating marine pile driving sound source
Liu et al. Jamming recognition based on feature fusion and convolutional neural network
Thottappilly OFDM for underwater acoustic communication
CN107426197A (en) A kind of secret communication method for realizing hidden voice call
CN101610117B (en) Surface mother ship and underwater user two-way acoustic communication method

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant