Keywords

1 Introduction

In this work, the massive multiple-input multiple-output (mMIMO) system model is adopted with secure point to point communication based on symbol decomposition. Massive MIMO has several advantages over the single input single output model due to its high spectral efficiency [1]. The idea of maximizing the spectral efficiency of mMIMO by using higher order modulation symbols can lead to high peak-to-average power ratio (PAPR), which can lead to distortion and reduce amplification efficiency [2]. Naturally, higher order modulation signals have amplitude difference in the multilevel constellation points to high envelope fluctuations. Alternatively, higher order modulation symbols can be decomposed into low order modulation symbols to experience low envelope fluctuations which results in improved amplification efficiency over signal with high envelope fluctuations. In this method, we decompose sixty four Quadrature Amplitude Modulation (64-QAM) symbol into six Binary Phase Shift Keying (BPSK) symbols and amplify all the six components in six transmitting antenna and eventually beamforms all the components simultaneously through a single channel. The receiver receives all the components as a single 16-QAM symbol [3]. In this paper, we decompose 64-QAM symbol into six BPSK symbols instead of three QPSK symbols, since using BPSK will have low envelop as compared to using QPSK. Also, based on the results from [3] in terms of security aspects, we use BPSK.

This method involves multi-layer system: a multi-branch amplification circuit forms the first layer; each branch is connected in parallel to an antenna array forming the second layer; the third layer is composed of a spatial multiplexing circuit. This multi-layer approach can create complex channel matrix and self-interference in the received signal. Therefore, this requires Iterative Block Decision Feedback Equalization (IB-DFE) to reduce the self-interference in the signal detection block [4]. IB-DFE is a low complexity receiving technique because it does not require matrix inversion unlike mMIMO conventional receivers, improving estimation and the information estimate with each iteration [5, 6].

We consider that both the transmitter and receiver positions are dynamic and, in this condition, the communication channel will also be dynamic, thus it is imperative to estimate channel information at the receiver. If receiver does not know accurate channel knowledge, then the symbol components will for distorted higher order modulation symbols. A robust channel estimation technique is employed in his system model, in which, we superimpose pilot symbol for every single data symbol and estimate at the receiver [7, 8]. This technique is robust and gives accurate estimates but consumes higher energy as compared to other techniques and this excess energy can be harvested from the superimposed signal at the receiver by using power splitting based simultaneous wireless power and information transmission technique (PS-SWIPT) with an power splitting circuit [9, 10]. The physical layer security aspects of the decomposition-based transmitter were already explained in [3]. Thus, in this paper we focus on energy efficiency and channel estimation technique with IB-DFE receiver.

The paper is organized as follows: Sect. 2 describes the system model. Section 3 describes joint information and channel estimation technique. The performance analysis is illustrated with simulation results in Sect. 4. Conclusion and future works are given in Sect. 5.

1.1 Technological Innovation and Contribution for Life Improvement

From autonomous cars to health monitoring systems, many real time mobile applications depend on the resilience of the communication system. The presented system combines energy efficiency, information security and improved channel estimation for non-static channels. It means that the overall communication will be more robust compared with classic schemes, thus providing a solid based for the development of more and better products for life improvement. This paper also adopts green communication technology SWIPT to improve energy efficiency and opportunistically use the available excess energy from RF signals and thereby this system model is suitable for sustainable wireless communication system.

2 System Model

In this system model, a 64-QAM signal is decomposed to six BPSK polar components and transmitted over Rayleigh frequency selective fading channel, which is denoted as \( H_{k} \). Single-carrier frequency-division multiple access (SC-FDMA) transmission scheme is used for transmitting all the symbol components over \( H_{k} \). The 64-QAM symbols are converted to its corresponding bits [5]. The modulated symbol set is given as \( {\mathcal{S}} \, = \,\{ s_{0} ,\,s_{1} , \ldots ,\,s_{N - 1} \} \), where \( N \) is the total number of symbols, where \( s_{n} \in {\mathbb{C}} \) and \( {\mathbb{C}} \) denotes constellation. To each constellation point i.e. \( s_{n} \) is associated to a set of bits, the total number bits in constellation points is denoted as \( \mu \), \( \mu = { \log }_{2} \left( M \right) \). The bits in polar form for each symbol is given as \( {\mathcal{B}}\, = \,\{ b_{n}^{0} ,\,b_{n}^{1} ,\, \ldots ,\,b_{n}^{{\left( {\mu - 1} \right)}} \} \), with \( b_{n}^{\left( i \right)} \, = \, \pm 1\, = \,2\beta_{n}^{\left( i \right)} - 1,\,\beta_{n}^{\left( i \right)} = 0\,{\text{or}}\,1 \). The set of six bits can be decomposed from each 64-QAM constellation points and the bit subsets of the transmitted signal is given as \( {\mathcal{B}}_{m} \), \( m\, = \,0,\,1, \ldots ,\,M - 1 \), where \( M \) is the total number of constellation points. The bits of 64-QAM signal can be written as

$$ s_{n} \, = \,\sum\nolimits_{m = 0}^{M - 1} {g_{m} \prod\nolimits_{{b_{n}^{\left( i \right)} \in {\mathcal{B}}_{m} }} {b_{n}^{\left( i \right)} ;\,n = 0,\,1, \ldots ,\,M - 1 ;} } $$
(1)

and this forms a system of \( M \) equations for each constellation points in \( {\mathcal{S}} \) and \( g_{m} \) is the unknown variable. Without loss of generality, \( m \) is associated with its corresponding \( \mu \) bits and then \( m \) is given as, \( m\, = \,\left( {\gamma_{{\left( {\mu - 1,m} \right)}} ,\,\gamma_{{\left( {\mu - 2,m} \right)}} , \ldots ,\,\gamma_{{\left( {1,m} \right)}} ,\,\gamma_{{\left( {0,m} \right)}} } \right) \), This defines \( {\mathcal{B}}_{m} \) as the set of bits, where \( b_{n}^{\left( i \right)} \) is included if \( \gamma_{{\left( {i,m} \right)}} \) is 1.

Then, \( s_{n} \) can be written as

$$ s_{n} \, = \,\sum\nolimits_{m = 0}^{M - 1} {g_{m} \prod\nolimits_{i = 0}^{\mu - 1} {\left( {b_{n}^{\left( i \right)} } \right)^{{\gamma_{{\left( {i,m} \right)}} }} } } . $$
(2)

The transmitter has \( 3 \) layers, in the \( 1^{st} \) layer, the symbols is decomposed into \( 6 \) BPSK polar components. In the \( 2^{nd} \) layer, all the \( 6 \) polar components amplified by \( 6 \) parallel non-linear amplifiers and these amplifiers are connected to \( 6 \) parallel antennas. These \( 6 \) parallel antennas denoted as \( N_{m} \), where \( m\, = \,1,\,2,\,3,\,4,\,5,\,6 \). In the \( 3^{rd} \) layer, to reduce the complexity in channel estimation, \( N_{m} \) is combined and connected to only \( 1 \) transmitting antenna as in [11] and therefore this system model is similar to single input and multiple output antenna model. Due to amplification of low envelop BPSK components, the power efficiency improves with the decrease in amplifier distortion and lower PAPR of smaller components, it is possible to use non-linear amplifier in transmitter with negligible distortion [12]. All the six components are beamformed by using an angle \( \Theta \) to a desired direction towards the receiver like conventional beamforming mMIMO and here all the components are transmitted in a single communication channel. The total number of receiving antennas are denoted as \( R \).

Similar to [3, Fig. 1] here, six antennas with the BPSK components have a uniform distance \( d \) between them, this \( d \) can create time delay of the transmitted components at the receiver. The difference between each component at the receiver due to \( d \) and angle \( \theta \) can be compensated by using phase shift in the transmitter. The antennas \( R \) receives this signal at a time delay \( \Delta t \) between each antenna due to the angle \( \theta \) and this time delay is assumed to be uniform, which is due to multiple receiving antennas. \( \theta \) is originally meant for transmitting all the components in a uni-direction, however due to multiple receiving antenna, there is space \( d \) between each antenna is similar to transmitting antenna. In practice, \( \Delta t \) is common in all beamforming mMIMO systems and it can be compensated.

Fig. 1.
figure 1

Block diagram of IB-DFE receiver with joint channel and information estimate.

3 Joint Channel and Information Estimation

The joint channel estimation and signal detection method by using IB-DFE receiver as shown in Fig. 1, which is presented in this section. The received signal at \( R \) is given as

$$ Y_{k,l,i} \, = \,\left( {1 - \alpha } \right)\left( {H_{k,l} \left( {\sqrt {P_{x} } X_{k,l} + \sqrt {P_{q} } Q_{k,l} } \right) + N_{l} } \right)\, + \,N_{e,l} , $$
(3)

where \( X_{k,l} \) and \( Q_{k,l} \) are the data and pilot symbols and \( k \) is the frequency of block \( l \), then \( k\, = \,0,\,1, \ldots ,\,K - 1 \) and \( l\, = \,0,\,1, \ldots ,\,L - 1 \). \( X_{k,l} = {\text{DFT}}\left\{ {X_{n,l} } \right\} \) and \( Q_{k,l} = DFT\left\{ {q_{n,l} } \right\} \) are the information and pilot signals converted from time to frequency domain, respectively. \( N_{l} \) is additive Gaussian white noise (AWGN) due to signal transmission and \( N_{e,l} \) is due to power allocation at power splitting circuit. \( \left( {1 - \alpha } \right) \) is power allocated for information and \( \alpha \) is the power allocated for energy harvesting.

To perform joint channel and information estimation; the average channel estimate \( \tilde{H}_{k,l}^{avg} \) without using IB-DFE receiver is estimated over a set of signal block. By using received \( \tilde{H}_{k,l}^{avg} \) in IB-DFE, the information estimate \( X_{k,l}^{\left( i \right)} \) and new channel estimates \( \tilde{H}_{k,l}^{\left( i \right)} \) are estimated, where \( i \) is the number of IB-DFE iterations. In each iteration, the accuracy of estimate improves, and the final information estimate is denoted as \( X_{k,l,F}^{\left( i \right)} \).

The channel can be estimated by using the received signal and pilot signal, which is given as

$$ \begin{aligned} \tilde{H}_{k,l} & \, = \,\frac{{Y_{k,l,i} }}{{P_{q} Q_{k,l} }} \\ & = \,\left( {1 - \alpha } \right)\left( {\frac{{H_{k,l} \sqrt {P_{X} } X_{k,l} + N_{l} }}{{\sqrt {P_{q} } Q_{k,l} }}} \right) + \left( {1 - \alpha } \right)H_{k,l} + \frac{{N_{e,l} }}{{P_{q} Q_{k,l} }} = \left( {1 - \alpha } \right)H_{k,l} + eH_{k,l} ; \\ \end{aligned} $$
(4)
$$ eH_{k,l} = \left( {1 - \alpha } \right)\left( {\frac{{H_{k,l} \sqrt {P_{X} } X_{k,l} + N_{l} }}{{\sqrt {P_{q} } Q_{k,l} }}} \right) + \frac{{N_{e,l} }}{{P_{q} Q_{k,l} }}, $$
(5)

where \( eH_{k,l} \) is the channel estimation error. \( eH_{k,l} \) depends on the power of the information, noise and pilot signal. The frame structure has \( N \) sub carriers per block and it is modelled as in [7]. The expected value of \( H_{k,l} \), \( N_{l} \), \( N_{e} \), \( X_{k,l} \) and \( Q_{k,l} \), respectively is given as

$$ \begin{array}{*{20}l} {{\mathbb{E}}\left[ {H_{k,l} } \right] = 2\sigma_{h,k,l}^{2} ;\;{\mathbb{E}}\left[ {N_{l} } \right] = 2\sigma_{n,l }^{2} ;\;{\mathbb{E}}\left[ {N_{e,l} } \right] = 2\sigma_{n,e }^{2} ;} \hfill \\ {{\mathbb{E}}\left[ {X_{k,l} } \right] = N{\mathbb{E}}\left[ {X_{n,l}^{2} } \right] = 2\sigma_{X,k,l}^{2} ;} \hfill \\ {{\mathbb{E}}\left[ {Q_{k,l} } \right] = N{\mathbb{E}}\left[ {q_{n,l}^{2} } \right] = 2\sigma_{Q,k,l }^{2} ,} \hfill \\ \end{array} $$
(6)

then applying (6) in (5) gives

$$ {\mathbb{E}}\left[ {eH_{k,l} } \right] = \left( {1 - \alpha } \right)\left( {\frac{{\sigma_{X,k,l}^{2} + \sigma_{n,l}^{2} }}{{\sigma_{Q,k,l}^{2} }}} \right) + \frac{{\sigma_{n,e}^{2} }}{{\sigma_{Q,k,l}^{2} }}. $$
(7)

Even though channel fades with the change in frequency but remains constant for the respective frequency of a set of transmitted signals. Therefore, (7) can be is written as

$$ eH_{k,l}^{avg} = \frac{1}{l}\sum\nolimits_{0}^{L - 1} {\left( {{\mathbb{E}}\left[ {eH_{k,l} } \right]} \right),\;\;l = 0,1, \ldots ,L - 1;} $$
(8)
$$ \tilde{H}_{k,l}^{avg} = \left( {1 - \alpha } \right)\left( {H_{k,l} + eH_{k,l}^{avg} } \right), $$
(9)

Where \( eH_{k,l}^{avg} \) and \( \tilde{H}_{k,l}^{avg} \) are the average of channel estimation error and channel estimate over \( l \) blocks, respectively. The information estimated from the received signal by using channel estimate \( \tilde{H}_{k,l}^{avg} \) and \( Q_{k,l} \) is given as

$$ \begin{aligned} X_{k,l}^{hd} & = \frac{{Y_{k,l,i} }}{{\tilde{H}_{k,l}^{avg} }} - \left( {1 - \alpha } \right)\sqrt {P_{q} } Q_{k,l} \\ & = \frac{{\left( {1 - \alpha } \right)\left( {H_{k,l} \left( {\sqrt {P_{X} } X_{k,l} + \sqrt {P_{q} } Q_{k,l} } \right) + N_{l} } \right) + N_{e,l} }}{{\left( {1 - \alpha } \right)\left( {H_{k,l} + eH_{k,l}^{avg} } \right)}}\, - \,\left( {1 - \alpha } \right)\sqrt {P_{q} } Q_{k,l} . \\ \end{aligned} $$
(10)

The information estimate by using IB-DFE is written as

$$ X_{k,l}^{\left( i \right)} = \left( {Y_{k,l,i} - \left( {1 - \alpha } \right)\sqrt {P_{q} } Q_{k,l} \tilde{H}_{k,l}^{\left( i \right)} } \right)F_{k,l}^{\left( i \right)} - X_{k,l}^{{\left( {i - 1} \right)}} B_{k,l}^{\left( i \right)} , $$
(11)

where \( i \) is the number of iterations with \( i = 0,1, \ldots ,n \). In the first iteration, \( \tilde{H}_{k,l}^{\left( 0 \right)} = \tilde{H}_{k,l}^{avg} \), then, the feed forward coefficient \( F_{k,l}^{\left( i \right)} \) is given as

$$ F_{k,l}^{\left( i \right)} = \frac{{F_{k,l}^{\left( i \right)} }}{{\frac{1}{N}\sum\nolimits_{k = 0}^{N - 1} {\left( {F_{k,l}^{\left( i \right)} \tilde{H}_{k,l}^{\left( i \right)} } \right)} }}, $$
(12)

where \( F_{k,l}^{\left( i \right)} \) is given as

$$ F_{k,l}^{\left( i \right)} = \frac{{\tilde{H}_{k,l}^{\left( i \right)} }}{{\left( {\frac{{\sigma_{n,k,l}^{2} }}{{\sigma_{x,k,l}^{2} }}} \right) + |\tilde{H}_{k,l}^{\left( i \right)} |^{2} (1 - \left( {Cr^{{\left( {i - 1} \right)}} )^{2} } \right)}} $$
(13)

and the correlation factor \( Cr^{{\left( {i - 1} \right)}} = \frac{{{\mathbb{E}}\left[ {X_{n,l}^{\left( i \right)} X_{n,l}^{ *} } \right]}}{{{\mathbb{E}}\left[ {\left| {X_{n,l}^{2} } \right|} \right]}} \). The feedback co-efficient \( B_{k,l}^{\left( i \right)} \) is given as

$$ B_{k,l}^{\left( i \right)} = F_{k,l}^{\left( i \right)} \tilde{H}_{k,l}^{\left( i \right)} - 1 . $$
(14)

By using \( X_{k,l}^{\left( i \right)} \) and \( Q_{k,l} \) in (4), the improved channel estimates is given as

$$ \tilde{H}_{k,l}^{\left( i \right)} = \left( {1 - \alpha } \right)\left( {\frac{{H_{k,l} \left( {\sqrt {P_{X} } \left( {X_{k,l}^{\left( i \right)} } \right) + \sqrt {P_{q} } Q_{k,l} } \right) + N_{l} }}{{\sqrt {P_{X} } X_{k,l}^{\left( i \right)} + \sqrt {P_{q} } Q_{k,l} }}} \right) + \frac{{N_{e,l} }}{{\sqrt {P_{X} } X_{k,l}^{\left( i \right)} + \sqrt {P_{q} } Q_{k,l} }}. $$
(15)

Here, two channel estimate obtained from this receiver, they are:

  • The channel estimate i.e. (9), obtained without by using IB-DFE.

  • The channel estimate i.e. (15), obtained by using IB-DFE block. Here, estimate accuracy improves with each IB-DFE iteration until a saturation point.

Applying (15) instead of (9) in (11) gives information estimates than the previous estimates, which is denoted as \( X_{k,l,F}^{\left( i \right)} \). Here, three information estimates obtained from this receiver, they are:

  • The information estimate obtained by using \( \tilde{H}_{k,l}^{avg} \).

  • The information estimate i.e. (11), obtained by using \( \tilde{H}_{k,l}^{avg} \) and by using IB-DFE.

  • The information estimate i.e. \( X_{k,l,F}^{\left( i \right)} \), obtained by using improved channel estimate obtained with the help of IB-DFE i.e. (15) and by using IB-DFE.

4 Results

The performance of the system model is analyzed by using Monte Carlo simulations. Bit error rate (BER) analysis is considered as the performance metric for this model. For all the simulation results \( R = 32 \) and \( K = 256 \), while \( l = 2;3 \). The power ratio between pilot signal and the data symbols at the transmitter is denoted as \( \beta \). To analyze the impact of channel estimation, energy harvesting power allocation \( \alpha \) is considered zero except Fig. 4.

Figure 2 illustrates the BER performance of the system model with IB-DFE, where irrespective of number of iterations, the BER does not improves with iteration due to high diversity provided by \( R \). The results demonstrate that, under high order diversity, BER reaches its saturation point and there is no scope of improvement. Therefore, this simulation model is can be used in imperfect channel model condition, to understand the impact of IB-DFE receiver on the channel estimation technique.

Fig. 2.
figure 2

BER performance of the system model with perfect channel, where \( i \) is the IB-DFE receiver iteration and \( \alpha = 0 \).

Figure 3 illustrates the BER performance of the system model with IB-DFE receiver for an imperfect channel. Here, we compare the results based on the block size \( L = 2,3 \) and \( \beta = 1,3,5. \) The results shows that error rate performance improves with the increase in block size, while there is not much difference with \( \beta \) value, but results on the basis that \( \beta = 1 \), shows that BER performance is slightly better at \( L = 3 \), because the information estimate suffers due high power of pilot signal. This pattern shows the results of [8]. Thus, we consider \( L = 3 \) and \( \beta = 1 \) as an optimum configuration for this simulation.

Fig. 3.
figure 3

BER performance of the system model with IB-DFE receiver for an imperfect channel. Here, IB-DFE receiver iteration \( i = 4 \) for all the curves and \( \alpha = 0 \).

Figure 4 illustrates BER performance under imperfect channel condition and, the BER performance improves with increase IB-DFE iterations. The impact of power loss due to energy harvesting increases with increase in \( \alpha \) value. Thereby, IB-DFE incorporation with this channel estimate technique helps to achieve better performance as similar to system with perfect channel condition and helps in improving BER performance under energy harvesting condition. This simulation analysis explicitly demonstrates the advantage of using IB-DFE receiver for allowing rectenna to harvest \( 10\% \) of the total energy of the received signal with the loss of 1 dBm SNR.

Fig. 4.
figure 4

BER performance the system model with imperfect channel condition, where \( \alpha = \{ 0,\,0.1,\,0.2\} \), \( L = 3 \) and \( \beta = 1 \).

5 Conclusions

In this paper, we have analyzed the efficiency of IB-DFE receiver for joint channel and information estimation. The results show the improvement of error rate performance under the imperfect channel condition proves that IB-DFE receiver is perfect choice for employing SWIPT, which results in less one dBm SNR loss for using up to \( 10\% \) of power for energy harvesting. Also employing efficient power amplifier with less distortion due to signal low envelop makes this system model energy efficient at both transmitter and receiver. For the future, this work can be extended to multiple channel estimation for multi-user scenarios.