# Patent

## Title
Method for OFDM signal identification and parameter estimation

## Abstract

Spectrum sensing is one of the most challenging problems in cognitive radio systems. The spectrum of interest needs to be characterized and unused frequencies should be identified for possible exploitation in a simple and fast way, allowing the radio to catch up with the changing transmission parameters. A sensing method is presented where primary users are identified by matching the features extracted from the received signal to the a priori information about primary users' transmission characteristics. For estimation of some signal parameters, the cyclostationarity of the transmission spectrum is explored by using a suboptimal maximum likelihood (ML) estimator. The proposed algorithms can be used in cognitive radio for identifying various transmissions and for electronic surveillance.

## Background
CROSS REFERENCE TO RELATED APPLICATION
This application claims priority to U.S. Provisional Patent Application No. 60/954,115, entitled “OFDM Signal Identification and Transmission Parameter Estimation for Cognitive Radio Applications”, filed Aug. 6, 2007. The contents of which are herein incorporated by reference.


FIELD OF INVENTION
This invention relates to spectrum sensing for cognitive radio applications. Specifically, the invention relates to a method for signal identification and parameter estimation.
BACKGROUND OF THE INVENTION
As wireless communication systems are making transition from voice-only communications to interactive Internet data and multi-media applications, the desire for higher data rate transmissions are increasing tremendously. With the increase in wireless devices, future technologies will face spectral crowding, and coexistence of wireless devices will be a major problem. Considering the limited bandwidth availability, accommodating the demand for higher capacity and data rates is a challenging task, requiring innovative technologies that can offer new ways of exploiting the available radio spectrum.
Cognitive radio is a viable solution for spectrum scarcity problem that arises with increased number of users and applications that require higher data rates such as video and Internet by enabling opportunistic spectrum usage (Mitola, J. and J. Maguire, G. Q., “Cognitive radio: making software radios more personal,” IEEE Personal Commun. Mag., vol. 6, no. 4, pp. 13-18, August 1999). One crucial aspect of cognitive radio is related to its ability to autonomously identify the available unused spectrum and to adapt to environmental and spectral characteristics of a geographical location.
In cognitive radio, primary users have higher priority or legacy rights than secondary users on the spectrum usage. Conversely, secondary users exploit the idle spectrum as to not interfere with primary users. Therefore, secondary users need to have the capability to sense the spectrum usage reliably. Spectrum sensing requires obtaining the spectrum characteristics across multiple dimensions such as time, space, frequency, and code.
Matched filtering is the optimum method for detection of primary users. However, traditional matched filtering requires the cognitive radio to demodulate the received signal and requires perfect knowledge of the primary users signaling features. Moreover, since cognitive radio needs receivers for all signal types, traditional matched-filtering is practically difficult to implement. Spectrum sensing can also be performed by correlating the received signal with a known copy of itself, known as waveform-based sensing (H. Tang, “Some physical layer issues of wide-band cognitive radio systems,” in Proc. IEEE mt. Symposium on New Frontiers in Dynamic Spectrum Access Networks, Baltimore, Md., USA, November 2005, pp. 151-159). However, this method is only applicable to systems with known repeating patterns such as wireless metropolitan area network (WMAN) signals (IEEE Standard for Local and Metropolitan area networks Part 16, The Institute of Electrical and Electronics Engineering, Inc. Std. IEEE 802.16E-2005, 2005). Another method for detection of primary user transmission is cyclostationarity feature detection, where cyclostationarity features of the received signal are exploited (D. Cabric and R. W. Brodersen, “Physical layer design issues unique to cognitive radio systems,” in Proc. IEEE mt. Symposium on Personal, Indoor and Mobile Radio Commun., vol. 2, Berlin, Germany, September 2005, pp. 759-763). These features are caused by the periodicity in the signal or in its statistics (mean, autocorrelation etc.). Instead of power spectral density (PSD), cyclic correlation function is used for detecting the signals present in a given spectrum. Cyclostationarity based methods can differentiate noise from the primary users. Finally, energy detector-based spectrum sening is common because of low hardware complexity, robust performance, and blind primary users' signal characterization. The presence of a signal can be detected by comparing the output of the energy detector with a threshold which depends on the noise floor. However, energy detection suffers from an inability to differentiate interference from primary users' signals and poor performance under low signal-to-noise ratios (SNRs) (H. Tang, “Proc. IEEE mt. Symposium on New Frontiers in Dynamic Spectrum Access Networks, Baltimore, Md., USA, November 2005, pp. 151-159).
Cognitive radio devices must identify unused signal spectrum in a fast and efficient way. Conventional algorithms sense the spectrum without exploiting all the properties of primary users. The a priori information about transmission properties of possible primary users can be leveraged for developing a framework for spectrum sensing.
Current spectrum estimation techniques do not utilize most of the a priori knowledge about primary users. This knowledge about center frequencies and bandwidth of this type of a signal along with other transmission characteristics can be exploited to identify the presence of a transmission and thus to improve spectrum sensing accuracy. The features that describe devices can be obtained by investigating the received signal. These features then can be matched to the a priori sets of known parameters. By finding the exact transmission characteristics of primary users, possible estimation errors due to sensing algorithm and noise can be removed.
SUMMARY OF THE INVENTION
Multi-carrier techniques, specifically orthogonal frequency division multiplexing (OFDM), are commonly used in modern communications systems. For identification of active transmissions, signals can be classified as single-carrier and multi-carrier, reducing candidate size. Moreover, transmission parameters of an OFDM based system can be detected blindly if the system is not known to cognitive radio. The capability of identifying transmission parameters using blind detection can be useful for spectrum survey for radio monitoring systems to discover enemy and illegal transmissions as well.
To realize a fully cognitive radio, the cognitive devices should be aware of not only spectral content but also temporal and spatial contents of the environment that they are operating in. Information about the transmission properties of possible primary users, such as transmission bandwidths and center frequencies, are used to develop a signal identification method where the parameters estimated from the received signal. These estimated parameters are matched to the possible transmission parameters for achieving a more robust and reliable sensing. The sensing might also be performed by correlating the received signal with a known copy of itself. This method is only applicable to systems with known signal patterns such as wireless metropolitan area network (WMAN) signals, and it is termed as waveform-based sensing.
Another method for detection of primary user transmission is cyclostationarity feature detection. This algorithm exploits the cyclostationarity features of the received signal caused by the periodicity in the signal or in its statistics. Cyclic correlation function is used for detecting the signals present in the given spectrum. This method can also differentiate the noise from the primary users which is a result of the fact that noise is a wide-sense stationary (WSS) signal with no correlation while modulated signals are cyclostationaiy with spectral correlation due to the redundancy of signal periodicities. An energy detector-based approach is used as an intermediate step for obtaining features about the transmissions in the frequencies of interest. Energy detector based approaches (also called radiometry or periodogram) are the most common way of spectrum sensing because of their low computational complexity and because the receiver does not need any knowledge on the primary users' signal.
In one embodiment, waveform spectrum sensing is performed by utilizing a method of extracting signal features from a test frequency. A series of useful algorithms include energy detector or cyclostationary feature extraction fast Fourier transform, time frequency transform, discrete Fourier transform, and inverse discrete Fourier transform. These extracted signal features are used to identify use of a transmission spectrum by a primary user. The radio exploits a primary user's transmission spectrum use in characterizing the multi-dimensional spectrum.
Feature extraction in the spectrum sensing utilizes estimates of symbol duration and cyclic prefix size using known OFDM transmission characteristics in some embodiments. Symbol length estimates are used to calculate cyclic prefix duration and timing offset. Alternatively, the cyclic prefix duration is estimated by assuming cyclic prefix duration is multiple of symbol duration, from 1⁄4, 1⁄8, 1/16, or 1/32. Once the cyclic prefix is known, the repetition in cyclic prefix is used to find data duration.
To determine primary use of the transmission spectrum, the square of the fast Fourier transform output is compared to a threshold. Frequency domain features and time domain features are estimated in some embodiments to reduce cognitive radio hardware requirements. These features include bandwidth, center frequency, single-carrier transmission, multi-carrier transmission, maximum signal duration, multiple accessing, duplexing, frame duration, and spreading codes. The feature sets may be classified using Baysean decisions. In some embodiments, the method is a function of the received signal evaluated at the maximum likelihood estimates of signal parameters.
Also presented is a method of estimating waveform parameters for transmission spectrum. A carrier transmission must first be identified as either multi-carrier or single carrier, followed by a determination of the radio access technology (RAT) classification. Symbol duration and cyclic prefix length are identified using estimations based on the radio access technology. Symbol length is estimated to calculate cyclic prefix duration and timing offset, which are alternatively assumed to be a multiple of symbol duration, from 1⁄4, 1⁄8, 1/16, or 1/32. The number and frequencies of subcarriers are then estimated. In some embodiments, the method relies on a modified likelihood function for obtaining performance results of the proposed method. A Bayesian classification is used to assign features of RAT to the method.



BRIEF DESCRIPTION OF THE DRAWINGS
For a fuller understanding of the invention, reference should be made to the following detailed description, taken in connection with the accompanying drawings, in which:
FIG. 1 is a block diagram indicating the three main steps required in to model a cognitive engine in the present invention.
FIG. 2 is an illustration of the cylic prefix extension for OFDM systems.
FIG. 3 is a graph comparing the numerical performance of ML, suboptimal ML, and two-step algorithms. Normalized MSE of OFDM symbol duration estimation versus SNR is also shown.
FIG. 4 is a graph indicating the probability of incorrect detection for OFDM signal for different numbers of symbols under AWGN and fading channels.
FIG. 5 is a graph showing normalized MSE of OFDM symbol duration as a function of SNR for differing numbers of symbols under AWGN and fading channels.
FIG. 6 is a block diagram of the energy detector-based partial match-filtering method of the present invention.
FIG. 7 is a graph of receiver operating characteristics curves for different SNR values of a single frequency sample.
FIG. 8 is an illustration of a DUDE algorithm for bandwidth and center frequency estimation.
FIG. 9 is a graph showing power spectral density of tested WLAN and Bluetooth signals.
FIG. 10 is a graph indicating detection error rates for WLAN and Bluetooth systems at different SNR values for Lt=15.
FIG. 11 is a graph indicating detection error rates for WLAN and Bluetooth systems at different SNR values for Lt=30.
FIG. 12 is a graph indicating detection error rates for WLAN and Bluetooth systems at different SNR values for Lt=15 under a multi-path fading channel.



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT
Spectrum sensing is realized in three main steps: 1) extraction of a predefined set of features from the received signal; 2) using the extracted features for making decisions on the presence of an anticipated transmission; and 3) exploiting the gained knowledge about the active primary users for multi-dimensional spectrum characterization.
The system model for cognitive radio is shown in FIG. 1, indicating the three main steps. The cognitive engine governs the selection of bands for transmission and communicates with radio, user, and open systems interconnection reference (OSI) layers. Radio identification is performed based on the extracted features of received signal, environmental information and cognitive radio's prior knowledge. The set of possible radio access technologies (RATs) and their transmission parameters can be collected by cognitive devices using previous decisions (blind) or they can be broadcasted by a central unit (assisted). Alternatively, these parameters can be pre-configured to cognitive radio during hardware design.
In feature extraction, parameters are extracted using algorithms such as energy detector or cyclostationary based feature extraction. These algorithms can be used individually or together for estimating different parameters. These parameters include frequency domain features, such as bandwidth, center frequency, and single carrier or multi-carrier, time domain features like maximum signal duration, multiple accessing, duplexing method, frame duration, spreading codes, and other features, such as cumulants, power moment matrices, and zero-crossing rate.
The dimension of the feature set is a design parameter. Some or all of the parameters can be extracted from the received data and used for classification. However, the larger the dimension of the feature set, the more complex the system, and the more reliable the classification results. Selection of the employed features is performed by considering the characteristics of potential systems. The feature set can be represented as a vector x which is a point in multidimensional feature space. Then, this vector is used for classifying the detected transmission into one of K candidate transmissions using a classifier.
In decision making, the measured signal is associated with a RAT, using classification methods such as pattern recognition, neural networks, or statistical classification. Selection of a classification method depends on the considered RATs and available features. The classification performed in this step is different from blind modulation classification task because not only the type of signaling but also some parameters of the transmission need to be identified.
Multi-dimensional spectrum characterization, the output of classification is used for obtaining a complete multidimensional spectrum awareness in cognitive radio. The knowledge of primary users helps identify the transmission opportunities across different dimensions. For example, if the identified signal is a cordless phone, the range is expected to be around 100 meters and for Bluetooth signals it is around 10 meters. This knowledge can be used in a cooperative sensing environment for gaining knowledge in the space dimension. The characterization in time, frequency, and code dimensions is straightforward once primary users are associated with a particular RAT.
The feature extraction algorithms are based on FFT operation which is used to transform time-domain signal into frequency domain, using FFT output samples to decide whether an FFT frequency bin is occupied by a primary user. However, the primary user's signal is usually spread over a group of FFT output samples as the bandwidth of primary user is expected to be larger than the considered bandwidth divided by the FFT size. Using this fact, the FFT output is filtered for noise averaging in order to obtain a better performance.
The proposed algorithms are especially suitable for cognitive devices using OFDM as their transmission technology. The availability of FFT circuitry in these systems eases the requirements on the hardware. Moreover, the computational requirements of the spectrum sensing algorithm is reduced as receiver needs to apply FFT for data detection.
OFDM signals can be considered as a composition of large number of independent random signals. Hence, sampled OFDM signals can be assumed to have a Gaussian distribution thanks to central limit theorem. This observation is used to classify signals as single-carrier or multi-carrier by applying normality tests to incoming signal. However, in the presence of multipath-fading and interference, single-carrier signals at the cognitive radio receiver can also have Gaussian distribution. This is especially true for low signal-to-noise ratio (SNR) scenarios. Furthermore, when power control is employed in time-multiplexed OFDM (such as IEEE 802.16), the Gaussian approximation does not hold anymore, and as a consequence, classification based on normality tests may fail.
Various parameter estimation algorithms for OFDM systems exist. Methods for finding the OFDM symbol and cyclic prefix (CP) durations are mostly based on exploitation of the cyclostationarity of OFDM signaling due to cyclic prefix extension. To estimate the length of CP, time-frequency transform of the received signal is found and entropy is calculated, which is then used to find the CP size. Auto-correlation of the received signal with its delayed copy is calculated and its magnitude is normalized with average signal power. As the auto-correlation has non-zero terms only for the part that falls into the guard time, power ratio of the guard time to the symbol duration can be used to find the CP length. However, SNR knowledge is required for a precise estimation.
Subcarrier or discrete Fourier transform (DFT) size is another key parameter for OFDM systems. An algorithm based on multiple signal classification (MUSIC) technique is used to determine DFT size using eigen-decomposition to diagonalize the correlation matrix of the received signal Akaike's criteria is applied to estimate the number of sub-spaces or sub-carriers. However, zeros are used for generating OFDM signal instead of cyclic prefix. Once CP and symbol durations are estimated, DFT size is determined by testing DFT size values. DFT outputs are tested for Gaussianity to find the correct DFT size. If not all of the subcarriers are employed for transmission, active subcarriers should be identified. This can be performed by analyzing the power level at the output of FFT. However, this method requires knowledge of FFT size and synchronization to the received signal.
OFDM converts serial data stream into parallel blocks of size NFFT and modulates these blocks using inverse discrete Fourier transform (IDFT). Thus, the OFDM signaling system was modeled. The continuous-time OFDM signal in the baseband can be written as:




  
    
      
        
          
            
              x
              m
            
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              
                
                  ∑
                  
                    k
                    =
                    1
                  
                  K
                
                ⁢
                
                  
                
                ⁢
                
                  
                    
                      S
                      m
                    
                    ⁡
                    
                      (
                      k
                      )
                    
                  
                  ⁢
                  
                    e
                    
                      j
                      ⁢
                      
                        
                          2
                          ⁢
                          
                            πΓ
                            ⁡
                            
                              (
                              k
                              )
                            
                          
                          ⁢
                          t
                        
                        
                          T
                          D
                        
                      
                    
                  
                
              
              ⁢
              
                
              
              -
              
                T
                G
              
            
            ≤
            t
            <
            
              T
              D
            
          
        
        ,
      
    
    
      
        (
        1
        )
      
    
  




where Sm(k) is the transmitted data symbol at kth (data) subcarrier of mth OFDM symbol and Γ denotes the set of K used subcarriers. Useful data duration and CP length are represented by TD and TG respectively and they make up the total duration of OFDM symbol TS, i.e. TS=TD+TG, as shown in FIG. 2. Transmitted signal can now be written as:





  
    
      
        
          x
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            ∑
            
              m
              =
              
                -
                ∞
              
            
            ∞
          
          ⁢
          
            
              
                x
                m
              
              ⁡
              
                (
                
                  t
                  -
                  
                    mT
                    S
                  
                
                )
              
            
            .
          
        
      
    
    
      
        (
        2
        )
      
    
  




After D/A conversion and radio frequencies (RF) modulation, signal is passed through the mobile radio channel which can be modeled as a time-variant linear filter





  
    
      
        
          
            h
            ⁡
            
              (
              t
              )
            
          
          =
          
            
              ∑
              
                i
                =
                1
              
              L
            
            ⁢
            
              
            
            ⁢
            
              
                
                  h
                  i
                
                ⁡
                
                  (
                  t
                  )
                
              
              ⁢
              
                δ
                ⁡
                
                  (
                  
                    t
                    -
                    
                      τ
                      i
                    
                  
                  )
                
              
            
          
        
        ,
      
    
    
      
        (
        3
        )
      
    
  




where L is the number of sample-spaced channel taps and τi is the delay of each tap. The channel is assumed to be constant over an OFDM symbol, but time-varying across OFDM symbols, which is a reasonable assumption for low and medium mobility, is possible.

At the receiver, the signal is received along with noise. After down-conversion, simplified baseband model of the received signal can be formulated as:




  
    
      
        
          y
          ⁡
          
            (
            t
            )
          
        
        =
        
          
            
              e
              
                j2πζ
                ⁢
                
                  
                
                ⁢
                t
              
            
            ⁡
            
              [
              
                
                  x
                  ⁡
                  
                    (
                    t
                    )
                  
                
                ⁢
                ★
                ⁢
                
                  
                
                ⁢
                
                  h
                  ⁡
                  
                    (
                    t
                    )
                  
                
              
              ]
            
          
          +
          
            ω
            ⁡
            
              (
              t
              )
            
          
        
      
    
    
      
        (
        4
        )
      
    
  
  
    
      
        
          
        
        ⁢
        
          =
          
            
              
                e
                
                  j2πζ
                  ⁢
                  
                    
                  
                  ⁢
                  t
                
              
              ⁢
              
                ∫
                
                  
                    x
                    ⁡
                    
                      (
                      τ
                      )
                    
                  
                  ⁢
                  
                    h
                    ⁡
                    
                      (
                      
                        t
                        -
                        τ
                      
                      )
                    
                  
                  ⁢
                  
                    d
                    τ
                  
                
              
            
            +
            
              ω
              ⁡
              
                (
                t
                )
              
            
          
        
      
    
    
      
        (
        5
        )
      
    
  
  
    
      
        
          
        
        ⁢
        
          
            =
            
              
                
                  e
                  
                    j2πζ
                    ⁢
                    
                      
                    
                    ⁢
                    t
                  
                
                ⁢
                
                  
                    ∑
                    
                      i
                      =
                      1
                    
                    L
                  
                  ⁢
                  
                    
                  
                  ⁢
                  
                    
                      x
                      ⁡
                      
                        (
                        
                          t
                          -
                          
                            τ
                            i
                          
                        
                        )
                      
                    
                    ⁢
                    
                      
                        h
                        i
                      
                      ⁡
                      
                        (
                        t
                        )
                      
                    
                  
                
              
              +
              
                ω
                ⁡
                
                  (
                  t
                  )
                
              
            
          
          ,
        
      
    
    
      
        (
        6
        )
      
    
  




Where ★ denotes convolution operation and ω(t) is the additive white Gaussian noise (AWGN) sample with zero mean and variance of σ2ω. ξ corresponds to the frequency offset due to inaccurate frequency synchronization with the received signal. The average received signal power is denoted by σ2ξ.

It is assumed that the receiver does not know the correct sampling rate and samples the received signal with a sampling time of Δt. Now, the discrete-time received signal can be obtained from continuous time signal as:

y[n]:=y(nΔt).  (7)

In discrete-time representation, ND, NG and NS are used to represent the number of samples corresponding to data, CP and total symbol durations respectively. The time durations can be calculated by multiplying these numbers with the sampling period, e.g. TS=NSΔt. The received signal is assumed to have D samples which contains an unknown single-carrier signal or unknown number of OFDM symbols. In vector notation the received signal can be represented as:

y=[y[1],y[2], . . . y[D]].  (8)

Mobile WiMax

The recently approved 802.16e standard (known as mobile WiMAX) uses scalable OFDMA as physical layer transmission technique. By changing the FFT size as a function of the transmission bandwidth, the subcarrier spacing is kept constant for all bandwidths aiming to reduce the inter-carrier interference (Id), due to mobility and frequency offsets, to a negligible level. The available FFT sizes and system bandwidths are given in the Table. Moreover, the fixed standard supports FFT sizes of 256 and 2048 with different bandwidths (multiple of 250 kHz and no less than 1.25 MHz). Hence the MSs should have support for various channel bandwidths. This requires dynamic detection of the FFT size and the channel bandwidth which is a similar problem to detection the primary users in a given frequency band and the algorithm proposed in the previous section might be used for this purpose. In both cases, some a priori information about the center frequencies and bandwidths of the expected signals is assumed to be available.







TABLE





Available bandwidth and FTT sizes for 802.16e.




















System bandwidth (MHz)
1.25
NA
5
10
20


Sampling frequency (MHz)
1.429
NA
5.714
11.429
22.857


FFT size
128
256
512
1024
2048









OFDM Cyclostationary Detection

OFDM signaling is shown to be cyclostationary with period TS in (R. Heath Jr and G. Giannakis, “Exploiting input cyclostationarity for blind channel identification in OFDM systems,” IEEE Trans. Signal Processing, vol. 47, no. 3, pp. 848-856, March 1999). It is also known that cyclostationarity of a signal holds when it is passed through a linear time-invariant channel. Hence, the received OFDM signal is also cyclostationary. Detection algorithms using the cyclostationarity characteristics of the OFDM signaling in AWGN channel were developed.
Symbol duration and CP size are estimated by assuming an OFDM transmission. For the observation vector y, OFDM symbol duration ND and CP size NG are not known as well as the start time of each symbol, denoted by θ. Let us denote the set of samples that belong to CP with Ω. The samples in CP and in the last part of OFDM symbol are pairwise correlated, as seen in FIG. 2. Using the structure of OFDM signaling, the correlation of the samples in observation vector can be written as:




  
    
      
        
          E
          ⁡
          
            [
            
              
                y
                ⁡
                
                  [
                  n
                  ]
                
              
              ⁢
              
                
                  y
                  *
                
                ⁡
                
                  [
                  
                    n
                    +
                    Δ
                  
                  ]
                
              
            
            ]
          
        
        =
        
          {
          
            
              
                
                  
                    σ
                    s
                    2
                  
                  +
                  
                    σ
                    ω
                    2
                  
                
              
              
                
                  
                    Δ
                    =
                    0
                  
                  ,
                
              
              
                
                  
                
              
            
            
              
                
                  
                    σ
                    s
                    2
                  
                  ⁢
                  
                    e
                    
                      
                        -
                        j2πξΔ
                      
                      ⁢
                      
                        
                      
                      ⁢
                      
                        tN
                        D
                      
                    
                  
                
              
              
                
                  
                    Δ
                    =
                    
                      N
                      D
                    
                  
                  ;
                
              
              
                
                  n
                  ∈
                  Ω
                
              
            
            
              
                0
              
              
                
                  otherwise
                  .
                
              
              
                
                  
                
              
            
          
        
      
    
    
      
        (
        9
        )
      
    
  




where E [·] denotes expectation operation.

Maximum likelihood estimates of the unknown parameters ND, NG, and θ can be obtained by maximizing the log-likelihood function of the observations y. Log-likelihood function be formulated as:




  
    
      
        
          Λ
          ⁡
          
            (
            
              
                y
                ;
                
                  N
                  D
                
              
              ,
              
                N
                G
              
              ,
              θ
            
            )
          
        
        =
        
          log
          (
          
            
              ∏
              
                n
                ∈
                Ω
              
            
            ⁢
            
              
            
            ⁢
            
              
                p
                ⁡
                
                  (
                  
                    
                      y
                      ⁡
                      
                        [
                        n
                        ]
                      
                    
                    ,
                    
                      y
                      ⁡
                      
                        [
                        
                          n
                          +
                          
                            N
                            D
                          
                        
                        ]
                      
                    
                  
                  )
                
              
              ⁢
              
                
                  ∏
                  
                    n
                    ∉
                    
                      Ω
                      ⋃
                      
                        Ω
                        ′
                      
                    
                  
                
                ⁢
                
                  
                
                ⁢
                
                  p
                  ⁡
                  
                    (
                    
                      y
                      ⁡
                      
                        [
                        n
                        ]
                      
                    
                    )
                  
                
              
            
          
          )
        
      
    
    
      
        
          
        
        ⁢
        
          (
          10
          )
        
      
    
  
  
    
      
        =
        
          log
          (
          
            
              ∏
              
                n
                ∈
                Ω
              
            
            ⁢
            
              
            
            ⁢
            
              
                
                  p
                  ⁡
                  
                    (
                    
                      
                        y
                        ⁡
                        
                          [
                          n
                          ]
                        
                      
                      ,
                      
                        y
                        ⁡
                        
                          [
                          
                            n
                            +
                            
                              N
                              D
                            
                          
                          ]
                        
                      
                    
                    )
                  
                
                
                  
                    p
                    ⁡
                    
                      (
                      
                        y
                        ⁡
                        
                          [
                          n
                          ]
                        
                      
                      )
                    
                  
                  ⁢
                  
                    p
                    ⁡
                    
                      (
                      
                        y
                        ⁡
                        
                          [
                          
                            n
                            +
                            
                              N
                              D
                            
                          
                          ]
                        
                      
                      )
                    
                  
                
              
              ⁢
              
                
                  ∏
                  n
                
                ⁢
                
                  
                
                ⁢
                
                  p
                  ⁡
                  
                    (
                    
                      y
                      ⁡
                      
                        [
                        n
                        ]
                      
                    
                    )
                  
                
              
            
          
          )
        
      
    
    
      
        (
        11
        )
      
    
  




where p (·) is the probability distribution function (PDF) of the observed samples and p (·, ·) is the joint PDF. The set Ω′ contains the samples of OFDM data that are copied to CP. Assuming that y has a jointly Gaussian distribution, it is shown that the log-likelihood function can be obtained as:





  
    
      
        
          Λ
          ⁢
          
            (
            
              
                y
                ;
                
                  N
                  D
                
              
              ,
              
                N
                G
              
              ,
              θ
            
            )
          
        
        =
        
          
            
            
              
                ∑
                
                  n
                  ∈
                  Ω
                
              
              ⁢
              
                
              
              ⁢
              
                
                  y
                  ⁡
                  
                    (
                    n
                    )
                  
                
                ⁢
                
                  
                    y
                    *
                  
                  ⁡
                  
                    (
                    
                      n
                      +
                      
                        N
                        D
                      
                    
                    )
                  
                
              
            
            
          
          -
          
            
          
          ⁢
          
            
              ∑
              
                n
                ∈
                Ω
              
            
            ⁢
            
              
                (
                
                  
                    
                      a
                      2
                    
                    ⁢
                    
                      (
                      
                        
                          
                            
                            
                              y
                              ⁡
                              
                                [
                                n
                                ]
                              
                            
                            
                          
                          2
                        
                        +
                        
                          
                            
                            
                              y
                              ⁡
                              
                                [
                                
                                  n
                                  +
                                  
                                    N
                                    D
                                  
                                
                                ]
                              
                            
                            
                          
                          2
                        
                      
                      )
                    
                  
                  +
                  
                    
                      1
                      
                        2
                        ⁢
                        a
                      
                    
                    ⁢
                    
                      (
                      
                        
                          σ
                          s
                          2
                        
                        +
                        
                          σ
                          ω
                          2
                        
                      
                      )
                    
                    ⁢
                    
                      (
                      
                        1
                        -
                        
                          a
                          2
                        
                      
                      )
                    
                    ⁢
                    
                      log
                      ⁡
                      
                        (
                        
                          1
                          -
                          
                            a
                            2
                          
                        
                        )
                      
                    
                  
                
                )
              
              .
            
          
        
      
    
    
      
        (
        12
        )
      
    
  




where a and Ω are as given in





  
    
      
        a
        =
        
          
            
              σ
              s
              2
            
            
              
                σ
                s
                2
              
              +
              
                σ
                ω
                2
              
            
          
          .
        
      
    
    
      
        (
        13
        )
      
    
  




Ω:={n|m(ND+NG)+p+θ, 





  
    
      
        
          m
          =
          
            {
            
              0
              ,
              1
              ,
              ...
              ⁢
              
                
              
              ,
              
                
                  ⌊
                  
                    D
                    
                      
                        N
                        D
                      
                      +
                      
                        N
                        G
                      
                    
                  
                  ⌋
                
                -
                1
              
            
            }
          
        
        ,
      
    
    
      
        (
        14
        )
      
    
  




p={1,2, . . . , NG}, θ={0,1, . . . , ND−1}}.

respectively. Please note that, ML estimation obtained by maximizing the likelihood function given by (12) is asymptotically optimal when the following assumptions are satisfied: 1) Noise is white, Gaussian, and uncorrelated with the transmission signal, 2) y has a joint Gaussian distribution with correlation function given by (9). Also note that, the ML estimation algorithm can be modified to estimate the frequency offset ξ as well. The likelihood function of the received signal is evaluated at ML estimates of signal parameters, providing an estimation of signal parameters via rotational invariance techniques (ESPRIT).
The value of a, σ2s and σ2n are difficult to estimate in a non-cooperative scenario. At low SNR values, the second and third terms of (12) approaches to zero and can be dropped. Therefore, a suboptimal ML estimation is possible with less complexity and no need for the knowledge of noise and signal variances. Modified likelihood function can now be obtained from (12) as:




  
    
      
        
          
            Λ
            ~
          
          ⁡
          
            (
            
              
                y
                ;
                
                  N
                  D
                
              
              ,
              
                N
                G
              
              ,
              θ
            
            )
          
        
        =
        
          
            
            
              
                ∑
                
                  m
                  =
                  0
                
                
                  ⌊
                  
                    D
                    
                      
                        N
                        D
                      
                      +
                      
                        N
                        G
                      
                    
                  
                  ⌋
                
              
              ⁢
              
                
                  ∑
                  
                    p
                    =
                    1
                  
                  
                    N
                    G
                  
                
                ⁢
                
                  
                    y
                    ⁡
                    
                      [
                      
                        
                          m
                          ⁡
                          
                            (
                            
                              
                                N
                                D
                              
                              +
                              
                                N
                                G
                              
                            
                            )
                          
                        
                        +
                        p
                        +
                        θ
                      
                      ]
                    
                  
                  ⁢
                  
                    
                      y
                      *
                    
                    ⁡
                    
                      [
                      
                        
                          m
                          ⁡
                          
                            (
                            
                              
                                N
                                D
                              
                              +
                              
                                N
                                S
                              
                            
                            )
                          
                        
                        +
                        
                          N
                          D
                        
                        +
                        p
                        +
                        θ
                      
                      ]
                    
                  
                
              
            
            
          
          .
        
      
    
    
      
        (
        15
        )
      
    
  




where |·| denoted the largest integer less than or equal to input argument. Note that the likelihood function is independent of the frequency offset. Hence, the proposed sub-obtimum ML estimation algorithm is





  
    
      
        
          
            N
            ^
          
          D
        
        ,
        
          
            N
            ^
          
          G
        
        ,
        
          θ
          =
          
            
              
                arg
                ⁢
                
                  
                
                ⁢
                max
              
              
                
                  N
                  D
                
                ,
                
                  N
                  G
                
                ,
                θ
              
            
            ⁢
            
              {
              
                
                  Λ
                  ~
                
                ⁡
                
                  (
                  
                    
                      y
                      ;
                      
                        N
                        D
                      
                    
                    ,
                    
                      N
                      G
                    
                    ,
                    θ
                  
                  )
                
              
              }
            
          
        
      
    
    
      
        (
        16
        )
      
    
  
  
    
      
        
          
        
        ⁢
        
          =
          
            
              
                arg
                ⁢
                
                  
                
                ⁢
                max
              
              
                
                  N
                  D
                
                ,
                
                  N
                  G
                
                ,
                θ
              
            
            ⁢
            
              
                {
                
                  
                  
                    
                      ∑
                      
                        m
                        =
                        0
                      
                      
                        ⌊
                        
                          D
                          
                            
                              N
                              D
                            
                            +
                            
                              N
                              G
                            
                          
                        
                        ⌋
                      
                    
                    ⁢
                    
                      
                        ∑
                        
                          p
                          =
                          1
                        
                        
                          N
                          G
                        
                      
                      ⁢
                      
                        
                          y
                          ⁡
                          
                            [
                            
                              
                                m
                                ⁡
                                
                                  (
                                  
                                    
                                      N
                                      D
                                    
                                    +
                                    
                                      N
                                      G
                                    
                                  
                                  )
                                
                              
                              +
                              p
                              +
                              θ
                            
                            ]
                          
                        
                        ⁢
                        
                          
                            y
                            *
                          
                          ⁡
                          
                            [
                            
                              
                                m
                                ⁡
                                
                                  (
                                  
                                    
                                      N
                                      D
                                    
                                    +
                                    
                                      N
                                      S
                                    
                                  
                                  )
                                
                              
                              +
                              
                                N
                                D
                              
                              +
                              p
                              +
                              θ
                            
                            ]
                          
                        
                      
                    
                  
                  
                
                }
              
              .
            
          
        
      
    
    
      
        (
        17
        )
      
    
  




Maximization given in (17) is a computationally demanding process. In order to reduce the complexity, maximization can be performed in two steps. First, OFDM symbol length ND is estimated, and second the estimated ND value is used to find CP duration NG and timing offset θ. The simplification on the computational complexity comes with the expense of performance degradation.
The repetition due to CP is used to find the OFDM data duration TD. As the receiver is assumed to have no prior information on the symbol duration, all expected values should be tested. The discrete correlation of the received signal can be written as:




  
    
      
        
          
            R
            y
          
          ⁡
          
            (
            Δ
            )
          
        
        =
        
          
            1
            
              D
              -
              Δ
            
          
          ⁢
          
            
              ∑
              
                n
                =
                1
              
              
                D
                -
                Δ
              
            
            ⁢
            
              
                y
                ⁡
                
                  [
                  n
                  ]
                
              
              ⁢
              
                
                  
                    y
                    *
                  
                  ⁡
                  
                    [
                    
                      n
                      +
                      Δ
                    
                    ]
                  
                
                .
              
            
          
        
      
    
    
      
        (
        18
        )
      
    
  




Assuming enough statistical averaging, the correlation can be obtained as:




  
    
      
        
          
            R
            y
          
          ⁡
          
            (
            Δ
            )
          
        
        =
        
          {
          
            
              
                
                  
                    σ
                    s
                    2
                  
                  +
                  
                    σ
                    ω
                    2
                  
                
              
              
                
                  
                    Δ
                    =
                    0
                  
                  ,
                
              
            
            
              
                
                  
                    
                      N
                      G
                    
                    
                      
                        N
                        D
                      
                      +
                      
                        N
                        G
                      
                    
                  
                  ⁢
                  
                    σ
                    s
                    2
                  
                  ⁢
                  
                    e
                    
                      
                        -
                        j
                      
                      ⁢
                      
                        
                      
                      ⁢
                      2
                      ⁢
                      
                        
                      
                      ⁢
                      πξΔ
                      ⁢
                      
                        
                      
                      ⁢
                      
                        tN
                        D
                      
                    
                  
                
              
              
                
                  
                    Δ
                    =
                    
                      N
                      D
                    
                  
                  ,
                
              
            
            
              
                0
              
              
                
                  otherwise
                  ,
                
              
            
          
        
      
    
    
      
        (
        19
        )
      
    
  




Hence, the auto-correlation based estimation of ND can be formulated as:




  
    
      
        
          
            
              N
              ^
            
            D
          
          =
          
            
              
                arg
                ⁢
                
                  
                
                ⁢
                max
              
              Δ
            
            ⁢
            
              {
              
                
                
                  
                    R
                    y
                  
                  ⁡
                  
                    (
                    Δ
                    )
                  
                
                
              
              }
            
          
        
        ,
        
          Δ
          >
          0.
        
      
    
    
      
        (
        20
        )
      
    
  




Once the OFDM symbol duration ND is found, the CP duration can be estimated by testing for different NG values as:




  
    
      
        
          
            N
            ^
          
          G
        
        ,
        
          
            θ
            ^
          
          =
          
            
              
                arg
                ⁢
                
                  
                
                ⁢
                max
              
              
                
                  N
                  G
                
                ,
                θ
              
            
            ⁢
            
              
                {
                
                  
                    Λ
                    ~
                  
                  ⁡
                  
                    (
                    
                      
                        y
                        ;
                        
                          
                            N
                            ^
                          
                          D
                        
                      
                      ,
                      
                        N
                        G
                      
                      ,
                      θ
                    
                    )
                  
                
                }
              
              .
            
          
        
      
    
    
      
        (
        21
        )
      
    
  




Usually, the CP length is chosen as a multiple of OFDM symbol duration. In IEEE 802.16 systems, for example, available sizes are 1⁄4, 1⁄8, 1/16 and 1/32 (IEEE Standard for Local and Metropolitan area networks Part 16: Air Interface for Fixed and Mobile Broadband Wireles Access Systems Amendment 2: Physical and Medium Access Control Layers for Combined Fixed and Mobile Operation in Licensed Bands and Corrigendum 1, The Institute of Electrical and Electronics Engineering, Inc. Std. IEEE 802.16E-2005, 2005). Therefore, in order to further reduce the complexity, the CP size hypotheses can be limited to these ratios of the estimated ND value from previous step.
A numerical comparison of ML, suboptimal ML, and two-step low-complexity algorithms is given in FIG. 3 which shows the normalized mean-squared-error (MSE) of OFDM data duration estimation for AWGN and multi-path fading channels for a WLAN system. As can be seen from this figure, the performance loss for the suboptimal method is around 1-1.5 dB as compared to the ML method. Please note that for ML algorithm, perfect noise and signal variance knowledge is assumed; and the performance of ML algorithm is expected to decrease when these parameters are not perfectly known. Another 2 dB performance loss is observed when low-complexity two step method is used instead of calculating (17).
The modified likelihood function (15), using the two-step estimation algorithm, evaluated at the estimated values can be written as:

γ={tilde over (∇)}(y;ND,NG,θ)|{circumflex over (N)}D,{circumflex over (N)}G,θ={tilde over (∇)}(y;{circumflex over (N)}D,{circumflex over (N)}G,θ).  (22)

Eqn. (22) gives us the likelihood of receiving y for an OFDM system which has the estimated parameters. For OFDM signals, γ converges to




  
    ⌊
    
      D
      
        
          
            N
            ^
          
          D
        
        +
        
          
            N
            ^
          
          G
        
      
    
    ⌋
  
  ⁢
  
    
      N
      ^
    
    G
  
  ⁢
  
    σ
    s
    2
  




and for non-OFDM signals, it becomes zero, i.e.





  
    
      
        γ
        =
        
          {
          
            
              
                
                  
                    ⌊
                    
                      D
                      
                        
                          
                            N
                            ^
                          
                          D
                        
                        +
                        
                          
                            N
                            ^
                          
                          G
                        
                      
                    
                    ⌋
                  
                  ⁢
                  
                    
                      N
                      ^
                    
                    G
                  
                  ⁢
                  
                    σ
                    s
                    2
                  
                
              
              
                
                  H
                  1
                
              
            
            
              
                0
              
              
                
                  
                    H
                    0
                  
                  ,
                
              
            
          
        
      
    
    
      
        (
        23
        )
      
    
  




where 0 and 1 represent the null hypothesis and alternate hypothesis, i.e. the received signal is not OFDM signal or it is an OFDM signal, respectively. This observation is used for classification of a signal as single-carrier or multi-carrier (OFDM) by comparing γ to a threshold ξ. The received signal is classified as OFDM signal if γ>ξ.

The value of decision threshold depends on the incoming (OFDM) signal power and desired probability of false alarm (or probability of miss). Signal power can be found using the total signal if SNR is known. In non-cooperative environments, SNR knowledge is usually unavailable. In such a case, a low SNR value can be assumed.
The proposed algorithms were tested using computer simulations. As for OFDM system, IEEE 802.11a based WLAN was utilized. For this system, the transmission bandwidth is 20 MHz and number of subcarriers is 64, out of which 11 subcarriers are used as guard subcarriers and one as DC subcarrier. Furthermore, 8 subcarriers are deactivated in order to model a cognitive radio system. CP contains 16 samples, hence, CP rate is 1⁄4. For this system, CP and total symbol durations correspond to 0.8 μs and 4 μs respectively. The sampling rate of received signal is assumed to be unknown to the cognitive radio and signal is assumed to be sampled at a rate 50% higher than Nyquist rate. An exponentially decaying power delay profile (PDP) with 90 ns root-mean-squared (RMS) delay spread is used for modeling wireless channel. Simulation results are obtained for 25, 50 and 100 consecutive OFDM symbols. The threshold for classifying signal as OFDM is chosen as




  ζ
  =
  
    
      
        
          σ
          s
          2
        
        +
        
          σ
          ω
          2
        
      
      2
    
    ⁢
    
      ⌊
      
        D
        
          
            
              N
              ^
            
            G
          
          +
          
            
              N
              ^
            
            D
          
        
      
      ⌋
    
    ⁢
    
      
        
          N
          ^
        
        G
      
      .
    
  




This value corresponds to a 0 dB SNR value. Normalized MSE is used as a performance measure of the estimation as it reflects both the bias and the variance of the estimation.

FIG. 4 shows the probability of incorrect classification of the simulated OFDM signal as single-carrier in AWGN and multi-path fading channels. As this figure shows, the classification is successful even at low SNR values and no error is observed for 10 dB or higher SNR values. Note that for AWGN channel, 0 dB is a decision boundary for classification. This boundary is the selected threshold ξ which corresponds to 0 dB SNR assumption for the received signal. If SNR information is available, the decision threshold 4 can be changed adaptively decreasing the false alarm and mis-detection rates. The normalized MSEs for TD is given in FIG. 5. Similar results were obtained for MSE of CP time TG. The estimation is more accurate for AWGN channel than multi-path fading channel as expected. It is noted that the number of available OFDM symbols used for estimation affects the performance of the estimation because of noise averaging.
Partial Match Filtering
Partial match-filtering (PMF) matches the primary users' signals parameters instead of the signals. The block diagram of the proposed algorithm is shown in FIG. 6. The signal that arrives to cognitive user y(t) is first filtered with a band pass filter (BPF) to extract the signal in the frequencies of interest. This filter may be adjustable and controlled by a control unit in order to scan a wider range. The output of the filter is sampled at Nyquist rate and N-point FFT is applied to obtain the frequency domain samples which can be modeled as




  
    
      
        
          
            Y
            m
          
          ⁡
          
            (
            k
            )
          
        
        =
        
          {
          
            
              
                
                  
                    
                      
                        
                          W
                          m
                        
                        ⁡
                        
                          (
                          k
                          )
                        
                      
                    
                    
                      
                        
                          H
                          0
                        
                        ,
                      
                    
                  
                  
                    
                      
                        
                          
                            S
                            m
                          
                          ⁡
                          
                            (
                            k
                            )
                          
                        
                        +
                        
                          
                            W
                            m
                          
                          ⁡
                          
                            (
                            k
                            )
                          
                        
                      
                    
                    
                      
                        
                          H
                          1
                        
                        ,
                      
                    
                  
                
                ⁢
                
                  
                
                ⁢
                k
              
              =
              1
            
            ,
            ...
            ⁢
            
              
            
            ,
            N
          
        
      
    
    
      
        (
        24
        )
      
    
  




where Sm(k) is the transmitted signal by primary users at the output of mth FFT operation, Wm(k) is the white noise sample at kth frequency sample, and N is the FFT size. 0 and 1 represent the null hypothesis and alternate hypothesis respectively. The white noise is modeled as a zero-mean Gaussian random variable with variance σ20, i.e. W(k)=N (0, σ20) The signal term is also modeled as a zero-mean Gaussian variable whose variance is a function of frequency, i.e. S(k)=N (0, σ2k), where σk is the local standard deviation. The variation of σk across frequency depends on the characteristics of primary users signals and operating frequency. Hence, by changing the signal variance across frequency, PSDs for various technologies can be represented.
Frequency Domain Filtering
The magnitude square of FFT output |Ym(k)|2 can be compared with a threshold value for detection of presence of transmission at a particular frequency. Filtering is applied before applying the threshold detector to smooth the FFT output. The optimum filter coefficients depend on the statistics of primary user's signal as well as the noise power. Minimum mean-square error (MMSE) filtering is applied for estimating the noise plus interference ratio for OFDM systems, using filter coefficients derived as a function of the statistics of interference. As an approximation to MMSE filtering, a sliding rectangular window can also be applied for smoothing the spectrum estimates Y(k). Use rectangular filter, the decision statistics at different frequencies can be calculated as




  
    
      
        
          
            
              Y
              ~
            
            m
          
          ⁡
          
            (
            k
            )
          
        
        =
        
          
            ∑
            
              l
              =
              
                k
                -
                
                  ⌊
                  
                    
                      L
                      f
                    
                    /
                    2
                  
                  ⌋
                
              
            
            
              k
              +
              
                ⌊
                
                  
                    L
                    f
                  
                  /
                  2
                
                ⌋
              
            
          
          ⁢
          
            
              ∑
              
                u
                =
                
                  m
                  -
                  
                    L
                    t
                  
                  +
                  1
                
              
              m
            
            ⁢
            
              
                
                
                  
                    Y
                    u
                  
                  ⁡
                  
                    (
                    l
                    )
                  
                
                
              
              2
            
          
        
      
    
    
      
        (
        25
        )
      
    
  




Where [·] represents the floor function, and Lt and Lf are the averaging filter lengths in time and frequency directions respectively.
B. Threshold Detector
The sensing time index m is dropped for simplicity. The output of the rectangular filter {tilde over (Y)}(k) is fed to a threshold device to identify the frequencies occupied by the primary users. This is equivalent to distinguishing between the following two hypotheses:

H0:Y(k)=W(k),  (26)

H1:Y(k)=S(k)+W(k).  (27)

As Y(k) has Gaussian distribution, the decision statistics {tilde over (Y)}(k) follows chi-square distribution with 2LfLt degrees freedom, χ2LfLt. Hence, {tilde over (Y)}(k) can be modeled as:




  
    
      
        
          
            Y
            ~
          
          ⁡
          
            (
            k
            )
          
        
        =
        
          {
          
            
              
                
                  
                    
                      
                        
                          
                            σ
                            0
                            2
                          
                          2
                        
                        ⁢
                        
                          χ
                          
                            2
                            ⁢
                            
                              
                            
                            ⁢
                            
                              L
                              f
                            
                            ⁢
                            
                              L
                              t
                            
                          
                          2
                        
                      
                    
                    
                      
                        
                          H
                          0
                        
                        ,
                      
                    
                  
                  
                    
                      
                        
                          
                            
                              σ
                              0
                              2
                            
                            +
                            
                              σ
                              k
                              2
                            
                          
                          2
                        
                        ⁢
                        
                          χ
                          
                            2
                            ⁢
                            
                              
                            
                            ⁢
                            
                              L
                              f
                            
                            ⁢
                            
                              L
                              t
                            
                          
                          2
                        
                      
                    
                    
                      
                        
                          H
                          1
                        
                        ,
                      
                    
                  
                
                ⁢
                
                  
                
                ⁢
                k
              
              =
              1
            
            ,
            ...
            ⁢
            
              
            
            ,
            N
          
        
      
    
    
      
        (
        28
        )
      
    
  




where σ20 and σ2k are as defined before. The performance of the detection algorithm can be summarized with two probabilities: probability of detection PD and probability of false alarm PF. PD is the probability of detecting a signal on the considered frequency when it is truly present, thus large detection probability is desired. It can be formulated as:

PD=PΓ({tilde over (Y)}>λ|H1),  (29)

where λ is the threshold value. PF is the probability that the test incorrectly decides that the considered frequency is occupied when it actually is not, and it can be written as:

PF=PΓ({tilde over (Y)}>λ|H0).  (30)

PF should be kept as small as possible. The probabilities PF and PD can be calculated as:




  
    
      
        
          
            P
            F
          
          =
          
            1
            -
            
              Γ
              ⁡
              
                (
                
                  
                    
                      L
                      f
                    
                    ⁢
                    
                      L
                      t
                    
                  
                  ,
                  
                    λ
                    
                      σ
                      0
                      2
                    
                  
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        31
        )
      
    
  
  
    
      
        
          
            P
            D
          
          =
          
            1
            -
            
              Γ
              ⁡
              
                (
                
                  
                    
                      L
                      f
                    
                    ⁢
                    
                      L
                      t
                    
                  
                  ,
                  
                    λ
                    
                      
                        σ
                        0
                        2
                      
                      +
                      
                        σ
                        k
                        2
                      
                    
                  
                
                )
              
            
          
        
        ,
      
    
    
      
        (
        32
        )
      
    
  




where Γ(a, x) is the incomplete gamma function provided in (M. Abramowitz and I. Stegun, Eds., Handbook of Mathematical Functions With Formulas, Graphs, and Mathematical Tables, 9th ed. National Bureau of Standards, 1970). In general, PD and PF are inversely proportional. FIG. 7 shows the receiver operating characteristics (ROCS) for a single frequency sample at the output of threshold device for different SNR values. SNR is defined as the ratio of the primary user's signal power to noise power, i.e. SNR=σ2k/σ20. The averaging size is set to 15 in this figure. As this figure clearly shows, the performance of the threshold detector increases at high SNR values. The noise power is normalized before it is fed into the threshold device.
The decision threshold λ can be selected for finding an optimum balance between PD and PF. However, this requires the knowledge of noise and detected signal powers as seen in (31) and (32). Noise power can be estimated, but the signal power is difficult to estimate as it changes depending on the ongoing transmission characteristics and the distance between the cognitive radio and primary user. It is assumed that the noise variance is estimated and the threshold is chosen to obtain a false alarm rate of 0.1 percent.
Feature Extraction
Two features, bandwidth and center frequencies of primary users, are extracted by using the threshold detector output. In order to achieve this, the bins exceeding the threshold value are grouped together. In order to reject the erroneous bins, a method similar to the initial stage of dynamic undersea detection extractor (DUDE) algorithm was applied by defining two parameters Bmin and Gnax. Bmin is the minimum assumed bandwidth for the primary users and Gmax is the maximum gap allowed between two frequency samples. The feature extraction algorithm searches for continuous frequencies which are marked by threshold detector as having signal, not Gmax samples away and with more adjacent frequency samples than Bmin as illustrated in FIG. 8. Hence, using these two parameters, the occupied frequency band can be identified. It is then straightforward to estimate the bandwidth and the center frequency of transmission which make up the feature set x. One drawback of this method, however, is that the cognitive radio may not differentiate between two (or more) superimposed primary user transmissions and threat them as a single transmission with a larger bandwidth.
Classification
Primary users are identified by using the a priori information about their transmission parameters. Feature sets are classified using Bayesian decision rule measurements in the previous section to systems or devices that have the highest a posteriori probability. The classifier can be represented in terms of a set of discriminant functions gi(x), i=1, . . . , K where K is the total number of systems. Note that the same RATs operating at different bands are regarded as different systems for the sake of classification. The transmission band of a system needs to be known for identifying the frequencies occupied by primary users. The discriminant functions are defined as:

gi(x)=log P(x|ωi)+log P(ωi).  (33)

where P(wi) is the prior probability of class wi, and P(x|wi) is the class conditional probability which is assumed to have multivariate normal distribution with mean vector μi and covariance matrix Σi. Under such an assumption, the discriminant functions can be obtained as:




  
    
      
        
          
            g
            i
          
          ⁡
          
            (
            x
            )
          
        
        =
        
          
            
              -
              
                1
                2
              
            
            ⁢
            
              
                (
                
                  x
                  -
                  
                    μ
                    i
                  
                
                )
              
              τ
            
            ⁢
            
              
                ∑
                i
                
                  -
                  1
                
              
              ⁢
              
                
              
              ⁢
              
                (
                
                  x
                  -
                  
                    μ
                    i
                  
                
                )
              
            
          
          -
          
            
              1
              2
            
            ⁢
            log
            ⁢
            
              
              
                ∑
                i
              
              ⁢
              
                
              
              
            
          
          +
          
            log
            ⁢
            
              
            
            ⁢
            
              
                P
                ⁡
                
                  (
                  
                    ω
                    i
                  
                  )
                
              
              .
            
          
        
      
    
    
      
        (
        34
        )
      
    
  




The classification regions defined by (34) are quadratic functions. Therefore, gi(x) is known as normal-based quadratic discriminant function. The classification rule can now be written as: assign a feature vector x to a system wi if gi(x)>gj(x) for all j≠i.
The discriminant function gi(x) is locally estimated. When collaborative sensing is used, classification results can be shared among cognitive radios. Alternatively, cognitive radios can share the discriminant function values, allowing a central unit to combine these values and make a more optimum decision based on the discriminant function values from all cognitive radios present in the environment.
The mean vector μi can be obtained by using the expected values of features, e.g. bandwidth and center frequency in this case. In practice, the covariance matrix Σi is unknown and can be replaced with estimates based on the training data as




  
    
      
        
          
            ∑
            i
          
          ^
        
        ⁢
        
          =
          
            
              1
              
                N
                -
                1
              
            
            ⁢
            
              
                ∑
                
                  j
                  =
                  1
                
                N
              
              ⁢
              
                
                  (
                  
                    
                      x
                      j
                    
                    -
                    
                      μ
                      i
                    
                  
                  )
                
                ⁢
                
                  
                    
                      (
                      
                        
                          x
                          j
                        
                        -
                        
                          μ
                          i
                        
                      
                      )
                    
                    τ
                  
                  .
                
              
            
          
        
      
    
    
      
        (
        35
        )
      
    
  




Furthermore, the covariance matrix can be assumed to be the same for all classes using the same RAT. For example, all the classes using WLAN are expected to have the same covariance matrix as only center frequency is changed. When the estimated features are not correlated to each other, the correlation matrix becomes a diagonal matrix. Different features will have different units and hence proper normalization of this features needs to be established. Moreover, the values of diagonal elements give the weights for each feature, and different weights can be assigned to different features.
The proposed algorithm was tested using two RATs in 2.4 GHz industrial, scientific and medical (ISM) bands, WLAN and Bluetooth signals. The ISM band is an unlicensed band, without primary users. It is assumed that Bluetooth signal is not changing hopping frequency during the detection phase, i.e., the detection process takes much less than 675 μs. Both WLAN and Bluetooth signals are assumed to have the same power and prior probabilities, and no overlap of spectrum between the two signal is allowed. The latter assumption can be justified by the adaptive frequency hopping (AFH) capability of Bluetooth. The spectrum of the ISM band is calculated using 1024-point FFT. Hence the frequency separation between two FFT output samples is 81.5 kHz, i.e. at least 12 frequency points will have Bluetooth signal and at least 246 frequency points will have WLAN signal. The length of the frequency domain filter is set to 5, i.e. Lf=5. FIG. 9 shows the PSD of a WLAN signal operating at channel 5 (Fc=2.432 GHz) and a Bluetooth signal at frequency Fc=2.462 GHz. WLAN and Bluetooth signals are generated by transmitting random data modulated using OFDM and Gaussian frequency shift keying (GFSK) modulations as defined in their corresponding standards.
The features, transmission bandwidth and center frequency, are calculated after passing the FFT output through a threshold device and finding continuous frequencies with power larger than the threshold. The threshold is chosen to have a low false alarm rate of 0.1 percent, i.e. PF=0.001. Note that while setting the threshold, the noise variance is assumed to be known. The feature extraction parameters are set as Gmax=30 and Bmin=6 which correspond to 2.44 MHz and 490 kHz respectively.
FIGS. 10 and 11 show the error rates in detecting the WLAN and Bluetooth signals for additive white Gaussian noise (AWGN) channels. The error rate for identifying a non-existent system (or a system in another band) is shown as well. The time averaging Lt for FIG. 10 and FIG. 11 are set to be 15 and 30 respectively. Hence, sensing durations correspond to 184 μs and 368 μs. The classification is quite successful even at very low SNR scenarios and performance increases with increasing averaging duration. Note that the error rates for WLAN is smaller for a given SNR value because of the large bandwidth and large separation of center frequencies of WLAN signals.
FIG. 12 shows the performance of energy detector based algorithm when the received signal is passed through a multipath channel with exponential delay power delay profile (PDP) with root-mean-squared (RMS) delay spread value of 50 ns. The performance of the Bluetooth signal detection is affected very little while the performance of WLAN detection is decreased. The power fluctuations within the large bandwidth of the WLAN signal due to frequency selectivity and fading causes the performance degradation.
When the received signal is weak (low SNR values), the signal power falls below the noise floor. In such cases, the comparator yields zero output mostly due to the low PF settings. Hence, no transmissions can be detected. As a result of this, the probability of missing ongoing transmissions becomes one and probability of detecting a non-existent signal approaches to zero.
It is also to be understood that the following claims are intended to cover all of the generic and specific features of the invention herein described, and all statements of the scope of the invention which, as a matter of language, might be said to fall therebetween. Now that the invention has been described,

## Claims

1. A method of sensing waveform spectrum, comprising the steps of:
obtaining a radio, further comprising
a cognitive radio engine adapted to examine a radiowave frequency band;
an input in electrical communication with the cognitive radio engine;

extracting signal features from a test frequency, further comprising:
collecting at least one signal from the test frequency, wherein the signal comprises radio access technology information and at least one transmission parameter of the signal;
inputting the at least one signal to the cognitive radio engine;
determining radio access technology from the radio access technology information using a classifier selected from the group consisting of pattern recognition, neural networks, and statistical classification;
extracting the at least one transmission parameter of the signal using at least one algorithm selected from the group consisting of energy detector or cyclostationary feature extraction fast Fourier transform time frequency transform, discrete Fourier transform, and inverse discrete Fourier transform;
wherein a square of the fast Fourier transform output is compared to a threshold value to determine primary use in the test frequency;

using the extracted signal features to identify primary use of a transmission spectrum; and
exploiting the primary use for multi-dimensional spectrum characterization.

2. The method of claim 1, wherein symbol duration and cyclic prefix size are extracted using known OFDM transmission characteristics.

3. The method of claim 1, wherein symbol length is used to calculate cyclic prefix duration and timing offset.

4. The method of claim 1, wherein repetition in cyclic prefix is used to find data duration.

5. The method of claim 4, wherein the cyclic prefix duration is a multiple of symbol duration selected from the group consisting of 1⁄4, 1⁄8, 1/16, and 1/32.

6. The method of claim 1, wherein a band pass filter is used to extract the signal features from the test frequency.

7. The method of claim 6, wherein the filtered signals are sampled at Nyquist rate.

8. The method of claim 1, wherein feature sets are classified using Baysean decisions.

9. The method of claim 1, wherein the discrete Fourier transform uses algorithms selected from the group consisting of multiple signal classification, akaike's criteria, and eigene decomposition.

10. The method of claim 1, wherein the transmission parameters are selected from the group consisting of bandwidth, center frequency, single carrier, multi-carrier, maximum signal duration, multiple accessing, duplexing method, frame duration, spreading codes, cumulants, power moment matrices, and zero-crossing rate.

11. The method of claim 1, wherein the at least one algorithm is energy detector and cyclostationary feature extraction.

12. The method of claim 1, further comprising identifying active subcarriers, wherein the active subcarriers are identified by analyzing the power level at the output of FFT.