Open Access Library Journal
Vol.04 No.05(2017), Article ID:75968,10 pages
10.4236/oalib.1103574

Interference Cancellation Method without Feedback Amount for Three Users Interference Channel

Xinji Tian, Haotian Zhang, Wenjie Jia

School of Physics & Electronic Information Engineering, Henan Polytechnic University, Jiaozuo, China

Copyright © 2017 by authors and Open Access Library Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: March 30, 2017; Accepted: May 1, 2017; Published: May 4, 2017

ABSTRACT

In this paper, we investigate how to cancel multi-user interference without feedback amount over multiple-input multiple-output (MIMO) interference channel using space-time code. Three transmitter-receiver pairs are considered and each node is equipped with four antennas. Space-time codeword with dimension 4 ´ 4 is designed for three transmitters, respectively. The number of interfering codeword is reduced through the design of codeword. Then, multi-user interference is cancelled by a linear operation on the received signals. Compared with the existing scheme for the same scene, the transmission efficiency is improved, and the unidirectional cooperative link is not required. Simulation results demonstrate the validity of theoretical analysis.

Subject Areas:

Information and Communication Theory and Algorithms

Keywords:

Interference Channel, Space-Time Codeword, Multi-User Interference, Interference Cancellation, Transmission Efficiency

1. Introduction

Reliability and effectiveness are important factors for wireless communications. Multiple-input multiple-output (MIMO) systems have the ability to improve both and hence, are immensely popular. In theory, MIMO systems can achieve extremely high spectral efficiency by transmitting independent data stream simultaneously. However, multiple users, transmitting different codewords to multiple receivers simultaneously using the same band, are not able to collaborate with each other [1] . There is serious interference at each receiver, and the effectiveness becomes limited. Therefore, it is quite important to design transmission techniques to eliminate multi-user interference [2] .

Several kinds of interference cancellation methods have been conducted, such as pre-coding at the transmitters, multi-user detection of the receivers, and interference alignment [3] [4] . Interference alignment methods align the interference caused at each receiver from all the transmitters along the same direction while the desired signals occupy linearly independent signal dimensions. In this case, the desired signals can be decoded free from interference by zero-forcing algorithm.

A popular way to approximating the sum-capacity for interference alignment is using the concept of degrees of freedom (DoF), in terms of transmission efficiency. Many works have been conducted to evaluate DoF [5] [6] . However, effectiveness is contract to reliability in communication systems. With high DoF, reliability in terms of the diversity gain is needed to be improved. Therefore, it is meaningful to combine interference alignment and space-time code, which can offer diversity gain [7] [8] . [9] proposes a new interference alignment scheme using Alamouti code for two users MIMO X channel, in which zero vector is introduced into each codeword. The diversity gain of [9] is 2, which is brought by Alamouti code. Space-time codeword with dimension 4 × 4 is designed and zero vectors are introduced into each codeword over two users MIMO X channel with four antennas at each user, in which the diversity gain is 4 [10] . The multi- user interference is mitigated through interference alignment and linear operation on the received signals. Motivated by [9] , interference alignment and Alamouti code are introduced into MIMO interference channel with three transmitter- receiver pairs [11] . The multi-user interference is mitigated through interference alignment and unidirectional cooperative link. However, channel state information (CSI) is required at the transmitters. It is very difficult to estimate the CSI in a mobility environment. Therefore, it is crucial to consider transmission scheme combined interference alignment and space-time code without CSI.

Under this motivation, interference cancellation method without feedback is proposed for MIMO interference channel with three transmitter-receiver pairs where each node is equipped with four antennas. Space-time codeword with dimension 4×4 is designed for three transmitters, respectively. Each codeword includes eight independent modulated symbols and zero matrices are introduced into each codeword, through which the number of interfering codeword is reduced. Then, multi-user interference is cancelled by a linear operation on the received signals. The transmission efficiency of the proposed scheme is higher than that of [11] . Theoretical analysis shows that the diversity of the proposed scheme is the same as that of [11] .

Notation: The operators ( ) and ( ) H stand for complex conjugation and conjugation transpose, respectively. denotes Frobenius norm of the enclosed term. E ( ) is the expectation operations.

2. System Model

The system model is shown in Figure 1. There are three dedicated transmitter- receiver pairs, where each transmitter has an individual message to be communicated to its respectively receiver. All the transmitters and receivers are assumed to have four antennas. We assume CSI is not available at each transmitter.

a k , b k and c k , which take values from a set such that E [ | a k | 2 ] = E [ | b k | 2 ] = E [ | c k | 2 ] = 1 , k = 1 , 2 , , 8 , represent the modulated signals at transmitter 1, transmitter 2 and transmitter 3, respectively. a k , b k and c k are encoded into 2 × 2 space-time codeword, denoted by A i , B i and C i , respectively, i = 1 , 2 . A i , B i and C i are expressed as

A i = [ a 4 i 3 + e j θ a 4 i 1 a 4 i 2 e j θ a 4 i a 4 i 2 + e j θ a 4 i a 4 i 3 + e j θ a 4 i 1 ] B i = [ b 4 i 3 + e j θ b 4 i 1 b 4 i 2 e j θ b 4 i b 4 i 2 + e j θ b 4 i b 4 i 3 + e j θ b 4 i 1 ] C i = [ c 4 i 3 + e j θ c 4 i 1 c 4 i 2 e j θ c 4 i c 4 i 2 + e j θ c 4 i c 4 i 3 + e j θ c 4 i 1 ] (1)

where, j = 1 , 0 < θ < 2 π . The elements of A i , the elements of B i and the elements of C i are nonzero by the proper value of θ . Each element of A i , each element of B i and each element of C i includes two independent symbols, with the purpose of improving transmission efficiency. If we see the elements of A i , the elements of B i and the elements of C i as the modulated signals, i = 1 , 2 , A i , B i and C i have an Alamouti structure. The motive for such design is that the effective channel matrix of Alamouti code has orthogonal columns, which can be used to reduce decoding complexity.

At transmitter 1, A 1 and A 2 are combined according to a certain rule to get codeword A . Similarly, B 1 and B 2 are combined to obtain codeword B while C 1 and C 2 are combined to obtain codeword C . A , B and C are given by

Figure 1. The system model of the proposed scheme.

A = [ A 1 0 2 A 2 0 2 ] , B = [ B 1 0 2 0 2 B 2 ] , C = [ 0 2 C 1 0 2 C 2 ] (2)

where, 0 2 represents 2 × 2 zero matrix. 0 2 is introduced because it can be used to mitigate partial multi-user interference, as shown later. A , B and C are the desired codewords for the three receivers, respectively.

As shown in Figure 1, three transmitters send A , B and C on the same frequency band at the same time. The received signals at receiver i, denoted by Y i , are expressed as

Y 1 = H 11 A + H 12 B + H 13 C + N 1 (3)

Y 2 = H 21 A + H 22 B + H 23 C + N 2 (4)

Y 3 = H 31 A + H 32 B + H 33 C + N 3 (5)

where, i = 1 , 2 , 3 , H i k is the channel from transmitter k to receiver i and N i is the complex additive white Gaussian noise at receiver i. The dimension of H i k and N i are both 4 × 4 . To avoid degenerate channel conditions, we assume that the channel coefficient values are drawn i.i.d. from a continuous distribution and the absolute value of all the channel coefficients is bounded between a nonzero minimum value and a finite maximum value.

The term H 11 A in Equation (3) is the desired signals for receiver 1 and the term H 12 B + H 13 C in Equation (3) is the interference signals for receiver 1. The term H 22 B in Equation (4) is the desired signals for receiver 2 and the term H 21 A + H 23 C in Equation (4) is the interference signals for receiver 2. The term H 33 C in Equation (5) is the desired signals for receiver 3 and the term H 31 A + H 32 B in Equation (5) is the interference signals for receiver 3.

3. Interference Cancellation Method

We describe the interference cancellation method at R 1 . It can be similarly used to mitigate interference at R 2 and to mitigate interference R 3 .

3.1. Remove the Undesired Codeword

Let H 11 = [ H 1 H 2 ] , H 12 = [ H 3 H 4 ] , H 13 = [ H 5 H 6 ] , Y 1 = [ Y 11 Y 21 ] and N 1 = [ N 11 N 21 ] . The dimension of H l , N i 1 and Y i 1 are both 4 × 2 , l = 1 , 2 , , 6 i = 1 , 2 . Then, (3) can be expressed as

[ Y 11 Y 21 ] = [ H 1 A 1 + H 2 A 2 + H 3 B 1 H 4 B 2 + H 5 C 1 + H 6 C 2 ] + [ N 11 N 21 ] (6)

(6) is equivalent to the following two equations,

Y 11 = H 1 A 1 + H 2 A 2 + H 3 B 1 + N 11 (7)

Y 21 = H 4 B 2 + H 5 C 1 + H 6 C 2 + N 21 (8)

The codeword of transmitter 1 and the partial codeword of transmitter 2 are encompassed by Y 11 . The codeword of transmitter 3 and the partial codeword of transmitter 2 are encompassed by Y 21 . So, the number of interfering codeword is reduced by the introduction of zero matrices. However, both the desired codeword for receiver 1 and the undesired codeword for receiver 1 in terms of B 1 are included by Y 11 . In what follows, we represent how to cancel the undesired codeword B 1 .

Let H 1 = [ h 11 1 h 12 1 h 21 1 h 22 1 h 31 1 h 32 1 h 41 1 h 42 1 ] , H 2 = [ h 11 2 h 12 2 h 21 2 h 22 2 h 31 2 h 32 2 h 41 2 h 42 2 ] and H 3 = [ h 11 3 h 12 3 h 21 3 h 22 3 h 31 3 h 32 3 h 41 3 h 42 3 ] . Use y i j and

n i j to denote the element of Y 11 and N 11 , i = 1 , 2 , 3 , 4 , j = 1 , 2 . The elements of Y 11 can be effectively written as

y 1 = G 11 [ a 1 + e j θ a 3 a 2 + e j θ a 4 ] + G 12 [ a 5 + e j θ a 7 a 6 + e j θ a 8 ] + G 13 [ b 1 + e j θ b 3 b 2 + e j θ b 4 ] + n 1 (9)

y 2 = G 21 [ a 1 + e j θ a 3 a 2 + e j θ a 4 ] + G 22 [ a 5 + e j θ a 7 a 6 + e j θ a 8 ] + G 23 [ b 1 + e j θ b 3 b 2 + e j θ b 4 ] + n 2 (10)

y 3 = G 31 [ a 1 + e j θ a 3 a 2 + e j θ a 4 ] + G 32 [ a 5 + e j θ a 7 a 6 + e j θ a 8 ] + G 33 [ b 1 + e j θ b 3 b 2 + e j θ b 4 ] + n 3 (11)

y 4 = G 41 [ a 1 + e j θ a 3 a 2 + e j θ a 4 ] + G 42 [ a 5 + e j θ a 7 a 6 + e j θ a 8 ] + G 43 [ b 1 + e j θ b 3 b 2 + e j θ b 4 ] + n 4 (12)

where, G i j = [ h i 1 j h i 2 j h i 2 j h i 1 j ] , y i = [ y i 1 y i 2 ] , n i = [ n i 1 n i 2 ] , i = 1 , 2 , 3 , 4 , j = 1 , 2 , 3 . y i

and n i are the effective received signals and the effective noise at the i-th an-

tenna, respectively. [ a 1 + e j θ a 3 a 2 + e j θ a 4 ] , [ a 5 + e j θ a 7 a 6 + e j θ a 8 ] and [ b 1 + e j θ b 3 b 2 + e j θ b 4 ] are the effec-

tive transmitted signals of the three transmitters. G i j is the effective channel matrix corresponding to the effective transmitted signals, i = 1 , 2 , 3 , 4 ,

j = 1 , 2 , 3 . G i j , having orthogonal columns, satisfies G i j H G i j G i j 2 = I 2 , where I 2 is

unit matrix with dimension 2 × 2 . The third terms at the right side of Equation (9-12) are the interference signals, which can be mitigated using interference cancellation method as follows

z 1 = G 13 H y 1 G 13 2 G 23 H y 2 G 23 2 = P 1 [ a 1 + e j θ a 3 a 2 + e j θ a 4 ] + P 2 [ a 5 + e j θ a 7 a 6 + e j θ a 8 ] + w 1 (13)

z 2 = G 33 H y 3 G 33 2 G 43 H y 4 G 43 2 = P 3 [ a 1 + e j θ a 3 a 2 + e j θ a 4 ] + P 4 [ a 5 + e j θ a 7 a 6 + e j θ a 8 ] + w 2 (14)

where, P i = G 13 H G 1 i G 13 2 G 23 H G 2 i G 23 2 , P j = G 33 H G 3 ( j 2 ) G 33 2 G 43 H G 4 ( j 2 ) G 43 2 ,

w 1 = G 13 H n 1 G 13 2 G 23 H n 2 G 23 2 , w 2 = G 33 H n 3 G 33 2 G 43 H n 4 G 43 2 , i = 1 , 2 , j = 1 , 2 . z i and w i are

the signals and the noise after linear operation, respectively, i = 1 , 2 . P i is effective channel matrix after linear operation, i = 1 , 2 , 3 , 4 . The dimension of z i and the dimension w i are both 2 × 1 , i = 1 , 2 . The dimension of P i is 2 × 2 , i = 1 , 2 , 3 , 4 . Symbols transmitted from transmitter 2 and symbols transmitted from transmitter 3 are not included in (13-14). So multi-user interference is mitigated through linear operation on the received signals.

3.2. Separate the Desired Codewords

Even if multi-user interference is cancelled, the desired codeword at receiver 1 includes 8 independent modulated symbols, namely a 1 , a 2 , , a 8 , as shown in (13-14). If we decode a 1 , a 2 , , a 8 based on (13-14), the complexity of maximum likelihood (ML) algorithm is proportional to the 8th power of modulation order, which is extremely high. So, it is necessary to divide these symbols into several groups. Take operation on z 1 and z 2 to get,

x 1 = P 2 H z 1 P 2 2 P 4 H z 2 P 4 2 = Q 1 [ a 1 + e j θ a 3 a 2 + e j θ a 4 ] + w 3 (15)

x 2 = P 1 H z 1 P 1 2 P 3 H z 2 P 3 2 = Q 2 [ a 5 + e j θ a 7 a 6 + e j θ a 8 ] + w 4 (16)

where, Q 1 = P 2 H P 1 P 2 2 P 4 H P 3 P 4 2 , Q 2 = P 1 H P 2 P 1 2 P 3 H P 4 P 3 2 , w 3 = P 2 H w 1 P 2 2 P 4 H w 2 P 4 2 ,

w 4 = P 1 H w 1 P 1 2 P 3 H w 2 P 3 2 . a 1 , a 2 , a 3 and a 4 are encompassed by A 1 whereas

a 5 , a 6 , a 7 and a 8 are encompassed by A 2 . So, the elements of A 1 and the elements of A 2 are separated, as shown in (15-16). The decoding complexity can be reduced.

From the above analysis, not only multi-user interference is cancelled, but also the two desired codewords A 1 and A 2 are separated. Similar operation can be performed at the other two receivers. Note that neither feedback information nor the unidirectional cooperative link is required.

4. Decoding Method

In this section, we will show how to decode, taking receiver 1 as example.

G i j has an Alamouti structure, i = 1 , 2 , 3 , 4 , j = 1 , 2 , 3 . Due to the completeness of matrix addition, matrix multiplication, and the scalar multiplication of the Alamouti matrix, Q 1 and Q 2 still have the Alamouti structure. Let

Q 1 = [ q 1 q 2 q 2 q 1 ] and Q 2 = [ q 3 q 4 q 4 q 3 ] . Multiply both sides of (15) with Q 1 H to

get

Q 1 H x 1 = [ x 1 x 2 ] = [ | q 1 | 2 + | q 2 | 2 0 0 | q 1 | 2 + | q 2 | 2 ] [ a 1 + e j θ a 3 a 2 + e j θ a 4 ] + Q 1 H w 3 (17)

Multiply both sides of (16) with Q 2 H to get

Q 2 H x 2 = [ x 3 x 4 ] = [ | q 3 | 2 + | q 4 | 2 0 0 | q 3 | 2 + | q 4 | 2 ] [ a 5 + e j θ a 7 a 6 + e j θ a 8 ] + Q 2 H w 4 (18)

In (17), the channel corresponding to a 1 + e j θ a 3 is orthogonal to the channel corresponding to a 2 + e j θ a 4 . So we can decode a 1 + e j θ a 3 and a 2 + e j θ a 4 separately. Similarly, we can decode a 5 + e j θ a 7 and a 6 + e j θ a 8 separately as well. According to the analysis mentioned above, the decoding method is as follows.

Step1, obtain G 11 , G 12 , , G 43 from the channel matrices, and then compute P 1 , P 2 , P 3 and P 4 from G 11 , G 12 , , G 43 ;

Step 2, compute Q 1 = P 2 H P 1 P 2 2 P 4 H P 3 P 4 2 and Q 2 = P 1 H P 2 P 1 2 P 3 H P 4 P 3 2 through di-

rect multiplications and subtraction ;

Step 3, obtain y 1 , y 2 , y 3 and y 4 according to the received signals, as in (9-12);

Step 4, take operation on y 1 , y 2 , y 3 and y 4 to get z 1 and z 2 ,

z 1 = G 13 H y 1 G 13 2 G 23 H y 2 G 23 2 , z 2 = G 33 H y 3 G 33 2 G 43 H y 4 G 43 2 ;

Step 5, take operation on z 1 and z 2 to get x 1 = P 2 H z 1 P 2 2 P 4 H z 2 P 4 2 and

x 2 = P 1 H z 1 P 1 2 P 3 H z 2 P 3 2 , furthermore, computer [ x 1 x 2 ] = Q 1 H x 1 and [ x 3 x 4 ] = Q 2 H x 2

through direct multiplications;

Step 6, estimate a 1 and a 3 by

( a ^ 1 , a ^ 3 ) = arg a 1 , a 3 min [ x 1 ( | q 1 | 2 + | q 2 | 2 ) ( a 1 + e j θ a 3 ) ] , where a ^ 1 and a ^ 3 are the es-

timation of a 1 and a 3 , respectively;

Step 7, estimate a 2 and a 4 by

( a ^ 2 , a ^ 4 ) = arg a 2 , a 4 min [ x 2 ( | q 1 | 2 + | q 2 | 2 ) ( a 2 + e j θ a 4 ) ] , where a ^ 2 and a ^ 4 are

the estimation of a 2 and a 4 , respectively;

Step 8, estimate a 5 and a 6 by

( a ^ 5 , a ^ 6 ) = arg a 5 , a 6 min [ x 3 ( | q 1 | 2 + | q 2 | 2 ) ( a 5 + e j θ a 6 ) ] , where a ^ 5 and a ^ 6 are the

estimation of a 5 and a 6 , respectively;

Step 9, estimate a 7 and a 8 by

( a ^ 7 , a ^ 8 ) = arg a 7 , a 8 min [ x 4 ( | q 1 | 2 + | q 2 | 2 ) ( a 7 + e j θ a 8 ) ] , where a ^ 7 and a ^ 8 are the

estimation of a 7 and a 8 , respectively;

The ML decoding complexity in step (6-9) is proportional to the square of modulation order.

It is easy to see from (17) that the channel corresponding to a 1 + e j θ a 3 is | q 1 | 2 + | q 2 | 2 , which includes two independent path with gain | q 1 | 2 and | q 2 | 2 . So a 1 + e j θ a 3 has a diversity gain of two, which is equal to that of [11] . Similar arguments can be used to prove that other symbols achieve a diversity of two as well.

5. Simulation Results

6 independent modulated signals are sent over two time slots in [11] . Feedback information, namely 16 complex numbers, are required at the transmitters, which occupy 16/3 time slot if these three receivers adopt the same transmission and the same code mode as the transmitters. So, 22/3 time slot is required to send 6 independent modulated signals in [11] . The transmission efficiency is 3/22 symbol per antenna per time slot. 24 independent modulated signals are sent over four time slots in the proposed scheme. The transmission efficiency is 1/2 symbol per antenna per time slot, which is 11/3 times as that of [11] .

Performance of the proposed scheme and performance of the scheme in [11] are compared in Table 1. Two schemes have the same diversity gain. The advantage of the proposed scheme is that neither feedback information nor unidirectional cooperative link is required, and the transmission efficiency is improved. The disadvantage is that the decoding complexity is proportional to the square of modulation order, higher than that of [11] , which is proportional to the modulation order.

In Figure 2, we simulate the average BER curves of the proposed scheme and the scheme in [11] . We consider uncoded systems, in which all the underlying links experience Rayleigh fading channels. The data symbols are assumed to be 4QAM symbols and 16QAM symbols, generated independently at each transmitter. The average BER is determined through monte-carlo simulations. The

Figure 2. The BER curves of the two schemes.

Table 1. Performance of the two schemes.

slope of BER curves represents diversity gain, which is 2 in these two schemes. It can be seen from this figure that the slope of BER curves of the proposed scheme is almost the same as that of [11] , which demonstrate the validity of theoretical analysis.

6. Conclusion

In this paper, the multi-user interference is cancelled by the design of the codeword and linear operation on the received signals. Compared with the existing scheme for the same scene, the proposed scheme doesn’t require feedback information and doesn’t require unidirectional cooperative link, while keeping the same diversity gain. However, the scheme is limited to three transmitter-receiver pairs, where each node is equipped with four antennas. Future work on this scheme includes extending the application scene.

Acknowledgements

The authors would like to thank the reviewers for their detailed reviews and constructive comments, which have helped improve the quality of this paper. This work was supported in part by National Natural Science Foundation of China under Grant No. 51174263.

Cite this paper

Tian, X.J., Zhang, H.T. and Jia, W.J. (2017) Interference Cancellation Method without Feedback Amount for Three Users Interference Channel. Open Access Library Journal, 4: e3574. https://doi.org/10.4236/oalib.1103574

References

  1. 1. Gesbert, D., Shafi, M. and Shiu, D. (2013) From Theory to Practice: An Overview of Mimo Space-Time Coded Wireless Systems. IEEE Journal on Selected Areas in Communications, 21, 281-302.
    https://doi.org/10.1109/JSAC.2003.809458

  2. 2. Chen, X., Song, S.H. and Letaief, K.B. (2014) Interference Alignment in Dual-Hop MIMO Interference Channel. IEEE Transactions on Wireless Communications, 13, 1274-1283.
    https://doi.org/10.1109/TWC.2014.012114.121952

  3. 3. Lu, L. and Zhang, W. (2014) Blind Interference Alignment with Diversity in K-User Interference Channels. IEEE Transactions on Communication, 62, 2850-2859.
    https://doi.org/10.1109/TCOMM.2014.2333516

  4. 4. Maddah-Ali, M., Motahari, A. and Khandani, A. (2008) Communication over MIMO X Channels: Interference Alignment, Decomposition, and Performance Analysis. IEEE Transactions on Information Theory, 54, 3457-3470.
    https://doi.org/10.1109/TIT.2008.926460

  5. 5. Xie, B., Li, Y., Minn, H. and Nosratinia, A. (2013) Adaptive Interference Alignment with CSI Uncertainty. IEEE Transactions on Communications, 61, 792-801.
    https://doi.org/10.1109/TCOMM.2012.121112.110589

  6. 6. Chen, X., Song, S.H. and Letaief, K.B. (2014) Interference Alignment in Dual-Hop MIMO Interference Channel. IEEE Transactions on Wireless Communications, 13, 1274-1283.
    https://doi.org/10.1109/TWC.2014.012114.121952

  7. 7. Ganesan, A. and Rajan, B.S. (2013) Interference Alignment with Diversity for the 2×2 X Network with Three Antennas. arXiv:1304.1432[cs.IT]

  8. 8. Shi, L., Zhang, W. and Xia, X.G. (2013) Space-Time Block Code Designs for Two-User MIMO X Channels. IEEE Transactions on Communications, 61, 3806-3815.
    https://doi.org/10.1109/TCOMM.2013.080613.120972

  9. 9. Li, L.B. and Jafarkhani, H. (2013) Maximum-Rate Transmission with Improved Diversity Gain for Interference Networks. IEEE Transactions on Information Theory, 59, 5313-5330.
    https://doi.org/10.1109/TIT.2013.2264479

  10. 10. Ganesan, A. and Rajan, B.S. (2014) Interference Alignment with Diversity for the 2 × 2 X Network with Four Antennas. IEEE Transactions on Information Theory, 60, 3576-3592.
    https://doi.org/10.1109/TIT.2014.2313614

  11. 11. Zaki, A., Wang, C. and Rasmussen, L.K. (2013) Combining Interference Alignment and Alamouti Codes for the 3-User Mimo Interference Channel. IEEE Wireless Communications and Networking Conference, Shanghai, 3563-3567.
    https://doi.org/10.1109/wcnc.2013.6555138