Next Article in Journal
Applying High-Speed Vision Sensing to an Industrial Robot for High-Performance Position Regulation under Uncertainties
Previous Article in Journal
Eigenspace-Based Generalized Sidelobe Canceler Beamforming Applied to Medical Ultrasound Imaging
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exact Distributions of Finite Random Matrices and Their Applications to Spectrum Sensing

by
Wensheng Zhang
1,*,
Cheng-Xiang Wang
2,
Xiaofeng Tao
3 and
Piya Patcharamaneepakorn
2
1
School of Information Science and Engineering, Shandong University, Jinan 250100, China
2
Institute of Sensors, Signals and Systems, School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh, EH14 4AS, UK
3
Wireless Technology Innovation Institute, Beijing University of Posts and Telecommunications, Beijing 100876, China
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(8), 1183; https://doi.org/10.3390/s16081183
Submission received: 12 June 2016 / Revised: 22 July 2016 / Accepted: 25 July 2016 / Published: 29 July 2016
(This article belongs to the Section Physical Sensors)

Abstract

:
The exact and simple distributions of finite random matrix theory (FRMT) are critically important for cognitive radio networks (CRNs). In this paper, we unify some existing distributions of the FRMT with the proposed coefficient matrices (vectors) and represent the distributions with the coefficient-based formulations. A coefficient reuse mechanism is studied, i.e., the same coefficient matrices (vectors) can be exploited to formulate different distributions. For instance, the same coefficient matrices can be used by the largest eigenvalue (LE) and the scaled largest eigenvalue (SLE); the same coefficient vectors can be used by the smallest eigenvalue (SE) and the Demmel condition number (DCN). A new and simple cumulative distribution function (CDF) of the DCN is also deduced. In particular, the dimension boundary between the infinite random matrix theory (IRMT) and the FRMT is initially defined. The dimension boundary provides a theoretical way to divide random matrices into infinite random matrices and finite random matrices. The FRMT-based spectrum sensing (SS) schemes are studied for CRNs. The SLE-based scheme can be considered as an asymptotically-optimal SS scheme when the dimension K is larger than two. Moreover, the standard condition number (SCN)-based scheme achieves the same sensing performance as the SLE-based scheme for dual covariance matrix K = 2 . The simulation results verify that the coefficient-based distributions can fit the empirical results very well, and the FRMT-based schemes outperform the IRMT-based schemes and the conventional SS schemes.

1. Introduction

Cognitive radio networks (CRNs), working as a sharp tool to deal with the spectrum scarcity problem, have been recognized as one of the most promising communication technologies in recent years [1,2,3]. As the key function of CRNs, spectrum sensing (SS) is used to check the accessibility of the target frequency bands, which belong to the primary user (PU) systems, but are not occupied [4,5]. The SS schemes for CRNs have been comprehensively summarized in [2,4,5]. Moreover, the random matrix theory (RMT) has also been used in wireless communication systems [6,7,8,9,10,11].
According to the size of matrix dimensions, the RMT can generally be divided into two types: infinite RMT (IRMT) with large dimensions and finite RMT (FRMT) with finite dimensions. The characteristics of the IRMT have been discussed by some mathematicians, e.g., Alan Edelman [12], Jinho Baik and Jack W. Silverstein [13] and Zhidong Bai [14]. The distributions of the eigenvalues and the condition number (CN) in the IRMT have been expressed as asymptotic formulations, such as the Marchenko–Pastur distributions [15] and the Tracy–Widom distributions [16]. The asymptotic results have been used to construct the SS systems for the CRNs. Based on the eigenvalue distributions of sample covariance matrices, some blind SS schemes have been discussed in [17,18,19]. The eigenvalue ratio in the RMT is a sharp tool to build blind SS schemes.
Compared to the IRMT with large dimensions, the FRMT with finite dimensions has the following advantages in the applications of the practical wireless communication systems. First, the distributions of the FRMT can be evaluated and formulated in the exact and tractable expressions. Second, the FRMT is more suitable for the practical wireless communication systems due to their finite dimensions. However, to the best of our knowledge, the dimension boundary between the IRMT and the FRMT has not yet been precisely defined. We provide the definition of the dimension boundary for Wishart matrices by comparing the theoretical eigenvalue distribution with the empirical eigenvalue distribution.
The FRMT distributions have been discussed in [20,21,22,23,24,25], in which the extreme eigenvalues (largest and smallest) distributions and scaled largest eigenvalue (SLE) distributions were given in [20,21], and [22,23,24], respectively. Our previous work in [25] also provided the extreme eigenvalue distributions based on coefficient matrices. We extend that work and derive the unified distributions of the SLE and the Demmel condition number (DCN) with proposed coefficient matrices (vectors). Note that the FRMT distributions have been given in different expressions, and some of them were not tractable due to high complexity. The unified and compact expressions of the FRMT distributions are required for the applications in the CRNs, because such expressions provide a fast and convenient way to calculate the exact thresholds.
The existing SS schemes can be generally divided into two categories. One is the single SS schemes, where energy detection, feature matched filter detection and cyclostationarity feature detection are usually exploited by the distributed sensors [4,26]; the other is the cooperative SS schemes, where multiple sensors send their sensing samples or sensing results to the fusion centre (FC), and the final sensing result can be determined [27,28]. The cooperative SS schemes can achieve superior sensing performance compared to the single SS schemes. Recently, a robust SS scheme based on crowd sensors, a cooperative SS scheme via Dirichlet process and a multi-agent-based SS framework for 5G networks have been discussed in [29,30,31], respectively. In particular, the random matrix theory (RMT) has been used to construct SS systems that can achieve better sensing performance compared to conventional SS schemes [32,33,34,35,36]. However, the RMT-based SS schemes have to compute the eigenvalues of the sample matrices, leading to extra computational complexity, especially for the IRMT-based schemes.
The FRMT-based SS systems in the CRNs have been discussed in [36], in which the exact SCN distributions of the dual Wishart matrix given in [9] was exploited. The results of [36] indicated that the FRMT-based SS systems can achieve superior sensing performance compared to the counterparts based on the IRMT due to the exact distributions calculated with the FRMT. However, the results of [9] were limited for the dual Wishart matrix, whose dimensions are 2 × 2 . The exact distributions for the SCN of the Wishart matrix with arbitrary dimensions is still unknown to us, although the sensing performance may be evaluated in an empirical way. Moreover, the sensing performances of other FRMT-based SS systems (e.g., the SLE-based scheme, the largest eigenvalue (LE)-based scheme and the DCN-based scheme) are not clear. This paper aims to address the above problems, and the key contributions are summarized as follows:
  • The exact distributions of the FRMT are unified with the proposed coefficient matrices (vectors). The coefficient reuse mechanism is studied, i.e., the LE distributions and the SLE distributions can be formulated with the same coefficient matrices. Moreover, the SE distributions and the DCN distributions share the same coefficient vectors. In particular, a new and simple CDF of the DCN is formulated with the coefficient vector. The dimension boundary between the IRMT and the FRMT is defined by evaluating the theoretical and empirical eigenvalue distributions.
  • The sensing performances of the FRMT-based SS schemes are analysed and evaluated. The asymptotical optimal SS schemes in the FRMT with varying dimensions are proposed. We demonstrate that the SLE-based scheme is asymptotically optimal when K > 2 , and the SCN-based scheme possesses identical sensing performance with the SLE-based scheme when K = 2 .
The remainder of the paper is organized as follows: Section 2 presents the system model. The IRMT and FRMT are discussed in Section 3 and Section 4, respectively. The SS systems based on the FRMT are provided in Section 5. The numerical results and analysis are provided in Section 6. This paper is concluded in Section 7.

2. System Model

In this section, we introduce the cooperative SS model based on the RMT, followed by some key definitions of the RMT. The system model is shown in Figure 1, in which the secondary user (SU) system working in the interweave CR paradigm [37] dynamically accesses the PU’s spectrum when the PU signals are absent in the target frequency bands. The SUs periodically sense PU signals and send PU samples to the FC. The FC in the SU system periodically gathers the samples from the distributed SUs and makes the final decisions of the availability of the spectrum. Let K and N denote the number of SUs and the number of the samples of PU signals per SU, respectively. At the FC, a K × K covariance matrix Y can be generated by Y = X · X , where X denotes the K × N sample matrix and denotes transpose-conjugate (Hermitian). Here, we assume that K N , and X is full row-rank.
Let λ 1 λ 2 λ K denote K ordered eigenvalues of Y . The standard condition number (SCN) denoted by ξ is defined as the ratio of the largest eigenvalue λ K to the smallest eigenvalue λ 1 , i.e.,
ξ λ K λ 1
The Demmel condition number (DCN) (the DCN in [38] was defined as the square-root of κ given in Equation (2)) denoted by κ is defined as the ratio of the trace to the smallest eigenvalue:
κ R λ 1
where R = Tr ( Y ) = i = 1 K λ i is the trace of Y The scaled largest eigenvalue (SLE) denoted by ψ is defined as the ratio of the largest eigenvalue to the trace, i.e.,
ψ λ K R
For the IRMT, the critical requirement is that the matrix dimensions should be very large, that is ( K , N ) , leading to the asymptotic distributions. Hence, we cannot get the exact distributions of the IRMT due to the large dimensions. The dimension requirements of the IRMT limit its further applications in wireless communication systems based on the following reasons. First, the asymptotic results cannot provide an exact evaluation to the practical systems. Second, the sizes of the practical systems cannot be large enough to satisfy the dimension requirements of the IRMT. Therefore, the FRMT with finite dimensions is more useful for the practical systems because both the number of sensors K and the number of the samples per sensor N cannot be very large for the practical CRNs. Moreover, the exact and closed-form expressions of the FRMT are required to calculate the exact thresholds in the practical CRNs.
The i-th PU signal sample gathered by the k-th sensor is:
x k i = h k i · s k i + n k i , ( k = 1 , 2 , , K , i = 1 , 2 , , N )
where h k i denotes the channel gain between the k-th sensor and the FC, the noise n k i denotes an independent and identical distributed (i.i.d.) circular complex Gaussian process with zero mean and variance σ n 2 , the PU signal s k i represents an i.i.d. circular complex Gaussian process with zero mean and variance σ s 2 and the signal-to-noise ratio (SNR) is defined as ρ = σ s 2 σ n 2 . Therefore, at the FC, a binary hypothesis test (HT) can be formulated as:
(5a) H 0 : X = X n (5b) H 1 : X = X s + X n
where H 0 and H 1 denote the absence and the presence of the PU signal, respectively, and the k-th row of the K × N sample matrix X is the N-length sample vector from the k-th sensor. The covariance matrix Y can be regarded as a Wishart matrix under the hypothesis H 0 .
The sensing performance is evaluated by the probability of detection ( P D ) under a given probability of false alarm ( P F ) denoted by δ. P D and P F can be formulated as:
P D = Pr { T γ | H 1 }
P F = Pr { T γ | H 0 }
where T denotes the proposed sensor and γ is the given threshold; Pr { · } denotes the probability operation. The threshold γ can be calculated with:
γ = P T - 1 ( 1 - δ )
where P T - 1 ( · ) denotes the inverse function of the cumulative distribution function (CDF) of the sensor T. The sensing performance is mainly affected by the threshold γ, which is calculated with the distributions of the sensor T. This is also the reason why the SS systems based on the FRMT can achieve superior sensing performance compared to those of the IRMT [36]. Based on the RMT, the sensor T is constructed by the characteristics of the covariance matrix Y . That is, it is possible to build the sensor with the SCN ξ, the DCN κ, the SLE ψ, and so on. We evaluate the sensors based on various parameters and try to obtain an optimal SS scheme in the FRMT paradigm. For example, if the SCN ξ is used to construct the sensor T, the formulation is defined as the ratio of the largest eigenvalue λ K to the smallest eigenvalue λ 1 :
T λ K λ 1 = ξ
For the FRMT, the distributions of the characteristics are critical in the construction of the SS systems because the corresponding thresholds are generated from the CDF of the sensors shown in Equation (8). We unify the distributions with the proposed coefficient matrices (vectors) and provide the simple and compact formulations. The corresponding inverse functions can be generated by the given closed-form expressions.

3. Asymptotic Distributions in the IRMT

Before the discussions of the FRMT, we investigate the distributions of the IRMT in this section for comparisons. These results are used to indicate the characteristics of the IRMT when the dimensions are very large, that is K , N . The dimension boundary between the IRMT and the FRMT is initially defined here to asymptotically divide the IRMT with large dimensions and the FRMT with finite dimensions.

3.1. Eigenvalue Distributions in the IRMT

3.1.1. General Eigenvalue Distributions in the IRMT

Lemma 1 (Marchenko-Pastur Distributions). 
For a central complex Wishart matrix W with the large dimensions ( K , N , K N = c ), the probability density function (PDF) [15] and the CDF [36] of eigenvalue λ can be formulated as:
f λ ( x ) = ( x - α ) ( β - x ) 2 π c x
F λ ( x ) = 1 2 + f λ ( x ) + ( 1 - c ) 2 π asin ( 1 + c ) x - ( 1 - c ) 2 2 x c + ( 1 + c ) 2 π asin 1 + c - x 2 c
where the upper bound and lower bound of λ are β ( 1 + c ) 2 and α ( 1 - c ) 2 , respectively.

3.1.2. Extreme (Largest or Smallest) Eigenvalue Distributions in the IRMT

Lemma 2 (Tracy-Widom Distributions). 
For a central complex Wishart matrix W with the large dimensions ( K , N , K N = c ), the PDF and CDF of the centralized and normalized extreme eigenvalue λ ¯ E can be formulated as:
f λ ¯ E ( x ) = d F TW 2 ( x ) d x
F λ ¯ E ( x ) = F TW 2 = exp - r ( x - r ) q 2 ( x ) d x
where q 2 ( x ) is the Hastings–McLeod solution of the Painlevé equation of Type II [16,39]. The centralized and normalized extreme eigenvalues are defined as:
λ ¯ E K λ K - β μ
λ ¯ E 1 λ 1 - α ν
where β and α are defined in Lemma 1; the normalized factors μ and ν are defined as:
μ ( K + N ) ( 1 / K - 1 / N ) 1 / 3
ν ( K - N ) ( 1 / K - 1 / N ) 1 / 3

3.2. Standard Condition Number Distributions in the IRMT

Lemma 3 (Tracy-Widom-Curtiss Distributions). 
For a central complex Wishart matrix W with the large dimensions ( K , N , K N = c ), the PDF [36,40] and the CDF [36] of ξ can be formulated as:
f ξ ( x ) = - 1 μ · ν 0 r · f λ ¯ E x · r - β ν · f λ ¯ E r - α μ d r
F ξ ( x ) = α μ f λ ¯ E ( - r ) · F λ ¯ E x · ( β - μ · r ) - α ν d r
where f λ ¯ E ( x ) and F λ ¯ E ( x ) are defined in Lemma 2.
Seeing the distributions of the IRMT from Lemma 1 to Lemma 3, the expressions are unclosed. This means that only the asymptotic formulations can be exploited to determine the thresholds, leading to low sensing performances.

3.3. Dimension Boundary between the IRMT and the FRMT

The dimension boundary is provided to asymptotically divide the IRMT and the FRMT. If the dimensions of a specific random matrix are larger than the dimension boundary, such a matrix can be considered as an infinite random matrix, and the IRMT is therefore applicable. Otherwise, the random matrix with the dimensions less than the dimension boundary is finite, and the FRMT can be used to analyse such a matrix.
Theorem 1 (Dimension Boundary). 
For a central complex Wishart matrix W with the dimensions K and N, the dimension boundary between the IRMT and the FRMT can be defined as:
( K b , N b ) = argmin λ e ( K , N ; λ ) - t ( K , N ; λ ) s.t. 0 < K N < , K / N = c , λ [ α , β ]
where ( K b , N b ) denotes the dimension boundary, e ( K , N ; λ ) and t ( K , N ; λ ) are the empirical and theoretical eigenvalue distributions, respectively, c 0 , 1 is a fixed value and all eigenvalues lie in the eigenvalue interval [ α , β ] .
Proof. 
Under the condition that K , N and K / N = c , the empirical eigenvalue distributions of the Wishart matrix can be precisely evaluated by the corresponding theoretical eigenvalue distributions. Based on the convergency of the theoretical eigenvalue distribution, the dimension boundary ( K b , N b ) can always be achieved. When the difference between the two distributions goes to the minimum, the dimension boundary is determined. ☐
Note that Theorem 1 is provided from the view of the IRMT, i.e., when the dimensions of the Wishart matrix are so large that the theoretical eigenvalue distributions of the IRMT can precisely fit the corresponding empirical eigenvalue distributions. Moreover, this theorem can also be defined from the view of the FRMT, and in that case, both the accuracy and the computational complexity of the theoretical eigenvalue distributions in the FRMT should be considered.
However, the dimension boundary determined by Theorem 1 only guarantees that the Wishart matrix with the dimensions larger than such a boundary can be analysed with the IRMT. In practice, there should be a large margin between the IRMT and the FRMT and the Wishart matrices with the dimensions lying in such a margin being able to be handled with both theories. The following corollary provides a scheme to numerically determine the dimension boundary using Theorem 1 and Marchenko–Pastur distributions in Lemma 1.
Corollary 1 (Dimension Boundary Decision Scheme). 
For an eigenvalue of the Wishart matrix with the dimensions K and N, the distance between the empirical and theoretical eigenvalue distributions is defined as:
d ( K , N ; λ ) = f e ( K , N ; λ ) - f t ( K , N ; λ )
where the theoretical eigenvalue PDF f t ( K , N ; λ ) can be implemented by the Marchenko-Pastur PDF in Equation (10), and the empirical eigenvalue PDF can be generated in a numerical way. For all of the eigenvalues lying in the eigenvalue interval [ α , β ] , the sum of the eigenvalue distance can be calculated with:
S ( K , N ) = λ [ α , β ] d ( K , N ; λ )
The dimension boundary ( K b , N b ) can be achieved if the following condition is satisfied:
S ( K b + Δ , N b + Δ / c ) S ( K b , N b ) ϕ
where Δ denotes the dimension step and ϕ is a given dimension threshold beyond which the corresponding Wishart matrix is considered as infinite.
Proof. 
For the theoretical eigenvalue distributions, such as Marchenko-Pastur distributions or Tracy-Widom distributions, when the dimensions get larger, the theoretical eigenvalue distributions converge to the empirical eigenvalue distributions. For the eigenvalues lying in the eigenvalue interval, the sum of the distance between two distributions can be larger than the given dimension threshold for the fixed dimension step. ☐
Based on Theorem 1 and Corollary 1, the dimension boundary between the IRMT and the FRMT can be determined. As shown in Figure 2, an example is provided to illustrate the boundary decision process. In the example, the dimensions are set to K = 2 : 6 : 26 and c = 0.1 , 0.2 , 0.5 . When the dimension step is Δ = 6 and the the dimension threshold is ϕ = 0.7 , the dimension boundaries for three cases are ( 20 , 200 ) , ( 20 , 100 ) and ( 20 , 40 ) . Note that for different Δ and ϕ , the corresponding boundary should be different.

4. Exact Distributions in the FRMT

In this section, we summarize the current research of the FRMT and unify the results with the proposed coefficient matrices (vectors). Based on the coefficient mechanism, the closed-form distributions of the FRMT are derived in exact and compact expressions.

4.1. Eigenvalue Distributions in the FRMT

Lemma 4 (James-Edelman Distribution). 
The joint PDF of K ordered eigenvalues of W is provided by [20,25,41]:
f Λ ( λ 1 , . . . , λ K ) = i = 1 K λ i N - K exp - λ i ( K - i ) ! ( N - i ) ! 1 i < j K λ i - λ j 2
where Λ denotes the eigenvalue set { λ 1 , . . . , λ K } .
The corresponding joint CDF of K ordered eigenvalues of W can be straightforwardly calculated with the definite integral of f Λ ( λ 1 , . . . , λ K ) :
F Λ ( λ 1 , . . . , λ K ) = 0 λ 1 0 λ K f Δ ( x 1 , , x K ) d x 1 d x K
However, to the best of our knowledge, there is still no closed-form solution of such a CDF.

4.2. Extreme (Largest or Smallest) Eigenvalue Distributions in the FRMT

Theorem 2 (Largest Eigenvalue Distributions). 
For a central complex Wishart matrix W with finite dimensions, the PDF of the largest eigenvalue λ K can be formulated in a closed-form formulation [20,25]:
f λ K ( x ) = k = 1 K exp - k x n = N - K k ( N + K - 2 k ) P k - 1 , n - N + K ( K , N ) x n
where P k , n ( K , N ) is the ( k , n ) -th entry of the coefficient matrix P ( K , N ) .
The corresponding CDF in a closed-form formulation can be given as [20,25]:
F λ K ( x ) = k = 0 K exp ( - k x ) n = 0 k ( N + K - 2 k ) C k , n ( K , N ) x n
where C k , n ( K , N ) is the ( k , n ) -th entry of a coefficient matrix C ( K , N ) associated with P ( K , N ) .
Proof. 
Based on the joint PDF of ordered K eigenvalues in Equation (24), the CDF of the largest eigenvalue λ K can be obtained by integration, that is,
F λ K ( x ) = 0 x 0 x i = 1 K λ i N - K exp - λ i det Ω K ! ( K - i ) ! ( N - i ) ! d λ 1 d λ K
where Ω is the square Vandermonde matrix of { λ 1 , . . . , λ K } , and its determinant can be calculated with ([20] Appendix A):
Sensors 16 01183 i001
Based on the determinant in Equation (29), the corresponding coefficient matrices P ( K , N ) for PDF and C ( K , N ) for CDF can be calculated for the specific K and N. More details about the generation algorithms for the coefficient matrices can be found in [25], and some typical coefficient matrices and heuristic examples are provided in Appendix A. More coefficient matrices can be calculated with the generation algorithms provided in [25]. ☐
Theorem 3 (Smallest Eigenvalue Distributions). 
For a central complex Wishart matrix W with finite dimensions, the PDF of the smallest eigenvalue λ 1 can be formulated in a closed-form formulation [21,25]:
f λ 1 ( x ) = - exp ( - K x ) m = N - K N K - K 2 p m - ( N - K - 1 ) ( K , N ) x m
where p m - ( N - K - 1 ) ( K , N ) is the m - N - K - 1 -th entry of the coefficient vector p ( K , N ) .
The corresponding CDF in a closed-form formulation can be given as [21,25]:
F λ 1 ( x ) = 1 - exp - K x m = 1 N K - K 2 + 1 c m ( K , N ) x m - 1
where c m ( K , N ) is the m-th entry of the coefficient vector c ( K , N ) . Different from the coefficient matrices P [ K , N ] and C [ K , N ] , the coefficients p [ K , N ] in Equation (30) and c [ K , N ] in Equation (31) are coefficient vectors.
Proof. 
The proofs for the PDF and CDF can be found in [21,25], and the coefficient vectors are generated in [25]. Some typical coefficient vectors and the corresponding examples are provided in Appendix B. ☐

4.3. Condition Number Distributions in the FRMT

4.3.1. Standard Condition Number

Lemma 5 (SCN Distributions). 
Let ξ denote the SCN of the central complex Wishart matrix W with finite dimensions. The PDF and the CDF of ξ can be formulated as: [9,36]
f ξ ( x ) = Φ ( N ) [ A ( x ; N , N - 1 ) + A ( x ; N - 2 , N + 1 ) - 2 A ( x ; N - 1 , N ) ]
F ξ ( x ) = Φ ( N ) [ B ( x ; N , N - 1 ) + B ( x ; N - 2 , N + 1 ) - 2 B ( x ; N - 1 , N ) ]
where:
Φ ( N ) = 1 2 ( N - 1 ) ! ( N - 2 ) !
A x ; M , N = ( N - 1 ) ! ( x + 1 ) 2 n = 0 N - 1 ( M + 1 - n x ) x n ( x + 1 ) n m = 1 M n + m x + 1
B ( x ; M , N ) = ( N - 1 ) ! M ! - 1 x + 1 n = 0 N - 1 x n ( x + 1 ) n m = 1 M n + m x + 1
Note that the above distributions are limited to the dual Wishart matrix W ( K = 2 ).

4.3.2. Demmel Condition Number

Theorem 4 (DCN Distributions). 
Let κ denote the DCN of the central complex Wishart matrix W with finite dimensions. The PDF of κ can be formulated as [42]:
f κ ( x ) = Γ ( m 2 i = N K ( N K ) K P i ( N K 1 ( K , N ) x m 2 Γ ( m 2 i 1 ) ( x K ) 2 + i m 2
where m = 2 K N and p ( K , N ) is the coefficient vector indicated in Theorem 3; Γ ( n ) = n ! denotes the Gamma function.
We here provide a new and simple CDF of κ in a closed-form solution using the hypergeometric function instead of the expression in ([42] Equation (18)), which is not tractable due to high complexity:
F κ ( x ) = 1 - 2 Γ m 2 x 1 - m 2 ( m - 2 ) i = N - K ( N - K ) K ( - 1 ) i - m 2 c i - ( N - K - 1 ) ( K , N ) F 1 2 2 - m 2 + i , 1 - m 2 ; 2 - m 2 ; x K Γ m 2 - i - 1 K 2 + i - m 2
Proof. 
The CDF of κ is the integral of the PDF in Equation (36), i.e.,
F κ ( x ) 0 x f κ ( r ) d r
Considering the key part of Equation (37) r - m 2 r - K i + 2 - m 2 , let - m 2 = μ - 1 , i + 2 + m 2 = ν and β = - 1 K
Based on the integral of ([43] Equation (3.194.1)):
0 x r μ - 1 1 + β r ν d r = x μ μ F 1 2 ν , μ ; 1 + μ ; - β x ,
the key part of Equation (38) can be calculated as:
0 x r - m 2 r - K i + 2 - m 2 d r = x 1 - m 2 1 - m 2 - K 2 - m 2 + i F 1 2 2 - m 2 + i , 1 - m 2 ; 2 - m 2 ; x K
Including the corresponding entry of the coefficient vector c i - ( N - K - 1 ) ( K , N ) and other coefficients of Equation (37), such as Γ m 2 and Γ m 2 - i - 1 , the final CDF of Equation (38) can be calculated after some algebra. Note that the coefficient vector c ( K , N ) is for the distribution of the smallest eigenvalue λ 1 ; the final CDF of the DCN κ can still use the corresponding coefficients accordingly. ☐

4.3.3. Scaled Largest Eigenvalue

The SLE distributions of the complex and central Wishart matrix W with any dimensions have been discussed in [23,24,44]. The centralized and normalized SLE asymptotically follows the Tracy–Widom distribution in the paradigm that the matrix dimensions should go to infinity [44]. The exact and closed-from distributions of the SLE have been independently given in [23] and [24]. However, the coefficients in the expressions of [24] were calculated with the methods given in [45], which were complicated and not tractable. We unify the results in [23] with the coefficient matrices C [ K , N ] and P [ K , N ] given in Theorem 1.
Theorem 5 (SLE Distributions). 
Let ψ denote the SLE of the central complex Wishart matrix W with finite dimensions. The PDF of ψ can be formulated as [23]:
f ψ ( x ) = m 2 - 1 ! K m 2 - 1 i = 1 K j = N - K ( N + K ) i - 2 i 2 ( K - i x ) m 2 - j - 2 ( m 2 - j - 2 ) ! P i , j - ( N - K - 1 ) ( K , N ) x j U 1 - i x K
where m = 2 K N , P ( K , N ) is the coefficient matrix and U x = 1 , x 0 0 , x < 0 denotes the step function. The corresponding CDF is formulated as:
F ψ ( x ) = ( m 2 - 1 ) ! K m 2 - 1 i = 1 K j = N - K ( N + K ) i - 2 i 2 i m 2 - j - 2 C i , j - ( N - K - 1 ) K , N Φ x
where C ( K , N ) is the coefficient matrix. The function Φ ( x ) is:
Φ ( x ) = B x U K i - x + B K i U x - K i - B ( 1 )
and the function B ( x ) is defined as:
B ( x ) = K i m 2 - j - 2 q = 0 m 2 - j - 2 - i K q x q + j + 1 m 2 - j - 2 - q ! ( q ) ! ( j + q + 1 )
Proof. 
The proof and the deductions can be found in [23], and the coefficient matrices are shown in Theorem 2. ☐

5. SS Schemes Based on the FRMT

In this section, the FRMT-based SS schemes are proposed for the CRNs. In the paradigm of the FRMT, the SLE-based scheme can be considered as an asymptotic generalized likelihood ratio test (GLRT) under the condition that the hypothesis H 0 is clear and the hypothesis H 1 is unknown.

5.1. SS Schemes Based on Asymptotic GLRT

In a sensing period, the covariance matrix Y can be generated from the sample X , that is Y = X X . The likelihood ratio (LR) of Y can be written as:
L = f Y ; ρ , σ n 2 | H 1 f ( Y ; σ n 2 | H 0 )
where the channel effect and the PU signal variance σ s 2 can be indicated by SNR ρ. Given a threshold γ, the corresponding GLRT can be expressed as:
L T = sup ρ , σ n 2 f ( Y ; ρ , σ n 2 | H 1 ) sup σ n 2 f ( Y ; σ n 2 | H 0 ) > H 1 < H 0 γ
where the symbol sup denotes the supreme. It is reasonable to assume that the noise variance σ n 2 can be estimated in advance for the practical CRNs, and its value is assumed to be unitary, that is σ n 2 = 1 . In particular, the PU signal SNR ρ is still unknown, and hence, the distributions of Y in the hypothesis H 1 cannot be obtained. The characteristics of the covariance Y can be indicated by its corresponding eigenvalues.
Loosing the requirements of the GLRT by replacing the distributions of Y with its eigenvalues, we can further revise the GLRT as:
L T = sup ρ , σ n 2 = 1 f ( λ 1 , , λ K ; ρ | H 1 ) sup σ n 2 = 1 f ( λ 1 , , λ K | H 0 ) > H 1 < H 0 γ
where the noise variance is assumed to be unitary. Seeing the right side of the formulation, the eigenvalue PDFs under H 0 can be obtained from the results given in Section 3 and Section 4. The distributions of the eigenvalues under H 1 can be asymptotically described by the spiked population model [13,46], in which the large eigenvalues liking the ‘spikes’ can be separated from the bulk of the eigenvalues. However, the number and amplitudes of the ‘spikes’ are hard to estimate. How to construct the ‘right’ sensor that can utilize more information of the covariance matrix Y is the key issue. In other words, how to use the eigenvalues to approximately denote the supremum of the PDF of all eigenvalues under H 1 , which is a tough problem, needs to be addressed.
The GLRT can be considered as the optimal test under the conditions that the PDFs of Y under H 0 and H 1 are clear and all of the parameters such as σ n 2 and ρ are exactly estimated [47]. The conditions are too strict to be satisfied; especially, we cannot obtain the exact distributions of Y under H 1 . In order to loose the rigid requirements and use the covariance matrix, the characteristics of the covariance matrix Y have been used to construct the SS schemes, which can asymptotically follow the GLRT in the IRMT paradigm [32,33,34,35]. In the IRMT, the largest eigenvalue λ K [32], the SCN ξ [33,34] and the SLE κ [35] have been exploited to construct the SS systems, which achieved relatively better sensing performance. Especially, the SLE-based scheme can be considered as the asymptotic GLR statistic ([35] Equation (5)). However, the SLE in the paradigm of IRMT is not suitable for the SS schemes because the SLE distributions in the IRMT cannot be determined.

5.2. FRMT-Based SS Schemes

In the FRMT, we have checked the sensing performance of the SCN-based SS scheme and found that the proposed scheme outperforms all of the IRMT-based schemes [36] due to the exact SCN distributions. We use the the exact distributions of the SLE and the SCN to construct the SS systems, and the corresponding sensing performances are evaluated.
The SS algorithm (taking the SCN-based scheme as an example) based on the FRMT can be described as follows:
  • Given the P F ( δ = 0.01 : 0.01 : 0.99 ), based on the CDF of ξ in Equation (33), the corresponding thresholds are generated by γ = F ξ - 1 ( 1 - δ ) .
  • Construct the sample matrix X = [ x k n ] k = { 1 , , K } , n = { 1 , , N } with K × N PU samples from K sensors, each gathering N samples.
  • Construct the K × K covariance matrix Y = X · X ; calculate K ordered eigenvalues λ 1 λ 2 λ K ; and let the SCN ξ denote the sensor T, that is T ξ = λ K λ 1 .
  • Compare T to the required threshold γ, and record the result C 1 = C 1 + 1 if ξ γ ; otherwise, C 0 = C 0 + 1 .
  • Repeat the above operations C times, C = C 0 + C 1 , and evaluate the sensing performance by P D = C 1 C .
The threshold is the key issue in the construction of the SS scheme based on the RMT. The distributions of the RMT characteristics are required to generate the corresponding threshold. We summarize the distributions of the IRMT and FRMT in Section 3 and Section 4. The results are indicated in the following table. For a given δ, the threshold is generated by γ = F T - 1 ( 1 - δ ) , where F T - 1 ( x ) is the inverse function of the CDF of T. If the CDF of T is not exact or too complicated, see the CDF of λ K in Theorem 3 (Tracy–Widom distributions), we cannot calculate the exact thresholds. Therefore, the sensing performance based on IRMT is not satisfied for the CRNs. However, in the FRMT paradigm, the thresholds can be exactly calculated based on the exact distributions of the sensor T; see Table 1. Based on the exact distributions of the FRMT, the corresponding thresholds can be calculated analytically. The eigenvalues are used to constructed the GLRT in Equation (48).
Compared to the SCN-based scheme and the DCN-based scheme, the SLE-based scheme can be considered as an asymptotically optimal SS scheme in the FRMT paradigm based on the following facts:
  • The SLE can be considered as the generalized likelihood ratio (GLR), leading to the fact that the SLE-based scheme is asymptotically GLRT [24,35], which can be considered as an asymptotically optimal SS scheme for the FRMT-based SS schemes;
  • Compared to the SS schemes based on the SCN or the DCN, all eigenvalue information under the hypothesis H 1 are included;
  • The compact and closed-form distributions of the SLE in the FRMT are available, and the exact thresholds can be generated.

6. Numerical Results and Analysis

In this section, the theoretical results, especially the distributions in the IRMT and FRMT, are verified. The sensing performances of the FRMT-based SS schemes are evaluated by the simulations. Note that only the PDFs are indicated to verify the proposed results. The CDFs can be theoretically or numerically generated by the corresponding PDFs, and the verifications of the CDFs are omitted for simplicity.

6.1. Theoretical Results’ Verifications

The main distributions in the IRMT and FRMT are verified in this subsection.

6.1.1. IRMT Verifications

The theoretical distributions in the IRMT are verified with the empirical distributions. In order to indicate the characteristics of the IRMT, the dimensions of the random matrices should be very large, that is, K , N . In Theorem 1 and Corollary 1, the dimension boundary between the IRMT and the FRMT is determined. For Wishart matrices with c = 0.1 , when the dimensions are larger than the boundary ( 20 , 200 ) , the theoretical results in the IRMT can work very well. In the following simulations, we set K = 30 and N = 300 .
The asymptotic results of the Marchenko–Pastur distributions and the Tracy–Widom distributions are indicated in Figure 3 and Figure 4, respectively. The results indicate that the theoretical distributions can fit the empirical distributions when the dimensions are very large. Note that the results in ([36], Figure 3) indicated that the Marchenko–Pastur distributions in Equation (10) cannot work well when the dimensions are not very large ( K = 5 , N = 25 ).

6.1.2. FRMT Verifications

The distributions in the FRMT, including extreme (largest and smallest) eigenvalues, the SCN, the SLE and the DCN, are verified in the following simulations. Figure 5 indicates the distributions of the extreme (largest and smallest) eigenvalues of the central and complex Wishart matrix with finite dimensions ( K = 3 , 4 , N = 4 , 5 ). The theoretical results match the empirical distributions very well, even if the dimensions are finite. The simulations verify the exact distributions in Theorem 2 and Theorem 3 and indicate that the coefficient matrices in the closed-form solutions can work well for the FRMT.
The DCN distributions provided in Theorem 4 are verified in Figure 6, in which the theoretical results match the empirical results very well. The simulations also indicate that the PDF expressions in Equation (37) based on the coefficient vector p ( K , N ) can precisely describe the DCN PDFs for the finite Wishart matrix. The PDFs of the SLE in Theorem 5 based on the coefficient matrix P ( K , N ) are verified in Figure 7. We can see that the unified expressions of the SLE match the empirical results very well; even the dimensions are very small ( K = 3 , N = 6 ). The numerical results in Figure 5, Figure 6 and Figure 7 indicate that the coefficient-based expressions with the proposed coefficient matrices (vectors) fit the empirical results very well, even though the matrix dimensions are finite.

6.2. Simulations Results for Sensing Performances

The sensing performances of the FRMT-based scheme, the IRMT-based scheme and the conventional cooperative SS scheme are shown in Figure 8, in which the number of the sensors is K = 3 , and the sample number per sensor is N = 6 . The SNR of PU signal is set to −5 dB. For the FRMT-based scheme, the SLE-based scheme is exploited. The largest eigenvalue calculated with the TW-law of the IRMT is used. For the conventional cooperative SS scheme, the FC gathers PU samples from the distributed sensors and determines the final decision with the energy detection method. The FRMT-based scheme and IRMT-based scheme can achieve better sensing performance compared to the conventional cooperative SS scheme, especially for low P F . Moreover, the FRMT-based scheme is superior to the IRMT-based scheme, because the precise sensing thresholds can be determined by the exact eigenvalue distributions. There is about a 20 % performance gain comparing the FRMT-based scheme with the IRMT-based scheme when P F = 0.1 .
However, the computational complexity of the RMT-based scheme is higher than that of the conventional cooperative SS scheme. The RMT-based schemes have to calculate the eigenvalues of sample matrices. For the SLE-based scheme, all K eigenvalues and the trace R should be calculated. The computational complexity of the eigenvalue calculation is about O K N 3 . For the IRMT-based scheme, the largest eigenvalue λ ¯ E K of Equation (14) has to be determined. As for the conventional cooperative SS scheme, the computational complexity is about O K N 2 under the assumption that the energy detection scheme is used at the FC.
The sensing performances of the proposed SS schemes based on the FRMT are evaluated in the following simulations. The receiver operating characteristic (ROC) performances, the probability of detection ( P D ) against the probability of false alarm ( P F ), are used to evaluate the SS schemes based on the FRMT. Figure 9 indicates that the SLE-based scheme outperforms the SCN-based scheme about 13 % when δ = 0.1 and the number of PU samples is 24. The thresholds in this simulation are precisely calculated with the exact distributions of the SCN and the SLE.
Figure 10 compares two SS schemes in the same sensing circumstance when P F = 0.05 and the total number of the samples is set to 18 and 24. When the SNR is 0 dB and K = 3 , N = 8 , the sensing performances of the SLE-based and SCN-based schemes are about 70 % and 45 % , respectively. There is about a 25 % performance gain. Moreover, the sensing performance of the SLE-based scheme with K = 4 , N = 6 and 0 dB is about 74 % , and there is only a 4 % performance gain compared to the case of K = 3 , N = 8 . From the results in Figure 10, we can see that the SLE-based scheme outperforms the SCN-based scheme when K > 3 under the FRMT paradigm.
In Figure 11, the sensing performance of the LE-based scheme is compared to those of the SLE-based scheme and the SCN-based scheme for K = 3 , N = 6 with P F = 0.05 . The LE-based scheme outperforms both the SLE-based scheme and SCN-based scheme for varying PU SNRs under the condition that the sensing thresholds of the LE-based scheme are predefined in an empirical way. However, like the energy detection scheme, the LE-based scheme also suffers from the noise uncertainty problem, which generates imprecise thresholds. The noise uncertainty problem of the LE-based scheme is illustrated by the simulations in Figure 12, in which the theoretical thresholds generated from the distributions of the LE cannot match the empirical thresholds very well. This means that we cannot calculate the exact thresholds using the distributions of the LE due to the noise uncertainty problem.

7. Conclusions

The coefficient-based distributions of the FRMT have been discussed in this paper. The coefficient reuse mechanism has been studied. The same coefficient matrices can be used by the distributions of the SE and the SLE. Moreover, the distributions of the SE and the DCN can share the same coefficient vectors. The FRMT-based SS schemes can achieve better sensing performance, compared to the IRMT-based schemes and the conventional SS schemes. The numerical results have verified that the coefficient-based distributions of the FRMT can fit the empirical results very well. The theoretical analysis and the simulation results have indicated that the SLE-based scheme is an asymptotically optimal scheme in the FRMT paradigm when the dimensions of the Wishart matrix are relatively large, i.e., K > 2 .

Acknowledgments

The authors acknowledge the support from the National Natural Science Foundation of China (Grant No. 61371110), the General and Special Program of China Postdoctoral Science Foundation (Grant Nos. 2012M521334 and 2013T60669), the Outstanding Young Scientist Research Award Foundation of Shandong Province (Grant No. BS2013DX004), the 863 project in 5G wireless networking, Ministry of Science and Technology of China (Grant No. 2014AA01A701), the EU FP7 QUICKproject (Grant No. PIRSES-GA-2013-612652) and the EU H2020 ITN5G Wireless project (Grant No. 641985).

Author Contributions

Wensheng Zhang, Cheng-Xiang Wang, Xiaofeng Tao and Piya Patcharamaneepakorn conceived of the main algorithms. Wensheng Zhang deduced the main mathematical results, and Xiaofeng Tao provided some suggestions. Wensheng Zhang wrote the paper, and Cheng-Xiang Wang and Piya Patcharamaneepakorn revised it.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The coefficient matrices P ( K , N ) and C ( K , N ) for finite dimensions K and N can be calculated with the generation algorithm given in [25]. The coefficient matrices P ( K , N ) and C ( K , N ) ( K = 3 , N K ) are indicated in Table A1 and Table A2, respectively. Note that only the cases of K = 3 , N = 3 , 4 , 5 , 6 are provided for simplicity.
Table A1. Coefficient matrices P ( K = 3 , N = 3 , 4 , 5 , 6 ) .
Table A1. Coefficient matrices P ( K = 3 , N = 3 , 4 , 5 , 6 ) .
N = 3N = 4
−3−66−2 1 4 6−8 9 2 −1 1 12 0
−66−3−1 - 1 2 −1241−1 - 1 12 - 1 12
3000064 1 2 000
N = 5
5−52 - 1 3 1 48 00
−10010 - 1 8 0 - 1 144
552 1 3 1 48 00
N = 6
5 2 −2 5 8 - 1 2 1 240 000
−5−1 1 4 1 12 - 1 120 - 1 120 1 2880 - 1 2880
5 2 3 13 8 1 2 7 80 1 120 1 2880 0
Table A2. Coefficient matrices C ( K = 3 , N = 3 , 4 , 5 , 6 ) .
Table A2. Coefficient matrices C ( K = 3 , N = 3 , 4 , 5 , 6 ) .
N= 3N= 4
100001000000
−30−31 - 1 4 −3−3 3 2 - 13 6 7 12 - 1 12 0
3031 1 4 360 4 3 11 12 1 6 1 24
−10000−1−3 - 3 2 - 1 6 000
N = 5
100000000
−3−3 - 3 2 7 6 - 23 24 5 24 - 1 48 00
366 2 3 1 3 1 3 1 9 1 72 1 288
−1−3 - 9 2 - 17 6 - 7 8 - 1 8 - 1 144 00
N = 6
10000000000
−3−3 - 3 2 - 1 2 1 2 - 3 10 13 240 - 1 240 000
3664 3 4 1 10 3 40 1 30 7 960 1 1440 1 5760
−1−3 - 9 2 - 9 2 - 11 4 - 21 20 - 61 240 - 3 80 - 1 320 - 1 8640 0
Based on the coefficient matrices in the above tables and the corresponding PDF in Equation (26) and CDF in Equation (27), the exact formulations can conveniently be generated. Let K = 3 and N = 4 , the coefficient matrices P ( 3 , 4 ) and C ( 3 , 4 ) can be given as:
P ( 3 , 4 ) = 6 - 8 9 2 - 1 1 12 0 - 12 4 1 - 1 - 1 12 - 1 12 6 4 1 2 0 0 0
C ( 3 , 4 ) = 1 0 0 0 0 0 0 - 3 - 3 3 2 - 13 6 7 12 - 1 12 0 3 6 0 4 3 11 12 1 6 1 24 - 1 - 3 - 3 2 - 1 6 0 0 0
and the corresponding PDF and CDF can be formulated as:
f λ K ( x ) = ( 6 e - 3 x - 12 e - 2 x + 6 e - x ) x + ( 4 e - 3 x + 4 e - 2 x - 8 e - x ) x 2 + ( 1 2 e - 3 x + e - 2 x + 9 2 e - x ) x 3 - ( e - 2 x + e - x ) x 4 ( 1 12 e 2 x 1 12 e x ) x 5 1 12 e 2 x x 6
F λ K ( x ) = 1 - e - 3 x + 3 e - 2 x - 3 e - x - ( 3 e - 3 x - 6 e - 2 x + 3 e - x ) x - ( 3 2 e - 3 x - 3 2 e - x ) x 2 - ( 1 6 e - 3 x - 4 3 e - 2 x + 13 6 e - x ) x 3 + ( 11 12 e 2 x + 7 12 e x ) x 4 + ( 1 6 e 2 x 1 12 e x ) x 5 + 1 24 e 2 x x 6

Appendix B

The coefficient vectors p ( K , N ) and c ( K , N ) for finite dimensions K and N are provided in Table B1 as examples. Let K = 3 and N = 5 , the coefficient vectors p ( 3 , 5 ) and c ( 3 , 5 ) can be given as:
p ( 3 , 5 ) = 5 5 2 1 3 1 48
c ( 3 , 5 ) = - 1 - 3 - 9 2 - 17 6 - 7 6 - 1 8 - 1 144
Based on the coefficient vectors and the corresponding PDF in Equation (30) and CDF in Equation (31), the corresponding formulations can conveniently be written as:
f λ 1 ( x ) = 5 e - 3 x x 2 + 5 e - 3 x x 3 + 2 e - 3 x x 4 + 1 3 e - 3 x x 5 + 1 48 e - 3 x x 6
F λ 1 ( x ) = 1 - e - 3 x - 3 e - 3 x x - 9 2 e - 3 x x 2 - 17 6 e - 3 x x 3 - 7 6 e - 3 x x 4 - 1 8 e - 3 x x 5 - 1 144 e - 3 x x 6
Table B1. Coefficient matrices p ( K = 3 , N = 3 , 4 , 5 , 6 ) and c ( K = 3 , N = 3 , 4 , 5 , 6 ) .
Table B1. Coefficient matrices p ( K = 3 , N = 3 , 4 , 5 , 6 ) and c ( K = 3 , N = 3 , 4 , 5 , 6 ) .
Coefficient Matrices p ( K = 3 , N = 3 , 4 , 5 , 6 )
N = 3N = 4N = 5
364 1 2 552 1 3 1 48
N = 6
5 2 3 13 8 1 2 7 80 1 120 1 2880
Coefficient Matrices c ( K = 3 , N = 3 , 4 , 5 , 6 )
N = 3N = 4
−1−1−3 - 3 2 - 1 6
N = 5
−1−3 - 9 2 - 17 6 - 7 6 - 1 8 - 1 144
N = 6
−1−3 - 9 2 - 9 2 - 11 4 - 21 20 - 61 240 - 3 80 - 1 320 - 1 8640

References

  1. Haykin, S. Cognitive radio: Brain-empowered wireless communications. IEEE J. Select. Areas Commun. 2005, 23, 201–220. [Google Scholar] [CrossRef]
  2. Akyildiz, I.F.; Lee, W.-Y.; Vuran, M.C.; Mohanty, S. Next generation/dynamic spectrum access/cognitive radio wireless networks: A survey. Comput. Netw. J. 2006, 50, 2127–2159. [Google Scholar] [CrossRef]
  3. Hong, X.; Wang, J.; Wang, C.-X.; Shi, J. Cognitive radio in 5g: A perspective on energy-spectral efficiency trade-off. IEEE Commun. Mag. 2014, 52, 46–53. [Google Scholar]
  4. Yücek, T.; Arslan, H. A survey of spectrum sensing algorithms for cognitive radio applications. IEEE Commun. Surv. Tutor. 2009, 11, 116–130. [Google Scholar]
  5. Sun, H.; Nallanathan, A.; Wang, C.-X.; Chen, Y. Wideband spectrum sensing for cognitive radio networks: A survey. IEEE Wirel. Commun. 2013, 20, 74–81. [Google Scholar]
  6. Tulino, A.M.; Verdu, S. Random Matrix Theory and Wireless Communications; Now Publishers Inc.: Hanover, MA, USA, 2004; pp. 1–182. [Google Scholar]
  7. Couillet, R.; Debbah, M. Random Matrix Methods for Wireless Communications; Cambridge University Press: Cambridge, UK, 2011. [Google Scholar]
  8. Kang, M.; Alouini, M.-S. Largest eigenvalue of complex Wishart matrices and performance analysis of MIMO MRC systems. IEEE J. Select. Areas Commun. 2003, 21, 418–426. [Google Scholar]
  9. Matthaiou, M.; McKay, M.; Smith, P.; Nossek, J. On the condition number distribution of complex wishart matrices. IEEE Trans. Commun. 2010, 58, 1705–1717. [Google Scholar]
  10. Zanella, A.; Chiani, M.; Win, M. On the marginal distribution of the eigenvalues of wishart matrices. IEEE Trans. Commun. 2010, 57, 1050–1060. [Google Scholar]
  11. Jin, S.; McKay, M.; Gao, X.; Collings, I. MIMO multichannel beamforming: SER and outage using new eigenvalue distributions of complex noncentral wishart matrices. IEEE Trans. Commun. 2010, 56, 424–434. [Google Scholar]
  12. Edelman, A. Eigenvalues and Condition Numbers of Random Matrices. Ph.D. Thesis, MIT, Cambridge, MA, USA, May 1989. [Google Scholar]
  13. Baik, J.; Silverstein, J.W. Eigenvalues of large sample covariance matrices of spiked population models. J. Multivar. Anal. 2006, 97, 1382–1408. [Google Scholar]
  14. Bai, Z.; Silverstein, J.W. Spectral Analysis of Large Dimensional Random Matrices; Springer: New York, NY, USA, 2010. [Google Scholar]
  15. Marchenko, V.A.; Pastur, L.A. Distributions of eigenvalues for some sets of random matrices. Math. USSR Sb. 1967, 1, 457–483. [Google Scholar]
  16. Tracy, C.; Widom, H. On orthogonal and sympletic matrix ensembles. Commun. Math. Phys. 1996, 177, 727–754. [Google Scholar]
  17. Abreu, G.; Zhang, W.; Sanada, Y. Finite random matrices for blind spectrum sensing. In Proceedings of the Forty Fourth Asilomar Conference on Signals, Systems and Computers (ASILOMAR), Pacific Grove, CA, USA, 7–10 November 2010; pp. 116–120.
  18. Kortun, A.; Ratnarajah, T.; Sellathurai, M.; Liang, Y.; Zeng, Y. On the eigenvalue-based spectrum sensing and secondary user throughput. IEEE Trans. Veh. Technol. 2014, 63, 1480–1486. [Google Scholar]
  19. Huang, L.; Fang, J.; Liu, K.; So, H.; Li, H. An eigenvalue-moment-ratio approach to blind spectrum sensing for cognitive radio under sample-starving environment. IEEE Trans. Veh. Technol. 2015, 64, 3465–3480. [Google Scholar]
  20. Dighe, P.; Mallik, R.; Jamuar, S. Analysis of transmit-receive diversity in rayleigh fading. IEEE Trans. Commun. 2003, 51, 694–703. [Google Scholar]
  21. Park, C.S.; Lee, K.B. Statistical multimode transmit antenna selection for limited feedback MIMO systems. IEEE Trans. Wirel. Commun. 2008, 7, 4432–4438. [Google Scholar]
  22. Wei, L.; Tirkkonen, O. Analysis of scaled largest eigenvalue based detection for spectrum sensing. In Proceedings of the IEEE International Conference on Communications (ICC), Kyoto, Japan, 5–9 June 2011; pp. 1–5.
  23. Wei, L.; Tirkkonen, O.; Dharmawansa, K.D.P.; McKay, M.R. On the exact distribution of the scaled largest eigenvalue. In Proceedings of the IEEE International Conference on Communications (ICC), Ottawa, ON, Canada, 10–15 June 2012; pp. 2422–2426.
  24. Kortun, A.; Sellathurai, M.; Ratnarajah, T.; Zhong, C. Distribution of the ratio of the largest eigenvalue to the trace of complex wishart matrices. IEEE Trans. Signal Process. 2012, 60, 5527–5532. [Google Scholar]
  25. Abreu, G.; Zhang, W. Extreme eigenvalue distributions of finite random wishart matrices with application to spectrum sensing. In Proceedings of the Forty Fifth Asilomar Conference on Signals, Systems and Computers (ASILOMAR), Pacific Grove, CA, USA, 6–9 November 2004; pp. 772–776.
  26. Cabric, D.; Mishra, S.; Brodersen, R. Implementation issues in spectrum sensing for cognitive radios. In Proceedings of the Thirty Eighth Asilomar Conference on Signals, Systems and Computers (ASILOMAR), Pacific Grove, CA, USA, 7–10 November 2004; pp. 772–776.
  27. Unnikrishnan, J.; Veeravalli, V. Cooperative sensing for primary detection in cognitive radio. IEEE J. Sel. Areas Commun. 2008, 21, 418–426. [Google Scholar]
  28. Li, Z.; Yu, F.; Huan, M. A cooperative spectrum sensing consensus scheme in cognitive radios. In Proceedings of the IEEE INFOCOM 2009, Rio de Janeiro, Brazil, 19–25 April 2009; pp. 2546–2550.
  29. Ding, G.; Wang, J.; Wu, Q.; Zhang, L.; Zou, Y.; Yao, Y.; Chen, Y. Robust spectrum sensing with crowd sensors. IEEE J. Select. Areas Commun. 2014, 62, 3129–3143. [Google Scholar]
  30. Huang, X.; Hu, F.; Chen, H.; Wang, G.; Jiang, T. Intelligent cooperative spectrum sensing via hierarchical Dirichlet process in cognitive radio networks. IEEE J. Select. Areas Commun. 2015, 33, 771–787. [Google Scholar]
  31. Zhang, Z.; Zhang, W.; Zeadally, S.; Wang, Y.; Liu, Y. Cognitive radio spectrum sensing framework based on multi-agent architecture for 5G networks. IEEE Wirel. Commun. 2015, 22, 34–39. [Google Scholar]
  32. Cardoso, L.S.; Debbah, M.; Bianchi, P.; Najim, J. Cooperative spectrum sensing using random matrix theory. In Proceedings of the 3rd International Symposium on Wireless Pervasive Computing (ISWPC), Santorini, Greece, 7–9 May 2008; pp. 334–338.
  33. Zeng, Y.; Liang, Y.-C. Eigenvalue-based sectrum sensing algorithms for cognitive radio. IEEE Trans. Wirel. Commun. 2009, 57, 1784–1793. [Google Scholar]
  34. Penna, F.; Garello, R.; Spirito, M. Cooperative spectrum sensing based on the limiting eigenvalue ratio distribution in wishart matrices. IEEE Commun. Lett. 2009, 13, 507–509. [Google Scholar]
  35. Bianchi, P.; Debbah, M.; Maida, M.; Najim, J. Performance of statistical tests for single-source detection using random matrix theory. IEEE Trans. Inform. Theory 2011, 57, 2400–2419. [Google Scholar]
  36. Zhang, W.; Abreu, G.; Inamori, M.; Sanada, Y. Spectrum sensing algorithms via finite random matrices. IEEE Trans. Commun. 2012, 60, 164–175. [Google Scholar]
  37. Goldsmith, A.; Jafar, S.; Maric, I.; Srinivasa, S. Breaking spectrum gridlock with cognitive radios: An information theoretic perspective. Proc. IEEE 2009, 97, 894–914. [Google Scholar]
  38. Edelman, A. On the distribution of a scaled condition number. Math. Comput. 1992, 58, 185–190. [Google Scholar]
  39. Dieng, M. Distribution functions for edge eigenvalues in orthogonal and symplectic ensembles: Painlevé representations. Ph.D. Thesis, University of California, Davis, CA, USA, 2005. [Google Scholar]
  40. Curtiss, J.H. On the distribution of the quotient of two chance variables. Ann. Math. Statist. 1941, 12, 409–421. [Google Scholar]
  41. James, A.T. Distributions of matrix variates and latent roots derived from normal samples. Ann. Math. Statist. 1964, 35, 475–501. [Google Scholar]
  42. Wei, L.; McKay, M.; Tirkkonen, O. Exact demmel condition number distribution of complex wishart matrices via the mellin transform. IEEE Commun. Lett. 2011, 15, 175–177. [Google Scholar]
  43. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products, 7th ed.; Academic Press: Burlington, MA, USA, 2007. [Google Scholar]
  44. Nadler, B. On the distribution of the ratio of the largest eigenvalue to the trace of a wishart matrix. J. Multivar. Anal. 2011, 102, 363–371. [Google Scholar]
  45. Maaref, A.; Aissa, S. Closed-form expressions for the outage and ergodic Shannon capacity of MIMO MRC systems. IEEE Trans. Commun. 2005, 53, 1092–1095. [Google Scholar]
  46. Johnstone, M. On the distribution of the largest eigenvalue in principal components analysis. Ann. Stat. 2001, 29, 295–327. [Google Scholar]
  47. Zeitouni, O.; Ziv, J.; Merhav, N. When is the generalized likelihood ratio test optimal? IEEE Trans. Inform. Theory 1992, 38, 1597–1602. [Google Scholar]
Figure 1. The system model with the PU system consisting of the PU base-station and some PUs (the number is larger than one) and the SU system consisting of the fusion centre (FC) and K SUs working as the sensors.
Figure 1. The system model with the PU system consisting of the PU base-station and some PUs (the number is larger than one) and the SU system consisting of the fusion centre (FC) and K SUs working as the sensors.
Sensors 16 01183 g001
Figure 2. Dimension boundary.
Figure 2. Dimension boundary.
Sensors 16 01183 g002
Figure 3. The Marchenko–Pastur distribution (PDF) for the eigenvalues of the central and complex Wishart matrix W with large dimensions.
Figure 3. The Marchenko–Pastur distribution (PDF) for the eigenvalues of the central and complex Wishart matrix W with large dimensions.
Sensors 16 01183 g003
Figure 4. The Tracy–Widom distribution (PDF) for the normalized and centralized extreme (largest) eigenvalue of the central and complex Wishart matrix W with large dimensions.
Figure 4. The Tracy–Widom distribution (PDF) for the normalized and centralized extreme (largest) eigenvalue of the central and complex Wishart matrix W with large dimensions.
Sensors 16 01183 g004
Figure 5. The PDFs of the extreme (largest and smallest) eigenvalues of finite Wishart matrices.
Figure 5. The PDFs of the extreme (largest and smallest) eigenvalues of finite Wishart matrices.
Sensors 16 01183 g005
Figure 6. The PDFs of the DCN for the central and complex Wishart matrix W with finite dimensions.
Figure 6. The PDFs of the DCN for the central and complex Wishart matrix W with finite dimensions.
Sensors 16 01183 g006
Figure 7. The PDFs of the SLE for the central and complex Wishart matrix W with finite dimensions.
Figure 7. The PDFs of the SLE for the central and complex Wishart matrix W with finite dimensions.
Sensors 16 01183 g007
Figure 8. The comparisons of the FRMT-based scheme, the IRMT-based scheme and the conventional cooperative spectrum sensing (SS) scheme. The SNR of PU signal is −5 dB.
Figure 8. The comparisons of the FRMT-based scheme, the IRMT-based scheme and the conventional cooperative spectrum sensing (SS) scheme. The SNR of PU signal is −5 dB.
Sensors 16 01183 g008
Figure 9. The comparisons of the SLE-based scheme and the SCN-based scheme with the same number of PU samples.
Figure 9. The comparisons of the SLE-based scheme and the SCN-based scheme with the same number of PU samples.
Sensors 16 01183 g009
Figure 10. The comparisons of the SLE-based scheme and the SCN-based scheme in the same sensing circumstance, P F = 5 % .
Figure 10. The comparisons of the SLE-based scheme and the SCN-based scheme in the same sensing circumstance, P F = 5 % .
Sensors 16 01183 g010
Figure 11. The sensing performances of the LE-based scheme, the SLE-based scheme and the SCN-based scheme for K = 3 , N = 6 . The P F is set to 0.05 .
Figure 11. The sensing performances of the LE-based scheme, the SLE-based scheme and the SCN-based scheme for K = 3 , N = 6 . The P F is set to 0.05 .
Sensors 16 01183 g011
Figure 12. The comparisons of the theoretical thresholds and the empirical thresholds of the LE-based scheme with K = 3 , N = 6 for varying P F . The number of the trials is 10,000.
Figure 12. The comparisons of the theoretical thresholds and the empirical thresholds of the LE-based scheme with K = 3 , N = 6 for varying P F . The number of the trials is 10,000.
Sensors 16 01183 g012
Table 1. Distributions of the random matrix theory (RMT). IRMT, infinite random matrix theory; FRMT, finite random matrix theory; SCN, standard condition number; DCN, Demmel condition number; SLE, scaled largest eigenvalue.
Table 1. Distributions of the random matrix theory (RMT). IRMT, infinite random matrix theory; FRMT, finite random matrix theory; SCN, standard condition number; DCN, Demmel condition number; SLE, scaled largest eigenvalue.
VariableIRMTFRMT
PDFCDFPDFCDF
Largest Eigenvalue: λ K Equation (12)Equation (13)Equation (26)Equation (27)
Smallest Eigenvalue: λ 1 Equation (12)Equation (13)Equation (30)Equation (31)
SCN: ξEquation (18)Equation (19)Equation (32) aEquation (33) a
DCN: κPDF of κ bEquation (37)Equation (38)
SLE: ψno results cEquation (42)Equation (43)
a, No exact results for K > 2 ; b, only for K × K matrices; c, to the best our knowledge, there are no corresponding results in the IRMT.

Share and Cite

MDPI and ACS Style

Zhang, W.; Wang, C.-X.; Tao, X.; Patcharamaneepakorn, P. Exact Distributions of Finite Random Matrices and Their Applications to Spectrum Sensing. Sensors 2016, 16, 1183. https://doi.org/10.3390/s16081183

AMA Style

Zhang W, Wang C-X, Tao X, Patcharamaneepakorn P. Exact Distributions of Finite Random Matrices and Their Applications to Spectrum Sensing. Sensors. 2016; 16(8):1183. https://doi.org/10.3390/s16081183

Chicago/Turabian Style

Zhang, Wensheng, Cheng-Xiang Wang, Xiaofeng Tao, and Piya Patcharamaneepakorn. 2016. "Exact Distributions of Finite Random Matrices and Their Applications to Spectrum Sensing" Sensors 16, no. 8: 1183. https://doi.org/10.3390/s16081183

APA Style

Zhang, W., Wang, C. -X., Tao, X., & Patcharamaneepakorn, P. (2016). Exact Distributions of Finite Random Matrices and Their Applications to Spectrum Sensing. Sensors, 16(8), 1183. https://doi.org/10.3390/s16081183

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

Article Metrics

Back to TopTop