1. Introduction
Underwater acoustic (UWA) communications are dominant technologies for data exchanges in broad sea areas, which have been successfully used to accomplish underwater tasks such as marine environmental monitoring, marine security surveillance and resource exploration. As the requirements for underwater network deployments and large-volume marine data acquisition increase, multicarrier UWA communications have attracted the attention of researchers due to high spectral efficiency (SE), where Orthogonal Frequency Division Multiplexing (OFDM) and index modulation with OFDM (OFDM-IM) are the outstanding and representative schemes. OFDM-IM has the robustness to multipath fading of UWA channels and exploits the idle space resources as compensation for the limited frequency bandwidth. OFDM-IM communications are highly desired to be further studied for emerging underwater applications.
The concept of index modulation (IM) originates from spatial modulation (SM) in multiple-input multiple-output (MIMO) systems. The difference between OFDM-IM from SM-MIMO is that OFDM-IM employs the indices of subcarriers in the frequency domain instead of the indices of antennas [
1]. IM separates the data into two parts, i.e., the constellation data and the index data, hence OFDM-IM conveys useful information using both constellation symbols and subcarrier indices. OFDM-IM is flexible to achieve the tradeoff between bit error rate (BER) performance and SE by designing different subcarrier activation patterns (SAPs), which determine the activated or inactivated subcarriers. In OFDM-IM, additional information requires special handling because the conventional detection for constellation symbols is only part of data detection modules. The detection of SAPs is the prerequisite for constellation data detection, which is the key to avoiding decoding error propagation. For OFDM-IM detection, maximum likelihood (ML) detection is the best, which searches the optimal indices and constellation points among all combinations. The complexity of ML detection increases exponentially with the number of index combinations as well as the order of constellation modulations. Even the searches have been employed for each subblock, the complexity is heavy, especially for the multicarrier systems with large-number subcarriers and high-order constellation modulations. A log-likelihood (LLR) detector can achieve near-optimal ML performance with reduced complexity [
1]. If the ratio of activated subcarriers is low, the LLR detector tends to find nonexistent index combinations with high probability. Subsequently, many sub-optimal detectors are developed [
2,
3,
4]. The work in [
3] proposes a low-complexity detector that conveys data using all possible SAPs to avoid errors from invalid SAPs. The mapping of non-fixed length bits makes the communication system complicated. The work in [
4] proposes a modified k-largest-value (klv) detector, which chooses
k active indices with the largest values according to the active likelihood metrics. This method provides a method for dealing with illegal SAPs and operates on each subcarrier.
As research deepens, the inherent characteristics of OFDM-IM symbols from IM concept and SAPs are further exploited. SAPs are constructed by zero subcarriers and activated subcarriers, which makes OFDM-IM symbols present inherent sparsity characteristics. The sparsity leverages powerful tools to solve the detection problem. The paper [
5] proposes to interpret data detection as a problem of convex optimization, on which the AP constraints are imposed, and then a semi-definite relaxation method is utilized. The convergence speed of this method is severely limited and the error floor exists in high SNR ranges. Exploiting the sparsity of symbols, compressed sensing (CS) based methods have recently emerged. The work [
6] proposes a CS-assisted signaling strategy, based on which an iterative residual check (IRC)-based detector is formulated. An AMP-based detector is proposed for OFDM-IM by exploiting the statistic of OFDM-IM symbols in the frequency domain [
7]. The work [
8] proposes to use the AMP framework to realize iterative channel estimation and data detection for OFDM-IM. Approximate Message Passing (AMP) algorithm belongs to the Bayesian estimation framework [
9,
10]. The simplifications of the Gaussian message according to the central limit theorem and Taylor expansions make the IM detector efficient depending on several posterior parameters. One constraint of the AMP method is the sensitivity to the non-Gaussianity of the dictionary matrix. The generalized AMP (GAMP) is incorporated as a step for joint phase-noise Estimation and decoding in OFDM-IM [
11]. The performance of these OFDM-IM detectors for UWA communications cannot be concluded and is to be verified and analyzed.
OFDM-IM was first introduced to UWA communications in [
12]. The work verifies the reduced PAPR effects due to power averaging by the inactivated subcarriers and better BER performance compared with OFDM. Since then, OFDM-IM-based UWA communications have been developing. Combining the IM with advanced OFDM schemes, many UWA OFDM-IM transceivers are proposed such as the fully quadrature subcarrier-index shift keying OFDM (FQSISK-OFDM) modulation scheme [
13], which put their emphasis on the design of activation strategy. The work [
14] combines the IM concept with Orthogonal Time Frequency Space (OTFS) and proposes a Hamming distance optimized model to modify the index combinations. Currently, there is a lack of a design scheme for a low-complexity OFDM-IM detector for UWA communications and a detailed analysis of the effects from the UWA channel for UWA OFDM-IM systems. To fill this gap, this paper considers the OFDM-IM receiver design and especially focuses on data detection, which is a major challenge for real-time and high SE data transmission schemes in the UWA physical layer.
In this paper, we involve the data detection for OFDM-IM in the Bayesian estimation framework. Instead of a loopy factor graph (FG), a vector approximate message passing (VAMP) detector based on a scalar FG is proposed. An appropriate statistical prior model is beneficial for the achievement of the optimal Bayesian solution. By exploiting the sparsity of symbols in the frequency domain, the SAP constraint is considered using the statistical prior and a novel minimum mean square error (MMSE)-optimal shrinkage function is derived. The data reconstruction performance is improved through the forward-back message passing scheduling, where the involved statistical parameters are learned automatically. Inherently from the robustness of VAMP for the deviations of Gaussian matrices [
15], the proposed detector is less sensitive to the non-Gaussianity of the measurement matrix which is composed of unknown channel components and presents good generality and convergence. Aiming to mitigate the possible invalid SAPs, a modified trick using the criterion of minimal Euclidean distance with the space of a look-up table is used to replace the possible non-existent results. The VAMP-based detection has relatively low complexity, which is advantageous for real-time data transmission. Simulation results verify the proposed receiver outperforms the benchmarks in terms of complexity, BER as well as robustness to the time-varying UWA channel.
The rest of this paper is organized as follows.
Section 2 introduces the UWA OFDM-IM communication systems. The proposed VAMP-based detector is presented in
Section 3. The computation complexity and numerical experiments are analyzed in
Section 4.
Section 5 makes the conclusive remarks.
Notations: Lower case boldface letter and upper case boldface letter denote vector and matrix, respectively. and stand for the transpose and conjugate transpose operation. is the diagonalized operation with as the diagonal element. and denote the expectation and variance operation. denotes the floor operation.
2. System Model
In this paper, the OFDM-IM-based UWA communication system is considered. The system framework is shown in
Figure 1. There are a total
B bits to be transmitted and input to the bit splitter. The system has
N subcarriers, and these subcarriers are split into
G groups. Each group includes
subcarriers. Correspondingly, the bits are also split as
G groups, and there are
p bits in each group. The
p bits include two parts, i.e.,
bits for index selection and
bits for amplitude and phase modulation (APM).
bits are mapped as the indices of
k activated subcarriers, which are extracted from the index combinations in a look-up table.
is determined by
which is the logarithm of the index combination number. The total number of data subcarriers is
. Considering the
g-th group of the OFDM-IM subblock, the indices of activate subcarriers in the look-up table
are given by
where
for
and
. The elements of
are arranged in an ascending order. The size of Table
is
.
bits are mapped as
M-ary constellation points from the constellation alphabet
.
where
is the modulation order. Mapping symbol
are appended on the
k activated subcarriers with normalized power, i.e.,
.
A coherent UWA communication system is required to track the channel effects. Besides data subcarriers, there are
subcarriers allocated as comb pilots for the tracking of time-varying UWA channels. The information of pilot subcarriers is perfectly known. The remaining
subcarriers are idle. The
gth OFDM subblock can be expressed as
. Concatenate
G OFDM-IM subblocks, and OFDM-IM data frame is given by
where
denotes the symbols to be modulated on
N subcarriers. The
mth subcarrier is with the frequency
, where
denotes the subcarrier spacing. The symbol duration is
. Before transmission, the OFDM-IM is transformed as the time-domain signal through the inverse fast Fourier transform (FFT) operation, then the
-length cyclic prefix (CP) is appended to mitigate the inter-symbol interference (ISI), and the time duration of CP
is larger than the maximum path delay. The SE is given by
The baseband signal is upshifted to the passband given by
where
is the pulse shaping filter, and
when
, otherwise
.
As shown in
Figure 1, the signal passes through the UWA channel. The UWA channel is generally expressed as a time-varying model given by
where
L is the total number of paths,
and
denote the path amplitude and path delay for the
lth path. Assuming all paths have the same Doppler scaling factor
a, then
[
16]. After passing through the channel, the received signal is given by
where
is the passband additive white Gaussian noise (AWGN). Due to the Doppler effects, the received signal suffers from compressing and broadening effects. The scalar coefficient
a is estimated by the resampling method as
, where
is the length of the transmitted signal and the length of the received signal
is calculated by cross-correlating the linear frequency modulated (LFM) preambles of neighboring frames. After the resampling and downshifting operation, the baseband received signal is given by
where
denotes the baseband AWGN and
is the residual carrier frequency offset (CFO). To discretize the baseband signal
at the baseband rate
, then the
pth received sample is given by
. Through expansion and union simplification operation, the received samples are denoted as
where
and
is the channel frequency response.
Based on Formula (
8), CFO matrix is defined as
. Part of inactive subcarriers is used for CFO mitigation. The process is operated in the frequency domain based on the goal of minimizing the leaked energy of the zero subcarriers after the signal passes through the channel. After CP removal and FFT transform, the received signal in the frequency domain is denoted as
, where
is the Discrete Fourier Transform (DFT) matrix and
is the received signal in the time domain. The selection matrix
extracts zero subcarriers from
, and the energy is expressed as
. Assuming CFO is perfectly known and compensated, the optimization problem is formulated as
One dimensional (1-D) search method is used to solve the problem (
9) [
17]. The detailed effects of CFOs are analyzed in the following experiments. With
, the received signal components are corrected through phase reversal. Then the inter-carrier interference (ICI) is mitigated and the ICI-free signal is obtained.
Channel state information (CSI) is obtained through the channel estimation module. pilot subcarriers are utilized as input, and the received pilot subcarriers are the output, then the channel estimation problem can be formed as a sparse signal recovery (SSR) problem. To design efficient channel estimation, representative CS-based methods can be involved.
Then frequency-domain channel
is fed to the data detection module. For the detection of OFDM-IM, symbol detection includes index detection and constellation symbol detection. To realize reliable data detection, the signal for the
gth group
in
in the frequency domain is expressed as
where
is a diagonal matrix with the components from the
gth group vector in
and
is the
gth group Gaussian noise vector. ML detection is the optimal method and searches all the combinations of indices and constellation points, which is generally defined as
where
and
,
, are the corresponding signal element and channel coefficient of the
gth OFDM-IM subblock. ML detector achieves the optimal error performance. However, it has high complexity which exponentially increases with the size of subblocks and the modulation order. It is impractical for cases with a large number of subcarriers or high-order modulation in UWA communications. Therefore, it is desired to design a low-complexity and effective detector for UWA communications. In the following contents, we propose a novel detector based on VAMP theory by exploiting the inherent sparsity of OFDM-IM symbols in the frequency domain.
3. Proposed Method
OFDM-IM symbols with different non-zero supports exhibit sparse structure, which is different from the OFDM frame. In the OFDM frame, almost all subcarriers are required for data transmission. To solve the problem of data detection in OFDM-IM, VAMP based framework is incorporated by considering the symbol sparsity generated from the SAPs. In this section, the concept of VAMP detection and detailed message scheduling is first introduced. Then the VAMP-based data detection with the designed prior aided shrinkage functions is presented.
3.1. VAMP Framework for IM Detection
According to the Equation (
10), the data detection problem is defined as
The problem means to recover the vector
from noisy linear observation
with noise
.
follows Gaussian distribution with the element
, and
is the noise precision. Once the problem is solved, the indices are detected and the symbols are equalized jointly. The detection is executed for each group, and the superscript
g is omitted for simplification.
In the field of digital communications, many problems such as channel estimation, data detection or user detection, are involved as Bayesian estimation problems. The key to the Bayesian message-passing graph is appropriate assumptions about the prior distribution and Gaussian noise. The standard linear regression problem corresponding to the problem in (
12) is given by
where
is the penalty function.
and
is the measurement matrix, which is the diagonal matrix with the channel components. Relate the measurement matrix
with the system model in
Section 2,
.
VAMP is first proposed in [
15] to solve the problem (
13), which has been proved robust to a broader class of large random matrices
compared with conventional AMP method [
9]. Assuming known prior function
and likelihood function
, the posterior function is calculated through Bayesian rule as
where
. According to different estimate criteria, the minimum mean square (MMSE) estimation is
and the maximum a posteriori (MAP) estimation is
.
The basic VAMP framework is given as in Algorithm 1. In Algorithm 1,
is the denoising function which is parameterized by
and
.
is the divergence at
and
. The Onsager term
cancels the correlation between
and
and guarantees the Gaussianity of estimated errors.
Algorithm 1 Vector AMP (SVD Version) |
Input: Measurement , dictionary matrix , maximum iteration number T, denoising function , noise precision Output: Recovered channel vector - 1:
Initialization: , , - 2:
Compute SVD of , and , , - 3:
Compute - 4:
while
do - 5:
- 6:
- 7:
- 8:
- 9:
- 10:
- 11:
- 12:
end while
|
To involve the problem (
12) into the VAMP framework, the factor graph and message scheduling of the VAMP-based detection are shown in
Figure 2.
In
Figure 2, the leftmost node denotes the likelihood function defined as
, where
is the
mth row of the measurement matrix
. The rightmost factor node denotes the prior distribution
. The original variable node
is divided as
and
, which are connected by the Dirac delta function
. Then the joint probability distribution is
. The message passing scheduling in the non-loopy graph includes three parts:
The expectation propagation (EP) is adopted for belief calculation using the Gaussian approximation and moment matching. For each variable node, the marginal function is the product of all impinged messages. We denote the approximate message for as , where , i.e., and . The approximate belief for the node is determined by expectation and covariance with respect to the marginal function , i.e., and .
The message from the variable node to the factor node is . Because all beliefs are approximated as the Gaussian messages, the messages at tth iteration are assumed as and .
Message from factor node to variable node is , where is to integrate over the other connected factors except . is the neighboring node of f.
According to the sum-product rules, the bidirectional messages for each edge are derived as , , and . Due to the equivalent relation in , and . Based on the above analysis, we can achieve the node message for tth iteration as follows:
(1) Calculate
with
, then
(2) Calculate with , then one key step is to define the prior . Without any assumptions, the prior is unknown; therefore, the shrinkage function and the divergence are unknown. The problem is solved in the next subsection.
(3) Calculate the edge message with , then its mean and precision are given by and .
With the forward-and-back message passing, the iterations are converged until the true belief is approximated. The non-loopy factor graph in
Figure 2 provides interpretable message scheduling with ease of implementation.
3.2. Prior-Aided VAMP Detection
Considering the explicit definition of the sparsity of OFDM-IM symbols in the frequency domain, it can be predicted that the posterior estimation of symbols based on the optimal MMSE criterion can be achieved.
has sparse characteristics and non-zero symbols are drawn from the constellation
. For each element in
, the prior of
x is given by
where
is the signal sparsity and
is the size of constellation map
. The constellation point
s has normalized power, i.e.,
.
As shown in step (2), the marginal belief
, which is given by
, where
with mean
and variance
. The MMSE denosing function is element-wise given by
. Expand the denominator as
The numerator is given by
The derivative of MMSE denoiser is required to be calculated, which needs the key term
.
, where the numerator is given by
With the equation from (
18) to (
20), the MMSE denoiser and its divergence are achieved, which are utilized to replace step 5 and step 6 in Algorithm 1.
3.3. Modification for Invalid SAPs and Constellation Detection
After each group is estimated, the invalid SAPs appear like the case in LLR detection. Instead of discarding the non-existent SAPs as wrong results, modifications are required to reduce the performance loss. Firstly, we sort the power amplitudes of in an ascending order and obtain the indices sequence , then the non-zero index group is guaranteed by the k indices with the largest values in . The step is important especially for low SNR cases due to the noise effects greatly blur the boundaries between noise and signal subspace.
However, the member
may be out of the scope of
. To minimize the effects of error detection, the value in the index group is checked to choose the closest one in
. The choice criterion is to minimize the Euclidean distance between
with the Table space, i.e.,
. It is worth noting that symbols
are simultaneously equalized and they have been located on the decision region of the constellation symbols. They are easily recovered using one-step ML symbol detection given by
Besides the VAMP with SVD transform in Algorithm 1, here we present the whole process considering the message passing of
Figure 2, and the VAMP with LMMSE step is shown in Algorithm 2.
Algorithm 2 Vector AMP-based OFDM-IM detection (LMMSE Version) |
Input: Received signal of each group , measurement matrix , Maximum iteration number T, denoising function , noise precision , activated subcarrier number k Output: Recovered index data and equalized symbols - 1:
Initialization: , , - 2:
while
do - 3:
- 4:
- 5:
- 6:
- 7:
- 8:
- 9:
- 10:
- 11:
- 12:
- 13:
end while - 14:
Sort as and estimated indices - 15:
if
then - 16:
Update as with the minimum ED in - 17:
end if - 18:
|