SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Open Badges Research

Capacity-maximization threshold design for wideband sensing with guaranteed minimum primary user rate

Peng Jia and Tho Le-Ngoc*

Author Affiliations

Department of Electrical and Computer Engineering, McGill University, Montreal, QC, H3A 0E9, Canada

For all author emails, please log on.

EURASIP Journal on Wireless Communications and Networking 2012, 2012:360  doi:10.1186/1687-1499-2012-360

The electronic version of this article is the complete one and can be found online at: http://jwcn.eurasipjournals.com/content/2012/1/360

Received:30 August 2011
Accepted:7 November 2012
Published:6 December 2012

© 2012 Jia and Le–Ngoc; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


We investigate the design of optimal threshold for energy detection of primary user (PU) beacon signals in spectrum sensing. The weighted sum of the PU and secondary user (SU) link rates is used as the objective function with constraint on the guaranteed minimum PU link rate. SU with wideband sensing capability selects the idle frequency-slot with the lowest sensed power to be the candidate to access. Numerical results show that the wideband sensing capability with the proposed optimal threshold and accessing strategy can help the SU to meet the stringent PU rate constraints while achieving much higher SU link rate and network capacity.

Cognitive radio; Spectrum sensing; Capacity; Wideband sensing

1. Introduction

Cognitive radio (CR) [1] is promising to leverage the spectrum holes which hide in the grossly underutilized licensed band due to the inflexible spectrum management [2,3]. To achieve this, the secondary user (SU) must be aware of their surrounding environments. Usually, SU performs spectrum sensing to obtain information of the licensed band usage. In this process, different techniques can be used while the simplest one is energy detection which requires no pre-knowledge of the primary user (PU) signal. However, one of the biggest challenges in energy detection is how to design the sensing threshold. Our earlier work [4] investigated the performance of different optimal sensing thresholds derived according to both classic Bayesian and capacity-maximization sensing objectives, indicating that in the communication systems, optimizing miss detection and false alarm probabilities is just the intermediate step while the network capacity should be optimized. Bayesian method is found to be more PU link rate protective while unconstrained maximization of network capacity may jeopardize the PU link rate to a much higher degree.

To maintain the performance of PU link, in this article,a we include the constraint on the guaranteed minimum PU rate in capacity-maximization threshold design for cognitive sensing. Furthermore, we consider SU with wideband sensing capability and PU transmission of imbedded pilot signals, so that SU can sense PU activity in a number of frequency slots of the licensed PU frequency band and select the suitable idle frequency slot for SU transmission. In this way, the SU may have more than one choice, and an appropriate selection algorithm can help to enhance the sensing reliability, which can in turn protect better the PU link rate while improving the achievable SU rate and sum capacity. From the viewpoint of reducing miss detection, we propose the SU to select the frequency slot with the lowest sensed power to be the candidate to access. Based on this wideband sensing approach, we further formulate the network capacity optimization problem with the guaranteed PU link rate, and derive the corresponding optimum threshold. The merits of the wide-band sensing and the corresponding access algorithm are firmly supported by the simulation results as compared to the narrow-band sensing case.

The remainder of this article is organized as follows: Section 2 introduces the single PU and single SU spatial network model and the wideband sensing channel model which includes both large- and small-scale fadings. Furthermore, the network capacity under the wideband operation is derived in Section 3. Section 4 defines the sensing objective function and the constraint function, and we also derive the optimal sensing threshold in this section as well. In Section 5, numerical results are illustrated to show the benefits of the wideband sensing and the spectrum access algorithm. Finally, Section 6 concludes this article.

2. Network and channel model

2.1. Network model

Figure 1 illustrates the network model under consideration. The SU transmitter (ST) location is used as the reference point, i.e., the origin of the polar coordinates (0, 0). The SU receiver (SR) is assumed to be uniformly located in a circular area centered at the ST with radius DS where ds≤ DS denotes the distance between ST and SR. ST can introduce noticeable interference to the PU receiver (PR) if the ST–PR distance, dSP, is within the impact distance DI. The PU transmitter (PT) lies in a circular area centered at the PR with radius DP and dp≤ DP denotes the PT–PR distance. dPS* is the distance between the PT and ST. We also assume a minimum allowable transmitter–receiver distance of ε.

thumbnailFigure 1. Network model.

2.2. Sensing channel model

Consider a segment of the licensed PU operation frequency band of B Hz. This segment can have contiguous or non-contiguous spectra with a total bandwidth of W Hz (where W < B), which can be divided into K equal frequency slots. PU transmission in each frequency slot of Δf = W/K Hz is assumed to be a statistically independent event with a priori probability (or PU transmission activity factor) λ. The channel between any Tx–Rx pair of distance d over the wide frequency band of B Hz, including large-scale path loss and small-scale multipath fading, can be represented by the following channel impulse response

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M1">View MathML</a>


where Δt = 1/B, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M2">View MathML</a>, d is the distance between the transmitter and receiver, α is the path loss exponent, L is the number of resolvable paths and the tap coefficient of the lth resolvable path, hFl, is assumed to be a complex circular Gaussian random variable having statistically independent real and imaginary components with zero mean and variance <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M3">View MathML</a>

2.3. Sensing structure

In each frequency slot, PU sends a beacon (e.g., a sinewave tone or a preamble known by SU), which can be regarded as deterministic signal, to assist SU to detect the spectrum occupation (e.g., [5]). We consider the SU having wideband sensing capability so that it can sense the PU beacon signals in a number of frequency slots to select one vacant slot for transmission as follows. To resolve the wideband channel into a number of narrowband frequency slots, the SU uses the sensing structure as shown in Figure 2. The wideband sample is first taken as input of the N-point Fast Fourier transform (FFT) to obtain the sensed beacon signals in K subcarriers of the interested frequency segment (where K < N and Δf = B/N). The power values, ρk’ s, of the sensed subcarrier beacons are then measured using energy detection and the minimum value, i.e., ρj = mink=1,2,.,Kρk is selected for further comparison with the optimal threshold to finally decide on the status (i.e., vacant or occupied) of frequency slot j.

thumbnailFigure 2. The sensing structure.

In the frequency domain, via discrete Fourier transform, the discrete frequency response in kth frequency slot is

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M4">View MathML</a>


From the characteristics of hls, it can be shown that

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M5">View MathML</a>


Then, for each frequency slot, the ST needs to perform a binary hypothesis testing to decide between the following two hypotheses:

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M6">View MathML</a>


<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M7">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M8">View MathML</a> and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M9">View MathML</a> denote the absence and presence of PU transmission in the kth frequency slot, respectively. Xk is the beacon signal transmitted by PU with power SPB and Vk ~ CN(0, σ02) is the complex Gaussian thermal noise in the kth frequency slot. It follows that the sensed power sample in the kth frequency slot, |Yk|2, is exponentially distributed with the probability density function

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M10">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M11">View MathML</a> and recall that <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M12">View MathML</a> is the distance between the PT and ST.

3. Network and capacity formulation

Consider the PU transmission in K frequency slots with equal probability λ and equal power SP. The SU keeps sensing every frequency slot and chooses the minimum sampled power to compare with the optimal threshold. If the sampled power is smaller than the predefined threshold, then the SU will perceive the corresponding frequency slot which is idle and initiate its communication. However, in this cognitive process, the SU will not be able to perfectly detect the idle/occupied frequency slot. Since miss detection will bring interference to the PU and false alarm will waste SU’s precious transmission opportunities, they will both have the impact on the network capacity. For a threshold γ used in the selected frequency slot, the conditional probabilities of miss detection Pm(γ) and false alarm Pf(γ) can be derived as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M13">View MathML</a>



<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M14">View MathML</a>


For the sensing structure and spectrum accessing strategy shown in Figure 2, and independently and exponentially distributed |Yk|2, k = 1, 2, …, K, the probability that the minimum power sample is in the frequency slot k* can be derived as <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M15">View MathML</a>. In a given time, the PU may occupy q frequency slots, where 0 ≤ qK with the probability of Pq = λq(1 − λ)(K − q). In the q occupied frequency slots, rj = (σ02 + σX2)−1, while in the (K – q) vacant frequency slots, rk = σ0−2. Therefore, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M16">View MathML</a>. It follows that, for a given situation with q occupied frequency slots and (K – q) vacant frequency slots, the conditional probability that the minimum power sample is in an occupied frequency slot is <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M17">View MathML</a>, and the conditional probability that the minimum power sample is in a vacant frequency slot is <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M18">View MathML</a>. Subsequently, the probability that the minimum power sample is in the area of occupied frequency slots can be expressed as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M19">View MathML</a>


Similarly, the probability that the minimum power sample is in the area of vacant frequency slots can be expressed as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M20">View MathML</a>


It is interesting to verify that Po + Pv = 1, and for K = 1, Po = λ, Pv = 1 − λ.

When miss detection occurs, the SU transmits and hence introduces interference to the PU in only one occupied frequency slot with the probability of Pm(γ)Po/K. In this event, the PU rate in the presence of SU interference is <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M21">View MathML</a> (in b/s/Hz) where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M22">View MathML</a> is the PU signal power received at the PR with SPK = KSP, <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M23">View MathML</a>, IS = SSh02dSPα|hFSP|2 represents the interference from SU to PR, and σ02 is the thermal noise power in a certain frequency slot. Otherwise, in the absence of SU interference, with the probability of [λPm(γ)Po/K], the PU rate in the presence of SU interference is <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M24">View MathML</a> (in b/s/Hz). In other words, the PU link rate of a certain frequency slot can be represented as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M25">View MathML</a>


where E{.} represents the expectation over fading and all PU and SU radio locations in the network.

When the false alarm occurs, SU will not transmit (while the frequency slot is actually vacant), effectively introducing capacity loss. This event occurs with the probability of Pf(γ)Pv. In other words, the event that frequency slot is actually vacant and the SU correctly detects the absence of the PU with the probability of [1 − Pf(γ)]Pv, and, in this case, the SU can transmit in the vacant frequency slot with the rate of <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M26">View MathML</a> (in b/s/Hz) where LS = SSh02dSα|hFS|2 is the SU signal power received at the SR. On the other hand, in the event of miss detection, SU transmits in an occupied frequency slot in the presence of interference from the PU with the probability of Pm(γ)Po and the rate of <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M27">View MathML</a> (in b/s/Hz) where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M28">View MathML</a> represents the interference from PU to SR with SPK = KSP, and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M29">View MathML</a>. Therefore, the SU link rate can be expressed as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M30">View MathML</a>


4. Objective function and optimal threshold

We aim at maximizing the weighted sum of the achieved PU and SU rates while guaranteeing the minimum PU rate, i.e., the optimization problem can be written as Maximize

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M31">View MathML</a>


subject to

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M32">View MathML</a>


where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M33">View MathML</a> (in b/s/Hz) represents the achieved PU rate in the absence of SU, and 0 ≤ μ ≤ 1 is selected to emphasize the relative importance between the PU and the SU link rates. Substituting Equations (11) and (12) into (13), the objective function becomes <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M34">View MathML</a> where

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M35">View MathML</a>

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M36">View MathML</a>

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M37">View MathML</a>


Taking the first derivative of Equation (13), <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M38">View MathML</a> and letting f'(γ) = 0, we get <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M39">View MathML</a>

The steady point is the maximum if it satisfies the following relation

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M40">View MathML</a>

Since b is non-positive, only a needs to be discussed. For a > 0 and f "(γ) < 0, the optimal threshold can be written as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M41">View MathML</a>


For a ≤ 0, f(γ) would be maximum as γ, corresponding to the case in which ST transmits without monitoring the PU transmission status.

Furthermore, the constraint (14) implies that <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M42">View MathML</a> where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M43">View MathML</a> (in b/s/Hz) denotes the actual PU capacity loss and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M44">View MathML</a> (in b/s/Hz) denotes the target PU capacity loss due to interference from the SU in K frequency slots. Since <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M45">View MathML</a> and ΔCPK are both non-negative, when ΔCPK φ ≤ 0, the constraint on the guaranteed target minimum PU rate is always satisfied, i.e., SU can use the optimal threshold derived in (15). Moreover, by solving the inequality above, the constraint on the threshold can be derived

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M46">View MathML</a>

In summary, the capacity maximization optimal threshold which guarantees a certain PU link rate can be written as

<a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M47">View MathML</a>


It is worth noting that, under the setting of φ = 1 (i.e., the guaranteed minimum PU rate is set exactly at its achievable rate in the absence of the SU), ΔCPKφ = ΔCPK, and γ = 0, indicating that the SU is never allowed to transmit.

5. Illustrative results

5.1. Simulation method

The network model in Figure 1 is used to set up the simulation. The locations of SU and PU’s transmitters and receivers are generated in the same way as the full location case introduced in [4]. However, the channel realization is more complicated for the wideband channel model. At first, a K × 3000K binary matrix (allocation matrix) is generated for each radio location set of K subcarriers, and 3000K channel realizations. In each row of this matrix, the proportion of ones and zeros is equal to λ/(1 − λ), but the positions of the ones and zeros are uniform randomly distributed from row to row. This binary matrix is used to represent the occupied/vacant situation of the licensed band.

For each radio location set, two L × 3000K matrices of samples of <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M48">View MathML</a> and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M49">View MathML</a> are generated first with 3000K realizations where L is the number of resolvable paths and is set to 4 in the simulation. The matrices of hFPl and hFPSl are processed by FFT to obtain samples of and <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M50">View MathML</a><a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M51">View MathML</a>, which are organized as two K × 3000K matrices. Two K × 3000 matrices of LS, Is samples are also generated. After performing E{.} and setting μ = 0.5, a and b can be calculated to derive the optimal threshold. Another K × 3000K random matrix is generated in the same way as LPK, IPK but with the channel <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M52">View MathML</a>, where <a onClick="popup('http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://www.jwcn.eurasipjournals.com/content/2012/1/360/mathml/M53">View MathML</a> for each row. By dot-multiplying HK by the allocation matrix then adding thermal noise, the detected signal sample matrix YK is obtained. Next, to calculate Pf(γ)Pv and Pm(γ)Po, two counters are used. For each column of YK, the minimum value is compared with the corresponding threshold γ. If it is larger than γ and corresponds to a vacant frequency slot (i.e., element in the allocation matrix with value of 0), then the false alarm counter is increased by one. Similarly, if it is smaller than γ and the corresponding value in the allocation matrix is 1, then the miss detection counter is increased by one. Eventually by dividing the false alarm counts by the number of zeros experienced across each column, Pf(γ)Pv is obtained. Pm(γ)Po is taken as the ratio of the miss detection counts to the number of 1’s experienced across each column. For each radio location set, the same procedure is performed. Eventually, CPU, CSU are calculated according to Equations (11) and (12).

5.2. Results and discussion

Figures 3, 4, and 5 plot the PU, SU link rates, and sum rate (in b/s/Hz) versus the PU activity factor λ, respectively, for constraint of 0.8CPUo where CPUo is also included in Figure 3 for reference. While the PU link rate is within the permitted range (as shown in Figure 3), the results plotted in Figure 4 indicate that the constraint on guaranteed minimum PU rate largely decreases the achievable SU link rate in the narrow-band case (i.e., one frequency slot), as compared to the wideband sensing scheme (i.e., four or eight frequency slots). As a result, the wideband sensing and spectrum access strategy maintain much higher SU link rate and eventually deliver a much more overall network capacity (sum rate) as shown in Figure 5.

thumbnailFigure 3. PU rates versus PU activity factor for different numbers of frequency slots with φ = 0.8.

thumbnailFigure 4. SU rates versus PU activity factor for different numbers of frequency slots with φ = 0.8.

thumbnailFigure 5. Sum rates versus PU activity factor for different numbers of frequency slots with φ = 0.8.

Moreover, it is interesting to notice that it exhibits a trend that the more candidate subcarriers SU can choose from, the better SU link rate can be achieved. This is confirmed numerically for a stringent requirement of guaranteed PU minimum rate with φ = 0.95. The corresponding PU, SU, and sum rates (in b/s/Hz) versus the PU activity factor λ are plotted in Figures 6, 7, and 8, respectively. The best achievable PU rate CPUo (in the absence of SU) is also plotted in Figure 6 for reference. Comparing the results in Figures 4 and 7 indicates that, when λ increases, the achieved SU rate is decreased to maintain the guaranteed PU minimum rate. Higher guaranteed PU minimum rate reduces the achieved SU rate. For the guaranteed PU rate increase of 15% (from 0.8CPUo in Figure 4 to 0.95CPUo in Figure 7), the achieved SU rate drops by more than 50% for narrowband sensing (i.e., single frequency slot), and less than 20% for wideband sensing (for four or eight frequency slots). The gain in the achieved SU rate due to wideband sensing is more pronounced at higher PU activity factor λ.

thumbnailFigure 6. PU rates versus PU activity factor for different numbers of frequency slots with φ = 0.95.

thumbnailFigure 7. SU rates versus PU activity factor for different numbers of frequency slots with φ = 0.95.

thumbnailFigure 8. Sum rates versus PU activity factor for different numbers of frequency slots with φ = 0.95.

Since the objective of our optimization problem is to maximize the weighted sum rate while guaranteeing the minimum PU rate (CPUφCPUo, 0 ≤ φ < 1), it can be expected that the sum rate is non-decreasing with λ as shown in Figures 5 and 8. However, at high value of λ (approaching 1), the achieved sum rate becomes saturated. This trend indicates that CR operation with SU sharing the PU spectrum band is more efficient at low PU activity.

6. Conclusion

In this article, the design of optimum sensing threshold for capacity maximization while guaranteeing the required minimum PU link rate is investigated for SU with wideband sensing capability. SU selects the idle frequency slot to access corresponding to the minimum measured PU beacon power exceeding the optimum sensing threshold. The wideband sensing capability along with the proposed access strategy helps the SU to reduce the miss detection probability and to increase its achievable link rate. Illustrative results show that, compared to the narrowband sensing, the wideband sensing approach offers a much better SU link rate and network sum capacity over the whole range of PU transmission activity factor and the resulting SU and sum rates increase monotonically with the number of sensed frequency slots while meeting stringent constraints on the guaranteed minimum PU rates.


aPart of this article was presented at the IEEE VTC2011-Fall Conference in San Francisco, September 5–8, 2011.

Competing interests

The authors declare that they have no competing interests.


  1. S Haykin, Cognitive radio: brain-empowered wireless communications. IEEE. J. Sel. Areas Commun. 23(2), 201–220 (2005)

  2. FCC, Spectrum policy task force report ET Docket. Technical Report, 02–155 (2002)

  3. FCC, Federal communications commission spectrum policy task force FCC Report of the Spectrum Efficiency Working Group (FCC, Technical Report, 2002)

  4. P Jia, M Vu, T Le-Ngoc, S Hong, V Tarokh, Capacity- and Bayesian-based cognitive sensing with location side information. IEEE. J. Sel. Areas Commun. 29(2), 276–289 (2010)

  5. N Maziar, A survey of cognitive radio access to tv/white spaces. Int. J. Digital Multimed. Broadcast 11 Article ID 236568 Publisher Full Text OpenURL