In particular, the estimation of doubly-dispersive channels
is a very important aspect for future multicarrier systems
especially when TF symbols are precoded. Since the provi-
sion of an accurate channel state information (CSI) and the
usage of an appropriate equalizer is essential to enable high
TF diversity gains. Vehicular channels are considered to be
doubly-dispersive, underspread, and often also to be sparse in
the continuous DD domain following the wide-sense stationary
uncorrelated scattering (WSSUS) model [18]. A channel is
underspread if all delay shifts and Doppler shifts are contained
within a small region, i.e., both are relatively small. The
channel is sparse when only a few point-like scatterers in the
continuous DD domain exist. For pulse-shaped multicarrier
filterbanks, the inherent sparsity of the channel cannot be
harnessed using any form of discrete Fourier transform (DFT)
for channel estimation in the discrete DD domain [14], [19],
[20]. A common way to estimate the channel is to get the least-
squares (LS) estimator from the pilot samples and to smooth
them by means of the Wiener filtering in the discrete DD
domain, which is commonly referred to as linear minimum
mean square error (LMMSE) estimator or Markov estimator
[21], [22]. This approach however suffers under leakage effects
[20], [23]. Leakage effects are caused by the presence of both
fractional Doppler shifts and fractional delay shifts which are
not consistent with the discrete nature of the DSFT [14].
To cope with leakage effects and to promote sparsity,
more complex estimation schemes are commonly followed.
In this scope, compressed sensing or even super resolution
are possible schemes, see for example [23], [24], respectively.
In [25], Rasheed et al. propose a compressed sensing based
algorithm using orthogonal matching and modified subspace
pursuit to estimate the time-varying channels. A framework
for sparse Bayesian learning with Laplace priors and a new
piloting scheme has been introduced by Zhao et al. in [26],
where they consider fractional Doppler shifts but not fractional
delay shifts. An off-grid sparse signal recovery to estimate
the original channel rather than the effective discrete channel
in the DD domain is proposed in [27]. In [28], an iterative
optimization method is presented by Liu et. al., where a
message passing signal recovery algorithm is utilized for
channel estimation which takes fractional Doppler shifts but
not fractional delay shifts into account. The listed schemes are
rather complex, require high computing power, and consider
longer time intervals, e.g., are computed adaptively over multi-
ple frames, which does not suite well to URLLC in the context
of rapidly changing vehicular channels, as it is known that
the WSSUS assumption only holds for a limited duration and
bandwidth [29]. This makes channel estimation challenging
and requires channel estimation on a per frame basis [10].
Computationally complex and iterative optimization methods
are therefore not considered in the presented paper.
Focusing on low-complexity estimators for URLLC, a com-
mon choice is the estimation of the channel main diagonal
(CMD) on a per frame basis. This can be done by using
an LMMSE estimator which however suffers from leakage
effects. In this paper, we propose a novel CMD estimator
in the TF domain in contrast to the estimation in the DD
domain used for OTFS and DFT based schemes for OFDM.
We place pilot symbols in the TF domain to enable higher
flexibility and reduced overhead for pilot signaling. However,
the pilot and data symbols still need to be properly arranged
within a rectangular frame. To apply fast orthogonal precoding
transformations, we typically require the input dimension to be
to the power of two which equals the number of data symbols.
Therefore, the placement of the pilot symbols is not obvious.
To control the number and position of the pilot symbols,
we propose an algorithm and a so called accordion pilot
placement to place pilots in between the precoded symbols
in the TF domain. The main contributions of this paper can
be summarized as follows:
•We study pulse-shaped multicarrier systems with linear
precoding for URLLC over doubly-dispersive channels,
•we numerically compare different linear precoding trans-
formations,
•we propose a novel smoothness optimized estimation
scheme of the CMD coefficients which minimizes the
energy of the discrete Hessian and takes the ratio between
the delay spread and Doppler spread, the self-interference
power, and receiver noise into account, and
•we introduce a pilot placement scheme, i.e., accordion
pilot placement, which enables a smooth control of the
number and position of the pilot symbols.
A. Paper Organization
In Section II, the Gabor signaling and doubly-dispersive
channel model is introduced. Linear precoding transforms and
their diversity gain are discussed in Section III. In Section
IV, we detail channel estimation, leakage effects, equalization,
data recovery, and the proposed channel estimation scheme.
The accordion pilot placement is presented in Section V.
In Section VI, we show our numerical results. Finally, we
summarize our conclusions in Section VII.
B. Notational Remarks
Random variable vectors, 2D-arrays and matrices are de-
noted with bold letters. Superscripts (·)∗and (·)Hdenote the
complex conjugate and the Hermitian transpose, respectively.
Let ∗denote the non-cyclic 2D convolution which only returns
the valid part. The column-wise vectorization operator, the
absolute value, the euclidean norm, and the Frobenius norm is
denoted as vec{·},|·|,k·k2, and k·kF, respectively. We denote
δ(·)as the Dirac distribution, as the Hadamard product,
E{·} as expectation operation, and j2=−1. We denote the
indices of down-converted received signal by (¯
·).
II. SYSTEM MODEL
In this section, we introduce the system model which
includes the doubly-dispersive channel and the input-output
mapping of the information resources. We use a time-
continuous Gabor (Weyl-Heisenberg) signaling to derive a dis-
crete system model for the pulse-shaped multicarrier scheme.
We define the Gabor grid Λ = FZM×TZNwith frequency