Optimal Joint Allocation of MultiSlot Spectrum Sensing and Transfer Power in MultiChannel Cognitive Radio

In multichannel, cognitive radio (CR), the secondary user (SU) is allowed to utilize multiple subaltern frequency bands of the primary user (PU), when these bands, namely, subchannels are not currently being used. To support this spectrum reuse functionality, the SU is required to sense each subchannel, and only the subchannels wherein the PU is inactive are available for the spectrum access of the SU. In this paper, a multislot spectrum sensing and transfer scheme for multichannel CR is proposed, whose sensing stage is divided into several time slots allocated to the subchannels for spectrum sensing. While guaranteeing the spectrum sensing performance on each subchannel and limiting the interference to the PU, we formulate an optimization problem thatmaximizes the SU’s aggregate throughput by jointly allocating the optimal number of sensing time slots and the optimal transfer power to each subchannel.Theoretical analysis is given to prove the feasibility of the proposed optimization problem and simulation results are presented to show the notable improvement on the SU’s throughput when the sensing time slots and the transfer power are both optimized by the proposed scheme.


Introduction
Cognitive radio (CR) has been proposed to overcome the shortage of spectrum resources that is introduced by the traditional fixed spectrum allocation method [1].In CR, the secondary user (SU) can transfer data through an idle channel that is licensed to the primary user (PU) but currently unoccupied, providing that the PU's normal communication is not disturbed [2,3].Hence, the spectrum sensing to detect the PU's activity is an essential method to reuse the unoccupied spectrum and has become one of the key technologies in CR [4].
Spectrum sensing requests the SU to efficiently and accurately detect the presence of the PU, in order to avoid causing harmful interference to the PU [5].Since, to the SU, the PU's signal is often with unknown location, structure, and strength, energy detection, which is widely used due to its simple implementation, serves as the optimal spectrum sensing method when the detector only receives the power of the received signal [6].The sensing time has deep impact on the performance of energy detection, and the sensing performance, reflected by the false alarm and detection probabilities, can be improved with the increase of the sensing time [7].In CR, the higher the detection probability is, the less interference the PU suffers, while the lower the false alarm probability becomes, the more opportunity the SU gets to access the spectrum [8].
In recent years, multichannel CR has been proposed to improve the throughput by allowing the SU to utilize multiple idle subchannels to transfer data synchronously; however, before communicating, the SU has to detect the PU's activity in each subchannel [9].A multiband joint detection method is introduced in [10], which jointly detects the PU's signals over multiple frequency bands rather than over one band in the equal sensing time.In [11], the authors consider an orthogonal frequency division multiplexing (OFDM) based multichannel CR system and present efficient algorithms to maximize the SU's aggregate throughput by jointly optimizing both the sensing threshold and the power allocation, taking into account the influence of the false alarm and detection probabilities.However, the sensing time is not considered by the authors, which may shorten the SU's transfer time and decrease its throughput.In [12,13], a sensing-throughput tradeoff problem is analyzed mathematically, and the energy detection is used to prove that the formulated problem indeed has one optimal sensing time that yields the highest throughput.However, in multichannel CR, how to obtain the optimal set of the sensing time to all the subchannels is very important.
In this paper, subject to the constraints of the false alarm and detection probabilities of each subchannel and the aggregate interference power and rate of the SU, a multislot spectrum sensing and transfer scheme for multichannel CR is formulated.The SU's aggregate throughput is maximized by jointly allocating the optimal number of sensing time slots and the optimal transfer power to each subchannel.The analytical and numerical results obtained in this paper clearly show that the maximal throughput can be achieved when the sensing time slots and the transfer power are both optimized.

Multichannel Cognitive Radio.
In multichannel CR, the SU may operate over multiple subchannels synchronously; therefore the SU's throughput can be improved greatly.Consider a PU system, based on multicarrier modulation, operating over a wideband spectrum that is divided into  nonoverlapping narrowband subchannels.Since the PU may not use all the subchannels to transfer in the specific time and area, the left unused subchannels are available for the SU's spectrum access, as shown in Figure 1.Before accessing some subchannel, the SU must estimate the PU's state (presence or absence) in the subchannel by performing spectrum sensing, and only when the PU's absence is detected, the SU can use this subchannel to transfer data [14].However, if the PU's presence is undetected, the SU may also use the subchannel by causing harmful interference to the PU.

SU's Spectrum
Sensing.In spectrum sensing, since the SU and the PU are two different kinds of communication systems, SU is difficult to obtain any communication information from the PU system.Hence, the energy detection is widely used in the SU's spectrum sensing, which can be directly implemented without obtaining any signal information.Energy detection firstly calculates the energy of the received signal and then compares it to a preset threshold [15].
If the PU is absent, only the noise is received by the SU's detector, while if the PU is present, the mixture of the noise and the PU's signal is received.Hence, according to the different activities of the PU, the SU's received signal   in subchannel  for  = 1, 2, . . .,  is given as follows: where the states  0 and  1 denote the PU's absence and presence, respectively,   () is the PU's signal with the power of    ,   () is the Gaussian white noise with the power of  2  , ℎ  is the subchannel gain between the PU and the SU, and   is the number of sampling nodes.By supposing that the sensing time and the sampling frequency are   and   , respectively,   is given by By averaging the energy of   sampling nodes, the energy Ω   of the received signal   is obtained as follows: where the sensing SNR   = ℎ 2     / 2  .Energy detection compares Ω   to the detection threshold .And the PU's presence is decided with Ω   ≥ , while the subchannel is determined to be idle and available with Ω   < .If the PU is really absent but the PU's presence is falsely detected, the false alarm is produced, while if the PU is really present but the PU's absence is falsely detected, the miss detection happens.According to (2) and (4), the false alarm probability    and the detection probability    are given as follows: where 2) d.According to (5), with a given    ,    can be denoted by    as follows: Sensing Frame structure of single-channel CR.

System Model
3.1.Related Work.In CR, the SU's frame structure is often divided into sensing and transfer stages.In the sensing stage, the SU detects each subchannel and judges which one can be used, while, in the transfer stage, the SU uses the detected idle subchannels to transfer data.
In the sensing stage of [12], the SU searches an idle subchannel from the  PU's subchannels one by one, and once an idle subchannel  is found, the SU will stop sensing and transfer data in subchannel  through the transfer stage, as shown in Figure 2.This scheme is named as single-channel CR.
In the sensing stage of [16], the bandwidth is divided into  subbands, and each subband is allocated to sense one subchannel.Hence, the SU can use all the sensing time and part of the frequency band to sense a subchannel, as shown in Figure 3.In the transfer stage, the SU may use all the detected idle subchannels to transfer data.This scheme is named as single-slot multichannel CR.
In the sensing stage of [10], each subchannel is allocated to the equal continuous time slot to sense spectrum, and in one time slot, only one subchannel is detected by the SU, as shown in Figure 4.In the transfer stage, all the detected idle subchannels are available for the spectrum access of the SU.This scheme is named as equal-slot multichannel CR.However, the gain of each subchannel is often different, and the sensing performance can be further improved if we allocate more sensing time to the subchannel with better sensing gain.

Proposed Scheme.
In this paper, the sensing stage is further divided into many small time slots, as shown in Figure 5.In each time slot, the SU can sense any subchannel with energy detection, and the number inside each time slot means the subchannel index to be sensed.Note that in the example, the SU senses spectrum firstly in channel 2, then in channels 1 and 3, and so forth, and finally in channels  − 1 and .Actually, the sensing order of the subchannels does not affect the final sensing performance of each subchannel.At the end of the sensing stage, the SU can obtain the used statuses of the  subchannels by combining the sensing information in each time slot.The number of the sensing time

Sensing Transfer
Figure 4: Frame structure of equal-slot multichannel CR. slots allocated to each subchannel can be different, which is chosen according to the subchannel gain.By supposing that the frame period is , the length of each time slot is , and the number of the time slots to sense subchannel  is   , the sensing time of subchannel  is given by and the SU's transfer time is given by By assuming that the subchannels are symmetrical, the SU's data transfer in subchannel  includes the following two cases.(a) If the PU's absence is detected accurately by the SU in the probability of 1−   , the SU can transfer in the subchannel effectively, with the rate of where   is the subchannel gain between the SU's transmitter and receiver.(b) If the PU's presence is undetected in the miss detection probability of 1 −    , the SU may use the same subchannel with the PU and cause severe interference to the PU, with the rate of Obviously,  0  >  1  .In case (b), if the PU is really present but the SU detects the PU's absence falsely in the probability of 1 −    , the SU will transfer data with the power of   that causes severe interference to the PU.Hence, the SU's aggregate interference power is given by should be constrained by letting   ≤  max , where  max is the maximal interference power suffered by the PU.Moreover in case (b), the SU transfers data with the rate of  1  at the presence of the PU; these data may also disturb the PU's normal communication, and therefore we can get the maximal aggregate interference rate of the SU as follows: It is also necessary to constrain the SU's interference rate as As log 2 () is a convex function, we get the Cauchy inequality as follows: By supposing  min = min{ 1 ,  2 , . . .,   }, from ( 13) and ( 14), we may derive the inequality as follows: where   =   /  .Hence, we can guarantee   ≤  max  by letting from which we further obtain that where  max = (2  max  / − 1)(1 +  min ).In each frame, the SU may transfer data in cases (a) and (b) with the rates of  0  and  1  at the absence and presence of the PU, respectively, during the transfer time   .The SU's aggregate throughput over  subchannels is given as follows: where   0 and   1 are the appearance probabilities of the states  0 and  1 , respectively.In this paper, our goal of jointly allocating sensing time slots and transfer power is to maximize the SU's aggregate throughput subject to the constraints on the false alarm probability, the detection probability and the transfer power of each subchannel, and the aggregate interference power and the aggregate interference rate of the SU, as follows: where the vectors k = { 1 ,  2 , . . .,   } and p = { 1 ,  2 , . . .,   },  is the upper bound of the false alarm probability,  is the lower bound of the detection probability,  max is the SU's maximal transfer power, and ⌈⌉ denotes the maximal integer no bigger than .In (19), we set  ≤ 0.5 and  ≥ 0.5, because the SU's sensing performance must be good enough in order to decrease the interference to the PU and improve the spectrum access of the SU, as indicated in CR. 19) is a multiparameter optimization problem about k and p, which can be solved by the alternating direction optimization (ADO) [16].In each iteration, one of the two parameters is fixed while the other one is optimized, and the two parameters are alternately optimized through some iterations until both of them tend to be constant.Firstly, we fix p (corresponding  0  and

Solution to Optimization Problem. Equation (
where  0  =  −1 ()(1 +   ) and  1  =   √  .From (18), since  improves with the increase of both    and    , it can achieve the maximum only when    = , and correspondingly, the equation of (20) is obtained.By substituting (20) into    ≤ , we get ( 0  +  1  √  ) ≤ , from which we obtain   ≥   , where   is given by where ⌊⌋ denotes the minimal integer no less than .
We suppose that the aggregate number of sensing time slots  = ∑  =1   , where  ≥ ∑  =1   and  ∈ .With the fixed p, the constraints of (19) about p can be ignored, while, with   ≥   , the constraints of false alarm and detection probabilities can be satisfied.Hence, (19) transforms the suboptimization problem about k as follows: In the following proposition, we will prove that (k) is a convex function about the vector k, which denotes that there is an optimal set of k that maximizes (k).
If  is given, since k is an integer vector, the optimal   for  = 1, 2, . . .,  can be obtained by the greedy algorithm, as shown in Figure 4.In this figure,   (  ) is denoted by As  is an integer within a specific range, it is not complicated to search through .With any given integer  ∈ [∑  =1   , ⌈/⌉], we calculate (k) with the optimal vector k obtained by the greedy algorithm, and then we choose the optimal  that maximizes (k), as follows:

Yes
No With the given , since the increment range of ∑  =1   is from ∑  =1   to , the time complexity of the greedy algorithm in Figure 6 is (( − ∑  =1   )log 2  + ).From (25), by enumerating  ∈ [∑  =1   , ⌈/⌉], the aggregate time complexity to search an optimal set of k is given by Secondly, we fix k (corresponding   ,    , and    are all fixed) and optimize p.With the fixed k, the constraints of false alarm and detection probabilities can be ignored.Noting    = , (19) transforms the suboptimization problem about p as follows: where  0  =    0    0 and  1  =    1    1 .Using the Lagrange multiplier algorithm to solve (27), we obtain the Lagrange function as follows: where  1 ,  2 , and  3 are the Lagrange multipliers.According to the Karush-Kuhn-Tucker (KKT) conditions, the optimal p is got by letting /p = 0, and if  1 , 2 , and  3 are not zero, they are got by letting the corresponding constraints acquire equations.Hence, we obtain Noting   ≥ 0, direct calculation /  = 0 yields where Using the subgradient method, we can further obtain  1 ,  2 , and  3 through the iterations as follows: where  is the iteration index and  () is a sequence of scalar step size.Once  1 ,  2 , and  3 are all obtained, we can get the optimal transfer power   for  = 1, 2, . . .,  through (30).
The time complexity of the Lagrange multiplier algorithm depends on the estimation accuracy  ∈ (0, 1).By supposing that the stopping conditions of the iteration in (31 3 | ≤ , the iteration complexity is (1/ 3 ), and the time complexity of each iteration is ( log 2 ).Hence, the aggregate time complexity of the Lagrange multiplier algorithm is given by The ADO method is used to obtain the optimal joint allocation of multislot spectrum sensing and transfer power by optimizing the two suboptimization problems ( 22) and ( 27) alternately and iteratively.The joint allocation algorithm based on ADO is described in Figure 7. Noting that the goal of optimizing ( 22) and ( 27) is to achieve the locally maximal value of the objective function, the value of the objective function is nondecreasing in each iteration as follows: where  is the iteration index.Equation (31) indicates that if (k, p) is convergent, both of k and p must be convergent.That is, if the stopping conditions of the iteration in ADO are ‖k () − k (−1) ‖ ≤  and ‖p () − p (−1) ‖ ≤ , the optimal (k () , p () ) can be obtained.
As the stopping conditions of the joint allocation algorithm are ‖p () − p (−1) ‖ ≤  and ‖k () − k (−1) ‖ ≤ , the iteration complexity of the algorithm is (1/ 2 ).In each iteration, the greedy algorithm with the complexity of  GE and the Lagrange multiplier algorithm with the complexity of  LM are both implemented once.Hence, the aggregate and ‖ ‖ ≤ ?

Fixing
, calculate = {k l } L l=1 by p (n)  k (n+1)   Initialize: the iteration index n = 0, the slot vector Choose the power vector with any p (n)   the greedy algorithm in Figure 4 Let n = n + 1

Yes
No the estimation accuracy , time complexity of the joint allocation algorithm is given as follows: However, since, in the conventional exhaustive searching method, we need to search through the vectors k and p including 2 parameters, the time complexity of the exhaustive searching method is (1/ 2 ) that greatly exceeds the time complexity of the proposed scheme if  ≥ 3.
Figure 8 shows the SU's aggregate throughput  versus the lower bound of detection probability , with different number of sensing time slots .From this figure, it is seen that there exists an optimal  (the aggregate sensing time is ) that maximizes , and the convex-shaped  is also consistent with Proposition 1, which verifies the correctness of our allocation scheme.In the traditional equalslot multichannel CR of [10], the sensing time allocated to each subchannel is /.With the same aggregate sensing time, the throughput of the proposed scheme is always larger than those of the equal-slot multichannel CR and the singleslot multichannel CR because, in our scheme, we allocate different sensing time to each subchannel according to their different gains.We also see that  decreases greatly with the increase of , because the false alarm probability increases with the increase of , which reduces the SU's spectrum access.Figure 9 shows  versus the maximal transfer power  max , with different .It is seen that the proposed scheme outperforms the other two schemes, with various  max .
Figure 10 indicates the uniform-time throughput / versus the frame period , with different .It is seen that with the increase of , the number of the available sensing time slots increases, and longer sensing time can improve the SU's sensing performance.With the improved sensing performance, the decrease of the false alarm probability may increase the SU's spectrum access, and thus the SU's throughput is improved.Figure 11 compares the achieved throughput of the proposed scheme and the single-channel CR.It is seen that the proposed scheme that uses multiple subchannels to transfer can get much more throughput than the single-channel CR does.
Figure 12 shows the transfer power and false alarm probability of each subchannel by the proposed scheme.It is seen that in our allocation scheme, with the given detection probability, the larger power is allocated to the subchannels with lower false alarm probabilities, because the SU has more opportunities to access these subchannels, and larger power can help the SU make full use of the better subchannels in order to achieve perfect throughput.However, no power is allocated to the subchannels with much higher false alarm probabilities (i.e., subchannels 1 and 9), because the SU often cannot access these subchannels, and the transfer power allocation in these worse subchannels will not achieve benefits.

Conclusions
In this paper, we propose a multislot spectrum sensing and transfer scheme of multichannel CR.The frame is divided into sensing and transfer stages, and the sensing stage is further divided into many small time slots that are allocated to each subchannel for energy detection.An optimization problem of jointly allocating sensing time slots and transfer  power to each subchannel is formulated to maximize the SU's aggregate throughput, subject to the constraints of the false alarm and detection probabilities of each subchannel and the aggregate interference power and rate of the SU.A joint allocation algorithm based on alternating directing optimization is proposed to solve the optimization problem.With numerical results, it is shown that there exists an optimal set of sensing time slots that maximizes the throughput

Figure 3 :
Figure 3: Frame structure of single-slot multichannel CR.

Figure 5 :
Figure 5: Proposed frame structure of multichannel cognitive radio.

Figure 6 :
Figure 6: Flowchart to solve   based on greedy algorithm.

Figure 7 :
Figure 7: Joint allocation algorithm based on ADO.

Figure 8 :
Figure 8: Throughput versus lower bound of detection probability.
multichannel CR Single-slot multichannel CR