# Joint Power Allocation at the Base Station and the Relay for Untrusted Relay Cooperation OFDMA Network

Authors: Weiheng Jiang; Wenjiang Feng
Journal: International Journal of Antennas and Propagation
 (2015)

---

## Abstract

The secure communication that multiple OFDMA-based cell-edge mobile stations (MS) can only transmit confidential messages to base station (BS) through an untrusted intermediate relay (UR) is discussed. Specifically, with the destination-based jamming (DBJ) scheme and fixed MS transmission power assumption, our focus is on the joint BS and US power allocation to maximize system sum secrecy rate. We
first analyze the challenges in solving this problem. The result indicates that our nonconvex joint power allocation is equivalent to a joint MS access control and power allocation. Then, by problem relaxation and the alternating optimization approach, two suboptimal joint MS access control and power allocation algorithms are proposed. These algorithms alternatively solve the subproblem of joint BS and UR power allocation and the subproblem of MS selection until system sum secrecy rate is nonincreasing. In addition, the convergence and computational complexity of the proposed algorithms are analyzed. Finally, simulations results are presented to demonstrate the performance of our proposed algorithms.


---

## Body

## 1. Introduction

Broadcast is a fundamental property of the wireless medium. This property makes wireless communication susceptible to eavesdropping. Traditionally, this problem is typically addressed via upper layer approaches, such as the cryptographic protocols in the application layer which relies on computational complexity. Nowadays, information-theoretic security which exploits the properties of wireless channel to secure communications has received considerable attention.In fact, the secret communication in the presence of an eavesdropper was first introduced and studied by Wyner who considered a wiretap channel model [1]. Wyner showed that when eavesdropper’s channel is a degraded version of the main channel, the source and destination can achieve a positive perfect information rate (defined as secrecy rate). The maximal secrecy rate from the source to the destination is defined as the secrecy capacity and for the degraded wiretap channel is given as the difference between the rate at the legitimate receiver and the rate at the eavesdropper. Inspired by the pioneering work of Wyner, Csiszar and Korner then considered the general broadcast wiretap channel [2]. Their research showed that even when the eavesdropper is not degraded with respect to the legitimate user, secure communication between the legitimate users is possible by exploiting the inherent randomness of the communication channel. Recently, the secret communication in wireless networks has been intensively studied for various scenarios [3–5].Although the property of broadcasting makes wireless communication susceptible to eavesdropping, it also provides us with the opportunity to improve the security of wireless transmission by cooperation and relay [6, 7]. In recent years, the relay channel without security constraints has been studied under various scenarios. In most of these works, cooperation strategies such as decode-and-forward (DF), compress-and-forward (CF), and amplify-and-forward (AF) have been constructed to increase the transmission rate or reliability function [7, 8]. By using cooperative relay to enhance communication security, the research can be grouped into two types. The first type corresponds to the classical sense of cooperation, where the cooperating nodes strengthen the main transmission (from legitimate transmitter to the legitimate receiver) by using common relaying techniques such as AF and DF [6, 9–12]. For the strategies of the second type, the cooperating parties improve the secrecy performance of the system by weakening the eavesdropping link (from legitimate transmitter to the eavesdropper), such as the noise-forwarding [13], cooperative jamming [10–12, 14], and artificial noise [15]. For these works of cooperation enhanced secure communication, an important assumption is that the cooperation node is trustworthy. In other words, the cooperation (or relay) node is not an eavesdropper and does not make malicious attack. As pointed out in [16–18], the assumption that the relay node is trustworthy is not always true in some sense, since the relay node may not belong to the same party as the source or the destination. Under this situation, we may wonder whether secure communication between the source and the destination is still possible with the cooperation of untrusted node, especially in the scenario that no direct link can be established between the source and the destination. In fact, this involves untrusted relay (UR) cooperation secure communication which has already been discussed in [16–28].The untrusted relay model was first studied in [16] for the general relay channel. The secrecy capacity of this system will be zero if the relay channel is degraded and will equal the wiretap channel capacity if the channel is reversely degraded. Although it is a pessimistic result, one might wonder whether there exists any situation where the cooperation of UR can enable a higher secrecy rate than simply treating it as an eavesdropper. The positive response came from He and Yener [17]. Their research showed that using an CF-based UR to relay information can surely achieve a higher secrecy rate than just treating the relay as an eavesdropper. Later in [18], with the assumption that there is no direct link between the source and the destination, the destination-based jamming (DBJ) scheme was proposed to obtain positive secrecy rate even though the relay is untrusted. Following these information-theoretical developments, the joint source and relay beamforming design problem for an untrusted MIMO relay channel was studied in [19] and the secrecy outage probability for untrusted AF relay channels was investigated in [20]. The work in [21] considered secrecy rate maximization problem in untrusted two-way relaying channels with friendly jammers. Then game theory-based distributed joint jammer selection and power allocation algorithm was presented therein. In [22], the authors studied joint transmit design and relay node selection for a relay network with a group of untrusted relay nodes. Reference [23] presented a comprehensive analysis of the secrecy capacity for two-hop AF cooperative systems with untrustworthy relay nodes. The work in [24] considered the optimal source and relay transmission design for the hybrid network with both internal and external wiretappers. The optimal signal transmission power is calculated therein. The untrusted relay cooperation scenario also had been discussed in the cognitive radio networks [25]. The authors considered that the secondary users (SUs) helped the primary users (PUs) to relay their confidential messages in reward for being allowed to share their spectrum bands. However, the PUs might be reluctant to accept the SUs’ help, since the SUs are untrustworthy and may try unauthorized decoding of the PUs’ messages. Further research about this topic could be found in [26–28].All the above-mentioned literatures [16–28] only considered the secure communication between one pair of nodes by the cooperation of untrusted relay. Our work in this paper extends these works to the multiuser case. In particular, multiple OFDMA-based cell-edge mobile stations (MS) transmit confidential messages to base station (BS) only through an untrusted intermediate relay (UR). To be specific, with the destination-based jamming (DBJ) scheme and fixed MS transmission power assumption, our focus is on the joint BS and US power allocation to maximize system sum secrecy rate. Firstly, we discuss the challenges in solving this problem. The result indicates that our nonconvex joint power allocation is equivalent to a joint MS access control and power allocation. Secondly, by problem relaxation and the alternating optimization approach, two suboptimal joint MS access control and power allocation algorithms are proposed. These algorithms alternatively solve the subproblem of joint BS and UR power allocation and the subproblem of MS selection until system sum secrecy rate is nonincreasing. Also, the convergence and computational complexity of the proposed algorithms are analyzed. Finally, simulations are conducted to evaluate the performance of our proposed algorithms. The following of the paper is organized as follows: system model and cooperation protocol are introduced in Section 2; the joint power allocation problem is discussed in Section 3 and suboptimal algorithms are also present therein; selected numerical results are shown in Section 4; and finally the conclusion is given.


## 2. System Model and Cooperation Protocol

### 2.1. System Model

As depicted in Figure1, the system investigated herein is about a single cell network with one BS and multiple cell-edge MS which have secure communication requirements with BS in the uplink. Define the MS set by N and N=|N|, and MS are indexed by i=1,…,N. Each of the MS in N has already been allocated one orthogonal channel (in frequency-domain with unit bandwidth, such as OFDMA-based broadband wireless network). Since these MS are all far away from BS, no direct link between MS and BS can be established. Suppose that another MS (not in N) or heterogeneous network node is happening to be in the right place and has the will of cooperating with these MS as a volunteer of relay. (We do not discuss the motivations of UR cooperation and only focus on the joint power optimization. As in [21], the incentives of UR cooperation can be analyzed by game theory.) Unfortunately, this (relay) node (denoted by UR) may not belong to the same operator as BS and MS and thus has lower security clearance. Therefore, the UR has a dual role and acts as both an essential helper and a potential eavesdropper but not making any malicious attack. Under this setup, effective scheme should be introduced to complete the confidential message transmission and also to protect them not to be wiretapped by the UR. In this paper, there are no extra jammers existing and the DBJ scheme [18] is adopted to achieve the goal; that is, the BS plays the role of jammer. In addition, the following assumptions are used in our paper: (i) the cooperation strategy of UR is AF-based and full channel state information (CSI) is available at all nodes; (ii) BS, MS, and UR are all equipped with single antenna and operated in half-duplex mode; (iii) all MS transmit with fixed power; that is, piS(>0), ∀i∈N. The system aims at maximizing system sum secrecy rate by joint UR and BS power allocation. (Although the joint optimization of BS, UR, and MS power allocation can further improve the system sum secrecy rate, for simplicity, we only consider the problem of joint BS and UR power allocation herein.)Figure 1
System model.


### 2.2. Cooperation Protocol

In our considered scenario, the MS-BS secure communication process is performed by two phases. During the first phase (Phase I, shown with solid lines in Figure1), MS send their confidential messages to UR over their allocated channels. Simultaneously, the BS injects interference signals in the form of Gaussian noise over multiple MS working channels. Thus, the received signal at UR over channel i (for MS-i) is(1)yiR=piShiSRxiS+piBgiBRxiB+niR,where piB is the BS jamming signal power over channel-i. xiS and xiB (both with unit energy) denote the transmit modulated symbol by MS-i and the jamming signal by BS over channel i, respectively. Suppose that all channels are reciprocal; we use hiSR and giRB (or giBR) to denote the channel coefficients from MS-i to UR and from UR to BS over channel i, respectively. niR stands for a complex valued additive white Gaussian noise (AWGN) at UR in channel i such that niR~CN(0,σi2).During the secondary phase (Phase II, shown with dashed line in Figure1), UR normalizes its received signals over different channels and forwards the scaled versions to BS. In channel i, the signal regenerated by UR is xiR=ζipiRyiR, where ζi=(piS|hiSR|2+piB|giBR|2+σi2)-1/2 is the normalized factor at UR for MS-i and piR is the power allocated by UR for MS-i. Thus, the signal received at BS for MS-i is yiB=ζipiRgiRByiR+niB, ∀i∈N, where niB is the AWGN at BS over the channel for MS-i and niB~CN(0,σi2). Then, after canceling the self-interference, that is, ζipiRpiBgiBRgiRBxiB, the instantaneous received signal at the BS over channel i can be written as(2)yiB=piSpiRhiSRgiRBpiShiSR2+piBgiBR2+σi2xiS+piRgiRBpiShiSR2+piBgiBR2+σi2niR+niB.Let |hiSR|2/σi2=γiSR, |giBR|2/σi2=γiBR, and |giRB|2/σi2=γiRB (since giRB=giBR, we have γiBR=γiRB). Then, the instantaneous received signal SNR (signal to noise ratio) at the UR (which is denoted by χiR) and at the BS (which is denoted by χiB) for MS-i can be expressed as, respectively,(3)χiR=piSγiSR1+piBγiBR,χiB=piSpiRγiSRγiRB1+piSγiSR+piRγiRB+piBγiBR.


## 2.1. System Model

As depicted in Figure1, the system investigated herein is about a single cell network with one BS and multiple cell-edge MS which have secure communication requirements with BS in the uplink. Define the MS set by N and N=|N|, and MS are indexed by i=1,…,N. Each of the MS in N has already been allocated one orthogonal channel (in frequency-domain with unit bandwidth, such as OFDMA-based broadband wireless network). Since these MS are all far away from BS, no direct link between MS and BS can be established. Suppose that another MS (not in N) or heterogeneous network node is happening to be in the right place and has the will of cooperating with these MS as a volunteer of relay. (We do not discuss the motivations of UR cooperation and only focus on the joint power optimization. As in [21], the incentives of UR cooperation can be analyzed by game theory.) Unfortunately, this (relay) node (denoted by UR) may not belong to the same operator as BS and MS and thus has lower security clearance. Therefore, the UR has a dual role and acts as both an essential helper and a potential eavesdropper but not making any malicious attack. Under this setup, effective scheme should be introduced to complete the confidential message transmission and also to protect them not to be wiretapped by the UR. In this paper, there are no extra jammers existing and the DBJ scheme [18] is adopted to achieve the goal; that is, the BS plays the role of jammer. In addition, the following assumptions are used in our paper: (i) the cooperation strategy of UR is AF-based and full channel state information (CSI) is available at all nodes; (ii) BS, MS, and UR are all equipped with single antenna and operated in half-duplex mode; (iii) all MS transmit with fixed power; that is, piS(>0), ∀i∈N. The system aims at maximizing system sum secrecy rate by joint UR and BS power allocation. (Although the joint optimization of BS, UR, and MS power allocation can further improve the system sum secrecy rate, for simplicity, we only consider the problem of joint BS and UR power allocation herein.)Figure 1
System model.


## 2.2. Cooperation Protocol

In our considered scenario, the MS-BS secure communication process is performed by two phases. During the first phase (Phase I, shown with solid lines in Figure1), MS send their confidential messages to UR over their allocated channels. Simultaneously, the BS injects interference signals in the form of Gaussian noise over multiple MS working channels. Thus, the received signal at UR over channel i (for MS-i) is(1)yiR=piShiSRxiS+piBgiBRxiB+niR,where piB is the BS jamming signal power over channel-i. xiS and xiB (both with unit energy) denote the transmit modulated symbol by MS-i and the jamming signal by BS over channel i, respectively. Suppose that all channels are reciprocal; we use hiSR and giRB (or giBR) to denote the channel coefficients from MS-i to UR and from UR to BS over channel i, respectively. niR stands for a complex valued additive white Gaussian noise (AWGN) at UR in channel i such that niR~CN(0,σi2).During the secondary phase (Phase II, shown with dashed line in Figure1), UR normalizes its received signals over different channels and forwards the scaled versions to BS. In channel i, the signal regenerated by UR is xiR=ζipiRyiR, where ζi=(piS|hiSR|2+piB|giBR|2+σi2)-1/2 is the normalized factor at UR for MS-i and piR is the power allocated by UR for MS-i. Thus, the signal received at BS for MS-i is yiB=ζipiRgiRByiR+niB, ∀i∈N, where niB is the AWGN at BS over the channel for MS-i and niB~CN(0,σi2). Then, after canceling the self-interference, that is, ζipiRpiBgiBRgiRBxiB, the instantaneous received signal at the BS over channel i can be written as(2)yiB=piSpiRhiSRgiRBpiShiSR2+piBgiBR2+σi2xiS+piRgiRBpiShiSR2+piBgiBR2+σi2niR+niB.Let |hiSR|2/σi2=γiSR, |giBR|2/σi2=γiBR, and |giRB|2/σi2=γiRB (since giRB=giBR, we have γiBR=γiRB). Then, the instantaneous received signal SNR (signal to noise ratio) at the UR (which is denoted by χiR) and at the BS (which is denoted by χiB) for MS-i can be expressed as, respectively,(3)χiR=piSγiSR1+piBγiBR,χiB=piSpiRγiSRγiRB1+piSγiSR+piRγiRB+piBγiBR.


## 3. Jammer and Relay Power Allocation under Individual Power Constraint

In this section, the problem of system sum secrecy rate maximization is first introduced and we analyze the challenges in solving this problem. Second, the alternating optimization (AO) approach is provided to handle our problem and then an AO algorithm is given. The convergence and computational complexity of the AO algorithm are also analyzed. Finally, another suboptimal algorithm with lower complexity is presented.


### 3.1. Secrecy Rate Maximization Problem

As mentioned earlier, in our considered system, BS plays the role of jammer and relay is untrusted. The system aims to maximize sum secrecy rate by joint BS and UR power allocation. From [18, 21], we know that the achievable secrecy rate for MS-i at the BS is(4)RipiB,piR=12RiBpiB,piR-RiRpiB+,where the factor 1/2 is due to the fact that the communication is divided into two phases and we will omit it below. [·]+≜max⁡{0,·}, RiB(piB,piR)=log⁡(1+χiB), and RiR(piB)=log⁡(1+χiR) represent the capacity between MS-i and BS and between MS-i and UR, without the secrecy constraint, respectively. Based on the formulation of χiR and χiB in (3), it is obvious that the obtained secrecy rate for MS-i depends on both relay power piR and jammer power piB. Therefore, the joint optimization of BS and UR power allocation to maximize sum secrecy rate is required. Since the UR and BS have independently power constraint, consequently, the optimization problem (OP1) that we seek to solve is(5a)OP1:max⁡(pR,pB)R(pR,pB)=∑i∈N‍Ri(piB,piR)(5b)s.t.0⩽∑i∈N‍piR⩽PR,0⩽∑i∈N‍piB⩽PB(5c)s.t.IpiR>0=IpiB>0,∀i∈N,where pR=(p1R,…,pNR) and pB=(p1B,…,pNB) denote the power allocation vectors at UR and BS for the MS in N, respectively. PR and PB are the corresponding power constraints at UR and BS, respectively. I{X} is an indicator function such that I{X}=1 if the event X is true and I{X}=0 otherwise. Constraint (5c) comes from the fact that the direct link between BS and MS does not exist. Therefore, zero jamming or relay power always leads to zero secrecy rate.With the problem definition in (5a), (5b), and (5c), it is noted that three factors make OP1 hard to deal with: (i) the nonsmooth property of Ri(piB,piR) which is caused by the nonnegative constraint of the secrecy rate definition; (ii) the nonconvex property of Ri(piB,piR) (with respect to piB and piR), even without considering the nonnegative constraint of the secrecy rate definition; (iii) the potential combinatorial property of OP1 which is caused by the coupled relationship between piB and piR defined in (5c). In order to handle OP1, an important lemma is introduced first as below.Lemma 1.
For MS-i, ∀i∈N, with transmit power piS and secrecy rate definition (4), the sufficient and necessary condition for MS-i to obtain positive secrecy rate is that the allocated jamming power piB and relay power piR satisfy(6)piR>pi,thR=1+1+piSγiSR/piBγiRBγiRB,piB>0.Also, the necessary condition for all MS in N to obtain positive secrecy rate is that the UR power PR should satisfy (7). (One can note that the lower bound of condition (7) is loose, because we have assumed that the allocated jamming power is piB=∞ for ∀i∈N. In fact, we can attain more tightly conditions by assuming that, for ∀i∈N, piB=PB. Actually, this is unnecessary. Our purpose is to declare that OP1 is equivalent to a problem of joint MS access control and power allocation.) Consider(7)PR>∑i∈N‍1γiRB.Proof.
Please see AppendixA.By Lemma1, the exact nonlinear coupled relationship between piB and piR which is caused by nonnegative secrecy rate definition is characterized by (6). For (7), it is worth strengthening that it is just the necessary condition for the system to promise all MS in N to achieve positive secrecy rate but not sufficient condition. Thus, if (7) is false, then surely some MS in N will be rejected by the system due to overload; while if it is true, then it is still possible for some MS in N being rejected. Therefore, in essence, OP1 is equivalent to the following joint MS access control and power allocation problem (which is denoted by OP2):(8a)OP2:max⁡N∘⊆N,(pR,pB)R∘(pR,pB)=∑i∈N∘‍Ri∘(piB,piR)(8b)s.t.∑i∈N∘‍piR⩽PR,∑i∈N∘‍piB⩽PB(8c)s.t.IpiR>pi,thR=IpiB>0,∀i∈N∘, where N∘ is a MS subset. This subset is determined by the MS access control so that all MS in this subset can achieve positive secrecy rate. Define N∘=|N∘|, and let pR=(p1R,…,pN∘R) and pB=(p1B,…,pN∘B) denote the power allocation vectors by UR and BS for MS in N∘, respectively. (Although in OP1 we have used pR=(p1R,…,pNR) and pB=(p1B,…,pNB) to denote the power allocation vectors at UR and BS for the MS in N, resp., for notational convenience, we use pR=(p1R,…,pN∘R) and pB=(p1B,…,pN∘B) below to denote the power allocation vectors at UR and BS for the MS in N∘, resp.) In (8a), we define Ri∘(piB,piR)=RiB(piB,piR)-RiR(piB). Thus, Ri∘(piB,piR) is equivalent to Ri(piB,piR) under the condition that (6) (or (8c)) is satisfied. Constraint (8c) comes from (6) and (5c), in which pi,thR is a function of piB.For OP2, constraint (8c) and MS access control requirements make OP2 still a mixed integer nonlinear programming (MINLP) and thus intractable. However, from OP2, a relaxed joint MS access control and power allocation problem can be attained as below (which is denoted by OP3): (9)OP3:max⁡N∘⊆Nmax⁡(pR,pB)R∘(pR,pB)=∑i∈N∘‍Ri∘(piB,piR)s.t.∑i∈N∘‍piR⩽PR,∑i∈N∘‍piB⩽PB,where the constraint (8c) is removed from OP2.By observing OP3, one can find that it has three important properties: (i) OP3 can be decomposed into two subproblems: the inner joint power allocation and the outer MS access selection; (ii) the inner subproblem has convex constraint set; (iii) the optimal solution of OP3 and OP2 has the following relationship.Proposition 2.
LetSOP3={NOP3⋆,pi,OP3⋆R,pi,OP3⋆B} be the optimal solution for OP3; then it is also the optimal solution for OP2. (Since both OP3 and OP2 are combinatorial-based problems, their optimal solution may be not unique. However, Proposition 2 is still valid for the case of more than one optimal solution.)Proof.
Please see AppendixB.Due to the properties of OP3, solving OP3 instead of OP2 would be more sensible. However, there are four challenges should be overcome: (i) OP3 is still a MINLP which is known to be nondeterministic polynomial-time hard (NP-hard); (ii) although OP3 could be handled by alternatively solving the outer and inner subproblem, the attained solution may be outside the feasible region of OP2 (because we have removed constraint (8c) to get OP3); (iii) the objective function of the inner subproblem is not jointly concave with decision variables; (iv) the outer subproblem of OP3 is combinatorial-based and we cannot construct rules about the optimal N∘. In this paper, we provide a suboptimal scheme to explore the solution for OP3. In addition, it is proved that the obtained solution is at least suboptimality for OP2. To be specific, our suboptimal scheme is as follows: (i) we alternatively solve the inner and outer subproblem of OP3; (ii) in handling the inner subproblem, alternative optimization could be used as well. Although the objective function of OP3 is nonconcave, BS power optimization problem has unique solution with previously given UR power allocation; also, the optimal UR power allocation problem is convex for fixed BS power allocation; (iii) the outer subproblem of OP3 can be handled by a heuristic scheme. Based on the above analysis, in the following, our focus is on handling OP3 by AO approach.


### 3.2. AO Approach and the Algorithm

In this subsection, we first assume that the MS access set has been predetermined and alternatively discuss the relay and jammer power allocation. Then, a suboptimal MS access control scheme is proposed.


#### 3.2.1. Optimal Relay Power Allocation under Fixed Jamming

We first discuss relay power allocation under fixed jamming powerp⋆B and predetermined MS set N∘. Then, we have the following optimization problem OP4:(10a)OP4:max⁡pR∑i∈N∘‍Ri∘(pi⋆B,piR)(10b)s.t.∑i∈N∘‍piR⩽PR,fixedp⋆B,N∘.Since both N∘ and p⋆B are fixed and also constraint (10b) is convex, it follows from [29] that OP4 is a degraded AF-based relay network power allocation problem with parallel orthogonal Gaussian relay channels. Therefore, the optimal solution of OP4 can be characterized by the following theorem.Theorem 3.
For OP4, its optimal solution is as follows:(11)pi⋆R=Bi2-4AiCi-Bi2Ai+,∀i∈N∘,where Ai=(γiBR)2(1+piSγiSR), Bi=γiRB(2+piSγiSR)(1+piSγiSR+pi⋆BγiBR), Ci=(1+piSγiSR+pi⋆BγiBR)2-piSγiSRγiRB(1+piSγiSR+pi⋆BγiBR)/λ⋆, and λ⋆ takes the value such that ∑i∈N∘‍pi⋆R(λ⋆,pi⋆B)=PR.Proof.
Please see AppendixC.At the moment, one may note that, in the above relay power allocation, we have not used the conditions presented in Lemma1 to restrict the relay power allocation or to exclude MS which cannot promise positive secrecy rate. The explanation is as follows: the purpose of introducing Lemma 1 is to uncover the fact that the joint relay and jammer power allocation is essentially equivalent to the joint MS access control and power allocation. In that case, the access control rather than power allocation is used to exclude MS which obtains negative secrecy rate.


#### 3.2.2. Optimal Jamming under Fixed Relay Power Allocation

Now, our attention turns to the jamming power allocation at BS for givenp⋆R and N∘. Thus, we have the following optimization problem OP5:(12a)OP5:max⁡pB∑i∈N∘‍Ri∘(piB,pi⋆R)(12b)s.t.∑i∈N∘‍piB⩽PB,fixedp⋆R,N∘.It is obvious that the jamming power constraint (12b) is convex. Therefore, the focus is on the objective function (12a). In order to verify the convexity of the objective function, we take the derivative of ∑i∈N∘‍Ri∘(piB,pi⋆R) with respect to piB, ∀i∈N∘, and then have(13)∂∑i∈N∘‍Ri∘(piB,pi⋆R)∂piB=piSγiSRγiBRDipiB2+EipiB+FiOi,where Di=(1-pi⋆RγiRB)(γiBR)2, Ei=2γiBR(1+pi⋆RγiBR)>0, Fi=(1+piSγiSR)[(1+pi⋆RγiBR)(piSγiSR+pi⋆RγiRB)+1]>0, and Oi=[(1+piSγiSR+pi⋆RγiRB+piBγiBR)2+piSpi⋆RγiSRγiRB(1+piSγiSR+pi⋆RγiRB+piBγiBR)]×[(1+piBγiBR)2+piSγiSR(1+piBγiBR)]>0. It can be noted that if pi⋆RγiRB>1 (from Lemma 1, one may notice that it is the necessary condition for MS-i to obtain positive secrecy rate), then Ri∘(piB,pi⋆R) is a quasi-concave function of piB and takes the maximum at(14)pi,1⋆B=-Ei-Ei2-4DiFi2Di.Namely, Ri∘(piB,pi⋆R) is an increasing function of piB for piB∈0,pi,1⋆B and a decreasing function of piB for piB∈pi,1⋆B,∞, while if pi⋆RγiRB⩽1 (from Theorem 3, one can find that this is possible after the relay power allocation), then we have ∂Ri∘(piB,pi⋆R)/∂piB>0, ∂2Ri∘(piB,pi⋆R)/∂(piB)2<0, and ∂2Ri∘(piB,pi⋆R)/∂(piB)∂(pjB)=0, i≠j. Therefore, the objective function is a monotonically increasing and concave function with respect to piB. With the BS power constraint PB, it is known that Ri∘(piB,pi⋆R) will be maximized at piB=PB. Following these analyses, three cases may appear for OP5: (i) for ∀i∈N∘ we have pi⋆RγiRB>1 and then the objective function is quasi-concave and OP5 becomes nonconvex (but quasi-concave); (ii) for ∀i∈N∘ we have pi⋆RγiRB⩽1 and then the objective function is concave and OP5 is a convex programming; (iii) ∃i∈N∘ let pi⋆RγiRB>1, and also ∃j∈N∘ but  j≠i let pj⋆RγiRB⩽1; then OP5 is a hybrid quasi-concave and concave problem. Nevertheless, it is still possible to analyze and derive the optimal solution of OP5 by the KKT conditions [30]. In other words, the KKT conditions are the necessary conditions (not sufficient conditions due to may be nonconcavity of Ri∘(piB,pi⋆R) with respect to piB, ∀i∈N∘). Among the KKT conditions, we have(15)∂Ri∘piB,pi⋆R∂piB-ξ=0,∀i∈N∘,where ξ is the Lagrange multiplier associated with the power constraint (12b) and ξ⩾0. Based on (15), we use the following theorem to characterize the properties of the solution for OP5.Theorem 4.
For OP5 with givenp⋆R and N∘ and ∀i∈N∘, if pi⋆RγiRB>1, let pi,1⋆B be defined as (14); if pi⋆RγiRB⩽1, let pi,1⋆B=PB. Also let pi,2⋆B(ξ) be the real positive root of (15) for given ξ and define Xi(ξ)={pi,2⋆Bξ∣pi,2⋆B(ξ)>0,∂Ri∘(piB,pi⋆R)/∂piBpi,2⋆B(ξ)-ξ=0}. The unique solution of OP5 is denoted by ξ⋆,{pi⋆B:∀i∈N∘}, then we have the following: (1)
If∑i∈N∘‍pi,1⋆B⩽PB, then pi⋆B=pi,1⋆B and ξ⋆=0; otherwise, pi⋆B=p~i⋆BI{Xi(ξ⋆)≠Φ}+0×I{Xi(ξ⋆)=Φ}, where p~i⋆B=argmax⁡pi,2⋆B(ξ⋆)∈Xi(ξ⋆)pi,2⋆B(ξ⋆) and Φ is null set; ξ⋆ takes the value such that ∑i∈N∘‍pi⋆B(pi⋆R,ξ⋆)=PB.(2)
If∑i∈N∘‍pi,1⋆B>PB, then the term ∑i∈N∘‍pi⋆B(pi⋆R,ξ) is a decreasing function of ξ.Proof.
Please see AppendixD.By the results in Theorem4, for jamming power allocation, we first calculate the jamming power pi,1⋆B, ∀i∈N∘, for each accessed MS. If ∑i∈N∘‍pi,1⋆B(pi⋆R)⩽PB, that is, the sum of jamming power demand is less than or equal to the available power at BS, then all MS in N∘ will be allocated the jamming power pi⋆B=pi,1⋆B, ∀i∈N∘, while if ∑i∈N∘‍pi,1⋆B(pi⋆R)>PB, we should perform a bisection search to obtain the optimal Lagrange multiplier ξ⋆ and calculate the corresponding jamming power pi⋆B, ∀i∈N∘, by (15).


#### 3.2.3. Access Control Scheme and Alternating Optimization Algorithm

Based on the alternative optimization result, a suboptimal access control scheme which can find a feasible suboptimal solution for OP3 (or OP2) is provided herein. The main idea of this scheme is as follows: when the alternative relay and jamming power optimization is convergent, we remove the MS inN∘ which has obtained the minimal (general) secrecy rate and then repeat the procedure until the system sum (general) secrecy rate does not increase. Therefore, the proposed suboptimal MS access control scheme and the AO-based power allocation together solve OP3. Hence, we have the following Algorithm 5.Algorithm 5 (alternative relay power optimization and jamming power optimization (OptROptJam)).
Step 1. Initialize t=0 and let N∘(t)=N, pB(t)=p0B.
Step 2. Based on Theorems 3 and 4, alternatively optimize the relay and jamming power until convergence and obtain piB(t) and piR(t), respectively; then go to Step 3.
Step 3. Calculate the system sum (general) secrecy rate ∑i∈N∘(t)‍Ri∘(piB(t),piR(t)); if ∑i∈N∘(t)‍Ri∘(piB(t),piR(t))>∑i∈N∘(t-1)‍Ri∘(piB(t-1),piR(t-1)), then let N∘(t)=N∘(t)∖i⋆ and t=t+1; go to Step 2, where i⋆=argmin⁡i∈N∘(t)Ri∘(t); otherwise, the algorithm is finished.In Algorithm5, t represents the iteration time of the algorithm, and N∘(t) denotes the accessed MS subset at t. In this paper, we set N∘(0)=N; that is, all MS are accessed by the system at the beginning. pB(t)=p0B is the start point of the algorithm; that is, we set a jamming power vector and then perform relay power iteration. The rule of start point setting is that it should satisfy the corresponding power constraint. The convergence condition in Step 2 is that the difference of two rounds of alternative relay and jamming power optimization is less than a threshold ɛ. In Step 3, if there are more than one MS feedback from argmin⁡i∈N∘(t)Ri∘(t), then one of them is randomly removed.Suboptimality and Convergence. For Algorithm 5, its optimality is affected by two factors: (i) the optimality of the access control scheme and (ii) the performance of alternative power optimization. For the former, the rules of constructing the optimal MS access set are unclear. Thus, it is difficult to judge the optimality of the proposed suboptimal access control scheme. In this paper, we use simulation result to demonstrate its performance. For the alternating-based power optimization approach, in general, it may fail to locate the stationary points, not to mention the global convergence to the optimal solution. However, we have the following conclusion for Algorithm 5.Proposition 6.
Algorithm5 always has the limit point and any limit point generated by Algorithm 5 is at least a suboptimal solution for OP2.Proof.
Please see AppendixE.Complexity. The complexity of Algorithm 5 is mainly affected by the alternative optimization and the access control. The former involves searching two optimal Lagrange multipliers which have the approximate complexity of O(Nlog⁡2(λ⋆/ɛ)) and O(Nlog⁡2(ξ⋆/ɛ)), respectively, where λ⋆ and ξ⋆ are provided in Theorems 3 and 4, respectively, and ɛ is the tolerable error presented in Algorithm 5. The approximate complexity of the proposed access control scheme is O(N), where N is the number of MS in the system.


### 3.3. Another Suboptimal Algorithm

For two reasons, we propose another suboptimal algorithm for our problem: (i) pursuing lower complexity algorithm and (ii) as the benchmark scheme of Algorithm5. This suboptimal algorithm is based on two suboptimal power allocation strategies at UR and BS, respectively, that is, equal relay power allocation and optimal proportional jamming power allocation. More precisely, at iteration t, the relay and jamming power allocation rules are as follows, respectively:(16)piR(t)=PRN∘t(17)piBt=PBpi,1⋆BpiRt∑i∈N∘t‍pi,1⋆BpiRtI∑i∈N∘t‍pi,1⋆BpiRt>PB+pi,1⋆BpiRtI∑i∈N∘t‍pi,1⋆BpiRt⩽PB,where pi,1⋆B(piR(t)) is defined in Theorem 4. From (16) and (17), one can find that UR equally allocates its power to the accessed MS at t and the BS power allocation is based on proportional rule. In particular, if ∑i∈N∘(t)‍pi,1⋆B(piR(t))⩽PB, then the obtained jamming power for MS-i (i∈N∘(t)) is pi,1⋆B; otherwise, its jamming power is in proportion to pi,1⋆B(piR(t))/∑i∈N∘(t)‍pi,1⋆B(piR(t)). The access control scheme for this new suboptimal algorithm is the same as Algorithm 5. To sum up, we have another suboptimal algorithm (Algorithm 7) which is summarized as below.Algorithm 7 (equal relay power allocation and optimal proportional jamming power allocation (EROptPJam)).
Step 1. Initialize t=0 and N∘(t)=N.
Step 2. Perform relay and jamming power allocation based on (16) and (17), respectively; then we have piB(t) and piR(t) and go to Step 3.
Step 3. Calculate the system sum (general) secrecy rate ∑i∈N∘(t)‍Ri∘(piB(t),piR(t)) and if ∑i∈N∘(t)‍Ri∘(piB(t),piR(t))>∑i∈N∘(t-1)‍Ri∘(piB(t-1),piR(t-1)), then let N∘(t)=N∘(t)∖i⋆ and t=t+1; go to Step 2, where i⋆=argmin⁡i∈N∘(t)Ri∘(t); otherwise, the algorithm is finished.One can note that the complexity of Algorithm7 is significantly lower. This is because we do not need to perform any power iteration but just calculate the jamming power pi,1⋆B by (14) for given relay power allocation. Therefore, it is of no doubt that Algorithm 7 has lower complexity but experiences worse performance than Algorithm 5.


## 3.1. Secrecy Rate Maximization Problem

As mentioned earlier, in our considered system, BS plays the role of jammer and relay is untrusted. The system aims to maximize sum secrecy rate by joint BS and UR power allocation. From [18, 21], we know that the achievable secrecy rate for MS-i at the BS is(4)RipiB,piR=12RiBpiB,piR-RiRpiB+,where the factor 1/2 is due to the fact that the communication is divided into two phases and we will omit it below. [·]+≜max⁡{0,·}, RiB(piB,piR)=log⁡(1+χiB), and RiR(piB)=log⁡(1+χiR) represent the capacity between MS-i and BS and between MS-i and UR, without the secrecy constraint, respectively. Based on the formulation of χiR and χiB in (3), it is obvious that the obtained secrecy rate for MS-i depends on both relay power piR and jammer power piB. Therefore, the joint optimization of BS and UR power allocation to maximize sum secrecy rate is required. Since the UR and BS have independently power constraint, consequently, the optimization problem (OP1) that we seek to solve is(5a)OP1:max⁡(pR,pB)R(pR,pB)=∑i∈N‍Ri(piB,piR)(5b)s.t.0⩽∑i∈N‍piR⩽PR,0⩽∑i∈N‍piB⩽PB(5c)s.t.IpiR>0=IpiB>0,∀i∈N,where pR=(p1R,…,pNR) and pB=(p1B,…,pNB) denote the power allocation vectors at UR and BS for the MS in N, respectively. PR and PB are the corresponding power constraints at UR and BS, respectively. I{X} is an indicator function such that I{X}=1 if the event X is true and I{X}=0 otherwise. Constraint (5c) comes from the fact that the direct link between BS and MS does not exist. Therefore, zero jamming or relay power always leads to zero secrecy rate.With the problem definition in (5a), (5b), and (5c), it is noted that three factors make OP1 hard to deal with: (i) the nonsmooth property of Ri(piB,piR) which is caused by the nonnegative constraint of the secrecy rate definition; (ii) the nonconvex property of Ri(piB,piR) (with respect to piB and piR), even without considering the nonnegative constraint of the secrecy rate definition; (iii) the potential combinatorial property of OP1 which is caused by the coupled relationship between piB and piR defined in (5c). In order to handle OP1, an important lemma is introduced first as below.Lemma 1.
For MS-i, ∀i∈N, with transmit power piS and secrecy rate definition (4), the sufficient and necessary condition for MS-i to obtain positive secrecy rate is that the allocated jamming power piB and relay power piR satisfy(6)piR>pi,thR=1+1+piSγiSR/piBγiRBγiRB,piB>0.Also, the necessary condition for all MS in N to obtain positive secrecy rate is that the UR power PR should satisfy (7). (One can note that the lower bound of condition (7) is loose, because we have assumed that the allocated jamming power is piB=∞ for ∀i∈N. In fact, we can attain more tightly conditions by assuming that, for ∀i∈N, piB=PB. Actually, this is unnecessary. Our purpose is to declare that OP1 is equivalent to a problem of joint MS access control and power allocation.) Consider(7)PR>∑i∈N‍1γiRB.Proof.
Please see AppendixA.By Lemma1, the exact nonlinear coupled relationship between piB and piR which is caused by nonnegative secrecy rate definition is characterized by (6). For (7), it is worth strengthening that it is just the necessary condition for the system to promise all MS in N to achieve positive secrecy rate but not sufficient condition. Thus, if (7) is false, then surely some MS in N will be rejected by the system due to overload; while if it is true, then it is still possible for some MS in N being rejected. Therefore, in essence, OP1 is equivalent to the following joint MS access control and power allocation problem (which is denoted by OP2):(8a)OP2:max⁡N∘⊆N,(pR,pB)R∘(pR,pB)=∑i∈N∘‍Ri∘(piB,piR)(8b)s.t.∑i∈N∘‍piR⩽PR,∑i∈N∘‍piB⩽PB(8c)s.t.IpiR>pi,thR=IpiB>0,∀i∈N∘, where N∘ is a MS subset. This subset is determined by the MS access control so that all MS in this subset can achieve positive secrecy rate. Define N∘=|N∘|, and let pR=(p1R,…,pN∘R) and pB=(p1B,…,pN∘B) denote the power allocation vectors by UR and BS for MS in N∘, respectively. (Although in OP1 we have used pR=(p1R,…,pNR) and pB=(p1B,…,pNB) to denote the power allocation vectors at UR and BS for the MS in N, resp., for notational convenience, we use pR=(p1R,…,pN∘R) and pB=(p1B,…,pN∘B) below to denote the power allocation vectors at UR and BS for the MS in N∘, resp.) In (8a), we define Ri∘(piB,piR)=RiB(piB,piR)-RiR(piB). Thus, Ri∘(piB,piR) is equivalent to Ri(piB,piR) under the condition that (6) (or (8c)) is satisfied. Constraint (8c) comes from (6) and (5c), in which pi,thR is a function of piB.For OP2, constraint (8c) and MS access control requirements make OP2 still a mixed integer nonlinear programming (MINLP) and thus intractable. However, from OP2, a relaxed joint MS access control and power allocation problem can be attained as below (which is denoted by OP3): (9)OP3:max⁡N∘⊆Nmax⁡(pR,pB)R∘(pR,pB)=∑i∈N∘‍Ri∘(piB,piR)s.t.∑i∈N∘‍piR⩽PR,∑i∈N∘‍piB⩽PB,where the constraint (8c) is removed from OP2.By observing OP3, one can find that it has three important properties: (i) OP3 can be decomposed into two subproblems: the inner joint power allocation and the outer MS access selection; (ii) the inner subproblem has convex constraint set; (iii) the optimal solution of OP3 and OP2 has the following relationship.Proposition 2.
LetSOP3={NOP3⋆,pi,OP3⋆R,pi,OP3⋆B} be the optimal solution for OP3; then it is also the optimal solution for OP2. (Since both OP3 and OP2 are combinatorial-based problems, their optimal solution may be not unique. However, Proposition 2 is still valid for the case of more than one optimal solution.)Proof.
Please see AppendixB.Due to the properties of OP3, solving OP3 instead of OP2 would be more sensible. However, there are four challenges should be overcome: (i) OP3 is still a MINLP which is known to be nondeterministic polynomial-time hard (NP-hard); (ii) although OP3 could be handled by alternatively solving the outer and inner subproblem, the attained solution may be outside the feasible region of OP2 (because we have removed constraint (8c) to get OP3); (iii) the objective function of the inner subproblem is not jointly concave with decision variables; (iv) the outer subproblem of OP3 is combinatorial-based and we cannot construct rules about the optimal N∘. In this paper, we provide a suboptimal scheme to explore the solution for OP3. In addition, it is proved that the obtained solution is at least suboptimality for OP2. To be specific, our suboptimal scheme is as follows: (i) we alternatively solve the inner and outer subproblem of OP3; (ii) in handling the inner subproblem, alternative optimization could be used as well. Although the objective function of OP3 is nonconcave, BS power optimization problem has unique solution with previously given UR power allocation; also, the optimal UR power allocation problem is convex for fixed BS power allocation; (iii) the outer subproblem of OP3 can be handled by a heuristic scheme. Based on the above analysis, in the following, our focus is on handling OP3 by AO approach.


## 3.2. AO Approach and the Algorithm

In this subsection, we first assume that the MS access set has been predetermined and alternatively discuss the relay and jammer power allocation. Then, a suboptimal MS access control scheme is proposed.


### 3.2.1. Optimal Relay Power Allocation under Fixed Jamming

We first discuss relay power allocation under fixed jamming powerp⋆B and predetermined MS set N∘. Then, we have the following optimization problem OP4:(10a)OP4:max⁡pR∑i∈N∘‍Ri∘(pi⋆B,piR)(10b)s.t.∑i∈N∘‍piR⩽PR,fixedp⋆B,N∘.Since both N∘ and p⋆B are fixed and also constraint (10b) is convex, it follows from [29] that OP4 is a degraded AF-based relay network power allocation problem with parallel orthogonal Gaussian relay channels. Therefore, the optimal solution of OP4 can be characterized by the following theorem.Theorem 3.
For OP4, its optimal solution is as follows:(11)pi⋆R=Bi2-4AiCi-Bi2Ai+,∀i∈N∘,where Ai=(γiBR)2(1+piSγiSR), Bi=γiRB(2+piSγiSR)(1+piSγiSR+pi⋆BγiBR), Ci=(1+piSγiSR+pi⋆BγiBR)2-piSγiSRγiRB(1+piSγiSR+pi⋆BγiBR)/λ⋆, and λ⋆ takes the value such that ∑i∈N∘‍pi⋆R(λ⋆,pi⋆B)=PR.Proof.
Please see AppendixC.At the moment, one may note that, in the above relay power allocation, we have not used the conditions presented in Lemma1 to restrict the relay power allocation or to exclude MS which cannot promise positive secrecy rate. The explanation is as follows: the purpose of introducing Lemma 1 is to uncover the fact that the joint relay and jammer power allocation is essentially equivalent to the joint MS access control and power allocation. In that case, the access control rather than power allocation is used to exclude MS which obtains negative secrecy rate.


### 3.2.2. Optimal Jamming under Fixed Relay Power Allocation

Now, our attention turns to the jamming power allocation at BS for givenp⋆R and N∘. Thus, we have the following optimization problem OP5:(12a)OP5:max⁡pB∑i∈N∘‍Ri∘(piB,pi⋆R)(12b)s.t.∑i∈N∘‍piB⩽PB,fixedp⋆R,N∘.It is obvious that the jamming power constraint (12b) is convex. Therefore, the focus is on the objective function (12a). In order to verify the convexity of the objective function, we take the derivative of ∑i∈N∘‍Ri∘(piB,pi⋆R) with respect to piB, ∀i∈N∘, and then have(13)∂∑i∈N∘‍Ri∘(piB,pi⋆R)∂piB=piSγiSRγiBRDipiB2+EipiB+FiOi,where Di=(1-pi⋆RγiRB)(γiBR)2, Ei=2γiBR(1+pi⋆RγiBR)>0, Fi=(1+piSγiSR)[(1+pi⋆RγiBR)(piSγiSR+pi⋆RγiRB)+1]>0, and Oi=[(1+piSγiSR+pi⋆RγiRB+piBγiBR)2+piSpi⋆RγiSRγiRB(1+piSγiSR+pi⋆RγiRB+piBγiBR)]×[(1+piBγiBR)2+piSγiSR(1+piBγiBR)]>0. It can be noted that if pi⋆RγiRB>1 (from Lemma 1, one may notice that it is the necessary condition for MS-i to obtain positive secrecy rate), then Ri∘(piB,pi⋆R) is a quasi-concave function of piB and takes the maximum at(14)pi,1⋆B=-Ei-Ei2-4DiFi2Di.Namely, Ri∘(piB,pi⋆R) is an increasing function of piB for piB∈0,pi,1⋆B and a decreasing function of piB for piB∈pi,1⋆B,∞, while if pi⋆RγiRB⩽1 (from Theorem 3, one can find that this is possible after the relay power allocation), then we have ∂Ri∘(piB,pi⋆R)/∂piB>0, ∂2Ri∘(piB,pi⋆R)/∂(piB)2<0, and ∂2Ri∘(piB,pi⋆R)/∂(piB)∂(pjB)=0, i≠j. Therefore, the objective function is a monotonically increasing and concave function with respect to piB. With the BS power constraint PB, it is known that Ri∘(piB,pi⋆R) will be maximized at piB=PB. Following these analyses, three cases may appear for OP5: (i) for ∀i∈N∘ we have pi⋆RγiRB>1 and then the objective function is quasi-concave and OP5 becomes nonconvex (but quasi-concave); (ii) for ∀i∈N∘ we have pi⋆RγiRB⩽1 and then the objective function is concave and OP5 is a convex programming; (iii) ∃i∈N∘ let pi⋆RγiRB>1, and also ∃j∈N∘ but  j≠i let pj⋆RγiRB⩽1; then OP5 is a hybrid quasi-concave and concave problem. Nevertheless, it is still possible to analyze and derive the optimal solution of OP5 by the KKT conditions [30]. In other words, the KKT conditions are the necessary conditions (not sufficient conditions due to may be nonconcavity of Ri∘(piB,pi⋆R) with respect to piB, ∀i∈N∘). Among the KKT conditions, we have(15)∂Ri∘piB,pi⋆R∂piB-ξ=0,∀i∈N∘,where ξ is the Lagrange multiplier associated with the power constraint (12b) and ξ⩾0. Based on (15), we use the following theorem to characterize the properties of the solution for OP5.Theorem 4.
For OP5 with givenp⋆R and N∘ and ∀i∈N∘, if pi⋆RγiRB>1, let pi,1⋆B be defined as (14); if pi⋆RγiRB⩽1, let pi,1⋆B=PB. Also let pi,2⋆B(ξ) be the real positive root of (15) for given ξ and define Xi(ξ)={pi,2⋆Bξ∣pi,2⋆B(ξ)>0,∂Ri∘(piB,pi⋆R)/∂piBpi,2⋆B(ξ)-ξ=0}. The unique solution of OP5 is denoted by ξ⋆,{pi⋆B:∀i∈N∘}, then we have the following: (1)
If∑i∈N∘‍pi,1⋆B⩽PB, then pi⋆B=pi,1⋆B and ξ⋆=0; otherwise, pi⋆B=p~i⋆BI{Xi(ξ⋆)≠Φ}+0×I{Xi(ξ⋆)=Φ}, where p~i⋆B=argmax⁡pi,2⋆B(ξ⋆)∈Xi(ξ⋆)pi,2⋆B(ξ⋆) and Φ is null set; ξ⋆ takes the value such that ∑i∈N∘‍pi⋆B(pi⋆R,ξ⋆)=PB.(2)
If∑i∈N∘‍pi,1⋆B>PB, then the term ∑i∈N∘‍pi⋆B(pi⋆R,ξ) is a decreasing function of ξ.Proof.
Please see AppendixD.By the results in Theorem4, for jamming power allocation, we first calculate the jamming power pi,1⋆B, ∀i∈N∘, for each accessed MS. If ∑i∈N∘‍pi,1⋆B(pi⋆R)⩽PB, that is, the sum of jamming power demand is less than or equal to the available power at BS, then all MS in N∘ will be allocated the jamming power pi⋆B=pi,1⋆B, ∀i∈N∘, while if ∑i∈N∘‍pi,1⋆B(pi⋆R)>PB, we should perform a bisection search to obtain the optimal Lagrange multiplier ξ⋆ and calculate the corresponding jamming power pi⋆B, ∀i∈N∘, by (15).


### 3.2.3. Access Control Scheme and Alternating Optimization Algorithm

Based on the alternative optimization result, a suboptimal access control scheme which can find a feasible suboptimal solution for OP3 (or OP2) is provided herein. The main idea of this scheme is as follows: when the alternative relay and jamming power optimization is convergent, we remove the MS inN∘ which has obtained the minimal (general) secrecy rate and then repeat the procedure until the system sum (general) secrecy rate does not increase. Therefore, the proposed suboptimal MS access control scheme and the AO-based power allocation together solve OP3. Hence, we have the following Algorithm 5.Algorithm 5 (alternative relay power optimization and jamming power optimization (OptROptJam)).
Step 1. Initialize t=0 and let N∘(t)=N, pB(t)=p0B.
Step 2. Based on Theorems 3 and 4, alternatively optimize the relay and jamming power until convergence and obtain piB(t) and piR(t), respectively; then go to Step 3.
Step 3. Calculate the system sum (general) secrecy rate ∑i∈N∘(t)‍Ri∘(piB(t),piR(t)); if ∑i∈N∘(t)‍Ri∘(piB(t),piR(t))>∑i∈N∘(t-1)‍Ri∘(piB(t-1),piR(t-1)), then let N∘(t)=N∘(t)∖i⋆ and t=t+1; go to Step 2, where i⋆=argmin⁡i∈N∘(t)Ri∘(t); otherwise, the algorithm is finished.In Algorithm5, t represents the iteration time of the algorithm, and N∘(t) denotes the accessed MS subset at t. In this paper, we set N∘(0)=N; that is, all MS are accessed by the system at the beginning. pB(t)=p0B is the start point of the algorithm; that is, we set a jamming power vector and then perform relay power iteration. The rule of start point setting is that it should satisfy the corresponding power constraint. The convergence condition in Step 2 is that the difference of two rounds of alternative relay and jamming power optimization is less than a threshold ɛ. In Step 3, if there are more than one MS feedback from argmin⁡i∈N∘(t)Ri∘(t), then one of them is randomly removed.Suboptimality and Convergence. For Algorithm 5, its optimality is affected by two factors: (i) the optimality of the access control scheme and (ii) the performance of alternative power optimization. For the former, the rules of constructing the optimal MS access set are unclear. Thus, it is difficult to judge the optimality of the proposed suboptimal access control scheme. In this paper, we use simulation result to demonstrate its performance. For the alternating-based power optimization approach, in general, it may fail to locate the stationary points, not to mention the global convergence to the optimal solution. However, we have the following conclusion for Algorithm 5.Proposition 6.
Algorithm5 always has the limit point and any limit point generated by Algorithm 5 is at least a suboptimal solution for OP2.Proof.
Please see AppendixE.Complexity. The complexity of Algorithm 5 is mainly affected by the alternative optimization and the access control. The former involves searching two optimal Lagrange multipliers which have the approximate complexity of O(Nlog⁡2(λ⋆/ɛ)) and O(Nlog⁡2(ξ⋆/ɛ)), respectively, where λ⋆ and ξ⋆ are provided in Theorems 3 and 4, respectively, and ɛ is the tolerable error presented in Algorithm 5. The approximate complexity of the proposed access control scheme is O(N), where N is the number of MS in the system.


## 3.2.1. Optimal Relay Power Allocation under Fixed Jamming

We first discuss relay power allocation under fixed jamming powerp⋆B and predetermined MS set N∘. Then, we have the following optimization problem OP4:(10a)OP4:max⁡pR∑i∈N∘‍Ri∘(pi⋆B,piR)(10b)s.t.∑i∈N∘‍piR⩽PR,fixedp⋆B,N∘.Since both N∘ and p⋆B are fixed and also constraint (10b) is convex, it follows from [29] that OP4 is a degraded AF-based relay network power allocation problem with parallel orthogonal Gaussian relay channels. Therefore, the optimal solution of OP4 can be characterized by the following theorem.Theorem 3.
For OP4, its optimal solution is as follows:(11)pi⋆R=Bi2-4AiCi-Bi2Ai+,∀i∈N∘,where Ai=(γiBR)2(1+piSγiSR), Bi=γiRB(2+piSγiSR)(1+piSγiSR+pi⋆BγiBR), Ci=(1+piSγiSR+pi⋆BγiBR)2-piSγiSRγiRB(1+piSγiSR+pi⋆BγiBR)/λ⋆, and λ⋆ takes the value such that ∑i∈N∘‍pi⋆R(λ⋆,pi⋆B)=PR.Proof.
Please see AppendixC.At the moment, one may note that, in the above relay power allocation, we have not used the conditions presented in Lemma1 to restrict the relay power allocation or to exclude MS which cannot promise positive secrecy rate. The explanation is as follows: the purpose of introducing Lemma 1 is to uncover the fact that the joint relay and jammer power allocation is essentially equivalent to the joint MS access control and power allocation. In that case, the access control rather than power allocation is used to exclude MS which obtains negative secrecy rate.


## 3.2.2. Optimal Jamming under Fixed Relay Power Allocation

Now, our attention turns to the jamming power allocation at BS for givenp⋆R and N∘. Thus, we have the following optimization problem OP5:(12a)OP5:max⁡pB∑i∈N∘‍Ri∘(piB,pi⋆R)(12b)s.t.∑i∈N∘‍piB⩽PB,fixedp⋆R,N∘.It is obvious that the jamming power constraint (12b) is convex. Therefore, the focus is on the objective function (12a). In order to verify the convexity of the objective function, we take the derivative of ∑i∈N∘‍Ri∘(piB,pi⋆R) with respect to piB, ∀i∈N∘, and then have(13)∂∑i∈N∘‍Ri∘(piB,pi⋆R)∂piB=piSγiSRγiBRDipiB2+EipiB+FiOi,where Di=(1-pi⋆RγiRB)(γiBR)2, Ei=2γiBR(1+pi⋆RγiBR)>0, Fi=(1+piSγiSR)[(1+pi⋆RγiBR)(piSγiSR+pi⋆RγiRB)+1]>0, and Oi=[(1+piSγiSR+pi⋆RγiRB+piBγiBR)2+piSpi⋆RγiSRγiRB(1+piSγiSR+pi⋆RγiRB+piBγiBR)]×[(1+piBγiBR)2+piSγiSR(1+piBγiBR)]>0. It can be noted that if pi⋆RγiRB>1 (from Lemma 1, one may notice that it is the necessary condition for MS-i to obtain positive secrecy rate), then Ri∘(piB,pi⋆R) is a quasi-concave function of piB and takes the maximum at(14)pi,1⋆B=-Ei-Ei2-4DiFi2Di.Namely, Ri∘(piB,pi⋆R) is an increasing function of piB for piB∈0,pi,1⋆B and a decreasing function of piB for piB∈pi,1⋆B,∞, while if pi⋆RγiRB⩽1 (from Theorem 3, one can find that this is possible after the relay power allocation), then we have ∂Ri∘(piB,pi⋆R)/∂piB>0, ∂2Ri∘(piB,pi⋆R)/∂(piB)2<0, and ∂2Ri∘(piB,pi⋆R)/∂(piB)∂(pjB)=0, i≠j. Therefore, the objective function is a monotonically increasing and concave function with respect to piB. With the BS power constraint PB, it is known that Ri∘(piB,pi⋆R) will be maximized at piB=PB. Following these analyses, three cases may appear for OP5: (i) for ∀i∈N∘ we have pi⋆RγiRB>1 and then the objective function is quasi-concave and OP5 becomes nonconvex (but quasi-concave); (ii) for ∀i∈N∘ we have pi⋆RγiRB⩽1 and then the objective function is concave and OP5 is a convex programming; (iii) ∃i∈N∘ let pi⋆RγiRB>1, and also ∃j∈N∘ but  j≠i let pj⋆RγiRB⩽1; then OP5 is a hybrid quasi-concave and concave problem. Nevertheless, it is still possible to analyze and derive the optimal solution of OP5 by the KKT conditions [30]. In other words, the KKT conditions are the necessary conditions (not sufficient conditions due to may be nonconcavity of Ri∘(piB,pi⋆R) with respect to piB, ∀i∈N∘). Among the KKT conditions, we have(15)∂Ri∘piB,pi⋆R∂piB-ξ=0,∀i∈N∘,where ξ is the Lagrange multiplier associated with the power constraint (12b) and ξ⩾0. Based on (15), we use the following theorem to characterize the properties of the solution for OP5.Theorem 4.
For OP5 with givenp⋆R and N∘ and ∀i∈N∘, if pi⋆RγiRB>1, let pi,1⋆B be defined as (14); if pi⋆RγiRB⩽1, let pi,1⋆B=PB. Also let pi,2⋆B(ξ) be the real positive root of (15) for given ξ and define Xi(ξ)={pi,2⋆Bξ∣pi,2⋆B(ξ)>0,∂Ri∘(piB,pi⋆R)/∂piBpi,2⋆B(ξ)-ξ=0}. The unique solution of OP5 is denoted by ξ⋆,{pi⋆B:∀i∈N∘}, then we have the following: (1)
If∑i∈N∘‍pi,1⋆B⩽PB, then pi⋆B=pi,1⋆B and ξ⋆=0; otherwise, pi⋆B=p~i⋆BI{Xi(ξ⋆)≠Φ}+0×I{Xi(ξ⋆)=Φ}, where p~i⋆B=argmax⁡pi,2⋆B(ξ⋆)∈Xi(ξ⋆)pi,2⋆B(ξ⋆) and Φ is null set; ξ⋆ takes the value such that ∑i∈N∘‍pi⋆B(pi⋆R,ξ⋆)=PB.(2)
If∑i∈N∘‍pi,1⋆B>PB, then the term ∑i∈N∘‍pi⋆B(pi⋆R,ξ) is a decreasing function of ξ.Proof.
Please see AppendixD.By the results in Theorem4, for jamming power allocation, we first calculate the jamming power pi,1⋆B, ∀i∈N∘, for each accessed MS. If ∑i∈N∘‍pi,1⋆B(pi⋆R)⩽PB, that is, the sum of jamming power demand is less than or equal to the available power at BS, then all MS in N∘ will be allocated the jamming power pi⋆B=pi,1⋆B, ∀i∈N∘, while if ∑i∈N∘‍pi,1⋆B(pi⋆R)>PB, we should perform a bisection search to obtain the optimal Lagrange multiplier ξ⋆ and calculate the corresponding jamming power pi⋆B, ∀i∈N∘, by (15).


## 3.2.3. Access Control Scheme and Alternating Optimization Algorithm

Based on the alternative optimization result, a suboptimal access control scheme which can find a feasible suboptimal solution for OP3 (or OP2) is provided herein. The main idea of this scheme is as follows: when the alternative relay and jamming power optimization is convergent, we remove the MS inN∘ which has obtained the minimal (general) secrecy rate and then repeat the procedure until the system sum (general) secrecy rate does not increase. Therefore, the proposed suboptimal MS access control scheme and the AO-based power allocation together solve OP3. Hence, we have the following Algorithm 5.Algorithm 5 (alternative relay power optimization and jamming power optimization (OptROptJam)).
Step 1. Initialize t=0 and let N∘(t)=N, pB(t)=p0B.
Step 2. Based on Theorems 3 and 4, alternatively optimize the relay and jamming power until convergence and obtain piB(t) and piR(t), respectively; then go to Step 3.
Step 3. Calculate the system sum (general) secrecy rate ∑i∈N∘(t)‍Ri∘(piB(t),piR(t)); if ∑i∈N∘(t)‍Ri∘(piB(t),piR(t))>∑i∈N∘(t-1)‍Ri∘(piB(t-1),piR(t-1)), then let N∘(t)=N∘(t)∖i⋆ and t=t+1; go to Step 2, where i⋆=argmin⁡i∈N∘(t)Ri∘(t); otherwise, the algorithm is finished.In Algorithm5, t represents the iteration time of the algorithm, and N∘(t) denotes the accessed MS subset at t. In this paper, we set N∘(0)=N; that is, all MS are accessed by the system at the beginning. pB(t)=p0B is the start point of the algorithm; that is, we set a jamming power vector and then perform relay power iteration. The rule of start point setting is that it should satisfy the corresponding power constraint. The convergence condition in Step 2 is that the difference of two rounds of alternative relay and jamming power optimization is less than a threshold ɛ. In Step 3, if there are more than one MS feedback from argmin⁡i∈N∘(t)Ri∘(t), then one of them is randomly removed.Suboptimality and Convergence. For Algorithm 5, its optimality is affected by two factors: (i) the optimality of the access control scheme and (ii) the performance of alternative power optimization. For the former, the rules of constructing the optimal MS access set are unclear. Thus, it is difficult to judge the optimality of the proposed suboptimal access control scheme. In this paper, we use simulation result to demonstrate its performance. For the alternating-based power optimization approach, in general, it may fail to locate the stationary points, not to mention the global convergence to the optimal solution. However, we have the following conclusion for Algorithm 5.Proposition 6.
Algorithm5 always has the limit point and any limit point generated by Algorithm 5 is at least a suboptimal solution for OP2.Proof.
Please see AppendixE.Complexity. The complexity of Algorithm 5 is mainly affected by the alternative optimization and the access control. The former involves searching two optimal Lagrange multipliers which have the approximate complexity of O(Nlog⁡2(λ⋆/ɛ)) and O(Nlog⁡2(ξ⋆/ɛ)), respectively, where λ⋆ and ξ⋆ are provided in Theorems 3 and 4, respectively, and ɛ is the tolerable error presented in Algorithm 5. The approximate complexity of the proposed access control scheme is O(N), where N is the number of MS in the system.


## 3.3. Another Suboptimal Algorithm

For two reasons, we propose another suboptimal algorithm for our problem: (i) pursuing lower complexity algorithm and (ii) as the benchmark scheme of Algorithm5. This suboptimal algorithm is based on two suboptimal power allocation strategies at UR and BS, respectively, that is, equal relay power allocation and optimal proportional jamming power allocation. More precisely, at iteration t, the relay and jamming power allocation rules are as follows, respectively:(16)piR(t)=PRN∘t(17)piBt=PBpi,1⋆BpiRt∑i∈N∘t‍pi,1⋆BpiRtI∑i∈N∘t‍pi,1⋆BpiRt>PB+pi,1⋆BpiRtI∑i∈N∘t‍pi,1⋆BpiRt⩽PB,where pi,1⋆B(piR(t)) is defined in Theorem 4. From (16) and (17), one can find that UR equally allocates its power to the accessed MS at t and the BS power allocation is based on proportional rule. In particular, if ∑i∈N∘(t)‍pi,1⋆B(piR(t))⩽PB, then the obtained jamming power for MS-i (i∈N∘(t)) is pi,1⋆B; otherwise, its jamming power is in proportion to pi,1⋆B(piR(t))/∑i∈N∘(t)‍pi,1⋆B(piR(t)). The access control scheme for this new suboptimal algorithm is the same as Algorithm 5. To sum up, we have another suboptimal algorithm (Algorithm 7) which is summarized as below.Algorithm 7 (equal relay power allocation and optimal proportional jamming power allocation (EROptPJam)).
Step 1. Initialize t=0 and N∘(t)=N.
Step 2. Perform relay and jamming power allocation based on (16) and (17), respectively; then we have piB(t) and piR(t) and go to Step 3.
Step 3. Calculate the system sum (general) secrecy rate ∑i∈N∘(t)‍Ri∘(piB(t),piR(t)) and if ∑i∈N∘(t)‍Ri∘(piB(t),piR(t))>∑i∈N∘(t-1)‍Ri∘(piB(t-1),piR(t-1)), then let N∘(t)=N∘(t)∖i⋆ and t=t+1; go to Step 2, where i⋆=argmin⁡i∈N∘(t)Ri∘(t); otherwise, the algorithm is finished.One can note that the complexity of Algorithm7 is significantly lower. This is because we do not need to perform any power iteration but just calculate the jamming power pi,1⋆B by (14) for given relay power allocation. Therefore, it is of no doubt that Algorithm 7 has lower complexity but experiences worse performance than Algorithm 5.


## 4. Simulations Results

In this section, we present simulation results to demonstrate the performance of all proposed algorithms. Our simulation scenario is shown in Figure2, where the BS is located at (100 m, 0 m), and multiple MS are randomly distributed over the dashed rectangle region where the X coordinate is limited from −100 m to 0 m, and the Y coordinate ranges from 50 m to −50 m. The Y coordinate of the UR is fixed at 0 m, while its X coordinate varies from 0 m to 100 m. The simulation scenario and parameters’ setting are similar to those in [31, 32]. The channel gains are η(0.0097/dτ)1/2, where d is the distance between any two nodes, for example, BS and UR, or MS and UR, and the path-loss exponent is τ=4; η represents fading and is taken here as η~CN(0,1). In the simulation, the noise variance is σ2=10-8 W and ɛ=10-5 is the convergence threshold of the proposed alternatively optimization algorithms. In order to better understand the performance of our proposed algorithms, two other benchmark schemes are introduced as follows.Figure 2
The simulated network configuration.JoPowerOpt (Joint Power Constraint-Based Alternating Optimization) [33]. This algorithm is based on joint BS and UR power constraint. In other words, the system has power constraint PR+PB=P. The procedure of this algorithm is similar to OptROptJam. It is obvious that, with joint power constraint, system can obtain better performance for its ability to balance power usage between jamming and relay [34].ExJoPowerOpt (Extension of Joint Power Constraint-Based Alternating Optimization) [33]. This algorithm extends JoPowerOpt to the scenario with individual power constraint. The main idea is that when JoPowerOpt is convergent, we separately check the conditions of BS power constraint PB and UR power constraint PR; if both of them are satisfied, then ExJoPowerOpt is over; otherwise, let P=P-ΔP and repeat JoPowerOpt, where ΔP is a small positive constant.


### 4.1. General Secrecy Rate

In this section, we plot the MS (e.g., MS-i) general secrecy rate Ri∘(piB,piR) as a function of relay power piR and jamming power piB, and the result is shown in Figure 3. BS is located at (100 m, 0 m), and UR is located at (10 m, 0 m), and MS-i is located at (−90 m, 5 m). The transmission power of MS-i is piS=0.1W. From Figure 3, we can observe that, for given jamming power piB, Ri∘(piB,piR) is an increasing and concave function of piR, while, for given piR, the convexity of Ri∘(piB,piR) with respect to piB depends on the value of piR. If piR is too small, then Ri∘(piB,piR) is a concave function of piB. When piR becomes larger, then Ri∘(piB,piR) is quasi-concave of piB. This is consistent with our theoretic analysis. In addition, from Figure 3, we can find that if the relay power is too small, that is, the necessary condition (6) in Lemma 1 is not satisfied, then MS-i obtains negative general secrecy rate. (In fact, the negative general secrecy rate has no physical meaning and it is introduced for convenient illustration. Herein, the negative general secrecy rate denotes the zero value of secrecy rate.)Figure 3
General secrecy rate as a function of relay power and jamming power.


### 4.2. Performance under Different System Parameters

In this section, we evaluate the performance of the proposed algorithms under different UR positions, number of MS, and the available BS and UR power.We first evaluate the system performance under different UR positions and the results are presented in Figures4–6. In the simulation, the UR is initially placed at (80 m, 0 m) and moves to (5 m, 0 m) along a line. Therefore, UR moves away from nearby the BS to far from the BS. The power constraints are PB=2.0 W at BS and PR=1.5 W at UR. The number of MS in the system is N=15, and the transmit power of each MS is pS=piS=0.2 W, ∀i∈N. Figure 4 illustrates the achievable system sum secrecy rate versus the UR-BS distance. It is clear that the attained secrecy rates are decreasing functions of UR-BS distance. This result can be attributed to three factors: (i) the decrease of multiple MS access gain (this can be explained by Lemma 1 and verified in Figure 5; we know that when UR moves away from BS, the number of successively accessed MS decrease); (ii) the decrease of BS jamming efficiency (if UR-BS distance increases, then the channel conditions between UR and BS become worse); (iii) the increase eavesdropping ability of UR (if UR moves away from BS, then the channel conditions between UR and MS become better).Figure 4
System sum secrecy rate versus UR-BS distance.Figure 5
Number of successively accessed MS versus UR-BS distance.Figure 6
Total power consumption at BS versus UR-BS distance.From Figure4, one can note that the algorithm JoPowerOpt outperforms other three algorithms in the performance of system sum secrecy rate. However, we must highlight that JoPowerOpt is based on the scenario with joint power constraint. Thus, it is impractical for realistic system, and it is introduced just for performance comparison. The algorithm OptROptJam takes the second place, while the performance gap between JoPowerOpt and OptROptJam is small when UR is far away from the BS and it becomes larger if UR is closer to BS. Finally, this performance gap decreases if UR is nearby the BS. The increase of this performance gap is due to the fact that joint power constraint-based algorithm has the ability to balance power usage between relay and jamming. Thus, it can result in higher efficiency in resource usage, while the decrease of performance gap owes to the channel degradation between MS and UR, especially if UR is very close to BS. Moreover, in Figure 4, one can observe that the performance gap between OptROptJam and EROptPJam/ExJoPowerOpt is large when UR is far away from BS and then is reduced when UR moves to BS. In fact, the reasons for performance gap reduction for EROptPJam and ExJoPowerOpt are similar: when UR moves towards BS, all channels between MS and UR are going to be bad and the channel gains between UR and BS are going to be the same. Therefore, with the movement of UR, the optimal proportional jamming is then equal to the optimal jamming, and the performance of optimal relay power allocation is approximated by equal relay power allocation. From this figure, one can note that the performance gain by optimal relay power allocation (or optimal jamming power allocation) is only little when UR is close to BS.Figure6 shows how BS power consumption is varying with UR-BS distance. One can observe that, for all proposed algorithms, the total jamming power consumption at BS increases if UR moves away from BS. This increase is consistent with intuitive and theoretic analysis; that is, the BS jamming efficiency is high when UR is close to BS (better channel conditions between UR and BS) and it is low if UR is far away from the BS (worse channel conditions between UR and BS). In Figure 6, we see that if UR-BS distance is larger than 50 m, the algorithms OptROptJam and EROptPJam render the BS to use its whole available power to jam, while if UR-BS distance keeps decreasing, then the total jamming power consumption at BS is sharply reduced. By Figure 6, it is not difficult to uncover the reasons why algorithm ExJoPowerOpt attains worse system sum secrecy rate performance. If the UR is far away from the BS, the jamming power consumption of ExJoPowerOpt is about 1.5 W (right now, the jamming power consumption of ExJoPowerOpt is the same as JoPowerOpt). Because we have P=PB+PR=3.5 W, about 2 W (>PR=1.5 W) power is allocated to UR to perform relay. We know that, except the joint power constraint, algorithm ExJoPowerOp should also satisfy the individual power constraint at BS and UR as well. This indicates that we need to perform P=P-Δ until the power allocated to UR is not larger than 1.5 W. This procedure finally renders the total used jamming power at BS less than 1.2 W. This power is smaller than that used by JoPowerOpt. Hence, it experiences some performance loss.We then analyze the system performance under different MS numbers, and the results are depicted in Figures7 and 8. The simulation is done with the parameters PB=2.0W, PR=1.5 W, and UR is located at (10 m, 0 m). From Figure 7, we see that the system sum secrecy rates are increasing if the number of MS becomes larger. The improvement of secrecy rate can be attributed to the multiple MS access gain (i.e., see Figure 8). In other words, the more MS in the system leads to the higher efficiency. If the number of MS is few, we note that the performance gaps among these four algorithms are small. If the number of MS increases, the performance gap between JoPowerOpt and OptROptJam is still not too much but JoPowerOpt always outperforms OptROptJam, while the performance gap between OptROptJam/JoPowerOpt and EROptPJam/ExJoPowerOpt becomes larger as the number of MS increases. For EROptPJam, the increase of performance gap is due to its equal relay power allocation rule. For ExJoPowerOpt, individual power constraints limit its performance. In addition, depending on the number of MS, EROptPJam obtains better secrecy rate performance when the MS number is small, and ExJoPowerOpt dominates EROptPJam if the number of MS becomes larger. From this simulation result, one can note that if the MS number is large, we should optimize both relay and jammer power to obtain better system performance; if the MS number is only a few, the suboptimal algorithm EROptPJam is enough. From Figure 8, although the number of successively accessed MS increases when the number of MS becomes larger, one can find that the number of MS be rejected by the system increases as well, while ExJoPowerOpt always accesses the least MS among these algorithms.Figure 7
System sum secrecy rate versus number of MS.Figure 8
Number of successively accessed MS number versus number of MS.At last, how system sum secrecy rate is affected by available UR and BS power is analyzed, and the result is shown in Figure9. The number of MS in the system is N=15, and UR is located at (10 m, 0 m). Available power at BS is varying from 0.8 W to 2.0 W and total power at UR is varying in [0.5W,1.5W]. From Figure 9, we can find that the secrecy rates of OptROptJam, EROptPJam, and JoPowerOpt are all increasing functions of available UR and BS power. In fact, the conclusion that the system sum secrecy rate is an increasing function of relay power is consistent with theoretic analysis (in the proof of Theorem 3). Also, since the UR is located at 90 m far away from the BS, from Figure 6, we know that the increase of jamming power still can improve the system secrecy rate performance. Although the increasing of UR power will significantly improve the obtained secrecy rate for ExJoPowerOpt, if the UR power is small, the secrecy rate for ExJoPowerOpt seems to remain unchanged with the increase of BS power. In other words, right now the system secrecy rate is restricted by the UR power (for its individual power constraint requirements). Additionally, in Figure 9, one can see that JoPowerOpt always attained the highest secrecy rate, and OptROptJam takes the second place. This performance difference has been explained in Figure 4. For EROptPJam and ExJoPowerOpt and the given power matrix plane (formulated by the UR and BS power), if we define a diagonal line (the blue dashed line in Figure 9) from point (0.5W,0.8W) to point (1.5W,2.0W), then ExJoPowerOpt obtained better performance on the right side and EROptPJam achieves higher secrecy rate on the other side. This result can be explained as that, due to the additional individual power constraint, the secrecy rate performance of ExJoPowerOpt is more sensitive to the unbalance of available UR and BS power.Figure 9
System sum secrecy rate versus UR and BS total power.


## 4.1. General Secrecy Rate

In this section, we plot the MS (e.g., MS-i) general secrecy rate Ri∘(piB,piR) as a function of relay power piR and jamming power piB, and the result is shown in Figure 3. BS is located at (100 m, 0 m), and UR is located at (10 m, 0 m), and MS-i is located at (−90 m, 5 m). The transmission power of MS-i is piS=0.1W. From Figure 3, we can observe that, for given jamming power piB, Ri∘(piB,piR) is an increasing and concave function of piR, while, for given piR, the convexity of Ri∘(piB,piR) with respect to piB depends on the value of piR. If piR is too small, then Ri∘(piB,piR) is a concave function of piB. When piR becomes larger, then Ri∘(piB,piR) is quasi-concave of piB. This is consistent with our theoretic analysis. In addition, from Figure 3, we can find that if the relay power is too small, that is, the necessary condition (6) in Lemma 1 is not satisfied, then MS-i obtains negative general secrecy rate. (In fact, the negative general secrecy rate has no physical meaning and it is introduced for convenient illustration. Herein, the negative general secrecy rate denotes the zero value of secrecy rate.)Figure 3
General secrecy rate as a function of relay power and jamming power.


## 4.2. Performance under Different System Parameters

In this section, we evaluate the performance of the proposed algorithms under different UR positions, number of MS, and the available BS and UR power.We first evaluate the system performance under different UR positions and the results are presented in Figures4–6. In the simulation, the UR is initially placed at (80 m, 0 m) and moves to (5 m, 0 m) along a line. Therefore, UR moves away from nearby the BS to far from the BS. The power constraints are PB=2.0 W at BS and PR=1.5 W at UR. The number of MS in the system is N=15, and the transmit power of each MS is pS=piS=0.2 W, ∀i∈N. Figure 4 illustrates the achievable system sum secrecy rate versus the UR-BS distance. It is clear that the attained secrecy rates are decreasing functions of UR-BS distance. This result can be attributed to three factors: (i) the decrease of multiple MS access gain (this can be explained by Lemma 1 and verified in Figure 5; we know that when UR moves away from BS, the number of successively accessed MS decrease); (ii) the decrease of BS jamming efficiency (if UR-BS distance increases, then the channel conditions between UR and BS become worse); (iii) the increase eavesdropping ability of UR (if UR moves away from BS, then the channel conditions between UR and MS become better).Figure 4
System sum secrecy rate versus UR-BS distance.Figure 5
Number of successively accessed MS versus UR-BS distance.Figure 6
Total power consumption at BS versus UR-BS distance.From Figure4, one can note that the algorithm JoPowerOpt outperforms other three algorithms in the performance of system sum secrecy rate. However, we must highlight that JoPowerOpt is based on the scenario with joint power constraint. Thus, it is impractical for realistic system, and it is introduced just for performance comparison. The algorithm OptROptJam takes the second place, while the performance gap between JoPowerOpt and OptROptJam is small when UR is far away from the BS and it becomes larger if UR is closer to BS. Finally, this performance gap decreases if UR is nearby the BS. The increase of this performance gap is due to the fact that joint power constraint-based algorithm has the ability to balance power usage between relay and jamming. Thus, it can result in higher efficiency in resource usage, while the decrease of performance gap owes to the channel degradation between MS and UR, especially if UR is very close to BS. Moreover, in Figure 4, one can observe that the performance gap between OptROptJam and EROptPJam/ExJoPowerOpt is large when UR is far away from BS and then is reduced when UR moves to BS. In fact, the reasons for performance gap reduction for EROptPJam and ExJoPowerOpt are similar: when UR moves towards BS, all channels between MS and UR are going to be bad and the channel gains between UR and BS are going to be the same. Therefore, with the movement of UR, the optimal proportional jamming is then equal to the optimal jamming, and the performance of optimal relay power allocation is approximated by equal relay power allocation. From this figure, one can note that the performance gain by optimal relay power allocation (or optimal jamming power allocation) is only little when UR is close to BS.Figure6 shows how BS power consumption is varying with UR-BS distance. One can observe that, for all proposed algorithms, the total jamming power consumption at BS increases if UR moves away from BS. This increase is consistent with intuitive and theoretic analysis; that is, the BS jamming efficiency is high when UR is close to BS (better channel conditions between UR and BS) and it is low if UR is far away from the BS (worse channel conditions between UR and BS). In Figure 6, we see that if UR-BS distance is larger than 50 m, the algorithms OptROptJam and EROptPJam render the BS to use its whole available power to jam, while if UR-BS distance keeps decreasing, then the total jamming power consumption at BS is sharply reduced. By Figure 6, it is not difficult to uncover the reasons why algorithm ExJoPowerOpt attains worse system sum secrecy rate performance. If the UR is far away from the BS, the jamming power consumption of ExJoPowerOpt is about 1.5 W (right now, the jamming power consumption of ExJoPowerOpt is the same as JoPowerOpt). Because we have P=PB+PR=3.5 W, about 2 W (>PR=1.5 W) power is allocated to UR to perform relay. We know that, except the joint power constraint, algorithm ExJoPowerOp should also satisfy the individual power constraint at BS and UR as well. This indicates that we need to perform P=P-Δ until the power allocated to UR is not larger than 1.5 W. This procedure finally renders the total used jamming power at BS less than 1.2 W. This power is smaller than that used by JoPowerOpt. Hence, it experiences some performance loss.We then analyze the system performance under different MS numbers, and the results are depicted in Figures7 and 8. The simulation is done with the parameters PB=2.0W, PR=1.5 W, and UR is located at (10 m, 0 m). From Figure 7, we see that the system sum secrecy rates are increasing if the number of MS becomes larger. The improvement of secrecy rate can be attributed to the multiple MS access gain (i.e., see Figure 8). In other words, the more MS in the system leads to the higher efficiency. If the number of MS is few, we note that the performance gaps among these four algorithms are small. If the number of MS increases, the performance gap between JoPowerOpt and OptROptJam is still not too much but JoPowerOpt always outperforms OptROptJam, while the performance gap between OptROptJam/JoPowerOpt and EROptPJam/ExJoPowerOpt becomes larger as the number of MS increases. For EROptPJam, the increase of performance gap is due to its equal relay power allocation rule. For ExJoPowerOpt, individual power constraints limit its performance. In addition, depending on the number of MS, EROptPJam obtains better secrecy rate performance when the MS number is small, and ExJoPowerOpt dominates EROptPJam if the number of MS becomes larger. From this simulation result, one can note that if the MS number is large, we should optimize both relay and jammer power to obtain better system performance; if the MS number is only a few, the suboptimal algorithm EROptPJam is enough. From Figure 8, although the number of successively accessed MS increases when the number of MS becomes larger, one can find that the number of MS be rejected by the system increases as well, while ExJoPowerOpt always accesses the least MS among these algorithms.Figure 7
System sum secrecy rate versus number of MS.Figure 8
Number of successively accessed MS number versus number of MS.At last, how system sum secrecy rate is affected by available UR and BS power is analyzed, and the result is shown in Figure9. The number of MS in the system is N=15, and UR is located at (10 m, 0 m). Available power at BS is varying from 0.8 W to 2.0 W and total power at UR is varying in [0.5W,1.5W]. From Figure 9, we can find that the secrecy rates of OptROptJam, EROptPJam, and JoPowerOpt are all increasing functions of available UR and BS power. In fact, the conclusion that the system sum secrecy rate is an increasing function of relay power is consistent with theoretic analysis (in the proof of Theorem 3). Also, since the UR is located at 90 m far away from the BS, from Figure 6, we know that the increase of jamming power still can improve the system secrecy rate performance. Although the increasing of UR power will significantly improve the obtained secrecy rate for ExJoPowerOpt, if the UR power is small, the secrecy rate for ExJoPowerOpt seems to remain unchanged with the increase of BS power. In other words, right now the system secrecy rate is restricted by the UR power (for its individual power constraint requirements). Additionally, in Figure 9, one can see that JoPowerOpt always attained the highest secrecy rate, and OptROptJam takes the second place. This performance difference has been explained in Figure 4. For EROptPJam and ExJoPowerOpt and the given power matrix plane (formulated by the UR and BS power), if we define a diagonal line (the blue dashed line in Figure 9) from point (0.5W,0.8W) to point (1.5W,2.0W), then ExJoPowerOpt obtained better performance on the right side and EROptPJam achieves higher secrecy rate on the other side. This result can be explained as that, due to the additional individual power constraint, the secrecy rate performance of ExJoPowerOpt is more sensitive to the unbalance of available UR and BS power.Figure 9
System sum secrecy rate versus UR and BS total power.


## 5. Conclusion

In this work, based on the DBJ secure communication protocol [18], joint BS and UR power allocation to maximize system sum secrecy rate is discussed for the untrusted relay cooperation OFDMA network. We develop the joint power optimization problem and our analysis indicates that the joint power allocation is equivalent to joint MS access control and power allocation. By problem relaxation and using the dual alternating optimization approach, two suboptimal MS access control and power allocation algorithms are proposed, that is, the algorithms OptROptJam and EROptPJam. The former is iteration-based algorithm and we have proved its convergence, while the complexity of this algorithm is O(N3log⁡2(λ⋆/ɛ)log⁡2(ξ⋆/ɛ)). The later can be implemented without any iteration and it has the complexity of O(N2). Two other benchmark schemes are introduced for performance comparison. These two benchmark algorithms both have the approximate complexity of O(N2log⁡2(λ⋆/ɛ)). Our simulation results indicate that (i) the benchmark scheme JoPowerOpt always has the best secrecy rate performance. However, it is based on joint BS and UR power constraint. Thus, it is useless for realistic system; (ii) our proposed algorithm OptROptJam takes the second place in secrecy rate performance. At some UR and BS power region, if the UR is far away from BS, the performance gap between JoPowerOpt and OptROptJam is only little; (iii) the secrecy rate performance of both ExJoPowerOpt and EROptPJam depends on particular system parameters. Additionally, if UR is close to BS or the number of MS is small, the performance gap between JoPowerOpt/OptROptJam and ExJoPowerOpt/EROptPJam is negligible.



---

*Source: 530632-2015-04-21.xml*