#### Abstract

The problem of channel estimation for multicarrier communications is addressed. We focus on systems employing the Discrete Cosine Transform Type-I (DCT1) even at both the transmitter and the receiver, presenting an algorithm which achieves an accurate estimation of symmetric channel filters using only a small number of training symbols. The solution is obtained by using either matrix inversion or compressed sensing algorithms. We provide the theoretical results which guarantee the validity of the proposed technique for the DCT1. Numerical simulations illustrate the good behaviour of the proposed algorithm.

#### 1. Introduction

In wireless communications, the channel filter is usually time-varying; for this reason, it is necessary to estimate the channel filter from time to time. To this aim, some training symbols (i.e., symbols known both by the transmitter and by the receiver) are typically used. In this way, when the training symbols are transmitted by the channel, the received signal is used to extract the information about the channel filter. Some well-known techniques for channel estimation are based on the Discrete Fourier Transform (DFT); in this case, the training symbols are OFDM waveforms.

Additionally, if the channel filter is sparse (i.e., containing only a small amount of nonzero coefficients), then compressed sensing techniques can be applied. Compressed sensing (CS) algorithms approximate the sparsest solution to a linear system [1]. This is very useful when the solution depends on a small number of degrees of freedom and only a few measurements of the vector are observed. For this reason, in the last few years CS algorithms have been applied to a wide variety of scenarios in communications: cognitive radio, radar, antenna arrays, multicarrier communications, and so forth. When CS is applied to channel estimation problems, it is usually denoted as compressed channel sensing (CCS). Several CCS algorithms have been proposed in the literature for different types of channels arising in communication problems, such as ultrawideband channels, underwater acoustic communications, or multipath channels [2–5]. Most of these techniques are based on DFTs or spread spectrum signals.

In this work, we consider a multicarrier communications system that is based on the Discrete Cosine Transform Type-I (DCT1) even instead of the standard DFT. Some Discrete Cosine Transforms have been widely used in the context of multicarrier modulation (MCM), as an alternative to the DFT, due to their good properties (e.g., good performance under carrier frequency offset) [6–13]. In particular, in a very recent work [14] the DCT1 is applied for MCM communications. The main advantages of the DCT1 are as follows:(i)The inverse of the DCT1 is the same transform DCT1, up to a scaling factor; so we can use the same transform at both the transmitter and the receiver [15].(ii)The convolution of two vectors is transformed by DCT1 into a pointwise product of their transforms (under some symmetry conditions on the vectors) [15, 16]. This is analogous to the circular convolution property of the DFT. This is a key property for signal reconstruction in MCM communications [14].

For these reasons, we investigate the use of DCT1 for channel estimation; in particular, we address the problem of estimation of whole-point symmetric (WS) channels by means of CS techniques in the DCT1 transform domain. The strategy consists of using only a few training symbols, which are transmitted through the channel, and reconstructing the impulse response of the filter in the receiver by using the same small number of measurements. Thus, the economy of the data can be exploited by CS algorithms, which are able to provide sparse filters.

In this work we will provide not only a new estimation procedure but also the training signals valid for our algorithm, and we will show that this technique is both simple and theoretically correct. These are the main contributions of this paper. Numerical simulations also illustrate the effectiveness of our results.

The paper is organized as follows. Firstly, in Section 2 we recall the general channel estimation problem. Secondly, in Section 3 the DCT1 is introduced and we obtain new important properties of this transform. Then, the proposed procedure is presented in Section 4, where its theoretical justification is also provided. Section 6 contains some numerical examples that illustrate the behaviour of our algorithm. Finally, we highlight the main contributions of this work in Section 7.

#### 2. The Channel Estimation Problem

Let us consider a multicarrier modulation communications system that performs an inverse transform in the transmitter and a direct transform in the receiver , as shown in Figure 1. Let us consider also a channel with the following impulse response: The transmission of an information symbol through this channel results in a received symbol , such that where is a term related to the additive noise.

In multicarrier systems, in order to eliminate interblock interference, we often add to the original symbol a left prefix and a right suffix , both of length :In matrix form, the received data symbol is given bywhere is the Toeplitz matrix of size defined by the filter: It is easy to see [8] that this received symbol can be written asTherefore, if we apply both the discrete transformations and in order to diagonalize , then and denoting , and , we get

Now, the question is if we know the training symbol and its corresponding received symbol , is it possible to estimate ? The answer is* yes*, whenever there is an invertible which transforms into the elements of the diagonal matrix . Hence, in the absence of noise, it suffices to compute ; but in the presence of noise we only obtain an estimation in the transform domain: Now we can recover the estimated filter as Of course, this estimation would be exact in absence of noise. See Figure 1 for a general diagram of the channel estimation problem.

The existence of such transform is a condition usually met in practice. For example, in OFDM systems the signal is extended by appending a cyclic prefix or suffix, so that the equivalent matrix is circulant and diagonalized by the DFT transform. The diagonal matrix contains the eigenvalues of , which in addition form the vector DFT(). Hence, is estimated simply by applying an inverse DFT.

As OFDM systems present poor behaviour under carrier frequency offset, other multicarrier modulation (MCM) techniques have been investigated, which are related to other transformations different from DFT. Among them, the eight types of Discrete Cosine Transforms (DCTs) have been studied in the literature, and for each one of them the corresponding extension technique has been proposed [6, 8, 9, 14]. These works focus on MCM signal reconstruction, and they provide good results due to the good properties of the DCTs. However, the channel estimation stage is essential in order to implement a DCT-MCM system in practice. For this reason, in this paper we apply the DCT1 to the channel estimation problem for the first time.

#### 3. The Discrete Cosine Transform Type-I (DCT1) Even

The DCT1 even of an -length signal is given by the matrix , whose th element is defined by where This is the definition of given in [15], except for the normalization factor ; it has been introduced here in order to ensure the involution property, , which simplifies the numerical calculations. In this way, the direct and inverse DCT1 transforms are identical.

The first contribution of this work is the demonstration of the following theorem regarding the invertibility of some submatrices of the DCT1 matrix. This is a key property which guarantees that the channel filter can be obtained by means of a small amount of received data; this will be applied in the following section, when using compressed sensing techniques. Let us now state and prove this important property.

Theorem 1. *Any submatrix of , whose columns have been extracted from the first columns of , is invertible.*

*Proof. *The submatrix formed by the first columns of is Let us consider any submatrix of this matrix; our aim is to show that is invertible. Notice that its rows can be indexed as , so the entries of areTo show that is invertible, it suffices to prove that the unique vector such that is Let be such vector; the condition is rewritten asBy defining we can rewrite the latter expression asOur aim is to prove that the numbers which fulfill (16) are necessarily null; it is equivalent to the fact that , , finishing the proof.

To this aim, let us now introduce the auxiliary self-reciprocal polynomial of degree : Notice that, for any , we have thatthus, if is a nonzero root of , then also is a root of

Our strategy is to prove that has roots, say, more than its degree ; if this occurs, then must be the null polynomial, and all its coefficients necessarily are 0, so , , and the claim follows. Let us find some roots of :(i)By denoting the complex numbers it is easy to see that where we have used (16). Hence, are roots of ; note that there are different numbers because their arguments lie in since(ii)As already mentioned, for any , also is a root of : So there are different roots of whose arguments lie in (iii)The union of the set of roots of (19) and (22) provide a total amount of different roots if and In this case, has more roots than its degree, so is the null polynomial and the claim holds. In case , the corresponding root has been counted twice; the same happens if , because the root would appear twice. In these cases, we can only guarantee that there are different roots but it is easy to see that any self-reciprocal polynomial of even degree satisfies the following property: if (or is a root of , then it is a root of multiplicity at least 2. This implies that, in our case, (or is a double root, so has at least roots, concluding the proof.

#### 4. Channel Estimation in DCT1 Multicarrier Systems

Let us assume that the channel filter presents whole-point (WS) symmetry: As we have already mentioned before, the interblock interference is eliminated by introducing as redundancy a left prefix and also a right suffix , both of length , into each data symbol to be transmitted. In order to apply DCT1, it is proved in [9, 14] that it suffices to consider the extended block in (3) with the choice of prefix and suffix as follows: which means that we apply a whole-point symmetry (WS) on the left and on the right sides of the original symbol. Figure 2 illustrates an example of the WS symmetric extension of . The received vector (4) is then

It is proved in [14] that the corresponding can be perfectly diagonalized via the DCT1: and the diagonal elements of eigenvalues of are themselves the DCT1 transform of the vector :where stands for the half-right filter of , padded with zeroes:

Thus, we have been able to find an easy solution to the channel estimation problem in DCT1 MCM communication systems. Following the general statement of the problem given in Section 2, we denote , , and and get the schemewhere .

As the components of the training signal are the symbols that can be stored in memory, from this equation we simply obtain an estimation of , for any component such that :In other words, if all the 1-tap filters are nonzero (, then we compute the -length vector , which gives a perfect estimation of in absence of noise, and we can straightforwardly obtain the symmetric channel filter

#### 5. Compressed Channel Sensing for DCT1

Now, the question is what can we do if a component of is null? In effect, this is a situation very common in practice. Indeed, we would like to have only a few training symbols and many null components Moreover, how can we obtain an estimated vector which is* sparse*? The answer to these two questions is given by two facts: on the one hand, we can apply compressed sensing techniques; on the other hand, the DCT1 matrix presents a key property which guarantees that the components of the sparse vector can be obtained by knowing only components of the vector

Let us explain this idea in detail. If there are only nonzero symbols , corresponding to the components , then only components of the vector are defined, by means of (29). As has length , it is impossible to recover its components, but at least we can try to estimate of them by means of the computed components of the estimated vector . Let us use the same notation to define the -length vector which contains the known components of (29). Then, we need to find the sparse vector that minimizes the norm where denotes the submatrix of formed by its corresponding rows .

Besides, we can exploit the structure of , which has final zeroes, so we can write where we have denoted , and stands for the submatrix of containing the* first * columns of and the corresponding rows In this way, we have the following minimization problem:where vectors and have length

Compressed sensing techniques show that it is possible to achieve the sparsest vector if its sparsity order is (the spark of a square matrix is its rank plus 1). As we have proved in Section 3, our Theorem 1 guarantees that the rank of any square submatrix of the first columns of is maximum, so has maximum rank. This means that we can reconstruct sparse filters of sparsity order In practice, it is possible to recover vector in two ways:(i)As is invertible, we can simply define so as to get null error ( This is true in absence of noise, but the drawback in practice is that we may obtain a nonsparse vector (ii)Alternatively, as has maximum rank , we can apply some well-known algorithms used in CS scenarios (e.g., Lasso techniques, OMP, and CoSamp [5]) in order to find the sparse solution of the problem. These algorithms converge when the matrix of the linear system satisfies either the Restricted Isometry Property (RIP) or a weaker property regarding the* coherence* of . Luckily, we can also guarantee that this good property is fulfilled by the DCT1, by placing symbols at equally spaced positions. This procedure has been applied in our simulations, and this clear advantage is the reason why the DCT1 matrix performs well for compressed sensing techniques in the simulations.

* Summary of the procedure is as follows*:(1)Choose a training signal of length , and compute .(2)Apply a whole-point symmetry of length at both edges of so as to obtain of length .(3)Transmit through the channel.(4)Take the central components of the received vector, which form .(5)Apply the DCT1 block: .(6)Compute .(7)In case for all , obtain which is the desired estimation of the half-right filter (8)In case some components of are null, and at least components of are nonzero, find the solution of the problem that can be solved via CS techniques or simply defining (9)In any case, from by WS symmetry we get the estimated filter channel

#### 6. Numerical Results

In this section, we analyse the behaviour of the proposed compressed channel sensing (CCS) scheme by testing it on three channels: a fixed simple channel of length (i.e., ), a more challenging fixed nonminimum phase channel of length (i.e., ), and a perturbed symmetric version of the ITU-T M.1225 pedestrian channel A [17]. In all cases, a sparse signal is constructed in the DCT1 domain by setting if (for ) and otherwise. Hence, is a -sparse vector containing only nonnull elements uniformly distributed, as stated in the previous section. For instance, when we have and the nonnull elements are only , , , and , whereas for we have and the nonnull elements are , , , , , and . In this way, we are truly performing a compressed sensing of the channel, since we are only exploring certain elements (which correspond to particular frequencies) in the transformed domain.

The inverse DCT1 is then performed and the time-domain transmitted vector is passed through the filter with symmetric impulse response . Then, zero-mean additive white Gaussian noise (AWGN) with variance is added, obtaining the received vector . The length DCT1 of the central elements of this vector is now computed, resulting in . Finally, the elements of corresponding to the nonnull positions of are extracted and a length inverse DCT1 is performed on them to estimate the right-half of the channel’s impulse response. The rest of the channel is reconstructed exploiting its symmetry. The performance measure used is the reconstruction signal to noise ratio (SNR), where , , and is the channel’s length.

As a first example, we select the following channel: We set the length of the DCT1 to and check the behaviour of the CCS scheme as the channel’s SNR increases from −10 dB to 30 dB using only 4 training pilots. simulations are performed for each SNR.

Figure 3, which displays the reconstruction SNR as a function of the channel’s SNR, shows that an increasingly accurate estimation of the channel can be obtained as the SNR increases, even by using only pilots. It can be seen that the reconstruction SNR increases linearly as the signal power to noise ratio increases. Indeed, the following relationship can be established: where , with , and in this case. This shows that an increase in SNR of 7.78 dB in the reconstruction is obtained with respect to the channel’s SNR.

Figure 4 shows three examples of the estimated channel’s impulse and frequency responses for three signal to noise ratios: SNR = 0 dB, SNR = 10 dB, and SNR = 20 dB. The channel’s impulse response is displayed on the left hand side (central dot with the true values and bar spanning the range between the minimum and maximum recovered values), whereas the right hand side shows the channel’s frequency response (true value in black line and shaded area showing the range between maximum and minimum values). Note the substantial decrease in the variation of the coefficients of the channel’s impulse response as the SNR increases (indeed, for SNR = 20 dB the bars cannot be appreciated, since the recovered coefficients are always virtually identical to the true coefficients) and the corresponding improvement in the estimation of the channel’s frequency response (with a decrease in the shaded area).

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

We have also tested the effect of the number of subcarriers, , by using for (i.e., ). The result, displayed in Figure 5, shows that the value of is irrelevant (in terms of accuracy of the reconstructed channel), as long as is an integer number and the pilot carriers can be uniformly distributed (as it happens in this case for ). When is not an integer number, the pilots cannot be uniformly distributed and an approximation error is obtained (as seen in the cases , , , and ). However, this error decreases as increases and can be completely avoided by zero-padding the channel’s impulse response until is integer.

As a second example, we consider a length non-minimum-phase channel: We set again the length of the DCT1 to and check the behaviour of the CCS scheme as the channel’s SNR increases from −10 dB to 30 dB using only 6 training pilots. simulations are performed for each SNR. However, even though this channel is much more challenging than the previous one, similar results are obtained in terms of the reconstruction error. Indeed, (39) is also valid in this case and the reconstruction SNR versus channel’s SNR curve for this channel (not shown) is virtually identical to Figure 3. In fact, we have also tested several other (both minimum and nonminimum phase) channels and this result seems to apply to all of them. Figure 6 shows three examples of the estimated channel’s impulse and frequency responses for three signal to noise ratios: SNR = 0 dB, SNR = 10 dB, and SNR = 20 dB.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

Finally, we test our approach on a perturbed symmetric version of the ITU-T M.1225 pedestrian channel A. The pedestrian channel A was generated using Matlab’s stdchan function using a carrier frequency GHz, a sampling period ns, and a length . The resulting channel’s impulse response, , is very sparse, since it typically has only 3 nonnull coefficients. The symmetric channel’s impulse response is constructed as , so its length is and we have . Then, the coefficients of are perturbed by adding independent white Gaussian noise samples with variance to each of them in order to analyze the effect of the lack of symmetry, typical of real-world channels. The results are shown in Figure 7: a small lack of symmetry (e.g., or ) practically has no effect; a moderate amount (e.g., ) lowers the performance but still provides a good estimate of the channel (with a reconstruction SNR around 30 dB); a large lack of symmetry (e.g., ) results in a low reconstruction SNR (around 10 dB), as the reconstructed channel is approximately equal to the symmetric part of the true channel. This highlights the limitations of our approach but also its potential in many approximately symmetric real-world channels (e.g., channels with a large central coefficient and small not completely symmetric coefficients around it).

#### 7. Conclusions

In this work, we have presented a general procedure for the estimation of any symmetric channel filter for multicarrier communication systems based on the Discrete Cosine Transform Type-I (DCT1) even. For any training signal transmitted through the channel, at the receiver, we show how to take into account the information of the training symbol so as to estimate the channel filter. The main contribution of this work is that it is possible to estimate the channel filter with a small amount of training signals, just knowing a small amount of the received samples, and regardless of the location of these samples. This is an important consequence of the good properties of the DCT1 matrix that have been also proved here for the first time. Thus, our proposed procedure with the DCT1 formulation meets the conditions that guarantee perfect estimation of the channel filter in absence of noise, whereas in noisy scenarios a very good estimation can also be achieved. We have also designed specific sparse training signals for our DCT1 procedure and showed that it can also be applied to channels whose impulse response is only approximately symmetric with good results. Future research lines include extending these procedure to nonsymmetric channels.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work has been partially supported by the Spanish Ministry of Economy and Competitiveness through Project TEC2012-38058-C03-01. David Luengo has also been supported by the BBVA Foundation through Project MG-FIAR (“I Convocatoria de Ayudas Fundación BBVA a Investigadores, Innovadores y Creadores Culturales”). All the authors are members of the UPM TACA Research Group and also thank the Universidad Politécnica de Madrid for its support.