Abstract
Accurate channel state information (CSI) is necessary at receiver for coherent detection in amplify-and-forward (AF) cooperative communication systems. To estimate the channel, traditional methods, that is, least squares (LS) and least absolute shrinkage and selection operator (LASSO), are based on assumptions of either dense channel or global sparse channel. However, LS-based linear method neglects the inherent sparse structure information while LASSO-based sparse channel method cannot take full advantage of the prior information. Based on the partial sparse assumption of the cooperative channel model, we propose an improved channel estimation method with partial sparse constraint. At first, by using sparse decomposition theory, channel estimation is formulated as a compressive sensing problem. Secondly, the cooperative channel is reconstructed by LASSO with partial sparse constraint. Finally, numerical simulations are carried out to confirm the superiority of proposed methods over global sparse channel estimation methods.
1. Introduction
Relay-based cooperative communication [1–6] has been studied in the last decade due to its capability of enhancing the transmission range and providing the spatial diversity for single-antenna receivers by employing the relay nodes as virtual antennas [7–9]. A typical example of cooperative communication system is shown in Figure 1, where source transmits signal to destination with the help of relay . It is well known that utilizing multiple-inputs multiple-outputs (MIMO) transmission can boost the channel capacity [10, 11] in broadband communication systems. In addition, diversity techniques in MIMO system could mitigate selective fading and hence improve the quality of service (QoS) [12, 13]. However, it poses a practical challenge to integrate multiple antennas onto a small handheld mobile terminal. To deal with the limitation, one could choose relay-based cooperation networks which have been investigated in last decade [1–5]. The main reason is that the diversity from relay nodes existing in the network could be exploited, where relay can either be provided by operators or be obtained from cooperating mobile terminals of other users.
In the relay-based cooperative communication system (CCS), data transmission is usually divided into two time slots. During the first time slot, the source broadcasts its information to both relay and destination. During the second time slot, the relay could select different protocols and then transmit the signal to the destination. Usually, there are two kinds of protocols in cooperative communication systems; one is to amplify the received signal at the relay and forward it to the destination, which is termed as amplify-and-forward (AF); the second is to decode the received signal, modulate it again, and then retransmit to destination, which is often termed as decode-and-forward (DF). Due to coherent detection in these systems, accurate channel state information (CSI) is required at the destination (for AF) or at both relay and destination (for DF). About DF cooperative communication systems, the channel estimation methods could be borrowed from point-to-point (P2P) communication systems directly. However, extra channel estimation will increase the computational complexity which is a burden at relay, and broadcasting the estimated channel information will result in further interference at the destination. On the other hand, AF cooperative communication technique can avoid this disadvantage and we focus on AF CCS in this study.
As shown in Figure 1, the cooperative channel consists of a direct link and a cascaded one . Based on the assumption of dense multipath, linear channel estimation for the relay-based AF cooperative networks has been proposed [4]. Even though the proposed method can achieve lower-bound performance, low spectrum efficiency is unavoidable since the training sequence takes a large amount of the bandwidth. Hence, one method to improve the spectral efficiency is by reducing the length of training sequence for channel estimation.
As channel measurement techniques improved in the last decade, broadband wireless channels have been confirmed to exhibit inherent sparse or cluster-sparse structure in delay spread. If we can take advantage of the sparse prior information, then the spectral efficiency can be improved. In allusion to point-to-point (P2P) communication systems, efficient sparse channel estimation methods [14–16] have been proposed. To improve the spectral efficient and/or estimation performance in CCS, we have studied channel estimation in CCS and proposed an effective sparse channel estimation method [17].
However, the cascaded channel () may be no longer sparse due to the linear convolutional operation [18]. Relay-based cooperative communication not only reduces the transmission range but also improves the received signal-to-noise ratio (SNR) when comparing with the P2P communication systems. Unlike the previous research under global sparse hypothesis, the cooperative channel consists of two parts: sparse part () and dense part (). A simple example is shown in Figure 2. As a result, our previous method will be degraded since it cannot fully take advantage of the prior information. Unlike the previous method, in this paper, we propose a partial sparse channel estimation method by using LASSO [19] (PEL) to improve the performance. Based on this idea, an improved partial sparse channel estimation by using LASSO (IEL) is proposed by utilizing both partial sparse constraint and global sparse constraint. On the one hand, the partial sparse constraint can improve estimation performance in the cooperative communication systems. On the other hand, the global sparse constraint can mitigate noise interference under low SNR regime. To confirm the effectiveness of the two proposed methods, we give various numerical simulation results in Section 4.
Section 2 introduces the system model and problem formulation. In Section 3, two improved channel estimation methods are proposed. The first method is the improved channel estimation method by using partial sparse constraint, and the second one is an improved channel estimation method by using both partial and global sparse constraint. In Section 4, we give some representative numerical simulation results and related discussions. Concluding remarks are presented in Section 5.
Notations. In this paper, we use boldface lower case letters to denote vectors and boldface capital letters to denote matrices. represents the complex Gaussian random variable. stands for the expectation operation. , denote transpose and conjugated transpose operations of , respectively. is the number of nonzero elements of , and is the Euclidean norm of . represents a diagonal matrix whose diagonal entries are from vector . denotes the convolution of two vectors and .
2. System Model
Consider a multipath fading AF CCS where the source sends data to destination with the help of relay as shown in Figure 1. The three terminals are assumed to be equipped with a single antenna each. , , and denote the impulse response of the frequency selective fading channel vectors between three links , , and , respectively. The three channel vectors are assumed to be zero-mean circularly symmetric complex Gaussian random variables with variance and are independent of each other. For the time being, we assume perfect synchronization among three terminals. Note that they differ from our previous research in [17], because the impulse response of the cooperative channels, and , is modeled as highly dense channel model due to the fact that the relay can reduce transmission range and improve channel quality. In other words, multipath taps arrive in a very short delay spread. The two channels are assumed to have length and , respectively. For simplicity, we assume that they have same length , and the channel model of can be written as where and represent the complex path gain with and symbol spaced time delay of the th path, respectively. The training sequence vector is denoted as where is the number of training length, and the transmit power is , where is the symbol power. According to the property of AF cooperative system, one full transmission can be divided into two time slots as shown in Figure 3.
At the first time slot, complex baseband received signal at and is given by where and denote two complex circulant channel matrices with their first columns and respectively [20]; and are a realization of a complex additive Gaussian white noise vector with zero mean and covariance matrix , and is the identity matrix. Then the relay amplifies the received signal and retransmits the signal during the second time slot. The received signal vector at the destination is given by where is a circulant channel matrix with first column ; is a composite noise with zero mean and covariance matrix , where is a realization of a complex additive Gaussian white noise (AWGN) vector with zero mean and covariance matrix . Considering long-time averaging, the amplification factor is given by where is the transmit power of relay. Using (2), the effective input-output relation in the AF cooperative communication system can be summarized as According to the matrix theory [21], all circulant matrices can share the same eigenvectors [20]. That is to say, the same unitary matrix can work for all circulant matrices. Hence, the matrices and in (5) are decomposited as and , respectively, where denotes a diagonal matrix and is the unitary discrete Fourier transform (DFT) matrix with entries , . At the same time, is the decomposition of a circulant matrix which is constructed from a cascaded channel impulse response . Here, both and are diagonal matrices which are given by respectively, where and are given by respectively, where denotes direct link from source to destination at the first time slot and represents cascaded channel from source to destination via help of the relay at the second time slot. Based on the above analysis, (5) can be rewritten as where denotes -length received signal vector; denotes equivalent training matrix, and it can be written as with dimension; represents -length cooperative channel vector; denotes -length complex AWGN vector; and are partial DFT matrices taking the first and columns of , respectively. And the is a realization of a complex Gaussian random vector with zero mean and covariance matrix .
3. Sparse Channel Estimation
In this section, we discuss the sparse channel estimation for AF CCS. Firstly, we review briefly CS theory and restricted isometry property (RIP) of training signal matrix. Then, sparse channel estimation method [17] is introduced. Finally, we propose improved sparse channel estimators by using partial sparse constraint to take full advantage of prior information in AF CCS.
3.1. Review of the CS
In a typical complex sparse identification system, one can use known matrix to estimate an -length unknown sparse signal vector based on the observation linear system model: where is a complex observation signal vector, is a noise vector, and is sparse vector which means the number of dominant entries is no more than , that is, . The position of dominant entries is randomly distributed. In addition, according to CS assumption. The optimal sparse solution can be obtained uniquely by solving minimization problem: where is regularized parameter which trades off the mean square error (MSE) and sparsity. However, solving norm is NP hard and cannot be utilized in practical applications [22].
Fortunately, alternative suboptimal sparse recovery methods have been studied if the known measurement matrix satisfies RIP [23]. Let , be the submatrix extracting those columns of that are indexed by the elements of . Then the -restricted isometry constant (RIC) of is defined as the smallest parameter such that for all with and all vector . Assume that is an random measurement matrix that satisfies the RIP of order with RIC , that is, . Consider an arbitrary sparse vector in observation model , where , by solving minimization problem, and suboptimal sparse solution is obtained by where and is a parameter which is decided by the noise level and RIC of . Hence, the estimator satisfies sparse recovery performance with where is a parameter which is also decided by noise level and RIC of . Let us recall the channel estimation problem for AF cooperative systems in (8); if the equivalent training matrix satisfies RIP, then accurate sparse channel estimation can be achieved. In the next, we will present improved sparse channel estimation methods by using LASSO algorithm [19].
3.2. Sparse Channel Estimation
Channel estimation is done on sparse channel by sending the training symbols. Conventional sparse channel estimation method using LASSO algorithm (SEL) has been proposed for deriving sparse impulse response for AF CCS [17]. According to the system model in (8), the global sparse channel estimator can be achieved by where is an identity matrix and is a regularization parameter which controls the tradeoff between square error and sparse constrained . However, the proposed method can only solve global sparse solution well while neglecting the inherent partial sparse structure. In the next, we propose a method to fully exploit the prior information in AF CCS.
3.3. Partial Sparse Channel Estimation
From signal processing perspective, extra prior information of partial sparse can be further utilized. In this situation, partial sparse channel estimation by using LASSO (PEL) could be achieved by where is a diagonal weighted matrix and is a regularization parameter which controls the tradeoff between square error and local partial sparse constrained .
Based on the partial sparse constraint on cooperative channel impulse response, we propose an improved PEL (IEL) estimator. On the one hand, the local sparse constrain can improve estimation performance. On the other hand, the global sparse constraint can mitigate noise interference in the low SNR regime. The IEL estimator can be obtained by where the regularization parameters and are given by the (15) and (18), respectively. In the following, we will give representative simulation results to confirm the effectiveness of the improved sparse channel estimation methods.
4. Numerical Simulations
In this section, we will compare the performance of the proposed estimators, that is, PEL and IEL, with SEL estimator and LS estimator. To achieve average estimation performance, 1000 independent Monte-Carlo runs are adopted. The length of training sequence is . The length of direct link is with the number of dominant channel taps . The two cooperative links and have fixed length of with the number of dominant channel taps . All of the nonzero channel taps are generated following Rayleigh distribution and set to . Transmit power and AF relay power are fixed as , where is the symbol power. The received SNR is defined as .
Channel estimator is evaluated by normalized mean square error (NMSE) which is defined by where and denote cooperative channel vector and its estimator, respectively. At first, we compare their estimation performance with different number of dominant channel taps and . As shown in Figure 4, when and , the direct link is sparse channel impulse response, and the cascaded link may not be sparse since the linear convolution between and . It can be observed from Figure 4 that the two proposed channel estimators are better than both SEL estimator and LS-based linear channel estimator. It is worth nothing that the IEL estimator has a better performance than PEL one, since the IEL takes not only advantage of partial sparse prior information but also utilizes global sparse constraint to mitigate noise interference. The same performance advantage can also be seen in other scenarios with different number of dominant channel taps as shown in Figures 5 and 6. When and , we can also find that IEL estimator has a better performance than PEL under low SNR (less than dB). On the other hand, if the direct link is highly sparse, for example, , while the cooperative links and are highly dense, for example, , the two proposed channel estimators have a more significant performance advantage over traditional methods. In addition, IEL estimator is worse than PEL estimator when the SNR is higher than dB. According to these results, we can conclude the following: if direct and cooperative links are highly sparse channel, then IEL can achieve obvious better estimation performance than PEL; if the direct link is highly sparse channel while cooperative link is highly dense channel, the estimation performance of the two proposed methods is very close. It is worth mentioning that estimation performance of IEL is better than PEL due to the fact that IEL utilizes global sparse constraint to mitigate noise interference and partial sparse constraint to take advantage of channel sparsity. However, the computational complexity of IEL is higher than PEL. That is to say, IEL uses higher computational complexity than PEL to obtain performance advantage. Hence, to use IEL or PEL will be decided by the requirement of practical communication systems. In addition, when the number of dominant channel taps is very small, for example, and , the estimation performance is close to the SEL estimator. However, if the direct link is sparse, for example, and cooperative link is highly dense, for example, , then the two proposed channel estimators are close to CRLB. According to the previous analysis, we can find that the proposed methods are generalized from both LS-based linear estimation method and SEL, since they are either based on dense or sparse channel assumption. Hence, our proposed methods can work well in different channel environments.
5. Conclusion
Accurate CSI is indispensable for coherent detection in AF CCS. Traditional channel estimation methods are based on assumptions of either dense channel model or sparse channel model in AF CCS. In this paper, the two kinds of channel models have been generalized as a partial sparse channel model. By means of compressive sensing and partial sparse constraint, we have proposed improved sparse channel estimation methods to fully exploit channel prior information. Numerical simulations have confirmed the performance superiority of the proposed method over the conventional global sparse channel estimation method and traditional linear LS method. The proposed method can also be extended to other cooperative communication systems such as MIMO AF CCS.
Acknowledgment
The authors would like to thank Professor Guan Yong Liang for helpful discussions and appreciate the constructive comments of the anonymous reviewers.