1. Introduction
Transmitted signals are distorted by multi-path propagation in wireless communication. Therefore, in wireless OFDM systems, equalization is popularly performed to compensate for the distorted signals. One of the most important jobs of an equalizer is CFR estimation. The communication specifications in a mobile environment are usually provided with pilots to estimate the CFR at the receiver. There are three common types of pilot signal arrangements: block arrangement, comb arrangement, and scattered arrangement [
1,
2,
3]. In comb-type pilot arrangement, the pilot signals are uniformly distributed in the frequency domain within each OFDM symbol, while in block-type pilot arrangement, the pilot signals are assigned to all subcarriers in a particular OFDM symbol, which is sent periodically in the time domain. In scattered-type pilot arrangement, pilot signals are scattered in the time and frequency domains. The channel estimation is generally divided into two steps at the receiver. First, the CFRs at the pilot subcarriers are estimated. Second, the CFRs at the data subcarriers are estimated using the estimated CFRs at the pilot subcarriers. In the second step, the common estimation methods include 1D or 2D linear interpolation [
4,
5], the DFT-based channel estimation [
6,
7,
8,
9], the LS method [
10,
11], and the compressed sensing method [
12,
13,
14,
15,
16,
17]. Linear interpolation requires the least computation, while the compressed-sensing method requires the most. In recent years, machine learning has also been used for channel estimation [
18,
19,
20,
21,
22]. However, these algorithms [
18,
19,
20,
21,
22] are usually computationally demanding and require a large number of training samples. Refs. [
23,
24,
25] are other related literatures we refer to.
If there is no guard band and the pilot spacing is the integer power of two for the comb-type pilot arrangement, the DFT-based channel estimation is a good estimation method. However, in the present communication specifications, there is a guard band in the spectrum to reduce the interference from neighboring bands. This leads to large errors in the estimated CFRs near the guard band based on the DFT-based channel estimation (i.e., the leakage effect [
6,
7,
8,
9]). In [
6,
7,
8], the CFRs at the virtual pilot subcarriers were estimated with different methods. The method proposed in [
6] inserts extra pilots at the specific subcarriers on both sides of the guard band. Then, the CFRs at the virtual pilot subcarriers are estimated using the inserted pilots. However, as the pilot signal arrangement used by this method is not applicable to any communication specification, it cannot be used in actual communication systems. The optimal linear estimator for leakage suppression is proposed in [
7]. The method [
7] must determine the statistical properties of a channel. However, in practical application, the statistical properties of a channel are unknown and unlikely to be estimated. In [
8], the estimated CFRs at the two pilot subcarriers at both ends of the guard band were used to estimate CFRs at the virtual pilot subcarriers by linear interpolation. However, linear interpolation may not provide a good estimate. The three methods can estimate the overall CFR using the conventional DFT-based channel estimation after estimating the CFRs at the virtual pilot subcarriers. The two estimating channel impulse response (CIR) methods of [
9,
11] are very similar, and their computational complexities are the same.
The method proposed in this paper does not require any statistical properties of a channel or the insertion of extra pilots. This method is applicable to the pilot signals in a comb arrangement. When the number of maximum channel delay points is smaller than or equal to half the number of pilots (including virtual pilots), the CFRs at the virtual pilot subcarriers can be estimated. The method proposed in this paper has lower computational complexity and stored complex numbers than the LS method [
11]. The channel estimation efficiency of our proposed method is still similar to that of the LS method, even if the number of maximum channel delay points is greater than half the number of pilots. Therefore, this paper performed simple system simulations in the time-varying channel to verify the channel estimation efficiency of the proposed method.
The system architecture is outlined in
Section 2 and the proposed method is introduced in
Section 3. The computational complexity is compared with the LS method in
Section 4. Furthermore,
Section 5 shows the simulation and discussion, while
Section 6 shows the conclusions of this paper.
2. Pilot Arrangement and CFR Estimation at Pilot Subcarrier
Figure 1 shows the locations of the pilot of the OFDM system with a comb-type pilot arrangement. Related parameters are described in
Table 1. When the time-varying channel is approximately invariant in an OFDM symbol, the digital CIR corresponding to some OFDM symbols will be set to
h[
n], where
n is the index for time. If the number of maximum channel delay points is
Np/2,
h[
n] = 0 when
n >
Np/2 − 1 or
n < 0. Let the frequency domain data transferred in an OFDM symbol be
X[
k], where
k is the index for frequency. The time-domain signal
x[
n] is derived from the
N-points IDFT of
X[
k] (i.e., the content free of CP (cyclic prefix) in this OFDM symbol). Let the OFDM symbol content received by the receiver after the CP is removed be
y[
n]. When the CP is long enough that two adjacent OFDM symbols will not disturb each other because of the channel,
y[
n] can be expressed as follows:
where
w[
n] is white Gaussian noise and ⊗ is the circular convolution operation. The CFR
H[
k] of this OFDM symbol can be derived from the
N-points DFT of
h[
n]; the mathematical expression is given as follows:
The
N-points DFTs of
x[
n],
y[
n], and
w[
n] are
X[
k],
Y[
k], and
W[
k], respectively. In addition:
where
W[
k] is also white Gaussian noise. The DFT of a circular convolution of two
N-points sequences is the product of their respective DFTs. In an OFDM symbol,
Y[
k] can be expressed as follows:
Afterward, CFRs at the pilot subcarriers are estimated by the LS method. The CFRs estimated at the pilot subcarriers can be written as follows:
All the above variables are described in
Table 2.
3. The Proposed Method Based on the LS Method
The number of maximum delay points of a CIR was assumed to be
Np/2. We divided the pilots into two groups according to their locations as even and odd multiples of ∆
p. From (2), two groups of CFR with adjacent pilots spaced by 2∆
p were first observed, i.e.,
The
Np/2-points of IDFTs of
H[2
k∆
p] and
H[(2
k + 1)∆
p] are shown as (7) and (8), respectively:
(7) and (8) are represented by the following matrices:
where
Gp is the inverse Fourier transform matrix of 0.5
Np × 0.5
Np with (
Gp)
i,k = 2
ej2πik/(Np/2)/
Np and
i,
k = 0, 1, 2, …,
Np/2 − 1. The CIR can be written as (9) and (10), so the two equations are equal as follows:
The CFRs at the virtual even pilot subcarriers were selected as the target of estimation. Equation (11) was changed as follows:
where (·)
H represents the conjugate transposed operation of a matrix; 0.5
Np,
; and
A = 0.5
Np(
Gp)
HDHGp is a 0.5
Np × 0.5
Np matrix. For convenient representation,
A is expressed as follows:
In this case, we assume that
and
are both even and equal and
is a multiple of 4. So
is also even. Otherwise, the matrix dimensions of
A1 to
A9 will vary according to different conditions. As the target of estimation is the CFRs at the virtual even pilot subcarriers, the unknown CFRs at the virtual subcarriers part in the left term
Hp,odd of (12) were removed, and the equation at the righthand side of (12) was divided into two parts: CFRs at virtual and non-virtual pilot subcarriers. Equation (12) can be changed to the following:
where
is the CFR vector at the virtual even pilot subcarriers that we aimed to estimate.
is a 0.5
× 0.5
matrix and
is a 0.5
× 0.5
matrix and
and
in Equation (14) can be estimated using pilot signals.
and
are two known matrices. Afterward, the LS method was used to estimate the CFRs at the virtual even pilot subcarriers, expressed as follows:
where
and
are the CFRs estimated by (4).
There is a problem in (15), where
approaches a singular matrix, so that the operation content value of the inverse matrix is abnormally large. The estimation error induced by noise will be amplified severely [
11], influencing the effectiveness of this method. The corresponding improvement method has been proposed in [
11], that is, to add the regularization parameter before the calculation of the inverse matrix. Thus, Equation (15) was changed to:
where
αe is the regularization parameter, which is a constant;
I is a
identity matrix; and
is the final result estimated with the regularization parameter. However, (16) is not directly used for operation in actual hardware realization. The right side of Equation (16) was expanded as follows:
and
can be worked out in advance and stored in the memory. Both of them are
matrices, and each has
complex numbers to be stored in the memory. The number of multiplication operations required by (17) is much smaller than that required by (16).
Afterward, the CFRs at the virtual odd pilot subcarriers were estimated. Equation (11) was changed as in Equation (12) to obtain the following equation:
where
B is a 0.5
Np × 0.5
Np matrix and
B = 0.5
Np(
Gp)
HDGp. For convenient representation,
B is expressed as follows:
As the target of estimation is the CFRs at the virtual odd pilot subcarriers, the unknown CFRs at the virtual subcarriers part in the left term
Hp,even of (18) were removed, and the equation at the right-hand side of (18) was divided into two parts: CFRs at virtual and non-virtual pilot subcarriers. Equation (18) can be changed to the following:
where
is the CFR vector at the virtual odd pilot subcarriers that we aimed to estimate.
is a 0.5
× 0.5
matrix and
is a 0.5
× 0.5
matrix and
and
are two known matrices. Afterward, the LS method was used to estimate the CFRs at virtual odd pilot subcarriers, expressed as follows:
where
and
are the CFRs estimated by (4). The regularization parameter
αo is used to overcome the problem that
approaches a singular matrix. Equation (20) was not directly used in actual hardware implementation. The right side of Equation (20) was expanded as follows:
and
can be calculated in advance and stored in the memory. Both of them are
matrices, and each has
complex numbers to be stored in the memory. The number of multiplication operations required by (21) is much smaller than that required by (20).
According to the results of (4), (17), and (21), the CFRs vector at the pilot and virtual pilot subcarriers can be obtained. The CIR of the Np points can be calculated by the Np-points IFFT of . According to the assumptions of this paper, the number of maximum delay points of a CIR is smaller than or equal to half the number of pilots (including virtual pilot) (i.e., 0.5Np points). Therefore, the values of 0.5Np points of the second half of can be regarded as the estimation error and discarded; that is, let the values of 0.5Np points of the second half of be zero to obtain . The complete CFR can be obtained by N-points DFT of .
The whole procedure of the proposed method in every OFDM symbol is organized as follows:
- Step 1.
The CFRs at pilot subcarriers are estimated by (4).
- Step 2.
The CFRs at virtual even and odd pilot subcarriers are estimated by (17) and (21) respectively.
- Step 3.
According to step 1 and step 2, the CFRs vector at pilot and virtual pilot subcarriers can be obtained.
- Step 4.
The CIR of the Np points can be calculated by the Np-points IFFT of .
- Step 5.
Set the values of 0.5Np points of the second half of be zero to obtain .
- Step 6.
The complete CFR can be obtained by N-points DFT of .
4. Comparison of Computational Complexity with the LS Method
The equation for the estimation of a CIR by the LS method in [
11] is shown as follows:
where
is an
partial Fourier transform matrix with
,
k ∈ pilot,
n = 0, 1, …,
− 1. Additionally,
I is an
identity matrix. When the relevant parameters are determined,
of (22) can be calculated in advance and stored in the memory. The acceptable channel maximum delay of the LS method is
. The channel maximum delay of our method is less than or equal to
Np/2. For a fair comparison of our method with the LS method [
11], let the channel maximum delay be equal to
Np/2. Therefore, the content of
in (22) will change, and the matrix dimension changes to 0.5
and has 0.5
complex numbers to be stored in the memory. The computational complexity of Equation (22) is
complex multiplications.
Consider the number of multiplications after CFRs estimation at the pilot subcarriers until the CIR estimation. The computation of our method was divided into two parts. Part 1 uses (17) and (21) to estimate the CFRs at the virtual pilot subcarriers. The number of multiplications of this part is
. Part 2 performs the
Np-points IFFT to obtain the CIR. The number of multiplications of this part is (
Np/2)log
2(
Np) [
26]. Therefore, the computation of the proposed method requires
+ (
Np/2)log
2(
Np) complex multiplications, which is much fewer than
in general cases. Comparisons with relevant parameters are described in subsequent simulations.
5. Simulation and Discussion
First, the simulation parameters were set up in
Table 3. The signal was QPSK modulation and the speed was 20 km/h.
Table 4 shows the typical urban time-varying channel model [
27]. The three power spectral density function models were Jakes, Gauss I, and Gauss II (see COST 207 specification [
27]). Forty channels were generated randomly with the channel model for simulation with different SNRs. Fifty OFDM symbols were transferred in one simulation. In terms of regularization parameters, the LS method [
11] selects 0.01, while the method proposed in this paper selects
αe =
αo = 0.02. The dimension of the estimated CIR for the LS method [
10] is also set to
Np/2. The problem of the matrix approaching singularity also exists in [
10], so we use a regularization parameter that is set to 0.00001 to avoid this issue. The regularization parameter selection process refers to [
28]. To ensure a fair comparison of the conventional DFT-based channel estimation with other methods, we only use first 0.5
Np points of the CIR estimated by conventional DFT-based channel estimation. All simulations were performed using MATLAB software.
Figure 2 is the simulated diagram of SNR to the bit error rate (BER). It is observed that the curves of the LS methods [
10,
11] and the method proposed in this paper are completely overlapped under the appropriate regularization parameter, meaning the BERs are very close. The BER curves of conventional DFT-based channel estimation and [
8] are very close. This shows that the leakage effect was not improved by the method proposed in [
8]. The conventional DFT-based channel estimation and [
8] have lower computational complexity than the proposed method in this paper. However, the BER simulation showed that the proposed method has better effectiveness than the conventional DFT-based channel estimation and [
8]. At low SNR, the effect of noise is higher than the leakage for the conventional DFT-based channel estimation. Therefore, the three curves of the conventional DFT, LS, and proposed methods are very close at low SNR. At high SNR, the effect of leakage is higher than the noise for the conventional DFT-based channel estimation. Therefore, the BER curve of the conventional DFT-based method is worse than the BER curves of the LS methods and the proposed method at high SNR. Linear interpolation has the lowest computational complexity but the worst BER curve.
Figure 3 and
Figure 4 show the simulation results with speeds equal to 100 km/h and 200 km/h, respectively. At low SNR, the six curves did not change from
Figure 2,
Figure 3 and
Figure 4 as the effect of noise was higher than speed. At high SNR, the six curves from
Figure 2,
Figure 3 and
Figure 4 all deteriorated with increasing speed as the effect of speed was higher than noise. Since the channel changes were extraordinarily slow, it can be assumed that the channel is approximately invariant within an OFDM symbol when speed is very low. The channel was no longer approximately invariant within an OFDM symbol when the speed is very high. Therefore, the simulation results at high speed are worse than the simulation results at low speed when SNR is high. The BER curves of the LS methods [
10,
11] and the method proposed in this paper are still overlapped in
Figure 3 and
Figure 4.
Consider the number of multiplications after CFRs at pilot subcarriers estimation until the CIR estimation. The number of complex multiplications of the LS method [
11] is
. There are 6848 complex multiplications after the relevant parameters are substituted in. Meanwhile, the number of complex multiplications of the method proposed in this paper is
+ (
Np/2)log
2(
Np). There are 2695 complex multiplications after the relevant parameters are substituted in. Both the LS method [
11] and the proposed method have some complex numbers that need to be stored in the memory. The numbers of complex numbers that need to be stored in the memory are
and
, respectively. In this example,
= 6848 and
= 2247.
Table 5 shows the comparisons of different methods on computational complexity and memory resource usage.
It is obvious that the method proposed in this paper has lower computational complexity and stored complex numbers than the LS method [
9,
11], while both have the same BER from simulations. Therefore, the method proposed in this paper is better than the LS method [
9,
11].
Let us consider a bad urban time-varying channel model [
27] as shown in
Table 6. The number of maximum channel delay points is equal to 100 and greater than
Np/2. Therefore, the dimension of the estimated CIR for both LS methods [
10,
11] is set to the maximum value of the dimension, which is
. But in this case, the matrices in both LS methods [
10,
11] are getting closer and closer to singular. The regularization parameters of [
10,
11] are set to 0.00005 and 0.05, respectively. For fair comparison, the proposed method, the DFT method, and [
8] also take the first
points of the estimated channel impulse response as the final estimated channel impulse response.
Figure 5 shows the simulation results with speed equal to 20 km/h. The six curves in
Figure 5 are all becoming worse than those in
Figure 2. The numbers of maximum channel delay points in bad and typical urban channel models are 100 and 50, respectively. The larger the maximum channel delay, the worse the channel estimation efficiency of linear interpolation will be. Therefore, the linear interpolation curve in
Figure 5 is obviously much worse than that in
Figure 2. Since the matrices are getting closer and closer to singular in the two LS methods, although the regularization parameter is used, it still makes the two LS curves in
Figure 5 worse than those in
Figure 2. The estimated CFRs at the virtual pilot subcarriers with our proposed method will be inaccurate when the number of maximum channel delay points is larger than
Np/2. Therefore, the curve of our proposed method in
Figure 5 becomes worse than that in
Figure 2. However, the curves of the two LS methods and our proposed method still overlap in
Figure 5. This shows the channel estimation efficiency of our proposed method is still similar to those of the LS methods [
10,
11], even if the number of maximum channel delay points is greater than
Np/2.