A spectrum sensing scheme for partially polarized waves

A spectrum holes sensing approach is developed for Cognitive Radio (CR) systems based on the estimation of the Stokes parameters of partially polarized electromagnetic waves. The performance of the proposed approach is compared with the energy detection method. An expression for the probability of detection based on the proposed method is derived using α - μ generalized gamma fading model.


Introduction
Cognitive radio (CR) technology has witnessed a growing interest over the past decade, as it promises more efficient use of the available spectrum [1,2]. A key stage in CR is spectrum sensing, in which the Secondary User (SU) must detect the presence of a Primary User (PU) in a certain channel, and thus, deems this part of the spectrum unused, and make the decision to share it. This entails a sequence of functions that the CR system should perform, such as power control [3] and spectrum management [4]. Several techniques were proposed to improve spectrum sensing such as energy detection [5], cyclostationary feature detection [6], sensing based on smart antennas [7,8], and wideband spectrum sensing [9,10]. These techniques primarily make use of the amplitude, frequency, and phase information of the PU signal.
It is possible, however, to improve the spectrum sensing process by exploiting the polarization state of the signal. In radar systems [11], the polarization state was used to improve the detection capability of the system. The new polarization-dependent detection statistics, which use the power and relative phase of the two orthogonal polarization components was proposed to enhance radar detection in homogenous channels [12]. The radar detection performance was enhanced based on the polarization difference between the clutter and the target [11]. The sine of the relative phase between two orthogonally polarized received signals has been proposed and tested as detection statistic in radar systems [13]. New statistics use two orthogonal polarization component powers and their relative phase to enhance target detection [12]. They are thus fundamentally different from the well-known Marcum-Swerling envelope detector [14], which operates on only a single polarization component of the received power, and the pseudocoherent detector [13], which operates on only the relative phase of the polarization components.
Recent research on polarization detection in CR was concerned with completely polarized waves. A Virtual Polarization Detection (VPD) method based on the vector signal processing was presented for effective spectrum sensing of cognitive radios [15]. Polarization spectrum hole sensing was proposed for cognitive radio to optimize the received polarization at the SU in order to protect the PU from interference caused by the SU and to reduce the interference from PU to SU [16]. Optimal Polarization Reception (OPR) was proposed for CR to improve the SINR [17]. A new blind spectrum sensing method based on the polarization characteristic of the received signal, which is completely represented by the orientation of a polarization vector, was proposed [18], and a closed-form expression for the probability of false alarm and probability of detection under Additive White Gaussian Noise (AWGN) and Rayleigh-fading channels was derived. The fading and noisy nature of a wireless communication channel places a major challenge for spectrum sensing. Since sensing decisions based on a single SU measurements may be unreliable, the idea of collaborative spectrum sensing has attracted a lot of research interest [19]. However, the polarization state is changed with spectrum sensing time which was not considered in this work.
In this paper, the above constraint is addressed by proposing a polarization-based spectrum sensing scheme where a statistical model for the PU signal polarization parameters is adopted. This model takes into account the channel backscatter and the partial polarization nature of the PU signal, with the assumption that the channel experiences slow fading. Cooperative polarization spectrum sensing is proposed to mitigate the effects of fading and shadowing, which can seriously degrade the sensing performance. A cluster-based cooperation scheme is proposed to decrease bit error probability. All SUs are grouped into few clusters and one cluster head is set for each cluster to collect the sensing results, make cluster decisions, and forward measurements to the central unit. Thus the bit error probability will be reduced greatly because most of SUs will be closer to the cluster heads than to the central unit. Analytical results show that significant improvement can be achieved with our proposed method.
The rest of this paper is organized as follows. In Section 2, the method of characterizing the signal polarization is provided. The spectrum sensing models are developed in Section 3. Section 4 describes the cooperative spectrum sensing mechanism and the cluster technique of the cooperative sensing is shown in Section 5. Simulation results are illustrated in Section 6, followed by the conclusion in Section 7.

The Method of Characterizing the Signal Polarization
The polarization of a monochromatic plane wave is completely specified by constant amplitude and relative phase of the two orthogonal electric-field components. In a cognitive radio system, the SU uses orthogonally dual polarized antennas to detect the PU signal ⃗ , which can be completely described in vector form as where and indicate the in-phase and quadrate phase components, respectively. Typically, the polarization of a monochromatic wave is determined by the Jones vector ⃗ of the signal ⃗ defined as [20] Alternatively, the polarization is determined by the geometrical parameters, namely the ellipticity angle , and orientation angle , which can be uniquely represented by a point on the Poincare sphere. The field component parameters and , and the geometrical parameters and , are described in the following set of equations [21] sin (2 ) cos ( ) = sin (2 ) cos (2 ) , However, if the amplitudes and phases encounter slow fluctuations with time, which is the case if the polarization of primary signal suffers from either noise or fading, the components of the Jones vector are said to be quasimonochromatic or narrowband. The polarization of a quasimonochromatic wave is quantified using an average polarization state vector, which may be defined in terms of four measurable components; namely, ⃗ = [ 0 1 2 3 ] . These components are known as the Stokes parameters (SPs) and are given by [12] where is the sampling time, and is the number of samples which is restricted by the bandwidth of a narrowband filter. 0 , 1 , 2 , and 3 are physically recognizable quantities as follows.
(1) 0 is the sum of the power in the ℎ and V electric-field components and thus represents the total power of the received signal. (2) 1 is the difference between the power in the ℎ and V electric-field components. (3) 2 is the difference between the power in the two orthogonal electric-field components whose axes are rotated 45 ∘ relative to the ℎ and V axes.
(4) 3 is the difference between the right-hand and the left-hand circularly polarized power.
The component 0 satisfies the relation Alternatively, the average polarization state can be obtained from the coherence matrix C of the received field, such that where ⟨⋅⟩ denotes averaging over and * represents the complex conjugate. The coherence matrix C is a linear combination of the SPs such that: where

Spectrum Sensing Based on the Polarization
Spectrum sensing is essentially a binary hypothesis testing problem, which indicates the PU's absence or presence, respectively, such that where ( ) is the observed signal at the CR, ( ) is the PU signal, ( ) is the zero mean Gaussian random process with identical autocorrelation and power spectral density 0 Watts/Hz, and ℎ is the amplitude gain of the channel having mean-square value Ω = ℎ 2 and Probability Density Function (PDF) ℎ (ℎ). The received instantaneous signal power is modulated by ℎ 2 and consequently the instantaneous Signalto-Noise Ratio (SNR) can be expressed as = ℎ 2 ( / 0 ) with an average = Ω( / 0 ) where is the signal energy accumulated over the observation period.
In this paper, we propose two schemes which use the degree of polarization and the axial ratio to detect the PU signal.

Spectrum Sensing Based on the Degree of Polarization.
The degree of polarization ( ) is a quantity used to describe the polarized portion of an electromagnetic wave. A perfectly polarized wave has equal to 1, whereas an unpolarized wave has equal to 0. A wave, which is partially polarized, can therefore be represented by superposition of a polarized and unpolarized component. This implies a somewhere in between 0 and 1. Alternatively, the degree of polarization is defined as the ratio of the polarized power to the total power in the wave; that is, The estimate of the ratio of the polarized power to the total power in the received signal is used as detection statistics in radar systems [12]. Figure 1 depicts a block diagram of the proposed spectrum sensing system where = . There are two orthogonal antennas, which detect the horizontal and vertical components ℎ and V of the signal , respectively. The Stokes vector estimation block delivers the Stokes vector ⃗ to the polarization degree estimator, which, in turn, produces the detection statistic . This statistic serves as the input to the threshold detector to decide whether the PU is present or not.
The distribution that governs the statistic can be used to estimate the probability of detection and the probability of false alarm. Here, the in-phase and quadrature components of the quasimonochromatic wave are assumed to be zeromean Gaussian random processes. The estimation of the SP and the elements of the sample correlation random process is equivalent, which is called Wishart distribution [22]. Thus, the probability density function of the detection statistic is given by [12] where ∞ and { 0 , 1 , 2 , 3 } are the actual degree of polarization, and the actual Stokes vector components, respectively. They are obtained in the limit as the number of independent samples approaches infinity. For a fixed threshold , the conditional probability of false alarm and detection can be expressed as [12] where is the Signal-to-Noise Ratio (SNR), 1− 2 ( + 1/2, ) is the incomplete beta function, = 2Γ( + 1/2)/√ , 4 ISRN Signal Processing Figure 1: The proposed spectrum sensing system model.
can be affected by fading that in turn affects and hence in order to incorporate its influence, must be averaged over all possible values of according to where (⋅) represents the PDF of the channel. For a fading signal with envelope ℎ, an arbitrary parameter > 0, and a -root mean value ℎ = √ (ℎ ), the -PDF ℎ (ℎ) is given by [23] where = E 2 (ℎ )/V(ℎ ) and E(⋅) and V(⋅) are the expectation and variance operators, respectively [23]. The PDF of the SNR is obtained by a change of variables as shown in [24] ( ) = ( /2) By substituting (15) in (13), can be written as where Using [25], one can get the average probability of detection as follows: where is the Hypergeometric function, (⋅, :) is the Meijer function, and is the Hypergeometric function [26].

Spectrum Sensing Based on the Axial Ratio of a Polarization Ellipse.
The polarized portion of the PU signal represents a net polarization ellipse traced by the electric field vector as a function of time. The ellipse has a magnitude ( ) such that The ellipticity is the ratio of the minor to the major axis of the corresponding electric field polarization ellipse and varies ISRN Signal Processing 5 from 0 for linearly polarized wave to 1 for circularly polarized wave. The polarization ellipse is alternatively described by its eccentricity, which is zero for a circularly polarized wave, and increases as the ellipse becomes thinner. It then becomes one for a linearly polarized wave. Alternatively is defined as the ratio of the polarized power to the total power in the wave. Figure 1 depicts a block diagram of the proposed spectrum sensing system, where = . This statistic serves as the input to the threshold detector to decide whether the PU is present or not. The distribution that governs the statistic is to be determined to estimate the probability of detection and the probability of false alarm [12], which is given by For a fixed threshold , the conditional probability of false alarm can be expressed as and the probability of detection can be given by must be averaged over all possible values of as follows: where ( ) represents the -probability density function. Following a similar approach to that of Section 3.1 and by substituting (15) in (24), can be written as where Using binomial theorem and [25], will be given by It is worth mentioning that the restriction for this scheme is that the polarization information of the primary signal must be known a priori. If the PU signal is linearly polarized, then will be close to 0 for high SNR, and if the PU signal is circularly polarized, then will be close to 1 for high SNR.

Spectrum Sensing Based on Cooperative Polarization Detection
The fading and noisy nature of a wireless communication channel places a major challenge on the accuracy of spectrum sensing. Sensing decisions that is based on measurements of a single SU may be unreliable. Cooperative spectrum sensing is one possible solution to overcome this unreliability. In Figure 2, a number of SUs, which are distributed in different locations independently, can detect the PU and make the decision whether the signal exists or not. According to the information received from various SUs, the central unit makes the final decision based on some rules. A general fusion rule is when a final decision of 1 is taken whenout-of -SUs report 1. When = 1, the -out-of -rule is equivalent to the OR rule. When = , the decision rule becomes the AND rule. By selecting different values of , different detection performances are obtained.
Two cases are considered, namely, cooperative spectrum sensing with perfect and imperfect reporting channels.

Perfect Reporting Channels.
If the channels between each SU and the central unit are noise free, then the overall probability of false alarm and the overall probability of detection of the cooperative spectrum sensing for -outof -rule of the cooperative spectrum sensing are given by [27] where the , and , are the probabilities of false alarm and detection of the degree of polarization and axial ratio as derived in (14), (18)

Imperfect Reporting Channels.
In practical systems, the reporting channels between the SUs and the central unit will experience fading. This, in turn, will degrade transmission reliability of the sensing results that are reported from the SUs to the central unit. Let , denote the probability of receiving 1 at the central unit when the th Secondary User sends 0 and the probability of receiving 0 at the central unit when the th Secondary User sends 1 which can be calculated over -fading channel as follows: where (1/2) erfc( √ ) is the probability of error for BPSK in an channel. The shown integral can be put in the form of Laplace transform and hence using [25] where Δ( , ) = / , ( + 1)/ , . . . , ( + − 1)/ .
Consequently, the overall probability of false alarm is obtained as where CU 1 , PU 1 , and SU 1 are the hypothesis in the central unit, PU, and SU, respectively.

Cluster-Based Cooperative Spectrum Sensing
The clustering method is proposed in the cooperative spectrum sensing scheme in order to improve the sensing performance by decreasing the reporting channel error, which is proved to exploit a selection diversity gain [28]. In this paper, we assume that the instantaneous channel state information of the channel is available for the Secondary Users, and the channel between any two Secondary Users in the same cluster is perfect because they are in the vicinity of each other. Figure 3 shows the system model of cluster-based cooperative spectrum sensing. We assume that there are SUs, who are divided into clusters. The th cluster has an integer number of SUs, which satisfies Firstly, all SUs are assumed to belong to few clusters by the distributed clustering algorithms [29]. Secondly, the Secondary User who has the largest instantaneous reporting channel gain will be selected as the cluster head. Then, the cooperative sensing is carried out through the following steps. Firstly, all SUs perform the local spectrum sensing. Every SU sends a decision to the cluster head. Secondly, the cluster head receives those local decisions from the SUs in the same cluster and then makes the decision according to certain fusion rule. All cluster heads send their decisions to the central unit. In the end, the central unit makes the final decision according to the fusion rule. The cluster head and central unit make the decision according to the OR-rule in order to limit the interference from the SUs to the PU. The -out-of -fusion rule is adopted in the cluster head and the OR fusion rule is adopted in the central unit [30]. If each cluster has / SUs, then the global false alarm probability and detection probability will be given by where (V, , ) ) . (45)

Results and Discussion
Results obtained using the proposed methods are presented. Monte Carlo simulation consisting of 100000 independent trials was performed. The degree of freedom is set to 3.5 to make fair comparison between the proposed methods and the energy detection method [31]. The equal and normalized polarization states of the PUs are randomly generated with ∈ (0, /2), ∈ (− /2, /2). It is assumed that SUs use dual polarized antennas, which is necessary to detect the horizontal and vertical components of the wave. Figure 4 compares the Receiver Operating Characteristic (ROC) of the DoP method with the directional method [18]. It can be noticed from this figure that the DoP method improves the at low SNR since the signal is polarized and the noise is unpolarized wave. This makes the signal detection, in particular at low SNR, more efficient, which is the range of interest for CR systems. Therefore, it cannot be used in high SNR but it offers a baseline for comparison. Figure 5 shows the ROC under AWGN compared with the energy detection (ED) method. It is clear from this figure that the proposed method improves the probability of detection over a wide range of SNR (from −5 dB to 5 dB). For example, in the case of SNR = 0 dB, improved from 22% to 88%. Moreover, it is observed that at SNR = 0 dB for the DoP method is better than for the ED method at SNR = 5 dB, meaning the improvement of more than +5 dB for the DoP method relative to the ED method. It is also noticed that at SNR = 10dB, the DoP method offers 100% detection probability for all values of . This improvement is due to the more accurate detection of the polarized portion.

The Degree of Polarization (DoP) Method.
The case of Rayleigh fading channel is considered in Figure 6, which illustrates that fading has a stronger impact on the ED than on the proposed DoP method. This is attributed to the fact that the DoP is based on a second order statistic; that is, the improvement in the probability of detection for SNR = 5 dB is approximately 30% at = 0.2. A better performance for the DoP is reported for Nakagami fading as shown in Figure 7, where, at SNR = 0dB, the improvement is approximately 20% when the = 0.2. Moreover, the at SNR = 5 dB for the proposed method is greater than the for the ED method at SNR = 10 dB, this means that the improvement is nearly +7 dB.
It can be noticed that at fading conditions, the depolarization effect increases, which means higher dispersion of polarization states. Therefore, the detection performance decreases. Thus, the smaller the depolarization effect on the primary signal, the more constant the polarization state and thus the better the detection performance.

The Axial Ratio (AR) Method.
To demonstrate the performance of the axial ratio (AR) method, Figure 8 shows the characteristics for different values of and . For = 2 and = 2, the probability of detection is 91% at = 0.2, which is a superior value relative to the ED method. This can be easily explained since at strong fading less energy is coupled between the cross-polarized channels, which results in small dispersion of the polarization state leading to better detection performance. Figure 9 shows the superior performance of the AR method relative to the ED method over a wide range of SNR under Rayleigh fading channel. In Figure 10, where Nakagami fading is considered, for SNR = −5 dB, an improvement of approximately 60% in the probability of detection is achieved for = 0.2 and 70% for SNR = 5dB relative to the ED method. Hence, the AR method is the most robust and applicable detection method in the case of fading condition and/or presence of noise power uncertainty.         A comparison between both approaches is drawn in Figure 11 showing the superiority of the AR method relative to the DoP, where for = 0.2, the AR method gives = 0.73 compared to = 0.58 for the DoP. It can be seen that the sensing performance of the AR is better than the DoP method due to the fact that the AR method resists the depolarization effect on the primary signal. Figure 12 shows the ROC of the cooperative scheme for the DoP method over a wide range of SNR under imperfect Nakagami fading channel. It is clear that the is very high at low especially at SNR = 5 dB and 10 dB. It can be observed that the increase rapidly with the increase of the SNR. When the SNR = 10 dB, the is almost 100% for all values of .

Conclusions
Spectrum sensing based on Stokes parameters was thoroughly analyzed using new detection statistics, namely, the degree of polarization and the axial ratio. The proposed approaches were studied under different fading scenarios and the obtained results demonstrated superior performance relative to the conventional energy detection method. An extensive study is reported on the two methods and their performance under -fading channels. The results demonstrated that the proposed algorithms are particularly applicable for the case of unknown primary polarization and/or presence of noise power uncertainty. In general, the proposed methods are better than energy detection in different ranges of the SNR and under different fading conditions. They are, however, more complex in terms of implementation due to the use of two orthogonally polarized antennas. It is also worth noting that the axial ratio method is generally better than the degree of polarization method.
Cooperative spectrum sensing was then considered and shown to be a powerful method for dealing with the hidden terminal problem. Simulation results show that the cluster rule gives superior performance for a wide range of SNR.