Research Article  Open Access
Jingke Dai, "BitErrorRate Analysis of Raptor Codes over Rician Fading Channels", Journal of Electrical and Computer Engineering, vol. 2020, Article ID 2685075, 8 pages, 2020. https://doi.org/10.1155/2020/2685075
BitErrorRate Analysis of Raptor Codes over Rician Fading Channels
Abstract
The biterrorrate expressions of nonsystematic Raptor (NR) codes and systematic Raptor (SR) codes over Rician fading channels are first derived using a Gaussian approximation (GA) approach. These BER expressions provide a significant reduction in computational complexity for analyzing system performance when compared with simulation and discretized density evolution (DDE). As shown by the asymptotic analysis, the NR codes originally designed for binary erasure channels still have good performance on Rician fading channels but SR codes do not. Therefore, the degree distributions of SR codes are specifically optimized on Rician channels which are superior to the existing SR codes and comparable to NR codes.
1. Introduction
Rateless codes are able to adapt the bit rate without channel state information at the transmitter and then are suitable for variant channels. Luby Transform (LT) codes [1] are the first class of practical rateless codes which are proposed for binary erasure channels (BEC). As an extension of LT codes, the Raptor codes with lower decoding complexity are constructed by serially concatenating inner LT codes with an outer highrate code [2]. If the original symbols are first transmitted followed by LTcoded symbols, the corresponding rateless codes could be called systematic LT (SLT) codes [3] or systematic Raptor (SR) codes [4].
LT codes and Raptor codes have been generalized to noisy channels, such as binary symmetric channels (BSC) [5], additive white Gaussian noise (AWGN) channels [5–8], multipleaccess channels [9], and fading channels [10], which have been shown to achieve good decoding performance. The asymptotic analyses of them are normally obtained by discretized density evolution (DDE) [4, 8] and some simplified onedimensional (1D) analysis methods, e.g., Gaussian approximation with message mean (GAmean) [3, 6] and extrinsic information transfer (EXIT) chart [7]. By combining those asymptotic analyses with differential evolution algorithm or linear programming (LP) method, the performance of LT codes and Raptor codes can be further improved [6–8]. As is well known, the DDE provides an accurate result but has high complexity, while those 1D methods sacrifice reasonable accuracy for much faster calculation. In addition, the lower bounds of biterrorrate (BER) for LT codes and lowdensity generatormatrix (LDGM) codes are analyzed based on GAmean [3], EXIT [7], and DDE [11], respectively, which match the corresponding asymptotic performance very well at high signaltonoise ratio (SNR) but significant discrepancy occurs at low SNR.
Meanwhile, most of the above literature studies with asymptotic analysis concentrate on BEC or AWGN channels, and the performance of LT/Raptor codes on fading channels is usually obtained by computer simulations [10, 12–15]. Among them, the study in [10] reveals that the rateless coding scheme has advantages in efficiency and reliability over fixedrate codes when channels information is not available and provides the simulated performance of Raptor codes on Rayleigh fading channels. The simulation results in [12] show that the Raptor codes can closely approach the capacity limits of both Rayleigh fading channels and Rician fading channels.
Recently, based on a modified GA method, the biterrorrate of systematic Raptor codes on Rayleigh/Nakagami fading channels has been derived in [16], providing a way of asymptotic analysis for Raptor codes on fading channels. We extend the work of [16] to nonsystematic Raptor codes and Rician channels. First, the BER expressions of both NR codes and SR codes are derived on Rician channels. Comparisons among simulation, DDE analysis, and the proposed analytical BER demonstrate that they are in good agreement and our analytical expression has a significant reduction in computation complexity compared to DDE analysis. In addition, the asymptotic results show that the decoding thresholds of NR codes are very close to the Shannonlimits though the employed degree distribution is originally designed for BEC. Second, the SR codes with the distribution designed for NR codes are prone to poor performance; hence, the proposed BER expression is combined with LP method to search good degree distributions for SR codes. The simulation results indicate that the optimized SR codes have BER performance as good as NR codes and outperform existing SR codes which are specifically designed on AWGN channels.
2. Encoding and Decoding
The coded BPSK symbols from transmitter are mapped as and then the received symbols are given bywhere is a Gaussian variable with zeromean and variance and the fading factor is Rician distributed with the probability density function (PDF) as [12]where is called a Rician factor and is the zerothorder modified Bessel function of the first kind. The channel side information is assumed as known and then the loglikelihood ratio (LLR) provided for decoder is calculated as
Raptor encoding is a twostage process. The original information sequence of length is precoded using a highrate (LDPC) code to produce the intermediate sequence of length . This intermediate sequence is further LTencoded to produce the output of Raptor code. Here LT code is the inner code, which can be described by its output degree distribution . The integer is the maximum degree and denotes the probability that degree is chosen in the random XORencoding process [3]. If the NR codes are taken, the generation matrix of inner nonsystematic LT (NLT) codes is , where is a matrix of size constructed according to and is the actual number of LTcoded symbols. If the SR codes are taken, the generation matrix is constructed by , where is an identity matrix of size and is of size . The actual rates of inner codes, outer codes, and whole codes are denoted as , , and , respectively.
Raptor decoding is also divided into two steps for inner and outer decoding, respectively. Firstly, the classic sumproduct algorithm (SPA) is employed to decode inner code, which is based on matrix , the transpose of , for both NLT codes and SLT codes. Each column and each row of correspond to the variable node (VN) and check node (CN) of the decoding bipartite graph. The messages are updated between VNs and CNs as follows [3, 12]:where is the message passed from the nth check node to the mth variable node at the th decoding iteration and is the message passed in the opposite direction, The set of all neighbors of nth node except for the mth node is denoted by .
When the iteration achieves the maximum number , the inner decoder outputs the soft information (NLT) or (SLT) as the initial message provided to outer decoder. The outer (LDPC) code is also decoded by SPA, the detail of which is found in [12] and is not presented here.
3. BER Analysis of NR Codes
To analyze the BER of NR codes, some notations are predefined as follows. The edge degree distribution of CNs is calculated by their node degree distribution as given by and the average degree of CNs is . The average degree of VNs is denoted by and the variable node and edge degrees are both assumed to be Poisson distributed as and , respectively. They can also be truncated to obtain polynomials as and , where is large enough to ensure and .
If allzero codewords are transmitted [3], the initial LLR has the following conditional PDF as [16]
Combining with (2), we have the unconditional PDFand the corresponding initial BER before decoding is calculated as
First, the inner NLT decoding is considered. We denote the error probability at iteration as the average probability that the VNs messages are negative and as the check nodes with degree . According to (4), the probability of at iteration is derived aswhere the VNs messages are assumed to be i.i.d. Note that is assumed to be symmetric Gaussian distributed. Denoting as the mean of , we get the following equation:where function . By averaging over all the possible check node degrees, the mean of message by check nodes is
The VN of degree outputs a Gaussian message with meanas well as variance ; thus, the conditional BER is
By averaging over all the possible variable node degrees, we have the BER as given by
Note that is the function of .
At the iteration , the message provided from inner decoder to outer decoder can be assumed to be Gaussian distributed [8] with the mean and variance as
Since the initial message for outer decoder is Gaussian distributed, the EXIT chart technique can be employed to track the message in decoding iterations. The mutual information from check nodes at the th iteration is updated by [16]where and are the edge degree distributions of variable nodes and check nodes of outer LDPC codes, and are the corresponding maximum degrees, and the function is defined in [17].
At the maximum iteration, the BER of Raptor code is calculated bywhere denotes the variable node degree distribution of outer LDPC codes.
4. BER Analysis of SR Codes
From equations (4)–(7), we know the SLT decoding is very similar to NLT decoding and then the BER analysis of SR codes is also like that of NR codes. However, in (5), the message from VNs is the sum of LLR from CNs, whereas, in (7), the message from VNs is the sum of LLR from CNs and initial LLR from channel. Hence, the process of VNs for SLT codes is a little different from NLT codes.
Since the message from check nodes is assumed to be Gaussian distributed, the sum of them at the variable nodes with degree is also Gaussian distributed (see (14)) with the PDF as
Therefore, the PDF of message from VNs with degree is the convolution of (9) and (20) as given byand the corresponding average BER is calculated as
As the outer decoding is the same as that in NR codes, the BER of SR codes is also calculated by (19), only needing to substitute (22) for (16) in advance. Note that the convolution in (22) is easily obtained by Fast Fourier Transform (FFT).
We define the decoding overhead as . The outer codes are set as the 0.98rate (4,200)regular LDPC code and the 0.80rate (4,20)regular LDPC code, respectively. The degree distributions of inner codes for both NR codes and SR codes take the optimized result in [2], which is rewritten as
In Figures 1 and 2, the BER of NR codes and SR codes obtained by (19) are compared with the simulation results with information length , which are indicated as “Analytical” and “Simulation,” respectively. Furthermore, decoding threshold achieved by DDE and Shannonlimit are also depicted in those figures. Note that the Shannonlimit here is calculated by , where is the capacity of Rician channel with fading factor and noise variance . It can be seen that the analytical curves closely match the simulated results at high overhead region. Generally, at the BER of , the difference between the simulated overhead and the analytical threshold of errorfree decoding is around 0.03–0.15 for all the BER curves. For example, we observe from Figure 1(a) that the simulated overhead to BER of is around 2.41, whereas the analytical threshold is 2.31, which has a difference of 0.10. The performance of NR codes is much better than the SR codes, since the employed degree distribution of inner codes is originally designed for NR codes, though on binary erasure channel. In fact, the decoding threshold of NR codes is very close to the corresponding channel capacity. For example, in Figure 2(a), the DDEthreshold of NR codes is 1.73, the gap of which to Shannonlimit is only 0.07. The above results show that the NR codes have good robustness on different classes of channels and then the performance improvement space is limited (also see [8]), whereas the SR codes require careful design for specific channels.
(a)
(b)
(a)
(b)
From Figures 1(b) and 2(b), the simulation results of Raptor codes with outer rate of 0.80 are similar to the ones with outer rate of 0.98; only the overheads of former codes are smaller than the corresponding later codes. But it does not mean the codes with outer rate of 0.80 are better than the codes with outer rate of 0.98 because their performance should be evaluated by the whole rates at this time. For example, the overhead thresholds of NR codes in Figures 1(a) and 1(b) are 2.22 and 2.17, respectively, but the whole rates of them are around 0.441 and 0.369, respectively. In fact, we observe that the codes with outer rate of 0.80 have gaps between Shannonlimit and decoding threshold which are much larger than the gaps of codes with outer rate of 0.98, indicating that the outer rate should be set higher (e.g., ).
5. Optimization of SR Codes
The degree distribution of inner codes of SR codes is normally optimized by linear programming, which is presented as
In this design model, the objective is minimizing inner code’s overhead, which is calculated by in (24), and equations (26) and (27) ensure the validity of the distributions. Note that is predefined. in (25) is obtained by (22) and is chosen from points in the interval , where is an integer indicating the number of sampling points, is the initial error probability given by (10), and is the target error probability, usually set as the desired BER for successful decoding of outer LDPC codes [16].
Setting different and , the optimized degree distributions of inner SLT codes are shown in Table 1. The BER performance of SR codes with those optimized degree distributions is compared with that of the NR codes and SR codes with existing degree distributions. Besides the distribution , the optimized distribution on AWGN channel for rate1/2 SR codes in [4] is also considered and denoted as . The 0.98rate (4,200)regular LDPC code is fixed as the outer code and the iteration number in SPA is set to 100 for both inner codes and outer codes. All the codes have the same information length as .

The simulation results of those Raptor codes on Rician channels are shown in Figures 3 and 4. It is seen that the optimized SR codes indeed outperform the SR codes with existing distributions and the performance is very close to that of NR codes on all the considered channels, which indicates the effectiveness of our analysis and optimization method. For example, on the channel with and , the required overheads for NR codes and optimized SR codes achieving BER below are both 2.45. The maximum difference between NR codes and optimized SR codes occurs on the channel with and , and the overhead gap is only 0.03. Although the performance of optimized SR codes does not exceed that of NR codes, the good features of systematic codes should be considered, such as visible information bit, convenience of error checking, and smaller decoding matrix. The SR codes with have moderate performance on Rician channel when but have poor performance when because of being originally designed for AWGN channel, where the distribution of initial LLR is more like the one in former case. In addition, it should be noted that all the optimized SR codes are designed under specific channel conditions ( and ), and their performance would be degraded if the channel changes; that is, the SR codes are prone to be affected by the channel when compared with NR codes.
(a)
(b)
(a)
(b)
6. Conclusions
The BER expressions of nonsystematic Raptor codes and systematic Raptor codes on Rician fading channels have been derived using GA approach. The simulation results are in good agreement with the analytical results; thus the BER expressions provide a computationally efficient way of analyzing BER performance and predicting the overhead threshold with reasonable accuracy. With the existing degree distributions, NR codes still have Shannonlimit approaching performance on Rician channels but SR codes have poor performance. The proposed analytical expression is combined with linear programming to search good degree distributions for SR codes and the simulation results indicate that the optimized SR codes have BER performance as good as NR codes.
Data Availability
The numerical and simulation data used to support the findings of this study are available upon request to the corresponding author Jingke Dai (email: djk029@163.com).
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Grant no. 61501469) and China Postdoctoral Science Foundation (Grant no. 2015M572782).
References
 M. Luby, “LT codes,” in Proceedings of 43rd Symposium on Foundations of Computer Science, pp. 271–280, Vancouver, Canada, November 2002. View at: Google Scholar
 A. Shokrollahi, “Raptor codes,” IEEE Transactions on Information Theory, vol. 52, no. 6, pp. 2551–2567, 2006. View at: Publisher Site  Google Scholar
 S. Xu and D. Xu, “Optimization design and asymptotic analysis of systematic Luby transform codes over BIAWGN channels,” IEEE Transactions on Communications, vol. 64, no. 8, pp. 3160–3168, 2016. View at: Publisher Site  Google Scholar
 A. Kharel and L. Cao, “Asymptotic analysis and optimization design of physical layer systematic rateless codes,” in Proceedings of 15th IEEE Annual Consumer Communications & Networking Conference, Las Vegas, NV, USA, January 2018. View at: Google Scholar
 R. Palanki and J. S. Yedidia, “Rateless codes on noisy channels,” in Proceedings of the IEEE International Symposium on Information Theory (ISIT), Chicago, IL, USA, June 2004. View at: Google Scholar
 O. Etesami and A. Shokrollahi, “Raptor codes on binary memoryless symmetric channels,” IEEE Transactions on Information Theory, vol. 52, no. 5, pp. 2033–2051, 2006. View at: Publisher Site  Google Scholar
 I. Hussain, X. Xiao, and L. K. Rasmussen, “Error floor analysis of LT codes over the additive white Gaussian noise channel,” in Proceedings of 2011 IEEE Global Telecommunication Conference, Houston, TX, USA, December 2011. View at: Google Scholar
 A. Kharel and L. Cao, “Analysis and design of physical layer raptor codes,” IEEE Communications Letters, vol. 22, no. 3, pp. 450–453, 2018. View at: Publisher Site  Google Scholar
 D. Deng, D. Xu, and S. Xu, “Optimisation design of systematic LT codes over AWGN multiple access channel,” IET Communications, vol. 12, no. 11, pp. 1351–1358, 2018. View at: Publisher Site  Google Scholar
 J. Castura and Y. Yongyi Mao, “Rateless coding over fading channels,” IEEE Communications Letters, vol. 10, no. 1, pp. 46–48, 2006. View at: Publisher Site  Google Scholar
 A. Kharel and L. Cao, “Analysis and design of serially concatenated LDGM codes,” 2018, https://arxiv.org/abs/1801.08270v2. View at: Google Scholar
 B. Sivasubramanian and H. Leib, “Fixedrate raptor codes over Rician fading channels,” IEEE Transactions on Vehicular Technology, vol. 57, no. 6, pp. 3905–3911, 2008. View at: Publisher Site  Google Scholar
 P. Palanisamy and T. V. S. Sreedhar, “Performance analysis of Raptor codes in WiMax systems over fading,” in Proceedings of 2008 IEEE Region 10 Conference, Hyderabad, India, November 2008. View at: Google Scholar
 W. Huang, H. Li, and J. Dill, “Digital fountain codes system model and performance over AWGN and Rayleigh fading channels,” in Proceedings of International Conference on Computing Communications and Control, Orlando, FL, USA, April 2010. View at: Google Scholar
 M. K. Sharma, A. Verma, and J. Kaur, “Performance analysis of fountain coded OFDM system over Rician fading channel,” in Proceedings of 2014 International Conference on Advances in Communication and Computing Technologies, Mumbai, India, August 2014. View at: Google Scholar
 J. Dai, X. Chen, F. Zhang, and K. Kang, “Optimisation design of systematic fountain codes on fading channels,” IET Communications, vol. 13, no. 20, pp. 3369–3376, 2019. View at: Publisher Site  Google Scholar
 S. T. Brink, G. Kramer, and A. Ashikhmin, “Design of lowdensity paritycheck codes for modulation and detection,” IEEE Transactions on Communications, vol. 52, no. 4, pp. 670–678, 2004. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Jingke Dai. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.