Disclosure of Invention
Aiming at the defects or improvement requirements of the prior art, the invention provides a frequency domain suppression method and a frequency domain suppression system for inter-carrier interference of underwater acoustic OFDM (orthogonal frequency division multiplexing), which are used for solving the technical problem that the data detection performance of the underwater acoustic OFDM system with a medium-high Doppler factor is limited due to the fact that ICI (inter-carrier interference) cannot be suppressed to the greatest extent by taking a subcarrier interval as a reference frequency offset in the prior art.
In order to achieve the above object, in a first aspect, the present invention provides a frequency domain suppression method for inter-carrier interference in an underwater acoustic OFDM system, including the following steps:
s1, constructing an ICI frequency domain suppression target optimization model based on the minimum mean square error principle, and estimating the optimal reference frequency offset of the current data frame by minimizing the objective function of the ICI frequency domain suppression target optimization model according to the pilot signals in the lead codes;
s2, after carrying out frequency spectrum shifting on the received signal for a plurality of times at preset frequency intervals, carrying out FFT operation on the signals subjected to frequency spectrum shifting by the same point number respectively to obtain a plurality of paths of FFT signals at each subcarrier; the preset frequency interval is a fraction times of the optimal reference frequency offset of the current data frame;
and S3, carrying out weighted combination on the generated multipath FFT signals subcarrier by subcarrier to obtain signals after the frequency domain suppression of the inter-carrier interference.
Further preferably, the step S1 includes:
s11, constructing an ICI frequency domain suppression target optimization model based on a minimum mean square error principle;
s12, dividing the pilot signals in the lead code into two parallel signals A, B; the channel A signal is a pilot signal in the lead code, and the channel B signal is a product of the pilot signal and the sampling moment;
s13, initializing the current reference frequency offset and the weight value of the sub-carrier corresponding to the pilot signal;
s14, carrying out frequency spectrum shifting on the A, B two paths of signals for a plurality of times at a frequency interval preset by the fractional multiple of the current reference frequency offset, and carrying out FFT operation on the signals subjected to frequency spectrum shifting with the same number of points to obtain a plurality of paths of FFT signals;
s15, weighting and combining the multiple paths of FFT signals corresponding to the path A signals to obtain demodulated signals on the path A pilot frequency; forming vectors by using multiple paths of FFT signals corresponding to the B path of signals, and performing weighted combination after a group of Hadamard product operation to obtain demodulated signals of the B path of pilot frequency;
s16, calculating a target value of the ICI frequency domain suppression target optimization model based on the demodulated signal on the pilot frequency of the path A, judging whether the ratio of the target values of the two adjacent times is smaller than a preset threshold value, and if so, turning to the step S17; otherwise, the operation is finished, and the current reference frequency offset is the optimal reference frequency offset of the current data frame;
s17, updating the current reference frequency offset and weight value based on the demodulated signal of the pilot frequency of the path A and the demodulated signal of the pilot frequency of the path B, and turning to the step S14;
the target function of the ICI frequency domain suppression target optimization model is:
wherein f is
eIs the current reference frequency offset;
the weight value at the p pilot signal subcarrier; b
kIs the original data symbol; k
P={k
1,...,k
P-1,k
P};k
pA subcarrier subscript corresponding to the p-th pilot signal in the preamble; p1, 2,. said, P; p is the number of pilot signals;
the demodulated signal on the p pilot frequency of the A path;
and a column vector formed by FFT signals corresponding to the p-th pilot signal of the A path.
Further preferably, in step S13, the A, B two signals are respectively spaced at intervals of frequency in frequency interval
Performing 2I-time frequency spectrum shifting to obtain (2I +1) paths of signals respectively;
wherein, the path A signal is v (t), and the path B signal is tv (t); the frequency spectrum shifted signal corresponding to the A path signal is
The frequency spectrum shifted signal corresponding to the B path signal is
I is a positive integer.
Further preferably, the ith FFT signal corresponding to the pth pilot signal of the a channel is:
vi(t) is the signal after the frequency spectrum shift corresponding to the A path signal, and Δ f is the subcarrier interval.
Further preferably, the demodulated signal on the p-th pilot in the B-path is:
wherein,
column vectors formed by FFT signals corresponding to the p-th pilot signal of the B path; β [ -I, · 1, 0, 1,. 1, I]
T;
Representing a Hadamard product operation;
the ith FFT signal corresponding to the p-th pilot signal of the B path is:
tvi(t) is a signal after the frequency spectrum corresponding to the B-path signal is shifted, and Δ f is a subcarrier interval.
Further preferably, in step S2, the multipath FFT signals at the k-th subcarrier are:
zk=[zk,-I,...,zk,-1,zk,0,zk,1,...,zk,I]T
wherein,
i is a positive integer, r (t) is a received signal, f
baseThe best reference frequency offset for the current data frame.
Further preferably, the multiple FFT signals z at the k-th sub-carrier are combined
kWeighting and combining to obtain a signal x after the inter-carrier interference frequency domain suppression on the kth path of sub-carrier
kWherein
w
kis the weight value at the k-th subcarrier.
Further preferably, the frequency domain suppression method for inter-carrier interference in the underwater acoustic OFDM system is used for symbol detection; the method specifically comprises the following steps: after steps S1-S2 are performed, the generated multiple FFT signals are weighted and combined subcarrier by subcarrier, and symbol detection is performed; and when the symbol detection is finished once, updating a weight value corresponding to the next subcarrier based on the current symbol detection error, and carrying out symbol detection after carrying out weighted combination on the multipath FFT signals of the next subcarrier by adopting the updated weight value.
Further preferably, the update formula for updating the weight value corresponding to the k +1 th subcarrier based on the detection error detected by the symbol at the k +1 th subcarrier is as follows:
wk+1=wk+μgk
wherein, w
k+1The weight value of the k +1 th path of subcarrier; mu is an iteration step length; g
kThe iterative gradient of the weight specifically includes:
"+" indicates a conjugate operation.
In a second aspect, the present invention provides a frequency domain suppression system for inter-carrier interference in an underwater acoustic OFDM system, including:
the optimal reference frequency offset estimation module is used for constructing an ICI frequency domain suppression target optimization model based on the minimum mean square error principle, estimating the optimal reference frequency offset of the current data frame by minimizing the objective function of the ICI frequency domain suppression target optimization model according to the pilot signals in the lead codes, and outputting the optimal reference frequency offset to the multi-path FFT signal acquisition module;
the multi-path FFT signal acquisition module is used for carrying out frequency spectrum shifting on the received signal for a plurality of times at preset frequency intervals, respectively carrying out FFT operation on the signals subjected to frequency spectrum shifting with the same number of points to obtain multi-path FFT signals at each subcarrier, and outputting the multi-path FFT signals to the symbol detection module; the preset frequency interval is a fraction times of the optimal reference frequency offset of the current data frame;
and the weighting combination module is used for carrying out weighting combination on the generated multipath FFT signals subcarrier by subcarrier to obtain signals after the inter-carrier interference frequency domain suppression.
In a third aspect, the present invention also provides a machine-readable storage medium storing machine-executable instructions which, when invoked and executed by a processor, cause the processor to implement any one of the above-mentioned methods for frequency-domain suppression of inter-carrier interference in an underwater acoustic OFDM system.
Generally, compared with the prior art, the above technical solution conceived by the present invention has the following beneficial effects:
1. the invention provides a frequency domain suppression method for inter-carrier interference in an underwater acoustic OFDM system, which is characterized in that a target optimization model of an ICI frequency domain suppression system is constructed based on a minimum mean square error criterion, the selection of the optimal reference frequency offset is converted into a multivariable joint optimization problem through a two-way parallel processing structure, and the multivariable joint optimization is realized under the condition of ensuring more ideal performance, namely, the optimal reference frequency offset of the current data frame is estimated by minimizing a target function of the ICI frequency domain suppression target optimization model, so that the ICI frequency domain suppression precision is improved to the maximum extent, and the system detection performance can be obviously improved under medium-high Doppler factors.
2. Compared with the existing F-FFT, the frequency domain suppression method for the inter-carrier interference in the underwater sound OFDM system carries out frequency spectrum shifting on the received signal for a plurality of times by taking the fractional multiple of the optimal reference frequency offset as the frequency interval without considering only [ F [ -F ] isk-Δf,fk+Δf]Sub-carrier interference in a frequency range, where fkAnd f is the center frequency of the kth carrier, and is the subcarrier interval. Under the condition of medium and high Doppler factors, the detection performance of the frequency ICI suppression method can be remarkably improved, and the application range of the frequency ICI suppression in different Doppler expansion scenes is widened.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention. In addition, the technical features involved in the embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
To achieve the above object, as shown in fig. 1, the present invention provides a frequency domain suppression method for inter-carrier interference in an underwater acoustic OFDM system, including the following steps:
s1, constructing an ICI frequency domain suppression target optimization model based on the minimum mean square error principle, and estimating the optimal reference frequency offset of the current data frame by minimizing the objective function of the ICI frequency domain suppression target optimization model according to the pilot signals in the lead codes;
s2, after carrying out frequency spectrum shifting on the received signal for a plurality of times at preset frequency intervals, carrying out FFT operation on the signals subjected to frequency spectrum shifting by the same point number respectively to obtain a plurality of paths of FFT signals at each subcarrier; the preset frequency interval is a fraction times of the optimal reference frequency offset of the current data frame;
and S3, carrying out weighted combination on the generated multipath FFT signals subcarrier by subcarrier to obtain signals after the frequency domain suppression of the inter-carrier interference.
To further illustrate the methods provided by the present invention, the following examples are given:
examples 1,
In this embodiment, differential coherent detection is taken as an example to describe the method flow of the present invention in detail, and the frequency domain suppression method for inter-carrier interference in an underwater acoustic OFDM system provided by the present invention is mainly applied to a receiving end.
Specifically, after applying the frequency domain suppression method for inter-carrier interference provided by the present invention to a receiving end, the complete communication process of the underwater acoustic OFDM system is as follows:
a sending end:
PSK symbol b to be transmittedkAfter serial-to-parallel conversion, differential coding in frequency domain is carried out in sequence to obtain symbol d to be modulatedkI.e. by
Wherein, bk∈{a0,a1,...,aQ-1},aq=ej2πq/QQ is 0, 1,., Q-1, Q is the modulation order, and K is the total number of subcarriers.
Will encode the symbol dkThe signal is transmitted to a channel through operations of IFFT modulation, cyclic prefix addition, parallel-to-serial conversion, and frequency up-conversion, and the transmission signal can be represented as
Wherein f isk=f0+ k Δ f is the kth subcarrier frequency, f0For the lowest frequency of the system, Δ f — B/K is the subcarrier spacing, T is the OFDM block duration, and B is the signal bandwidth.
Receiving end:
the received signal is processed by down-conversion, rough sampling, cyclic prefix removal and other operations to obtain a signal to be demodulated;
specifically, as shown in fig. 2, the frequency domain suppression method for inter-carrier interference in an underwater acoustic OFDM system provided by the present invention includes the following steps:
s1, constructing an ICI frequency domain suppression target optimization model based on the minimum mean square error principle, and estimating the optimal reference frequency offset of the current data frame by minimizing the objective function of the ICI frequency domain suppression target optimization model according to the pilot signals in the lead codes;
specifically, a two-way parallel processing structure is adopted, and the optimal reference frequency offset in the data frame is estimated according to the pilot frequency in the preamble signal v (t). The double-path parallel processing structure comprises a signal preprocessing module, a multi-path FFT operation module, a Hadamard product operation module, a combiner and a gradient updating module; the signal preprocessing module is used for generating A, B two-path signals which are processed in parallel; the multi-path FFT operation module is used for carrying out frequency spectrum shifting on the preprocessed signals at preset frequency intervals and then respectively carrying out FFT operation on the same points; the Hadamard product operation module consists of a group of Hadamard product devices and is used for Hadamard product operation of one path of signals; and the combiner is used for carrying out weighted combination on the multipath FFT-processed signals to obtain modulated signals on each subcarrier. And the gradient updating module is used for updating the reference frequency offset value according to the signal processing information at the pilot frequency until a certain cutoff condition is met.
The specific process is as follows:
s11, constructing an ICI frequency domain suppression target optimization model based on a minimum mean square error principle;
s12, dividing the pilot signals in the lead code into two parallel signals A, B; the path A signal is a pilot signal v (t) in the lead code, and the path B signal is a product tv (t) of the pilot signal and the sampling time;
s13, initializing the current reference frequency offset and the weight value of the sub-carrier corresponding to the pilot signal; specifically, the initial value of the current reference frequency offset may be set to Δ f or an integer multiple of the doppler frequency offset; where Δ f is the inter-subcarrier spacing.
S14, carrying out frequency spectrum shifting on the A, B two paths of signals for a plurality of times at a frequency interval preset by the fractional multiple of the current reference frequency offset, and carrying out FFT operation on the signals subjected to frequency spectrum shifting with the same number of points to obtain a plurality of paths of FFT signals;
specifically, the A path signal v (t) is separated by a preset frequency interval
In 2I frequency spectrum shifts, a (2I +1) channel signal (including the original preamble signal v (t)) can be generated, and the mathematical expression is as follows:
wherein f is
e(j) The reference frequency offset is the reference frequency offset under the jth iteration number;
for generated multipath signal vi(t) respectively performing the same Fourier operation to obtain the ith path of FFT signal corresponding to the pth pilot signal:
wherein k ispA subcarrier subscript corresponding to the p-th pilot signal in the preamble; k is a radical ofp∈{k1,...,kP-1,kP}; p is the number of pilot signals.
The path B signal is the product of the pilot signal and the sampling time, i.e. the signal tv (t), and the same processing is performed as above, so that the ith path FFT signal corresponding to the p-th pilot signal of the path B is obtained as follows:
s15, weighting and combining the multiple paths of FFT signals corresponding to the path A signals to obtain demodulated signals on the path A pilot frequency; forming vectors by using multiple paths of FFT signals corresponding to the B path of signals, and performing weighted combination after a group of Hadamard product operation to obtain demodulated signals of the B path of pilot frequency;
specifically, for the p-th pilot, the (2I +1) -way output on the A-way
Forming a column vector
The demodulated signal on the p pilot frequency can be obtained after weighted combination by the combiner
Wherein,
the weight value of the p pilot signal subcarrier under the jth iteration number is obtained; similarly, the p-th pilot signal subcarrier is constituted by (2I +1) channels, and the weight value here is a vector constituted by weight values of the p-th pilot signal subcarrier.
(2I +1) path output on B path
Forming a column vector
Then obtaining the demodulation signal on the p pilot frequency after Hadamard product operation and weighted combination:
wherein β [ -I, · 1, 0, 1,. 1, I]
T,
Representing a hadamard product operation.
S16, calculating a target value of the ICI frequency domain suppression target optimization model based on the demodulated signal on the pilot frequency of the path A, judging whether the ratio of the target values of the two adjacent times is smaller than a preset threshold value, and if so, turning to the step S17; otherwise, the operation is finished, and the current reference frequency offset is the optimal reference frequency offset of the current data frame;
specifically, the objective function of the ICI frequency domain suppression objective optimization model is:
wherein f is
eIs the current reference frequency offset;
the weight value at the p pilot signal subcarrier; b
kThe method comprises the steps of obtaining original data symbols, namely PSK symbols to be transmitted; k
P={k
1,...,k
P-1,k
P};k
pA subcarrier subscript corresponding to the p-th pilot signal in the preamble; p1, 2,. said, P; p is the number of pilot signals;
the demodulated signal on the p pilot frequency of the A path;
and a column vector formed by FFT signals corresponding to the p-th pilot signal of the A path.
The invention judges whether the cut-off condition is met by judging whether the ratio of the target values of two adjacent times is smaller than a preset threshold value, wherein the cut-off condition in the embodiment is as follows: judging whether the requirements are met
Wherein E (j) is a target value under the j iteration number; e (j-1) is the target value for the j-1 th iteration number.
S17, based on the demodulated signal of the pilot frequency in the A path and the demodulated signal of the pilot frequency in the B path, the current reference frequency offset and the weight value are updated, and the step goes to S14.
Specifically, the embodiment is based on the minimum mean square error principle, and adopts a coordinate descent method to perform iterative update to solve the optimal reference frequency offset of the current data frame; the specific updating process is as follows:
……
for each sub-variable, the reference frequency offset iterative formula is:
fe(j+1)=fe(j)+μγ(j)
wherein,
is the k-th
pSymbol error at each pilot, "+" indicates conjugate operation.
The weight iterative formula is
S2, after carrying out frequency spectrum shifting on the received signal for a plurality of times at preset frequency intervals, carrying out FFT operation on the signals subjected to frequency spectrum shifting by the same point number respectively to obtain a plurality of paths of FFT signals at each subcarrier; the preset frequency interval is a fraction times of the optimal reference frequency offset of the current data frame;
specifically, the multipath FFT signals at the kth subcarrier in the received signal are:
zk=[zk,-I,...,zk,-1,zk,0,zk,1,...,zk,I]T
wherein,
i is a positive integer, r (t) is a received signal, f
baseThe best reference frequency offset for the current data frame.
And S3, carrying out weighted combination on the generated multipath FFT signals subcarrier by subcarrier to obtain signals after the frequency domain suppression of the inter-carrier interference.
Preferably, the frequency domain suppression method for inter-carrier interference in the underwater acoustic OFDM system can be used for symbol detection; at this time, the above-mentioned two-way parallel processing structure further includes a data detection module for recovering the original data symbols from the demodulated signal, as shown in fig. 1.
The method specifically comprises the following steps: after steps S1-S2 are performed, the generated multiple FFT signals are weighted and combined subcarrier by subcarrier, and symbol detection is performed (i.e., symbol detection is performed after the multiple FFT signals at the subcarriers are weighted and combined in sequence for each subcarrier in the received signal); and when the symbol detection is finished once, updating a weight value corresponding to the next subcarrier based on the current symbol detection error, and carrying out symbol detection after carrying out weighted combination on the multipath FFT signals of the next subcarrier by adopting the updated weight value. The specific process is as follows:
s31, making k equal to 0;
s32, receiving the multipath FFT signal z at the k path subcarrier in the signal
kWeighted combination to obtain the demodulated signal x on the k path sub-carrier
kWherein
w
kthe weight value at the k-th path of subcarrier, i.e. each path of weight value of the combiner.
S33, demodulating signal x on k path sub carrier
kCarrying out data detection to obtain detection symbols
And for the detected symbol
Performing symbol decision to obtain decision symbol
Specifically, the data detection in the present embodiment is differential coherent detection; data symbols on the kth subcarrier
The expression of (a) is as follows:
after symbol decision, the decision symbol is obtained as:
wherein,
indicating symbol decision, i.e. judgement
The Hamming distance from each point in the Q-order PSK symbol set is the final one when the Hamming distance is the shortest
And (4) taking values.
S34, calculating decision symbol
And detecting the symbol
The difference is obtained to obtain the detection error e of symbol detection at the k path subcarrier
kAnd calculating an iterative gradient g of the weights
k(ii) a Judging whether e is satisfied
kLess than a predetermined error threshold e
thAnd is
Less than a predetermined gradient threshold g
th(ii) a If yes, go to step S35; otherwise, updating the weight value corresponding to the k +1 th path of subcarriers based on the detection error of the symbol detection at the k path of subcarriers according to the minimum mean square error principle;
specifically, the detection error e of symbol detection at the k-th subcarrierkComprises the following steps:
according to the minimum mean square error principle, the update formula of the iterative gradient of the weight is as follows:
if | e
k|<e
thAnd is
The weights at the next carrier are updated
wk+1=wk+μgk
Wherein, w
k+1The weight value of the k +1 th path of subcarrier; mu is an iteration step length and is a constant close to 0; g
kThe iterative gradient of the weight specifically includes:
"+" indicates a conjugate operation.
Otherwise, no update is performed.
In this embodiment, the iterative gradient g of the weightthA value of 2, a predetermined error threshold ethThe value is 1.
S35, let k equal to k + 1;
and S36, repeating the steps S32-S35 to iterate until the symbol detection at each subcarrier is completed.
In order to further explain the intercarrier interference suppression effect of the method provided by the invention, the intercarrier interference suppression precision of the method (marked as A-FFT) provided by the invention and the classical P-FFT and F-FFT in the underwater sound OFDM is respectively tested. Fig. 3 is a time delay-path gain diagram of an underwater acoustic simulation channel for testing, and based on the simulation environment setting, the inter-carrier interference suppression precision of the classic P-FFT, F-FFT and the method a-FFT of the present invention for underwater acoustic OFDM is analyzed and compared. The specific experimental results are as follows:
(1) detection results under different Doppler factors
As shown in fig. 4, which is a graph of the MSE performance of the method provided by the present invention and the existing method according to the test result of Doppler factor (Doppler a) variation, it can be seen from fig. 3 that, under the condition of SNR 30dB, the a-FFT provided by the present invention has the performance equivalent to the classical P-FFT and F-FFT at low Doppler factor, while in the differential underwater acoustic OFDM system with medium and high Doppler factor, the a-FFT provided by the present invention exhibits the performance advantage of Mean Square Error (MSE) which is more significant.
(2) Detection results at different SNR of medium Doppler factor
FIG. 5 is a graph showing the MSE performance of the method of the present invention and the prior art method as a function of the signal-to-noise ratio, and FIG. 4 shows that the MSE performance is 2.5 × 10-4Compared with P-FFT and F-FFT, the MSE performance of the A-FFT method provided by the invention is reduced by 16.61-52.91%, and the performance is better.
Examples 2,
A frequency domain suppression system for inter-carrier interference in an underwater acoustic OFDM system, comprising:
the optimal reference frequency offset estimation module is used for constructing an ICI frequency domain suppression target optimization model based on the minimum mean square error principle, estimating the optimal reference frequency offset of the current data frame by minimizing the objective function of the ICI frequency domain suppression target optimization model according to the pilot signals in the lead codes, and outputting the optimal reference frequency offset to the multi-path FFT signal acquisition module;
the multi-path FFT signal acquisition module is used for carrying out frequency spectrum shifting on the received signal for a plurality of times at preset frequency intervals, respectively carrying out FFT operation on the signals subjected to frequency spectrum shifting with the same number of points to obtain multi-path FFT signals at each subcarrier, and outputting the multi-path FFT signals to the symbol detection module; the preset frequency interval is a fraction times of the optimal reference frequency offset of the current data frame;
and the weighting combination module is used for carrying out weighting combination on the generated multipath FFT signals subcarrier by subcarrier to obtain signals after the inter-carrier interference frequency domain suppression.
The related technical scheme is the same as embodiment 1, and is not described herein.
In a third aspect, the present invention also provides a machine-readable storage medium storing machine-executable instructions, which when invoked and executed by a processor, cause the processor to implement the method for frequency-domain suppression of inter-carrier interference in an underwater acoustic OFDM system according to embodiment 1.
The related technical scheme is the same as embodiment 1, and is not described herein.
It will be appreciated by those skilled in the art that the foregoing is only a preferred embodiment of the invention, and is not intended to limit the invention, such that various modifications, equivalents and improvements may be made without departing from the spirit and scope of the invention.