In this article, adaptive resource allocation (ARA) is investigated for multiple primary networks based-cognitive radio networks under a more practical system model, where the bandwidth of each secondary user is assumed to be limited and the maximum allowable interference for each primary network is different. We first formulate the ARA as a constrained optimization problem with the objective function of maximizing the proportional fairness-based ergodic sum capacity. The multiple constraints optimization problem is NP-hard and therefore, we propose a scheme to decompose the optimization problem into two unconstrained optimization problems by designing alternative objective functions and penalty functions. Then, a suboptimal heuristic solution framework based on particle swarm optimization is proposed to solve the unconstrained optimization problems. Computation simulations are carried out and the results show that the proposed scheme outperforms traditional ARA schemes.
Keywords:Cognitive radio networks; Resource allocation; Multiple primary networks; Particle swarm optimization
Recently, cognitive radio networks (CRNs)  get extensive attentions to alleviate the contradiction between the scarcity and low-utilization of the spectrum resources. In CRNs, two approaches are widely studied for dynamic spectrum access (DSA), termed as underlay mode and overlay mode. In underlay mode, secondary users (SU) transmit over all the frequency bands as long as the serviced SUs do not cause excessive interference to primary users (PU). One possible drawback of this mode is that it requires exact information of all primary receivers’ positions, which usually expenses a high complexity of hardware and computation. On the other hand, the basic idea of overlay is to detect the absence/presence of licensed primary radios, and opportunistically use the idle band for transmissions without causing harmful interference to the authorized signals [2,3]. Compared to the underlay mode, the overlay mode is more practical and easily accepted by primary systems. Also, the overlay spectrum access can further improve spectrum efficiency by exploiting multiuser diversity and adaptive resource allocation (ARA) algorithm. In this article, we consider the problem of ARA for overlay CRNs.
ARA is one of the most effective methods to improve the spectral efficiency and has widely been investigated in the last decade. Unlike existing wireless communication networks, the ARA problem in CRNs is more challenging. First, the detection errors of spectrum sensing greatly deteriorates the performance of ARA algorithm [4,5]. The authors of  proposed an iterative algorithm to obtain the optimal sensing time and the corresponding power allocation strategy with imperfect sensing information. The authors of  presented a primal-dual decomposition-based cross-layer scheduling for power allocation and subchannel assignment with raw sensing information. Furthermore, cognitive radios are usually required to be self-regulating, to control the interference to primary systems in a tolerable level [5-7]. The authors of  considered the average interference constraint on the nearest PU receiver. In , a total interference constraint to a certain PU who uses the same channel with SU is imposed. In [8,9], the interference control is considered as a total power constraint of all SUs.
Besides, because of the absence of the license, the available spectrum resources for SUs are usually temporary and unstable, which results in transmission interruption with high probability. One of considerable solutions is to enable SUs to search for spectrum resources from multiple primary networks (PRNs) simultaneously. The authors of  considered the mobile users operating on different radio access technologies (RATs), and the single-user case is analyzed. In the multiple PRNs, however, the available spectrum resources are usually quite wide as well as discontinuous. This leads to that the cognitive radios operating on such spectrum resources are price prohibitive due to the increase of sampling rate. In addition, different PRNs restrict different power values of CRNs in one resource allocation, which greatly increases the computational complexity of ARA.
This article aims to design an effective ARA algorithm in multiple PRNs environment. Specifically, we assume that the spectrum resources available for SUs are distributed in multiple PRNs, and discontinuous in frequency domain. Furthermore, the sampling rates of SUs are limited for practical consideration. Other constraints concerned in traditional multiuser communications networks like proportional fairness and bit error rate (BER) are also involved in this study.
We formulate these concerns as a constrained optimization problem with the objective function of maximizing sum capacity. The formulation is a mixed integer nonlinear programming (MINP) that is a typical NP-hard problem. Therefore, the most intractable problem is to find a suboptimal alternative with polynomial time complexity. In this article, we proposed an integral solution to the optimization problem. The main works and contributions of this article can be summarized as follows.
(1) The ARA problem for CRN in multiple PRNs is formulated as a constrained optimization problem.
(2) A hybrid alternative objective function instead of sum capacity is proposed to balance the tradeoff between the complexity of optimization and the speed of convergence.
(3) The optimization problem is transformed into two unconstrained problems through the design of penalty functions. In particular, we propose multiplicative penalty functions that are independent with channel gains and power limits, contributing to the algorithm more robust.
(4) A particle swarm optimization (PSO)-based heuristic intelligence algorithm is proposed to solve the two unconstrained problems. Specifically, we first present an initial topology of PSO according to the fact that an SU suffers channel fading of the same level for all frequencies (because large scale attenuation plays the key role for wireless channel). The proposed topology can greatly decrease the computation complexity from O(KN) to O(K3). Based on this, we design different particle structures for subchannels and power allocation, respectively.
In order to verify the proposed ARA algorithm, three computer simulations are carried out. The first one illustrates the efficiency of the proposed hybrid objective function in terms of convergence and quality; for comparison, we put a classical ARA algorithm  into the same scenario. The simulation results reflect that the proposed outperforms the algorithm in  both in capacity and proportional fairness; at last, we discuss the effect of imperfect sensing information, and a possible optimization of sensing strategy is proposed.
The rest of this article is organized as follows. The following Section “Related studies” introduces the state of the art in ARA for CRNs. In Section “System model and problem formulation”, we describe the system model and formulate the ARA as a constrained optimization problem. In Section “Objective function and problem decomposition”, we decompose the optimization into two unconstrained optimization problems. In Section “Particle swarm optimization based framework”, we discuss the solution to the optimization problems by a PSO based algorithm in detail. The computational complexity of proposed is analyzed in Section “Computational complexity analysis”. To verify the proposed scheme, computer simulations are carried out in Section “Results and discussion”, and the brief conclusions are drawn at last.
The main foundation of ARA is based on that different users suffer from different wireless channel fading. ARA algorithms can adaptively allocate a dimension to the users, and achieve higher capacity. Two classes of resource allocation techniques have widely been addressed in traditional OFDM-based wireless communication systems, namely, (1) margin adaptive (MA)  and (2) rate adaptive (RA) . The MA aims to minimize the overall transmit power, given the constraints on the users transmit rates or BER. The RA aims to maximize users’ capacity sum, given a total power constraint. Both optimization problems are MINP problem, and computationally prohibitive. The authors of  presented a suboptimal resource allocation algorithm by transforming the MA problem into a convex optimization problem. The authors of [11,14] discussed RA problem with the constraint of proportional rate, and proposed an optimal power allocation scheme by the Newton–Raphson method.
Right now, a lot of research works is currently ongoing to deal with the optimization of the spectrum utilization in CRNs. The authors of  presented a primal-dual decomposition approach for the RA problem with raw sensing information. However, they ignored proportional rate and fairness, which ensure that all users suffering different channel fading can experience the similar quality of service [15,16]. In this article, we assure proportional fairness among SUs by imposing a set of nonlinear constraints into the optimization problem and formulate a new optimization problem. Regarding to interference control, an important concept named “interference temperature”  is widely used to evaluate the interference on the primary system. Wang et al.  considered the average interference constraint on the nearest PU receiver. Nguyen and Lee  imposed a total interference constraint on PU, who uses the same channel with SU. Although “interference temperature” is an exact measurement, it is difficult to implement in realistic systems. To calculate the aggregation interference on each PU, the sensing algorithm based on receivers is needed, which usually requires information of receivers’ positions and channel gains from each SU transmitter to the corresponding PU receiver. In this article, we simplify the constraint with an only total power limit on specific PRN. This simplification is reasonable since the value of power constraint could be negotiated between PRN and CRN through a spectrum broker . Choi et al.  considered the mobile users who can access different RATs, and the single-user case was analyzed. In this article, we equip the SU as the same capability of operating on crossing bandwidth and extend the scenario to multiuser case. Further, we consider that the sampling capability of SU is limited and only part of the spectrum resources can be exploited at any time.
Because of extra constraints in CRNs, it is nontrivial to obtain optimal solution for ARA problem. Bio-inspired swarm intelligence-based optimization algorithms [19,20], as an effective alternative solution, have widely been adopted in CRNs. A successful example of exploring the swarm intelligence algorithms in CRN is the testbed designed by Virginia Tech group [21,22], which shown the flexibility and stability of CRN with the genetic algorithm (GA)-embedded cognitive engine. Newman et al.  proposed a GA-based cognitive radio suitable for Emergency (minimize BER) and Low Power (minimize power consumption, like MA in a traditional network scenarios). Since then, lots of alternative schemes such as the quantum genetic algorithm , cross entropy , and PSO  is exploited to address the problems of resource allocation. PSO is proposed as a new swarm intelligence algorithm in recent years due to its capabilities of convergence rapidity, optima finding, and matching problem easily. In this article, we consider PSO as a tool to solve the problem of resource allocation in CRNs.
System model and problem formulation
We consider an infrastructure-based CRN system under multiple PRNs, as depicted in Figure 1. PRN i is licensed Ninonoverlapping orthogonal frequency subchannels. The bandwidth of each subchannel is W Hz. For simplicity, we assume that there is no guard band and then the spectrum spacing between two adjacent subchannels is also W Hz. The coverage radius of cognitive base station is d and the distance between each primary base station (PBS) and cognitive base station (CBS) is 2d. Within one resource allocation, there exist K active SUs in the CRN.
Figure 1. The cognitive OFDM system with multiple primary networks.
Before performing resource allocation, the CBS first senses the primary signals from its neighboring PBSs. According to the sensing results, the CBS decides whether the corresponding spectrum band is busy or idle. Then, based on the spectrum states, the CBS dynamically allocates the available spectrum resources to SUs in service. In practice, the spectrum sensing is not always perfect and accurate. Two sensing errors, miss detection (MD) and false alarm (FA) are inevitable for spectrum sensing. An MD means that an SU identifies a channel unoccupied while in fact a PU is transmitting on which. An FA means that an SU identify a channel BUSY while that is actually IDLE. We denote the FA probability and MD probability as Pf and Pm, respectively. Both of them are the signal-to-noise ratio (SNR) as well as sensing algorithm dependent. If there are busy subchannels and idle subchannels in PRN i, the number of available subchannels for CRN is
Hence, the total available subchannels for CRN is . We assume that the resource allocation period is updated as fast as CSI feedback between CBS and SUs, and perfect instantaneous CSI is assumed available for the CBS. It is also assumed that information of both resource allocation results and CSI are transmitted through a cognitive pilot channel .
where N0 is the two-sided noise spectral density. We assume that the noise is additive white gaussian noise, and all spectrum has the same value of N0. is the SINR gap due to modulation, and BER is the BER for SUs. represents the interference on subchannel m of SU k, from the PRN i. is the channel gain for SU k in subchannel m of the PRN i and is the power allocated in subchannel m of the PRN i, can either be 1 or 0, indicating whether subchannel m of PRN i is allocated to SU k or not.
Hence, the total sum capacity for all SUs can be written as:
To avoid interference on PRNs or control interference at a tolerable level, a power control scheme is necessary in CRNs. Traditionally, the “interference temperature” is an evaluation indicator representing the aggregation interference of SUs in terms of the power spectrum . In order to evaluate the interference temperature from SU transmitter to each PU receiver, we need the position information of PUs and the channel gains between the SU transmitter and all PU receivers. This requirement is too expensive for implementation. Therefore, in this article, we adopt a simple scheme which restrains the total power of CRN in each PRN [8,9]. Denote as the total available power in PRN i, and the power constraint can be written as:
In addition, proportional fairness is also an important performance metric in multiuser wireless communication systems. Under the constraint of proportional fairness, we can explicitly control the capacity ratios among users, and generally ensure that each user can achieve its target data rate, especially for cell-edge users. In this article, we adopt the definition of proportional fairness index the same as in .
where is the sum capacity of SU k in one allocation period and γ kis the expected transmit rate of SU k. The maximum value of 1 to be achieved when R1/γ1=Rk/γk,∀k. To ensure each SU is served in one allocation period, we define a noncontinuous point Φ=0, if ∀Rk=0.
Notably, to guarantee the total transmit rate, cognitive radios have to collect the spectrum resources from multiple PRNs. This leads that the available spectrum for CRN is usually band-crossing as well as discontinuous, as illustrated in Figure 2. If we assume that the cognitive radio can transmit over all spectrum at the same time, then the resource allocation is as trivial as that in a traditional network. However, this assumption seems too strong in practice. According to the theory of bandpass sampling, the sampling rate for bandpass signal of bandwidth B must satisfy the following equation :
Figure 2. Spectrum resources distribution in multiple primary networks.
where m is an arbitrary integer, fsis the sampling rate, fcand B are the central frequency and bandwidth of bandpass signal, respectively. That is to say, the maximal bandwidth a radio can operate on is subject to the A/D capability, which is price prohibitive with the increase of fs. In this article, we consider SU is limited by bandpass capability, and denote as the width of maximal bandpass capability of SU k. Without loss of generality, we index the subchannel by lexicographical order, and denote fnas the central frequency of subchannel n
where the subchannel n belongs to the PRN m and . ψm is the starting frequency of PRN m. We denote Ωkas the subchannels set allocated to SU k in one allocation and define and . The constraint of bandpass capability can be expressed as follows:
Mathematically, the optimization objective (3), and constraints (4), (5), (8) can be modeled as the following mixed-integer nonlinear programming:
where and indicate the assignments of subchannel and power, respectively. The constraint C1 implies all power value allocated to subchannels are positive; C2 ensures the interference to PRNs at a tolerable level; C3 indicates that each subchannel can only be used by one SU at any time; C4 is the bandpass constraint of SU; C5 is the proportional fairness constraint.
Objective function and problem decomposition
For a CRN system consisting of N available subchannels and K SUs that can crossover at most m subchannels, there are at most possible subchannel allocations. The feasible set is too huge to exhaust. Besides, the nonlinear constraints C4 and C5 in (9) increase the complexity to obtain the optimization solution. Therefore, it is prohibitive to find an optimizer in terms of computational complexity. In this article, we attempt to find a suboptimal alternative to decrease the complexity significantly while still delivering performance close to the global optimum.
Before designing of the suboptimal algorithm, we attempt to simplify the problem in (9) first. In Section “Objective function”, we analyze two equivalent objective functions to delete the constraint C5, and then a hybrid objective functions is proposed. In Section “Problem decomposition”, the problem in (9) is decomposed into subchannels allocation and power allocation, and each of them has only two constraints. Section “Penalty function” discusses the design of penalty functions, by which transforming the two sub-problems into unconstrained optimization problems.
The global objective function is the maximum of ergodic sum capacity of SUs. It is difficult to achieve due to the constraint of proportional fairness. We need an equivalent objective function that is maximized simultaneously with global objective function. The maximum of proportional fairness and maximum of minimal capacity SU are two alternatives. For convenience, we denote the objective function of maximum of proportional fairness as MPF and maximum of minimal capacity SU as Max–Min–SU.
We relax the constraint of proportional fairness Φ∈(ε,1], where ε is a rational number close to 1. Then we extract the constraint C5 in (9) as an objective and remodel the problem as a multi-objective problem that:
We assume that for ∀ε,Φ∈(ε,1], , and the element number of is denoted as . It is trivial to verify that for ∀ε∗≥ε, , and . Thus, . It means is the asymptotical minimal cover for problem (9). If , there is no optimal solution for problem (1); if , then (ρ,p) is the optima we search for; if , is Pareto optima for problem (9) and the global optima must be included in . Therefore, we can get an approximately equivalent optimization problem as follows:
Due to the noncontinuous of feasible set, generally can only find a set , such that Φ∈(ε,1]. It is nontrivial to find out the Pareto-optimal front of set . The question becomes if there exist a equivalent function such that if and only if . Max–Min–User is one of them.
where the constraints are the same as in (10). □
Hybrid objective function
As discuss above, we proposed two alternative objective functions MPF and Max–Min–SU. For MPF, we adopt method narrowing the feasible set to search for the approximate solution. This model has fast velocity of convergence since its aim is Pareto optima for problem (9), at the expense of quality of solution; for Max–Min–SU in a contrary, we derive that its solution is equivalent to the optimization problem (9). However, due to the uniqueness, the searching process is more different than that in MPF. Therefore, we propose a hybrid objective function that possesses advantages of both MPF and Max–Min–SU. The format expression is as follows:
where the constraints are the same as in (10). And we will discuss the performance for the three fitness functions in Section “Particle swarm optimization based framework”.
Ideally, subchannel and power should be allocated jointly to achieve the optimal solution in (12). However, due to the mixed binary integer programming, it is prohibitive computational burden at CBS, which leads the increase of computation cost and allocation delay. Hence, we separates the problem (10) into subchannel allocation and power allocation to reduce the complexity, because the continuous variable and binary variable can be handled independently.
We first formulate subchannel allocation as an optimization problem with the assumption that all subchannels are allocated equal power distribution in each PRN, i.e.,
Thus the problem of subchannel allocation can be formulated as follows:
For a determined subchannel allocation, the optimization problem is formulated as follows:
The penalty function technique is an effective method to solve the constrained optimization problems . By penalizing the constraints and building a single objective function, the constrained problem can be transformed into an unconstrained one, which is in turn maximized using an unconstrained optimization algorithm. A penalty function is generally defined as :
where f(x) is the original objective function of the optimization problem. h(k) is a dynamically modified penalty value, k is the algorithm’s current iteration number; is the feasible set and H(x) is a penalty factor, which is always problem dependent.
The definition of penalty function in (17) is additive. In this article, two penalty functions based on (17) to solve the constraint C2 in problem (15) and (16) can be formulated as follows.
Additive penalty function
where IA(x) is indicative function. IA(x)=0 if x∉A and IA(x)=1 if x∈A. λ1 and λ2are the Lagrangian multipliers. In our problem, they are dependant on channel gains and power allocation results , resulting in changing in every resource allocation. This will introduce great burden in CBS. To simplify this problem, we propose multiplicative penalty functions, which is independent with and . Formally, the expressions of our proposed are written as follows.
Multiplicative penalty function
where β1and β2 are normalized penalty factors and enable penalty value to be adjusted dynamically. The expressions can be written as
By analysis above, we transform the problem (9) into an alternative with objective functions (20) (21) and constraint C1 in (15) (16). In the next section, we will propose a heuristic algorithm to solve the alterative problem.
Particle swarm optimization based framework
The PSO algorithm was introduced by James Kennedy and Eberhart (1995) as an effective batch of heuristic algorithms . Compared with other algorithms, PSO has better global searching ability at the beginning of the run and a local searching ability near the end of the run . It is also effective for integer programming .
The standard PSO algorithm basically involves the following steps :
(1) Construct the particle to map the solution of interest problem;
(2) Create the initial topology for swarm and parameters to initial the optimization;
(3) Calculate fitness value for each particle;
(4) Renew particle position;
(5) Return to (2) until the solution satisfies the requirement of interest problem.
Following the work, in this section we design a PSO-based framework to solve the subchannel allocation and power allocation for CR system. At the beginning of Section “PSO structure”, we discuss the design of particle structure to establish the mapping between PSO and resource allocation problem. Then, in order to improve the velocity of convergence, we propose an initial topology of swarm exploiting the channel characteristics in Section “Suboptimal power distribution for a fixed sub-channel allocation”; Section “Particle renew” discusses the particle renew and the process of the proposed algorithm.
Suboptimal subchannel allocation
Due to the constraint of proportional fairness, each SU should be allocated at least one subchannel in an resource allocation duration. Thus, we can exclude all of particles who do not traverse the set of SU. Furthermore, the subchannels allocated to one SU should be approximately bunching since the sampling ability is limited. So we can arrange SU sequentially and user k is preallocated |Ωk| subchannels as shown in Figure 3.
Figure 3. The base particle for subchannel allocation.
We name this particle as “base particle” because it is the start of particle structure for subchannel allocation and denote it as a 1×Nsensed dimension vector particlebase. Later we will explain the relationship between “base particle” and “particle” who, in fact, executes the solution searching work.
To establish particlebase, a method is needed to estimate variable |Ωk| for each SU. The capacity equation for SU k on subchannel n can be rewritten as follows:
For the purposes of analysis, we define the virtual power P as:
Thus we can get the expression as:
as the virtual channel-to-noise ratio for user k in subchannel n. Then (22) can be rewritten as:
According to channel fading theory, large-scale fading is the main factor that affects the value of Hk,n. Hence, for the same SU, Hk,n for all subchannel is approximate or fluctuate in the same level. Therefore, we can approximate function log(1 + x) as linear if the variation of x is small, as shown in Figure 4.
Figure 4. Illustration of approximation for Equations (24) and (25).
On the basis of the approximation, we can estimate the capacity sum for SU k as:
To satisfy the requirement of proportional fairness, we should make
Therefore, we can get the relationships as follows
and the total sensed subchannels Nsensecan be expressed as
Substituting (23) into (24) and we can get the approximate |Ωk| for each SU as:
particlebase is not optimal because we just arrange SU order in a sequential way. Intuitively, we can use sorting algorithm to rearrange SU order to achieve higher objective (fitness) value. Bubble Sorting is a classical sorting algorithm with complexity O(n2). In this article, we use a modified bubble sorting algorithm to rearrange SU order. We denote the SU order as a 1×K dimension vector K and the element of K indicates the SU index. fitnessK represents the fitness value with the SU order K and can be calculated by (20). The process is depicted in the following:
Algorithm 1 SU sorting algorithm
Generate K in a random permutation; for k=1 to K−1 for l=1 to K−k−1 Calculate fitnessKby (20); Ktemp←Ktemp(l)←K(l + 1)Ktemp(l + 1)←K(l) Calculate fitnessKtempby (20); if fitnessKtemp>fitnessKK←Ktempelse continue; end for end for
We denote the result of sorting as . The feasible set in this initial topology is only K, which is much less than the whole feasible set . The subchannel allocation and the corresponding elements are illustrated in Figure 5.
Figure 5. The optimized base particle for subchannel allocation.
It is an intuitive idea to generate particle based on . However, since variable K is usually much smaller than Nsense, the activity of particle with variable SU index is limited extremely. Furthermore, with the renew of such particle, it is difficult to control the solution under the constraint C4. Therefore, we design a particle with variable space [0,Nsense]. First, we generate a 1×Nsense dimension particle, and denoted as particleinit:
Based on particleinit, we can generate a set of 1×Nsensedimension vector particle as follows:
for i=1 to Npso
for n=1 to Nsense
where Npso is the number of particles, particlei refer to the ith particle. ω follows uniform distribution in [0,Nsense].
It is noticed that the elements of particle are no physical significance. To calculate the fitness value for each particle, we should transform the vector particle into a vector with the subchannel allocation information. Therefore, we define a vector, namely particleshadow corresponding to particle. The mapping relationship between particleshadowand particle can be written as:
(x,y) represents the operation of modulus after division and we define special point mod(x,y)=x if x/y=m, m is an arbitrary integer. particleshadow,i(n) indicates the nth element of the ith particleshadow.
The reason why we call the particle “shadow particle” lies in that each shadow particle maps a particle who is responsible for searching optimum. To illustrate the relationships among these particles, Figure 6 shows the transformations by taking an example with 8 SUs and 8 subchannels case.
Figure 6. Illustration for particle transformations.
The Rk in the particlei can be obtained as:
Substituting (32) into (20) and we can get the fitness(i) of the ith particle.
Suboptimal power distribution for a fixed sub-channel allocation
Once suboptimal subchannel allocation is determined, we shouldrestructure the particle, since the requirement has been changed in power allocation. Each element in vector particleshadow(n) is changed to represent a power value. For example, the fourth subchannel is allocated to power 0.011 W and the transmission rate for each SU can be calculated by joint subchannel and power allocation results, as shown in Figure 7.
Figure 7. The particle structure for power allocation.
According to the objective function, the search space we focus is Nsense-dimensional, the ith particle of the swarm is represented by a Nsense-dimensional vector particlei. The best particle of the swarm, i.e., the particle with the highest fitness value, is denoted by Gparticle, and its fitness value is denoted as Gbest. The best previous fitness value and corresponding position of the particleiis recorded and represented as and Pparticlei and the position change of the ith particle is denoted as velocityi. The particles are manipulated according to the following equations:
where χ is inertia coefficient in PSO algorithm and the variables ω1and ω2 are random positive numbers, drawn from a uniform distribution within interval [0,1]. An extensively employed function of χ, c1, and c2 is presented in :
where c=4.1,χc1=χc2=1.149445, and so χ≈0.729 are recommended.
Therefore, the PSO-based resource allocation can be described as follows:
Algorithm 2 PSO-based adaptive resource allocation algorithm
Step 1: Initialization for i=1 to S sensingSpectrum(i) end for calculate |Ωk| for all k according to (28); generate base particle according to Algorithm 1 generate particle and shadow particle according to (30) (31) Step 2: Subchannel Allocation for i=1 to Npso calculate the fitness(i) according to (32) end for renew Gbest, Gparticle, and Pparticlei∀irenew velocityiand particlei∀iaccording to (33) (34) if satisfy the stop conditions go to Step 3 else go to the top of Step 2; end if Step 3: Power Allocation Transform particle according to Section “Max–Min–SU”; for i=1 to Npso calculate the fitness(i) according to (32) end for renew Gbest, Gparticle , and Pparticlei∀irenew velocityiand particlei∀iaccording to (33) (34) if satisfy the stop conditions stop algorithm else go to the top of Step 3; end if
Computational complexity analysis
According to the analysis above, the computational complexity of proposed algorithm can be segmented into four parts: initialization, initialization topology of swarm, subchannel allocation, and power allocation, where subchannel allocation and power allocation can be further segmented into adapting and updating parts.
For initial stage, the algorithm should calculate the achievable rate of each SU in every subchannel. This cost complexity of O(K×Nsense). After that, sorting algorithm expenses O(K2) complexity, and in each sorting, the algorithm also needs to calculate O(K×Nsense) computation unit. Hence, the total complexity in initial topology is O(K2×K×Nsense). In the stage of adapting of subchannel, the proposed algorithm needs to calculate fitness value of each particle, which cost totally O(K×Nsense×Npso) for all particles. Then, in order to renew the Gbest and Pbest for all particles, comparison operations with complexity O(Npso) is needed. At the end of one searching in subchannel, a particle renewing with the complexity of O(Nsense×Npso) is executed. The complexity analysis for power allocation is the same as subchannel allocation and the corresponding results for each step is listed in Table 1.
Table 1. Computational complexity upper bounds of four algorithms
Simulation results and discussion
In this section, we present the simulation results to show the performances of the proposed resource allocation algorithm.
To evaluate the performance of proposed algorithm in heterogeneous environment, we consider a system consisting of one CRN and three ambient PRNs as shown in Figure 8. A CBS is located in the center of the map with the coverage radius d=1000 m. Three PBSs are arranged around the CBS as vertex of isosceles triangle, like a typical cell’s structure. The distance between CBS and some PBS is 2d=2000 m. We select the long-term evolution (LTE) as the background PRNs. According to the proposal of IEEE 1900, LTE network will adopt the DSA technique, and all radios are equipped with cognitive functionality in the near future [35,36].
Figure 8. Simulation topology for the infrastructure-based CRN with multiple primary networks.
In all simulations presented in this section, the COST-231 Hata path loss model is considered .
where hBSand hMSare the height of base station and mobile station, respectively. Without loss of generality, we assume hBS=32 m and hMS=1.5 m; Fi is the central transmission frequency. In this article, three PRNs are assumed to be licensed in LTE potential band, with central frequency F1=1805 MHz, F2=1930 MHz and F3=2110 MHz, respectively. The shadowing is implemented by lognormal distribution with standard deviation values of 6 dB. Each PRN share 20 MHz spectrum with CRN in overlay mode. We assume that the subchannel spacing is approximate to 180 kHz, which is equal to one resource block (RB) in LTE network and hence the maximal available subchannels in each PRN is about 110(19.8 MHz). For small-scale fading, we adopt Clarke’s flat fading model with six independent Rayleigh multipaths the same as in . The power delay profile is assumed exponentially decaying with e−2l, where l is the multipath index. Hence, the relative power of the six multipath components are [0,−0.869,−17.37,−26.06,−34.74,−43.43].
We assume that the maximal passband for SU is 200 MHz, which usually needs approximate 1000 MSPS or higher A/D converter in engineering implementation. The total transmission power of CRN on different PRNs is dependent on the results of negotiation between PRN and CRN. How to evaluate the utility of spectrum trading is out of our scope, but in order to reflect differences, we assume that the power limits in PRNs is different. Specifically, we adopt W, W and W, respectively. All other simulation parameters are list in Table 2.
Table 2. Simulation parameter list
Results and discussion
Comparisons with different objective functions
Figure 9 shows the performance comparisons of three objective functions discussed in Section “Objective function”, with K=20 and Nsense=200. Figure 9a,b illustrates the convergence process in subchannel allocation and power allocation respectively, and Figure 9c indicates the final capacity allocation of each SU. In Figure 9, we can see that the MAX–MIN–SU scheme performs worst in terms of proportional fairness, minimum user’s capacity and sum capacity. The reason lies in that the objective of the MAX–MIN–SU is unique, which results in that the search process is sensitive to the initial conditions and prone to be trapped in a local minimum. The MPF scheme has the best proportional fairness (Φ≥0.98) and the convergence process is smooth, both in subchannel allocation and power allocation. As we analyze in Section “Objective function”, the MPF objective is the Pareto optima for problem (9), hence the search process converges faster.
Figure 9. Performance comparisons with three objective functions.
As shown in Figure 9, the proposed hybrid scheme outperforms the other two schemes in terms of both minimum user’s capacity and sum capacity. This is because that the hybrid scheme combines the advantages of the MPF and the MAX–MIN–SU schemes. Specifically, the proposed scheme possesses a medium convergence ability with slight fluctuations (as shown in Figure 9a), which enables the search process jumping out the local traps with high probability; furthermore, with the consideration of maxminRk, the hybrid scheme can also keep the quality of solution in a high level. Figure 9c shows the final resource allocation results with three schemes. We can see that the hybrid scheme has a higher capacity than the other two schemes for most SUs. Therefore, the proposed hybrid objective function can provide better system performance and user experience.
In addition to the comparisons of different objective functions, two other results should be noticed as well. First, we can see that the MPF scheme outperforms the hybrid scheme in terms of proportional fairness while the hybrid scheme outperforms the MPF schemes in terms of capacity. This result reflects the tradeoff between the proportional fairness and the sum capacity, which is consistent with the conclusion in . Second, as shown in Figure 9b, we can see that the improvement of system performance is limited by power allocation (less than 10%). Therefore, we can remove the power allocation process with little performance loss, if the system is subject to the computational complexity constraint.
Comparisons with the conventional algorithm
To the best of the authors’ knowledge, there is no research jointly considering such constraints simultaneously. One similar study was presented in . The authors of  proposed an algorithm to solve the resource allocation problem considering proportional fairness. However, they ignored the bandpass limit of users in . To ensure the algorithm in  suitable for broadband spectrum environment, we have to modify it with the constraint of bandpass limit of SU. We assume that each SU in the algorithm in  can operate simultaneously over two PRNs, in other words, the maximal frequency distance for each SU is 220 MHz in our simulation scenario. Then, we divide all K SUs into two parts, SUs from 1 to K/2 are grouped as “first half” and SUs from K/2 + 1 to K are grouped as “second half”. Each PRN includes SUs as shown in Figure 10. Based on the modification, we can allocate the subchannels and power in each PRN using Zukang’s algorithm independently.
In Figure 11, we show the minimum SU’s capacity versus the number of SU with the proposed algorithm and the algorithm in . The spectrum utilization refers to the spectrum occupancy percentage of primary systems. Specifically, the spectrum utilization 50% means the available idle subchannels is (1−0.5)N≈165 and the spectrum utilization 40% means the available idle subchannels is (1−0.4)N≈198. As shown in Figure 11, the proposed algorithm achieves significant capacity gain over the algorithm in  under the constraint of bandpass limit. The reason is that in the proposed algorithm, SUs can be allocated to subchannels more flexibly than that in the algorithm in . In addition, we can observe that the minimum capacity decreases with the number of available subchannels. This is because low power is allocated to each subchannel under a given maximum power as the number of subchannel increases.
Figure 11. Minimum SU’s capacity versus K and spectrum utilization.
However, as shown in Figure 12, we note that achievable transmission rate actually increases with the number of available subchannel for SUs. For example, when K=8, the total transmission rate with Nsense=198 (spectrum utilization is 0.4) is 12.16 Mbits, while the total rate is 11.01 Mbits with Nsense=165 (spectrum utilization is 0.5).
Figure 12. Minimum SU’s achievable rate versus K and spectrum utilization.
Figure 13 shows the comparisons of proportional fairness with the two algorithms. It can be seen that the proposed algorithm outperforms the algorithm in  in terms of proportional fairness, with Φ>0.98 for all K. This is because the proposed algorithm can allocate subchannels and power to SUs in a global viewpoint, whereas Zukang’s algorithm is constrained by the division of SU groups. Figure 14 shows the capacity distribution among SUs for one channel realization. It also shows the phenomenon that the capacity is unbalanced between the first-half group and the second-half group. In addition, we can see that with the increase in the number of SUs, the proportional fairness of CRN system decrease sightly. In the light of (5), the more the number of SUs is ongoing in the system, the harder the condition Φ=1 is achieved.
Figure 13. Proportional fairness versus K and spectrum utility.
Figure 14. Capacity distributions for one channel realization.
Impact by spectrum sensing errors
The imperfect performance of a spectrum sensing algorithm is also a significant effect on the usage of the spectrum resources for CRNs. In this section, we take an energy detector for example to illustrate the impact of the spectrum sensing errors on the transmission rate for SUs. The Pf and Pd can be expressed as 
where m is the number of sampling points, γis the detector threshold and Q(·) is the standard Gaussian complementary CDF. P is the received average signal power and σ2 represents the noise variance, and thus the average SNR is P/σ2.
Figure 15 shows the achievable transmission sum rate (R·Nsense·W) versus Pf, where “dash line” represents the sum rate calculated by the ARA algorithm. In practice, the CBS cannot judge which subchannels have been occupied by PRNs because of sensing errors. Therefore, the sum rate is calculated without considering the interference from PRNs. “solid line” represents the allocation results with the interference from PRNs. As shown in Figure 15, the gap between the sum rate calculated by the CBS and the actual achievable sum rate decreases with Pf. The reason is when Pf increases, we can achieve lower Pm, resulting in the calculation error caused by PRNs interference reduces. This means that the ARA algorithm is more credible with lower Pm.
Figure 15. Achievable transmission sum rate versus SNR and Pf.
In addition, for a given SNR, with the increase in Pf, the number of reduces. Hence, when Pfis small, the CRN gains higher achievable transmission sum rate. However, it is at the expense that more subchannels occupied by PRNs are miss detected. Therefore, there exists an optimal Pf to balance the tradeoff between the system performance of CRN and interference on PRNs. How to obtain the solution of the optimal Pfis out of our scope of this article, and we just define a simple utility function U as an example to illustrate the existence of the tradeoff.
The utility function U consists of two parts, where ωsuNidle(1−Pf) is reward part and ωpu·NbusyPmis penalty part. This means that when CRN senses correctly an idle subchannels, he would be rewarded and when CRN senses wrongly a busy subchannel, he would be penalized. ωsu and ωpu represent the importance of CRN and PRN, respectively.
Figure 16 shows U versus Pf with different SNR. It is seen that for a given SNR, there exists an optimal Pf such that U is maximal. When SNR is at a low level, the CRN has to set a higher Pf to ensure Pm is low enough, whereas if SNR is high, CRN can select a lower Pf to obtain more transmission opportunities.
Figure 16. Utility function U versus SNR and Pf.
In this article, we develop an adaptive resources allocation scheme for CRNs in multiple PRNs environment. An MINP formulation is established to describe the scenario where a CRN exists around multiple PRNs. To obtain an adaptive resources allocations algorithm with polynomial time complexity, we design a hybrid objective function and two penalty functions, transforming the MINP into two unconstrained optimizations. Then, a PSO-based heuristic algorithm with low complexity initial topology is presented. Three computer simulations are carried out to verify the performance of our proposal. The results reflects that the proposed algorithm is efficient in terms of convergence and quality. Through comparisons, we show that the proposed algorithm outperforms traditional ARA algorithm both in sum capacity and proportional fairness.
The authors declare that they have no competing interests.
This study was supported by the National Basic Research Program of China (973 Program) (Grant No. 2009CB320402), and the National Natural Science Foundation of China (Grant No. 61032003 and 61001092).
J Mitola, GQ Maguire, Cognitive radio: making software radios more personal. IEEE Personal Commun. Mag 6(4), 13–18 (1999). Publisher Full Text
K Cumanan, R Krishna, Z Xiong, S Lambotharan, Multiuser spatial multiplexing techniques with constraints on interference temperature for cognitive radio networks. IET Signal Process 4(6), 666–672 (2010). Publisher Full Text
J Jang, KB Lee, Transmit power adaptation for multiuser ofdm systems. IEEE J. Sel. Areas Commun 21(2), 171–178 (2003). Publisher Full Text
CY Wong, R Cheng, K Lataief, R Murch, Multiuser ofdm with adaptive subcarrier, bit, and power allocation. IEEE J. Sel. Areas Commun 17(10), 1747–1758 (1999). Publisher Full Text
X Dong, J Wang, Y Zhang, M Song, R Feng, End-to-end qos provisioning in future cognitive heterogeneous networks. Proceeding of IEEE Conference on Communications Technology and Applications (ICCTA), 425–429 (2009)
D Kumar, S Mahalaxmi, J Kumar, R Ramya, Adaptive resource allocation for real-time services in ofdma based cognitive radio systems. Kaleidoscope: Beyond the Internet-Innovations for Future Networks and Services ITU-T, 1–5 (2010)
IF Akyildiz, W-Y Lee, MC Vuran, S Mohanty, Next generation/dynamic spectrum access/cognitive radio wireless networks: a survey 50(13), 2127–2159 ([Online], 2006), . Available at http://www.sciencedirect.com/science/article/pii/S1389128606001009 webcite
WJ Tang, QH Wu, Biologically inspired optimization: a review. Trans. Inst. Meas. Control 31(6), 495–515 (2009). Publisher Full Text
T Renk, C Kloeck, D Burgkhardt, FK Jondral, D Grandblaise, S Gault, JC Dunat, Bio-inspired algorithms for dynamic resource allocation in cognitive wireless networks. Proceeding of International Conference on 2nd Cognitive Radio Oriented Wireless Networks and Communications (CrownCom), 351–356 (2007)
C Rieser, T Rondeau, C Bostian, T Gallagher, Cognitive radio testbed: further details and testing of a distributed genetic algorithm based cognitive engine for programmable radios. Proceeding of IEEE Conference on Military Communications (MILCOM), vol. 3, 1437–1443 (2004)
D Maldonado, B Le, A Hugine, T Rondeau, C Bostian, Cognitive radio applications to dynamic spectrum allocation: a discussion and an illustrative example. Proceeding of IEEE Symposium on New Frontiers in Dynamic Spectrum Access Networks (DySPAN), 597–600 (2005)
TR Newman, R Rajbanshi, AM Wyglinski, JB Evans, GJ Minden, Population adaptation for genetic algorithm-based cognitive radios. Proceeding of International Conference on 2nd Cognitive Radio Oriented Wireless Networks and Communications (CrownCom), 279–284 (2007)
S Udgata, K Kumar, S Sabat, Swarm intelligence based resource allocation algorithm for cognitive radio network. Proceeding of International Conference on 1st Parallel Distributed and Grid Computing (PDGC) (Solan (HP), 2010), pp. 324–329
J Perez-Romero, O Salient, R Agusti, L Giupponi, A novel on-demand cognitive pilot channel enabling dynamic spectrum allocation. Proceeding of IEEE Symposium on New Frontiers in Dynamic Spectrum Access Networks (DySPAN), 46–54 (2007)
R Fletcher, An exact penalty function for nonlinear programming with inequalities. Math. Program 5, 129–150 (10, 1973), . 1007/BF01580117. [Online]. Available at http://dx.doi.org/10.1007/BF01580117 webcite Publisher Full Text
J-M Yang, Y-P Chen, J-T Horng, C-Y Kao, Applying family competition to evolution strategies for constrained optimization. in Evolutionary Programming VI, ser, ed. by . Lecture Notes in Computer Science, vol. 1213 (Springer, 1997), pp. 201–211 (10, 1997), . 1007/BFb0014812. [Online]. Available at http://dx.doi.org/10.1007/BFb0014812 webcite
A Khosla, S Kumar, K Ghosh, A comparison of computational efforts between particle swarm optimization and genetic algorithm for identification of fuzzy models. Fuzzy Information Processing Society (NAFIPS), 245–250 (2007)
M Clerc, J Kennedy, The particle swarm—explosion, stability, and convergence in a multidimensional complex space. IEEE Trans. Evol. Comput 6(1), 58–73 (2002). Publisher Full Text
S Buljore, H Harada, S Filin, P Houze, K Tsagkaris, O Holland, K Nolte, T Farnham, V Ivanov, Architecture and enablers for optimized radio resource usage in heterogeneous wireless access networks: the IEEE 1900.4 working group. IEEE Commun. Mag 47(1), 122–129 (2009)
D Gomez-Barquero, A Bria, J Monserrat, N Cardona, Minimal cost planning of dvb-h networks on existing wireless infrastructure. Proceeding of IEEE Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC), 1–5 (2006)