Next Article in Journal
An Anti-Electromagnetic Attack PUF Based on a Configurable Ring Oscillator for Wireless Sensor Networks
Next Article in Special Issue
Achieving Congestion Mitigation Using Distributed Power Control for Spectrum Sensor Nodes in Sensor Network-Aided Cognitive Radio Ad Hoc Networks
Previous Article in Journal
Ultra-Sensitive NT-proBNP Quantification for Early Detection of Risk Factors Leading to Heart Failure
Previous Article in Special Issue
Joint Resource Optimization for Cognitive Sensor Networks with SWIPT-Enabled Relay
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Optimal Cooperative Sensing with Energy Detection in Cognitive Radio

Department of Electronic Engineering, Sogang University, Seoul 04107, Korea
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(9), 2111; https://doi.org/10.3390/s17092111
Submission received: 2 August 2017 / Revised: 10 September 2017 / Accepted: 11 September 2017 / Published: 15 September 2017
(This article belongs to the Special Issue Cognitive Radio Sensing and Sensor Networks)

Abstract

:
In this paper, we propose an optimal cooperative sensing technique for cognitive radio to maximize sensing performance based on energy detection. In most spectrum sensing research, many cooperation methods have been proposed to overcome the sensitivity of energy detection so that both primary and secondary users are better off in terms of spectral efficiency. However, without assigning a proper sensing threshold to each sensing node, cooperation may not be effective unless the received average primary user signal-to-noise ratio (SNR) is identical. We show that equal threshold energy detection severely degrades sensing performance when cooperative sensing nodes experience diverse average SNRs, and it is not unusual for even single-node sensing to be better than cooperative sensing. To this end, based on the Neyman–Pearson criterion we formulate an optimization problem to maximize sensing performance by using optimized thresholds. Since this is a non-convex optimization problem, we provide a condition that convexifies the problem and thus serves as a sufficient optimality condition. We find that, perhaps surprisingly, in all practical cases one may consider this condition satisfied, and thus optimal sensing performance can be obtained. Through extensive simulations, we demonstrate that the proposed technique achieves a globally optimal solution, i.e., it maximizes the probability of detection under practical operating parameters such as the target probability of false alarm, different SNRs, and the number of cooperative sensing nodes.

1. Introduction

We are entering the era of a fourth industrial revolution where both hyper-connectivity and digital automation enable a large number of different types of economic entities to share information quickly for optimal decision-making based on industrial IoT sensor networks. This wave has driven 5G New Radio (NR) development to provide greater capacity and higher reliability for tactile internet, artificial intelligence (AI), and smart factory applications, etc. Since intelligent and efficient ways of utilizing scarce spectrum are required to meet the huge demand for increasing capacity, new, emerging technologies for dynamic spectrum management have been developed to improve spectral efficiency.
Cognitive radio (CR) is a promising form of technology for combating spectrum scarcity in next-generation wireless networks. A group of secondary users (SUs) can exploit the unused spectrum bands of a primary user (PU) in an opportunistic manner as long as no harmful interference is guaranteed [1]. Recent worldwide migration from analog to digital television broadcasting has expedited sharing of underutilized spectrum bands called TV white space (TVWS) [2]. For example, standards such as IEEE 802.22 and IEEE 802.11af opportunistically utilize TVWS, and regulatory bodies such as the FCC and OFCOM have updated rules for TVWS and other unlicensed bands [3]. TVWS-based CR was also introduced to enable smart grid communications in rural areas because of the good propagation properties of TVWS [4]. LTE-CR was proposed as a cellular application to extend the LTE time-division duplex (TDD) to support TVWS [5]. In addition to TVWS, 3GPP has developed LTE in an unlicensed spectrum so that LTE can harmoniously coexist with other wireless systems in a 5-GHz spectrum [6].
To support CR, one of the key roles of SUs is to sense the PU with high detection probability and evacuate the PU’s band immediately. There are several detection techniques, such as energy detection, matched filter detection, cyclostationary detection, and compressive sensing, [7]. Among these techniques, energy detection of [8] has been widely adopted in practice because it has low complexity and does not require prior knowledge of the PU [9]. However, since sensing performance of energy detection can deteriorate when the signal-to-noise ratio (SNR) is low [10], many cooperative sensing techniques have been proposed to overcome performance degradation [11,12,13,14,15,16,17,18,19,20,21,22,23,24].
Specifically, to maximize sensing performance, a series of cooperation algorithms have been investigated in terms of sensing duration, the number of cooperative sensing nodes, and thresholds of energy detection [11,12,13,14,15,16,17,18]. In [11], sensing-throughput tradeoff depending on sensing duration was investigated. The author of [12] jointly optimized sensing time and reporting time to maximize throughput based on maximal ratio combination (MRC)-based soft decision. The fundamental performance limits of the cooperative sensing using energy detection were analyzed by considering the unlimited number of sensing nodes, and it has been shown that the OR rule achieves zero false alarm probability for any given target detection probability irrespective of the received PU SNR [13]. In [14], the authors proposed optimal counting rules of hard decision under the Neyman–Pearson criterion and Bayesian criterion, respectively. In [15], both sensing duration and fusion rule were jointly optimized to maximize the throughput of SUs under the constraint of the probability of detection. In [16,17], an optimal k was investigated for the so called k-out-of-N rule, and then the same threshold for all sensing nodes is optimized by a given k. The authors of [18] proposed an on/off reporting mechanism to achieve robust cooperative sensing over an imperfect reporting channel as well as graceful degradation against sensor failures.
However, the existing work [11,12,13,14,15,16,17,18] assumed that the average SNR is identical for all cooperative sensing nodes, and so is the sensing threshold. This assumption may not be realistic in practical systems. In addition to multipath fading and shadowing effects, there are several practical factors for heterogeneous cooperative sensing, such as geographically different positions (indoor or outdoor), inherent RF front-end sensitivity, obstacles, temperature variations, and so on. Thus, cooperative sensing performance highly depends on the heterogeneity of the received PU SNRs at individual SUs. Although previous work, such as [19], considered heterogeneous average SNRs for cooperative sensing nodes, all sensing nodes still used the same sensing threshold.
In addition, the authors in [20,21,22,23,24] considered optimizing thresholds of cooperative sensing nodes based on different average SNRs. In [20], an iterative threshold selection algorithm was proposed to determine whether a specific sensing node can participate in cooperation or not, based on its contribution to detection performance. The work of [21] proposed a weighted decision rule by jointly optimizing the test threshold of likelihood-ratio and the local threshold of SUs. Although heuristic methods of optimizing two thresholds at a time did not guarantee the achievement of global optimal solutions, the proposed algorithm provided the highest throughput results over the non-jointly optimized schemes [21]. In [22,23], a heuristic algorithm was proposed to achieve suboptimal thresholds at cooperative sensing nodes by minimizing false alarm probability subject to target missed detection probability. Different fusion rules with iterative algorithms were proposed to decrease the total error detection probability for both single channel and multichannel cooperative sensing [24].
However, even though the existing literature [20,21,22,23,24] has investigated distinct thresholds for cooperative sensing nodes according to their own average SNRs, they are suboptimal and theoretically achieving the optimal detection performance still remains unsolved. Moreover, previous heuristic approaches would suffer from high computational complexity as the number of cooperative sensing nodes increases [20,21,22,23,24]. In this regard, we are motivated to theoretically develop a globally optimal solution for energy detection under a heterogeneous PU SNR environment.
We summarize our key contributions as follows. Based on the Neyman–Pearson criterion we formulate an optimization problem to maximize the sensing performance when a hard decision is used under different average SNRs. It is found that this is a non-convex optimization problem, and thus we provide a condition that convexifies the original problem. We then verify that all practically meaningful parameters satisfy the convexifying condition. Thus, our solution guarantees the maximized detection performance of cooperative sensing. Considering that using a single common threshold severely degrades the sensing performance when cooperative sensing nodes experience different average SNRs, our solution is promising. Our extensive simulations show that a globally optimized solution achieves an improvement of several orders of magnitude improvement under practical operating parameters such as the target probability of false alarm, different average SNRs, and the number of cooperative sensing nodes.
The rest of this paper is organized as follows. In Section 2, the CR system model and problem formulation are introduced. Cooperative sensing with optimal thresholds is presented with a sufficient optimality condition in Section 3. We present numerical results under practical operating parameters in Section 4. Finally, we draw our conclusions in Section 5.

2. System Model

In this section, we provide a scenario of cooperative sensing operation in the CR system. Using energy detection and its sensing performance metrics, we formulate an optimization problem.

2.1. Cognitive Radio System

Figure 1 illustrates a scenario for CR system with cooperative sensing in wireless networks. A SU can temporally and spatially utilize the licensed or unlicensed spectrum band of the PU. To avoid possible interference to the PU, the CR system requires SUs to simultaneously perform in-band sensing, which means SUs periodically sense the operating spectrum band to verify whether the PU appears or not. In addition to in-band sensing, out-of-band sensing implies that SUs sense the non-operating bands and maintain a list of available bands. Even though our work is applicable to both in-band and out-of-band sensing purposes, hereafter we consider in-band sensing. When SUs finish sensing, they send their own sensing results to a fusion center so that it combines all local results and makes a final decision about the presence of PUs. We assume that SUs honestly report their local decisions to the fusion center via common control channels [11,12,13,14,15,16,17], e.g., an individual reporting channel is reliable and error-free for simplicity [19,20,21,22,23,24].

2.2. Energy Detection and Sensing Performance Metrics

Let τ be a continuous value denoting sensing duration. Under the absence of PU (hypothesis H 0 ), the received signal at a sensing node is given by
r ( t ) = n ( t ) ,
where n ( t ) is the noise power at time t. Under the presence of PU (hypothesis H 1 ), we have
r ( t ) = h ( t ) · s ( t ) + n ( t ) ,
where s ( t ) is the received PU signal and h ( t ) is the channel gain at time t. We assume that n ( t ) is white Gaussian noise with two-sided power spectral density N 0 (AS I), and sensing duration τ is short so that the channel can be assumed to be static (AS II). According to [8], the test statistic using energy detection is given by
T = 1 N 0 0 τ | r ( t ) | 2 d t .
Let γ be the received SNR given by γ = P N 0 B , where P denotes the received power of PU at a sensing node, and B is the channel bandwidth. Under the assumptions of AS I and AS II, the received energy of PU signal is simply P τ . When PU is absent, T has the central chi-square distribution with the degrees of freedom equal to 2 τ B . When PU is present, T has the non-central chi-square distribution with 2 τ B degrees of freedom and non-centrality parameter δ = P τ N 0 equal to γ τ B . Based on the central limit theorem (CLT), when 2 τ B is more than 250, T under two conditions can be considered as Gaussian random variables [8].
T N 2 τ B , 4 τ B , under   H 0 , N 2 τ B + γ τ B , 4 τ B + 4 γ τ B , under   H 1 ,
where N m , σ 2 represents the normal distribution with the mean m and the variance σ 2 . From now on, we assume that 2 τ B is large enough so that the CLT condition is met, and thus T is a Gaussian random variable.
Suppose that we have N cooperative sensing nodes using energy detection. Let ϵ i be the detection threshold and P i d be the probability of detection, i.e., P i d = P Y ϵ i | H 1 at a sensing node i 1 , , N  [11,13]. Then, we have
P i d ( ϵ i , γ i , τ , B ) = Q ϵ i 2 τ B + γ i τ B 4 τ B + 4 γ i τ B ,
where
Q ( t ) = 1 2 π t exp u 2 2 d u .
Let P i f denote the probability of false alarm, i.e., P i f = P Y ϵ i | H 0 at a sensing node i 1 , , N  [11,13], and we have
P i f ( ϵ i , τ , B ) = Q ϵ i 2 τ B 4 τ B .
We notice that there is no PU signal present under H 0 , and then P i f is independent of γ i . However, when h ( t ) is varying due to shadowing/fading, P i d gives the detection probability conditioned on the instantaneous SNR γ i . In this case, the average detection probability can be derived by averaging (2) over fading statistics.
P ˜ i d = Q ϵ i 2 τ B + γ i τ B 4 τ B + 4 γ i τ B f γ ( γ i ) d γ i ,
where f γ ( γ i ) is the probability distribution function (PDF) of SNR under fading/shadowing and its mean value is γ ¯ i  [25]. In addition to the aforementioned probabilities, the probability of miss detection is P i m = 1 P i d , and the probability of correct rejection is P i c = 1 P i f at a sensing node i 1 , , N . During cooperative sensing, a PU is assumed to be either present or absent. Although multiple PUs can arrive or depart randomly, this assumption is still reasonable based on the long time average scenario of observed PU behavior. One may refer to [26,27] to optimize both the sensing time and sensing period (duty cycle) so that instantaneous sensing performance improves under random arrivals of multiple PUs in cooperative sensing.

2.3. Problem Formulation of Cooperative Sensing

We investigate optimal sensing thresholds for sensing nodes using hard decision. Even though use of the soft decision achieves better performance than the hard decision, the burden of reporting overhead hinders the use of soft decision [28]. One may want to directly consider the k-out-of-N rule, but finding an globally optimal solution for the k-out-of-N rule is not very tractable unless the same average SNR value is assumed for cooperative sensing scenarios [16]. Moreover, [21,22,23] pointed out the very high computational complexity and impracticality of finding suboptimal k. Thus, we first focus on AND and OR rules in this paper. The fusion center using the AND rule decides that PU is present when all cooperative sensing nodes report PU presence to the fusion center [9]. Then, the final detection and the final false alarm probabilities at the fusion center are simply given by
P d = i = 1 N P i d ,
P f = i = 1 N P i f .
Based on the Neyman–Pearson criterion, our aim is to maximize the final detection probability P d subject to the target final false alarm probability P f ¯ , which is formulated as follows:
max ϵ i , i 1 , , N P d = i = 1 N P i d , s . t . P f = i = 1 N P i f P f ¯ .
The fusion center using the OR rule decides that the PU is present when any cooperative sensing nodes report PU presence to the fusion center [9]. Then, the final detection and the final false alarm probabilities at the fusion center are simply
P d = 1 i = 1 N 1 P i d ,
P f = 1 i = 1 N 1 P i f .
Then, we can also use the OR rule to maximize the final detection probability P d subject to the target final false alarm probability P f ¯ as follows:
max ϵ i , i 1 , , N P d = 1 i = 1 N 1 P i d , s . t . P f = 1 i = 1 N 1 P i f P f ¯ .
Considering the symmetric mathematical structure of P i m and P i c , an alternative optimization problem is formulated as follows:
min ϵ i , i 1 , , N P m = i = 1 N P i m , s . t . P c = i = 1 N P i c P c ¯ ,
where P c ¯ is according given by 1 P f ¯ .
If the globally optimal solution can be derived from the AND rule of (6), our optimization approach is also applicable to the OR rule because of its symmetric mathematical structure. AND and OR rules are two opposite extremes. For comparison, the existing literature [7,11,13,28] have shown that dominant performance regions of both rules highly depend on system parameters such as τ , γ ¯ , B , and P f ¯ for different purposes. Herein we mainly focus on the existence of globally optimized sensing thresholds. Therefore, in the rest of paper, our problem formulation of cooperative sensing mainly exploits the mathematical structure of the AND rule so that we can provide a globally optimal solution to maximize the detection performance of cooperative sensing under different average SNRs.

3. Cooperative Sensing with Optimal Thresholds

3.1. Limitation of Equal Threshold—A Motivating Example with Two Sensing Nodes

To better understand the effect of sensing thresholds, in this section we first verify the detection performance with both equal and optimal thresholds when only two cooperative sensing nodes experience different average SNRs.
Example 1.
Since the same average SNR for all sensing nodes in [11,13,14,15,16,17] may not be realistic in practice because of spatial diversity, herein we use different average SNRs in cooperative sensing [19,20,21,22,23,24]. As a motivating example, Figure 2 shows the receiver operating characteristic (ROC) curves under four different scenarios: the sensing performance of a single node sensing with the worst and the best average SNRs (CASE I and CASE II) as well as that of both equal and optimal thresholds under different average SNRs (CASE III and CASE IV). To obtain the equal threshold ϵ for two cooperative sensing nodes [19], the P i f at individual sensing node is simply given by P f ¯ from (6), and then ϵ is accordingly given by 4 τ B Q 1 P f ¯ + 2 τ B from (4). Interestingly, compared to a single node sensing with a higher SNR ( 12  dB), the sensing performance of cooperative sensing with equal threshold is even worse, and thus cooperation is not beneficial at all. To overcome the degradation of the detection performance, we assign the optimal threshold to each sensing node based on an exhaustive search for illustrative purposes. By doing so, optimal thresholds can improve the sensing performance. Since an exhaustive search is not suitable for a large number of sensing nodes, in Section 3.2, we further develop a convex optimization algorithm to achieve optimal thresholds. When many sensing nodes cooperate, we are also interested in knowing how much cooperative gain is expected by using globally optimized sensing thresholds.
Example 2.
In Figure 3, we verify how the difference of average SNRs affects the detection performance of cooperative sensing. In this example with two sensing nodes, we vary the average SNR of one sensing node and fix that of the other sensing node. For the equal threshold case, we notice that the effective range of the average SNR is very limited. However, using optimal thresholds alway outperforms other cases regardless of the average SNR difference. To this end, along with the previous work [19,20,21,22,23,24], we are motivated to find out the globally optimal solution rather than suboptimal solutions. We will see that cooperative sensing dramatically (by several orders of magnitudes) improves sensing performance under practical operating parameters such as the target probability of false alarm, different average SNRs, and the number of cooperative sensing nodes.

3.2. Cooperative Sensing with Optimal Thresholds

In solving (6), the individual sensing threshold for sensing node i is ϵ i = 4 τ B Q 1 P i f + 2 τ B by (4). Then, we can use P i f as optimization variables instead of ϵ i , and rewrite the optimization problem (6) by using (2) for analytical tractability:
max P i f , i 1 , , N i = 1 N Q 1 1 + γ i Q 1 P i f γ i τ B 2 , s . t . i = 1 N P i f P f ¯ .
With one further step, we define the normalized threshold x i ϵ i 2 τ B 4 τ B from H 0 of (1). Then, it can be easily drawn from (4) such that x i = Q 1 P i f .
After taking the logarithm of (7), the problem is equivalent to the following optimization problem:
max x i , y i , i 1 , , N i = 1 N log Q y i , s . t . y i = 1 1 + γ i x i γ i τ B 2 , i = 1 N log Q x i log P f ¯ ,
where y i is introduced as an auxiliary variable. The convex problem has three requirements: the objective function must be concave for the maximization problem; the equality constraint functions must be affine; and the inequality constraint convex functions must be non-positive [29]. Note that although the objective function of (8) is concave, (8) still remains a non-convex optimization problem. This is because when we make the inequality constraint functions convex, it becomes non-negative which cannot satisfy a standard form of convex optimization. Thus, the problem cannot be posed as convex optimization [29]. Hence, we now try to find the condition that can convexify the problem as presented in the following proposition.
Proposition 1.
Suppose that observations of cooperative sensing nodes follow the normal distribution as in (1) with different SNR γ i . Then, based on the Neyman–Pearson criterion, (8) can be convexified if
R x i + x i 1 1 + γ i R y i + y i , i 1 , , N ,
where a function R ( t ) is defined as R ( t ) Q ( t ) Q ( t ) .
Proof. 
See Appendix A. ☐
Remark 1.
Note that Proposition 1 serves as a sufficient optimality condition. As long as the sufficient optimality condition is satisfied, the problem can be convexified so that a globally optimal solution can be obtained.
We notice that the sufficient optimality condition may not always hold for x i R . However, by considering all the practically meaningful scenarios, we rigorously verify that as long as P i f is on [ ε , 1 ε ] where ε is infinitesimally small and sensing observations follow the normal distribution, the sufficient optimality condition holds so that a globally optimal solution is guaranteed in cooperative sensing, as described in Examples 3 and 4.
Example 3.
To investigate the practicality of Remark 1, in Figure 4 we plot both the left-hand side (LHS) and the right-hand side (RHS) of (9) with τ = 1  ms and B = 6  MHz for various PU SNR γ i . As can be seen, LHS is greater than RHS, and thus (9) holds for x i [ 10 , 10 ] . Recall that x i is a normalized version of ϵ i that determines the false alarm probability, i.e., P i f = Q x i and x i has a region of interest depending on P i f . Ideally, P i f should be any value on [ 0 , 1 ] . We notice that when (9) holds, the corresponding P i f covers almost all values on ( 0 , 1 ) . In this example, P i f is in the interval [ ε , 1 ε ] where ε 7 . 620 × 10 24 , and thus it is obvious that a set of x i satisfying (9) is a super set of the practical region of interest. In Figure 5, we also show that (LHS-RHS) is positive in (9). Note that although the gap between LHS and RHS depends on the range of SNR, the inequality condition itself still holds as long as P i f is on [ ε , 1 ε ] regardless of SNR. For example, even for the worst case of zero SNR, the condition (9) holds because LHS becomes equal to RHS.
Example 4.
One may want to verify practicality of Remark 1 with other values of τ and B. Thus we also consider the worst case condition such that 2 τ B is small and just satisfies the minimum CLT condition of (1), e.g., 2 τ B 250. We confirm that when x i 5 , (9) also holds and the corresponding P i f spans in the interval ( 0 , 1 ε ] where ε 2 . 867 × 10 7 . Thus, although the problem is a non-convex optimization, it can be convexified with Proposition 1 in practice, and a globally optimal solution can be derived as in the following subsection.

3.3. Convex Optimization

Under the region of satisfying (9), now we can obtain the globally optimal solution by applying the standard convex optimization procedure. Let v = ( v 1 , , v N ) , and the Lagrangian function of (A1) is
L ( v , λ ) = i = 1 N log Q a i ( log Q ) 1 ( v i ) γ ˜ i + λ log P f ¯ i = 1 N v i ,
where λ is the Lagrange multiplier. The KKT condition for this problem consists of four parts, which are a primal feasible, dual feasible, complementary slackness, and zero gradient:
i = 1 N v i * log P f ¯ 0 , λ * 0 , λ * i = 1 N v i * log P f ¯ = 0 , Q a i ( log Q ) 1 ( v i * ) γ ˜ i Q ( log Q ) 1 ( v i * ) a i Q a i ( log Q ) 1 ( v i * ) γ ˜ i Q ( log Q ) 1 ( v i * ) λ * = 0 , i = 1 , , N .
Since there is one Lagrange multiplier λ and also v i can be expressed by a function of λ according to (10), one can obtain the optimal ( v * , λ * ) by adjusting λ based on some iterative methods such as bisection method and Newton’s method. The details of the convex optimization algorithm are shown in Algorithm 1 based on the bisection method.
Algorithm 1: Find the globally optimal v * and λ * for (10).
Input: Given N, τ , B, P f ¯ , and γ i for all i.
  1:  Set initial λ min and λ max .
  2:  Calculate a i and γ ˜ i for all i.
  3:  while | λ min λ max | > ε do
  4:  λ = ( λ min + λ max ) / 2 .
  5: Given λ , solve (10) to compute v i for all i.
  6: If i = 1 N v i log P f ¯ > 0 , then set λ min = λ , else set  λ max = λ .
  7:  end while.
Output: Calculate P i f (and also ε i ), and obtain the maximum P d in (6) and (7).

4. Numerical Results

In this section we present the numerical results with optimal cooperative sensing. For numerical results, the sensing duration τ is set at 1 msec and the channel bandwidth B is set at 6 MHz, both of which satisfy the test statistic of the sensing node in (1). The extensive numerical results are mainly based on the AND rule to verify the effect of optimal thresholds on both miss detection and false alarm performances so that no harmful interference as well as high spectral efficiency are guaranteed for PU and SUs, respectively. As aforementioned in Section 2.3, one may easily verify that our approach can be extended to the OR rule in a similar manner because of its symmetric mathematical structure as with the AND rule. Liu et al. [24] proposed a local optimal solution using iterative-threshold method, which is effective when SNRs are different. When average received SNRs are same for all cooperative sensing nodes, their sensing thresholds should be same [24]. However, when SNRs are same, the method in [24] cannot achieve the performance of the uniform threshold method because of its inherent limitation of iteration procedure. By contrast, our proposed algorithm can achieve a globally optimal solution whether SNRs are same or not. For the identical SNR case, we also verify that the uniform threshold method exactly matches our globally optimal solution.
In Figure 6, ROC curves are presented for both a single node sensing and cooperative sensing with seven nodes, i.e., N = 7 . Each cooperative sensing node experiences its own average received SNR γ ¯ i , which is randomly selected in the interval [ 10 , 16 ] dB. We also set the γ ¯ i of a single node sensing to be either 10 or 16 dB to represent the best and the worst case scenarios, respectively. As can be seen, although cooperative sensing with equal threshold (-⋄-) is obviously better than that of a single node sensing with γ ¯ i = 16 dB (-□-), it is interesting to observe that cooperative sensing (-⋄-) can be worse than single node sensing with γ ¯ i = 10 dB (-∘-). Namely, cooperation does not necessarily guarantee better sensing performance unless the detection thresholds are properly set. However, cooperative sensing with optimal thresholds (-✡-) significantly improves the sensing performance. Specifically, when compared to cooperative sensing with equal thresholds, under the target P f ¯ is 0.01 , and P m is reduced from 30% to 2% (and thus the corresponding P d is improved from 70% to 98%).
Next, to investigate the impact of N on optimal cooperative sensing performance, we compare the miss detection probabilities P m of the equal threshold method and the proposed optimal threshold in Figure 7. As can be seen, the more sensing nodes participate in cooperation, the lower the achieved P m is, in both cases. However, the slope of optimal threshold is much steeper than that of equal threshold. Furthermore, the miss detection probability is reduced by more than 100 times when N is 30 or more. Since the improvement is roughly proportional to N, we conclude that using optimized thresholds is essential for a large scale cooperative sensing, e.g., without optimal thresholds, many cooperative sensing nodes possibly lead to stagnant sensing performance.
In Figure 8, we investigate the impact of the mean of different average SNRs E [ γ ¯ i ] on P m . When it is relatively low (below 17 dB), we do not observe a noticeable difference in P m between equal and optimal thresholds. However, as it grows, optimal thresholds can achieve much lower P m than that of equal thresholds. For example, a miss detection probability as low as 10 5 times can be achieved when E [ γ ¯ i ] is 10 dB.

5. Conclusions

In this paper we proposed an optimal cooperative sensing solution using energy detection when sensing nodes experience different PU SNRs. Specifically, the fusion center and its cooperative sensing nodes used the hard decision under the Neyman–Pearson criterion. Despite the non-convexity of the problem, we provided a sufficient optimality condition that interestingly holds in the practical region of interest, and thus a globally optimal solution was obtained. Furthermore, utilizing optimal thresholds provides a gain of several orders of magnitude of cooperation in terms of miss detection probability compared with the conventional equal threshold under various scenarios of the target probability of false alarm, the number of cooperative sensing nodes, and the different average SNRs. Some practical extensions of the proposed scheme may include topics such as imperfect reporting channels, reduced reporting overheads, and energy-efficient sensing techniques in more general hard and soft decisions.

Acknowledgments

This research was supported in part by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science and ICT under Grant (NRF-2017R1A1A1A05001377).

Author Contributions

Sunghwan Bae designed the algorithm, performed the simulations, and prepared the manuscript. Jaewoo So contributed to the design of the system and supervised the activities. Hongseok Kim led the reasearch project and supervised the activities as the corresponding author. All authors discussed the results and approved the publication.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof of Proposition 1.
For notational simplicity we define a i 1 1 + γ i and γ ˜ i γ i τ B 2 . Since (8) is not convex optimization, we again define v i log Q x i to make the inequality constraint of v i form a convex set. Using x i = ( log Q ) 1 ( v i ) , the optimization in (8) becomes
max v i , i 1 , , N i = 1 N log Q a i ( log Q ) 1 ( v i ) γ ˜ i , s . t . i = 1 N v i log P f ¯ .
Here, we define a function f ( t ) as
f ( t ) Q a i ( log Q ) 1 ( t ) γ ˜ i .
Then, for the objective function of (A1), the pointwise sum i = 1 N log f ( v i ) is concave if f ( v i ) is a log-concave function of v i for all i [29]. The log-concavity property of f ( v i ) needs to hold for all i, so we now omit the subscript i for simplicity. To check whether f ( v ) is log-concave or not, we use the log-concave condition, e.g., the inequality is f ( v ) f " ( v ) f 2 ( v ) [29].
Although one may directly obtain the first and second derivatives of f ( v ) with respect to v, it is quite difficult to deal with ( log Q ) 1 ( v ) in f ( v ) . Interestingly, we notice that it is possible to express both f ( v ) and f " ( v ) as a function of x. Once we use a function defined as z ( t ) Q ( a ( t γ ˜ ) ) , then f ( v ) = z ( x ( v ) ) = Q ( a ( x ( v ) γ ˜ ) ) , where x ( v ) = ( log Q ) 1 ( v ) . To obtain f ( v ) and f " ( v ) , based on the chain rule d f ( v ) d v = d z ( x ( v ) ) d v = d z ( x ) d x · d x d v as well as the derivative formula ( u 1 ) ( v ) = 1 u ( u 1 ( v ) ) for any analytic function u : R R [30], we have
d x d v = Q ( x ) Q ( x ) ,
which is only a function of x. Thus, we have f ( v ) as follows:
f ( v ) = a Q a x γ ˜ Q ( x ) Q ( x ) .
In a similar way, we also have f " ( v ) as follows:
f " ( v ) = a 2 Q " a x γ ˜ Q ( x ) 2 Q ( x ) 2 + a Q a x γ ˜ Q ( x ) Q ( x ) a Q a x γ ˜ Q " ( x ) Q ( x ) 2 Q ( x ) 3 .
Applying the log-concave condition f " ( v ) f ( v ) f 2 ( v ) with (A2) and (A3), the following long inequality condition is given by
a 2 Q " a x γ ˜ Q ( x ) 2 Q ( x ) 2 + a Q a x γ ˜ Q ( x ) Q ( x ) a Q a x γ ˜ Q " ( x ) Q ( x ) 2 Q ( x ) 3 · Q a x γ ˜ a 2 Q a x γ ˜ 2 Q ( x ) 2 Q ( x ) 2 .
After several steps of manipulations with Q " ( t ) = t 2 π exp t 2 2 = t Q ( t ) by (3), the inequality condition in (A4) boils down to:
x ( 1 a 2 ) + a 2 γ ˜ + Q ( x ) Q ( x ) a Q a x γ ˜ Q a x γ ˜ .
Finally, using a function R ( t ) = Q ( t ) Q ( t ) , the inequality condition is given by
R ( x ) + x a R a x γ ˜ + a ( x γ ˜ ) .
With notations of x i and y i , we obtain
R x i + x i 1 1 + γ i R y i + y i , i 1 , , N ,
which completes the proof. ☐

References

  1. Haykin, S. Cognitive radio: Brain-empowered wireless communications. IEEE J. Select. Areas Commun. 2005, 23, 201–220. [Google Scholar] [CrossRef]
  2. Forde, T.; Doyle, L. A TV whitespace ecosystem for licensed cognitive radio. Telecommun. Policy 2013, 37, 130–139. [Google Scholar] [CrossRef]
  3. Sum, C.-S.; Villardi, G.P.; Rahman, M.A.; Baykas, T.; Tran, H.N.; Lan, Z.; Sun, C.; Alemseged, Y.; Wang, J.; Song, C. Cognitive communication in TV white spaces: An overview of regulations, standards, and technology. IEEE Commun. Mag. 2013, 51, 138–145. [Google Scholar] [CrossRef]
  4. Ghassemi, A.; Bavarian, S.; Lampe, L. Cognitive radio for smart grid communications. In Proceedings of the 2010 First IEEE International Conference on Smart Grid Communications (SmartGridComm), Gaithersburg, MD, USA, 4–6 October 2010. [Google Scholar]
  5. Xiao, J.; Hu, R.Q.; Qian, Y.; Gong, L.; Wang, B. Expanding LTE network spectrum with cognitive radios: From concept to implementation. IEEE Wirel. Commun. 2013, 20, 12–19. [Google Scholar] [CrossRef]
  6. Ratasuk, R.; Mangalvedhe, N.; Ghosh, A. LTE in unlicensed spectrum using licensed-assisted access. In Proceedings of the 2014 Globecom Workshops (GC Wkshps), Austin, TX, USA, 8–12 December 2014. [Google Scholar]
  7. Akyildiz, I.F.; Lo, B.F.; Balakrishnan, R. Cooperative spectrum sensing in cognitive radio networks: A survey. Phys. Commun. 2011, 4, 772–776. [Google Scholar] [CrossRef]
  8. Urkowitz, H. Energy detection of unknown deterministic signals. Proc. IEEE 1967, 55, 523–531. [Google Scholar] [CrossRef]
  9. Yucek, T.; Arslan, H. A survey of spectrum sensing algorithms for cognitive radio applications. IEEE Commun. Surv. Tutor. 2009, 11, 116–130. [Google Scholar] [CrossRef]
  10. Axell, E.; Leus, G.; Larsson, E.G.; Poor, H.V. Spectrum sensing for cognitive radio: State-of-the-art and recent advances. IEEE Signal Process. Mag. 2012, 29, 101–116. [Google Scholar] [CrossRef]
  11. Liang, Y.-C.; Zeng, Y.; Peh, E.C.Y.; Hoang, A.T. Sensing-throughput tradeoff for cognitive radio networks. IEEE Trans. Wirel. Commun. 2008, 7, 1326–1337. [Google Scholar] [CrossRef]
  12. So, J. Opportunistic reporting-based sensing-reporting-throughput optimization scheme for cooperative cognitive radio networks. KSII Trans. Internet Inf. Syst. 2017, 11, 1319–1355. [Google Scholar]
  13. Bae, S.; Kim, H. Unlimited cooperative sensing with energy detection for cognitive radio. J. Commun. Netw. 2014, 16, 172–182. [Google Scholar] [CrossRef]
  14. Shen, J.; Liu, S.; Zeng, L.; Xie, G.; Gao, J.; Liu, Y. Optimisation of cooperative spectrum sensing in cognitive radio network. IET Commun. 2009, 3, 1170–1178. [Google Scholar] [CrossRef]
  15. Peh, E.C.Y.; Liang, Y.-C.; Guan, Y.L.; Zeng, Y. Optimization of cooperative sensing in cognitive radio networks: A sensing-throughput tradeoff view. IEEE Trans. Veh. Technol. 2009, 58, 5294–5299. [Google Scholar] [CrossRef]
  16. Zhang, W.; Mallik, R.; Letaief, K. Optimization of cooperative spectrum sensing with energy detection in cognitive radio networks. IEEE Trans. Wirel. Commun. 2009, 8, 5761–5766. [Google Scholar] [CrossRef]
  17. Liu, Y.; Zeng, C.; Wang, H.; Wei, G. Energy detection threshold optimization for cooperative spectrum sensing. In Proceedings of the 2010 2nd International Conference on Advanced Computer Control (ICACC), Shenyang, China, 27–29 March 2010. [Google Scholar]
  18. Bae, S.; Kim, H. Robust cooperative sensing with on/off signaling over imperfect reporting channels. IEEE Trans. Ind. Informat. 2016, 12, 2196–2205. [Google Scholar] [CrossRef]
  19. Peh, E.C.Y.; Liang, Y.-C. Optimization for cooperative sensing incognitive radio networks. In Proceedings of the IEEE Wireless Communications and Networking Conference, Kowloon, Hongkong, China, 11–15 March 2007. [Google Scholar]
  20. Teo, D.; Zhong, K.; Ng, B.-C. An iterative threshold selection algorithm for cooperative sensing in a cognitive radio network. In Proceedings of the 2010 IEEE Symposium on New Frontiers in Dynamic Spectrum, Singapore, 6–9 April 2010. [Google Scholar]
  21. Peh, E.C.Y.; Liang, Y.-C.; Guan, Y.L.; Zeng, Y. Cooperative spectrum sensing in cognitive radio networks with weighted decision fusion schemes. IEEE Trans. Wirel. Commun. 2010, 9, 3838–3847. [Google Scholar] [CrossRef]
  22. Lai, J.; Dutkiewicz, E.; Liu, R.-P.; Vesilo, R. Cooperative sensing with detection threshold optimization in cognitive radio networks. In Proceedings of the 2012 International Symposium on Communications and Information Technologies (ISCIT), Gold Coast, Australia, 2–5 October 2012. [Google Scholar]
  23. Lai, J.; Dutkiewicz, E.; Liu, R.-P.; Vesilo, R. Performance optimization of cooperative spectrum sensing in cognitive radio networks. In Proceedings of the 2013 IEEE Wireless Communications and Networking Conference (WCNC), Shanghai, China, 7–10 April 2013. [Google Scholar]
  24. Liu, X.; Jia, M.; Tan, X. Threshold optimization of cooperative spectrum sensing in cognitive radio networks. Radio Sci. 2013, 48, 23–32. [Google Scholar] [CrossRef]
  25. Ghasemi, A.; Sousa, E.S. Collaborative spectrum sensing for opportunistic access in fading environments. In Proceedings of the 2005 First IEEE International Symposium on New Frontiers in Dynamic Spectrum Access Networks (DySPAN 2005), Baltimore, MD, USA, 8–11 November 2005. [Google Scholar]
  26. Tang, L.; Chen, Y.; Hines, E.L.; Alouini, M.-S. Effect of primary user traffic on sensing-throughput tradeoff for cognitive radios. IEEE Trans. Wirel. Commun. 2011, 10, 1063–1068. [Google Scholar] [CrossRef]
  27. Düzenli, T.; Akay, O. A new spectrum sensing strategy for dynamic primary users in cognitive radio. IEEE Commun. Lett. 2016, 20, 752–755. [Google Scholar] [CrossRef]
  28. Chaudhari, S.; Lunden, J.; Koivunen, V.; Poor, H.V. Cooperative sensing with imperfect reporting channels: Hard decisions or soft decisions? IEEE Trans. Signal Process. 2012, 60, 18–28. [Google Scholar] [CrossRef]
  29. Boyd, S.; Vandenberghe, L. Convex Optimization; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  30. Larson, R.; Edwards, B.H. Calculus, 9th ed.; Cengage Learning, Inc.: Boston, MA, USA, 2010. [Google Scholar]
Figure 1. A scenario for cooperative sensing of cognitive radio system in wireless networks; there is a primary user, cooperative secondary users, and a fusion center. SU: secondary user.
Figure 1. A scenario for cooperative sensing of cognitive radio system in wireless networks; there is a primary user, cooperative secondary users, and a fusion center. SU: secondary user.
Sensors 17 02111 g001
Figure 2. Receiver operating characteristic (ROC) curves of cooperative sensing with two sensing nodes; node 1 with 15 dB average signal-to-noise ratio (SNR) and node 2 with 12 dB SNR.
Figure 2. Receiver operating characteristic (ROC) curves of cooperative sensing with two sensing nodes; node 1 with 15 dB average signal-to-noise ratio (SNR) and node 2 with 12 dB SNR.
Sensors 17 02111 g002
Figure 3. The detection performance for a single and two cooperative sensing cases when the average SNR of the sensing node 1 varies from 20 dB to 10 dB and that of the sensing node 2 is fixed at 15 dB. We notice that the effective range of the average SNR is very limited for the equal threshold case.
Figure 3. The detection performance for a single and two cooperative sensing cases when the average SNR of the sensing node 1 varies from 20 dB to 10 dB and that of the sensing node 2 is fixed at 15 dB. We notice that the effective range of the average SNR is very limited for the equal threshold case.
Sensors 17 02111 g003
Figure 4. Under practical operating parameters, the y-axis is used to depict value for both left- and right-hand sides of inequality (9) when x i = Q 1 P i f varies. Correspondingly, the table is used to represent value of P i f . LHS: left-hand side; RHS: right-hand side.
Figure 4. Under practical operating parameters, the y-axis is used to depict value for both left- and right-hand sides of inequality (9) when x i = Q 1 P i f varies. Correspondingly, the table is used to represent value of P i f . LHS: left-hand side; RHS: right-hand side.
Sensors 17 02111 g004
Figure 5. Inequality condition (9) in Proposition 1 holds for practical operating parameters.
Figure 5. Inequality condition (9) in Proposition 1 holds for practical operating parameters.
Sensors 17 02111 g005
Figure 6. Receiver operating characteristic (ROC) curves for a single node sensing and cooperative sensing; different SNRs are considered in cooperation with both an equal threshold and optimal thresholds.
Figure 6. Receiver operating characteristic (ROC) curves for a single node sensing and cooperative sensing; different SNRs are considered in cooperation with both an equal threshold and optimal thresholds.
Sensors 17 02111 g006
Figure 7. The probability of miss detection for two different thresholds methods are shown with the increasing number of cooperative sensing nodes; the target P f ¯ = 1 % .
Figure 7. The probability of miss detection for two different thresholds methods are shown with the increasing number of cooperative sensing nodes; the target P f ¯ = 1 % .
Sensors 17 02111 g007
Figure 8. The probability of miss detection for two different thresholds methods are shown when the mean of different average SNRs of cooperative sensing nodes increases; the target P f ¯ = 1 % and N = 7 .
Figure 8. The probability of miss detection for two different thresholds methods are shown when the mean of different average SNRs of cooperative sensing nodes increases; the target P f ¯ = 1 % and N = 7 .
Sensors 17 02111 g008

Share and Cite

MDPI and ACS Style

Bae, S.; So, J.; Kim, H. On Optimal Cooperative Sensing with Energy Detection in Cognitive Radio. Sensors 2017, 17, 2111. https://doi.org/10.3390/s17092111

AMA Style

Bae S, So J, Kim H. On Optimal Cooperative Sensing with Energy Detection in Cognitive Radio. Sensors. 2017; 17(9):2111. https://doi.org/10.3390/s17092111

Chicago/Turabian Style

Bae, Sunghwan, Jaewoo So, and Hongseok Kim. 2017. "On Optimal Cooperative Sensing with Energy Detection in Cognitive Radio" Sensors 17, no. 9: 2111. https://doi.org/10.3390/s17092111

APA Style

Bae, S., So, J., & Kim, H. (2017). On Optimal Cooperative Sensing with Energy Detection in Cognitive Radio. Sensors, 17(9), 2111. https://doi.org/10.3390/s17092111

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop