Massive-MIMO Sparse Uplink Channel Estimation Using Implicit Training and Compressed Sensing

Massive multiple-input multiple-output (massive-MIMO) is foreseen as a potential technology for future 5G cellular communication networks due to its substantial benefits in terms of increased spectral and energy efficiency. These advantages of massive-MIMO are a consequence of equipping the base station (BS) with quite a large number of antenna elements, thus resulting in an aggressive spatial multiplexing. In order to effectively reap the benefits of massive-MIMO, an adequate estimate of the channel impulse response (CIR) between each transmit–receive link is of utmost importance. It has been established in the literature that certain specific multipath propagation environments lead to a sparse structured CIR in spatial and/or delay domains. In this paper, implicit training and compressed sensing based CIR estimation techniques are proposed for the case of massive-MIMO sparse uplink channels. In the proposed superimposed training (SiT) based techniques, a periodic and low power training sequence is superimposed (arithmetically added) over the information sequence, thus avoiding any dedicated time/frequency slots for the training sequence. For the estimation of such massive-MIMO sparse uplink channels, two greedy pursuits based compressed sensing approaches are proposed, viz: SiT based stage-wise orthogonal matching pursuit (SiT-StOMP) and gradient pursuit (SiT-GP). In order to demonstrate the validity of proposed techniques, a performance comparison in terms of normalized mean square error (NCMSE) and bit error rate (BER) is performed with a notable SiT based least squares (SiT-LS) channel estimation technique. The effect of channels’ sparsity, training-to-information power ratio (TIR) and signal-to-noise ratio (SNR) on BER and NCMSE performance of proposed schemes is thoroughly studied. For a simulation scenario of: 4× 64 massive-MIMO with a channel sparsity level of 80% and signal-to-noise ratio (SNR) of 10 dB, a performance gain of 18 dB and 13 dB in terms of NCMSE over SiT-LS is observed for the proposed SiT-StOMP and SiT-GP techniques, respectively. Moreover, a performance gain of about 3 dB and 2.5 dB in SNR is achieved by the proposed SiT-StOMP and SiT-GP, respectively, for a BER of 10−2, as compared to SiT-LS. This performance gain NCME and BER is observed to further increase with an increase in channels’ sparsity.


Introduction
Today's modern cellular communication networks are witnessing an ever increasing demand for higher data rates and link reliability. A major hindrance in meeting these demands is scarcity of the electromagnetic spectrum. One of the key solutions to effectively deal with these challenges is to spatially reuse the available spectrum. This is achieved by employing multiple antennas at the access point and/or user equipment (UE), thus enabling the propagation channel to provide more degrees of freedom. Such multiple-input multiple-output (MIMO) communication systems are are computationally complex and have slow convergence rate. Semi-blind methods use some training symbols along with the statistics of received signal to estimate the CIR. Semi-blind channel estimation techniques share the benefits and trade-offs of the blind and training based estimation techniques. Recently, SiT based techniques have gained significant attention for the purpose of channel estimation. In an SiT based technique, a low power and periodic training sequence is superimposed over the information sequence and transmitted. Consequently, SiT based techniques avoid the allocation of any dedicated time/frequency slots for the training sequence.
Over the past few years, compressed sensing (CS) has emerged as a new paradigm for the recovery of sparse signals. The authors in [27,28] have established the fact that a finite-dimensional sparse signal can be exactly reconstructed from fewer, linear and nonadaptive measurements. The CS approach has been established as an efficient solution to estimate sparse multipath channels-see e.g., [14,29]. Computing the sparse solution requires solving a 0 -minimization problem, which is computationally non-deterministic polynomial-time (NP) hard. An alternative approach is to relax the 0 -minimization to 1 -minimization and solve the problem for a sparse solution. This approach is termed as Basis Pursuits (BP) [30]. A BP based approach achieves more accurate solutions but requires higher computational complexity. Another approach used to find the sparse solution of an underdetermined system of linear equations is that of greedy algorithms, which iteratively approximate the sparse signal by suitably choosing the columns from the sparse measurement matrix. In literature, a wide variety of the greedy algorithms have been proposed to solve the CS problem such as matching pursuit (MP) [31], orthogonal matching pursuit (OMP) [32], regularized orthogonal matching pursuit (ROMP) [33], stagewise orthogonal matching pursuit (StOMP) [34], orthogonal complementary matching pursuit (OCMP) [35], compressive sampling matching pursuit (CoSaMP) [36], and gradient pursuit (GP) [37].
Several massive-MIMO channel estimation techniques have been proposed in the literature-see e.g., [23,38] and references therein. In order to estimate CIR in massive-MIMO systems, linear minimum mean-squared error (LMMSE) estimation techniques are usually used. Several CS based techniques have also been proposed to estimate the CIR in massive-MIMO (see e.g., [39,40]). In [41], authors have proposed a low complexity polynomial channel estimation based on Bayesian channel estimators to estimate CIR in massive-MIMO. An estimation technique for an uplink channel of massive-MIMO has been proposed in [42] that exploits the joint sparsity of the channels in a massive-MIMO system. In [43], authors have devised a Gaussian-mixture Bayesian learning based channel estimation technique for massive-MIMO. In [44], authors have devised a channel estimation technique that requires a small number of pilot sequences and exploits the sparsity and common support property of massive-MIMO communication channels. In [45], authors have exploited CS along with random linear network coding (RLNC) to devise an energy efficient scheme for vital signal telemonitoring in wireless body area networks (WBAN).
By exploiting the sparsity of wireless multipath channels, SiT sequence based compressive channel sensing methods have been studied in various contexts such as single-input single-output (SISO) systems [14,46], sparse MIMO channels [47,48], and underwater acoustic channels [49]. In [46], a genetic algorithm (GA) based channel estimation method is proposed using an SiT sequence for SISO systems. In [14], a Dantzig selector (DS) algorithm based method is proposed for estimation of SISO sparse multipath channels using a known SiT sequence. This study is further extended in [47,48] for the case of multiuser MIMO systems, where SiT based DS and MP algorithms are proposed.
To the best of the authors' knowledge, no study has considered the use of CS based greedy algorithms in conjunction with an SiT sequence for estimation of massive-MIMO sparse uplink channels. The prime motivation for using a superimposed training based channel estimation approach is its improved spectral efficiency as compared to the conventional training based approaches. In this regard, the main contributions of this paper is to devise an SiT sequence based estimation of massive-MIMO sparse uplink channels by exploiting StOMP and GP as CS greedy algorithms. For large-scale sparse reconstruction problems, it has been established in [34] that the StOMP algorithm achieves better performance as compared to MP and OMP. Similarly, the GP algorithm enhances the performance (in terms of faster computations) of the OMP algorithm by adopting a directional gradient pursuit based approach [37]. As the CIR estimation of massive-MIMO sparse uplink channels involves a large number of channel coefficients, the choice of StOMP and GP as sparse reconstruction algorithms is more favorable for this purpose.
The rest of the paper is structured as follows. In Section 2, the system model under consideration for multiuser massive-MIMO uplink communications is presented. Section 3 presents first-order statistics of the received signals along with the proposed channel estimation techniques. In Section 3.1, an SiT least squares (SiT-LS) based channel estimation technique available in literature [26] is presented. Section 3.2 presents the proposed SiT-StOMP and SiT-GP channel estimation techniques for the case of sparse uplink channels in massive-MIMO. Section 4 discusses a minimum mean square error (MMSE) equalizer for a massive-MIMO uplink communication scenario. Section 5 presents the simulation results and a performance analysis for the proposed channel estimation techniques. In Section 6, conclusions are presented.
Notations: Matrices are denoted by boldface uppercase letters, e.g., X, while boldface lowercase letters are used to represent vectors, e.g., x. Small case letters, e.g., x are used to denote scalar quantities. Transpose and Hermitian transpose are represented by the superscripts T and H, respectively. Γ i represents the set of indices of the elements selected up to and including iteration i. The matrix X Γ i represents a sub-matrix of X consisting of only those columns of X whose indices belong to the set Γ i . The same convention is also followed for the case of vectors.

Massive-MIMO System Model for Uplink Communications
The uplink communications system model for massive-MIMO is assumed to consist of N t single antenna UEs communicating with N r antennas at the BS, as depicted in Figure 1. The information sequence to be transmitted by n t -th UE is represented by b n t = [b n t (0), b n t (1), · · · , b n t (M − 1)] T , such that b n t is zero-mean, statistically independent of other UEs, and with E{| b n t (k) | 2 } = 1. Let c n t = [c n t (0), c n t (1), · · · , c n t (M − 1)] T represent a low power and periodic training sequence for n t -th UE. The training sequence c n t is periodic with period P such that c n t (k) = c n t (k + aP), for k and a being any integers, and is assumed to also be known at the BS side. After superimposing the training sequence c n t over the information sequence, b n t , the resultant transmit signal for n t -th UE is given as The sequence x n t = [x n t (0), x n t (1), · · · , x n t (M − 1)] T is then transmitted by the n t -th UE over the massive-MIMO sparse uplink channel. Between n t -th UE and n r -th receive antenna at the BS, a frequency selective and time invariant channel is assumed with CIR denoted by a sparse vector h n r n t = [h 0 n r n t , h 1 n r n t , · · · , h L−1 n r n t ] T . The number of resolvable mutipaths is represented by L. Moreover, the CIR vector h n r n t is assumed to be Q sparse i.e., {Q = h nm 0 } L, and has support p = [p 0 ,p 1 , · · · ,p Q−1 ] such that h n r n t = = 0 ; ∈p; = 0 ; otherwise.
The signal received at time instant k by the n r -th antenna element of receiver array at the BS is given below, where z n r (k) denotes kth sample of zero-mean, additive white Gaussian noise (AWGN) with variance σ 2 z . At time instant k, the overall combined received signal by all antennas at the BS is denoted by y(k) = [y 1 (k), y 2 (k), · · · , y N t (k)] T and given by, , · · · , z N r (k)] T and the N r × N t channel matrix H for th tap delay is given by, . . .
The channel estimator (CE) block, shown in Figure 1, is implemented by using the SiT and first-order statistics based least squares (SiT-LS) technique presented in [50] and the proposed SiT-StOMP and SiT-GP techniques. Once the CIR is estimated, the effect of superimposed training sequence is removed from the information sequence at the receiver side by the training effect remover (TER) block, as shown in Figure 1. After eliminating the superimposed training sequence effect, the output of TER is fed as an input to the LMMSE equalizer.

Massive-MIMO Sparse Uplink Channel Estimation Using a First-Order Statistics Based Approach
For the uplink scenario of a massive-MIMO communication system, if each mobile user is assigned with a specific pilot sequence that is superimposed on the information sequence, a first order statistics of the received signal can be used to estimate the CIR as outlined in [50]. In this section, an SiT based channel estimation technique of [26] is extended for the estimation of massive-MIMO sparse uplink channels. For n t -th transmitter, the training sequence c n t (k) is periodic with period P =PN t , whereP is a positive integer. The training sequence c n t (k) is given as, where j = √ −1 and The pilot sequence c n t (k) is selected such that onlyP coefficients among P are non-zero. Hence, the training sequence c n t (k) for n t -th UE is given by, where α i,n t = 2π(iN t + n t − 1)/P, and c i,n t are suitably chosen coefficients for 1 ≤ n t ≤ N t and 0 ≤ i ≤P − 1. The training sequence c n t (k) can be designed by first choosing a periodic base sequencē c o (k) with a period ofP [26] such that, The training sequencec 1 (n t ) with period P is then generated by replicatingc o (k) for N t times. Therefore, for n t -th UE, the training sequence is obtained as [26] c n t (k) = σ cn tc1 (k) e j(2π/P)(n t −1)k for n t = 1, 2, · · · , N t .
Taking the expected value of the received signal y n r (k) at n r -th receive antenna gives For n 1 = n 2 , we have α i 1 ,n 1 = α i 2 ,n 2 for any {i 1 , i 2 } ∈ 0, 1, ...,P − 1. Let d n r n t = [d n r n t ,0 , d n r n t ,1 , · · · , d n r n t ,(P−1) ] T , where d n r n t ,i is given by The mean square consistent estimated n r n t = [d n r n t ,0 ,d n r n t ,1 , · · · ,d n r n t ,(P−1) ] T of d n r n t is computed as in [26] and is given byd where M is the total number of received symbols, as M → ∞,d n r n t ,i → d n r n t ,i . The vector form of (13) is given below,d n r n t = C n t h n r n t , where the matrix C n t is computed as where V n t is the Vandermonde matrix given by For distinct α i,n t , the rank of matrix V n t is L ifP ≥ L + 1 ([51], p. 257). Since all α i,n t are distinct and c i,n t = 0 ∀i, rank(C n t ) = rank(diag {c 0,n t , c 1,n t , · · · , c (P−1),n t }V n t ) = L, h n r n t can be uniquely determined from Equation (14).

SiT Based Least Squares (SiT-LS) Channel Estimation Approach
The least squares estimate of CIR between n t -th transmitter and n r -th receiver can be obtained from the linear model in Equation (14), as proposed in [50], and is given beloŵ The above estimate can also be obtained aŝ h SiT−LS n r n t = (C H n t C n t ) −1 C H n td n r n t .
To obtain the channel estimate for non-zero mean noise, setP ≥ L + 1, omit the first row from C n t andd n r n t ,0 fromd n r n t .

Proposed SiT Based Massive-MIMO Sparse Uplink Channel Estimation Techniques
The estimation error betweend n r n t ,i and d n r n t ,i can be computed by using Equation (3) in Equation (13), and is given byd n r n t ,i = d n r n t ,i + ε n r n t ,i , where ε n r n t ,i denotes the estimation error of d n r n t ,i . This estimation error constitutes contributions from interference due to superimposed information sequence of all UEs (b n r n t ,i ), additive noise (z n r n t ,i ), and interference due to training sequences of cross channels (c n rñt ,i ). The estimation error is thus given by, ε n r n t ,i =c n rñt ,i +b n r n t ,i +z n r n t ,i , wherẽ b n r n t ,i = 1 M Taking into account this inherent error ε n r n t ,i results in more adequate estimate of the CIR. Moreover, the first-order statistics based technique discussed previously does not take into consideration the sparse nature of uplink massive-MIMO channels. This section, thus, presents two extensions of the first-order statistics based channel estimation techniques for sparse uplink massive-MIMO channels by using the model in Equation (14).

SiT Based Stage-Wise Orthogonal Matching Pursuit (SiT-StOMP)
In [34], StOMP was proposed with an objective of performance enhancement in the reconstruction of sparse signals for large-scale CS problems while keeping the computational cost low. This performance enhancement in StOMP is achieved by allowing for selecting multiple columns per iteration as opposed to single column selection based strategy of MP and OMP algorithms. It has been established in [34] that for the reconstruction of large-scale sparse signals, StOMP performs faster than MP and OMP.
At each iteration of the StOMP algorithm, several columns from the measurement matrix are added to the active set as compared to the single column selection of the OMP algorithm. In StOMP, the selection of columns takes place according to a certain pre-determined threshold value. Only those columns are selected whose absolute correlations with the current residual exceed the threshold value. After selection of columns, StOMP then solves for a least squares problem to update the residual vector. Thus, StOMP converges faster than OMP since it requires less number of iterations to reconstruct the sparse solution.
Due to faster convergence and better performance for large-scale systems, StOMP has been considered as a viable solution for the estimation of massive-MIMO sparse uplink channels. Therefore, in order to reconstruct the sparse channel vector h n r n t from the model presented in Equation (14), we incorporate the StOMP algorithm. The proposed SiT-StOMP algorithm for the estimation of channel vector between n r -th receiver and n t -th transmitter is outlined as below: Input: Matrix C n t , vectord n r n t , and threshold i . Output: Channel estimate vectorĥ StOMP n r n t .
1. Initialize residual r 0 =d n r n t , index set Φ 0 = ∅, and iteration counter i = 1. 2. Create a set Ω i consisting of the indices of all elements in the vector, ψ i = C H n t r i−1 , which are above the threshold i 3. Update the index set by Φ i = Φ i−1 ∪ Ω i and residual by h n r n t = arg min h nr n t ∈R Φ i d n r n t − C n th n r n t 2 2 , r i =d n r n t − C n tȟ n r n t . 4. Check stopping criteria; if it is not met then update index i = i + 1, and go to step 2; if stopping criteria is met, set the final output vector asĥ StOMP n r n t =ȟ n r n t .
The stopping criteria in this case is a fixed number of maximum iterations, Λ. In order to avoid false alarms and missed detection, as proposed in [34], the threshold is set as i = t i r i 2 / √ P, where 2 ≤ t i ≤ 3. When the channel is sufficiently sparse, after the algorithm exits, Φ i is expected to have no more than P entries and all the non-zeros inĥ StOMP n r n t are selected in Φ i .

SiT Based Gradient Pursuit (SiT-GP)
The GP [37] is a greedy algorithm for the reconstruction of sparse signals. It utilizes the steepest descent methodology to compute the step-size for each iteration, i.e., the sparse solution vector is updated at each iteration with a directional update computed on the basis of gradient or conjugate gradient. The only additional computational cost compared to MP is that of the evaluation of the step size. The SiT-GP algorithm is given below: 1. Initialize the residual vector r 0 =d n r n t , the estimate of the channel coefficients vectorĥ 0 n r n t = 0, and Γ 0 = ∅; 2. for i = 1; i = i + 1 until stopping criteria is met, do The adoption of directional gradient pursuit in the GP algorithm results in fast approximations to OMP but with the same computational complexity as that of MP algorithm.

Minimum Mean Square Error (MMSE) Based Equalizer
Once an adequate estimate of the channel vector h n r n t has been obtained, the effect of the training sequence must be removed from the received signal before equalization. This is achieved by using TER block, as shown in Figure 1. Since the training sequence for each UE is also known at the BS, the effect of this superimposed training sequence can be removed as below n r n t c n t (k − ), (23) whereĥ n r n t is the estimate of th tap CIR from n t -th UE to n r -th receive antenna. The CIR estimateĥ n r n t may be taken from any of the estimation techniques discussed previously, i.e.,ĥ SiT−LS n r n t ,ĥ SiT−StOMP n r n t orĥ SiT−GP n r n t . The optimal weights,w n r , for the equalizer at n r -th receive antenna, can be computed as in [52], and are given below w n r n t = ĤĤ * + 2σ 2 where n r is the receiver index such that 1 ≤ n r ≤ N r , L e is the length of equalizer, τ d represents the symbol mappers' decision delay at equalizer's output, I denotes the (N t × L e ) × (N t × L e ) identity matrix andĤ| i is the ith column ofĤ. The estimate of the noise variance σ 2 n r is obtained as in [26]. The L e × (L e + L − 1) convolutional matrixĤ is given bŷ . . . 0ĥ 0 n r n tĥ 1 n r n t · · ·ĥ L−1 The overall convolutional matrixĤ is given beloŵ The information symbols transmitted from the n t -th transmitter can thus be estimated as below, followed by a symbol mapper, as shown in Figure 1: The estimated symbol vector,b n , is then fed as an input to a symbol mapper, as shown in Figure 1. The symbol mapper performs mapping of the estimated symbols as per used modulation scheme at the UEs. The output of the symbol mapper is the decoded symbol vectorb n .

Results and Discussion
This section presents the computer based simulation results of the proposed techniques and analysis of the obtained results. The normalized channel mean square error (NCMSE) and bit error rate (BER) have been used as the performance metrics for this purpose. The NCMSE for the proposed estimation techniques is defined as A massive-MIMO system consisting of 64-antenna BS serving four UEs i.e., N r = 64 and N t = 4, is considered for the purpose of simulations. The underlying massive-MIMO sparse uplink channels are assumed to be time-invariant and frequency-selective. The channel vectors h n r n t , having a fixed sparsity level Q/L, are independently generated for each Monte Carlo run. The non-zero channel coefficients of h n r n t are drawn from a zero-mean Gaussian distribution with variance 1/(N r (L + 1)). Moreover, a fixed channel length of L = 14 is used for all of the channels. The positions of non-zero channel taps of all the channels from a certain UE to all of the receive antennas at the BS is taken as the same because of the small separation distance of the antennas as compared to the line-of-sight (LoS) path length. The periodic and low power training sequence for each UE is generated by using the m− sequence based approach, as presented in [26]. For this purpose, a base sequence {−1, −1, −1, 1, 1, 1, 1, −1, 1, −1, 1, 1, −1, −1, 1} with periodP = 15 is used for all of the simulation results. AWGN is independently generated at each receiver for a specific signal-to-noise (SNR) ratio. The SNR at n r -th receiver is defined as the ratio of power of received signal σ 2 y nr to noise power σ 2 n r , i.e., SNR n r = σ 2 y nr /σ 2 n r . The zero mean binary phase-shift keying (BPSK) modulated information sequences (b n t ∈ {1, −1}) are generated mutually independent for each transmitter.
Performance comparison of the proposed SiT-StOMP and SiT-GP with that of the SiT-LS is presented in Figure 2 for NCMSE and BER against SNR, respectively. For this purpose, length of the information sequence is set as M = 1500 bits, sparsity level of the channel vectors as Q/L = 3/14, and training-to-information ratio (TIR) as σ 2 c n t /σ 2 b n t = 0.2, respectively. From Figure 2a, it can be seen that the proposed schemes perform much better as compared to SiT-LS. In terms of NCMSE at SNR = 10 dB, the proposed techniques SiT-StOMP and SiT-GP give a performance gain of 18 dB and 13 dB, respectively, over the SiT-LS. The BER based performance comparison is shown in Figure 2b. For a BER of 10 −2 , a performance gain of about 3 dB and 2.5 dB in SNR is achieved by proposed SiT-StOMP and SiT-GP, respectively, over SiT-LS. In order to demonstrate the performance of proposed schemes for different variants of channel's sparsity parameter, i.e., Q/L, the NCMSE and BER are plotted in Figures 3 and 4, respectively, for several values of Q/L (i.e., Q/L = 3/14, 5/14 and 7/14) by keeping the value of TIR equal to 0.2.
It can be observed from Figures 3 and 4 that NCMSE and BER decrease with increasing channel sparsity and vice versa. However, SiT-StOMP performs better than SiT-GP, in terms of NCMSE, even if the channel is less sparse. To determine the effect of variations in training power sequence, NCMSE is plotted in Figure 5 for the proposed SiT-StOMP and SiT-GP for several values of TIR i.e., σ 2 c n t /σ 2 b n t = 0.2, 0.5 and 1. For this purpose, the channel sparsity parameter Q/L is set equal to 3/14. It can be observed from Figure 5 that, for a fixed information sequence power, σ 2 b n t , increasing the training sequence power σ 2 c n t results in an improved NCMSE performance for both SiT-StOMP and SiT-GP. However, increasing the TIR leads to a poor BER performance, as this increased training power could have been utilized for the information sequence. Therefore, a suitable value of TIR must be chosen in order to relieve the NCMSE and BER trade-off. A detailed performance comparison (in terms of NCMSE versus channel sparsity and TIR) of SiT-StOMP, SiT-GP, and SiT-LS algorithms is presented in Figure 6a,b, respectively. In Figure 6, it can be observed that, for all levels of TIR and channel sparsity, both SiT-StOMP and SiT-GP outperform SiT-LS. Moreover, it is evident from Figure 6a that for a TIR of 0.3, SiT-StOMP outperforms SiT-GP for high levels of channel sparsity; however, for a less sparse channel, a converse behaviour is observed. This performance gain (for highly sparse channels) of SiT-StOMP compared to SiT-GP is also influenced by the TIR level (which is also depicted in Figure 5). Therefore, the effect of TIR levels on NCMSE performance of the proposed algorithms is demonstrated in Figure 6b. It can be observed that SiT-StOMP performs equally well both for lower as well as higher values of TIR in comparison with SiT-GP and SiT-LS (plotted for a highly sparse channel). However, there is a performance trade-off (as discussed earlier) between the accuracy of channel estimates and retrieved information sequence for different values of TIR. Therefore, it is desirable to keep TIR at a certain low value (e.g., TIR = 0.2) that promises an adequate BER performance. With such setting of TIR for an environment exhibiting highly sparse CIR, it is realized that SiT-StOMP is a favorable choice for obtaining channel estimates.

Conclusions
Two channel estimation techniques based on implicit training and compressed sensing have been proposed for massive-MIMO sparse uplink channels. A comprehensive analysis based on the results obtained from computer based simulations of these techniques has been presented. Performance of the proposed schemes is evaluated on the basis of NCMSE and BER as the performance criteria. In order to prove the validity of simulation results, NCMSE and BER based performance comparison of the proposed schemes with that of a notable SiT-LS scheme has been presented. Moreover, effect of variations in channel sparsity parameter and training-to-information ratio has also been presented for the proposed techniques. It has been established that the proposed SiT-StOMP and SiT-GP techniques outperform the first-order statistics based SiT-LS in terms of NCMSE and BER for the case of sparse multipath channels. It has been shown that the proposed SiT-StOMP and SiT-GP can provide a performance gain of 18 dB and 13 dB, respectively, in terms of NCMSE at an SNR of 12 dB and channel sparsity of 80%, over SiT-LS. Similarly, for a BER of 10 −2 , a gain of about 3 dB and 2.5 dB in SNR is achieved by SiT-StOMP and SiT-GP, respectively, over SiT-LS, for a channel sparsity level of 80%. Furthermore, in order to compare the performance of proposed SiT-StOMP and SiT-GP techniques, the NCMSE has been plotted against variations in channel sparsity and TIR. It has been demonstrated that, for sparser channels, SiT-StOMP achieves better performance in terms of NCMSE as compared to SiT-GP. However, as the channel becomes less sparse, the NCMSE for SiT-StOMP increases in comparison to that of SiT-GP. Furthermore, both of the proposed techniques perform better than SiT-LS even if the channel is less sparse. Moreover, it has been shown that SiT-StOMP performs better as compared to SiT-GP for lower as well as higher values of TIR.