DOI QR코드

DOI QR Code

Joint Antenna Selection and Multicast Precoding in Spatial Modulation Systems

  • Wei Liu (College of Physical Science and Technology, Central China Normal University) ;
  • Xinxin Ma (College of Physical Science and Technology, Central China Normal University) ;
  • Haoting Yan (The Research Institute for Frontier Science, Beihang University) ;
  • Zhongnian Li (College of Physical Science and Technology, Central China Normal University) ;
  • Shouyin Liu (College of Physical Science and Technology, Central China Normal University)
  • Received : 2023.07.07
  • Accepted : 2023.09.25
  • Published : 2023.11.30

Abstract

In this paper, the downlink of the multicast based spatial modulation systems is investigated. Specifically, physical layer multicasting is introduced to increase the number of access users and to improve the communication rate of the spatial modulation system in which only single radio frequency chain is activated in each transmission. To minimize the bit error rate (BER) of the multicast based spatial modulation system, a joint optimizing algorithm of antenna selection and multicast precoding is proposed. Firstly, the joint optimization is transformed into a mixed-integer non-linear program based on single-stage reformulation. Then, a novel iterative algorithm based on the idea of branch and bound is proposed to obtain the quasioptimal solution. Furthermore, in order to balance the performance and time complexity, a low-complexity deflation algorithm based on the successive convex approximation is proposed which can obtain a sub-optimal solution. Finally, numerical results are showed that the convergence of our proposed iterative algorithm is between 10 and 15 iterations and the signal-to-noise-ratio (SNR) of the iterative algorithm is 1-2dB lower than the exhaustive search based algorithm under the same BER accuracy conditions.

Keywords

1. Introduction

Spatial modulation is an energy efficient transmission scheme, in which only one antenna is activated during each transmission. Hence only one radio frequency (RF) chain is needed in the spatial modulation systems (SMSs). In SMSs, the input signals are transmitted at both the transmit antenna (TA) domain and amplitude-phase modulation (APM) domain [1]. Although SMSs utilize the spatial dimension to transmit information bits, they lack transmit diversity due to limited number of active transmit antennas. Therefore, the achievable transmission rate is imposed a few limitations at the APM domain. On the other hand, physical layer multicasting is a spectral efficient transmission scheme, in which a common message is multicasted to a couple of users [2]. The integration of multicast and spatial modulation is highly desired. The multicast based SMSs can achieve good balance between the energy efficiency and spectral efficiency.

When channel state information is available at the base station, transmit antenna selection and transmit precoding can be employed to enhance the performance of SMSs. Specifically, the diversity order can be obtained by employing transmit antenna selection [3], [4]. Moreover, for SMSs, when the total number of transmit antennas is not the power of 2 in practice, the antenna selection process is also needed [4]. On the other hand, the precoding can be optimized to shape the received APM signal constellation and boost the bit error rate (BER) performance [5]–[7]. Note that in SMSs, only one antenna is activated in each transmission period and the precoding matrix should be diagonal [7].

In the references of multicast spatial modulation, reference [9] firstly analyzes the BER performance of the multicast SMSs and provide a closed-form asymptotic BER upper by exploiting the system statistics. Finally, reference [9] studies the impact of the receiver number on BER in the multicast spatial modulation. Furthermore, reference [10] proposes two heuristic transmit mode selection algorithms, including tree pruning technique and approximating the signal constellation. Then reference [10] also design a precoder approach to improve the system. Similarly, reference [8] considers both antenna selection and the precoding coefficients, and gives the optimization problem. However, the optimization problem was not solved directly and reference [8] uses exhaustive search algorithm and convex optimization tools in turn to solve antenna selection and precoding design problems.

In this paper, antenna selection and multicast precoding in multicast based SMSs are design jointly. To facilitate the joint optimization of antenna selection and precoding, a single-stage reformulation is proposed to reconstruct the problem of joint optimization into a mix-integer non-linear program. Secondly, a branch and bound based algorithm is designed to obtain the quasi-optimal solution, whose convergence behavior and performance superiority are verified by the numerical results. Specifically, when the number of iterations reaches 10 to 15, the proposed branch and bound based algorithm converges with the accuracy of 0.0001, which is defined as the gap between the obtained upper bounds and lower bounds. This shows that the proposed algorithm converges fast, thus proving the effectiveness of joint antenna selection and multicast precoding in SMS. Moreover, a low-complexity deflation algorithm based on successive convex approximation (SCA) is proposed to obtain a suboptimal solution, which can achieve BER performance close to the exhaustive search based algorithm. To the best of our knowledge, this is the first work studying the single-stage reformulation for jointly solving the antenna selection and precoding problems. Meanwhile, under the same BER condition, the SNR corresponding to the deflation algorithm achieve 1dB lower than the result of the exhaustive search based algorithm.

Notations: We employ uppercase boldface letters for matrices and lowercase boldface for vectors. E(⦁), Tr(⦁), ║⦁║p, (⦁)T, (⦁)H and Re(⦁) return the expectation, trace, pl -norm, transpose, conjugate transpose and the real part of the input respectively.

2. System Description

As shown in Fig. 1, we consider the downlink of a multicasting SMS where a base station (BS) multicasts a common message to K users simultaneously using both the TA domain and APM domain signal transmission [1]. The BS is equipped with Ntot antennas whereas each user equipment (UE) has NR antennas. For brevity, we denote K = {1, 2,...,K} and Ntot = {1, 2,..., Ntot} as the user set and transmit antenna set, respectively.

E1KOBZ_2023_v17n11_3204_f0001.png 이미지

Fig. 1. Downlink of a multicasting SMS.

2.1 Signal Model

In each transmission, the BS selects NSM antennas from Ntot to perform spatial modulation, where NSM is power of 2 [3]. There are total CNSMNtot possible combinations of antenna selection. I is denoted as the set of selected NSM antennas with |I| = NSM, and L as the set of enumerations of all possible antenna selections with |I| = CNSMNtot. Let Hk ∈ CNRxNtot be the channel matrix from the BS to the k -th user. Then after antenna selection, the channel of user k used for spatial modulation becomes Hk,I ∈ CNSMxNtot where I ∈ L. In each transmission, the first log2 NSM bits of the transmitted signal are mapped to one point in the spatial constellation set defined as

S = {e1, ..., ei, ..., eNSM}       (1)

where ei ∈ CNSMx1 is an all-zero vector except for the i-th entry, which is 1. The selection of ei indicates that the i-th antenna among all the selected NSM antennas is activated. Next, the last log2 M bits are transmitted from the activated antenna by M-QAM or M-PSK signal constellation defined by

M = {s1, ..., sm, ..., sM}       (2)

where sm ∈ C with E{|sm|2} = 1

The spectral efficiency in the described system can be computed as log2 NSM M bits per channel use. The transmitted spatial modulation constellation SM with cardinality NSMM is the Cartesian product of S and M which is given by

SM = {e1s1, ..., eism, ..., eNSMsM}       (3)

Next, denote the transmitted spatial modulation symbol as xl = eism ∈ CNSMx1 with xl ∈ SM. Then xl is precoded by a diagonal matrix W = diag(w) ∈ CNSMxNSM with w = [wI1, ...,WISM]T ∈ CNSMx1 given I. Then the received signal at user k given antenna selection I ∈ L can be written as

yk = Hk,IWxl + nk = hk,IiWIiSm + nk       (4)

Where, Ii is the i-th element in I and hk,Ii is the i-th column in Hk,I which also is the Ii-th column in the total channel matrix Hk. And wIi is the precoding coefficient for antenna Ii in Ntot. Moreover, nk ~ CN(0,σ2INR) denotes the additive noise at user k.

At the receiver side, user k decodes the TA domain signal (i.e., index of the activated antenna) and APM domain signal jointly with maximum likelihood (ML) detector as follows

\(\begin{align}\hat{x}_{l}^{\mathrm{ML}}=\arg \min _{\hat{x}_{l} \in \mathcal{S} M}\left\|\mathbf{y}_{k}-\mathbf{H}_{k, I} \mathbf{W} \hat{x}_{l}\right\|_{2}^{2}, \forall k \in \mathcal{K}\end{align}\)       (5)

The error performance of the ML detector at user k can be approximated by the sum of the pairwise error probability (PEP), which is given by [5]

\(\begin{align}P_{k}^{\mathrm{Errror}}(I, \mathbf{w}) \leq \sum_{l=1}^{N_{\mathrm{SM}} M} \sum_{t=1, t \neq l}^{N_{\mathrm{sM}} M} Q\left(\sqrt{\frac{1}{2 \sigma^{2}} d_{k, l t}(I, \mathbf{w})}\right)\end{align}\)       (6)

where, \(\begin{align}Q(x)=\frac{1}{\sqrt{2 \pi}} \int_{x}^{\infty} e^{-\frac{y^{2}}{2}} d y\end{align}\) denotes the Gaussian tail probability [1], and \(\begin{align}d_{k, l t}(I, \mathbf{w})=\left\|\mathbf{H}_{k, I} \mathbf{W}\left(x_{l}-x_{t}\right)\right\|_{2}^{2}\end{align}\) is the squared Euclidian distance between two different spatial modulation symbols xl ∈ SM and xt ∈ SM observed at receiver k. Similar to [1] and [5], we employ the nearest neighbor approximation for PEP in (6) and we have

\(\begin{align}\tilde{p}_{k}^{\mathrm{e}}(I, \mathbf{w})=\lambda Q\left(\sqrt{\frac{1}{2 \sigma^{2}} d_{\min }^{k}(I, \mathbf{w})}\right)\end{align}\)       (7)

where, λ is the number of neighboring constellation points [1] and

\(\begin{align}\begin{array}{l}d_{\text {min }}^{k}(I, \mathbf{w})=\min _{\substack{x_{l}, x_{x} \in \mathcal{S} \\ x_{l} \neq x_{i}}}^{k}\left\|\mathbf{H}_{k, I} \mathbf{W}\left(x_{l}-x_{t}\right)\right\|_{2}^{2} \\ =\min _{\substack{i, j \in \mathcal{N}_{\mathrm{S}}, i \leq j \\ S_{m}, S_{n}, \mathcal{M} \\(i, m) \neq(j, n)}}\left\|\mathbf{h}_{k, I_{i}} \mathrm{~W}_{I_{i}} S_{m}-\mathbf{h}_{k, I_{j}} \mathrm{~W}_{I_{j}} S_{n}\right\|_{2}^{2} \\\end{array}\end{align}\)       (8)

In multicast system, the system performance is dominated by the worst user. Similar to (9), the BER performance of the considered multicasting SMS is defined by the worst BER among all the users, and can be written by

\(\begin{align}\tilde{P}^{\mathrm{e}}(I, \mathbf{w})=\max _{k \in \mathcal{K}} \tilde{P}_{k}^{\mathrm{e}}(I, \mathbf{w})\end{align}\)       (9)

2.2 Problem Formulation

In this paper, the BER of the multicasting SMS by jointly optimizing the transmit antenna selection and multicast precoding subject is minimized to the power constraint at the BS, which can be formulated as

\(\begin{align}\begin{array}{c}\min _{I, \mathbf{w}} \max _{k \in \mathcal{K}} \tilde{P}_{k}^{\mathrm{e}}(I, \mathbf{w}) \end{array} \end{align}\)       (10a)

s.t. ||w||22 ≤ PT       (10b)

where, PT is the total power constraint at the BS.

Since \(\begin{align}\tilde{P}_{k}^{\mathrm{e}}(I, \mathbf{w})\end{align}\) in (7) is a monotonically decreasing function of dkmin(I,w), the objective (10a) is equivalent to maxI,w mink∈Kdkmin(I,w). Hence, problem (10) can be equivalently written as

\(\begin{align}\max _{I, \mathbf{W}} \min _{k \in \mathcal{K}} \min _{\substack{i, j \in \mathcal{N}_{\mathrm{SM}}, i \leq j \\ S_{m}, S_{n} \in \mathcal{M} \\(i, m) \neq(j, n)}}\left\|\mathbf{h}_{k, I_i} \text {w}_{I_i} S_{m}-\mathbf{h}_{k, I_j} \text {w}_{I_j} s_{n}\right\|_{2}^{2}\end{align}\)       (11)

s.t. ||w||22 ≤ PT       (11b)

3. Single-stage Reformulation

In the formulation (11), the antenna selection variable I is a set of antenna indices, which cannot be designed jointly with precoding coefficients w unless enumerating all possible antenna combinations as in formulation (8). In what follows, we integrate the antenna selection and multicast precoding into an equivalent single-stage reformulation of (11), by which we can leverage the tools of optimization to design them jointly. To this end, we introduce a binary vector b = [b1,...,bNtot] to denote the antenna selection, i.e., bi =1 if antenna i is selected; otherwise, bi = 0. Since there are NSM antennas selected, we have \(\begin{align}\sum_{i=1}^{N_{\mathrm{tot}}} b_{i}=N_{\mathrm{SM}}\end{align}\). Note that there is an one-to-one mapping between b and I. Next, using the idea of Big-M formulation (11), formulation (11) is reformulated as follows

\(\begin{align}\max _{I, \mathbf{W}} \min _{k \in \mathcal{K}} \min _{\substack{i, j \in \mathcal{N}_{\mathrm{SM}}, i \leq j \\ s_{m}, S_{m} \in \mathcal{M} \\(i, m) \neq(j, n)}}\left\|\mathbf{h}_{k, i} v_{i} S_{m}-\mathbf{h}_{k, j} v_{j} s_{n}\right\|_{2}^{2}+\Phi_{i j}^{k}\end{align}\)       (12a)

s.t. bi∈{0,1}, ∀i∈Ntot       (12b)

\(\begin{align}\sum_{i=1}^{N_{\text {tot }}} b_{i}=N_{S M}\end{align}\)       (12c)

|vi|2 ≤ biPT, ∀i ∈ Ntot       (12d)

||v||22≤PT       (12e)

where Φkij = (1 - bi)Ωk + (1 - bj)Ωk. Constant Ωk has a sufficiently large value which should be chosen such that when bi = 0 or bj = 0 the term ||hk,ivism - hk,jvjsm||22 + Φkij is inactive in the second min function of (12a) for every feasible channel realization, precoding coefficients, and APM symbols. A simple choice could be \(\begin{align}\Omega_{k}=\max _{i \in \mathcal{N}_{\text {tot }}} 2 P_{T}\left\|\mathbf{h}_{k, i}\right\|_{2}^{2}\end{align}\). The variable vi is the precoding coefficient. Besides, v = [v1,...,vNtot]T ∈ CNtotx1. Since there are NSM antennas selected from Ntot, we have ||v||0 = NSM which is guaranteed by constraints (12c) and (12d). We note that the indices i, j in (12) take values from 1 to Ntot, which is different from i, j in (11) taking values from 1 to NSM.

Proposition 1: Problems (11) and (12) are equivalent.

Proof: First, based on the definition of binary variable bi =1 and constraint (12c), we can see that there is an one-to-one mapping between and I ∈ L in (11) and b in (12). Also, constraints (12c) and (12d) ensure ||v||0 = NSM and only the selected antennas have non-zero precoding coefficients. Hence, there is an one-to-one mapping between w in (11) and v in (12).

On the other hand, for the objective in (12a), we denote

\(\begin{align}\tilde{d}_{\min }^{k}(\mathbf{b}, \mathbf{v})=\min _{\substack{i, j \in \mathcal{N}_{\mathrm{SM}}, i \leq j \\ s_{m}, s_{n} \in \mathcal{N} \\(i, m) \neq(j, n)}}\left\|\mathbf{h}_{k, i} v_{i} s_{m}-\mathbf{h}_{k, j} v_{j} s_{n}\right\|_{2}^{2}+\Phi_{i j}^{k}\end{align}\)       (13)

Then to prove that (11) and (12) are equivalent, we only need to prove \(\begin{align}\tilde{d}_{\text {min }}^{k}(\mathbf{b}, \mathbf{v})=\tilde{d}_{\text {min }}^{k}(I, \mathbf{w})\end{align}\) for any given antenna selection, precoding coefficients and user index. Suppose I = {I1,...,INSM} are the selected NSM antennas and the corresponding precoding coefficients are w = {wI1,..., wINSM}. Let tot Ic = Ntot\I, i.e., the relative complement of set I with respect to set Ntot, denote the set of unused antennas. Hence we have bi =1 and vi = wi if i ∈ I; otherwise, if i ∈ Ic, we have bi = 0 and vi = 0. Then it follows that \(\begin{align}\tilde{d}_{\text {min }}^{k}(\mathbf{b}, \mathbf{v})=\min \left\{\tilde{d}_{\min }^{k, I}(\mathbf{b}, \mathbf{v}), \tilde{d}_{\text {min }}^{k, I^{c}}(\mathbf{b}, \mathbf{v})\right\}\end{align}\), where

\(\begin{align}\tilde{d}_{\min }^{k, I}(\mathbf{b}, \mathbf{v})=\min _{\substack{i, j \in I, i \leq j \\ s_{m}, s_{n} \in \mathcal{M} \\(i, m) \neq(j, n)}}\left\|\mathbf{h}_{k, i} v_{i} s_{m}-\mathbf{h}_{k, j} v_{j} s_{n}\right\|_{2}^{2}\end{align}\)       (14)

\(\begin{align}\tilde{d}_{\min }^{k, I^{c}}(\mathbf{b}, \mathbf{v})=\min _{\substack{i \in I^{c} \| j \in I^{c}, i \leq j, j \\ s_{m}, s_{n} \in \mathcal{M} \\(i, m) \neq(j, n)}}\left\|\mathbf{h}_{k, i} v_{i} s_{m}-\mathbf{h}_{k, j} v_{j} s_{n}\right\|_{2}^{2}+\Phi_{i j}^{k}\end{align}\)       (15)

where i ∈ IC || j ∈ IC means that at least one element in {i, j} belongs to Ic, i.e., at least one element in {bi, bj} is zero.

Note that in (14), i and j are the indices of selected transmit antennas and bi = bj = 1. Hence, (1 - bi)Ωk + (1 - bj)Ωk = 0.

Recalling (8), we can obtain that \(\begin{align}\tilde{d}_{\text {min }}^{k, I}(\mathbf{b}, \mathbf{v})=\tilde{d}_{\text {min }}^{k}(\mathbf{b}, \mathbf{v})\end{align}\). Since \(\begin{align}\Omega_{k}=\max _{i \in \mathcal{N}_{\text {tot }}} 2 P_{T}\left\|\mathbf{h}_{k, i}\right\|_{2}^{2}\end{align}\), we have \(\begin{align}\tilde{d}_{\text {min }}^{k, I}(\mathbf{b}, \mathbf{v})<\tilde{d}_{\min }^{k, I^{c}}(\mathbf{b}, \mathbf{v})\end{align}\) which results in \(\begin{align}\tilde{d}_{\min }^{k}(\mathbf{b}, \mathbf{v})=\tilde{d}_{\min }^{k, I}(\mathbf{b}, \mathbf{v})=\tilde{d}_{\min }^{k}(I, \mathbf{w})\end{align}\). Hence, (11) and (12) are equivalent.

4. Joint Antenna Selection and Multicast Precoding Design

The single-stage reformulation in problem (12) is a mixed-integer non-linear program (MINLP), which is NP-hard to solve. Besides, the objective in (12a) is non-smooth and non-convex. To proceed, we first show that problem (12) can be further reformulated into a mixed-integer quadratically constrained quadratic programming (MI-QCQP) problem.

4.1 Equivalent formulation

Proposition 2: The MINLP problem in (12) is equivalent to the following MI-QCQP problem

\(\begin{align}\max _{\mathbf{b}, \mathbf{v}} \min _{k \in \mathcal{K}} \min _{\substack{i, j \in \mathcal{N}_{\text {tot }} \\ i \leq j \\ s_{m}, s_{n} \in \mathcal{M} \\(i, m) \neq(j, n)}} \mathbf{v}^{H} \mathbf{R}_{i m, j n}^{k} \mathbf{v}+\Phi_{i j}^{k}\end{align}\)       (16a)

s.t. bi ∈ {0,1}, ∀i ∈ Ntot       (16b)

\(\begin{align}\sum_{i=1}^{N} b_{i}=N_{S M}\end{align}\)       (16c)

|vi|2 ≤ biPT, ∀i ∈ Ntot       (16d)

||v||22 ≤ PT       (16e)

where, Rkim,jn = HHkHkTim,jn is positive semidefinite matrix in which the matrix Θim,jn = (eism - ejsn)(eism - ejsn)H. The operator o denotes the Hadamard product. Moreover, ei ∈ CNtotx1 is an all-zero vector except for the i-th entry, which is 1.

Proof: In order to see the equivalence between the objective functions in (16a) and (12a), we denote the distance as dkim,jn = ||hk,ivism - hk,jvjsn||2. Next, introduce a diagonal matrix V = diag(v) ∈ CNtotxNtot where diag(v) represents the diagonal matrix which its main diagonal elements are the components of the vector v. Then we have

\(\begin{align}\begin{aligned} d_{i m, j n}^{k} & =\left\|\mathbf{H}_{k} \mathbf{V}\left(\mathbf{e}_{i} s_{m}-\mathbf{e}_{j} s_{n}\right)\right\|^{2} \\ & =\operatorname{Tr}\left\{\mathbf{V}^{H} \mathbf{H}_{k}^{H} \mathbf{H}_{k} \mathbf{V}\left(\mathbf{e}_{i} s_{m}-\mathbf{e}_{j} s_{n}\right)\left(\mathbf{e}_{i} s_{m}-\mathbf{e}_{j} s_{n}\right)^{H}\right\} \\ & =\operatorname{Tr}\left\{\mathbf{V}^{H} \mathbf{H}_{k}^{H} \mathbf{H}_{k} \mathbf{V} \boldsymbol{\Theta}_{i m, j n}\right\} \\ & \stackrel{\text { (a) }}{=} \mathbf{v}^{H}\left(\mathbf{H}_{k}^{H} \mathbf{H}_{k} \circ \boldsymbol{\Theta}_{i m, j n}^{T}\right) \mathbf{v} \\ & =\mathbf{v}^{H} \mathbf{R}_{i m, j n}^{k} \mathbf{V}\end{aligned}\end{align}\)      (17)

where, (a) follows from Eq. (1.10.6) in (12). Hence, the equivalence of problems (16) and (12) has been proved.

Next, the MI-QCQP problem (16) can be further converted (16) to the following equivalent problem

\(\begin{align}\max _{\mathbf{b}, \mathbf{v}, t} t\end{align}\)       (18a)

s.t. vHRkim,jnv + Φkij ≥ t, ∀k ∈ K, i, j ∈N, i ≤ j       (18b)

\(\begin{align}\begin{array}{r}s_{m}, s_{n} \in \mathcal{M},(i, m) \neq(j, n) \\ (16 \mathrm{~b})-(16 \mathrm{e})\end{array}\end{align}\)       (18c)

4.2 Quasi-Optimal Solution Based on Branch and Bound

In this subsection, we design a quasi-optimal algorithm for solving (18) by employing the branch and bound framework, in which branching and bounding are main tasks. In general, the branching step divides the feasible region into subsets and constructs sub-problems with those subsets. The bounding step finds the upper and lower bounds for those sub-problems within the corresponding subset. Specifically, we denote Q = − {b, v|(16b) - (16e)} as the feasible region. Then we branch the feasible region Q according to integer variables b , and divide the original problem into sub-problems accordingly.

Branching Step: Pick any index i , and we can divide the origin set Q into two subsets Q1 = − {b, v|bi = 0, (16b) - (16e)} and Q2 = {b, v|bi = 1, (16b) - (16e)}. Then the bounding steps described below are performed to obtain the upper and lower bounds of the original problem in these subsets respectively. Next, the subset with respect to the current overall upper bound is chosen to be further branched.

Bounding Step: In this step, we calculate the upper and lower bounds for the objective of (18) in different subsets. For convenience, we introduce ξ = [ξ1,...,ξNtot] ∈ RNtotx1 as the state vector for the integer variables b and the corresponding subset Qξ in the branching step, the elements of which take values from {0,1, -1}. Specifically, ξi =1 or 0 indicates that bi =1 or 0. If ξi = −1, bi has not been branched yet. Denote Aξ = {i ∈ Ntoti = 0or1}. Then subset Qξ which is divided by Aξ can be written as Qξ = {b, v|bi = ξi(i ∈ Aξ), (16b) - (16e)}.

Next, to obtain the upper bound in subset Qξ, we relax the undetermined binary variables bi(i∈Ntot\Aξ) into 0 ≤ bi ≤ 1, and then solve problem (18) using the SCA framework [13], which convexity the nonconvex constraint (18b) and iteratively approximate the nonconvex problem to convex ones. Specifically, at SCA iteration l, the left-hand-side of (18b) can be lower bounded by its first-order Taylor series expansion around vl, and (18b) can approximated by

2Re(Vl,HRkim,jnv) - Vl,HRkim,jnvl + Φki,j≥ t       (19)

Then the approximated convex problem at SCA iteration l can be written as

\(\begin{align}\max _{\mathbf{b}, \mathbf{v}, t} t\end{align}\)       (20a)

s.t. bi = ξi, i ∈ Aξ       (20b)

0 ≤ bi≤ 1, i ∈Ntot\Aξ       (20c)

(16c) - (16e), (19)       (20d)

which is convex and can be solved by modern solvers. The SCA algorithm is summarized for obtaining v giving ξ in Algorithm 1. The value of t at convergence is the upper bound.

Algorithm 1 SCA Algorithm for Obtaining (v,t) v t Given ξ

JAKO202305234914567_algor1.png 이미지

Function [U*,b*] = UB(ξ) is defined to represent the processing for computing the upper bound, in which U* is the obtained upper bound and b* = [b*1,...,b*Ntot] is not integral. In order to get the lower bound for the objective of (18) in Qξ, round b*i(i ∈ Aξ) to 0 or 1, and denote the obtained integer variables as be \(\begin{align}\tilde{\mathbf{b}}\end{align}\). Then similar to Algorithm 1, we apply SCA algorithm to solve (18) given be \(\begin{align}\tilde{\mathbf{b}}\end{align}\), and the obtained lower bound is denoted as L* = LB(ξ, \(\begin{align}\tilde{\mathbf{b}}\end{align}\)). If it is infeasible, the lower bound in subset Qξ is L* = ∞.

The branching and bounding steps are performed iteratively until the overall upper bound and overall lower bound converge. The details of the proposed branch and bound algorithm for solving problem (18) is summarized in Algorithm 2.

Algorithm 2 Algorithm Based on Branch and Bound

JAKO202305234914567_algor2.png 이미지

4.3 Low-complexity Suboptimal Deflation Algorithm Based on SCA

In this subsection, we propose a low-complexity algorithm for solving (18) in which the integer variables b are treated by the deflation process [14]. At first, all Ntot antennas are assumed to be selected. Next, we drop some poor communication link one by one until there remains NSM antennas, i.e., ║b║0 = NSM. In each step of the deflation process, after the remaining set of antennas is determined, i.e., giving b , the precoding vector v is obtained by solving (18) using the SCA framework similar to Algorithm 1.

Then the approximated convex problem at SCA iteration l can be formulated as

\(\begin{align}\max _{\mathbf{v}, t} t\end{align}\)       (21a)

s.t. (16c) - (16e), (19)       (21b)

The overall procedure of deflation SCA algorithm for solving problem (18) is summarized in Algorithm 3.

Algorithm 3 Deflation Algorithm Based on SCA

JAKO202305234914567_algor3.png 이미지

4.4 Complexity Analysis

In this section, the complexity of the proposed algorithms is discussed. Firstly, in Algorithm 3, the number of the subsets are needed to be considered increases exponentially with the problem dimension. Problem (20) approximately consist of 2Ntot + 1 variables and 1 second-order cone constraint of dimension Ntot . Hence, In the worst case, we need to traverse all 2Ntot − 2 subsets and the worst-case complexity is given by O(wNtot+1N2tot). Secondly, the iteration times in Algorithm 3 is Ntot - NSM, and the constraints in problem (21) consist of Ntot +1 variables and 1 second-order cone constraint of dimension Ntot approximately. Thus, the complexity of Algorithm 3 can be given by O(N3tot).

5. Numerical results

In this section, we present numerical results to evaluate the performance of the proposed algorithms, where K = 2, NR = 2, NSM = 2, ε = 10−4 and PT = NSM Watt. Moreover, QPSK constellations are employed to modulate the source information bits.

Firstly, we investigate the convergence behaviors of the branch and bound based algorithm in Algorithm 2 and the SCA algorithm in Algorithm 3 as depicted in Fig. 2 and Fig. 3, respectively. The parameters are set as Ntot = 6 and SNR = 3dB. The results in Fig. 2 corroborate the fact that the proposed branch and bound algorithm converges to its optimal value exactly. Note that in the first three iterations, the value of the lower bound is ∞, so it is not drawn in Fig. 2. Then, as the number of iterations increases, the upper bound decreases and the lower bound increases until convergence. Moreover, as shown in Fig. 3, we can see that the SCA algorithm converges fast toward the optimal value illustrated by the dotted line therein.

Next, we investigate the BER comparison versus the signal-to-noise-ratio (SNR) in Fig. 4 and Fig. 5 with Ntot = 4 or 6. As shown in Fig. 4 and Fig. 5, the branch and bound algorithm can achieve the best BER performance among all the algorithms since it can get the optimal solution, and the superiority becomes more pronounced as Ntot increases. Moreover, the performance of the deflation SCA algorithm is close to the exhaustive search based algorithm while avoiding enumerating all the antenna combinations. Furthermore, the BER performance of these algorithms with Ntot = 6 in Fig. 4 is better than that with Ntot = 4 in Fig. 5.

E1KOBZ_2023_v17n11_3204_f0002.png 이미지

Fig. 2. Convergence behavior of Algorithm 2 giving Ntot = 6, SNR = 3dB with QPSK.

E1KOBZ_2023_v17n11_3204_f0003.png 이미지

Fig. 3. Convergence behavior of SCA in Algorithm 3 giving Ntot = 6 , SNR = 3dB with QPSK.

E1KOBZ_2023_v17n11_3204_f0004.png 이미지

Fig. 4. BER for SMS giving NSM = 2 with QPSK, Ntot = 4 .

E1KOBZ_2023_v17n11_3204_f0005.png 이미지

Fig. 5. BER for SMS giving NSM = 2 with QPSK, Ntot = 4 .

6. Conclusion

In this paper, the transmit antenna selection and multicast precoding problem has been treated in multicast based SMSs. To this end, an equivalent single-stage reformulation of the original problem has been proposed to leverage the tools of optimization to design antenna selection and precoding coefficients jointly and simultaneously. Then, two algorithms to solve the reformulated mixed-integer non-linear problem have been proposed, where a branch and bound based iterative algorithm that can achieve quasi-optimal solution, and a deflation SCA algorithm that achieves elegant performance with low-complexity. The outcomes of this study highlight the significance of our approaches in improving the performance of transmit antenna selection and multicast precoding in SMSs. Numerical results demonstrate the convergence and efficiency of the proposed schemes, the convergence of our proposed iterative algorithm is between 10 and 15 iterations. Moreover, when the iterative algorithm and the exhaustive search based algorithm take the same BER accuracy conditions, the signal-to-noise-ratio (SNR) of the iterative algorithm is 1-2dB lower than the exhaustive search based algorithm.

References

  1. P. Yang, M. Di Renzo, Y. Xiao, S. Li and L. Hanzo, "Design Guidelines for Spatial Modulation," IEEE Communications Surveys & Tutorials, vol. 17, no. 1, pp. 6-26, Firstquarter 2015. https://doi.org/10.1109/COMST.2014.2327066
  2. N. D. Sidiropoulos, T. N. Davidson and Zhi-Quan Luo, "Transmit beamforming for physical-layer multicasting," IEEE Transactions on Signal Processing, vol. 54, no. 6, pp. 2239-2251, June 2006. https://doi.org/10.1109/TSP.2006.872578
  3. R. Rajashekar, K. V. S. Hari and L. Hanzo, "Quantifying the Transmit Diversity Order of Euclidean Distance Based Antenna Selection in Spatial Modulation," IEEE Signal Processing Letters, vol. 22, no. 9, pp. 1434-1437, Sept. 2015. https://doi.org/10.1109/LSP.2015.2408574
  4. R. Rajashekar, K. V. S. Hari and L. Hanzo, "Antenna Selection in Spatial Modulation Systems," IEEE Communications Letters, vol. 17, no. 3, pp. 521-524, March 2013. https://doi.org/10.1109/LCOMM.2013.012213.122650
  5. P. Cheng, Z. Chen, J. A. Zhang, Y. Li and B. Vucetic, "A Unified Precoding Scheme for Generalized Spatial Modulation," IEEE Transactions on Communications, vol. 66, no. 6, pp. 2502- 2514, June 2018. https://doi.org/10.1109/TCOMM.2018.2796605
  6. M. Lee, W. Chung, and T. Lee, "Generalized Precoder Design Formulation and Iterative Algorithm for Spatial Modulation in MIMO Systems With CSIT," IEEE Transactions on Communications, vol. 63, no. 4, pp. 1230-1244, April 2015. https://doi.org/10.1109/TCOMM.2015.2396521
  7. P. Yang, Y. L. Guan, Y. Xiao, M. D. Renzo, S. Li and L. Hanzo, "Transmit Precoded Spatial Modulation: Maximizing the Minimum Euclidean Distance Versus Minimizing the Bit Error Ratio," IEEE Transactions on Wireless Communications, vol. 15, no. 3, pp. 2054-2068, March 2016. https://doi.org/10.1109/TWC.2015.2497692
  8. M. Carosino and J. A. Ritcey, "Optimizing joint precoding and transmit antenna selection for spatial modulation," in Proc. of 2016 50th Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, pp. 915-919, 2016.
  9. M. -C. Lee, W. -H. Chung and T. -S. Lee, "BER Analysis for Spatial Modulation in Multicast MIMO Systems," IEEE Transactions on Communications, vol. 64, no. 7, pp. 2939-2951, July 2016. https://doi.org/10.1109/TCOMM.2016.2555912
  10. M. -C Lee and W. -H. Chung, "Configuration selection and precoder design for spatial modulation in multicast MIMO systems," in Proc. of 2015 IEEE 26th Annual International Symposium on Personal, Indoor, and Mobile Radio Communications (PIMRC), Hong Kong, China, pp. 45-50, 2015.
  11. P. K. Korrai, E. Lagunas, A. Bandi, S. K. Sharma and S. Chatzinotas, "Joint Power and Resource Block Allocation for Mixed-Numerology-Based 5G Downlink Under Imperfect CSI," IEEE Open Journal of the Communications Society, vol. 1, pp. 1583-1601, 2020. https://doi.org/10.1109/OJCOMS.2020.3029553
  12. L. B. Stotts, "Mathematical Preliminaries," in Free Space Optical Systems Engineering: Design and Analysis, Wiley, pp.1-49, 2017.
  13. A. Philipp, S. Ulbrich, Y. Cheng and M. Pesavento, "Multiuser downlink beamforming with interference cancellation using a SDP-based branch-and-bound algorithm," in Proc. of 2014 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Florence, Italy, pp. 7724-7728, 2014.
  14. Y.-F. Liu, Y.-H. Dai, and Z.-Q. Luo, "Joint Power and Admission Control via Linear Programming Deflation," IEEE Transactions on Signal Processing, vol. 61, no. 6, pp. 1327-1338, March15, 2013. https://doi.org/10.1109/TSP.2012.2236319