Outage Analysis of OFDM-Based Cognitive AF Relay Network in the Presence of Narrowband Interference

ETRI Journal.
2015.
May,
37(3):
460-470

- Received : March 10, 2014
- Accepted : February 03, 2015
- Published : May 01, 2015

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

Orthogonal frequency-division multiplexing (OFDM) is one of the most widely used technologies in current wireless communication systems and standards. Cognitive radio (CR) provides a robust solution to the problem of spectrum congestion as it offers opportunistic usage of frequency bands that are not occupied by primary users. Due to the underlying sensing, spectrum shaping, scaling, and interoperable capabilities of OFDM, it has been adapted as a best transmission technology for CR wireless systems. However, the performance of an OFDM-based CR wireless system is affected by the existence of narrowband interference (NBI) from other users. Further, due to carrier frequency offset in NBI sources, NBI energy may spread over all subcarriers of an OFDM signal. In this paper, a fixed Amplify-and-Forward (AF) relay that operates at a frequency band that is different from that of direct mode is introduced to suppress the effect of NBI. Analytical expressions are derived for outage probability in direct, AF-relay, and incremental relaying modes. The outage performance of the proposed AF relay–based CR network is proven to be better than that of direct mode.
Notations
. Time-domain and frequency-domain vectors and matrices are denoted by
a
^{TD}
,
a
^{FD}
and
A
^{TD}
,
A
^{FD}
, respectively, and scalars are indicated by lowercase letters. An
N
×
N
discrete Fourier transform (DFT) matrix and an
N
×
N
identity matrix are denoted by
F
_{N}
and
I
_{N}
, respectively. An all-zero matrix of size
m
×
n
is denoted by 0
_{m×n}
. The matrix complex conjugate transpose, complex conjugate, transpose, and inverse operations are denoted by ( )
^{H}
, ( )
^{*}
, ( )
^{T}
, and ( )
^{−1}
, respectively.
N
subcarriers and uses
v
zero-padded guard subcarriers to suppress intersymbol interference
[18]
.
AF relay–based CR network.
L
_{1}
× 1 channel impulse response (CIR) vector between SU1 and SU3 is defined as
h
_{13}
= [
h
_{13}
(0),
h
_{13}
(1), ... ,
h
_{13}
(
L
_{1}
−1)]
^{T}
. The
L
_{2}
× 1 CIR vector of NBI source two at SU3 is defined as
h
_{b3}
= [
h
_{b3}
(0),
h
_{b3}
(1), ... ,
h
_{b3}
(
L
_{2}
−1)]
^{T}
.
The
P
× 1 time-domain received vector at SU3 is given by
(1) $${\text{y}}_{13}^{\text{TD}}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}{\text{H}}_{13}^{\text{TD}}{\text{F}}_{\text{zp}}{\text{x}}_{1}^{\text{FD}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\text{A}}_{\text{fo},3}^{(1),\text{TD}}{\text{H}}_{b3}^{\text{TD}}{\text{b}}_{3}^{\text{TD}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\text{z}}_{3}^{\text{TD}}\text{\hspace{0.17em}},$$
where
P
=
N
+
v
;
P
×
P
circulant matrices of an CIR vectors
h
_{13}
and
h
_{b3}
, respectively;
F
_{zp}
is a
P
×
N
zero-padded inverse DFT matrix at SU1 (it is defined as
N
× 1 data vector at
N
subcarriers of an OFDM symbol; and
P
× 1 NBI signal vector at SU3. The interference power per subcarrier (IPPS) of an NBI signal is denoted by
P
× 1 zero-mean complex white Gaussian noise vector with elements that are independent and identically distributed with variance
P
×
P
diagonal matrix
β
is a uniformly distributed random variable over the interval [−1/2, 1/2]
[19]
. By taking a DFT, the frequency-domain representation of the received vector
(2) $${\text{y}}_{13}^{\text{FD}}\text{\hspace{0.17em}}={\text{F}}_{p}{\text{y}}_{13}^{\text{TD}}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}{\text{F}}_{p}\left({\text{H}}_{13}^{\text{TD}}{\text{F}}_{\text{zp}}{\text{x}}_{1}^{\text{FD}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\text{A}}_{\text{fo},3}^{(1),\text{TD}}{\text{H}}_{b3}^{\text{TD}}{\text{b}}_{3}^{\text{TD}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\text{z}}_{3}^{\text{TD}}\right).$$
Since
λ
_{13}
= [
λ
_{13}
(0),
λ
_{13}
(1), ... ,
λ
_{13}
(
P
−1)]
^{T}
and
λ
_{b3}
= [
λ
_{b3}
(0),
λ
_{b3}
(1), ... ,
λ
_{b3}
(
P
−1)]
^{T}
be
P
-point DFTs of
h
_{13}
and
h
_{b3}
, respectively. Then, the matrices
λ
_{13}
(0),
λ
_{13}
(1), ... ,
λ
_{13}
(
P
−1)} and
λ
_{b3}
(0),
λ
_{b3}
(1), ... ,
λ
_{b3}
(
P
−1)}. Now, (2) can be simply written as
(3) $${\text{y}}_{13}^{\text{FD}}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}{\text{H}}_{13}^{\text{FD}}\text{V}{\text{x}}_{1}^{\text{FD}}+\text{\hspace{0.17em}}{\text{A}}_{3}^{(1),\text{FD}}{\text{H}}_{b3}^{\text{FD}}{\text{b}}_{3}^{\text{FD}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\text{z}}_{3}^{\text{FD}},$$
where
V
=
F
_{P}
F
_{zp}
is the precoding matrix;
P
× 1 NBI signal vector the vector
k
th subcarrier element
b
_{3}
(
k
) ≠ 0 if
k
∈
I
_{b3}
, where
I
_{b3}
is a set that contains an index of the non-zero elements which is NBI affected subcarriers.
L
_{3}
× 1 CIR vector between SU1 and SU2 is defined as
h
_{12}
= [
h
_{12}
(0),
h
_{12}
(1), ... ,
h
_{12}
(
L
_{3}
−1)]
^{T}
. The
L
_{4}
× 1 CIR vector between NBI source one and SU2 is defined as
h
_{b2}
= [
h
_{b2}
(0),
h
_{b2}
(1), ... ,
h
_{b2}
(
L
_{4}
−1)]
^{T}
. The
P
× 1 time-domain received vector
(4) $${\text{y}}_{12}^{\text{TD}}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}{\text{H}}_{12}^{\text{TD}}{\text{F}}_{\text{zp}}{\text{x}}_{1}^{\text{FD}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\text{A}}_{\text{fo},2}^{(2),\text{TD}}{\text{H}}_{b2}^{\text{TD}}{\text{b}}_{2}^{\text{TD}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\text{z}}_{2}^{\text{TD}},$$
where
P
×
P
circulant matrices of the CIR vectors
h
_{12}
and
h
_{b2}
, respectively. In addition,
(5) $${\text{y}}_{12}^{\text{FD}}={\text{F}}_{p}{\text{y}}_{12}^{\text{TD}}=\text{\hspace{0.17em}}{\text{F}}_{p}\left({\text{H}}_{12}^{\text{TD}}{\text{F}}_{\text{zp}}{\text{x}}_{1}^{\text{FD}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\text{A}}_{\text{fo},2}^{(2),\text{TD}}{\text{H}}_{b2}^{\text{TD}}{\text{b}}_{2}^{\text{TD}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\text{z}}_{2}^{\text{TD}}\right).$$
Since
(6) $${\text{y}}_{12}^{\text{FD}}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}{\text{H}}_{12}^{\text{FD}}\text{V}{\text{x}}_{1}^{\text{FD}}+\text{\hspace{0.17em}}{\text{A}}_{2}^{(2),\text{FD}}{\text{H}}_{b2}^{\text{FD}}{\text{b}}_{2}^{\text{FD}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}{\text{z}}_{2}^{\text{FD}},$$
where
λ
_{12}
and
λ
_{b2}
, respectively, and
L
_{5}
× 1 CIR vector between SU2 and SU3 is defined as
h
_{23}
= [
h
_{23}
(0),
h
_{23}
(1), ... ,
h
_{23}
(
L
_{5}
−1)]
^{T}
The
L
_{6}
× 1 CIR vector of NBI source two at SU3 is defined as
h
_{b3}
= [
h
_{b3}
(0),
h
_{b3}
(1), ... ,
h
_{b3}
(
L
_{6}
−1)]. The time-domain received vector
(7) $${\text{y}}_{23}^{\text{TD}}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}\rho {\text{H}}_{23}^{\text{TD}}\left({\text{y}}_{12}^{\text{TD}}\right)+{\text{A}}_{\text{fo},3}^{(3),\text{TD}}{\text{H}}_{b3}^{\text{TD}}{\text{b}}_{3}^{\text{TD}}+{\text{n}}_{3}^{\text{TD}},$$
where
h
_{23}
and
h
_{b3}
, respectively. In addition,
(8) $$\begin{array}{l}{\text{y}}_{23}^{\text{TD}}\text{\hspace{0.17em}\hspace{0.17em}}=\text{\hspace{0.17em}}\rho {\text{H}}_{23}^{\text{TD}}{\text{H}}_{12}^{\text{TD}}{\text{F}}_{\text{zp}}{\text{x}}_{1}^{\text{FD}}+\text{\hspace{0.17em}}\rho {\text{H}}_{23}^{\text{TD}}{\text{A}}_{\text{fo},2}^{(2),\text{TD}}{\text{H}}_{b2}^{\text{TD}}{\text{b}}_{2}^{\text{TD}}\\ \text{\hspace{1em}\hspace{1em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}}+{\text{A}}_{\text{fo},3}^{(3),\text{TD}}{\text{H}}_{b3}^{\text{TD}}{\text{b}}_{3}^{\text{TD}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}\rho {\text{H}}_{23}^{\text{TD}}{\text{z}}_{2}^{\text{TD}}+{\text{n}}_{3}^{\text{TD}}.\end{array}$$
By taking a DFT, the frequency-domain representation of the received vector
(9) $$\begin{array}{l}{\text{y}}_{23}^{\text{FD}}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}{\text{F}}_{p}{\text{y}}_{23}^{\text{TD}}\text{\hspace{0.17em}}=\text{\hspace{0.17em}}{\text{F}}_{p}\{\rho {\text{H}}_{23}^{\text{TD}}{\text{H}}_{12}^{\text{TD}}{\text{F}}_{\text{zp}}{\text{x}}_{1}^{\text{FD}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}\rho {\text{H}}_{23}^{\text{TD}}{\text{A}}_{\text{fo},2}^{(2),\text{TD}}{\text{H}}_{b2}^{\text{TD}}{\text{b}}_{2}^{\text{TD}}\\ \text{\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}}+{\text{A}}_{\text{fo},3}^{(3),\text{TD}}{\text{H}}_{b3}^{\text{TD}}{\text{b}}_{3}^{\text{TD}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}\rho {\text{H}}_{23}^{\text{TD}}{\text{z}}_{2}^{\text{TD}}+{\text{n}}_{3}^{\text{TD}}\},\end{array}$$
where
(10) $$\begin{array}{l}{\text{y}}_{23}^{\text{FD}}=\rho {\text{H}}_{23}^{\text{FD}}{\text{H}}_{12}^{\text{FD}}\text{V}{\text{x}}_{1}^{\text{FD}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}\rho {\text{H}}_{23}^{\text{FD}}{\text{A}}_{2}^{(2),\text{FD}}{\text{H}}_{b2}^{\text{FD}}{\text{b}}_{2}^{\text{FD}}\\ \text{\hspace{1em}\hspace{1em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}}+{\text{A}}_{3}^{(3),\text{FD}}{\text{H}}_{b3}^{\text{FD}}{\text{b}}_{3}^{\text{FD}}+\rho {\text{H}}_{23}^{\text{FD}}{\text{z}}_{2}^{\text{FD}}+{\text{n}}_{3}^{\text{FD}}.\end{array}$$
Let
λ
_{23}
= [
λ
_{23}
(0),
λ
_{23}
(1), ... ,
λ
_{23}
(
P
−1)]
^{T}
and
λ
_{b3}
= [
λ
_{b3}
(0),
λ
_{b3}
(1), ... ,
λ
_{b3}
(
P
−1)]
^{T}
be
P
-point DFTs of
h
_{23}
and
h
_{b3}
, respectively. The elements of the diagonal matrices
λ
_{23}
and
λ
_{b3}
, respectively;
k
th subcarrier is given by
(11) $$\begin{array}{l}{y}_{13}(k)={\lambda}_{13}(k)v(k,k){x}_{1}(k)+{\lambda}_{13}(k){\displaystyle \sum _{q=0,q\ne k}^{N-1}v}(k,q){x}_{1}(q)\\ \text{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}}+\text{}{\displaystyle \sum _{m=0,\text{\hspace{0.17em}}m\in {I}_{b3}}^{P-1}{\lambda}_{b3}}(m){a}_{\text{3}}^{\text{1}}\left(k,{(m-k)}_{P}\right){b}_{3}(m)+{z}_{3}(k)\\ \text{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}for}k=0,\text{\hspace{0.17em}}1,\text{\hspace{0.17em}}2,\text{\hspace{0.17em}\hspace{0.17em}}\mathrm{...}\text{\hspace{0.17em}\hspace{0.17em}},P-1.\end{array}$$
In this expression, the first term is the desired signal, where
x
_{1}
(
k
) is the
k
th subcarrier source data. The second and third terms denote the residual zero-padded interference and NBI signal, respectively. Since
k
th subcarrier.
Assuming that the CIR vectors
λ
_{23}
and
λ
_{b3}
are perfectly known, the instantaneous SINR of the direct-mode receive signal
k
th subcarrier is given by
(12) $${\text{\Gamma}}_{13}(k)=\frac{{\left|{\lambda}_{13}(k)\right|}^{2}{\left|v(k,k)\right|}^{2}\text{E}\left({\left|{x}_{1}(k)\right|}^{2}\right)}{\begin{array}{l}\{{\left|{\lambda}_{13}(k)\right|}^{2}\text{E}\left({\left|{\displaystyle \sum _{q=0,q\ne k}^{N-1}v(k,q){x}_{1}}(q)\right|}^{2}\right)\\ +\text{E}\left({\left|{\displaystyle \sum _{m=0,\text{\hspace{0.17em}}m\in {I}_{b3}}^{P-1}{\lambda}_{b3}}(m){a}_{\text{3}}^{\text{1}}\left(k,{(m-k)}_{P}\right){b}_{3}(m)\right|}^{2}\right)\\ +\text{E}\left(\left|{z}_{3}(k)\right|{}^{2}\right)\}.\end{array}}$$
Using the Cauchy–Schwartz inequality, the first term in the denominator of (12) can be expressed as
(13) $$\text{E}\left({\left|{\displaystyle \sum _{q=0,\text{\hspace{0.17em}\hspace{0.17em}}q\ne k}^{N-1}v(k,q){x}_{1}(q)}\right|}^{2}\right)\le \text{}{\displaystyle \sum _{q=0,\text{\hspace{0.17em}\hspace{0.17em}}q\ne k}^{N-1}{\left|v(k,q)\right|}^{2}}\text{E}\left({\displaystyle \sum _{q=0,\text{\hspace{0.17em}\hspace{0.17em}}q\ne k}^{N-1}{\left|{x}_{1}(q)\right|}^{2}}\right).$$
Since
v
(
k
,
k
)|
^{2}
; then (13) can be written as
(14) $$\text{E}\left({\left|{\displaystyle \sum _{q=0,q\ne k}^{N-1}v(k,q){x}_{1}(q)}\right|}^{2}\right)\le \left(1-{\left|v(k,k)\right|}^{2}\right){\text{P}}_{\text{s}}.$$
Similarly, the second term in the denominator of (12) can be expressed as
(15) $$\text{E}\left({\left|{\displaystyle \sum _{m=0,\text{\hspace{0.17em}}m\in {I}_{b3}}^{P-1}\text{}{\lambda}_{b3}}(m){a}_{\text{3}}^{\text{1}}\left(k,{(m-k)}_{P}\right){b}_{3}(m)\right|}^{2}\right)\le \text{}{\displaystyle \sum _{m=0,\text{\hspace{0.17em}}m\in {I}_{b3}}^{P-1}{\left|{\lambda}_{b3}(m)\right|}^{2}}{\varsigma}_{3}{\sigma}_{b3}^{2},$$
where
_{13}
(
k
) is simply written as
(16) $${\text{\Gamma}}_{13}(k)\ge \frac{{\left|{\lambda}_{13}(k)\right|}^{2}{\left|v(k,k)\right|}^{2}{\text{P}}_{\text{s}}}{{\left|{\lambda}_{13}(k)\right|}^{2}\left(1-{\left|v(k,k)\right|}^{2}\right){\text{P}}_{\text{s}}+\text{}{\displaystyle \sum _{m=0,\text{\hspace{0.17em}}m\in {I}_{b3}}^{P-1}{\left|{\lambda}_{b3}(m)\right|}^{2}}{\varsigma}_{3}{\sigma}_{b3}^{2}+{\sigma}_{z3}^{2}},$$
where
z
_{3}
(
k
)|
^{2}
) is the noise variance. Dividing the numerator and denominator by the source power P
_{s}
, we get
(17) $${\text{\Gamma}}_{13}(k)\text{\hspace{0.17em}}\ge \text{\hspace{0.17em}}\frac{{\left|{\lambda}_{13}(k)\right|}^{2}{\left|v(k,k)\right|}^{2}}{{\left|{\lambda}_{13}(k)\right|}^{2}\left(1-{\left|v(k,k)\right|}^{2}\right)\text{\hspace{0.17em}}+\text{}{\displaystyle \sum _{m=0,\text{\hspace{0.17em}}m\in {I}_{b3}}^{P-1}{\left|{\lambda}_{b3}(m)\right|}^{2}}{\varsigma}_{3}\frac{1}{\text{SIR1}}\text{\hspace{0.17em}}+\text{\hspace{0.17em}}\frac{1}{\text{SNR}}},$$
where SIR 1 = P
_{s}
/
_{s}
/
h
_{12}
and
h
_{23}
with amplification factor
ρ
at SU2. The SINR between SU1 and SU3 at the
k
th subcarrier in AF relay mode, Γ
_{123}
(
k
), is given by
[21]
as
(18) $${\text{\Gamma}}_{123}(k)=\left[{\text{\Gamma}}_{12}(k){\text{\Gamma}}_{23}(k)\right]/\left[{\text{\Gamma}}_{12}(k)+{\text{\Gamma}}_{23}(k)+1\right],$$
where Γ
_{IJ}
(
k
) is the SINR of the link joining SUI and SUJ,
I
,
J
∈ (1, 2) . The SINR of the link joining SU1 and SU2 is determined as follows. From (6), the received signal
k
th subcarrier is given by
(19) $$\begin{array}{ll}{y}_{12}\left(k\right)=\hfill & {\lambda}_{12}\left(k\right)v\left(k,k\right){x}_{1}\left(k\right)+{\lambda}_{12}\left(k\right){\displaystyle \sum _{q=0,q\ne k}^{N-1}v\left(k,q\right)}{x}_{1}\left(q\right)\hfill \\ \hfill & +{\displaystyle \sum _{m=0,m\in {I}_{b2}}^{P-1}{\lambda}_{b2}\left(m\right)}{a}_{2}^{2}\left(k,{\left(m-k\right)}_{p}\right){b}_{2}\left(m\right)+{z}_{2}\left(k\right),\hfill \\ \hfill & \text{}k=0,1,\mathrm{...}\text{,}P-1.\hfill \end{array}$$
In this expression, the first term is the desired signal. The second and third terms denote the residual zero-padded interference and the NBI signal, respectively. The fourth term is the noise. The instantaneous SINR of the received signal
y
_{12}
(
k
) is given by
(20) $$\begin{array}{cc}{\text{\Gamma}}_{12}(k)\hfill & =\frac{\text{\hspace{0.17em}}{\left|{\lambda}_{12}(k)\right|}^{2}{\left|v(k,k)\right|}^{2}\text{E}\left(\left|{x}_{1}(k)\right|{}^{2}\right)}{\text{\hspace{0.17em}\hspace{0.17em}}\{{\left|{\lambda}_{12}(k)\right|}^{2}\text{E}\left(\left|{\displaystyle \sum _{q=0,q\ne k}^{N-1}v(k,q){x}_{1}}(q)\right|{}^{2}\right)}\hfill \\ \hfill & \text{}+\text{E}\left(\left|{\displaystyle \sum _{m=0,m\in {I}_{b2}}^{P-1}{\lambda}_{b2}}(m){a}_{\text{2}}^{\text{2}}\left(k,{(m-k)}_{P}\right){b}_{2}(m)\right|{}^{2}\right)\hfill \\ \hfill & \text{}+\text{E}\left(\left|{z}_{2}(k)\right|{}^{2}\right)\}.\hfill \end{array}$$
The first and second terms in the denominator of (20) can be simplified as
(21) $$\text{E}\left(\left|{\displaystyle \sum _{q=0,q\ne k}^{N-1}v(k,q){x}_{1}(q)}\right|{}^{2}\right)\le \left(1-{\left|v(k,k)\right|}^{2}\right){\text{P}}_{\text{s}}$$
and
(22) $$\text{E}\left(\left|\text{}{\displaystyle \sum _{m=0,\text{\hspace{0.17em}}m\in {I}_{b2}}^{P-1}\text{}}{\lambda}_{b2}(m){a}_{\text{2}}^{\text{2}}\left(k,{(m\text{\hspace{0.17em}}-\text{\hspace{0.17em}}k)}_{P}\right){b}_{2}(m)\right|{}^{2}\right)\le \text{}{\displaystyle \sum _{m=0,\text{\hspace{0.17em}}m\in {I}_{b2}}^{P-1}\text{}}{\left|{\lambda}_{b2}(m)\right|}^{2}\text{\hspace{0.05em}}{\varsigma}_{2}{\sigma}_{b2}^{2},$$
respectively, where
(23) $${\text{\Gamma}}_{12}(k)\text{\hspace{0.17em}}\ge \text{\hspace{0.17em}}\frac{{\left|{\lambda}_{12}(k)\right|}^{2}{\left|v(k,k)\right|}^{2}}{{\left|{\lambda}_{12}(k)\right|}^{2}\left(1\text{\hspace{0.17em}}-{\left|v(k,k)\right|}^{2}\right)\text{\hspace{0.05em}}+\text{}{\displaystyle \sum _{m=0,\text{\hspace{0.17em}}m\in {I}_{b2}}^{P-1}{\left|{\lambda}_{b2}(m)\right|}^{2}}{\varsigma}_{2}\frac{1}{\text{SIR2}}\text{\hspace{0.17em}}+\frac{1}{\text{SNR}}},$$
where SIR2 = P
_{s}
/
_{s}
/
(24) $$\begin{array}{ll}{\text{\Gamma}}_{23}(k)\hfill & =\frac{{\left|{\lambda}_{23}(k)\right|}^{2}{\left|v(k,k)\right|}^{2}\text{E}\left(\left|{x}_{1}(k)\right|{}^{2}\right)}{\{{\left|{\lambda}_{23}(k)\right|}^{2}\text{E}\left(\left|{\displaystyle \sum _{q=0,q\ne k}^{N-1}v(k,q){x}_{1}}(q)\right|{}^{2}\right)}\hfill \\ \hfill & \text{\hspace{0.17em}}+\text{E}\left(\left|{\displaystyle \sum _{m=0,m\in {I}_{b3}}^{P-1}{\lambda}_{b3}}(m){a}_{\text{3}}^{\text{3}}\left(k,{(m-k)}_{P}\right){b}_{3}(m)\right|{}^{2}\right)\hfill \\ \hfill & \text{\hspace{0.17em}}+\text{E}\left(\left|{n}_{3}(k)\right|{}^{2}\right)\}.\hfill \end{array}$$
The first and second terms in the denominator of (24) can be simplified as
(25) $$\text{E}\left(\left|{\displaystyle \sum _{q=0,q\ne k}^{N-1}v(k,q){x}_{1}(q)}\right|{}^{2}\right)\le \left(1-{\left|v(k,k)\right|}^{2}\right){\text{P}}_{\text{s}}$$
and
(26) $$\text{E}\left(\left|{\displaystyle \sum _{m=0,\text{\hspace{0.17em}}m\in {I}_{b3}}^{P-1}{\lambda}_{b3}}(m){a}_{\text{3}}^{\text{3}}\left(k,{(m-k)}_{P}\right){b}_{3}(m)\right|{}^{2}\right)\le \text{}{\displaystyle \sum _{m=0,\text{\hspace{0.17em}}m\in {I}_{b3}}^{P-1}{\left|{\lambda}_{b3}(m)\right|}^{2}{\varsigma}_{3}^{\u2033}}{\sigma}_{b3}^{2},$$
respectively, where
(27) $${\text{\Gamma}}_{23}(k)\ge \frac{{\left|{\lambda}_{23}(k,k)\right|}^{2}{\left|v(k,k)\right|}^{2}}{{\left|{\lambda}_{23}(k,k)\right|}^{2}\left(1-{\left|v(k,k)\right|}^{2}\right)+\text{}{\displaystyle \sum _{m=0,\text{\hspace{0.17em}}m\in {I}_{b3}}^{P-1}{\left|{\lambda}_{b3}(m)\right|}^{2}}{\varsigma}_{3}^{\u2033}\frac{1}{\text{SIR3}}+\frac{1}{\text{SNR}}},$$
where
_{12}
(
k
) and Γ
_{23}
(
k
). Hence, the SINR Γ
_{123}
(
k
) in AF mode is written as
(28) $${\text{\Gamma}}_{123}(k)=\frac{{\text{\Gamma}}_{12}(k){\text{\Gamma}}_{23}(k)}{{\text{\Gamma}}_{12}(k)+{\text{\Gamma}}_{23}(k)+1}\le \mathrm{min}\left({\text{\Gamma}}_{12}(k),{\text{\Gamma}}_{23}(k)\right).$$
R
b/s/Hz. In addition, outage of the incremental AF relay is analyzed.
_{13}
at the
k
th subcarrier for the given data rate of
R
b/s/Hz is defined as
(29) $${P}_{\text{out},k}^{\text{dir}}(R)\le P\left({\mathrm{log}}_{2}\left(1+{\text{\Gamma}}_{13}(k)\right)<R\right).$$
Let the threshold SINR
γ
_{1}
= 2
^{R}
− 1. Then, the upper bound of
(30) $${P}_{\text{out},k}^{\text{dir}}(R)\le P({\text{\Gamma}}_{13}(k)<{\gamma}_{1}).$$
The probability density function (PDF) of Γ
_{13}
(
k
) is determined as follows. Let (17) be defined as
(31) $${\text{\Gamma}}_{13}(k)=\frac{aX}{bY+cZ+d},$$
where
d
= 1/SNR. The random variables
X
and
Y
are exponentially distributed, and
Z
is chi-square distributed with 2
r
degrees of freedom
[22]
; the cardinality of the NBI signal vector is
r
. The PDF of the variables
X
,
Y
, and
Z
are given, respectively, by
(32) $${f}_{X}(x)=(1/2{\sigma}_{x}^{2})\mathrm{exp}(-x/2{\sigma}_{x}^{2}),x\ge 0,$$
(33) $${f}_{Y}(y)=(1/2{\sigma}_{y}^{2})\mathrm{exp}(-y/2{\sigma}_{y}^{2}),y\ge 0,$$
(34) $${f}_{Z}(z)=\frac{1}{{\sigma}_{z}^{{}^{2r}}{2}^{r}\text{\Gamma}(r)}{z}^{r-1}\mathrm{exp}\left(\raisebox{1ex}{$-z$}\!\left/ \!\raisebox{-1ex}{$2{\sigma}_{z}^{2}$}\right.\right),z\ge 0.$$
Since the desired signal channel and the interference channel are negatively quadrant dependent
[23]
, the outage probability of
X
for the given random variables
Y
and
Z
is given by
(35) $${P}_{\text{out},k}^{\text{dir}}(R)\text{\hspace{0.05em}}\le \text{\hspace{0.05em}}{\displaystyle \underset{0}{\overset{\infty}{\int}}{\displaystyle \underset{0}{\overset{\infty}{\int}}P\left(X\text{\hspace{0.05em}}<{\gamma}_{1}\frac{(by+cz+d)}{a}|{}_{Y=y,Z=z}\right)}}\text{\hspace{0.17em}}{f}_{Y}(y){f}_{Z}(z)\text{\hspace{0.17em}d}y\text{d}z.$$
Since |
v
(
k
,
k
)| ≈ 1 and
a
≈1, (35) can be simplified as
(36) $${P}_{\text{out},k}^{\text{dir}}(R)\le {\displaystyle \underset{0}{\overset{\infty}{\int}}{\displaystyle \underset{0}{\overset{\infty}{\int}}{\displaystyle \underset{0}{\overset{\left(by+cz+d\right){\gamma}_{1}}{\int}}{f}_{X}(x)}}}\text{\hspace{0.17em}}{f}_{Y}(y){f}_{Z}(z)\text{\hspace{0.17em}d}x\text{d}y\text{d}z.$$
Evaluating the inner integral by substituting (32), it becomes
(37) $${P}_{\text{out},k}^{\text{dir}}(R)\text{\hspace{0.17em}}\le \text{\hspace{0.17em}}{\displaystyle \underset{0}{\overset{\infty}{\int}}{\displaystyle \underset{0}{\overset{\infty}{\int}}\left[1-\mathrm{exp}\left(-{\gamma}_{1}(by\text{\hspace{0.17em}}+\text{\hspace{0.17em}}cz\text{\hspace{0.17em}}+\text{\hspace{0.17em}}d)\right)\right]\text{\hspace{0.17em}}}}{f}_{Y}(y){f}_{Z}(z)\text{\hspace{0.17em}d}y\text{d}z.$$
Now, (37) can be rewritten as
(38) $${P}_{\text{out},k}^{\text{dir}}(R)\text{\hspace{0.17em}}\le \text{\hspace{0.17em}}1-{\text{e}}^{-{\gamma}_{1}d}{\displaystyle \underset{0}{\overset{\infty}{\int}}{\text{e}}^{-{\gamma}_{1}by}}\text{\hspace{0.17em}}{f}_{Y}(y)\text{\hspace{0.17em}d}y{\displaystyle \underset{0}{\overset{\infty}{\int}}{\text{e}}^{-{\gamma}_{1}cz}}{f}_{Z}(z)\text{\hspace{0.17em}d}z.$$
Substituting (33) and (34) in (38), then integrating, the outage probability is determined as
(39) $${P}_{\text{out},k}^{\text{dir}}(R)\le 1-\mathrm{exp}(-{\gamma}_{1}d)\frac{{\left[1+({\gamma}_{1}c/r)\right]}^{-r}}{(1+{\gamma}_{1}b)}.$$
Substituting
b
,
c
,
d
, and
γ
_{1}
, the outage probability of the direct-mode SINR in terms of the data rate is given by
(40) $$\begin{array}{c}{P}_{\text{out},k}^{\text{dir}}(R)\le 1-\mathrm{exp}\left(-\left({2}^{R}-1\right)/\text{SNR}\right)\\ \text{}\times \frac{{\left[1+\left(\left({2}^{R}-1\right){\varsigma}_{3}/r\text{SIR1}\right)\right]}^{-r}}{\left[1+\left({2}^{R}-1\right)\left(1-{\left|v(k,k)\right|}^{2}\right)\right]}.\end{array}$$
_{AF}
, at the
k
th subcarrier is defined as
(41) $${P}_{\text{out},k}^{\text{AF}}(R)=P\left((1/2){\mathrm{log}}_{2}\left(1+{\text{\Gamma}}_{\text{AF}}(k)\right)<R\right).$$
Since maximal ratio combining is applied at the destination in the second time slot, the AF relay–mode SINR Γ
_{AF}
is defined as Γ
_{AF}
(
k
) = Γ
_{13}
(
k
) + Γ
_{123}
(
k
) . Now, (41) can be written as
(42) $${P}_{\text{out},k}^{\text{AF}}(R)=P\left({\text{\Gamma}}_{13}(k)+{\text{\Gamma}}_{123}(k)<{2}^{2R}-1\right).$$
Given the threshold SINR
γ
_{1}
for Γ
_{13}
(
k
) , the lower bound of the AF outage probability is given by
(43) $${P}_{\text{out},k}^{\text{AF}}(R)\ge {\displaystyle \underset{0}{\overset{{2}^{2R}-1}{\int}}P\left({\text{\Gamma}}_{123}(k)<{2}^{2R}-1-{\gamma}_{1}\right)\text{\hspace{0.17em}}}{f}_{{\text{\Gamma}}_{13}}({\gamma}_{1})\text{\hspace{0.17em}d}{\gamma}_{1}.$$
Let
γ
_{2}
= 2
^{2R}
−1−
γ
_{1}
. Then, (43) can be written as
(44) $${P}_{\text{out},k}^{\text{AF}}(R)\ge {\displaystyle \underset{0}{\overset{{2}^{2R}-1}{\int}}{\displaystyle \underset{0}{\overset{{2}^{2R}-1-{\gamma}_{1}}{\int}}{f}_{{\text{\Gamma}}_{123}}({\gamma}_{2}){f}_{{\text{\Gamma}}_{13}}({\gamma}_{1})\text{\hspace{0.17em}d}{\gamma}_{1}}}\text{d}{\gamma}_{2}.$$
where
f
_{Γ123}
(
γ
_{2}
) and
f
_{Γ13}
(
γ
_{1}
) are PDF functions of Γ
_{123}
(
k
) and Γ
_{13}
(
k
), respectively.
Since Γ
_{123}
(
k
) is the minimum of Γ
_{12}
(
k
) and Γ
_{23}
(
k
), the cumulative distribution function (CDF) of Γ
_{123}
(
k
) can be written as
[21]
(45) $${F}_{{\text{\Gamma}}_{123}}({\gamma}_{2})=1-\left[\left(1-{F}_{{\text{\Gamma}}_{12}}({\gamma}_{2})\right)\left(1-{F}_{{\text{\Gamma}}_{23}}({\gamma}_{2})\right)\right].$$
Using (39), the CDF of Γ
_{12}
(
k
) is computed as
(46) $${F}_{{\text{\Gamma}}_{12}}({\gamma}_{2})=1-\mathrm{exp}\left(-\frac{{\gamma}_{2}}{\text{SNR}}\right)\frac{{\left[1+\left({\gamma}_{2}{\varsigma}_{2}/r\text{SIR2}\right)\right]}^{-r}}{(1+{\gamma}_{2}b)}.$$
By using the approximations (1 + (
a
/
r
))
^{−r}
≈ exp(−
a
) and |
v
(
k
,
k
)|
^{2}
≈ 1, (46) can be simplified as
(47) $${F}_{{\text{\Gamma}}_{12}}({\gamma}_{2})=1-\mathrm{exp}\left(-\left(1/\text{SNR}\right)-\left({\varsigma}_{2}/\text{SIR2}\right)\right){\gamma}_{2}.$$
Similarly, the CDF of Γ
_{23}
(
k
) can be expressed as
(48) $${F}_{{\text{\Gamma}}_{23}}({\gamma}_{2})=1-\text{exp}\left(-\left(1/\text{SNR}\right)-\left({\varsigma}_{3}^{\u2033}/\text{SIR}3\right)\right){\gamma}_{2}.$$
Substituting (47) and (48) in (45), the CDF of of Γ
_{123}
(
k
) is obtained as
(49) $${F}_{{\text{\Gamma}}_{123}}\left({\gamma}_{2}\right)=1-\text{exp}\left(-\left(\frac{2}{\text{SNR}}+\frac{{\varsigma}_{2}}{\text{SIR}2}+\frac{{\varsigma}_{3}^{\u2033}}{\text{SIR}3}\right)\right){\gamma}_{2}.$$
Let
(50) $${F}_{{\text{\Gamma}}_{123}}({\gamma}_{2})=1-\mathrm{exp}(-{g}_{2}{\gamma}_{2}).$$
Differentiating (50) with respect to
γ
_{2}
, the PDF of Γ
_{123}
(
k
) is determined as
(51) $${f}_{{\text{\Gamma}}_{123}}({\gamma}_{2})=\text{d}\left({F}_{{\text{\Gamma}}_{123}}({\gamma}_{2})\right)/\text{d}{\gamma}_{2}={g}_{2}\mathrm{exp}(-{g}_{2}{\gamma}_{2}).$$
Using (39), the CDF of Γ
_{13}
(
k
) is given by
(52) $${F}_{{\text{\Gamma}}_{13}}({\gamma}_{1})=1-\left[\mathrm{exp}\left(-\frac{{\gamma}_{1}}{\text{SNR}}\right)\frac{{\left[1+\left({\gamma}_{1}{\varsigma}_{3}/\text{SIR1}r\right)\right]}^{-r}}{\left[1+{\gamma}_{1}\left(1-{\left|v(k,k)\right|}^{2}\right)\right]}\right].$$
It can be simplified as
(53) $${F}_{{\text{\Gamma}}_{13}}({\gamma}_{1})=1-\mathrm{exp}(-{g}_{1}{\gamma}_{1}),$$
where
g
_{1}
= (1/SNR +
ς
_{3}
/SIR1). Then, the PDF of Γ
_{13}
(
k
) is determined as
(54) $${f}_{{\text{\Gamma}}_{13}}({\gamma}_{1})=\text{d}\left({F}_{{\text{\Gamma}}_{13}}({\gamma}_{1})\right){\displaystyle /}\text{d}{\gamma}_{1}={g}_{1}\mathrm{exp}(-{g}_{1}{\gamma}_{1}).$$
Substituting the PDF of Γ
_{13}
(
k
) in (44) and evaluating the inner integral, we get
(55) $${P}_{\text{out},k}^{\text{AF}}(R)\ge {\displaystyle \underset{0}{\overset{{2}^{2R}-1}{\int}}\left[1-\mathrm{exp}\left(-{g}_{2}\left({2}^{2R}-1-{\gamma}_{1}\right)\right)\right]}\text{\hspace{0.17em}}{f}_{{\text{\Gamma}}_{13}}({\gamma}_{1})\text{\hspace{0.17em}d}{\gamma}_{1}.$$
Substituting the PDF of Γ
_{13}
(
k
) in (55), the lower bound of the outage probability for AF relay mode is determined as
(56) $$\begin{array}{c}{P}_{\text{out},k}^{\text{AF}}(R)\ge \left[1-\mathrm{exp}\left(-{g}_{1}\left({2}^{2R}\text{\hspace{0.05em}}-1\right)\right)\right]-[{\left({g}_{1}/\left({g}_{1}-{g}_{2}\right)\right)}_{{}_{}^{}}^{{}^{}}\\ \text{}\times \mathrm{exp}\left(-{g}_{2}\left({2}^{2R}\text{\hspace{0.05em}}-1\right)\right)]\left[1-\mathrm{exp}\left(-\left({g}_{1}-{g}_{2}\right)\left({2}^{2R}\text{\hspace{0.05em}}-1\right)\right)\right].\end{array}$$
R
b/s/Hz through a specified channel for the duration of the signal transmission. Indirectly, it reduces the number of degrees of freedom of the channel. But, in incremental AF relaying protocols, data is transmitted at either
R
/2 b/s/Hz or
R
b/s/Hz. The AF relay transmits data at the rate of
R
/2 b/s/Hz only if the direct-mode transmission fails to send the data at the data rate of
R
b/s/Hz, thereby an incremental AF relay network utilizes the relay channel effectively and improves the system performance
[17]
. The outage probability of the incremental AF relay mode is given by
(57) $$\begin{array}{ll}{P}_{\text{out},k}^{\text{IAF}}\left(R\right)\hfill & =P\left({\text{\Gamma}}_{13}\left(k\right)<R\right)\hfill \\ \hfill & \text{}\times P\left(\left(1/2\right){\mathrm{log}}_{2}\left(1+{\text{\Gamma}}_{\text{AF}}\left(k\right)\right)\left(R/2\right){|}_{{\text{\Gamma}}_{13}\left(k\right)R}\right).\hfill \end{array}$$
In an incremental AF relaying protocol, the intersection of the direct mode and AF relay mode outage probabilities is exactly the AF relay mode outage probability at half the rate of the data rate. Therefore, it can be written as
(58) $${P}_{\text{out},k}^{\text{IAF}}(R)=P\left((1/2){\mathrm{log}}_{2}\left(1+{\text{\Gamma}}_{\text{AF}}(k)\right)<(R/2)\right).$$
Using (56), the lower bound of the outage probability of the incremental AF relay mode is given by
(59) $$\begin{array}{c}{P}_{\text{out},k}^{\text{IAF}}(R)\ge \left[1-\mathrm{exp}\left(-{g}_{1}\left({2}^{R}-1\right)\right)\right]-\left({g}_{1}/({g}_{1}-{g}_{2})\right)\\ \text{}\times \mathrm{exp}\left(-{g}_{2}\left({2}^{R}-1\right)\right)\left[1-\mathrm{exp}\left(-({g}_{1}-{g}_{2})\left({2}^{R}-1\right)\right)\right].\end{array}$$

The outage performance of the proposed cognitive AF relay network in direct mode is shown in
Fig. 2
for the data rate of 1 b/s/Hz. The SIR is varied from 10 dB to 25 dB. It is observed that the analytical results are very close to the simulation results. As the noise dominates NBI at low SNR, the effect of NBI is more significant at high SNR. In direct mode, as the SIR increases from 20 dB to 25 dB, the outage probability decreases from 1 × 10
^{−2}
to 3 × 10
^{−3}
at the data rate of 1 b/s/Hz and SNR per subcarrier of 30 dB.
Outage performance of proposed cognitive AF relay network in direct mode at 1 b/s/Hz.
Figure 3
shows the outage performance of the OFDM-based cognitive AF relay network in direct mode for the data rates 1 b/s/Hz and 2 b/s/Hz. It is observed that when the data rate is increased from 1 b/s/Hz to 2 b/s/Hz, the minimum SNR requirement increases from 20 dB to 25 dB in the absence of NBI to operate the network at the maximum outage probability of 1 × 10
^{−2}
. The minimum SNR requirement rises to 30 dB SNR with a data rate of 1 b/s/Hz at an SIR of 20 dB.
Outage performance of proposed cognitive AF relay network in direct mode at 1 b/s/Hz and 2 b/s/Hz.
The outage performance for the proposed OFDM-based CR network in AF relay mode transmission is shown in
Fig. 4
at the data rate of 1 b/s/Hz for different SIR values. Since the maximal ratio combining provides diversity gain at the destination node, the minimum SNR requirement decreases to 15 dB from 20 dB in the direct mode transmission at the outage probability of 1 × 10
^{−2}
. In AF relay mode at 30 dB SNR, as the SIR increases from 20 dB to 25 dB, the outage probability decreases from 1 × 10
^{−3}
to 1 × 10
^{−4}
at the data rate of 1 b/s/Hz.
Outage performance of AF mode transmission at 1 b/s/Hz.
The outage performance of the proposed cognitive network in AF relay mode for the data rates of 1 b/s/Hz and 2 b/s/Hz is shown in
Fig. 5
. It is observed that when the data rate moves from 1 b/s/Hz to 2 b/s/Hz, the minimum SNR requirement to operate the network at the maximum outage probability of 1 × 10
^{−2}
increases from 15 dB to 22 dB in the absence of NBI, and for the SIR level of 20 dB, it further rises to 17 dB SNR with a data rate of 1 b/s/Hz and exceeds the limit at a data rate of 2 b/s/Hz.
Outage performance in AF relay mode at 1 b/s/Hz and 2 b/s/Hz.
The outage performance at the direct and AF relay modes is compared in
Fig. 6
. It is observed that the AF mode provides 5 dB SNR gain compared to the direct mode in the absence of NBI. At 30 dB SNR, when the mode is switched from direct to AF relay, the outage probability decreases from 1 × 10
^{−2}
to 1 × 10
^{−3}
at a data rate of 1 bit/s/Hz with 20 dB SIR.
Comparison of outage performance in AF relay mode and direct mode at 1 b/s/Hz.
In
Fig. 7
, the outage performance of the proposed OFDM-based CR network in AF relay mode is compared with direct mode for various relay interference values. The data rate is 1 b/s/Hz. In this Figure, RY and DS denote the relay and destination nodes, respectively. It is noted that the outage probability decreases from 1 × 10
^{−3}
to 4 × 10
^{−4}
when relay node operates with no NBI, at 30 dB SNR per subcarrier. In direct mode, the outage probability of the network is 1 × 10
^{−2}
for an SIR of 20 dB and an SNR of 30 dB. It indicates that the proposed network with limited interference at the relay performs better than the direct mode.
Comparison of outage performance in AF relay and direct modes at 1 b/s/Hz, with and without NBI at relay.
In
Fig. 8
,
r
represents the number of subcarriers that are directly affected by NBI. The outage performance of the proposed CR network in direct mode is shown in
Fig. 8
at a data rate of 1 b/s/Hz for different
r
values. At 30 dB SNR, when
r
decreases from five to three, the outage probability decreases from 4.4 × 10
^{−2}
to 2.7 × 10
^{−2}
; furthermore, it reduces to 1 × 10
^{−2}
when
r
= 1.
Outage performance of proposed cognitive AF relay network in direct mode at 1 b/s/Hz for various r values.
The outage performance of the AF relay and incremental AF relay modes is compared in
Fig. 9
. It is observed that for the outage probability of 1 × 10
^{−4}
, the incremental AF relay mode provides 5 dB SNR gain compared to the AF relay mode, in the absence of NBI. This is due to the additional degree of freedom in the incremental AF relay mode.
Comparison of outage performance in AF relay and incremental AF relay modes at 1 b/s/Hz.
The outage performance of the AF relay mode for the data rates between 0.5 b/s/Hz and 3 b/s/Hz is shown in
Fig. 10
. At a low data rate, the outage is exactly reduced by 1 × 10
^{−1}
for every 5 dB increase in SIR level. For the maximum outage probability of 1 × 10
^{−2}
, the CR network supports a data rate of up to 2 b/s/Hz at 25 dB SIR.
Outage performance of AF relay mode vs. transmission rates for fixed SNR = 40 dB.
Corresponding Author ssraj.tce@gmail.com
Samikkannu Rajkumar received his BE degree in electronics and communication engineering from the Angala Amman College of Engineering and Technology, Tamil Nadu, India, in 2000 and his ME degree in communication systems from Thiagarajar College of Engineering, Tamil Nadu, India, in 2006. At present, he is working as a full-time research scholar with the Department of Electronics and Communication Engineering, Thiagarajar College of Engineering, Tamil Nadu, India. His current areas of research include developing algorithms in cooperative relay networks and multiple antenna systems.
vnsenthilkumaran@tce.edu
V.N. Senthilkumaran received his BE degree in electronics and communication engineering from VelTech Engineering College, Tamil Nadu, India, in 2002. He obtained his ME degree in communication systems from S.R.M. Engineering College, Tamil Nadu, India, in 2004 and his PhD degree in the area of wireless relay network design from Anna University, Tamil Nadu, India, in 2014. At present, he is working as an assistant professor with the Department of Electronics and Communication Engineering, Thiagarajar College of Engineering, Tamil Nadu, India. His current areas of research include physical-layer network coding and capacity analysis of wireless networks.
sjtece@tce.edu
S.J. Thiruvengadam received his BE degree in electronics and communication engineering from Thiagarajar College of Engineering, Tamil Nadu, India, in 1991. He obtained his ME degree in applied electronics from the College of Engineering, Tamil Nadu, India, in 1994 and his PhD degree in the area of signal detection algorithms from Madurai Kamaraj University, Tamil Nadu, India, in 2005. At present, he is working as a professor with the Department of Electronics and Communication Engineering, Thiagarajar College of Engineering. Previously, he was a visiting associate professor with the Department of Electrical Engineering, Stanford University, CA, USA, under a postdoctoral fellowship sponsored by the Department of Science and Technology, Government of India, from January to December 2008. His current areas of research include statistical signal processing and multiple-input and multiple-output wireless communications.

I. Introduction

Cognitive radio (CR) uses the available radio spectrum opportunistically to cope with the growing number of bandwidth-limited wireless services
[1]
. In a CR network, secondary users effectively utilize the unused spectrums of the primary users. Orthogonal frequency-division multiplexing (OFDM) is employed in a CR system due to its inherent characteristics of spectrum sensing, spectrum shaping, scaling, and interoperable capabilities
[2]
. Further, OFDM reduces fading caused by multiple receptions of a signal and hence improves spectral efficiency.
However, the performance of an OFDM-based CR system is affected by narrowband interference (NBI) from other wireless users. The following interference cancellation techniques are widely used in wireless communication systems: filter-based approach, cyclostationarity-based approach, higher-order statistics–based approach, transform domain approach, joint detection/multiuser detection (MUD), and spatial processing (for example, beamforming)
[3]
. Filter-based approaches were used to suppress NBI in OFDM-based conventional CR networks in
[4]
–
[6]
. In a coded OFDM system, NBI is suppressed by whitening the NBI spectrum using a prediction-error filter
[4]
. In
[5]
, an interfering carrier frequency is estimated and cancelled using a notch filter. The success of this is dependent upon the accuracy of the estimation of the interfering carrier frequency. In
[6]
, the active interfering tones in a multiband OFDM CR network are turned off and the resultant inter-carrier interference is suppressed using a notch filter. A higher-order statistics–based Bayesian decision theoretic approach is used to suppress partial-band interference in a frequency hop communication system
[7]
. However, this cannot be directly extended to an OFDM-based high-data-rate communication system due to its computational complexity. In
[8]
, NBI is estimated by measuring the output signal energy, and receiver windowing is proposed to reduce the effect of NBI in an OFDM system.
In
[9]
, using a linear minimum mean square error–based estimator, unmodulated subcarriers are employed to measure NBI power. As the energy of NBI is different at each subcarrier, it is considered as a colored noise, which is then compensated for using an iterative decoder based on an expectation–maximization algorithm
[10]
. A constrained minimum mean-output-energy–based algorithm is developed for NBI estimation under the assumption that the second-order statistics of a received signal are known at the receiver
[11]
. The major issue with this method is its high computational complexity. Compressive sensing is used to estimate an NBI signal by exploiting its sparsity in the frequency domain
[12]
. The effects of NBI are analyzed in an ultra-wideband cooperative relay network in
[13]
.
Joint detection and MUD-based approaches were used for detecting secondary user symbols in the presence of a primary user signal, at the cost of high complexity
[14]
. In
[15]
and
[16]
, interference is minimized in a CR system using the methods of opportunistic interference cancellation and phase adjustment, respectively. Amplify-and-Forward (AF) and Decode-and-Forward (DF) relaying protocols are used in cooperative communication networks to improve throughput and quality of service by exploiting spatial diversity
[17]
. The success of such a spatial processing technique has inspired us to use a cooperative AF half-duplex relay for interference cancellation. In this paper, the first major contribution is the introduction of a half-duplex AF relaying protocol in an OFDM-based CR network to suppress NBI. The use of such a relaying protocol provides spatial diversity and maximizes the output signal-to-interference-plus-noise ratio (SINR) of the CR network at each subcarrier of the OFDM signal. The second major contribution of this work is the derivation of the analytical expressions for analyzing the outage performance of the proposed OFDM-based cognitive AF relay network in the presence of NBI.
The rest of the paper is organized as follows. A system model for an OFDM-based cognitive AF relay network is described in both the time domain and the frequency domain for direct and AF relay modes, in Section II. The analytical expressions are derived for the SINR of the direct and AF relay modes, in Section III. The outage performance of an OFDM-based CR network is analyzed for direct, AF-relay, and incremental relaying modes, in Section IV. Numerical results are presented in Section V, and concluding remarks are given in Section VI.
II. System Model

Consider the OFDM-based cognitive AF relay network with three secondary users shown in
Fig. 1
. Such a network can operate in two modes; namely, direct mode and AF relay mode. In direct mode, SU1 transmits data to SU3 in the first time slot. In AF relay mode, SU1 transmits data to SU2, which acts as a relay, during the first time slot. The signal at SU2 is amplified and forwarded to SU3 in the second time slot. The frequency bands FB1, FB2, and FB3, which are available for transmission/reception at each secondary user, are also indicated in
Fig. 1
. The received signals at SU2 in relay mode and SU3 in direct/relay modes are affected by NBI from other sources in the scenario. It is assumed that the first secondary user, SU1, has to transmit an OFDM signal to the third secondary user, SU3, through a frequency-selective fading channel. The OFDM signal has
PPT Slide

Lager Image

- 1. Direct Mode

In direct mode, SU1 transmits the OFDM signal to SU3 using the band FB1. The
H 13 TD

and
H b3 TD

are
F zp = [ F N 0 N×v ] Η ); x 1 FD

is a
b 3 TD

is a
σ b3 2

;
z 3 TD

is a
σ z3 2

. The carrier frequency offset (CFO) for an NBI signal in band FB1 is represented by a
A fo,3 (1),TD =diag{ 1,exp( j 2πβ(1)/P ) , ... , exp( j 2πβ(P−1)/P ) },

where
y 13 TD

in (1) is given by
H 13 TD

and
H b3 TD

are circulant matrices, they can be decomposed as
H 13 TD = F P H H 13 FD F P

and
H b3 TD = F P H H b3 FD F P

[20]
, where
H 13 FD

and
H b3 FD

are diagonal matrices. Let
H 13 FD

and
H b3 FD

are defined as
H 13 FD

= diag {
H b3 FD

= diag {
A 3 (1),FD =F P A3 (1),TD F P H

is a circulant matrix of CFO at SU3; and
b 3 FD

and
z 3 FD

are the NBI signal and noise vectors in the frequency domain, respectively. The
b 3 FD

is sparse in nature and indicates that
b 3 FD

contains only a few non-zero elements. The
- 2. AF Relay Mode

In AF relay mode, SU1 transmits an OFDM signal to SU2 using the band FB2 in the first time slot. The
y 12 TD

at SU2 is given by
H 12 TD

and
H b2 TD

are the
b 2 TD

and
z 2 TD

are the NBI signal and noise vectors at SU2, respectively, and
A fo,2 (2),TD

is the CFO of the NBI signal in band FB2. By taking a DFT, the frequency-domain representation of the received vector
y 12 TD

is given by
H 12 TD

and
H b2 TD

are circulant matrices, (5) can be simplified as
H 12 FD

and
H b2 FD

are the diagonal matrices of the vectors
A 2 (2),FD =F P A2 (2),TD F P H

is a circulant matrix.
In the second time slot, the node of SU2 plays as an AF relay. It amplifies the received signal
y 12 TD

and forwards it to SU3 using the band FB3. The
y 23 TD

at SU3 is given by
ρ= Ρ s / ( Ρ s σ h 2 + σ b2 2 + σ z2 2 )

is the amplification factor at SU2;
σ h 2 =E[ | h 12 | 2 ]

is the instantaneous squared-norm of the CIR vector
h 12 ; P s =E[ | x 1 FD | 2 ]

is the transmitted source power, which is the same as in SU2;
σ b2 2

is the IPPS of the NBI signal vector
b2 TD

;
σ z2 2

is the variance of the noise vector
z 2 TD

; and
H 23 TD

and
H b3TD

are the circulant matrices of the CIR vectors
A fo,3 (3),TD

,
b 3TD

, and
n 3TD

represent the CFO for NBI in band FB3 and the NBI signal and noise vectors at SU3, respectively. Substituting (4) in (7), the received vector
y 23TD

is written as
y 23TD

is given by
H 23 TD

,
H 12 TD

,
H b2 TD

, and
H b3TD

are circulant matrices. By decomposing them, (9) can be simplified as
H 23FD

and
H b3FD

are the vectors
A 3 (3),FD =F P A3 (3),TD F P H

is a circulant matrix.
III. SINR Analysis

In this section, the analytical expressions for the SINR of the received signals at relay and destination nodes are derived at each subcarrier in both direct and AF relay modes.
- 1. Direct Mode

From (3), the direct-mode receive signal
y 13FD

at the
A 3 (1),FD

is a circulant matrix, the third term is written as a circular convolution. The fourth term represents the noise at the
y 13FD

at the
∑ q=0 N−1 | v(k,q) | 2 =1,

the term
∑ q=0,q≠k N−1 | v(k,q) | 2

can be simplified as 1 − |
ς 3 = ∑ m=0,m∈ I b3 P−1 | a 3 1 ( k, (m−k) P ) | 2

is the squared-norm of the CFO term and
σ b3 2 =E( ∑ m=0,m∈ I b3 P−1 | b 3 (m) | 2 )

is the IPPS of the NBI term. Now, substituting (14) and (15) in (12), the lower bound of Γ
σ z3 2

= E(|
σ b3 2

and SNR = P
σ z3 2

.
- 2. AF Relay Mode

In AF relay mode, a transmit signal passes through the two different channels
y 12FD

from SU1 to SU2 at the
ς 2 = ∑ m=0,m∈ I b2 P−1 | a 2 2 ( k, (m−k) P ) | 2

is the squared-norm of the CFO term and
σ b2 2 =( ∑ m=0,m∈ I b2 P−1 | b 2 (m) | 2 )

is the IPPS of the NBI term. Substituting (21) and (22) in (20), we get
σ b2 2

and SNR = P
σ z2 2

, and
σ z2 2

is the noise variance. The instantaneous SINR of the received signal
y 23FD

can be expressed as
ς 3 ″ = ∑ m=0,m∈ I b3 P−1 | a 3 3 ( k, (m−k) P ) | 2

is the squared-norm of the CFO and
σ b3 2 =( ∑ m=0,m∈ I b3 P−1 | b 3 (m) | 2 )

is the IPPS of the NBI term. Substituting (25) and (26) in (24), we get
SIR3= P s / σ b3 2 , SNR=P s / σ n3 2

and
σ n3 2 =E( | n 3 ( k ) | 2 )

is the noise variance. In AF relay mode, the maximum data transfer rate depends on the minimum value between Γ
IV. Performance Analysis

In this section, the outage performance of the OFDM-based cognitive AF relay network is analyzed at the subcarrier level in direct and AF relay modes for the data rate of
- 1. Outage Probability Analysis in Direct Mode

In direct mode, the upper bound of the outage probability of Γ
P out,k dir (R)

can be written as
X = Y = | λ 13 (k) | 2 ,Z= ∑ m=0,m∈ I b3 P−1 | λ b3 (m) | 2 ,a= | v(k,k) | 2 ,b=( 1− | v(k,k) | 2 ) , c= ς 3 ( 1/ SIR1 )

, and
- 2. Outage Probability Analysis in AF Relay Mode

The outage probability for AF relay–mode SINR, Γ
g 2 =( 2 SNR + ς 2 SIR2 + ς 3 ″ SIR3 ).

Then, it can be written as
- 3. Outage Probability for Incremental AF Relaying

In AF relay mode, a relay transmits the amplified received signal at the data rate of
V. Numerical Results

The outage performance of the proposed OFDM-based CR relay network is analyzed in the presence of an NBI signal for direct, AF relay, and incremental AF relay modes. The simulation parameters are listed in
Table 1
.
Simulation parameters of cognitive AF relay network.

S. No | Parameters | Value |
---|---|---|

1 | Number of subcarriers in OFDM signal, | 64 |

2 | Guard sequence length, | 16 |

3 | Number of channel taps, | 10 |

4 | Number of subcarriers with NBI, | 1, 3, 5 |

5 | CFO parameter of NBI signal, | 0.5 |

6 | Data rates (in b/s/Hz), | 0.5, 1, 2, 3 |

7 | Signal to NBI interference ratio (in dB) | 10, 15, 20, 25 |

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

VI. Conclusion

In this paper, a fixed AF relaying protocol is introduced in an OFDM-based CR network to suppress the effect of NBI. The performance of the CR network with an AF relay is analyzed in terms of outage probability. It is shown that the use of the AF relay suppresses the effect of NBI and that the outage performance is improved by 5 dB per subcarrier compared to the direct-mode transmission at a data rate of 1 b/s/Hz. Further, the outage performance is improved by the incremental AF relaying, where the system is operated at variable rate nature. It is found that for no NBI at the relay node, the outage of the system is significantly reduced. Finally, we conclude that the effect of the NBI can be reduced by the AF relay protocol in an OFDM-based cognitive radio system at constant data rate.
BIO

Haykin S.
2005
“Cognitive Radio: Brain-Empowered Wireless Communications,”
IEEE J. Sel. Areas Commun.
23
(2)
201 -
220
** DOI : 10.1109/JSAC.2004.839380**

Mahmoud H.A.
,
Yücek T.
,
Arslan H.
2009
“OFDM for Cognitive Radio: Merits and Challenges,”
IEEE Wireless Commun.
16
(2)
6 -
15

Hong X.
2009
“Cognitive Radio Networks: Interference Cancellation and Management Techniques,”
IEEE Veh. Technol. Mag.
4
(4)
76 -
84

Batra A.
,
Zeidler J.R.
“Narrowband Interference Mitigation in OFDM Systems,”
IEEE Military Commun. Conf.
San Diego, CA, USA
1 -
7

Coulson A.J.
2006
“Bit Error Rate Performance of OFDM in Narrowband Interference with Excision Filtering,”
IEEE Trans. Wireless Commun.
5
(9)
2484 -
2492
** DOI : 10.1109/TWC.2006.1687772**

Yamaguchi H.
“Active Interference Cancellation Technique for MB-OFDM Cognitive Radio,”
IEEE European Microw. Conf.
Amsterdam, Netherlands
Oct. 12–14, 2004
1105 -
1108

Baum C.W.
,
Pursley M.B.
1992
“Bayesian Methods for Erasure Insertion in Frequency-Hop Communication Systems with Partial Band Interference,”
IEEE Trans. Commun.
40
(7)
1231 -
1238
** DOI : 10.1109/26.153368**

Redfern A.J.
2002
“Receiver Window Design for Multicarrier Communication Systems,”
IEEE J. Sel. Areas Commun.
20
(5)
1029 -
1036
** DOI : 10.1109/JSAC.2002.1007383**

Nilsson R.
,
Sjoberg F.
,
LeBlanc J.
2003
“A Rank-Reduced LMMSE Canceller for Narrowband Interference Suppression in OFDM-Based Systems,”
IEEE Trans. Commun.
51
(12)
2126 -
2140
** DOI : 10.1109/TCOMM.2003.820761**

Sanguinetti L.
,
Morelli M.
,
Poor H.V.
“Robust EM-Based Detection of BICM-OFDM Transmissions in the Presence of Narrowband Interference,”
IEEE European Wireless Conf.
Lucca, Italy
Apr. 12–15, 2010
696 -
700

Darsena D.
2007
“A Constrained Maximum-SINR NBI-Resistant Receiver for OFDM Systems,”
IEEE Trans. Signal Process.
55
(6)
3032 -
3047
** DOI : 10.1109/TSP.2007.893946**

Gomaa A.
,
Al-Dhahir N.
2011
“A Sparsity-Aware Approach for NBI Estimation in MIMO-OFDM,”
IEEE Trans. Wireless Commun.
10
(6)
1854 -
1862
** DOI : 10.1109/TWC.2011.040411.101118**

Maichalernnukul K.
,
Kaiser T.
,
Zheng F.
“Performance of Ultra-wideband Cooperative Relay Systems in the presence of Narrowband Interference,”
IEEE Int. Conf. Ultra-wideband
Bologna, Italy
Sept. 14–16, 2011
405 -
409

Hassan M.H.
,
Hossain Md. J.
2013
“Cooperative Beamforming for Cognitive Radio Systems with Asynchronous Interference to Primary User,”
IEEE Trans. Wireless Commun.
12
(11)
5468 -
5479
** DOI : 10.1109/TWC.2013.092013.121275**

Popovski P.
“Opportunistic Interference Cancellation in Cognitive Radio Systems,”
IEEE Dynamic Spectr. Access Netw.
Dublin, Ireland
Apr. 17–20, 2007
472 -
475

Alian E.H.M.
,
Mitran P.
2013
“A Phase Adjustment Approach for Interference Reduction in OFDM-Based Cognitive Radios,”
IEEE Trans. Wireless Commun.
12
(9)
4668 -
4679
** DOI : 10.1109/TW.2013.072213.121994**

Laneman J.N.
,
Tse D.N.C.
,
Wornell G.W.
2004
“Cooperative Diversity in Wireless Networks: Efficient Protocols and Outage Behavior,”
IEEE Trans. Inf. Theory
50
(12)
3062 -
3080
** DOI : 10.1109/TIT.2004.838089**

Muquet B.
2002
“Cyclic Prefixing or Zero Padding for Wireless Multicarrier Transmissions?”
IEEE Trans. Commun.
50
(12)
2136 -
2148
** DOI : 10.1109/TCOMM.2002.806518**

Batra A.
,
Zeidler J.R.
“Narrowband Interference Mitigation in BICM OFDM Systems,”
IEEE Acoust. Speech Signal Process.
Taipei, Taiwan
Apr. 19–24, 2009
2605 -
2608

Gray R.M.
2006
Toeplitz and Circulant Matrices: A Review
Now Publishers Inc.
Hanover, MA, USA
155 -
239

Rabiei P.
,
Namgoong W.
,
Al-Dhahir N.
2011
“On the Performance of OFDM-Based Amplify-and-Forward Relay Networks in the Presence of Phase Noise,”
IEEE Trans. Commun.
59
(5)
1458 -
1466
** DOI : 10.1109/TCOMM.2011.030411.100206**

Proakis J.G.
,
Salehi M.
2008
“Digital Communication,”
Mc-Graw Hill
New York, USA
1 -
1150

Lehmann E.L.
1966
“Some Concepts of Dependence,”
Ann. Mathematical Statistics
37
(5)
1137 -
1153
** DOI : 10.1214/aoms/1177699260**

Citing 'Outage Analysis of OFDM-Based Cognitive AF Relay Network in the Presence of Narrowband Interference
'

@article{ HJTODO_2015_v37n3_460}
,title={Outage Analysis of OFDM-Based Cognitive AF Relay Network in the Presence of Narrowband Interference}
,volume={3}
, url={http://dx.doi.org/10.4218/etrij.15.0114.0300}, DOI={10.4218/etrij.15.0114.0300}
, number= {3}
, journal={ETRI Journal}
, publisher={Electronics and Telecommunications Research Institute}
, author={Rajkumar, Samikkannu
and
Senthilkumaran, V.N.
and
T, hiruvengadam}
, year={2015}
, month={May}