# Goodness-of-Fit Based Secure Cooperative Spectrum Sensing for Cognitive Radio Network

Authors: Hiep Vu-Van; Insoo Koo
Journal: The Scientific World Journal
 (2014)

---

## Abstract

Cognitive radio (CR) is a promising technology for improving usage of frequency band. Cognitive radio users (CUs) are allowed to use the bands without interference in operation of licensed users. Reliable sensing information about status of licensed band is a prerequirement for CR network. Cooperative spectrum sensing (CSS) is able to offer an improved sensing reliability compared to individual sensing. However, the sensing performance of CSS can be destroyed due to the appearance of some malicious users. In this paper, we propose a goodness-of-fit (GOF) based cooperative spectrum sensing scheme to detect the dissimilarity between sensing information of normal CUs and that of malicious users, and reject their harmful effect to CSS. The empirical CDF will be used in GOF test to determine the measured distance between distributions of observation sample set according to each hypothesis of licensed user signal. Further, the DS theory is used to combine results of multi-GOF tests. The simulation results demonstrate that the proposed scheme can protect the sensing process against the attack from malicious users.


---

## Body

## 1. Introduction

Nowadays, more bandwidth and higher bit-rates have been required to meet usage demands due to an explosion in wireless communication technology. According to the Federal Communications Commission’s spectrum policy task force report [1], the actual utilization of the licensed spectrum varies from 15% to 80%. In some cases, the utilization is only a small percentage of the total capacity. Cognitive radio (CR) technology [2] has been proposed to solve the problem of ineffective utilization of spectrum bands. Both unlicensed and licensed users, termed the cognitive radio user (CU) and primary user (PU), respectively, operate in CR networks. In CR network, CUs are allowed to access the frequency assigned to PU when it is free. But CU must vacate the occupied frequency when the presence of PU is detected. Therefore, reliable detection of the PU’s signal is a requirement of CR networks.In order to ascertain the presence of a PU, CUs can use one of several common detection methods, such as matched filter, feature, and energy detection [2, 3]. Energy detection is the optimal sensing method if the CU has limited information about PU’s signal (e.g., only the local noise power is known) [3]. In energy detection, frequency energy in the sensing channel is received in a fixed bandwidth W over an observation time window T to compare with the energy threshold and determine whether or not the channel is utilized. However, the received signal power may fluctuate severely due to multipath fading and shadowing effects. Therefore, it is difficult to obtain reliable detection with only one CU. Better sensing performance can be obtained by allowing some CUs to perform cooperative spectrum sensing [4–6].CSS can use some combination methods such as equal gain combination (EGC) and maximum gain combination (MGC) [7] to combine sensing information of all CUs in the network and make a global decision about status of PU signal. Since EGC gives the same weight for all CUs in the network, it is easy to execute but with limited performance. MGC is known as the optimal combination rule. However, it requires information about the SNRs of the sensing channel, which is difficult to obtain in practice. In addition, MGC is sensitive to attack by malicious users who send false sensing data to the fusion center (FC) [8]. The research presented in [8, 9] determined that the presence of a few malicious users can severely reduce the performance of a CSS scheme. Algorithms used to identify the malicious users have been proposed in the studies of [8, 9]. In previous research, a simple technique (i.e., outlier-detection) is used to detect less damage malicious CUs such asalways No oralways Yes CU. In addition, the technique is unable to protect the CSS in the event of a large number of malicious users in the network.In this paper, we utilize multi-goodness-of-fit (GOF) tests to design a robust CSS, in which the event detection technique [10, 11] will be used to provide the combination of different evidence of each type of GOF test which are supported by a particular hypothesis of PU signal. The proposed scheme considers two types of GOF tests, Kolmogorov-Smirnov (KS) and Cramer-von Mises (CM) tests. The proposed scheme can distinguish the sensing information of normal CUs and that of malicious users and reject the harmful effect of malicious user to sensing combination process. Three common types of malicious users includingalways Yes,always No, andopposite are considered in this paper.


## 2. Background

### 2.1. Goodness-of-Fit Test

The GOF test summarizes the discrepancy between the observed samples with theoretical distributions or empirical distributions and the reference distribution. For then independent and identical distributed observation, the sample is first arranged in ascending order such that s
1
≤
s
2
≤
⋯
≤
s
n. The GOF test is used to determine whether or not the samples set was drawn from the same distribution with a cumulative distribution function (CDF) F
0. The testing hypothesis can be formulated as follows:
(1)
F
(
s
)
=
F
0
(
s
)
:
H
o
,
F
(
s
)
≠
F
0
(
s
)
:
H
1
,
where F
(
s
) is the empirical CDF of the sample. It can be calculated as follows:
(2)
F
(
s
)
=
1
n
∑
i
=
1
n
‍
I
{
s
i
≤
s
}
,
where I
{
·
} is the indicator of event {
·
}.There are many types of GOF tests, for instance, Cramer-von Mises (CM), Kolmogorov-Smirnov (KS), AndersonDarling (AD), and Hosmer-Lemeshow (HL) tests. In this paper, we consider two types of GOF tests,CM and KS tests, which can run well with a low number of samples.(1)
Kolmogorov-Smirnov (KS) test: the KS test, which is based on the empirical CDF of the samples set and the reference CDF, can be calculated according to the largest difference of two distributions as follows:(3)
D
KS
=
sup
⁡
{
|
F
(
s
i
)
-
F
0
(
s
i
)
|
:
i
=
1
,
…
,
n
}
,
wheresup
⁡
{
·
} is supremum function, which indicates the greatest element of the set. If the sample comes from distribution F
0
(
x
), then D
KS will converge to 0.
(2)
Cramer-von Mises (CM) test: CM test is used for judging the goodness-of-fit of the sample set’s CDFF
(
s
) and reference distribution’s CDF F
0
(
s
). The test statistic is given by
(4)
D
CM
=
n
∫
-
∞
+
∞
‍
[
F
(
s
)
-
F
0
(
s
)
]
2
d
F
0
(
s
)
and can be approximated as(5)
D
CM
=
1
12
n
+
∑
i
=
1
n
‍
[
2
i
-
1
2
n
-
F
0
(
s
i
)
]
2
.
If this value,D
CM, is larger than the threshold, the hypothesis that the sample data come from the reference distribution F
0 can be rejected.


### 2.2. Combination of Evidence in Dempster-Shafer Theory

Dempster-Shafer (DS) theory was first introduced by Demperster and was later extended by Shafer. This is a potentially valuable tool for the evaluation of risk and reliability in engineering applications when it is not possible to obtain a precise measurement from experiments or when knowledge is obtained from expert elicitation. An important aspect of this theory is the combination of evidence obtained from multiple sources and the modeling of conflict between them.In DS theory [12], a representation of ignorance is provided by assigning a nonzero mass function to hypothesis m, also called the basic probability assignment (BPA), and is defined for every hypothesis A such that the mass value m
(
A
) belongs to the interval [
0,1
] and satisfies the following conditions:
(6)
m
(
ϕ
)
=
0
,
∑
‍
m
(
A
)
=
1
,
A
⊆
Ω
,
where Ω is the frame of discernment, which is a fixed set of q mutually exclusive and exhaustive elements.By assigning a nonzero mass in a compound hypothesis,A
∪
B means that there exists the option to not make a decision between A and B but to leave the formulation in the A
∩
B class. In DS theory, two functions, belief (Bel) and plausibility (Pls), are defined to characterize the uncertainty and support of certain hypotheses. Bel measures the minimum or necessary support, whereas Pls reflects the maximum or potential support for that hypothesis [13]. These two measures, derived from mass values, are defined as a map from a set of hypotheses to interval [
0,1
] as follows:
(7)
Bel
(
A
)
=
∑
B
⊆
A
‍
m
(
B
)
,
Pls
(
A
)
=
∑
A
∩
B
≠
0
‍
m
(
B
)
.The sum of mass functions from different information source,m
j
(
j
=
1,2
,
…
M
), combined with the DS rule is known as the orthogonal sum, which is commutative and associative. The result is a new mass function, m
(
A
k
)
=
(
m
1
⊕
m
2
⊕
⋯
m
d
)
(
A
k
), which incorporates the joints information provided by the sources as follows:
(8)
m
(
A
k
)
=
1
1
-
K
∑
A
1
∩
A
2
⋯
A
d
=
A
k
‍
(
∏
1
≤
j
≤
M
‍
m
j
(
A
j
)
)
,
K
=
∑
A
1
∩
A
2
⋯
A
d
=
ϕ
‍
(
∏
1
≤
j
≤
M
‍
m
j
(
A
j
)
)
,
where K is the measure of conflict between the different sources and is introduced as a normalization factor.


## 2.1. Goodness-of-Fit Test

The GOF test summarizes the discrepancy between the observed samples with theoretical distributions or empirical distributions and the reference distribution. For then independent and identical distributed observation, the sample is first arranged in ascending order such that s
1
≤
s
2
≤
⋯
≤
s
n. The GOF test is used to determine whether or not the samples set was drawn from the same distribution with a cumulative distribution function (CDF) F
0. The testing hypothesis can be formulated as follows:
(1)
F
(
s
)
=
F
0
(
s
)
:
H
o
,
F
(
s
)
≠
F
0
(
s
)
:
H
1
,
where F
(
s
) is the empirical CDF of the sample. It can be calculated as follows:
(2)
F
(
s
)
=
1
n
∑
i
=
1
n
‍
I
{
s
i
≤
s
}
,
where I
{
·
} is the indicator of event {
·
}.There are many types of GOF tests, for instance, Cramer-von Mises (CM), Kolmogorov-Smirnov (KS), AndersonDarling (AD), and Hosmer-Lemeshow (HL) tests. In this paper, we consider two types of GOF tests,CM and KS tests, which can run well with a low number of samples.(1)
Kolmogorov-Smirnov (KS) test: the KS test, which is based on the empirical CDF of the samples set and the reference CDF, can be calculated according to the largest difference of two distributions as follows:(3)
D
KS
=
sup
⁡
{
|
F
(
s
i
)
-
F
0
(
s
i
)
|
:
i
=
1
,
…
,
n
}
,
wheresup
⁡
{
·
} is supremum function, which indicates the greatest element of the set. If the sample comes from distribution F
0
(
x
), then D
KS will converge to 0.
(2)
Cramer-von Mises (CM) test: CM test is used for judging the goodness-of-fit of the sample set’s CDFF
(
s
) and reference distribution’s CDF F
0
(
s
). The test statistic is given by
(4)
D
CM
=
n
∫
-
∞
+
∞
‍
[
F
(
s
)
-
F
0
(
s
)
]
2
d
F
0
(
s
)
and can be approximated as(5)
D
CM
=
1
12
n
+
∑
i
=
1
n
‍
[
2
i
-
1
2
n
-
F
0
(
s
i
)
]
2
.
If this value,D
CM, is larger than the threshold, the hypothesis that the sample data come from the reference distribution F
0 can be rejected.


## 2.2. Combination of Evidence in Dempster-Shafer Theory

Dempster-Shafer (DS) theory was first introduced by Demperster and was later extended by Shafer. This is a potentially valuable tool for the evaluation of risk and reliability in engineering applications when it is not possible to obtain a precise measurement from experiments or when knowledge is obtained from expert elicitation. An important aspect of this theory is the combination of evidence obtained from multiple sources and the modeling of conflict between them.In DS theory [12], a representation of ignorance is provided by assigning a nonzero mass function to hypothesis m, also called the basic probability assignment (BPA), and is defined for every hypothesis A such that the mass value m
(
A
) belongs to the interval [
0,1
] and satisfies the following conditions:
(6)
m
(
ϕ
)
=
0
,
∑
‍
m
(
A
)
=
1
,
A
⊆
Ω
,
where Ω is the frame of discernment, which is a fixed set of q mutually exclusive and exhaustive elements.By assigning a nonzero mass in a compound hypothesis,A
∪
B means that there exists the option to not make a decision between A and B but to leave the formulation in the A
∩
B class. In DS theory, two functions, belief (Bel) and plausibility (Pls), are defined to characterize the uncertainty and support of certain hypotheses. Bel measures the minimum or necessary support, whereas Pls reflects the maximum or potential support for that hypothesis [13]. These two measures, derived from mass values, are defined as a map from a set of hypotheses to interval [
0,1
] as follows:
(7)
Bel
(
A
)
=
∑
B
⊆
A
‍
m
(
B
)
,
Pls
(
A
)
=
∑
A
∩
B
≠
0
‍
m
(
B
)
.The sum of mass functions from different information source,m
j
(
j
=
1,2
,
…
M
), combined with the DS rule is known as the orthogonal sum, which is commutative and associative. The result is a new mass function, m
(
A
k
)
=
(
m
1
⊕
m
2
⊕
⋯
m
d
)
(
A
k
), which incorporates the joints information provided by the sources as follows:
(8)
m
(
A
k
)
=
1
1
-
K
∑
A
1
∩
A
2
⋯
A
d
=
A
k
‍
(
∏
1
≤
j
≤
M
‍
m
j
(
A
j
)
)
,
K
=
∑
A
1
∩
A
2
⋯
A
d
=
ϕ
‍
(
∏
1
≤
j
≤
M
‍
m
j
(
A
j
)
)
,
where K is the measure of conflict between the different sources and is introduced as a normalization factor.


## 3. The Proposed Secure Cooperative Spectrum Sensing Based on GOF Test

There is a definite difference between the CDF of received signal energy of normal CU and that of the malicious users as shown in Figure1. The CDF of the received signal energy of normal CU corresponding to the presence of the PU is always “under” that one corresponding to the absence of the PU. On the contrary, theopposite malicious CU has the CDF corresponding to the presence of PU to be “above” the CDF corresponding to the absence of PU. Thealways Yes andalways No malicious CUs have a similar CDF corresponding to presence and absence of PU. Due to the difference between CDF of normal and malicious CUs, we utilize GOF test to detect the appearance of malicious users in the network, so that their harmful effect can be rejected out of CSS process. Multi-GOF tests including KS and CM tests will be applied for adaptive robust CSS. The DS theory will be used to combine results of multi-GOF tests.The CDF of received signal energy at CU under absence and presence hypothesis of PU signal for (a) normal CU, (b)opposite malicious CU, (c)always Yes malicious CU, and (d)always No malicious CU.
(a)
(b)
(c)
(d)In this paper, we consider a CR network includingN CUs who cooperate to sense the signal from a PU. There are p
<
N malicious CUs appearing in the network which can be classified as three common types:always Yes,always No, andopposite malicious CUs. All CUs use energy detectors to perform spectrum sensing and send their sensing data to the FC through a control channel. Based on the sensing data obtained from the CUs, the FC makes a global decision concerning the presence or absence of the PU signal by using the proposed data fusion scheme. The proposed scheme has 3 steps as follows.Step 1.
All CUs perform spectrum sensing by using energy detection method to determine received signal energyE
j
=
{
e
1
,
j
,
e
2
,
j
,
…
,
e
M
,
j
}, where M is the number of sensing samples that the jth CU takes in the sensing interval.Step 2.
At the FC, GOF test statistics of each CU will be computed according to hypothesis of the PU as given in (11). After that, BPA and final BPA for current sensing data will be estimated based on the “reputation level” of each CU, which is updated from previous sensing interval. Based on final BPA, a global decision rule will be proposed to make global decision about status of PU signal.Step 3.
Update “reputation level” of each CU according to the global decision.The detailed description of each step will be given in the following subsections.


### 3.1. Energy Detection

At the sensing interval for thejth CU, the local spectrum sensing is to decide between the two following hypotheses:
(9)
H
0
:
s
j
(
k
)
=
n
j
(
k
)
,
H
1
:
s
j
(
k
)
=
h
j
p
(
k
)
+
n
j
(
k
)
,
where H
0 and H
1 correspond to the hypothesis of the absence and presence of the PU signal, respectively, h
j denotes the amplitude gain of the channel, s
(
k
) is the signal transmitted from the PU, n
j
(
k
) is the additive white Gaussian noise, and k is index of sensing sample at each sensing interval.A received signal energy of a sensing sample,e
k
,
j, is given as
(10)
e
k
,
j
=
{
|
n
j
(
k
)
|
2
,
H
0
|
h
j
s
(
k
)
+
n
j
(
k
)
|
2
,
H
1
.


### 3.2. BPA Estimation

The GOF test statistics of the current sensing datae
k
,
j
(
k
=
1
,
…
,
M
) of the j
th CU will be calculated according to each hypothesis of PU signal based on (3) and (5) as follows, respectively:
(11)
D
h
,
j
KS
=
sup
⁡
{
|
F
(
e
k
,
j
)
-
F
h
(
e
k
,
j
)
|
:
k
=
1,2
,
…
,
M
}
,
D
h
,
j
CM
=
1
12
M
+
∑
i
=
1
M
‍
[
2
i
-
1
2
M
-
F
h
(
e
k
,
j
)
]
2
,
where h
=
{
0,1
} is index of hypothesis H
h of PU signal, e
k
,
j is the received signal energy of kth sensing sample of the jth CU, F
(
·
) and F
h
(
·
) are empirical CDF of observed sensing sample and CDF of H
h hypothesis of PU, and M is the number of samples for each sensing interval.It is noteworthy that normal CU and malicious CU have different characteristics ofD
1
,
j
t and D
0
,
j
t as shown in Table 1, where t indexes types of GOF tests: KS and CM.Table 1
Reputation ranges according to each type of CUs.
Status of PU
Normal CU
A
l
w
a
y
s
Y
e
s CU
A
l
w
a
y
s
N
o CU
O
p
p
o
s
i
t
e CU
H
1
D
0
,
j
t
≫
D
1
,
j
t
D
0
,
j
t
≫
D
1
,
j
t
D
0
,
j
t
≪
D
1
,
j
t
D
0
,
j
t
≪
D
1
,
j
t
D
1
,
j
t
≈
0
D
1
,
j
t
≈
0
D
0
,
j
t
≈
0
D
0
,
j
t
≈
0
r
0
,
j
t
(
i
)
≫
0
r
0
,
j
t
(
i
)
≫
0
r
0
,
j
t
(
i
)
≪
0
r
0
,
j
t
(
i
)
≪
0
H
0
D
0
,
j
t
≪
D
1
,
j
t
D
0
,
j
t
≫
D
1
,
j
t
D
0
,
j
t
≪
D
1
,
j
t
D
0
,
j
t
≫
D
1
,
j
t
D
0
,
j
t
≈
0
D
1
,
j
t
≈
0
D
0
,
j
t
≈
0
D
1
,
j
t
≈
0
r
1
,
j
t
(
i
)
≫
0
r
1
,
j
t
(
i
)
≪
0
r
1
,
j
t
(
i
)
≫
0
r
1
,
j
t
(
i
)
≪
0Based on the values ofD
1
,
j
t and D
0
,
j
t, we will estimate BPA of current sensing data of each CU and their “reputation level” for robust CSS as follows:
(12)
Δ
1
,
j
t
=
D
0
,
j
t
D
1
,
j
t
+
D
0
,
j
t
R
0
,
j
t
,
Δ
0
,
j
t
=
D
1
,
j
t
D
1
,
j
t
+
D
0
,
j
t
R
1
,
j
t
,
where R
h
,
j
t is “reputation level” of the jth CU according to hypothesis h, and it can be determined based on history observation of the jth CU as follows:
(13)
R
h
,
j
t
(
i
)
=
r
h
,
j
t
(
i
-
1
)
∑
j
r
h
,
j
t
(
i
-
1
)
,
where i is the index of current sensing interval and r
1
,
j
t
(
i
-
1
) and r
0
,
j
t
(
i
-
1
) are updated from the previous sensing interval according to global decision:
(14)
r
1
,
j
t
(
i
-
1
)
=
r
1
,
j
t
(
i
-
2
)
+
(
D
1
,
j
t
(
i
-
1
)
-
D
0
,
j
t
(
i
-
1
)
)
,
(15)
r
0
,
j
t
(
i
-
1
)
=
r
0
,
j
t
(
i
-
2
)
+
(
D
0
,
j
t
(
i
-
1
)
-
D
1
,
j
t
(
i
-
1
)
)
.By usingr
1
,
j
t and r
0
,
j
t, types of CUs will be easily distinguished. The normal CU has positive value of both r
1
,
j
t and r
0
,
j
t that will be increased after updating step. r
1
,
j
t ofalways Yes and r
0
,
j
t ofalways No malicious CUs are almost negative and tend to decrease after updating step. On the other hand, both values of opposite CU are negative and have a tendency to decrease. We define “malicious threshold” as ρ to reject the attack of malicious CR in CSS, so that the CU, which has either r
1
,
j
t
<
0 or r
0
,
j
t
<
0, will be determined as malicious CU. The sensing data of malicious CUs will not be considered to make global decision by giving them r
1
,
j
t
=
0 and r
0
,
j
t
=
0.The BPA of all CUs will be combined with their reputation levels as(16)
Δ
1
t
=
1
n
Ω
∑
j
∈
Ω
‍
D
0
,
j
t
D
1
,
j
t
+
D
0
,
j
t
R
0
,
j
t
,
Δ
0
t
=
1
n
Ω
∑
j
∈
Ω
‍
D
1
,
j
t
D
1
,
j
t
+
D
0
,
j
t
R
1
,
j
t
,
where Ω and n
Ω are set of normal CUs and number of members of the set, respectively.Because the error in estimatingΔ
0
t and Δ
1
t,  Δ
0
t
+
Δ
1
t can be bigger than 1, we need to normalize those values as
(17)
Δ
0
t
*
=
Δ
0
t
Δ
0
t
+
Δ
1
t
,
Δ
1
t
*
=
Δ
1
t
Δ
0
t
+
Δ
1
t
.


### 3.3. DS Theory Combination

The DS theory will be used to combine the BPA of both GOF tests according to each hypothesis as follows:(18)
Δ
1
=
Δ
1
KS
*
⊕
Δ
1
CM
*
=
Δ
1
KS
*
Δ
1
CM
*
1
-
(
Δ
1
KS
*
Δ
0
CM
*
+
Δ
0
KS
*
Δ
1
CM
*
)
,
Δ
0
=
Δ
0
KS
*
⊕
Δ
0
CM
*
=
Δ
0
KS
*
Δ
0
CM
*
1
-
(
Δ
1
KS
*
Δ
0
CM
*
+
Δ
0
KS
*
Δ
1
CM
*
)
.Finally, the global decision will be made as follows:(19)
G
=
H
1
,
if
Δ
1
Δ
0
≥
η
,
G
=
H
0
,
otherwise
,
where η is the threshold for global decision.According to the global decision,r
1
,
j
t or r
0
,
j
t will be updated for the next sensing interval as follows, respectively.(i)
If the global decision isG
(
i
)
=
0, we update r
1
,
j
t
(
i
) by using (14).
(ii)
Otherwise, we updater
0
,
j
t
(
i
) by using (15).


## 3.1. Energy Detection

At the sensing interval for thejth CU, the local spectrum sensing is to decide between the two following hypotheses:
(9)
H
0
:
s
j
(
k
)
=
n
j
(
k
)
,
H
1
:
s
j
(
k
)
=
h
j
p
(
k
)
+
n
j
(
k
)
,
where H
0 and H
1 correspond to the hypothesis of the absence and presence of the PU signal, respectively, h
j denotes the amplitude gain of the channel, s
(
k
) is the signal transmitted from the PU, n
j
(
k
) is the additive white Gaussian noise, and k is index of sensing sample at each sensing interval.A received signal energy of a sensing sample,e
k
,
j, is given as
(10)
e
k
,
j
=
{
|
n
j
(
k
)
|
2
,
H
0
|
h
j
s
(
k
)
+
n
j
(
k
)
|
2
,
H
1
.


## 3.2. BPA Estimation

The GOF test statistics of the current sensing datae
k
,
j
(
k
=
1
,
…
,
M
) of the j
th CU will be calculated according to each hypothesis of PU signal based on (3) and (5) as follows, respectively:
(11)
D
h
,
j
KS
=
sup
⁡
{
|
F
(
e
k
,
j
)
-
F
h
(
e
k
,
j
)
|
:
k
=
1,2
,
…
,
M
}
,
D
h
,
j
CM
=
1
12
M
+
∑
i
=
1
M
‍
[
2
i
-
1
2
M
-
F
h
(
e
k
,
j
)
]
2
,
where h
=
{
0,1
} is index of hypothesis H
h of PU signal, e
k
,
j is the received signal energy of kth sensing sample of the jth CU, F
(
·
) and F
h
(
·
) are empirical CDF of observed sensing sample and CDF of H
h hypothesis of PU, and M is the number of samples for each sensing interval.It is noteworthy that normal CU and malicious CU have different characteristics ofD
1
,
j
t and D
0
,
j
t as shown in Table 1, where t indexes types of GOF tests: KS and CM.Table 1
Reputation ranges according to each type of CUs.
Status of PU
Normal CU
A
l
w
a
y
s
Y
e
s CU
A
l
w
a
y
s
N
o CU
O
p
p
o
s
i
t
e CU
H
1
D
0
,
j
t
≫
D
1
,
j
t
D
0
,
j
t
≫
D
1
,
j
t
D
0
,
j
t
≪
D
1
,
j
t
D
0
,
j
t
≪
D
1
,
j
t
D
1
,
j
t
≈
0
D
1
,
j
t
≈
0
D
0
,
j
t
≈
0
D
0
,
j
t
≈
0
r
0
,
j
t
(
i
)
≫
0
r
0
,
j
t
(
i
)
≫
0
r
0
,
j
t
(
i
)
≪
0
r
0
,
j
t
(
i
)
≪
0
H
0
D
0
,
j
t
≪
D
1
,
j
t
D
0
,
j
t
≫
D
1
,
j
t
D
0
,
j
t
≪
D
1
,
j
t
D
0
,
j
t
≫
D
1
,
j
t
D
0
,
j
t
≈
0
D
1
,
j
t
≈
0
D
0
,
j
t
≈
0
D
1
,
j
t
≈
0
r
1
,
j
t
(
i
)
≫
0
r
1
,
j
t
(
i
)
≪
0
r
1
,
j
t
(
i
)
≫
0
r
1
,
j
t
(
i
)
≪
0Based on the values ofD
1
,
j
t and D
0
,
j
t, we will estimate BPA of current sensing data of each CU and their “reputation level” for robust CSS as follows:
(12)
Δ
1
,
j
t
=
D
0
,
j
t
D
1
,
j
t
+
D
0
,
j
t
R
0
,
j
t
,
Δ
0
,
j
t
=
D
1
,
j
t
D
1
,
j
t
+
D
0
,
j
t
R
1
,
j
t
,
where R
h
,
j
t is “reputation level” of the jth CU according to hypothesis h, and it can be determined based on history observation of the jth CU as follows:
(13)
R
h
,
j
t
(
i
)
=
r
h
,
j
t
(
i
-
1
)
∑
j
r
h
,
j
t
(
i
-
1
)
,
where i is the index of current sensing interval and r
1
,
j
t
(
i
-
1
) and r
0
,
j
t
(
i
-
1
) are updated from the previous sensing interval according to global decision:
(14)
r
1
,
j
t
(
i
-
1
)
=
r
1
,
j
t
(
i
-
2
)
+
(
D
1
,
j
t
(
i
-
1
)
-
D
0
,
j
t
(
i
-
1
)
)
,
(15)
r
0
,
j
t
(
i
-
1
)
=
r
0
,
j
t
(
i
-
2
)
+
(
D
0
,
j
t
(
i
-
1
)
-
D
1
,
j
t
(
i
-
1
)
)
.By usingr
1
,
j
t and r
0
,
j
t, types of CUs will be easily distinguished. The normal CU has positive value of both r
1
,
j
t and r
0
,
j
t that will be increased after updating step. r
1
,
j
t ofalways Yes and r
0
,
j
t ofalways No malicious CUs are almost negative and tend to decrease after updating step. On the other hand, both values of opposite CU are negative and have a tendency to decrease. We define “malicious threshold” as ρ to reject the attack of malicious CR in CSS, so that the CU, which has either r
1
,
j
t
<
0 or r
0
,
j
t
<
0, will be determined as malicious CU. The sensing data of malicious CUs will not be considered to make global decision by giving them r
1
,
j
t
=
0 and r
0
,
j
t
=
0.The BPA of all CUs will be combined with their reputation levels as(16)
Δ
1
t
=
1
n
Ω
∑
j
∈
Ω
‍
D
0
,
j
t
D
1
,
j
t
+
D
0
,
j
t
R
0
,
j
t
,
Δ
0
t
=
1
n
Ω
∑
j
∈
Ω
‍
D
1
,
j
t
D
1
,
j
t
+
D
0
,
j
t
R
1
,
j
t
,
where Ω and n
Ω are set of normal CUs and number of members of the set, respectively.Because the error in estimatingΔ
0
t and Δ
1
t,  Δ
0
t
+
Δ
1
t can be bigger than 1, we need to normalize those values as
(17)
Δ
0
t
*
=
Δ
0
t
Δ
0
t
+
Δ
1
t
,
Δ
1
t
*
=
Δ
1
t
Δ
0
t
+
Δ
1
t
.


## 3.3. DS Theory Combination

The DS theory will be used to combine the BPA of both GOF tests according to each hypothesis as follows:(18)
Δ
1
=
Δ
1
KS
*
⊕
Δ
1
CM
*
=
Δ
1
KS
*
Δ
1
CM
*
1
-
(
Δ
1
KS
*
Δ
0
CM
*
+
Δ
0
KS
*
Δ
1
CM
*
)
,
Δ
0
=
Δ
0
KS
*
⊕
Δ
0
CM
*
=
Δ
0
KS
*
Δ
0
CM
*
1
-
(
Δ
1
KS
*
Δ
0
CM
*
+
Δ
0
KS
*
Δ
1
CM
*
)
.Finally, the global decision will be made as follows:(19)
G
=
H
1
,
if
Δ
1
Δ
0
≥
η
,
G
=
H
0
,
otherwise
,
where η is the threshold for global decision.According to the global decision,r
1
,
j
t or r
0
,
j
t will be updated for the next sensing interval as follows, respectively.(i)
If the global decision isG
(
i
)
=
0, we update r
1
,
j
t
(
i
) by using (14).
(ii)
Otherwise, we updater
0
,
j
t
(
i
) by using (15).


## 4. Simulation Results

In this section, simulation results of the proposed scheme and other soft combination schemes such as maximum gain combination (MGC) and equal gain combination (EGC) are provided. The network is considered in which5 CUs exist and some of them can be malicious CUs.In order to verify the reliability of the proposed combination scheme, we perform a simulation without considering malicious CU. The sensing results in Figure2 show that the proposed scheme can obtain better sensing performance in comparison with EGC scheme and obtain a similar sensing performance to that of the MGC scheme when no malicious CU is considered.Figure 2
ROC of the proposed scheme and reference schemes when no malicious CU is considered.The robustness of the proposed scheme will be investigated in the network with the appearance ofalways No,always Yes, andopposite malicious CUs in the network. Figures 3 and 4 show performance of the proposed scheme when 4 CUs arealways No oralways Yes malicious CUs among 5 CUs in the network. The results show that the proposed scheme with all CUs can achieve much better sensing performance than that one of the MGC and EGC schemes. This means that, by applying GOF test to CSS, the proposed scheme can detect the presence of those types of malicious CUs and reject their harmful effects to sensing process.Figure 3
ROC of the proposed scheme and reference schemes when4
always No malicious CUs are considered.Figure 4
ROC of the proposed scheme and reference schemes when4
always Yes malicious CUs are considered.Opposite malicious CU causes the most damage to sensing performance. However, the proposed scheme is expected to protect CSS against this type of malicious CU. Figure5 shows the sensing performance of the network when 4 CUs areopposite malicious among 5 CUs. MGC and EGC with all CUs provide very low performance due to the attack ofopposite malicious CU. However, the proposed scheme can defend their attacks and achieve high sensing performance.Figure 5
ROC of the proposed scheme and reference schemes when4
opposite malicious CUs are considered.


## 5. Conclusion

In this paper, multi-GOF tests are proposed to measure the difference between sensing data of normal CU and that of malicious CU. Further, the DS theory is used to combine results of multi-GOF tests. The proposed scheme considers the appearance of the most common types of malicious CU:always Yes,always No, andopposite types. The simulation results prove that the proposed scheme can reject almost harmful effect from those malicious CUs to protect CSS.



---

*Source: 752507-2014-05-18.xml*