Sub-band selection activation-based multi-band hyperbolic frequency modulation spread spectrum underwater acoustic communication method

Abstract

Disclosed by the present invention is a sub-band selection activation-based multi-band hyperbolic frequency modulation spread spectrum underwater acoustic communication method. The present invention discloses: dividing the available bandwidth of an underwater acoustic system into a plurality of sub-bands, performing hyperbolic frequency modulation on each of the sub-bands respectively, and performing spread spectrum modulation on the plurality of sub-bands within the same frequency modulation period, thus implementing multi-band parallel transmission. Hence, within each frequency modulation period, the divided plurality of sub-bands is grouped, and each sub-band group activates different sub-bands for transmission according to different options for transmitting data. Compared to other underwater acoustic hyperbolic frequency modulation communication solutions, the present invention further improves the frequency band utilization of the system, and the energy efficiency is also improved.

Claims

1. A sub-band selection activation-based multi-band hyperbolic frequency modulation spread spectrum underwater acoustic communication method, comprising the following steps: S1, dividing a bandwidth of a communication system into a plurality of sub-bands, and performing hyperbolic frequency modulation respectively; S2, grouping data and sub-bands, determining a sub-band activation series of steps, and performing modulation on a signal; S3, adding a frame header to the signal after modulation; S4, performing, at a receiving end, preprocessing, synchronization, and channel estimation on a received signal; S5, calculating correlation values for sub-band detection and demodulation; S6, detecting an activated sub-band; and S7, demodulating the activated sub-band.

2. The sub-band selection activation-based multi-band hyperbolic frequency modulation spread spectrum underwater acoustic communication method of claim 1, wherein in step S1, an available bandwidth B=f.sub.1−f.sub.0 of the system, which is a difference between a highest frequency f.sub.1 and a lowest frequency f.sub.0, is divided into K sub-bands, and a value of K is an integer multiple of 2; the bandwidth of each of the sub-bands is $\frac{B}{K},$ hyperbolic frequency modulation is performed on each of the sub-bands respectively with a frequency modulation period being T, for subsequent use in spread spectrum modulation; and a hyperbolic frequency modulation signal corresponding to a k(1≤k≤K)th sub-band is: $x_k\left[n\right]=\exp \left(j\frac{2\pi }{{\alpha }_k}\ln \left(1+{\alpha }_k{f}_{k0}\frac{n}{f_s}\right)\right),n=0,\mathrm{.Math.},N-1$ where f.sub.s is a sampling rate, $N=f_{s}T,{\alpha }_k=-\frac{B}{f_{k0}\left({\mathrm{Kf}}_{k0}+B\right)T}\mathrm{and}{f}_{k0}=f_0+\left(k-1\right)\frac{B}{K}.$

3. The sub-band selection activation-based multi-band hyperbolic frequency modulation spread spectrum underwater acoustic communication method of claim 1, wherein in step S2, the data and sub-bands are grouped, and the grouping is specifically: grouping every K bits into Q data groups, a value of Q being a number of the sub-bands divided by 2, i.e. $\frac{K}{2},$ and meanwhile, also grouping the K sub-bands into Q sub-band groups, each of the data groups and each of the sub-band groups being corresponding to each other and forming as one group, wherein a result of the grouping is that two sub-bands in each group transmit two bits.

4. The sub-band selection activation-based multi-band hyperbolic frequency modulation spread spectrum underwater acoustic communication method of claim 1, wherein in step S2, determining the sub-band activation series of steps, and performing modulation on the signal are specifically: activating, when a first bit of the two bits in each group is 0, a first sub-band in each group for transmission, and activating, when that bit is 1, a second sub-band in the group for transmission; using BPSK symbol mapping for a second bit of the two bits in each group to modulate a carrier signal of the activated sub-band; and denoting two bits in a j(j=1, . . . , Q)th group as b.sub.j1b.sub.j2, then a sending signal of a current multi-band hyperbolic frequency modulation spread spectrum symbol frame being expressed as $s\left[n\right]=\underset{j=1}{\overset{Q}{.Math.}}\left(2b_{j2}-1\right)x_{k_j}\left[n\right]n=0,\mathrm{.Math.},N-1$ where k.sub.j=2j+b.sub.j1−1 is a number of an activated sub-band in the jth group, and 2b.sub.j2−1 represents that BPSK symbol mapping is used for the second bit in the jth group.

5. The sub-band selection activation-based multi-band hyperbolic frequency modulation spread spectrum underwater acoustic communication method of claim 1, wherein in step S5, implementing an operation of correlation values for sub-band detection and demodulation includes: S51, calculating a correlation peak offset, and assuming that a result of the channel estimation in step S4 is that: a number of multipaths is L, and parameters corresponding to each path are a delay {circumflex over (τ)}.sub.1 and a Doppler factor {circumflex over (β)}.sub.1; and calculating a peak offset for a k th sub-band, an m th symbol, and an i th path $\Delta n_{k,m,l}=.Math.{\hat{\tau }}_l*f_s+\left(\frac{1-{\hat{\beta }}_l}{{\hat{\beta }}_l{\alpha }_k{f}_{k0}}\right)*f_s+\left(m-1\right)*\frac{N}{{\hat{\beta }}_l}.Math.$ where N is a number of points of a hyperbolic frequency modulation symbol; and f.sub.s is a signal sampling rate, $\frac{1-{\hat{\beta }}_l}{{\hat{\beta }}_l{\alpha }_k{f}_{k0}}$ is a shift caused by Doppler spread, f.sub.k0 represents a starting frequency of the kth sub-band, and $\frac{N}{{\hat{\beta }}_l}$ is a number of points of a hyperbolic frequency modulation symbol after Doppler spread compensation; and S52, performing a correlation operation to obtain a peak value; based on Δn.sub.k, m, l, respectively using each sub-band and a corresponding position of a received signal to perform the correlation operation so as to obtain the peak value; and the peak value I.sub.k, m, l for the kth sub-band, the mth symbol, and the ith path being expressed as $I_{k,m,l}=\underset{n=0}{\overset{N-1}{.Math.}}{x}_k\left[n\right]r\left[n+\Delta n_{k,m,l}\right]$ where x.sub.k[n] is a spread spectrum symbol of the kth sub-band, and r[n] represents a signal at the receiving end.

6. The sub-band selection activation-based multi-band hyperbolic frequency modulation spread spectrum underwater acoustic communication method of claim 5, wherein in S6, letting j=1, . . . , Q, a detection of a jth group of activated sub-bands includes: if $\underset{l=1,\mathrm{.Math.}L}{.Math.}.Math.I_{2j-1,m,l}.Math.>\underset{l=1,\mathrm{.Math.}L}{.Math.}.Math.I_{2j,m,l}.Math.,$ a first sub-band in the sub-band group is considered as the activated sub-band, and meanwhile, according to the modulation of step S2, it may be determined that {circumflex over (β)}.sub.jl=0 in symbols correspondingly transmitted by the sub-band group, and the number of the activated sub-band is k.sub.j=2j−1; and if $\underset{l=1,\mathrm{.Math.}L}{.Math.}.Math.I_{2j-1,m,l}.Math.\le \underset{l=1,\mathrm{.Math.}L}{.Math.}.Math.I_{2j,m,l}.Math.,$ a second sub-band in the sub-band group is considered as the activated sub-band, and meanwhile, according to a modulation regulation in S2, it can be determined that {circumflex over (β)}.sub.jl=1 in symbols correspondingly transmitted by the sub-band group, and the number of the activated sub-band is k.sub.j=2j.

7. The sub-band selection activation-based multi-band hyperbolic frequency modulation spread spectrum underwater acoustic communication method of claim 1, wherein in S7, letting j=1, . . . , Q, and the demodulation of symbols transmitted on a jth group of activated sub-bands is specifically implemented as follows: a maximum merge ratio criterion is used to perform merge processing on correlation peak values of L paths and make a sign decision, $S_{j2}=\mathrm{sign}\left(\mathrm{Re}\left\{\underset{l=1}{\overset{L}{.Math.}}{\hat{A}}_l^{*}{I}_{k_j,m,l}\right\}\right)$ where Â.sub.l is an estimated amplitude of an lth path, and then BPSK symbol demapping is performed, that is, ${\hat{b}}_{J^2}=\frac{S_{j2}+1}{2}.$

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a schematic structural diagram of an underwater acoustic hyperbolic frequency modulation spread spectrum communication solution combining multi-band transmission and carrier combination according to the present invention.

(2) FIG. 2 is a flowchart of sub-band group spread spectrum modulation with a first sub-band group as an example.

(3) FIG. 3 is a flowchart of sub-band group de-spreading with a first sub-band group as an example.

(4) FIG. 4 is a structural diagram of a signal frame.

(5) FIG. 5 is a simulation diagram of a bit error rate of a specific embodiment.

DETAILED DESCRIPTION OF EMBODIMENTS

(6) In order to make the objectives, technical solutions, and advantages of the present invention clearer and more specific, the present invention is further described in detail below with reference to the accompanying drawings and embodiments. It is to be understood that the specific embodiments described herein are merely illustrative of the present invention and are not intended to limit the present invention.

(7) With reference to FIG. 1, FIG. 1 is a schematic diagram of a system structure of the present invention.

(8) Meanings of various reference signs are as follows:

(9) B: bandwidth, where B=6 KHz in this example.

(10) f.sub.0: the starting frequency of the bandwidth, where f.sub.0=9 KHz in this example.

(11) f.sub.1: the ending frequency of the bandwidth, where f.sub.1=15 KHz in this example.

(12) T.sub.H: the symbol period of a hyperbolic frequency modulation signal, where T.sub.H=10 ms in this example.

(13) T.sub.G: nil symbol guard interval cycle, where T.sub.G=20 ms in this example.

(14) T.sub.L: LFM frame header duration, where T.sub.L=20 ms in this example.

(15) f.sub.s: the sampling rate of the hyperbolic frequency modulation signal, where f.sub.s=100 KHz in this example.

(16) N: the number of sampling points of the hyperbolic frequency modulation signal within a period, N=T.sub.H*f.sub.s, where N=1000 in this example.

(17) α.sub.k: the frequency modulation rate of the hyperbolic frequency modulation signal.

(18) f.sub.k0: representing the starting frequency of the k th sub-band.

(19) K: the number of sub-bands, where K=8 in this example.

(20) Q: the sub-bands being grouped into Q groups, where Q=4 in this example.

(21) M: the number of multi-band hyperbolic frequency modulation symbols contained in each frame of signal, where M=10 in this example.

(22) In this embodiment, an underwater acoustic hyperbolic frequency modulation spread spectrum communication solution combining multi-band transmission and carrier combination comprises the following steps:

(23) S1, generating and grouping hyperbolic frequency modulation spread spectrum sub-band carrier signals. B=f.sub.1−f.sub.0 is divided into K=8 sub-bands, the bandwidth for each of the sub-bands is

(24) $\frac{B}{K},$
and hyperbolic frequency modulation is performed on each of the sub-bands respectively for spread spectrum modulation. A hyperbolic frequency modulation carrier signal corresponding to the k(k=1, . . . , K) th sub-band is

(25) 0$\begin{array}{cc}{x}_k\left[n\right]=\exp \left(j\frac{2\pi }{{\alpha }_k}\ln \left(1+{\alpha }_k{f}_{k0}\frac{n}{f_s}\right)\right),n=0,\mathrm{.Math.},N-1,k=1,\mathrm{.Math.},K,& \left(1\right)\end{array}$
where

(26) ${\alpha }_k=-\frac{B}{f_{k0}\left(K{f}_{k0}+B\right)T_H}$
is the hyperbolic frequency modulation rate of the kth sub-band, and

(27) $f_{k0}=f_0+\left(k-1\right)\frac{B}{K}$
is the starting frequency of the k th sub-band.

(28) Meanwhile, K sub-bands are grouped into Q=4 groups, each group having two sub-bands.

(29) S2, signal modulation. Every K bits in sent data are used to modulate one multi-band hyperbolic frequency modulation spread spectrum symbol. Every K bits are grouped into Q groups to correspond to Q sub-band groups, each group having two bits. The two bits in each group are used to control the activation of the group of sub-band carriers and the modulation of each of the activated sub-band carriers, and a specific implementation is as follows:

(30) activating, when the first bit of the two bits in each group is 0, the first sub-band in each group for transmission, and activating, when that bit is 1, the second sub-band in the group for transmission;

(31) using BPSK symbol mapping for the second bit of the two bits in each group to modulate a carrier signal of the activated sub-band; and

(32) denoting two bits in the j (j=1, . . . , Q) th group as b.sub.j1b.sub.j2, then a sending signal of a current multi-band hyperbolic frequency modulation spread spectrum symbol frame being expressed as

(33) $\begin{array}{cc}s\left[n\right]=\underset{j=1}{\overset{Q}{.Math.}}\left(2b_{j2}-1\right)x_{k_j}\left[n\right]n=0,\mathrm{.Math.},N-1& \left(2\right)\end{array}$
where k.sub.j=2j+b.sub.j1−1 is the number of an activated sub-band in the jth group, and in the above formula, 2b.sub.j2−1 represents that BPSK symbol mapping is used for the second bit in the jth group.

(34) The above-mentioned signal modulation manner is used to modulate M multi-band hyperbolic frequency modulation symbols to form one frame of signal.

(35) S3, adding a frame header. A frame header is added to the signal after spread spectrum modulation. The frame header uses an LFM signal, and the LFM signal can be expressed as

(36) $\begin{array}{cc}c\left[n\right]=\cos \left(2\pi f_{0}n+\pi \frac{f_1-f_0}{T}{n}^2\right)n=0,\mathrm{.Math.},N_L-1& \left(3\right)\end{array}$
where N.sub.L=T.sub.L*f.sub.s.

(37) The structure of the signal frame is shown in FIG. 4, and T.sub.L, T.sub.G, and T.sub.H all can be adjusted according to actual requirements. After DA conversion and analog transmission front-end power amplification, the above-mentioned data frame is transmitted in an underwater acoustic channel.

(38) S4, receiving preprocessing. After the signal undergoes multipath and Doppler effects in the underwater acoustic channel, at the receiving end, band pass filtering is performed first, and then the frame header is used for synchronization and channel estimation; the present invention does not specifically discuss the synchronization and channel estimation, and general synchronization and channel estimation algorithms may be used. The received signal that has been synchronized is denoted as r[n], and it is assumed that a result after the channel estimation is that the number of multipaths is L, and parameters corresponding to each path are an amplitude Â.sub.l, a delay {circumflex over (τ)}.sub.l, and a Doppler factor {circumflex over (β)}.sub.l.

(39) S5, calculating correlation values for sub-band detection and demodulation. Here, the first sub-band group is taken as an example.

(40) S51, calculating a correlation peak offset. The peak offset for each path is calculated respectively based on the channel estimation parameters in S4. For the m(1, . . . , M) th symbol in a data frame, the peak offset for the l(1, . . . , L) th path thereof should be

(41) $\begin{array}{cc}\Delta n_{k,m,l}=.Math.{\hat{\tau }}_l*f_s+\left(\frac{1-{\hat{\beta }}_l}{{\hat{\beta }}_l{\alpha }_k{f}_{k0}}\right)*f_s+\left(m-1\right)*\frac{N}{{\hat{\beta }}_l}.Math.k=1,\mathrm{.Math.},Q& \left(4\right)\end{array}$
where

(42) $\frac{1-{\hat{\beta }}_l}{{\hat{\beta }}_l{\alpha }_k{f}_{k0}}$
is a shift caused by Doppler spread, f.sub.k0 represents the starting frequency of the k th sub-band, and

(43) $\frac{N}{{\hat{\beta }}_l}$
is the number of points of a hyperbolic frequency modulation symbol after Doppler spread compensation.

(44) S52, performing a correlation operation to obtain a peak value. Based on Δn.sub.k,m,l, the hyperbolic frequency modulation carrier signal x.sub.k[n] of each sub-band and a corresponding position of a received signal r[n] are respectively used to perform a correlation operation so as to obtain a peak value. For each data frame, the peak value I.sub.k,m,l for the k th sub-band, the m th symbol, and the l th path can be expressed as

(45) $\begin{array}{cc}{I}_{k,m,l}=\underset{n=0}{\overset{N-1}{.Math.}}{x}_k\left[n\right]r\left[n+\Delta n_{k,m,l}\right]& \left(5\right)\end{array}$

(46) S6, detecting an activated sub-band. A decision is made on the activated sub-band according to the correlation peak value I.sub.k,m,l obtained in S5, and for any m th symbol in the data frame, the decision criterion for the activated sub-band is as follows (here, the first sub-band group is taken as an example).

(47) If

(48) $\underset{l=1,\mathrm{.Math.},L}{\overset{}{.Math.}}.Math.I_{1,m,l}.Math.>\underset{l=1,\mathrm{.Math.},L}{\overset{}{.Math.}}.Math.I_{2,m,l}.Math.,$
the activated sub-band in the sub-band group is considered as the first sub-band, and meanwhile, according to a modulation regulation in S2, it may be determined that {circumflex over (b)}.sub.11=0 in symbols correspondingly transmitted by the sub-band group.

(49) If

(50) 0$\underset{l=1,\mathrm{.Math.}L}{\overset{}{.Math.}}.Math.I_{1,m,l}.Math.\le \underset{l=1,\mathrm{.Math.}L}{\overset{}{.Math.}}.Math.I_{2,m,l}.Math.,$
the activated sub-band in the sub-band group is considered as the second sub-band, and meanwhile, according to a modulation regulation in S2, it may be determined that {circumflex over (b)}.sub.11=1 in symbols correspondingly transmitted by the sub-band group.

(51) S7, demodulating the activated sub-band. The sequence number k of the activated sub-band is determined according to the result of the detection in S6, and symbol demodulation is performed on the sub-band. Here, a maximum merge ratio criterion is used to perform merge processing on the peak values of L paths and make a sign decision

(52) $\begin{array}{cc}{S}_{12}=\mathrm{sign}\left(\mathrm{Re}\left\{\underset{l=1}{\overset{L}{.Math.}}{\hat{A}}_l^{*}{I}_{k,m,l}\right\}\right)& \left(6\right)\end{array}$
then performing BPSK symbol demapping

(53) $\begin{array}{cc}{\hat{b}}_{12}=\frac{S_{12}+1}{2}& \left(7\right)\end{array}$

(54) Since then, the de-spreading of a sub-band group is completed, and the demodulation of the two symbols b.sub.11b.sub.12 transmitted by the sub-band group is achieved. The remaining sub-band groups are processed accordingly.

(55) In the present invention, the available bandwidth of the system is divided into a plurality of sub-bands for data transmission, which improves the frequency band utilization of the entire system. With the manner of selecting a sub-band to be activated, the energy consumption of the system is reduced on the premise of transmitting the same amount of data. It can also be seen from a simulation diagram of a bit error rate in FIG. 5 that, in this embodiment, after undergoing Doppler spread, the signal can still maintain a relatively good bit error rate (BER) performance at a relatively low signal-to-noise ratio (SNR), and with the manner of selecting a sub-carrier to be activated, an energy loss and interference between adjacent carriers are reduced.