Journal of Computer Networks and Communications
Volume 2016 (2016), Article ID 7918456, 7 pages
http://dx.doi.org/10.1155/2016/7918456

Research Article

ASEP of MIMO System with MMSE-OSIC Detection over Weibull-Gamma Fading Channel Subject to AWGGN

Keerti Tiwari,¹ Davinder S. Saini,² and Sunil V. Bhooshan¹

¹Department of Electronics and Communication Engineering, Jaypee University of Information Technology, Waknaghat, Solan 173234, India
²Department of Electronics and Communication Engineering, Chandigarh College of Engineering and Technology, Chandigarh 160019, India

Received 14 November 2015; Revised 28 January 2016; Accepted 15 February 2016

Academic Editor: Pantelis-Daniel Arapoglou

Copyright © 2016 Keerti Tiwari et al. 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.

Abstract.

Ordered successive interference cancellation (OSIC) is adopted with minimum mean square error (MMSE) detection to enhance the multiple-input multiple-output (MIMO) system performance. The optimum detection technique improves the error rate performance but increases system complexity. Therefore, MMSE-OSIC detection is used which reduces error rate compared to traditional MMSE with low complexity. The system performance is analyzed in composite fading environment that includes multipath and shadowing effects known as Weibull-Gamma (WG) fading. Along with the composite fading, a generalized noise that is additive white generalized Gaussian noise (AWGGN) is considered to show the impact of wireless scenario. This noise model includes various forms of noise as special cases such as impulsive, Gamma, Laplacian, Gaussian, and uniform. Consequently, generalized -function is used to model noise. The average symbol error probability (ASEP) of MIMO system is computed for 16-quadrature amplitude modulation (16-QAM) using MMSE-OSIC detection in WG fading perturbed by AWGGN. Analytical expressions are given in terms of Fox-H function (FHF). These expressions demonstrate the best fit to simulation results.

1. Introduction

In multiple-input multiple-output (MIMO) systems, the role of spatial multiplexing (SM) and spatial diversity is to provide high data rate and reliable communication, respectively. However, a tradeoff occurs between diversity and multiplexing for multiple access channels [1]. The demand to achieve high data rate is increasing for next generation communication systems. Using diversity, a simple detection can be formed at the cost of capacity reduction. Thus, spatial multiplexing is recommended for achieving significant capacity gain with improved transmission rate. However, demultiplexing at the receiver is still an issue for efficient system design.

A variety of detection techniques exist in literature to improve the system performance [2–7]. Zero forcing (ZF) and minimum mean square error (MMSE) are simple detection techniques although they lead to the reduced error performance [3, 4]. Maximum likelihood (ML) detection provides an optimal error rate performance. However, the hardware complexity increases with the increase in transmitter antennas and modulation order. Therefore, ordered successive interference cancellation (OSIC) is recommended with MMSE for improving the error rate performance [5]. In [5–8], efficient approaches have been investigated to develop OSIC and to reduce the complexity of receiver for Vertical-Bell Laboratories Layered Space-Time (V-BLAST) system. In [9], antenna selection has been used for OSIC detection to enhance the error rate performance. Also, a near-optimal selection technique has been proposed to decrease the complexity without substantial performance reduction. Different detection ordering such as signal-to-noise ratio (SNR) based, column norm based and, signal-to-interference plus noise ratio (SINR) based ordering has been investigated. It is examined that postdetection SINR based ordering achieves the best performance among all three ordering techniques [10]. Hence, MMSC-OSIC detection is used in this paper to achieve the improved error rate performance with adaptable complexity level in composite fading scenario. Identifying an efficient detection technique is still under investigation in MIMO systems.

Various small scale multipath channel models including Rayleigh, Rician, Nakagami-m, and Weibull have been recommended to analyze the wireless system performance. Large scale fading can be modeled using log-normal (LN) or Gamma distribution. Both the small scale and large scale fading effects can be observed simultaneously due to rapidly changing wireless environment. Therefore, channel modeling of composite fading which comprises multipath as well as shadowing effects is imperative to figure out. Also, it is important to resolve a number of practical problems with interference effects in MIMO wireless communications. The composite channel models such as Rayleigh-LN, shadowed-Rician, Gamma-Gamma, K, generalized-K, correlated shadowed- [11–17], and comparatively new-fangled Weibull-Gamma (WG) [18, 19] are widely used to incorporate both small scale fading and shadowing. Gamma models are simpler and more accurate than LN models and hence preferred in recent approaches. The WG composite distribution is appropriate for MIMO system design in the present wireless scenario due to its extensive flexibility, experimental efficiency, and analytical conformity. Weibull fading model is adaptive for modeling severe and nonsevere multipath fading conditions, Gamma model is used for shadowing, and collectively WG composite fading model is formed to quantify both the effects. WG fading model is a generic model and it includes Rayleigh-Gamma or K and exponential-Gamma distributions as its special cases. Consequently, this model approximates several other fading models [11, 18].

Both higher and lower modulation orders have been considered in multipath fading with the existence of AWGN [20–22]. However, the deviation of actual noise is possible; therefore, generic noise model is required. In power line communication (PLC), the system performance is extensively affected by additive and multiplicative power line noises. Further, the additive noise is categorized into background noise and impulsive noise. The background noise and impulsive noise models follow Nakagami-m and Middleton class A distribution, respectively, although the multiplicative PLC noise induces fading in the received signal power. Also, the system error rate performance has been evaluated in the presence of such noise scenarios [23–25]. The generalized Gaussian distribution (GGD) is the emerging research interest for modeling different noise effects. This generic noise model considers various forms of noise such as impulsive, Gamma, Gaussian, Laplacian [26]. Rectangular quadrature amplitude modulation (QAM) is defined by a combination of in-phase and quadrature phase pulse amplitude modulation (PAM) signals. This modulation technique has been used to compute average symbol error probability (ASEP) using Gaussian -function in composite fading scenario perturbed by additive white generalized Gaussian noise (AWGGN) [27].

To the best of our knowledge, the MIMO system performance has not been evaluated for composite WG fading channel with the consideration of generalized noise. To improve the MIMO system performance spatial multiplexing (SM) is used with efficient detection technique, that is, MMSE-OSIC. The expressions derived for ASEP in [27] are analyzed again which were previously limited to single-input single-output (SISO) system. To achieve the high data rate of the wireless link, higher order modulation techniques are preferred, although they are less flexible to noise and interference. Exact analytical expressions are computed in terms of Fox-H function (FHF) for 16-QAM in SM-MIMO WG fading subject to AWGGN. Two special cases of AWGGN, namely, AWGN and Laplacian noise, are considered. AWGN is popularly known and Laplacian noise has also achieved attention in the signal processing and wireless systems to model impulsive noise. In addition, the variations of fading and shadowing parameters are also illustrated.

The rest of the paper is organized as follows. In Section 2, system model and proposed detection technique are described. Section 3 provides the simulation results and analysis for ASEP of SM-MIMO in WG fading subject to generic noise along with specified detection technique. Finally, the paper is concluded in Section 4.

2. System and Channel Model

Consider a MIMO system having antennas, where and denote the number of transmit and receive antennas, respectively. is the transmit signal modulated by 16-QAM and multiplied by a composite flat fading channel envelope . Previously, the noise is generally considered as AWGN. However, the AWGGN noise is assumed here with zero mean and variance . MIMO system model is represented as

The probability density function (PDF) of the AWGGN noise is described in [28, Equation ()] over as

where and denote shaping parameter of and mean, respectively . Furthermore, the coefficient can be defined by normalizing the noise power using normalizing coefficient with respect to as

where denotes Gamma function, , , is the power per positive frequency of AWGGN, and is expectation operator.

The random variable of AWGGN distribution rigorously depends on its shaping parameter . This distribution reveals a superior fit to the quantified noise statistics with the changing physical channel conditions and forms various noise categories as special cases of AWGGN. For , , , and , it characterizes the prominent Gaussian, Laplacian, Gamma, and impulsive noise, respectively. Consequently, for and , the statistical properties and accurate simulation technique have been developed in the presence of AWGGN [29].

The PDF of the received signal envelope given in [30] is simplified to form WG distribution defined over and represented as

where is extended incomplete Gamma function, which is given as , , [28, Equation ()]. and denote fading figure and shadowing shaping parameter, respectively, where , , and () is the average power of the received signal envelope.

The SNR for received symbols in the presence of AWGGN follows WG PDF, which is defined over ). The average SNR per symbol is , where is the energy of transmitted symbols. By exchanging variables, PDF is represented in the form of SNR as

Equation (5) can be expressed in terms of the FHF using [28, Equation ()], [31, Equations (), () and ()] and simplifying [32, Equation ()] as

where

where is FHF defined in [31, Equation ()], [33] and represents Mellin-Barnes contour.

3. MMSE-OSIC Detection

MMSE detection technique maximizes the postdetection SINR by minimizing mean-square error (MSE). In OSIC, SINR based ordering improves the performance of linear detection technique by maximizing SINR. This technique maintains the low complexity for designing hardware. It holds number of linear receivers in which each receiver classifies one of the parallel data streams with detected signal components. These signal components are successively canceled from the received signal at each stage [34].

The MMSE detection technique explained in [10] offers the 1st estimated stream with the 1st row vector of the MMSE weight matrix which is given by

where denotes the Hermitian operator. In MMSE detection, to determine the required statistical information of , the th row vector of is obtained by

where is the th column vector of the channel matrix and is the Frobenius norm of matrix. denotes the th order detected symbol which depends on the order of detection; hence, this symbol may be different from the transmit signal at the th antenna. The sliced value of is given by . The remaining signal is represented by

Assume that , , , is accurately generated. When it means the interference is canceled and can be estimated. When then error propagation takes place. The entire performance of OSIC technique is affected by the order of detection. The erroneous outcomes in the previous stages cause the error propagation.

Primarily, signals containing a higher postdetection SINR are detected in SINR based ordering. The linear MMSE detection with the postdetection SINR is represented by

Once the selection of the SINR values using is completed, the corresponding layer with the highest SINR is chosen. The second detected symbol is selected by canceling the interference due to first detected symbol from the received signals. If th symbol is canceled first, then of (8) is transformed into (12) by deleting the channel gain vector as per th symbol:

Again, is calculated after substituting of (8) by (12). Then, SINR values (i.e., ) are computed by selecting the symbol containing highest SINR. After canceling the next symbol with the highest SINR, the same process is continued with the remaining signal. The total number of calculated SINR values is generated by .

The OSIC technique can offer diversity order greater than for all symbols. Following the ordering approach, the diversity order of the first detected symbol is also greater than . Nevertheless, the diversity order of remaining symbols depends on whether the previously detected symbols are exact or not. Suppose all the symbols are exact; then, the diversity order of the th detected symbol is . The th detected symbol is different from the one transmitted from the th transmit antenna. Since the ordering is based on SINR for MMSE detection, therefore, (11) is used to improve the ASEP performance.

4. Average Symbol Error Probability for -QAM

The symbol error probability (SEP) has been given in [35, Equation ()] for QAM in the presence of AWGN. The formation of -QAM signal constellation is given by two independent in-phase and quadrature -ary PAM signals, where -ary PAM and -ary PAM are in-phase and quadrature signals, respectively, and . Given that, GGD and Gaussian distribution demonstrate the identical symmetry properties. According to [27], the identical symmetry properties can be used to define SEP of -QAM, given as

where , , and and denote the decision distances for in-phase and quadrature phase components, respectively. is generalized- function for defined in [32] as

In [32, A.5], the representation of (14) in the form of FHF using [33, Equation (), ()] is given as

The ASEP is obtained by averaging the conditional SEPs in (13) under slow fading conditions over the PDF of . Then, is represented by

where

It is difficult to formulate and using the conventional expressions of WG distribution and GGD. Therefore, alternative expressions (6) and (15) are used to compute simplified analytical expressions for and and then resulting expression for the ASEP. In (17), consists of an integral including the product of two FHFs which is comparable to that of [32] considering the normalized value of fading shaping factor and severity of shadowing. Unlike [32], we prefer an efficient SINR based ordering for MMSE detection to improve the error rate performance of MIMO system. Using [27] and [31, Equation ()], can be represented in the form of FHF by a closed form expression given as

Substituting (6) and (15) in (18), an integral which includes the product of three FHFs is used to describe (20). Then, using [36, Equation ()], is expressed in terms of the FHF of two variables known as the bivariate Fox-H function (BFHF).

Substituting (19) and (20) in (16), the ASEP of -QAM is computed. This ASEP expression is given for rectangular (), square () QAM in arbitrary WG fading with AWGGN. Consequently, it maintains substantial range of noise and fading parameters. The commonly considered noise cases of AWGGN in composite fading scenario are as follows.

Case 1 (WG fading with Laplacian noise). The first special case of AWGGN emerges when , and the noise is considered Laplacian. Taking , is represented as

Using [31, 36], FHF and BFHF functions are well explored and utilized to make a simplified form of (21) by reducing number of terms in as

Similarly, can be written as

Using [36, Equation ()], for the description of BHFH and [28, Equations () and ()], is represented as

Equations (22) and (24) are used to calculate ASEP when Laplacian noise is present.

Case 2 (WG fading with AWGN). For , (19) and (20) can be rearranged to find the ASEP in AWGN environment. Again, the expressions for and are reduced as

Here, eliminates the effect of shadowing. When , WG distribution follows Weibull distribution and , converted WG distribution into Rayleigh distribution [30, Table ]. Thus, the fading scenario can be changed by setting the parameters and .

5. Simulation Results and Analysis

To evaluate the MIMO system performance, 16-QAM modulation is used as a function of SNR for the generalized case of noise. Therefore, distinct values of and arbitrary values of and are taken into consideration. The in-phase-to-quadrature phase decision distance ratio is represented as . For this case, the average total energy per symbol is given as and hence [27]. Taking a fixed , the identical average energies of the in-phase and quadrature signals are obtained.

The impact of the parameter is emerged to observe the system performance. When , the most favorable case occurs; this implies that the in-phase and quadrature distance are identical for both the Laplacian and Gaussian noise. For , same energy is obtained between the in-phase and quadrature signal; thus the system performance is reduced with a small amount, that is, approximately 1 dB SNR reduction for large SNRs. When the quadrature signal contains 10.5 times the average energy of the in-phase signal, for this instant loss is more essential as it gets approximately 4 dB SNR loss for large SNRs, comparative to the aforementioned case, where .

Firstly, composite WG fading is considered in Laplacian noise environment. To obtain Weibull and Rayleigh fading, the parameters are settled to and , , respectively. Figure 1 depicts the ASEP as a function of average SNR per symbol for both Gaussian and Laplacian cases of noise. In addition, distinct values of and are chosen to determine the severity of fading. Analytical results presented in this paper by (19) and (20) demonstrate the perfect match of the simulation results. The performance of the system is improved by increasing both the parameters and . Results shown in Figure 1 illustrate that the ASEP performance in Laplacian noise is superior to that of Gaussian noise for lower SNR or less than 15 dB SNR. However, for high SNR, less fading (), the situation is upturned and ASEP performance improves in the Gaussian noise compared with Laplacian noise. For severe fading (), Laplacian noise offers better results than Gaussian noise.

Figure 1: ASEP of 2 × 2 MIMO system using MMSE-OSIC over WG fading channel subject to Laplacian and Gaussian noise.

Afterward, Weibull and Rayleigh fading which are the special cases of WG fading are considered. In Figure 2, Rayleigh fading case is taken into account. In this case, system gives superior performance by diminishing , which validates the previous result in which the Laplacian noise gives better performance than the Gaussian noise in severe fading. It is previously mentioned that large fading parameter refers to less fading. In Figure 3, the ASEP is demonstrated as a function of the SNR per QAM symbol in Weibull fading environment () with AWGGN for which . In this case, when the less fading condition occurs, the different two regions are investigated. At low SNR, the ASEP decreases with increasing and at high SNR it improves by increasing .

Figure 2: ASEP of 2 × 2 MIMO system for 16-QAM using MMSE-OSIC over Rayleigh fading channel with arbitrary values of .

Figure 3: ASEP of 2 × 2 MIMO system for 16-QAM using MMSE-OSIC detection over Weibull fading channel () with arbitrary values of .

6. Conclusion

This paper evaluates the ASEP performance of MIMO system in composite WG fading environment subject to AWGGN. Analytical expressions for ASEP are derived using 16-QAM consisting of two independent in-phase and quadrature signals of PAM. The MMSE-OSIC detection is used to improve the error rate performance of MIMO system. It is concluded from the results that the ASEP performance in Laplacian noise is better than that of Gaussian noise for low SNR. However, in less fading, performance is degraded for high SNR and improved error performance is obtained in Gaussian noise compared with Laplacian noise. In severe fading, improved error rate performance can be achieved in the presence Laplacian noise compared with Gaussian noise. In Rayleigh fading case, the system gives superior performance for low noise shaping parameter . This result again proves that the Laplacian noise gives better performance than the Gaussian noise in severe fading. In Weibull fading, the different two regions are inspected for lower amount of fading or large fading parameter. Moreover, the ASEP reduces with at low SNR and it increases by at high SNR. Simulation results validate the analytical results.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

References

E. Biglieri, R. Calderbank, A. Constantinides, A. Goldsmith, A. Paulraj, and H. V. Poor, MIMO Wireless Communications, Cambridge University Press, Cambridge, UK, 2007.

Y. Lee, H.-C. Shih, C.-S. Huang, and J.-Y. Li, “Low-complexity MIMO detection: a mixture of ZF, ML and SIC,” in Proceedings of the 19th International Conference on Digital Signal Processing (DSP '14), pp. 263–268, Hong Kong, August 2014.

C. Siriteanu, S. D. Blostein, A. Takemura, H. Shin, S. Yousefi, and S. Kuriki, “Exact MIMO zero-forcing detection analysis for transmit-correlated rician fading,” IEEE Transactions on Wireless Communications, vol. 13, no. 3, pp. 1514–1527, 2014.

Y. Jiang, M. K. Varanasi, and J. Li, “Performance analysis of ZF and MMSE equalizers for MIMO systems: an in-depth study of the high SNR regime,” IEEE Transactions on Information Theory, vol. 57, no. 4, pp. 2008–2026, 2011.

H. Zhang, H. Dai, and B. L. Hughes, “Analysis on the diversity-multiplexing tradeoff for ordered MIMO SIC receivers,” IEEE Transactions on Communications, vol. 57, no. 1, pp. 125–133, 2009.

T.-H. Liu, “Some results for the fast MMSE-SIC detection in spatially multiplexed MIMO systems,” IEEE Transactions on Wireless Communications, vol. 8, no. 11, pp. 5443–5448, 2009.

Y. Song, C. Liu, and F. Lu, “Lattice reduction-ordered successive interference cancellation detection algorithm for multiple-input–multiple-output system,” IET Signal Processing, vol. 9, no. 7, pp. 553–561, 2015.

Y. Song, C. Liu, and F. Lu, “The optimal MMSE-based OSIC detector for MIMO system,” IEICE Transactions on Communications, vol. E99-B, no. 1, pp. 232–239, 2016.

Y. Bae and J. Lee, “Low complexity antenna selection for V-BLAST systems with OSIC detection,” EURASIP Journal on Wireless Communications and Networking, vol. 2011, article 6, 7 pages, 2011.

Y. S. Cho, J. Kim, W. Y. Yang, and C. G. Kang, MIMO-OFDM Wireless Communications with MATLAB, John Wiley & Sons, 2010.

P. M. Shankar, Fading and Shadowing in Wireless Systems, Springer, New York, NY, USA, 2012.

T. Hashem and M. I. Islam, “Performance analysis of MIMO link under fading channels,” in Proceedings of the 17th International Conference on Computer and Information Technology (ICCIT '14), pp. 498–503, Dhaka, Bangladesh, December 2014.

A. Abdi and M. Kaveh, “K distribution: an appropriate substitute for Rayleigh-lognormal distribution in fading-shadowing wireless channels,” Electronics Letters, vol. 34, no. 9, pp. 851–852, 1998.

M. Matthaiou, N. D. Chatzidiamantis, G. K. Karagiannidis, and J. A. Nossek, “On the capacity of generalized-K fading MIMO channels,” IEEE Transactions on Signal Processing, vol. 58, no. 11, pp. 5939–5944, 2010.

M. R. Bhatnagar and M. K. Arti, “On the closed-form performance analysis of maximal ratio combining in shadowed-rician fading LMS channels,” IEEE Communications Letters, vol. 18, no. 1, pp. 54–57, 2014.

M. R. Bhatnagar, “On the sum of correlated squared shadowed random variables and its application to performance analysis of MRC,” IEEE Transactions on Vehicular Technology, vol. 64, no. 6, pp. 2678–2684, 2015.

M. R. Bhatnagar and M. K. Arti, “Performance analysis of AF based hybrid satellite-terrestrial cooperative network over generalized fading channels,” IEEE Communications Letters, vol. 17, no. 10, pp. 1912–1915, 2013.

P. S. Bithas, “Weibull-gamma composite distribution: alternative multipath/shadowing fading model,” Electronics Letters, vol. 45, no. 14, pp. 749–751, 2009.

K. Tiwari and D. S. Saini, “SER improvisation of MIMO-MRC system over Weibull-Gamma fading channel,” in Proceedings of the International Conference on Signal Processing and Communication (ICSC '15), pp. 70–73, IEEE, Noida, India, March 2015.

N. C. Beaulieu, “A useful integral for wireless communication theory and its application to rectangular signaling constellation error rates,” IEEE Transactions on Communications, vol. 54, no. 5, pp. 802–805, 2006.

G. K. Karagiannidis, “On the symbol error probability of general order rectangular QAM in nakagami-m Fading,” IEEE Communications Letters, vol. 10, no. 11, pp. 745–747, 2006.

H. A. Suraweera and J. Armstrong, “A simple and accurate approximation to the SEP of rectangular QAM in arbitrary Nakagami-m fading channels,” IEEE Communications Letters, vol. 11, no. 5, pp. 426–428, 2007.

A. Mathur and M. R. Bhatnagar, “PLC performance analysis assuming bpsk modulation over nakagami-m additive noise,” IEEE Communications Letters, vol. 18, no. 6, pp. 909–912, 2014.

A. Mathur, M. R. Bhatnagar, and B. K. Panigrahi, “PLC Performance analysis over Rayleigh fading channel under Nakagami-m additive noise,” IEEE Communications Letters, vol. 18, no. 12, pp. 2101–2104, 2014.

A. Mathur, M. R. Bhatnagar, and B. K. Panigrahi, “Performance evaluation of PLC under the combined effect of background and impulsive noises,” IEEE Communications Letters, vol. 19, no. 7, pp. 1117–1120, 2015.

R. Viswanathan and A. Ansari, “Distributed detection of a signal in generalized Gaussian noise,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 37, no. 5, pp. 775–778, 1989.

H. Soury, F. Yilmaz, and M.-S. Alouini, “Error rates of M-PAM and M-QAM in generalized fading and generalized gaussian noise environments,” IEEE Communications Letters, vol. 17, no. 10, pp. 1932–1935, 2013.

M. A. Chaudhry and S. M. Zubair, On a Class of Incomplete Gamma Functions with Applications, Chapman & Hall/CRC, Boca Raton, Fla, USA, 1st edition, 2002.

M. Nardon and P. Pianca, “Simulation techniques for generalized Gaussian densities,” Working Papers 145, Departement of Applied Mathematics, Università Ca' Foscari Venezia, 2006.

F. Yilmaz and M.-S. Alouini, “A new simple model for composite fading channels: second order statistics and channel capacity,” in Proceedings of the 7th International Symposium on Wireless Communication Systems (ISWCS '10), pp. 676–680, IEEE, York, UK, September 2010.

A. A. Kilbas and M. Saigo, H-Transforms: Theory and Applications (Analytical Methods and Special Functions), CRC Press, New York, NY, USA, 1st edition, 2004.

H. Soury, F. Yilmaz, and M.-S. Alouini, “Average bit error probability of binary coherent signaling over generalized fading channels subject to additive generalized Gaussian noise,” IEEE Communications Letters, vol. 16, no. 6, pp. 785–788, 2012.

A. P. Prudnikov, Y. A. Brychkov, and O. I. Marichev, Integrals and Series: More Special Functions, vol. 3, Gordon and Breach Science Publishers, New York, NY, USA, 1990.

G. D. Golden, C. J. Foschini, R. A. Valenzuela, and P. W. Wolniansky, “Detection algorithm and initial laboratory results using V-BLAST space-time communication architecture,” Electronics Letters, vol. 35, no. 1, pp. 14–15, 1999.

N. C. Beaulieu, “A useful integral for wireless communication theory and its application to rectangular signaling constellation error rates,” IEEE Transactions on Communications, vol. 54, no. 5, pp. 802–805, 2006.

P. K. Mittal and K. C. Gupta, “An integral involving generalized function of two variables,” Indian Academic Science, vol. 75, pp. 117–123, 1972.