Spatial modulation based on reconfigurable antennas: performance evaluation by using the prototype of a reconfigurable antenna
Dung Nguyen Viet1,
Marco Di Renzo ORCID: orcid.org/0000-0003-0772-87931,
Vedaprabhu Basavarajappa2,
Beatriz Bedia Exposito2,
Jose Basterrechea3 &
Dinh-Thuy Phan-Huy4
EURASIP Journal on Wireless Communications and Networking volume 2019, Article number: 149 (2019) Cite this article
In this paper, we study the performance of spatial modulation based on reconfigurable antennas. Two main contributions are provided. We introduce an analytical framework to compute the error probability, which is shown to be accurate and useful for system optimization. We design and implement the prototype of a reconfigurable antenna that is specifically designed for application to spatial modulation and that provides multiple radiation patterns that are used to encode the information bits. By using the measured antenna radiation patterns, we show that spatial modulation based on reconfigurable antennas works in practice and that its performance can be optimized by appropriately selecting the radiation patterns to use for a given data rate.
Methods/experimental
The methods used in the paper are based on mathematical tools and theories. A new analytical framework for performance analysis is introduced. The theoretical framework is validated against Monte Carlo simulations, by using empirical data.
Spatial modulation (SM) [1] is a promising low-complexity [2] and energy-efficient [3] multiple-antenna modulation scheme, which is considered to be especially suitable for application to the Internet of Things (IoT) [4]. For these reasons, SM has attracted the attention of several academic and industrial researchers [5]. A comprehensive description of the main achievements and latest developments on SM research can be found in [6–10]. Some pioneering and recent experimental activities can be found in [11–14]. The research literature on SM is vast, and various issues have been tackled during the last few years, which include the analysis of the error probability [15], the design and analysis of transmit-diversity schemes [16, 17], and the analysis and optimization of the achievable rate [18, 19]. A recent comprehensive literature survey on SM and its generalizations is available in [10].
Among the many SM schemes that have been proposed in the literature, an implementation that is suitable for IoT applications is SM based on reconfigurable antennas (RectAnt-SM) [20]. In SM, the information bits are encoded onto the indices of the antenna elements of a given antenna array. In RectAnt-SM, by contrast, the information bits are encoded onto the radiation patterns (RPs) of a single-RF and reconfigurable antenna. This implementation has several advantages, especially for IoT applications [4].
In spite of the potential applications of RectAnt-SM in future wireless networks, to the best of the authors' knowledge, no analytical framework for computing the error probability of this emerging transmission technology is available. In the present paper, motivated by these considerations, we introduce an analytical framework that allows us to estimate the performance and to optimize the operation of RectAnt-SM. We prove, in particular, that the diversity order of RectAnt-SM is the same as the diversity order of SM, which coincides with the number of antennas at the receiver.
In order to substantiate the practical implementation and performance of RectAnt-SM, we design a single-RF and reconfigurable antenna that provides us with eight different RPs for encoding the information bits at a low-complexity and high energy efficiency. The proposed antenna is designed, and a prototype is implemented and measured in an anechoic chamber. Based on the manufactured prototype, we employ the measured RPs to evaluate the performance of RectAnt-SM. With the aid of our proposed analytical framework, in particular, we show that the error probability can be improved by appropriately choosing the best RPs, among those available, that minimize the average bit error probability.
Together with [4], the results contained in the present paper constitute the first validation of the performance of RectAnt-SM by using a realistic reconfigurable antenna that is capable of generating multiple RPs with adequate spatial characteristics for modulating information bits.
The remainder of the present paper is organized as follows. In Section 3, the system model is introduced. In Section 4, the analytical framework of the error probability is described. In Section 5, the prototype of the single-RF and reconfigurable antenna is presented. In Section 6, numerical results are illustrated, and the performance of RectAnt-SM is analyzed. Finally, Section 7 concludes the paper.
Notation: We adopt the following notation. Matrices, vectors, and scalars are denoted by boldface uppercase (e.g., A), boldface lowercase (e.g., a), and lowercase respectively (e.g., a). The element (u,v) of a matrix A is denoted by Au,v, and the uth entry of a vector a is denoted by au. The transpose, complex conjugate, and complex conjugate transpose of A are denoted by AT, A∗, and AH respectively. The absolute value of a complex number a is defined by \({\left | a \right |}\). \({\mathbb {E}_{a}}\left \{ \cdot \right \}\) and \({\mathbb {E}_{\mathbf {A}}}\left \{ \cdot \right \}\) denote the expectation operator of random variable a and matrix A, respectively. \(j = \sqrt { - 1} \) is the imaginary unit. \(\Pr \left \{ \cdot \right \}\) denotes probability. \({\cdot \choose \cdot }\) denotes the binomial coefficient. The Q-function is defined as \(Q(x) = \left ({1/\sqrt {2\pi } } \right)\int _{x}^{\infty } {\exp } \left ({ - \frac {{{u^{2}}}}{2}} \right){\mkern 1mu} du\). The moment generating function (MGF) of random variable X is defined as \({M_{X}}\left (s \right) = {\mathbb {E}_{X}}\left \{ {\exp \left ({ - sX} \right)} \right \}\). The gamma function is defined as \(\Gamma (z) = \int _{0}^{\infty } {{x^{z - 1}}} {e^{- x}}\,dx\). The modified Bessel function of order zero is defined as \({I_{0}}\left (\cdot \right)\).
In this section, we introduce the signal model, the channel model, and the demodulator.
Signal model
Let us consider a Nr×Nt multiple-input-multiple-output (MIMO) system that use a M-ary signal constellation diagram. In a conventional SM transmission scheme [1], the data stream is divided into two blocks, where the first block of \(\log _{2}(N_{t})\) bits is used to identify the index of the transmitted antenna and \(\log _{2}(M)\) bits are used to identify a symbol of the signal constellation diagram. By assuming that channel state information (CSI) is known at the receiver, the objective of the detector is to jointly estimate the active antenna and the data symbol that is transmitted in order to retrieve the entire transmitted bitstream.
Let us now consider the RecAnt-SM transmission scheme. In this case, we consider that the transmitter is equipped with a reconfigurable antenna that is capable of generating P different RPs. In this case, the RP that is used for transmission is also used to encode the information, in addition to the signal constellation diagram [20]. More precisely, the joint combination of antenna's RP and channel constitutes the physical resource that is used to encode the information bits. If P RPs are available, then \(\log _{2}(P)+\log _{2}(M)\) bits of information can be transmitted.
Compared with conventional SM, RecAnt-SM has the advantage of not requiring an array of antennas to enhance the data rate. In addition, RecAnt-SM can be implemented at a low cost, by using simple and compact antennas, where different RPs can be obtained by realizing appropriate circuits that modify the current flowing through the physical antenna. These specific features of RectAnt-SM make it useful for IoT applications. In the sequel, we will discuss an antenna that we have designed and fabricated and that allows us to implement RectAnt-SM in practice. As an example, let us consider P=M=4. Then, \(\log _{2}(P) = 2\) bits are used to identify the RP and \(\log _{2}(M) = 2\) are used to identify a symbol of the signal constellation diagram. Figure 1 illustrates a simple example of the implementation of RectAnt-SM, by focusing only on the RPs, e.g., RectAnt-SSK (space shift keying).
a–d Illustration of RecAnt-SSK with P=4. The RPs are obtained from a designed and manufactured antenna that is described in the sequel
Let us assume that the pth RP and the symbol xm are selected based on information bits to be transmitted. The Nr×1 received signal vector can be formulated as follows:
$$\begin{array}{*{20}l} \mathbf{y} &= \sqrt{\rho}\mathbf{He}_{p} x_{m} + \mathbf{w} \end{array} $$
where ρ is the average signal-to-noise-ratio (SNR) at each receive antenna; w is the Nr×1 complex additive white Gaussian noise vector of zero mean and unit variance; ep for p=1,…,P is a P×1 vector whose pth entry is equal to one and the other entries are equal to zero; and H is the Nr×P channel matrix that accounts for the antenna RPs as well. The vector ep allows one to select the specific RP given the bits to be transmitted. The channel matrix H is introduced in the next section.
RecAnt-SSK constitutes a special case of RecAnt-SM, where the information bits are encoded only into the RPs. In this case, the signal model simplifies as follows:
$$\begin{array}{*{20}l} \mathbf{y} = \sqrt{\rho}\mathbf{H}\mathbf{e}_{p} + \mathbf{n} \end{array} $$
Based on the signal model in (1), the maximum likelihood (ML-) optimum demodulator, by assuming full CSI available at the receiver, can be formulated as follows [2]:
$$\begin{array}{*{20}l} \left({\overset{\smallfrown}{q}},{x_{\overset{\smallfrown}{n}}} \right) = \underset{{\text{for}}\:q = 1, \ldots,P\:{\text{and}}\:n = 1, \ldots,M}{\arg \max} \left\{{D\left({q,{x_{n}}}\right)}\right\} \end{array} $$
$$\begin{array}{*{20}l} D\left({q,{x_{n}}} \right) = \sum\limits_{{n_{r}} = 1}^{{N_{r}}} {\left[ {\mathbf{y}_{{n_{r}}}^{*}\left(\mathbf{H}_{{n_{r}},q}{x_{n}}\right) - \frac{1}{2}{{\left| {{\mathbf{H}_{{n_{r}},q}}{x_{n}}} \right|}^{2}}} \right]} \end{array} $$
where \(\mathbf {y}_{{n_{r}}}\) is the nrth entry of y and \(\mathbf {H}_{{n_{r}},q}\) the entry in the row nr and column q of H. The demodulator of RectAnt-SSK can be obtained in a similar way.
Channel model
In this section, we introduce the channel model H that we briefly mentioned in the previous section. As far as RectAnt-SM is concerned, the channel model plays an important role, since the combined effect of RP and channel determines the system performance. Since we consider three-dimensional RPs (see Fig. 1), the considered channel model is chosen appropriately. More precisely, the channel model is based on a ray-based and cluster approach similar to [21]. A sketched representation of the channel model is given in Fig. 2.
Sketched representation of the channel model
For simplicity and without loss of generality, we consider a channel model with a single cluster and with multiple rays. The number of rays is denoted by K. The Nr×P channel matrix, which accounts for the RPs of the reconfigurable antenna as well, can be formulated as follows:
$$\begin{array}{*{20}l} {\mathbf{H}} = \frac{1}{{\sqrt K }}\sum\limits_{k = 1}^{K} {{\beta_{k}}} {\mathbf{a}^{r}}\left({{\theta_{k}^{r}},{\phi_{k}^{r}}} \right){\left({{\mathbf{a}^{t}}\left({{\theta_{k}^{t}},{\phi_{k}^{t}}} \right)} \right)^{T}} \end{array} $$
where the normalization factor \(1/{\sqrt K }\) preserves the average unit energy of the channel, and the following notation is used:
βk is the fading coefficient of the kth ray;
ar is the Nr×1 array response vector of the receiver;
at is the Nt×1 array response vector of the transmitter;
\(\left ({{\theta _{k}^{t}},{\phi _{k}^{t}}} \right)\) are the azimuth and elevation angles of departure (AoD) of the kth ray; and
\(\left ({{\theta _{k}^{r}},{\phi _{k}^{r}}} \right)\) are the azimuth and elevation angles of arrival (AoA) of the kth ray.
As far as the statistical distributions of βk, θk, and ϕk as concerned, Table 1 summarizes the most commonly used models. In particular, θk is often modeled as a truncated Laplacian random variable, and ϕk is often modeled as a Von-Mises, truncated Gaussian, or an uniform random variable. For each ray, the random variables are assumed to be independent and identically distributed.
Table 1 Distribution of variables for the considered channel model
The specific characteristics of the receiver and transmitter are determined by ar and at, respectively. As far as the receiver is concerned, we consider, as an example, a uniform linear array with an onmi-directional RP of unit gain. As far as the transmitter is concerned, we consider a reconfigurable antenna with P different RPs. The RPs are denoted as follows:
$$\begin{array}{*{20}l} {{\mathcal{G}}_{p}}\left({\theta,\phi} \right) = \sqrt {{G_{p}}\left({\theta,\phi} \right)} \exp \left({j{\Omega_{p} }\left({\theta,\phi} \right)} \right) \end{array} $$
for p=1,…,P, and Gp(θ,ϕ) and Ωp(θ,ϕ) are the amplitude and phase of the pth RP, respectively.
Based on these assumptions, \(\mathbf {H}_{{n_{r}},p}\) can be formulated as follows:
$$\begin{array}{*{20}l} {} {\mathbf{H}}_{n_{r},p} = \frac{1}{{\sqrt{K}}} {\sum\limits_{k = 1}^{K} {{\beta_{k}}}} &\underbrace{\exp (j\parallel{\mathbf{k}}\parallel {d}(n_{r} - 1)\sin ({\theta_{k}^{r}})\sin ({\phi_{k}^{r}}))}_{{\text{Receiver\,part}}} \times \\ &\underbrace {\sqrt {{G_{p}}({\theta_{k}^{t}},{\phi_{k}^{t}})} \exp (j{\Omega_{p}}({\theta_{k}^{t}},{\phi_{k}^{t}}))}_{{\text{Transmitter\,part}}} \end{array} $$
where nr=1,…,Nr, p=1,…,P, d is the distance between adjacent antennas, and k(θ,ϕ) is the wavevector defined as follows:
$$ {}{\mathbf{k}}(\theta,\phi) = \frac{{2\pi }}{\lambda }{\left[ {\sin (\theta)\cos (\phi),\sin (\theta)\sin (\phi),\cos (\theta)} \right]^{T}} $$
where λ is the wavelength.
Average bit error probability
From [15], it is known that the average bit error probability (ABEP) of the system model under analysis can be formulated as follows:
$$ {\begin{aligned} \text{ABEP} &\le \frac{1}{{PM}}\frac{1}{{{{\log }_{2}}\left({PM} \right)}}\sum\limits_{p = 1}^{P} \sum\limits_{m = 1}^{M} \sum\limits_{q = 1}^{P} \sum\limits_{n = 1}^{M}\\&\quad {{N_{H}}\text{APEP}\left({\left({p,{x_{m}}} \right) \to \left({q,{x_{n}}} \right)} \right)} \end{aligned}} $$
where APEP denotes the average pairwise error probability, which is the probability of demodulating the RP q and the symbol xn if the RP p and the symbol xm have been transmitted and are the only two possible options, and NH denotes the number of bits that the latter two constellation points different from each other.
The pairwise error probability (PEP), \({\Pr }_{E}\), of deciding for the qth radiation pattern while the pth radiation pattern is transmitted and deciding for the symbol xn while the symbol xm is transmitted can be written as follows:
$$\begin{array}{*{20}l} {\Pr}_{E}\left({\left({p,{x_{m}}} \right) \to \left({q,{x_{n}}} \right)\left| \mathbf{H} \right.} \right) &= \Pr \left({D\left({p,{x_{m}}} \right) < D\left({q,{x_{n}}} \right)} \right) \notag\\ &=Q\left({\sqrt {\frac{\rho }{2}{\gamma_{p,q,{x_{m}},{x_{n}}}}\left(\mathbf{H} \right)}} \right) \end{array} $$
$$\begin{array}{*{20}l} {\gamma_{p,q,{x_{m}},{x_{n}}}\left(\mathbf{H} \right)} &= \sum\limits_{{n_{r}} = 1}^{{N_{r}}} {{{\left| {{\mathbf{H}_{{n_{r}},q}}{x_{n}} - {\mathbf{H}_{{n_{r}},p}}{x_{m}}} \right|}^{2}}} \end{array} $$
As a result, the APEP can be written as follows:
$$\begin{array}{*{20}l} \text{APEP}&\left(\left({p,{x_{m}}} \right) \to \left({q,{x_{n}}} \right) \right) \notag\\ &= {\mathbb{E}_{\mathbf{H}}}\left\{ {Q\left({\sqrt {\frac{\rho }{2}{\gamma_{p,q,{x_{m}},{x_{n}}}}\left(\mathbf{H} \right)}} \right)} \right\} \notag\\ & {\stackrel{(a)}{=}}{\mathbb{E}_{\mathbf{H}}}\left\{ {\frac{1}{\pi }\int_{0}^{\pi /2} {\exp} \left({\frac{{ - \rho{\gamma_{p,q,{x_{m}},{x_{n}}}} }}{{4{{\sin }^{2}}(\vartheta)}}} \right)d\vartheta} \right\}\notag\\ &{\stackrel{(b)}{=}}\frac{1}{\pi }\int_{0}^{\pi /2} { {{M_{_{{\gamma_{p,q,{x_{m}},{x_{n}}}}}}}\left({\frac{{ - \rho }}{{4{{\sin }^{2}}(\vartheta)}}} \right)}} d\vartheta \end{array} $$
where (a) and (b) follow by applying Craig's formula [15] and from the definition of MGF of \({\gamma _{p,q,{x_{m}},{x_{n}}}}\), respectively.
Setup with N r=1
We start by considering the system setup with Nr=1. The APEP is formulated in the following proposition.
Let \(\bar \zeta \left ({K,\vartheta }\right) = \frac {\rho }{{4K{{\sin }^{2}}(\vartheta)}}\). If Nr=1, the APEP of RecAnt-SM is as follows:
$$ {\begin{aligned} &{\text{APEP}}=\\ & \frac{1}{\pi }\int_{0}^{\pi /2} {\int_{0}^{\infty} {\exp}} (- z)\left(\int_{0}^{\pi} {\int_{- \pi }^{\pi} {\exp} } \left(- z\bar \zeta \left({K,\vartheta} \right)\right.\right.\\& \left. \left.\psi\left(\ {p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \right) \right){f_{\theta} }\left({\theta^{t}}\right){f_{\phi} }\left({\phi^{t}}\right)\,d{\theta^{t}}\,d{\phi^{t}} \vphantom{{\int_{0}^{\pi}}}\!\right)^{K}\!dz\,d\vartheta \end{aligned}} $$
$$ {{} \begin{aligned} &\!\psi \left({p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \right) \,=\, \!\left| \sqrt {{G_{q}}\left({\theta^{t}},{\phi^{t}}\right)} \!\exp \left(j{\Omega_{q}}\left({\theta^{t}},{\phi^{t}}\right)\right){x_{n}} \right.\\& \left. - \sqrt {{G_{p}}({\theta^{t}},{\phi^{t}})} \exp (j{\Omega_{p}}({\theta^{t}},{\phi^{t}})){x_{m}} \right|^{2} \end{aligned}} $$
See the Appendix. □
In the high-SNR regime, the APEP is given in the following proposition.
If Nr=1, the APEP of RecAnt-SM in the high-SNR regime is as follows:
$$ {{} \begin{aligned} {\text{APEP}} &\!\le\! \frac{1}{\rho}\int_{0}^{\infty} \!\left(\int_{0}^{\pi} {\int_{- \pi }^{\pi} {\exp} } \left(\,-\, z\frac{1}{K}\psi \left({p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \!\right) \right) \right.\\ &\left. {f_{\theta}}(\theta^{t}){f_{\phi} }(\phi^{t})\,d\theta^{t} \,d\phi^{t} \right)^{K}\,dz \end{aligned}} $$
where ψ is defined in (14).
Equations (13) and (15) provide one with accurate performance predictions of the APEP. The analytical expression is, however, quite complex due to the discrete number of rays that are considered in the system model. In the following two propositions, we provide an asymptotic expression of the APEP under the assumption \(K \to \infty \). In the sequel, we will show that the simplified analytical expression of the APEP is accurate even for moderate values of K.
Assume \(K \to \infty \) and Nr=1. The APEP can be simplified as follows:
$$\begin{array}{*{20}l} &{\text{APEP}} \leq \frac{1}{2}\left({1 - \sqrt {\frac{{\rho \Theta }}{{\rho \Theta + 4}}}} \right) \end{array} $$
$$\begin{array}{*{20}l} \Theta = \int_{0}^{\pi} {\int_{- \pi }^{\pi} {\psi \left({p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \right)}} {f_{\theta} }({\theta^{t}}){f_{\phi} }({\phi^{t}})d{\theta^{t}}d{\phi^{t}} \end{array} $$
Assume \(K \to \infty \) and Nr=1. The APEP in the high-SNR regime can be simplified as follows:
$$\begin{array}{*{20}l} {\text{APEP}} \le \frac{1}{\rho\Theta } \end{array} $$
where Θ is defined in (17).
By direct inspection of the obtained expression of the APEP, we evince that the error probability is minimized as Θ increases. As expected, therefore, we evince that the error probability decreases as the difference between the RPs increases. It is worth mentioning, however, that the difference between the RPs is weighted by the distribution of the AoD of the rays. Thus, the APEP depends on both the RPs themselves and the specific characteristics of the channel.
In this section, we study the system setup where two antennas are available at the receiver. The following proposition generalizes Proposition 3
Assume \(K \to \infty \) and Nr=2. The APEP in the high-SNR regime can be formulated as follows:
$$ {{}\begin{aligned} &{\text{APEP}} \le \\ & \frac{3}{{{\rho^{2}}{{\left({\int_{- \pi }^{\pi} {\int_{0}^{\pi}\! {\psi \left({p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \right)}} {f_{\theta} }({\theta^{t}}){f_{\phi} }({\phi^{t}})d{\theta^{t}}d{\phi^{t}}} \right)}^{2}}\Psi }}\\ &\begin{array}{cc} \Psi = 1 - \left(\int_{- \pi }^{\pi} \int_{0}^{\pi} \cos \left(\left\| \mathbf{k} \right\|{d}\left({\sin \left({{\theta^{r}}} \right)\sin \left({{\phi^{r}}} \right)} \right) \right)\right. \\ \left.{f_{\theta} }({\theta^{r}}){f_{\phi} }({\phi^{r}})d{\theta^{r}}d{\phi^{r}} \right)^{2} - \\ \qquad \:\:\:\:\left(\int_{- \pi }^{\pi} \int_{0}^{\pi} \sin \left({\left\| \mathbf{k} \right\|{d}\left({\sin \left({{\theta^{r}}} \right)\sin \left({{\phi^{r}}} \right)} \right)} \right)\right.\\ \left.{f_{\theta} }({\theta^{r}}){f_{\phi} }({\phi^{r}})d{\theta^{r}}d{\phi^{r}} \right)^{2} \end{array} \end{aligned}} $$
By direct inspection of the obtained APEP, we evince that the diversity order is equal to two if two antennas are available at the receiver. This is consistent with conventional SM [15].
Setup with generic N r
In this section, we generalize the previous analytical frameworks for Nr=1 and Nr=2, by considering a generic value of Nr.
Assume \(K \to \infty \). The APEP in the high-SNR regime can be formulated as follows:
$$ {{}\begin{aligned} &{\text{APEP}} \le \\ & \frac{{{\alpha_{{N_{r}}}}}}{{{\rho^{{N_{r}}}}{{\left[ {\psi \left({p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \right)} \right]}^{{N_{r}}}}{\mathbb{E}_{{\theta^{r}},{\phi^{r}}}}\left\{ {F\left({{\theta^{1,r}}, \ldots,{\theta^{{N_{r}},r}},{\phi^{1,r}}, \ldots,{\phi^{{N_{r}},r}}} \right)} \right\}}}\\ &{\alpha_{{N_{r}}}} = \frac{1}{2}\left({\begin{array}{*{20}{c}} {2{N_{r}}}\\ {{N_{r}}} \end{array}} \right) + \sum\limits_{k = 0}^{{N_{r}} - 1} {{{\left({ - 1} \right)}^{{N_{r}} - k}}2\left({\begin{array}{*{20}{c}} {2{N_{r}}}\\ k \end{array}} \right)\frac{{\sin \left({\pi \left({{N_{r}} - k} \right)} \right)}}{{2\pi \left({{N_{r}} - k} \right)}}} \end{aligned}} $$
By direct inspection of the obtained APEP, we observe that the diversity order is equal to Nr, as in conventional SM [15].
The main limitation of (20) is that \(\mathbb {E}_{{\theta ^{r}},{\phi ^{r}}}\left \{ {F\left ({{\theta ^{1,r}}, \ldots,{\theta ^{{N_{r}},r}},{\phi ^{1,r}}, \ldots,{\phi ^{{N_{r}},r}}} \right)} \right \}\) is not analytically tractable and cannot, in general, be formulated in closed-form. In some special cases, however, this is possible, notably if Nr=3, as reported in the following proposition.
$$ {\begin{aligned} &{\text{APEP}} \le\\ & \frac{{{10}}}{{{\rho^{{3}}}{{\left[ {\psi \left({p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \right)} \right]}^{{3}}}\mathbb{E}_{{\theta^{r}},{\phi^{r}}}\left\{ {F\left({{\theta^{1,r}}, {\theta^{2,r}}, {\theta^{{3},r}},{\phi^{1,r}}, {\phi^{2,r}},{\phi^{{3},r}}} \right)} \right\}}} \end{aligned}} $$
where the following definition holds true:
$$ {\begin{aligned} &\mathbb{E}_{{\theta^{r}},{\phi^{r}}}\left\{ {F\left({{\theta^{1,r}}, {\theta^{2,r}}, {\theta^{{3},r}},{\phi^{1,r}}, {\phi^{2,r}},{\phi^{{3},r}}} \right)} \right\} = \\ &1 + 2{\left({{E_{1}}} \right)^{2}}{E_{3}} - 2{\left({{E_{2}}} \right)^{2}}{E_{3}} + 4{E_{1}}{E_{2}}{E_{4}} - {\left({{E_{3}}} \right)^{2}} \\&\quad- {\left({{E_{4}}} \right)^{2}} - 2{\left({{E_{1}}} \right)^{2}} - 2{\left({{E_{2}}} \right)^{2}} \end{aligned}} $$
$$\begin{array}{*{20}l} {E_{1}} &= \int_{- \pi }^{\pi} \int_{0}^{\pi} {\cos \left({\left\| \mathbf{k} \right\|{d}\left({\sin \left({{\theta^{r}}} \right)\sin \left({{\phi^{r}}} \right)} \right)} \right)}\\ &\quad {f_{\theta} }\left({{\theta^{r}}} \right){f_{\phi} }\left({{\phi^{r}}} \right)d{\theta^{r}}d{\phi^{r}} \end{array} $$
$$\begin{array}{*{20}l} {E_{2}} &= \int_{- \pi }^{\pi} \int_{0}^{\pi} \sin \left({\left\| \mathbf{k} \right\|{d}\left({\sin \left({{\theta^{r}}} \right)\sin \left({{\phi^{r}}} \right)} \right)} \right)\\ &\quad{f_{\theta} }\left({{\theta^{r}}} \right){f_{\phi} }\left({{\phi^{r}}} \right)d{\theta^{r}}d{\phi^{r}} \end{array} $$
$$\begin{array}{*{20}l} {E_{3}} &= \int_{- \pi }^{\pi} \int_{0}^{\pi} {\cos \left({2\left\| \mathbf{k} \right\|{d}\left({\sin \left({{\theta^{r}}} \right)\sin \left({{\phi^{r}}} \right)} \right)} \right)}\\ &\quad {f_{\theta} }\left({{\theta^{r}}} \right){f_{\phi} }\left({{\phi^{r}}} \right)d{\theta^{r}}d{\phi^{r}} \end{array} $$
$$\begin{array}{*{20}l} {E_{4}} &= \int_{- \pi }^{\pi} \int_{0}^{\pi} \sin \left({2\left\| \mathbf{k} \right\|{d}\left({\sin \left({{\theta^{r}}} \right)\sin \left({{\phi^{r}}} \right)} \right)} \right)\\ &\quad {f_{\theta} }\left({{\theta^{r}}} \right){f_{\phi} }\left({{\phi^{r}}} \right)d{\theta^{r}}d{\phi^{r}} \end{array} $$
Antenna prototype with reconfigurable radiation patterns
In order to test the performance of RectAnt-SM and to assess its practical feasibility, we have designed and fabricated a reconfigurable antenna that yields multiple radiation patterns. A photo of the antenna prototype specifically designed to implement RectAnt-SM is given in Fig. 3. The proposed antenna needs a single radio frequency chain and is capable of generating eight RPs. Four of them are illustrated in Fig. 1. The eight different RPs are obtained by considering a 4×4 array, as illustrated in Fig. 4, and by using a different excitation matrix. Five excitation matrices are reported as follows:
$$ {\begin{aligned} & {\stackrel{\left[ {\begin{array}{*{20}{c}} 0& + & + &0\\ 0& + & + &0\\ 0& - & - &0\\ 0& - & - &0 \end{array}} \right],\:}{\text{A}}} {\stackrel{\left[ {\begin{array}{*{20}{c}} 0&0& + & + \\ 0&0& + & + \\ - & - &0&0\\ - & - &0&0 \end{array}} \right],\:}{\text{B}}} {\stackrel{\left[ {\begin{array}{*{20}{c}} 0&0&0&0\\ - & - & + & + \\ - & - & + & + \\ 0&0&0&0 \end{array}} \right],}{\text{C}}} \\&\quad {\stackrel{\left[ {\begin{array}{*{20}{c}} - & - &0&0\\ - & - &0&0\\ 0&0& + & + \\ 0&0& + & + \end{array}} \right],\:}{\text{D}}} {\stackrel{\left[ {\begin{array}{*{20}{c}} 0& - & - &0\\ 0& - & - &0\\ 0& + & + &0\\ 0& + & + &0 \end{array}} \right]}{\text{E}}} \end{aligned}} $$
Fabricated antenna for implementing RectAnt-SM
Phasing network of the designed antenna
Each entry of the 4×4 matrix represents the excitation that is fed into the corresponding antenna of the array. The excitation has unit amplitude, and its phase is either 0 or 180, which is denoted by the sign "+" and " −", respectively, in the excitation matrices.
The designed prototype has been studied and optimized by simulating the 4×4 antenna array with the full-wave solver Ansys HFSS. Based on the optimized design, a prototype antenna has been fabricated and its RPs have been measured in an anechoic chamber, as illustrated in Fig. 5. In Fig. 6, we report, as an example, the simulated and measured RP RP1 that is reported in Fig. 1 as well.
Measurement of the antenna prototype in an anechoic chamber
Comparison of simulated (Sim) and measured (Meas) RP RP1. The beam cuts are measured at differently oriented angles (0 and 90 ∘) in space. "Co" stands for "Copolar gain," and "XP" stands for "Cross Polar gain"
Numerical results and discussion
In this section, we provide some numerical results in order to validate the analytical derivation of the bit error probability and in order to asses the achievable performance by using the RPs that are obtained from the fabricated antenna prototype. We show, in particular, that by appropriately selecting the RPs that minimize the analytical framework of the error probability, the error probability can be greatly decreased.
In Fig. 7, we compare Monte Carlo simulations against the proposed analytical framework of the bit error probability. In Fig. 8, we study the impact of K on the error probability with the aim of assessing the accuracy of the asymptotic framework for \(K \to \infty \). Both figures confirm that our analytical frameworks are accurate and in agreement with Monte Carlo simulations.
ABEP of RecAnt-SM (P=2, Nr=1, M=2 (BPSK)). The results are obtained by using RP1 and RP2
ABEP of RecAnt-SM (P=2, Nr=1, M=2 (BPSK)), as a function of K. The results are obtained by using RP1 and RP2
In Figs. 9, 10, and 11, we report the bit error probability as a function of the number of antennas at the receiver. We observe the good accuracy of the proposed analytical frameworks, and we note, in particular, that RectAnt-SM provides one with a diversity order equal to the number of antennas at the receiver.
ABEP of RecAnt-SM (P=2, Nr=1,2,3, M=2 (BPSK)). The results are obtained by using RP1 and RP2
ABEP of RecAnt-SM (P=4, Nr=1,2,3, M=2 (BPSK)). The results are obtained by using RP1- RP4
ABEP of RecAnt-SM (P=4, Nr=1,2,3, M=4 (QPSK)). The results are obtained by using RP1- RP4
In Fig. 12, we report the APEP as a function of the RPs. In particular, with the aid of the analytical framework of the error probability, we compute the ABEP for all possible combinations of four out of eight RPs in order to identify the impact of the RPs on the error performance. We observe that an appropriate choice of the RPs can yield a significant performance gain. The reason is that, due to practical design constraints, it may not be possible to design several RPs that are very different from each other. Therefore, identifying the best of them that provide good performance as a function of the channel model is an important optimization problem. The proposed analytical framework allows us to solve this optimization problem at a low complexity and high efficiency.
ABEP of RecAnt-SSK (P=4, Nr=3). The different markers show the ABEP by choosing four RPs out of eight RPs in order to identify those that offer the best performance
In this paper, we have studied the performance of spatial modulation based on reconfigurable antennas. We have introduced an analytical framework to compute the error probability and have described the prototype of a reconfigurable antenna that is specifically designed for application to spatial modulation. By using the radiation patterns obtained from the manufactured antenna prototype, we have shown that spatial modulation based on reconfigurable antennas works in practice and that its performance can be optimized by appropriately choosing the radiation patterns that minimize the proposed analytical framework of the error probability.
Let \({\nu _{{n_{r}}}} = {\mathbf {H}_{{n_{r}},q}}{x_{n}} - {\mathbf {H}_{{n_{r}},p}}{x_{m}}\) for nr=1,…,Nr. Recall that Nr=1. We note that ν1 is a zero mean complex Gaussian variable with variance:
$$ {\begin{aligned} \sigma_{{\nu_{1}}}^{2} &= \frac{1}{K}\sum\limits_{k = 1}^{K} \left| \sqrt {{G_{q}}\left({\theta_{k}^{t}},{\phi_{k}^{t}}\right)} \exp \left(j{\Omega_{q}}\left({\theta_{k}^{t}},{\phi_{k}^{t}}\right)\right){x_{n}}\right.\\& \left.- \sqrt {{G_{p}}\left({\theta_{k}^{t}},{\phi_{k}^{t}}\right)} \exp \left(j{\Omega_{p}}\left({\theta_{k}^{t}},{\phi_{k}^{t}}\right)\right){x_{m}} \right|^{2} \end{aligned}} $$
Thus, \({\gamma _{p,q,{x_{m}},{x_{n}}}} = \mid {\nu _{1}}{\mid ^{2}}\) is an exponential random variable whose probability density function is:
$$\begin{array}{*{20}l} f_{x}(x;\lambda) = \lambda\exp(-\lambda x), ~x\geqslant 0 \end{array} $$
with \(\lambda = 1/\sigma _{{\nu _{1}}}^{2}\).
The APEP can then be formulated as follows:
$$\begin{array}{*{20}l} \text{APEP} = \frac{1}{\pi }\int_{0}^{\pi /2} {{\mathbb{E}_{\sigma_{{\nu_{1}}}^{2}}}} \left\{ {{{\left({1 + \frac{{\rho \sigma_{{\nu_{1}}}^{2}}}{{4{{\sin }^{2}}(\vartheta)}}} \right)}^{- 1}}} \right\}{\mkern 1mu} d\vartheta \end{array} $$
We note that \(\sigma _{{\nu _{1}}}^{2}\) is a random variable that depends on \(\left \{ {{\theta _{k}^{t}},{\phi _{k}^{t}}} \right \}_{k = 1}^{K}\). The expectation can be computed by using the approach introduced in [24], as follows:
$$\begin{array}{*{20}l} \text{APEP} = \frac{1}{\pi }\int_{0}^{\pi /2} {\left({\int_{0}^{\infty} {{M_{S}}} (z){\mkern 1mu} dz} \right)} {\mkern 1mu} d\vartheta \end{array} $$
where \(S = 1 + \frac {{\rho \sigma _{{\nu _{1}}}^{2}}}{{4{{\sin }^{2}}(\vartheta)}}\), and:
$$\begin{array}{*{20}l} {M_{S}}(z) = \exp (- z){M_{\sigma_{{\nu_{1}}}^{2}}}\left({z\frac{\rho }{{4{{\sin }^{2}}(\vartheta)}}} \right) \end{array} $$
With the aid of some algebraic manipulations, we obtain the following:
$$ {{}\begin{aligned} &{M_{\sigma_{{\nu_{1}}}^{2}}}\left({z\frac{\rho }{{4{{\sin }^{2}}(\vartheta)}}} \right)\\ &\,=\, {\mathbb{E}_{{\theta_{k}^{t}},{\phi_{k}^{t}}}}\left\{ {\exp \left({ \!- z\frac{\rho }{{4K{{\sin }^{2}}(\vartheta)}}\sum\limits_{k = 1}^{K} {\psi \left({p,q,{x_{m}},{x_{n}},{\theta_{k}^{t}},{\phi_{k}^{t}}} \right)}} \right)} \right\} \end{aligned}} $$
$$ {{}\begin{aligned} &={\mathbb{E}_{{\theta_{k}^{t}},{\phi_{k}^{t}}}}\left\{ {\prod\limits_{k = 1}^{K} {\exp \left({ \,-\, z\frac{\rho }{{4K{{\sin }^{2}}(\vartheta)}}\psi \left({p,q,{x_{m}},{x_{n}},{\theta_{k}^{t}},{\phi_{k}^{t}}} \right)} \right)}} \right\}\\ &\,=\,\prod\limits_{k = 1}^{K} \left(\int_{- \pi }^{\pi} \int_{0}^{\pi} \exp \left(\!- z\frac{\rho }{{4K{{\sin }^{2}}(\vartheta)}}\psi \left(p,q,{x_{m}},{x_{n}},{\theta_{k}^{t}},{\phi_{k}^{t}} \right) \right)\right.\\ & \left. \quad {f_{\theta} }({\theta_{k}^{t}}){f_{\phi} }({\phi_{k}^{t}})d{\theta_{k}^{t}}d{\phi_{k}^{t}} {\vphantom{\int_{- \pi }^{\pi}}}\right) \end{aligned}} $$
$${{} \begin{aligned} &\psi \left({p,q,{x_{m}},{x_{n}},{\theta_{k}^{t}},{\phi_{k}^{t}}} \right) = \left| \sqrt {{G_{q}}({\theta_{k}^{t}},{\phi_{k}^{t}})} \exp \left(j{\Omega_{q}}\left({\theta_{k}^{t}},{\phi_{k}^{t}}\right)\right){x_{n}} \right. \\ & \left.- \sqrt {{G_{p}}({\theta_{k}^{t}},{\phi_{k}^{t}})} \exp (j{\Omega_{p}}({\theta_{k}^{t}},{\phi_{k}^{t}})){x_{m}} \right|^{2} \end{aligned}} $$
The proof follows from the following identity:
$$ {\begin{aligned} &\int_{- \pi }^{\pi} \int_{0}^{\pi} {f_{\theta} }\left({\theta^{t}}\right){f_{\phi}}({\phi^{t}})d{\theta_{k}^{t}}d{\phi_{k}^{t}} \\& \quad= \int_{- \pi }^{\pi} {f_{\theta} }\left({\theta^{t}}\right)d{\theta_{k}^{t}}\int_{0}^{\pi} {f_{\phi} }\left({\phi^{t}}\right)d{\phi_{k}^{t}} = 1 \end{aligned}} $$
From (29), we obtain, in the high-SNR regime, the following:
$${\begin{aligned} \text{APEP}_{p,q} &= \mathbb{E}_{\sigma^{2}_{\nu_{1}}}\left\{\frac{1}{\pi} \int_{0}^{\pi/2}\left(1+\frac{\rho\sigma^{2}_{\nu_{1}}}{4\sin^{2}(\vartheta)}\right)^{-1}\,d\vartheta \right\} \end{aligned}} $$
$${\begin{aligned} &\overset{\rho \gg 1}{\leq} \mathbb{E}_{\sigma^{2}_{\nu_{1}}}\left\{ \frac{1}{\pi} \int_{0}^{\pi/2}\left(\frac{\rho\sigma^{2}_{\nu_{1}}}{4\sin^{2}(\vartheta)}\right)^{-1}\,d\vartheta \right\} \end{aligned}} $$
$${\begin{aligned} &\,=\, \frac{1}{\rho}\mathbb{E}_{\sigma^{2}_{\nu_{1}}}\left\{\!\frac{1}{\sigma^{2}_{\nu_{1}}} \left(\underbrace{\frac{1}{\pi} \int_{0}^{\pi/2}\left(\frac{1}{4\sin^{2}(\vartheta)}\right)^{-1}\,d\vartheta}_{=1} \right) \right\} \end{aligned}} $$
$$ {\begin{aligned} &= \frac{1}{\rho}\mathbb{E}_{\sigma^{2}_{\nu_{1}}}\left\{\frac{1}{\sigma^{2}_{\nu_{1}}}\right\} \end{aligned}} $$
The rest of the proof follows by using similar steps as for the proof of Proposition 1.
By using the Maclaurin series expansion and keeping the first two dominant terms, we have the following:
$$\begin{array}{*{20}l} &\exp \left(- z\bar \zeta \left({K,\vartheta} \right)\psi \left({p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \right) \right) &\\&\quad= 1 - z\bar \zeta \left({K,\vartheta} \right)\psi \left({p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \right) + \\ &\qquad\:\:\:\:{\mathcal{O}}\left({z\bar \zeta \left({K,\vartheta} \right)\psi \left({p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \right)} \right)\\ &\quad\approx 1 - z\bar \zeta \left({K,\vartheta} \right)\psi \left({p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \right) \end{array} $$
By using this approximation, we have the following:
$$ \begin{aligned} &\int_{0}^{\pi} {\int_{- \pi }^{\pi} {\exp} } \left(- z\bar \zeta \left({K,\vartheta} \right)\psi \left({p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \right) \right)\\&\quad{f_{\theta} }\left({\theta^{t}}\right){f_{\phi} }\left({\phi^{t}}\right){\mkern 1mu} d{\theta^{t}}{\mkern 1mu} d{\phi^{t}}\\ &\approx \int_{0}^{\pi} \int_{- \pi }^{\pi} \left({1 - z\bar \zeta \left({K,\vartheta} \right)\psi \left({p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \right)} \right)\\&\quad {f_{\theta} }\left({\theta^{t}}\right){f_{\phi} }\left({\phi^{t}}\right){\mkern 1mu} d{\theta^{t}}{\mkern 1mu} d{\phi^{t}} \\ & = 1 - z\bar \zeta \left({K,\vartheta} \right)\int_{0}^{\pi} {\int_{- \pi }^{\pi} {\psi \left({p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \right)}}\\&\quad {f_{\theta} }\left({\theta^{t}}\right){f_{\phi} }\left({\phi^{t}}\right){\mkern 1mu} d{\theta^{t}}{\mkern 1mu} d{\phi^{t}} \end{aligned} $$
Let us consider the following integral:
$$\begin{array}{*{20}l} {} \Theta = \int_{0}^{\pi} {\int_{- \pi }^{\pi} {\psi \left({p,q,{x_{m}},{x_{n}},{\theta^{t}},{\phi^{t}}} \right)}} {f_{\theta} }\left({\theta^{t}}\right){f_{\phi} }\left({\phi^{t}}\right)d{\theta^{t}}d{\phi^{t}} \end{array} $$
In addition, the following holds true:
$$\begin{array}{*{20}l} {\underset{K \to + \infty}{\lim }} {\left({1 - \frac{{z\frac{\rho }{{4{{\sin }^{2}}(\vartheta)}}\Theta }}{K}} \right)^{K}} = \exp \left({ - z\frac{\rho }{{4{{\sin }^{2}}(\vartheta)}}\Theta} \right) \end{array} $$
where we used the following notable limit:
$$\begin{array}{*{20}l} {\underset{x \to + \infty}{\lim }} {\left({1 + \frac{k}{x}} \right)^{x}} = {e^{k}}, \end{array} $$
Therefore, the asymptotic APEP is:
$$\begin{array}{*{20}l} \text{APEP} &\!\leq\! \frac{1}{\pi }\int_{0}^{\pi /2} {\int_{0}^{\infty} {\exp \left({ \,-\, z\left({1 \,+\, \frac{\rho }{{4{{\sin }^{2}}(\vartheta)}}\Theta} \right)} \right)}} dzd\vartheta \notag\\ &{{\stackrel{(a)}{=}} \frac{1}{\pi }\int_{0}^{\pi /2} {{{\left({1 + \frac{\rho }{{{4{\sin }^{2}}(\vartheta)}}\Theta} \right)}^{- 1}}} d\vartheta} \notag\\ &{ {\stackrel{(b)}{=}} \frac{1}{2}\left({1 - \sqrt {\frac{{\rho \Theta }}{{\rho \Theta + 1}}}} \right)} \end{array} $$
where the first equality (a) comes from the fact that \(\int _{0}^{\infty } {\exp \left ({ - zc} \right) =} \frac {1}{c},c > 0\) and the second equality (b) follows from (5A.4a) in [25]. This concludes the proof.
By using steps similar to Proposition 3, the asymptotic APEP can be written as follows:
$$\begin{array}{*{20}l} \text{APEP} &\leq \frac{1}{\pi }\int_{0}^{\pi /2} {{{\left({1 + \frac{\rho }{{{4{\sin }^{2}}(\vartheta)}}\Theta} \right)}^{- 1}}} d\vartheta \notag\\ & {\stackrel{\rho \gg 1}{\le}} \frac{1}{{\rho \Theta }}\left({\underbrace {\frac{1}{\pi }\int_{0}^{\pi /2} {{{\left({\frac{1}{{4{{\sin }^{2}}(\vartheta)}}} \right)}^{- 1}}} {\mkern 1mu} d\vartheta }_{= 1}} \right), \end{array} $$
which concludes the proof.
We first introduce the following lemma [26] for application to hermitian quadratic forms in complex normal variables.
Lemma 1
Let vn(n=1,…,N) be a set of complex Gaussian random variables having zero mean. Let κ, with v=[v1,⋯,vN]T, be an Hermitian quadratic form:
$$\begin{array}{*{20}l} \kappa = \mathbf{v}^{H}\mathbf{I}_{N}\mathbf{v} \end{array} $$
Its MGF is as follows:
$$\begin{array}{*{20}l} M_{\kappa}(s) = \prod_{n = 1}^{N} (1-s\lambda_{n})^{-1} \end{array} $$
where λn is the nth eigenvalue of the covariance matrix \(\mathbf {R}_{\mathbf {v}} = \mathbb {E}\{\mathbf {v}\mathbf {v}^{H}\}\).
See [26]. □
The proof of Proposition 5 can be split in three steps.
Step 1: By using Lemma 1 with N=Nr=2, we have:
$$\begin{array}{*{20}l} \text{APEP} \approx \frac{1}{{{\rho^{2}}}}{\mathbb{E}_{{\theta_{k}^{t}},{\phi_{k}^{t}},{\theta_{k}^{r}},{\phi_{k}^{r}}}}\left\{ {\left({\frac{3}{{{\lambda_{1}}{\lambda_{2}}}}} \right)} \right\} \end{array} $$
where λ1 and λ2 are eigenvectors of the covariance matrix:
$$\begin{array}{*{20}l} \mathbf{R} = \left[ {\begin{array}{*{20}{c}} {\mathbb{E}\left\{ {\left. {{\nu_{1}}{\nu^{*}}_{1}} \right\}} \right.}&{\mathbb{E}\left\{ {\left. {{\nu_{1}}{\nu^{*}}_{2}} \right\}} \right.}\\ {\mathbb{E}\left\{ {\left. {{\nu_{2}}{\nu^{*}}_{1}} \right\}} \right.}&{\mathbb{E}\left\{ {\left. {{\nu_{2}}{\nu^{*}}_{2}} \right\}} \right.} \end{array}} \right] \end{array} $$
and R is assumed to be full rank.
It is worth mentioning that λ1 and λ2 depend on the random variables \({{\theta _{k}^{t}},{\phi _{k}^{t}},{\theta _{k}^{r}},{\phi _{k}^{r}}}\). In particular, from Lemma 1, we have:
$$\begin{array}{*{20}l} {} \text{APEP} = \mathbb{E}_{{\theta_{k}^{t}},{\phi_{k}^{t}},{\theta_{k}^{r}},{\phi_{k}^{r}}}\left\{ {\left. {\frac{1}{\pi }\int\limits_{0}^{\pi /2} {\prod\limits_{{n_{r}} = 1}^{2} {{{(1 + \frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{n_{r}}})}^{- 1}}} d\vartheta} } \right\}} \right. \end{array} $$
Also, we have the following:
$$\begin{array}{*{20}l} &\frac{1}{{\left({1 + \frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{1}}} \right)\left({1 + \frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{2}}} \right)}} \\&= \frac{{{A_{1}}}}{{\left({1 + \frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{1}}} \right)}} + \frac{{{A_{2}}}}{{\left({1 + \frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{2}}} \right)}} \end{array} $$
which yields:
$$\begin{array}{*{20}l} {A_{1}} = \frac{{{\lambda_{1}}}}{{{\lambda_{1}} - {\lambda_{2}}}}, \qquad {A_{2}} = - \frac{{{\lambda_{2}}}}{{{\lambda_{1}} - {\lambda_{2}}}} \end{array} $$
For high SNR, we have:
$$\begin{array}{*{20}l} {\left({1 + \frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{n_{r}}}} \right)^{- 1}} &\approx {\left({\frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{n_{r}}}} \right)^{- 1}}\\&\quad - {\left({\frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{n_{r}}}} \right)^{- 2}} \end{array} $$
where we have used the second-order Taylor approximation.
Thus, by substituting (48) and (49) in (46), we obtain:
$$ {} \begin{aligned} \text{APEP} &= {\mathbb{E}_{{\theta_{k}^{t}},{\phi_{k}^{t}},{\theta_{k}^{r}},{\phi_{k}^{r}}}}\left\{ \left. \sum\limits_{i = 1}^{2} {A_{i}}\left[ \frac{1}{\pi }\int\limits_{0}^{\pi /2} {{\left({\frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{i}}} \right)}^{- 1}}d\vartheta\right. \right. \right.\\ &\quad \left. \left. - \frac{1}{\pi }\int\limits_{0}^{\pi /2} {{{\left({\frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{i}}} \right)}^{- 2}}d\vartheta} \right] \right\} \\ &= {\mathbb{E}_{{\theta_{k}^{t}},{\phi_{k}^{t}},{\theta_{k}^{r}},{\phi_{k}^{r}}}}\left\{ \frac{1}{\rho }\left[ \frac{1}{{{\lambda_{i}}}}\frac{1}{\pi }\int\limits_{0}^{\pi /2} {{{\left({\frac{1}{{4{{\sin }^{2}}\vartheta }}} \right)}^{- 1}}d\omega} \right. \right.\\ &\quad \left.\left. - \frac{1}{\rho }\frac{1}{{{\lambda^{2}}_{i}}}\frac{1}{\pi }\int\limits_{0}^{\pi /2} {{{\left({\frac{1}{{4{{\sin }^{2}}\left(\vartheta \right)}}} \right)}^{- 2}}d\vartheta} \right] \right\} \\ & {\stackrel{(a)}{=}} { \mathbb{E}_{{\theta_{k}^{t}},{\phi_{k}^{t}},{\theta_{k}^{r}},{\phi_{k}^{r}}}}\left\{ \frac{1}{\rho }\left(\frac{{{\lambda_{1}}}}{{{\lambda_{1}} - {\lambda_{2}}}}\left[ {\frac{1}{{{\lambda_{1}}}} - \frac{3}{{\rho {\lambda^{2}}_{1}}}} \right] \right.\right.\\ &\quad \left.\left. - \frac{{{\lambda_{2}}}}{{{\lambda_{1}} - {\lambda_{2}}}}\left[ {\frac{1}{{{\lambda_{2}}}} - \frac{3}{{\rho {\lambda^{2}}_{2}}}} \right] \right) \right\} \\ &= {\mathbb{E}_{{\theta_{k}^{t}},{\phi_{k}^{t}},{\theta_{k}^{r}},{\phi_{k}^{r}}}}\left\{ {\frac{1}{{{\rho^{2}}}}\left({\frac{3}{{{\lambda_{1}}{\lambda_{2}}}}} \right)} \right\} \end{aligned} $$
where (a) follows from:
$$\begin{array}{*{20}l}{} \frac{1}{\pi }\int\limits_{0}^{\pi /2} {{{\left(\! {\frac{1}{{4{{\sin }^{2}}\left(\omega \right)}}}\! \right)}^{- 1}}d\omega} =\! 1, \;\; \frac{1}{\pi }\int\limits_{0}^{\pi /2} {{{\left(\! {\frac{1}{{4{{\sin }^{2}}\left(\omega \right)}}}\! \right)}^{- 2}}d\omega} = 3 \end{array} $$
Step 2: We compute the explicit expression of the product λ1λ2. To this end, we introduce the following lemma.
Assume that R, defined in (45), is full rank and has two distinct eigenvalues λ1 and λ2. The product of the two eigenvectors λ1 and λ2 is as follows:
$$ {\begin{aligned} &{\lambda_{1} \lambda_{2}} \,=\, {\mathbb{E}_{{\theta^{t}},{\phi^{t}},{\theta_{k}^{r}},{\phi_{k}^{r}}}}\left({\frac{1}{{{K^{2}}}}\left[ {\sum\limits_{k = 1}^{K} {\sum\limits_{k^{\prime} = 1}^{K} {\!\chi \left({{\theta_{k}}^{t},{\phi_{k}^{t}}} \right)\chi \left({{\theta_{k}}^{\prime{t}},\phi_{k}^{\prime{t}}} \right)}} \Upsilon} \right]} \right) \end{aligned}} $$
$$\begin{array}{*{20}l} \chi \left({{\theta_{k}}^{t},{\phi_{k}^{t}}} \right) = &\sqrt {{G_{q}}\left({{\theta_{k}}^{t},{\phi_{k}^{t}}} \right)} \exp \left({j{\Omega_{q}}\left({{\theta_{k}}^{t},{\phi_{k}^{t}}} \right)} \right) \\&- \sqrt {{G_{p}}\left({{\theta_{k}}^{t},{\phi_{k}^{t}}} \right)} \exp \left({j{\Omega_{p}}\left({{\theta_{k}}^{t},{\phi_{k}^{t}}} \right)} \right) \end{array} $$
$$\begin{array}{*{20}l} \Upsilon =1 - \cos \left(\left\| \mathbf{k} \right\|{d}\left(\sin \left({{\theta^{r}}} \right)\sin \left({{\phi^{r}}} \right) \right.\right. \\ -\left.\left. \sin \left({\theta^{\prime{r}}} \right)\sin \left({\phi^{\prime{r}}} \right) \right) \right) \end{array} $$
By definition:
$$\begin{array}{*{20}l} \nu_{1} &= \frac{1}{{\sqrt K }}\sum\limits_{k = 1}^{K} {{\beta_{k}}\chi \left({{\theta_{k}}^{t},{\phi_{k}^{t}}} \right)} \end{array} $$
$$\begin{array}{*{20}l} {\nu_{2}} &\,=\, \frac{1}{{\sqrt K }}\sum\limits_{k = 1}^{K} {\beta_{k}}{\chi \left({{\theta_{k}}^{t},{\phi_{k}^{t}}} \right)\!\exp \left({j\left\| k \right\|{d}\sin \left({{\theta_{k}^{r}}} \right)\sin \left({{\phi_{k}^{r}}} \right)} \right)} \end{array} $$
Moreover, it is known that the product of the eigenvalues is equal to the determinant of the covariance matrix (i.e., \(\det (\mathbf {R}) = {\lambda _{1}}{\lambda _{2}}\)). Thus, λ1λ2 can be computed directly from \(\det (\mathbf {R})\) as follows:
$$\begin{array}{*{20}l} \det \left(\mathbf{R} \right) \,=\, \mathbb{E}\left\{ {\left. {{\nu_{1}}{\nu^{*}}_{1}} \right\}} \right.\mathbb{E}\left\{ {\left. {{\nu_{2}}{\nu^{*}}_{2}} \right\}} \right. \,-\, \mathbb{E}\left\{ {\left. {{\nu_{2}}{\nu^{*}}_{1}} \right\}} \right.\mathbb{E}\left\{ {\left. {{\nu_{1}}{\nu^{*}}_{2}} \right\}} \right. \end{array} $$
By substituting (55) and (56) into (57), we obtain (52) where we have used the identity:
$$\frac{{\exp \left({j\phi} \right) + \exp \left({ - j\phi} \right)}}{2} = \cos \left(\phi \right)$$
and if \(k = k^{\prime }\) we have \( 1 - \exp \left (j\left \| {{\mathbf {k}}} \right \|{d}\left (\sin \left ({{\theta _{k}^{r}}} \right)\sin \left ({{\phi _{k}^{r}}} \right) - \sin \left ({\theta _{k}^{\prime {r}}} \right)\sin \left ({\phi _{k}^{\prime {r}}} \right) \right) \right) = 0\). □
Step 3: We exploit the asymptotic analysis introduced in Proposition 3 to derive a closed-form expression of the APEP. We apply the following formula that is a special case of (40) in [24] for X=λ1λ2:
$$\begin{array}{*{20}l} \mathbb{E}\left({{X^{- 1}}} \right) = \int_{0}^{\infty} {{M_{X}}\left({ - z} \right)dz} \end{array} $$
where \({M_{X}}\left (z \right) = \mathbb {E}_{X}\left \{\exp \left ({Xz} \right) \right \}\).
Then, we obtain:
$$ {\begin{aligned} {\mathbb{E}_{X}}\left({\exp \left({Xz} \right)} \right) &= {\mathbb{E}_{{\theta^{t}},{\phi^{t}},{\theta^{r}},{\phi^{r}}}}\left(\prod\limits_{k = 1}^{K} {\prod\limits_{k^{\prime} = 1}^{K} {\exp} } \left(\frac{1}{{{K^{2}}}}{\chi}\left({{\theta_{k}}^{t},{\phi_{k}^{t}}} \right)\right.\right.\\ &\quad \left.\left. {\chi}\left({{\theta_{k}}^{\prime{t}},\phi_{k}^{\prime{t}}} \right) \Upsilon z {\vphantom{\prod\limits_{k = 1}}}\right) \right)\\ &=\prod\limits_{k = 1}^{K} \prod\limits_{k^{\prime} = 1}^{K} \left(\int_{- \pi }^{\pi}\! \int_{- \pi }^{\pi}\! \int_{- \pi }^{\pi}\! \int_{- \pi }^{\pi}\! \int_{0}^{\pi}\! \int_{0}^{\pi}\!\! \int_{0}^{\pi}\! \int_{0}^{\pi} \exp\! \left(\! \frac{1}{{{K^{2}}}}\!{\chi}\left(\! {{\theta_{k}}^{t},{\phi_{k}^{t}}} \!\right)\!{\chi} \right.\right.\\ & \left.\left({{\theta_{k}}^{\prime{t}},\phi_{k}^{\prime{t}}} \right) \Upsilon z {\vphantom{\int_{- \pi }^{\pi}}}\right){f_{{\phi^{t}}_{k^{\prime}}}}{f_{{\phi^{t}}_{k}}}{f_{{\phi^{r}}_{k^{\prime}}}}{f_{{\phi^{r}}_{k}}}{f_{{\theta^{t}}_{k^{\prime}}}}{f_{{\theta^{t}}_{k}}}{f_{{\theta^{r}}_{k^{\prime}}}}{f_{{\theta^{r}}_{k}}}\\ &\left.\quad{\phi_{k^{\prime}}}^{t}d{\phi_{k}}^{t}d{\phi_{r^{\prime}}}^{r}d{\phi_{k}}^{r} d{\theta_{k}}^{\prime{t}}d{\theta_{k}}^{t}d{\theta_{r}}^{\prime{r}}d{\theta_{k}}^{r} {\vphantom{\int_{- \pi }^{\pi}}}\right)\\ &\,=\,\left(\! \int_{- \pi }^{\pi}\! \int_{- \pi }^{\pi}\! \int_{- \pi }^{\pi}\! \int_{- \pi }^{\pi}\! \int_{0}^{\pi}\!\! \int_{0}^{\pi}\!\! \int_{0}^{\pi}\!\! \int_{0}^{\pi} \!\exp\! \left(\! {\frac{1}{{{K^{2}}}}{\chi}\!\left(\! {{\theta^{t}},{\phi^{t}}}\! \right)\chi\!\left(\! {\theta^{\prime{t}}}\!,\phi^{\prime{t}} \right)\!\Upsilon z}\! \right)\right.\\ &\left.\quad{f_{{\phi^{t}}}}{f_{\phi^{\prime{t}}}}{f_{{\phi^{r}}}}{f_{\phi^{\prime{r}}}}{f_{{\theta^{t}}}}{f_{\theta^{\prime{t}}}}{f_{{\theta^{r}}}}{f_{\theta^{\prime{r}}}}d{\phi^{t}}d\phi^{\prime{t}}d{\phi^{r}}d\phi^{\prime{r}} d{\theta^{t}}d\theta^{\prime{t}}d{\theta^{r}}d\theta^{\prime{r}} {\vphantom{\int_{- \pi }^{\pi}}}\right)^{{K^{2}}} \end{aligned}} $$
where the last equality follows from the assumption of independent and identically distributed random variables.
Moreover, by using the following approximation:
$$\begin{array}{*{20}l} &\exp \left({\frac{1}{{{K^{2}}}}{\chi}\left({{\theta^{t}},{\phi^{t}}} \right)\chi\left({\theta^{\prime{t}},\phi^{\prime{t}}} \right)\Upsilon z} \right) \approx 1 \\ & \qquad+ \frac{1}{{{K^{2}}}}{\chi}\left({{\theta^{t}},{\phi^{t}}} \right)\chi\left({\theta^{\prime{t}},\phi^{\prime{t}}} \right)\Upsilon z \end{array} $$
in (59) and using the notable limit in (39), we obtain the following:
$$ {\begin{aligned} &{\mathbb{E}_{X}}\left({\exp \left({Xz} \right)} \right) {\stackrel{K \to \infty}{=}} \\ &\exp \left(z\int_{- \pi }^{\pi} \int_{- \pi }^{\pi} \int_{- \pi }^{\pi} \int_{- \pi }^{\pi} \int_{0}^{\pi} \int_{0}^{\pi} \int_{0}^{\pi} \int_{0}^{\pi} \left({{\chi}\left({{\theta^{t}},{\phi^{t}}} \right)\chi\left({\theta^{\prime{t}},\phi^{\prime{t}}} \right)\Upsilon} \right)\right.\\ &\quad\left.{f_{{\phi^{t}}}}{f_{\phi{\prime{t}}}}{f_{{\phi^{r}}}}{f_{\phi^{\prime{r}}}}{f_{{\theta^{t}}}}{f_{\theta^{\prime{t}}}}{f_{{\theta^{r}}}}{f_{\theta^{\prime{r}}}}d{\phi^{t}}d \phi^{\prime{t}}d{\phi^{r}}d\phi^{\prime{r}} d{\theta^{t}}d\theta^{\prime{t}}d{\theta^{r}}d\theta^{\prime{r}} {\vphantom{\int_{- \pi }^{\pi}}}\right) \end{aligned}} $$
Finally, substituting (61) in (58), and then (58) in (44), we conclude the proof by using the following result:
$$ {}\begin{aligned} & \int_{- \pi }^{\pi} \int_{0}^{\pi} \int_{- \pi }^{\pi} \int_{0}^{\pi} \left(1 - \cos \left(\left\| \mathbf{k} \right\|{d}\left(\sin \left({{\theta^{r}}} \right)\sin \left({{\phi^{r}}} \right)\right.\right.\right.\\ &\, \left.\left.\left.- \sin \left({\theta^{\prime{r}}} \right)\sin \left({\phi^{\prime{r}}} \right) \right) \right) \right) {f_{{\phi^{r}}}}{f_{\phi^{\prime{r}}}}{f_{{\theta^{r}}}}{f_{\theta^{\prime{r}}}}d{\phi^{r}}d{\theta^{r}}d \phi^{\prime{r}}d \theta^{\prime{r}} = \\ &1 - {\left({\int_{- \pi }^{\pi} {\int_{0}^{\pi} {\cos \left({\left\| \mathbf{k} \right\|{d}\left({\sin \left({{\theta^{r}}} \right)\sin \left({{\phi^{r}}} \right)} \right)} \right)} {f_{{\phi^{r}}}}{f_{{\theta^{r}}}}d{\phi^{r}}d{\theta^{r}}}} \right)^{2}}\\ &\!\!\!\!\,\quad - {\left({\int_{- \pi }^{\pi} {\int_{0}^{\pi}\! {\sin \left({\left\| \mathbf{k} \right\|{d}\left({\sin \left({{\theta^{r}}} \right)\sin \left({{\phi^{r}}} \right)} \right)} \right){f_{{\phi^{r}}}}{f_{{\theta^{r}}}}d{\phi^{r}}d{\theta^{r}}}} } \right)^{2}} \end{aligned} $$
Proof of Theorem 1
The proof of Theorem 1 generalizes the steps of Proposition 5 as follows.
Step 1: From Lemma 1, we have:
$$\begin{array}{*{20}l}{} \text{APEP} = {\mathbb{E}_{{\theta_{k}^{t}},{\phi_{k}^{t}},{\theta_{k}^{r}},{\phi_{k}^{r}}}}\left\{ {\left. {\frac{1}{\pi }\int\limits_{0}^{\pi /2} {\prod\limits_{n = 1}^{{N_{r}}} {{{(1 + \frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{n}})}^{- 1}}} d\vartheta} } \right\}} \right. \end{array} $$
Then, we need to identify the coefficients \(\left \{ {{A_{{n_{r}}}}} \right \}_{{n_{r}} = 1}^{{N_{r}}}\) so that the following is satisfied:
$$\begin{array}{*{20}l} \frac{1}{{\prod\limits_{{n_{r}} = 1}^{{N_{r}}} {\left({1 + \frac{\rho }{{4{{\sin }^{2}}\left(\omega \right)}}{\lambda_{n_{r}}}} \right)} }} = \sum\limits_{{n_{r}} = 1}^{{N_{r}}} {\frac{{{A_{n_{r}}}}}{{\left({1 + \frac{\rho }{{4{{\sin }^{2}}\left(\omega \right)}}{\lambda_{n_{r}}}} \right)}}} \end{array} $$
To this end, we can exploit the general formula of partial fraction decomposition [pp. 66–67], [27]:
$$\begin{array}{*{20}l} {A_{n_{r}}} = \frac{{{\lambda_{n_{r}}}^{{N_{r}} - 1}}}{{\prod\limits_{{n_{r}} \ne l = 1}^{{N_{r}}} {\left({{\lambda_{n_{r}}} - {\lambda_{l}}} \right)} }} \end{array} $$
for l,nr=1,…,Nr.
Thus, we have:
$$ {\begin{aligned} \text{APEP} &\!\approx\! {\mathbb{E}_{_{{\theta^{t}},{\phi^{t}},{\theta^{r}},{\phi^{r}}}}}\!\left\{\!\! {\left. {\frac{1}{\pi }\!\int\limits_{0}^{\pi /2} {\sum\limits_{{n_{r}} = 1}^{{N_{r}}}\! {{A_{n}}\!\left[\! {\sum\limits_{m = 1}^{{N_{r}}}\!\! {{{\left(\! { - 1} \right)}^{m + 1}}{{\left(\!\! {\frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{{n_{r}}}}}\!\! \right)}^{\! - m}}}}\!\! \right]} d\vartheta} }\! \right\}} \right. \\ &\!= {\mathbb{E}_{{\theta^{t}},{\phi^{t}},{\theta^{r}},{\phi^{r}}}}\left\{ {\sum\limits_{{n_{r}} = 1}^{{N_{r}}} {{A_{n}}\left[ {\sum\limits_{m = 1}^{{N_{r}}} {{{\left({ - 1} \right)}^{m + 1}}\frac{1}{\pi }\int_{0}^{\pi /2}}}\right.}}\right.\\ &\quad\left.{{\left.{{{{\left({\frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{{n_{r}}}}}\! \right)}^{\! - m}}d\vartheta} }\! {\vphantom{\sum\limits_{m = 1}^{{N_{r}}}}}\right]}}\! {\vphantom{\sum\limits_{{n_{r}} = 1}^{{N_{r}}}}}\right\} \\ &\!= {\mathbb{E}_{_{{\theta^{t}},{\phi^{t}},{\theta^{r}},{\phi^{r}}}}}\left\{ {\sum\limits_{{n_{r}} = 1}^{{N_{r}}} {A_{n}}\left[ {\sum\limits_{m = 1}^{{N_{r}}} {{{\left({ - 1} \right)}^{m + 1}}{{\left({\rho {\lambda_{{n_{r}}}}} \right)}^{- m}}}}\right.}\right.\\ &\quad\left.{{\left.{{\left({\frac{1}{\pi }\int_{0}^{\pi /2} {{{\left({\frac{1}{{4{{\sin }^{2}}\left(\vartheta \right)}}} \right)}^{- m}}d\vartheta} } \right)}} \right]}} \right\} \end{aligned}} $$
In the high SNR regime, we can use the Nrth-order Taylor approximation as follows:
$$ {\begin{aligned} {\left({1 + \frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{{n_{r}}}}}\! \right)^{- {N_{r}}}} &\approx \sum\limits_{m = 1}^{{N_{r}}} {{{\left({ - 1} \right)}^{m + 1}}{{\left({\frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{{n_{r}}}}}\! \right)}^{- m}}}\\ &\quad + \mathcal{O}\left({{{\left({ - 1} \right)}^{{N_{r}} + 1}}{{\left({\frac{\rho }{{4{{\sin }^{2}}\left(\vartheta \right)}}{\lambda_{{n_{r}}}}} \right)}^{- ({N_{r}} + 1)}}} \right) \end{aligned}} $$
We note that the following holds true:
$$ {\begin{aligned} \frac{1}{\pi }\!\! \int\limits_{0}^{\pi /2}\! {{{\left(\! {\frac{1}{{4{{\sin }^{2}}\!\left(\vartheta \right)}}}\! \right)}^{\! - m}}\! d\vartheta} &= \frac{1}{2}\! \left(\! {\!\begin{array}{*{20}{c}} {2m}\\ m \end{array}}\! \right)\\ &\quad +\! \sum\limits_{k = 0}^{m - 1}\! {{{\left(\! { - 1}\! \right)}^{m - k}}2\!\left(\! {\begin{array}{*{20}{c}} {2n}\\ k \end{array}}\! \right)\!\frac{{\sin \left({\pi\! \left({m - k} \right)} \right)}}{{2\pi \left({m - k} \right)}}} \end{aligned}} $$
where we have used the following identity [27, p. 31, eq. (1320,1)]:
$$ {}\begin{aligned} {\sin^{2m}}\left(\vartheta \right) &= \frac{1}{{{2^{2m}}}}\left\{\! {\sum\limits_{k = 0}^{m - 1} {{\left({ - 1} \right)}^{m - k}}2\left({\begin{array}{*{20}{c}} {2m}\\ k \end{array}} \right)\!\cos \left({2\left({m - k} \right)\vartheta} \right)}\right.\\ &\quad\left.{+ \left({\begin{array}{*{20}{c}} {2m}\\ m \end{array}}\! \right)} \right\} \end{aligned} $$
Let us define:
$$\begin{array}{*{20}l} {\alpha_{m}} = \left({\frac{1}{\pi }\int_{0}^{\pi /2} {{{\left({\frac{1}{{4{{\sin }^{2}}\left(\vartheta \right)}}} \right)}^{- m}}d\vartheta} } \right) \end{array} $$
Then, we have the following:
$$ {\begin{aligned} \text{APEP} &\approx {\mathbb{E}_{{\theta_{k}^{t}},{\phi_{k}^{t}},{\theta_{k}^{r}},{\phi_{k}^{r}}}}\left\{ {\sum\limits_{{n_{r}} = 1}^{{N_{r}}} {{A_{{n_{r}}}}\left[ {\sum\limits_{m = 1}^{{N_{r}}} {{{\left({ - 1} \right)}^{m + 1}}{{\left({\rho {\lambda_{{n_{r}}}}} \right)}^{- m}}{\alpha_{m}}}}\! \right]}}\! \right\} \\ &= {\mathbb{E}_{{\theta_{k}^{t}},{\phi_{k}^{t}},{\theta_{k}^{r}},{\phi_{k}^{r}}}}\left\{ {\sum\limits_{{n_{r}} = 1}^{{N_{r}}} {{A_{{n_{r}}}}\left[ {{{\left({ - 1} \right)}^{{N_{r}} + 1}}{{\left({\rho {\lambda_{{n_{r}}}}} \right)}^{- {N_{r}}}}{\alpha_{{N_{r}}}}}\right.}}\right.\\ &\quad\left.{\left.{ + \sum\limits_{m = 1}^{{N_{r}} - 1} {{{\left({ - 1} \right)}^{m + 1}}{{\left({\rho {\lambda_{{n_{r}}}}} \right)}^{- m}}{\alpha_{m}}}} \right]} \right\} \\ &= {\mathbb{E}_{{\theta_{k}^{t}},{\phi_{k}^{t}},{\theta_{k}^{r}},{\phi_{k}^{r}}}}\left\{ {\sum\limits_{{n_{r}} = 1}^{{N_{r}}} {{A_{{n_{r}}}}{{\left({ - 1} \right)}^{{N_{r}} + 1}}{{\left({\rho {\lambda_{{n_{r}}}}} \right)}^{- {N_{r}}}}{\alpha_{{N_{r}}}}}}\right.\\ &\quad\left.{+ \sum\limits_{n = 1}^{{N_{r}}} {{A_{n}}\left[ {\sum\limits_{m = 1}^{{N_{r}} - 1} {{{\left({ - 1} \right)}^{m + 1}}{{\left({\rho {\lambda_{{n_{r}}}}} \right)}^{- m}}{\alpha_{m}}}} \right]}} \right\} \end{aligned}} $$
By induction, the following can be proved:
$$\begin{array}{*{20}l} \sum\limits_{{n_{r}} = 1}^{{N_{r}}} {{A_{{n_{r}}}}\left[ {\sum\limits_{m = 1}^{{N_{r}} - 1} {{{\left({ - 1} \right)}^{m + 1}}{{\left({\rho {\lambda_{{n_{r}}}}} \right)}^{- m}}{\alpha_{m}}}} \right]} = 0 \end{array} $$
$$\begin{array}{*{20}l} \sum\limits_{{n_{r}} = 1}^{{N_{r}}} {{A_{{n_{r}}}}{{\left({ - 1} \right)}^{{N_{r}} + 1}}{{\left({\rho {\lambda_{{n_{r}}}}} \right)}^{- {N_{r}}}}{\alpha_{{N_{r}}}}} = \frac{{{\alpha_{{N_{r}}}}}}{{{\rho^{{N_{r}}}}}}\frac{1}{{\prod\limits_{n = 1}^{{N_{r}}} {{\lambda_{n}}} }} \end{array} $$
Step 2: Let \(X = \prod \limits _{{n_{r}} = 1}^{{N_{r}}} {{\lambda _{{n_{r}}}}}\). By using steps similar to those of Lemma 2, we obtain the following:
$$\begin{array}{*{20}l} \det \left({{\mathbf{R}_{\mathbf{v}}}} \right) &= \frac{1}{{{K^{{N_{r}}}}}}\left(\sum\limits_{{k_{1}} = 1}^{K} \ldots \sum\limits_{{k_{{N_{r}}}} = 1}^{K} \left(\prod\limits_{i = 1}^{{N_{r}}} {\chi}\left({{\theta_{{k_{i}}}}^{t},\phi_{{k_{i}}}^{t}} \right)\right. \right. \\& \left. \left. F\left({\theta_{{k_{1}}}^{r}, \ldots,\theta_{{k_{{N_{r}}}}}^{r},\phi_{{k_{1}}}^{r}, \ldots,\phi_{{k_{{N_{r}}}}}^{r}} {\vphantom{\sum\limits_{{k_{1}} = 1}}}\right) \right) \right) \end{array} $$
where \(F\left ({\theta _{{k_{1}}}^{r}, \ldots,\theta _{{k_{{N_{r}}}}}^{r},\phi _{{k_{1}}}^{r}, \ldots,\phi _{{k_{{N_{r}}}}}^{r}} \right)\).
Step 3: The APEP is then the following:
$$ {}{\begin{aligned} &\mathbb{E}_{X}\left\{ {\exp \left({Xz} \right)} \right\}\\ &={\mathbb{E}_{{\theta^{t}},{\phi^{t}},{\theta^{r}},{\phi^{r}}}}\left\{ \exp\!\! \left(\!\! \frac{z}{{{K^{{N_{r}}}}}}\!\!\left(\sum\limits_{{k_{1}} = 1}^{K} \ldots\! \sum\limits_{{k_{{N_{r}}}} = 1}^{K} \left[\! \left(\prod\limits_{i = 1}^{{N_{r}}} {{\chi}\left({{\theta_{{k_{i}}}}^{t},\phi_{{k_{i}}}^{t}} \right)} \!\right)\right.\right. \right.\right.\\& \left.\left.\left.\left. F\left({\theta_{{k_{1}}}^{r}, \ldots,\theta_{{k_{{N_{r}}}}}^{r},\phi_{{k_{1}}}^{r}, \ldots,\phi_{{k_{{N_{r}}}}}^{r}} \right) \vphantom{{\sum\limits_{{k_{1}} = 1}^{K}}}\right] \right) \right) \right\}\\ &={\mathbb{E}_{{\theta^{t}},{\phi^{t}},{\theta^{r}},{\phi^{r}}}}\left\{ \prod\limits_{{k_{1}} = 1}^{K} \ldots \prod\limits_{{k_{{N_{r}}}} = 1}^{K} \exp \left(\frac{z}{{{K^{{N_{r}}}}}}\left[ {\prod\limits_{i = 1}^{{N_{r}}} {{\chi}\left({{\theta_{{k_{i}}}}^{t},\phi_{{k_{i}}}^{t}} \right)}} \right]\right.\right. \\& \left. \left. F\left({\theta_{{k_{1}}}^{r}, \ldots,\theta_{{k_{{N_{r}}}}}^{r},\phi_{{k_{1}}}^{r}, \ldots,\phi_{{k_{{N_{r}}}}}^{r}} \right) {\vphantom{\prod\limits_{{k_{{N_{r}}}} = 1}^{K}}}\right) \right\} \\ &=\prod\limits_{{k_{1}} = 1}^{K} \ldots \prod\limits_{{k_{{N_{r}}}} = 1}^{K} {\mathbb{E}_{{\theta^{t}},{\phi^{t}},{\theta^{r}},{\phi^{r}}}}\left\{\exp \left(\frac{z}{{{K^{{N_{r}}}}}}\left[ {\prod\limits_{i = 1}^{{N_{r}}} {{\chi}\left({{\theta_{{k_{i}}}}^{t},\phi_{{k_{i}}}^{t}} \right)}} \right]\right. \right. \\& \left. \left. F\left({\theta_{{k_{1}}}^{r}, \ldots,\theta_{{k_{{N_{r}}}}}^{r},\phi_{{k_{1}}}^{r}, \ldots,\phi_{{k_{{N_{r}}}}}^{r}} \right) {\vphantom{\prod\limits_{i = 1}^{{N_{r}}}}}\right) \right\} \\ &=\left[{\mathbb{E}_{{\theta^{t}},{\phi^{t}},{\theta^{r}},{\phi^{r}}}}\left\{ \exp \left(\frac{z}{{{K^{{N_{r}}}}}}\left[ {\prod\limits_{i = 1}^{{N_{r}}} {{\chi^{i}}\left({{\theta^{i,t}},{\phi^{i,t}}} \right)}} \right] \right.\right.\right. \\& \left.\left. \left. F\left({{\theta^{1,r}}, \ldots,{\theta^{{N_{r}},r}},{\phi^{1,r}}, \ldots,{\phi^{{N_{r}},r}}} \right) {\vphantom{\prod\limits_{i = 1}^{{N_{r}}}}}\right) \right\} \right]^{{N_{r}}} \end{aligned}} $$
By using the following approximation:
$$ {\begin{aligned} {}&\exp\!\left(\! {\frac{z}{{{K^{{N_{r}}}}}}\!\left[ {\prod\limits_{i = 1}^{{N_{r}}} \!{{\chi}\!\left({{\theta^{i,t}},{\phi^{i,t}}} \right)}} \!\!\right]\!F\!\left({{\theta^{1,r}}, \!\ldots,{\theta^{{N_{r}},r}},{\phi^{1,r}}, \ldots,{\phi^{{N_{r}},r}}} \right)}\! \right) \approx \end{aligned}} $$
$$ \begin{aligned} {}&1 + \frac{z}{{{K^{{N_{r}}}}}}\left[ {\prod\limits_{i = 1}^{{N_{r}}} {{\chi}\left({{\theta^{i,t}},{\phi^{i,t}}} \right)}}\! \right]\!F\left({{\theta^{1,r}}, \ldots,{\theta^{{N_{r}},r}},{\phi^{1,r}}, \ldots,{\phi^{{N_{r}},r}}} \right) \end{aligned} $$
$$ {}{\begin{aligned} \mathbb{E}\left\{ {\exp \left({Xz} \right)} \right\}&=\left[ {\mathbb{E}_{{\theta^{t}},{\phi^{t}},{\theta^{r}},{\phi^{r}}}}\left\{1 + \frac{z}{{{K^{{N_{r}}}}}}\left[ {\prod\limits_{i = 1}^{{N_{r}}} {{\chi}\left({{\theta^{i,t}},{\phi^{i,t}}} \right)}} \right]\right. \right. \\& \left. \left. F\left({{\theta^{1,r}}, \ldots,{\theta^{{N_{r}},r}},{\phi^{1,r}}, \ldots,{\phi^{{N_{r}},r}}} \right) {\vphantom{\prod\limits_{i = 1}^{{N_{r}}}}}\right\} \right]^{{N_{r}}}\\ &=\left[ 1 + \frac{z}{{{K^{{N_{r}}}}}}{\mathbb{E}_{{\theta^{t}},{\phi^{t}},{\theta^{r}},{\phi^{r}}}}\left\{ \left[ {\prod\limits_{i = 1}^{{N_{r}}} {{\chi}\left({{\theta^{i,t}},{\phi^{i,t}}} \right)}} \right] \right. \right. \\& \left. \left. F\left({{\theta^{1,r}}, \ldots,{\theta^{{N_{r}},r}},{\phi^{1,r}}, \ldots,{\phi^{{N_{r}},r}}} \right) \prod\limits_{i = 1}^{{N_{r}}}\right\} \right]^{{N_{r}}}\\ &\approx \exp \left(z{\mathbb{E}_{{\theta^{t}},{\phi^{t}},{\theta^{r}},{\phi^{r}}}}\left\{ \left[ {\prod\limits_{i = 1}^{{N_{r}}} {{\chi}\left({{\theta^{i,t}},{\phi^{i,t}}} \right)}} \right]\right. \right. \\& \left. \left. F\left({{\theta^{1,r}}, \ldots,{\theta^{{N_{r}},r}},{\phi^{1,r}}, \ldots,{\phi^{{N_{r}},r}}} \right) \prod\limits_{i = 1}^{{N_{r}}}\right\} \right) \end{aligned}} $$
Thus, we obtain:
$$ {\begin{aligned} &\frac{{{\alpha_{{N_{r}}}}}}{{{\rho^{{N_{r}}}}}}\int_{0}^{\infty} {{M_{X}}\left({ - z} \right)dz} = \\ &\frac{{{\alpha_{{N_{r}}}}}}{{{\rho^{{N_{r}}}}}}\int_{0}^{\infty} \exp \left(z{\mathbb{E}_{{\theta^{t}},{\phi^{t}},{\theta^{r}},{\phi^{r}}}}\left\{\left[ \prod\limits_{i = 1}^{{N_{r}}} {{\chi}\left({{\theta^{i,t}},{\phi^{i,t}}} \right)} \right]\right. \right. \\ & \left. \left. F\left({{\theta^{1,r}}, \ldots,{\theta^{{N_{r}},r}},{\phi^{1,r}}, \ldots,{\phi^{{N_{r}},r}}} \right) {\vphantom{\prod\limits_{i =1}^{{N_{r}}}}}\right\} \right)dz \\ &= \frac{{{\alpha_{{N_{r}}}}}}{{{\rho^{{N_{r}}}}{\mathbb{E}_{{\theta^{t}},{\phi^{t}},{\theta^{r}},{\phi^{r}}}}\left\{ {\left[ {\prod\limits_{i = 1}^{{N_{r}}} {{\chi}\left({{\theta^{i,t}},{\phi^{i,t}}} \right)}} \right]F\left({{\theta^{1,r}}, \ldots,{\theta^{{N_{r}},r}},{\phi^{1,r}}, \ldots,{\phi^{{N_{r}},r}}} \right)} \right\}}} \\ &= \frac{{{\alpha_{{N_{r}}}}}}{{{\rho^{{N_{r}}}}{\mathbb{E}_{{\theta^{t}},{\phi^{t}}}}\left\{ {\left[ {\prod\limits_{i = 1}^{{N_{r}}} {{\chi}\left({{\theta^{i,t}},{\phi^{i,t}}} \right)}} \right]} \right\}{\mathbb{E}_{{\theta^{r}},{\phi^{r}}}}\left\{{F\left({{\theta^{1,r}}, \ldots,{\theta^{{N_{r}},r}},{\phi^{1,r}}, \ldots,{\phi^{{N_{r}},r}}} \right)} \right\}}}\\ &= \frac{{{\alpha_{{N_{r}}}}}}{{{\rho^{{N_{r}}}}{{\left[ {{\mathbb{E}_{{\theta^{t}},{\phi^{t}}}}\left\{ {\chi\left({{\theta^{t}},{\phi^{t}}} \right)} \right\}} \right]}^{{N_{r}}}}{\mathbb{E}_{{\theta^{r}},{\phi^{r}}}}\left\{ {F\left({{\theta^{1,r}}, \ldots,{\theta^{{N_{r}},r}},{\phi^{1,r}}, \ldots,{\phi^{{N_{r}},r}}} \right)} \right\}}} \end{aligned}} $$
This concludes the proof.
From Theorem 1, we can obtain (21). We need to prove how to obtain the explicit form of the function F(·). If Nr=3, the determinant of the covariance matrix R can be computed as follows:
$$ {} \begin{aligned} &\det \left({{\mathbf{R}}} \right) = \\ &\mathbb{E}\left\{ {\left. {{\nu_{1}}{\nu^{*}}_{1}} \right\}} \right.\mathbb{E}\left\{ {\left. {{\nu_{2}}{\nu^{*}}_{2}} \right\}} \right.\mathbb{E}\left\{ {\left. {{\nu_{3}}{\nu^{*}}_{3}} \right\}} \right. + \mathbb{E}\left\{ {\left. {{\nu_{1}}{\nu^{*}}_{2}} \right\}} \right.\mathbb{E}\left\{ {\left. {{\nu_{2}}{\nu^{*}}_{3}} \right\}} \right.\\& \mathbb{E}\left\{ {\left. {{\nu_{3}}{\nu^{*}}_{1}} \right\}} \right. + \mathbb{E}\left\{ {{\nu_{1}}{\nu^{*}}_{3}} \right\}\mathbb{E}\left\{ {{\nu_{2}}{\nu^{*}}_{1}} \right\}\mathbb{E}\left\{ {{\nu_{3}}{\nu^{*}}_{2}} \right\} - \mathbb{E}\left\{ {{\nu_{1}}{\nu^{*}}_{3}} \right\}\\& \mathbb{E}\left\{ {{\nu_{2}}{\nu^{*}}_{2}} \right\} \mathbb{E}\left\{ {{\nu_{3}}{\nu^{*}}_{1}} \right\} - \mathbb{E}\left\{ {\left. {{\nu_{1}}{\nu^{*}}_{2}} \right\}} \right.\mathbb{E}\left\{ {\left. {{\nu_{2}}{\nu^{*}}_{1}} \right\}} \right.\mathbb{E}\left\{ {\left. {{\nu_{3}}{\nu^{*}}_{3}} \right\}} \right.\\ & - \mathbb{E}\left\{ {\left. {{\nu_{1}}{\nu^{*}}_{1}} \right\}} \right.\mathbb{E}\left\{ {\left. {{\nu_{2}}{\nu^{*}}_{3}} \right\}} \right.\mathbb{E}\left\{ {\left. {{\nu_{3}}{\nu^{*}}_{2}} \right\}} \right. \end{aligned} $$
We note that:
$$\begin{array}{*{20}l} &\mathbb{E}\left\{ {\left. {{\nu_{u}}{\nu^{*}_{v}}} \right\}} \right. =\notag\\&\quad\!\! { \frac{1}{K}\sum\limits_{k = 1}^{K} {{\chi\left({{{\theta^{t}_{k}}},{{\phi^{t}_{k}}}} \right)}} \exp \left({j\left({u - v} \right)\left\| \mathbf{k} \right\|{d}\sin \left({{\theta_{k}^{r}}} \right)\sin \left({{\phi_{k}^{r}}} \right)} \right)} \end{array} $$
for u,v=1,…,Nr. Then, from (80) and (81), we have:
$$\begin{array}{*{20}l} &\det(\mathbf{R})=\frac{1}{{{K^{3}}}}\sum\limits_{{k_{1}} = 1}^{K} {\sum\limits_{{k_{2}} = 1}^{K} {\sum\limits_{{k_{3}} = 1}^{K} {{\chi}\left({{\theta_{{k_{1}}}}^{t},\phi_{{k_{1}}}^{t}} \right)}} } {\chi}\left({{\theta_{{k_{2}}}}^{t},\phi_{{k_{2}}}^{t}} \right)\notag\\&\qquad{\chi}\left({{\theta_{{k_{3}}}}^{t},\phi_{{k_{3}}}^{t}} \right)F\left({\theta_{{k_{1}}}^{r},\theta_{{k_{2}}}^{r},\theta_{{k_{3}}}^{r},\phi_{{k_{1}}}^{r},\phi_{{k_{2}}}^{r},\phi_{{k_{3}}}^{r}} \right) \end{array} $$
$$ {\begin{aligned} &F\left({\theta_{{k_{1}}}^{r},\theta_{{k_{2}}}^{r},\theta_{{k_{3}}}^{r},\phi_{{k_{1}}}^{r},\phi_{{k_{2}}}^{r},\phi_{{k_{3}}}^{r}} \right) =\\ &\left[\! \begin{array}{l} 1\\ +\exp\! \left({ \!- j\left\| k \right\|{d}\left[ \!{\sin \left({\theta_{{k_{1}}}^{r}} \right)\sin\! \left({\phi_{{k_{1}}}^{r}}\! \right) \,+\, \sin\! \left({\theta_{{k_{2}}}^{r}} \!\right)\sin\! \left({\phi_{{k_{2}}}^{r}}\! \right)\! -\! 2\sin\! \left({\theta_{{k_{3}}}^{r}} \right)\!\sin\! \left({\phi_{{k_{3}}}^{r}} \!\right)} \!\right]} \right)\\ + \exp\! \left({ \!- j\left\| k \right\|{d}\left[ \!{2\sin \left({\theta_{{k_{1}}}^{r}} \!\right)\sin\! \left({\phi_{{k_{1}}}^{r}} \!\right) \,-\, \sin\! \left({\theta_{{k_{2}}}^{r}} \!\right)\sin\! \left({\phi_{{k_{2}}}^{r}} \!\right) \,-\, \sin\! \left({\theta_{{k_{3}}}^{r}}\! \right)\sin\! \left({\phi_{{k_{3}}}^{r}} \right)}\! \right]} \right)\\ - \exp \left({ - j\left\| k \right\|{d}\left[ \!{2\sin\! \left({\theta_{{k_{1}}}^{r}} \right)\sin\! \left({\phi_{{k_{1}}}^{r}} \right) - 2\sin\! \left({\theta_{{k_{3}}}^{r}} \right)\sin\! \left({\phi_{{k_{3}}}^{r}} \right)} \right]} \right)\\ - \exp \left({ - j\left\| k \right\|{d}\left[ \!{\sin\! \left({\theta_{{k_{1}}}^{r}} \right)\sin\! \left({\phi_{{k_{1}}}^{r}} \right) - \sin\! \left({\theta_{{k_{2}}}^{r}} \right)\sin\! \left({\phi_{{k_{2}}}^{r}} \right)} \right]} \right)\\ - \exp \left({ - j\left\| k \right\|{d}\left[ \!{\sin\! \left({\theta_{{k_{2}}}^{r}} \right)\sin\! \left({\phi_{{k_{2}}}^{r}} \right) - \sin\! \left({\theta_{{k_{3}}}^{r}} \right)\sin\! \left({\phi_{{k_{3}}}^{r}} \right)} \right]} \right) \end{array} \right] \end{aligned}} $$
which can be further simplified as follows:
$$ {\begin{aligned} &F\left({\theta_{{k_{1}}}^{r},\theta_{{k_{2}}}^{r},\theta_{{k_{3}}}^{r},\phi_{{k_{1}}}^{r},\phi_{{k_{2}}}^{r},\phi_{{k_{3}}}^{r}} \right) =\\ &\left[\! \begin{array}{l} 1\\ + \cos \left({ \,-\, \left\| k \right\|{d}\left[ \!{\sin\! \left({\theta_{{k_{1}}}^{r}} \!\right)\sin\! \left({\phi_{{k_{1}}}^{r}} \right) \,+\, \sin\! \left({\theta_{{k_{2}}}^{r}} \!\right)\sin \!\left({\phi_{{k_{2}}}^{r}}\! \right)\! -\! 2\sin \!\left({\theta_{{k_{3}}}^{r}} \right)\sin \left({\phi_{{k_{3}}}^{r}} \right)} \right]} \right)\\ + \cos \left({ - \left\| k \right\|{d}\left[ {2\sin\! \left({\theta_{{k_{1}}}^{r}} \right)\!\sin\! \left({\phi_{{k_{1}}}^{r}}\! \right)\! -\! \sin \left({\theta_{{k_{2}}}^{r}} \!\right)\!\sin \!\left({\phi_{{k_{2}}}^{r}} \right) \,-\, \sin\! \left({\theta_{{k_{3}}}^{r}}\! \right)\!\sin\! \left({\phi_{{k_{3}}}^{r}}\! \right)} \right]} \right)\\ - \cos \left({ - \left\| k \right\|{d}\left[ {2\sin\! \left({\theta_{{k_{1}}}^{r}} \right)\sin\! \left({\phi_{{k_{1}}}^{r}} \right) - 2\sin\! \left({\theta_{{k_{3}}}^{r}} \right)\sin\! \left({\phi_{{k_{3}}}^{r}} \right)} \right]} \right)\\ - \cos \left({ - \left\| k \right\|{d}\left[ {\sin\! \left({\theta_{{k_{1}}}^{r}} \right)\sin\! \left({\phi_{{k_{1}}}^{r}} \right) - \sin\! \left({\theta_{{k_{2}}}^{r}} \right)\sin\! \left({\phi_{{k_{2}}}^{r}} \right)} \right]} \right)\\ - \cos \left({ - \left\| k \right\|{d}\left[ {\sin\! \left({\theta_{{k_{2}}}^{r}} \right)\sin\! \left({\phi_{{k_{2}}}^{r}} \right) - \sin\! \left({\theta_{{k_{3}}}^{r}} \right)\sin\! \left({\phi_{{k_{3}}}^{r}} \right)} \right]} \right) \end{array}\! \right] \end{aligned}} $$
where we have used similar observations as for Nr=2.
Now, we need to compute \(\mathbb {E}_{{\theta ^{r}},{\phi ^{r}}}\left ({F\left ({{\theta ^{1,r}},{\theta ^{2,r}},{\theta ^{3,r}},{\phi ^{1,r}},{\phi ^{2,r}},{\phi ^{3,r}}} \right)} \right)\). Let us define the following:
$$\begin{array}{*{20}l} \begin{array}{l} a = \left\| k \right\|{d}\sin \left({\theta_{{k_{1}}}^{r}} \right)\sin \left({\phi_{{k_{1}}}^{r}} \right)\\ b = \left\| k \right\|{d}\sin \left({\theta_{{k_{2}}}^{r}} \right)\sin \left({\phi_{{k_{2}}}^{r}} \right)\\ c = \left\| k \right\|{d}\sin \left({\theta_{{k_{3}}}^{r}} \right)\sin \left({\phi_{{k_{3}}}^{r}} \right) \end{array} \end{array} $$
By applying the following trigonometric identities:
$$ {} \begin{aligned} &\cos \left({a + b - 2c} \right) = \cos \left(a \right)\cos \left(b \right)\cos \left({2c} \right)\\& \quad - \sin \left(a \right)\sin \left(b \right)\cos \left({2c} \right) + \sin \left(a \right)\cos \left(b \right)\sin \left({2c} \right)\\&\quad + \sin \left(b \right)\cos \left(a \right)\sin \left({2c} \right) \\ &\cos \left({c + b - 2a} \right) = \cos \left(c \right)\cos \left(b \right)\cos \left({2a} \right)\\&\quad - \sin \left(c \right)\sin \left(b \right)\cos \left({2a} \right) + \sin \left(c \right)\cos \left(b \right)\sin \left({2a} \right)\\ &\quad + \sin \left(b \right)\cos \left(c \right)\sin \left({2a} \right) \end{aligned} $$
and noting that a,b, and c are independent, we eventually obtain the desired result after some algebraic manipulations.
M. Di Renzo, H. Haas, P. M. Grant, Spatial modulation for multiple-antenna wireless systems: a survey. IEEE Commun. Mag.49(12), 182–191 (2011).
A. Younis, S. Sinanovic, M. Di Renzo, R. Y. Mesleh, H. Haas, Generalised sphere decoding for spatial modulation. IEEE Trans. Commun.61(7), 2805–2815 (2013).
A. Stavridis, S. Sinanovic, M. Di Renzo, H. Haas, in 2013 IEEE 78th Vehicular Technology Conference (VTC Fall). Energy evaluation of spatial modulation at a multi-antenna base station (IEEELas Vegas, 2013), pp. 1–5. https://doi.org/10.1109/VTCFall.2013.6692187.
D. -T. Phan-Huy, et al., Single-carrier spatial modulation for the Internet of Things: design and performance evaluation by using real compact and reconfigurable antennas, submitted (2018). https://arxiv.org/abs/1812.07514. https://doi.org/10.1109/ACCESS.2019.2895754.
P. Patcharamaneepakorn, S. Wu, C. -X. Wang, E. -H. M. Aggoune, M. M. Alwakeel, X. Ge, M. Di Renzo, Spectral, energy, and economic efficiency of 5G multicell massive MIMO systems with generalized spatial modulation. IEEE Trans. Veh. Technol.65(12), 9715–9731 (2016).
M. Di Renzo, H. Haas, A. Ghrayeb, S. Sugiura, L. Hanzo, Spatial modulation for generalized MIMO: challenges, opportunities and implementation. Proc. IEEE. 102(1), 56–103 (2014).
P. Yang, M. Di Renzo, Y. Xiao, S. Li, L. Hanzo, Design guidelines for spatial modulation. IEEE Commun. Surv. Tuts.16(1), 6–26 (2015).
P. Yang, Y. Xiao, Y. L. Guan, K. V. S. Hari, A. Chockalingam, S. Sugiura, H. Haas, M. Di Renzo, C. Masouros, Z. Liu, L. Xiao, S. Li, L. Hanzo, Single-carrier SM-MIMO: a promising design for broadband large-scale antenna systems. IEEE Commun. Surv. Tuts.18(3), 1687–1716 (2016).
M. Di Renzo, H. Haas, A. Ghrayeb, L. Hanzo, S. Sugiura, in Wiley Encyclopedia of Electrical and Electronics Engineering. Spatial Modulation for Multiple-Antenna Communication, (2016). https://doi.org/10.1002/047134608X.W8327.
E. Basar, M. Wen, R. Mesleh, M. Di Renzo, Y. Xiao, H. Haas, Index modulation techniques for next-generation wireless networks. IEEE Access. 5:, 16693–16746 (2017). https://doi.org/10.1109/ACCESS.2017.2737528.
A. Younis, W. Thompson, M. Di Renzo, C. -X. Wang, M. A. Beach, H. Haas, P. M. Grant, Performance of spatial modulation using measured real-world channels. IEEE Veh. Technol. Conf. - Fall, 1–5 (2013).
N. Serafimovski, A. Younis, R. Mesleh, P. Chambers, M. Di Renzo, C. -X. Wang, P. M. Grant, M. A. Beach, H. Haas, Practical implementation of spatial modulation. IEEE Trans. Veh. Technol.62(9), 4511–4523 (2013).
P. Liu, M. Di Renzo, A. Springer, Line-of-sight spatial modulation for indoor mmWave communication at 60 GHz. IEEE Trans. Wirel. Commun.15(11), 7373–7389 (2016).
P. Liu, J. Blumenstein, N. S. Perovic, M. Di Renzo, A. Springer, Performance of generalized spatial modulation MIMO over measured 60 GHz indoor channels. IEEE Trans. Commun.66(1), 133–148 (2018).
M. Di Renzo, H. Haas, Bit error probability of SM-MIMO over generalized fading channels. IEEE Trans. Veh. Technol.61(3), 1124–1144 (2012).
M. Di Renzo, H. Haas, On transmit-diversity for spatial modulation MIMO: impact of spatial-constellation diagram and shaping filters at the transmitter. IEEE Trans. Veh. Technol.62(6), 2507–2531 (2013).
M. -T. Le, V. -D. Ngo, H. -A. Mai, X. -N. Tran, M. Di Renzo, Spatially modulated orthogonal space-time block codes with nonvanishing determinants. IEEE Trans. Commun.62(1), 85–99 (2014).
D. A. Basnayaka, M. Di Renzo, H. Haas, Massive but few active MIMO. IEEE Trans. Veh. Technol.65(9), 6861–6877 (2016).
N. S. Perovic, P. Liu, J. Blumenstein, M. Di Renzo, A. Springer, Optimization of the cut-off rate of generalized spatial modulation with transmit precoding. IEEE Trans. Commun.66(10), 4578–4595 (2018).
M. Di Renzo, Spatial modulation based on reconfigurable antennas - a new air interface for the IoT. IEEE Mil. Commun. Conf., 1–6 (2017).
Q. U. A. Nadeem, A. Kammoun, M. Debbahn, M. S. Alouini, A generalized spatial correlation model for 3D MIMO channels based on the Fourier coefficients of power spectrums. IEEE Trans. Sig. Proc.63(14), 3671–3686 (2015).
J. Zhang, C. Pan, F. Pei, G. Liu, X. Cheng, Three-dimensional fading channel models: a survey of elevation angle research. IEEE Commun. Mag.52(6), 218–226 (2014).
T. S. Pollock, T. D. Abhayapala, R. A. Kennedy, Three-dimensional fading channel models: a survey of elevation angle research. Telecommun. Syst.24(2), 415–436 (2003).
M. Di Renzo, C. Merola, A. Guidotti, F. Santucci, G. E. Corazza, Error performance of multi-antenna receivers in a Poisson field of interferers: a stochastic geometry approach. IEEE Trans. Commun.61(5), 2025–2047 (2013).
MK Simon, M-S Alouini, Digital Communication over Fading Channels (Wiley, 2005).
GL Turin, The characteristic function of Hermitian quadratic forms in complex normal variables. Biometrika. 67(1/2), 199–201 (1960).
A Jeffrey, D Zwillinger, Table of Integrals, Series, and Products (Elsevier, 2007).
Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study. The paper is built upon mathematical analysis and experimental validation.
Laboratoire des Signaux et Systèmes, CNRS, CentraleSupelec, Univ Paris-Sud, Université Paris-Saclay, Plateau du Moulon, Gif-sur-Yvette, 91192, France
Dung Nguyen Viet & Marco Di Renzo
TTI Norte, Parque Cientifico y Tecnologico de Cantabria, Albert Einstein 14, Santander, 39011, Cantabria, Spain
Vedaprabhu Basavarajappa & Beatriz Bedia Exposito
University of Cantabria, Av. de los Castros, s/n, Santander, 39005, Cantabria, Spain
Jose Basterrechea
Orange Labs, Orange Gardens, 44 avenue de la Republique, CS 50010, Chatillon Cedex, 92326, France
Dinh-Thuy Phan-Huy
Dung Nguyen Viet
Marco Di Renzo
Vedaprabhu Basavarajappa
Beatriz Bedia Exposito
The authors declare that they have equally contributed to the paper. All authors read and approved the final manuscript.
Correspondence to Marco Di Renzo.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Nguyen Viet, D., Di Renzo, M., Basavarajappa, V. et al. Spatial modulation based on reconfigurable antennas: performance evaluation by using the prototype of a reconfigurable antenna. J Wireless Com Network 2019, 149 (2019). https://doi.org/10.1186/s13638-019-1427-z
Spatial modulation
Reconfigurable antennas
Antenna prototype
Error probability
5Gwireless - Innovative Architectures, Wireless Technologies and Tools for High Capacity and Sustainable 5G Ultra-Dense Cellular Networks