Eigenbased Multi-Antenna Spectrum Sensing: Experimental Validation on a Software-Defined Radio Testbed

TL;DR

This paper experimentally validates multi-antenna spectrum sensing techniques for 6G networks using a software-defined radio testbed.

Contribution

The paper provides experimental validation and computation of GLRT false-alarm probability for four-antenna systems.

Findings

RLRT outperforms other detection schemes in terms of detection capability.

Experimental results align with simulation outcomes for SS techniques.

GLRT false-alarm probability is computed for the first time in explicit form for four-antenna systems.

Abstract

Spectrum Sensing (SS) is expected to play a crucial role in forthcoming 6G Cognitive Radio Networks (CRNs), where unlicensed users will be able to dynamically access the spectrum and perform opportunistic transmissions without generating interference for licensed users. In this work, we investigate multiple-antenna SS techniques by analyzing the performance of several widely used detection schemes—namely, Roy’s Largest Root Test (RLRT), the Generalized Likelihood Ratio Test (GLRT), the Eigenvalue Ratio Detector (ERD), and the Energy Detector (ED)—under varying false-alarm probabilities and signal-to-noise ratios (SNRs). We assume there are a fixed number of sensors at the secondary-user receiver, namely, four. To evaluate the behavior of these detectors in realistic conditions, we developed a software-defined radio (SDR) testbed using Universal Software Radio Peripherals (USRPs), enabling both primary-user signal transmission and secondary-user data acquisition. The experimental results, illustrated through Receiver Operating Characteristic (ROC) and performance curves, are compared with simulation outcomes. The analysis is complemented by a detailed state-of-the-art listing of the available analytical characterizations of the false-alarm probabilities for the considered SS schemes. In particular, the GLRT false-alarm probability, previously unavailable in explicit form for a four-antenna equipped receiver, is computed as well. These results validate the superior detection capability of RLRT over the other schemes tested, confirming its effectiveness not only in theoretical analysis but also in practical SDR-based implementations.

Figures

Full text

1. Introduction

Cognitive Radio Networking (CRN) [1] was conceived as a paradigm through which to counteract the growing shortage of available spectrum resources by promoting the intelligent coexistence of licensed and unlicensed users [2]. Through dynamic sensing and adaptation, CRN allows secondary systems to opportunistically exploit underused frequency bands without causing harmful interference for primary users. Owing to these capabilities, CRN is expected to become a key enabling technology for future 6G communication infrastructures, which will demand unprecedented flexibility, spectrum efficiency, and context awareness.

At the heart of CRN’s functionality lies Spectrum Sensing (SS), a critical process responsible for identifying whether a given frequency band is occupied by a primary transmitter (the busy condition) or currently vacant (in an idle state). The reliability of this sensing process directly impacts spectrum utilization efficiency and interference avoidance. As such, the development of accurate, fast, and robust SS algorithms has remained one of the cornerstone challenges in cognitive radio research over the past two decades [3,4,5]. These efforts have led to a wide range of techniques capable of coping with low signal-to-noise ratios (SNRs), channel variability, and practical hardware constraints.

The present work focuses on a measurement-based evaluation of several well-established SS strategies that operate based on the statistical properties of the received signal covariance matrix [6,7,8,9]. In systems equipped with multiple antennas, such eigenvalue-based approaches exploit spatial diversity to enhance detection reliability. The algorithms considered throughout this work include both statistically grounded procedures, such as the Generalized Likelihood Ratio Test (GLRT) [10] and Roy’s Largest Root Test (RLRT) [11], and more heuristic alternatives, such as the Energy Detector (ED) [12] and the Eigenvalue Ratio Detector (ERD) [6]. While the former rely on rigorous estimation-theoretic principles, the latter are often favored for their implementation simplicity in practical scenarios. In particular, in this work, we extend the content in [13] in two main directions:

• With reference to statistically grounded procedures, we provide a detailed analytical characterization of the false-alarm probabilities for all SS schemes presented both at finite size as well as in the asymptotic regime, and, especially for GLRT, we elaborate on the findings of [14] to provide a more explicit performance evaluation for the finite-dimensional case;
• We add a further heuristic scheme, the above mentioned ERD, to the investigation and, after providing a semi-asymptotic analysis of its performance, carry out a performance evaluation of this SS algorithm both experimentally (i.e., using measurements) as well as via simulations.

The practical implementation of CRN principles has been greatly facilitated by the emergence of software-defined radio (SDR) technology [1,15], which offers a highly flexible and reconfigurable platform for real-time signal acquisition, processing, and control (for a recent account on the topic, see [16]). By shifting signal-processing functionalities from rigid hardware components to programmable software modules [17], SDR devices enable dynamic reconfiguration of the radio front end in response to variations in the surrounding electromagnetic environment—an essential requirement for effective CRN [18]. This inherent adaptability allows radios to modify transmission parameters, waveforms, and protocols on the fly, thereby supporting intelligent spectrum awareness and decision making. Among the available SDR platforms, the Universal Software Radio Peripheral (USRP) has gained particular prominence due to its cost-effectiveness, modular hardware architecture, and seamless compatibility with widely used open-source software frameworks [19,20]. These features make the USRP a popular and reliable choice for the prototyping, implementation, and experimental evaluation of CRN functionalities. The close synergy between CRN concepts and SDR platforms further enables the experimental validation of theoretical spectrum sensing and adaptation strategies, effectively bridging the gap between analytical modeling and real-world performance assessment. This convergence marks a significant step toward the realization of more adaptive, autonomous, and intelligent radio systems capable of continuous self-optimization in dynamic wireless environments.

Extensive experimental research has been reported in the literature to validate SS algorithms using SDR-based setups, particularly focusing on well-known implementations of the ED [12,21,22,23] and cyclostationary feature detectors [24]. Building upon such foundational studies, the present contribution offers a comparative performance analysis of eigenvalue-based detection schemes, including the GLRT, the RLRT, and the ERD, alongside the classical ED, using both simulated and experimentally collected datasets. The experimental framework, designed around a USRP-based receiver, facilitates reproducible real-world testing while capturing the stochastic nature of wireless propagation and hardware-induced imperfections.

From an analytical perspective, the statistical performance of the above detectors has been rigorously studied under various assumptions regarding the number of antennas, observation samples, and channel conditions. Random Matrix Theory has emerged as a powerful instrument for deriving closed-form or asymptotic characterizations of key performance indicators, notably the probabilities of false alarms and detection. These analyses typically fall into two complementary regimes. The first corresponds to the finite-sample setting (see, e.g., [25]), where both the number of antennas and the number of time-domain snapshots are limited, allowing for exact, but often cumbersome, derivations. The second regime considers the large-system limit [26], in which the number of sensors and samples grows simultaneously while a constant ratio is maintained. This asymptotic framework provides insightful scaling laws and compact expressions for detection thresholds, offering valuable guidance for system design and parameter tuning.

Both the finite-sized and asymptotic approaches offer insight into the dependence of sensing performance on the SNR, the number of antennas, and the observation interval length. However, due to their simplicity with respect to the finite sample settings, asymptotic formulations have become an indispensable analytical tool, complementing experimental investigations and serving as the foundation for robust CRN algorithm design. The present study contributes to this ongoing line of research by combining theoretical understanding, numerical simulation, and real-world measurements to provide a comprehensive assessment of eigenvalue-based spectrum sensing within a practical multiple-antenna SDR framework.

This work is structured as follows: Section 2 relates the development of the system model and describes the sensing schemes. It also provides an outline of the corresponding analytical performance computation. The experimental setup and numerical results are reported in Section 4, while Section 5 concludes the work.

Notation

Throughout the paper, we denote vectors and matrices with boldface lowercase and uppercase letters, respectively. A hat on the variable denotes an estimated quantity (either numerically or measurement-based). We use $[eqn]$ to denote the statistical expectation. $[eqn]$ denotes the Euclidean norm of a vector, and $[eqn]$ denotes the Frobenius norm of a matrix, while we use $[eqn]$ to denote the trace of a square matrix. For matrices and vectors, $[eqn]$ indicates the conjugate transpose. Finally, $[eqn]$ is the determinant of a square matrix.

2. System Model

We focus on an SS setting where a receiver equipped with K antennas collects N time samples from each antenna. We use $[eqn]$ to denote the $[eqn]$ received vector at time $[eqn]$ , with generic entry $[eqn]$ denoting the discrete baseband complex sample at the k-th receiving antenna.

Under $[eqn]$ , namely, in the absence of a primary signal, $[eqn]$ is a vector of complex Gaussian noise samples with zero mean and variance $[eqn]$

$$[eqn]$$

where $[eqn]$ . Under $[eqn]$ , instead, both the primary signal and noise are present; therefore,

$$[eqn]$$

where $[eqn]$ is the transmitted signal sample, modeled, without a loss of generality, as Gaussian-distributed with zero mean and variance $[eqn]$ , while $[eqn]$ is the $[eqn]$ unknown complex channel vector, affected by Rayleigh flat fading. The channel is assumed to be memoryless and constant during the detection period; i.e., block fading is assumed. Under $[eqn]$ , we define the signal-to-noise ratio (SNR) at the receiver as

$$[eqn]$$

The detector collects the received samples in the $[eqn]$ matrix:

$$[eqn]$$

where $[eqn]$ is a $[eqn]$ signal vector, and $[eqn]$ is a $[eqn]$ noise matrix. The sample covariance matrix $[eqn]$ of the received samples is therefore given by

$$[eqn]$$

In our next analysis, we will make reference to the spectral decomposition of $[eqn]$ , with $[eqn]$ denoting the random unitary matrix of the eigenvectors and $[eqn]$ corresponding to the diagonal matrix of the random eigenvalues, whose largest value is denoted as $[eqn]$ and whose smallest value is denoted as $[eqn]$ , throughout the work.

Crucial to our signal detection task is, in particular, the exploitation of the largest eigenvalue. Indeed, with reference to the scenario where a primary signal is present, and upon defining, for given system parameters, the quantities

$$[eqn]$$ $$[eqn]$$

and

$$[eqn]$$

one can describe a phase transition phenomenon for the largest the eigenvalue of $[eqn]$ (see [27]), depending on the relative strength of the primary signal, such that, almost surely (a.s.)

$$[eqn]$$

3. Spectrum-Sensing Schemes and Performance Computation

The test statistic employed by the detector to discriminate between the null hypothesis $[eqn]$ and the presence of the primary signal $[eqn]$ is denoted as T; to make a decision, the detector compares it against a pre-defined threshold t: if $[eqn]$ , it chooses $[eqn]$ ; otherwise, it chooses $[eqn]$ . As a consequence, the probability of a false alarm

$$[eqn]$$

and the probability of detection is defined as

$$[eqn]$$

Hereinafter, we detail the SS schemes exploited in the experimental setups. The first one, the ED, exploits the difference in energy levels corresponding to the two alternative scenarios of the presence vs. absence of a primary signal. In turn, the ERD relies on the increased value of the ratio between the largest and smallest eigenvalue under $[eqn]$ w.r.t. $[eqn]$ . The remaining two schemes are based on classical estimation-theoretic criteria and depend on the statistical behavior of the largest eigenvalue of the sample covariance matrix of the collected data.

3.1. Energy Detector

The test statistic $[eqn]$ computes the average energy of the received signal over the N collected samples, normalized by the noise variance $[eqn]$ , namely,

$$[eqn]$$

Note that, since, according to the definition of the Frobenius norm, $[eqn]$ , $[eqn]$ can be expressed in terms of the spectrum of the sample covariance, as follows:

$$[eqn]$$

For this test, since the law of the normalized trace $[eqn]$ can be determined, e.g., as in Equation (8) of [14], by solving

$$[eqn]$$

the corresponding $[eqn]$ can be written as

$$[eqn]$$

Receiver operating characteristic (ROC) curves are computed based on both values of $[eqn]$ and $[eqn]$ , while the inversion of (15) leads to a proper threshold setting for performance computation. For eigenvalue-based schemes, as anticipated in the Introduction, there will be a need to resort to an asymptotic approach to carry out a performance analysis.

3.2. ERD

The idea behind the ERD is that the presence of a primary-user signal has a greater impact on the largest eigenvalue of the covariance of the collected samples than on the corresponding smallest one [6].

As a consequence, the test statistic of this sensing scheme is expressed as

$$[eqn]$$

Therefore, its probability law coincides with that of the so-called standard condition number of the matrix $[eqn]$ .

While the statistical distribution of (16) for a dual-antenna receiver (i.e., when $[eqn]$ ) is known thanks to [28], its closed-form expression for arbitrary values of K can be found in Theorem 1 of [29].

The corresponding false-alarm probability, for a given threshold $[eqn]$ , can be written as

$$[eqn]$$

where

$$[eqn]$$

and the $[eqn]$ matrix $[eqn]$ has the following entries:

$$[eqn]$$

In order to assess the analytical performance of the ERD, we resort to an asymptotic approach, where both the number of samples N and the number of receiving antennas K are assumed to grow large at the same rate; that is, the so-called aspect ratio of the random matrix $[eqn]$ , e.g., $[eqn]$ , is a number between 0 and 1 (for our purposes, indeed, $[eqn]$ ).

Under this assumption, we can exploit the convergence results for the extreme eigenvalues of $[eqn]$ , stating that

$$[eqn]$$ $$[eqn]$$

while, in distribution, upon defining the scaling parameters

$$[eqn]$$

and

$$[eqn]$$

both the random variables

$$[eqn]$$

and

$$[eqn]$$

converge to a random variable distributed according to a second order Tracy–Widom distribution (TW2), whose Cumulative Distribution Function (CDF) for the Gaussian Unitary Ensemble (GUE) is explicitly given by

$$[eqn]$$

where $[eqn]$ is the unique solution to the Painlevé II differential equation

$$[eqn]$$

satisfying the boundary condition

$$[eqn]$$

where $[eqn]$ denotes the Airy special function, one of the two linearly independent solutions to the following differential equation:

$$[eqn]$$

As a consequence of the above analysis, the Probability Density Function (PDF) of $[eqn]$ can be written as per Equations (12) and (13) of [6], relying on proper linear transformations of the $[eqn]$ law for the involved eigenvalues, namely,

$$[eqn]$$ $$[eqn]$$

and approximating, for the sake of simplicity, the behavior of $[eqn]$ and $[eqn]$ as if they were independent random variables in order to obtain the PDF and, therefore, the CDF of $[eqn]$ , e.g., $[eqn]$ .

Indeed, the PDF of $[eqn]$ , e.g., $[eqn]$ , can be obtained via Equation (14) of [6] as:

$$[eqn]$$

with $[eqn]$ being the indicator function of the inner interval. As a consequence of the above computation, threshold setting for ERD is performed upon computing $[eqn]$ .

3.3. RLRT

Roy’s test exploits as test statistics the largest eigenvalue of the sample covariance, normalized by the noise variance. Therefore,

$$[eqn]$$

$[eqn]$ can be immediately computed, for this sensing scheme, as

$$[eqn]$$

with $[eqn]$ as in (18), and where

$$[eqn]$$

Let us recall that, under $[eqn]$ , the CDF of $[eqn]$ is provided in Corollary 2 of [30], while, under $[eqn]$ , it assumes the expression derived in [31].

As discussed, a phase transition threshold exists for the identifiability of a signal with RLRT or GRLT, which can take place only if

$$[eqn]$$

After a phase transition, i.e., in case signals are actually identifiable, the scaled and normalized largest eigenvalue converges in distribution to a standard real Gaussian random variable [32], i.e.,

$$[eqn]$$

with

$$[eqn]$$

and

$$[eqn]$$

In the remainder of this section, for ease of notation, we forgo the dependence on $[eqn]$ and $[eqn]$ in both the expressions of $[eqn]$ and $[eqn]$ .

As a consequence of the convergence results for the largest eigenvalue, it turns out that an approximate expression for the threshold of RLRT, for a given tolerated false-alarm rate of $[eqn]$ , is

$$[eqn]$$

where $[eqn]$ is the inverse of the TW2 CDF (24). We also remark that, for single-signal detection, the presence of a phase transition phenomenon, for detectability under the alternate hypothesis, can be explicitly formulated in terms of a critical detection threshold, depending on the SNR, as follows:

$$[eqn]$$

3.4. GLRT

In the case of GLRT-based SS, the test statistic, derived by applying the GLR criterion to our detection problem, reads as

$$[eqn]$$

Its analytical characterization turns out to be challenging; however, even though the correlation between $[eqn]$ and the normalized trace in the denominator of (39) is intractable, as observed in [33], (39) and (13) are statistically independent. Therefore, the computation of an exact expression for both the PDF as well as the CDF for finite values of K and N has been provided in [14,34] (a thorough analysis of the literature in the asymptotic case can be found, e.g., in [10,26]). With reference to Equation (14) of [14], the false-alarm probability of the GLRT corresponding to a threshold $[eqn]$ can be expressed as

$$[eqn]$$

where

$$[eqn]$$

and the ways of expressing the values of $[eqn]$ are reported in the Appendix A. Under $[eqn]$ , a finite-sample statistical analysis of the GLRT can be performed by exploiting the results reported in [35], where the random variable $[eqn]$ is characterized. This is the subject of ongoing work. The setting of the threshold for GLRT exploits asymptotic arguments, as done for the RLRT in the previous sub-paragraph, since they are both affected by the phase transition phenomenon, impacting signal detectability. An effective approximation of the $[eqn]$ of the GLRT scheme can be found in [26]:

$$[eqn]$$

The definitions of $[eqn]$ and $[eqn]$ are given in the RLRT paragraph. Notice that (41) can be numerically inverted in an efficient way to find the required threshold $[eqn]$ .

3.5. Analytical vs. Simulated Statistics

The degree of matching between analytically computed and numerically simulated statistics of the considered SS schemes has been thoroughly investigated in the literature (see, e.g., [14,26,29]), discussing, in turn, issues regarding computing exact formulae obtained in finite-sized settings or the suitability of asymptotic characterizations for quantifying performance in real scenarios, involving a moderate number of sensors and/or processed samples. For the sake of ensuring the self-consistency of our analysis, in Figure 1, we depict both the analytical PDFs of all the sensing strategies considered (obtained via the formerly described asymptotic procedures), as well as their simulated counterparts, for different values of the ratio $[eqn]$ , ranging from $[eqn]$ (which is already suitable for resorting to random-matrix-theory-grounded analysis) to $[eqn]$ , a reduced-sample-size scenario that, while theoretically of interest, does not match usual SS settings, where the number of samples is quite a bit larger than the number of collecting sensors. The sensitivity of the ERD, RLRT, and GLRT to the behavior of the largest eigenvalue of the sample covariance of the channel can be easily remarked, whereas for the ED, only the variance changes, decreasing with increasing values of $[eqn]$ , while the expected value of the received energy remains constant.

4. Experimental Setup and Numerical Results

Figure 2 and Figure 3 show the SDR experimental setup, in which the multi-antenna SS algorithms were tested and validated through over-the-air real-time signal transmission and acquisition and subsequent data processing and performance evaluation. The setup is composed of three main components:

• The first component was a host computer, a SiComputer Extrema Workstations W200, equipped with an Intel Core i9-14900K CPU, with 3200 MHz and 24 cores, 64 GB of DDR4 RAM, and a Mellanox ConnectX-5 100 Gbps Ethernet adapter. The SDR software installed was MATLAB R2024b together with the Wireless Testbench support package for NI USRP radios.
• The second was an Ettus Research USRP X310, equipped with a Xilinx Kintex-7 FPGA, with two independent TX/RX channels, capable of streaming up to 200 MSamples/s per channel. The USRP X310 was used as single-antenna transmitter and connected to the workstation through a 10 Gbps Ethernet adapter. Notice that in order to obtain sufficiently low SNR levels suitable for testing the SS schemes at hand, we did not connect an external antenna to the transmitter.
• The third component was an NI Ettus USRP X410, equipped with a Xilinx Zynq Ultrascale+ ZU28DR RF System on Chip (RFSoC), with four independent TX/RX channels, capable of streaming up to 250 MSamples/s per channel. The USRP X410 was used as a multi-antenna receiver and was connected to the workstation through the QSFP28 interface using a 100 Gbps Ethernet adapter.

In order to assess the performance of the RLRT, GRLT, ERD and ED SS schemes, we assumed, as a primary user signal, there was a single carrier QPSK-modulated signal at a carrier frequency of $[eqn]$ GHz with a baud rate $[eqn]$ of 1 MBd/s. We used four samples per symbol; therefore, the USRP radios were set to a sample rate of $[eqn]$ Msamples/s. The narrow band was chosen for the primary signal to ensure there was a flat fading channel. Given the four channels of the USRP X410, $[eqn]$ , and we performed the experiment for both $[eqn]$ and $[eqn]$ samples.

Thresholds were computed previously for several values of $[eqn]$ and stored offline in a lookup table. While for the ED, the situation is quite straightforward, for RLRT, we used Momar Dieng’s MATLAB package “RMLab” [36], which provides the CDF $[eqn]$ of the TW2 distribution. We computed the centering and scaling parameters as in (20) and (22), and then we used the bisection method to retrieve the correct threshold $[eqn]$ value for any input $[eqn]$ . For ERD, we first used the same “RMLab” package to compute the scaled and centered PDFs $[eqn]$ and $[eqn]$ of the TW2 distribution as defined in (28) and (29), respectively, for different values of t; secondly, we obtained the PDF of the ERD as in (30) by means of numerical integration for each value of t; lastly, we obtained the CDF of ERD by numerically integrating the pdf, and we used the bisection method again to retrieve the correct threshold $[eqn]$ value for any input $[eqn]$ . Finally, for GLRT, we used the routines in [26], providing the inverted modified distribution of (41).

The SDR experiment can be condensed into three phases:

• 1.Noise variance estimation under $[eqn]$ : In this phase, no transmission occurred, and the USRP X410 acquired samples over its four channels for a total duration of $[eqn]$ s. Under the null hypothesis, the Maximum Likelihood estimation of the noise variance $[eqn]$ can be written as follows:

$$[eqn]$$

• 2.Primary user signal transmission and multi-antenna acquisition under $[eqn]$ : The USRP X310 was set up to continuously transmit the QPSK-modulated signal, and, simultaneously, the USRP X410 acquired samples under $[eqn]$ again for $[eqn]$ s. We can estimate the SNR $[eqn]$ at the receiver using the following approach:

$$[eqn]$$

• 3.Performance evaluation and comparison with simulated data: In the last phase, we processed the acquired samples into slots of N samples, and we computed the RLRT, ERD, GLRT, and ED test statistics. The results were generated through error counting by means of ROC curves for a specific SNR value and by means of a $[eqn]$ vs. SNR curve for a specific $[eqn]$ value. In particular, for a fixed $[eqn]$ , the points of each curve were obtained from the experimental setup by progressively increasing the transmit gain by $[eqn]$ dB at each step and repeating the procedure detailed in these three steps. Experimental results were compared against simulated ones; i.e., we performed a parallel Monte Carlo simulation by generating Gaussian samples with signal variance $[eqn]$ , and we applied the same test statistics and metric evaluations.

ROC curves for $[eqn]$ are reported in Figure 4, for all the considered SS schemes. Nonidealities in the experimental setup (such as interference during noise estimation, uneven noise floor levels among the four antennas, quantization noise, mulitpath fading, etc.) clearly impact the corresponding performance, which results in poorer performance than the simulated counterparts. As expected, RLRT outperformed all the other schemes, i.e., GLRT, in which there is no prior knowledge of the noise variance; the ED, which cannot make use of any information about the primary signal subspace; and, finally, ERD, which suffers from high noise sensitivity. Another relevant issue affecting the ERD is the underestimation of the denominator with respect to the case of test statistics based on estimation-theoretic criteria, such as the best-performing RLRT and the GLRT as well.

Among all nonidealities, multipath fading should be mitigated by the narrowband signal, while the uneven noise floor levels and received signal powers among the four independent channels play a pivotal role in the performance gap between the experimental and simulated results. A quantitative analysis of the unevenness of the received signal power between USRP antennas is reported in Figure 5, in which the relative full-scale received power levels are plotted against the input transmit radio gain for each USRP channel.

The performance curves representing the behavior of the probability of detection vs. SNR, at a target false-alarm probability of $[eqn]$ , are depicted in Figure 6 for $[eqn]$ . A very similar behavior was observed for all tests, as shown in Figure 4, as the performance ranking among the analyzed algorithms was retained.

A complementary perspective, i.e., the probability of misdetection ( $[eqn]$ ) vs. SNR, with $[eqn]$ , is reported in Figure 7, wherein, for $[eqn]$ , the curves do not hold to the already observed behavior over the whole SNR range due to the highly demanding target false alarm, combined with the processing of a small number of samples. However, RLRT still outperforms the remaining schemes, and ERD offers the poorest performance in this case as well.

The analysis is repeated for an increased sample size of $[eqn]$ ; in particular, ROC curves are reported in Figure 8, while the probability of detection vs. SNR is depicted in Figure 9. Performance ranking is unaffected by the increase of the number of collected samples, whereas the absolute performance benefits of the availability of a higher number of data. Remarkably, Figure 10 confirms that the performance ranking is fully retained for $[eqn]$ .

Overall, the experimental test statistics, with data acquired through SDR platforms, require approximately 0.5–1 more dB to achieve the same performance as the corresponding simulated ones.

As a final remark, we evaluated, based on simulated data only, the impact of employing a Gaussian primary signal on the performance of our tests and contrasted this, for each SS scheme, with our main working assumption, i.e., transmitting a QPSK-modulated signal. The curves, reported in Figure 11, confirm that, for the range of values of interest regarding the probability of detection, QPSK transmission is beneficial w.r.t. the use of a Gaussian signal.

5. Conclusions

An SDR-based experimental setup was developed to evaluate the performance of multi-antenna SS algorithms under realistic conditions. In this framework, a secondary user terminal equipped with $[eqn]$ receiving antennas was employed to collect measurement data. The acquired signals were processed to estimate the false-alarm and detection probabilities associated with four SS schemes, namely, ED, ERD, the GLRT, and RLRT. A thorough review of the mathematical literature was carried out to individuate, whenever available, explicit expressions for the false-alarm performance of each considered scheme, both in the finite-size as well as in an asymptotic regime, where the number of collected samples and that of the receiving sensors both tend toward infinity at the same rate. On top of this, an explicit characterization of the GLRT test statistics in the case of a four-sensor-equipped receiver is provided. The measurement-based performance curves were then compared with the results obtained from numerical simulations. The comparison demonstrates a very close match between experimental and simulated data, confirming the reliability of the proposed setup and the validity of the adopted simulation models. Future research plans include the investigation of finite-size statistical analysis of the presented eigenbased SS schemes under an alternate hypothesis (in the presence of the primary signal); moreover, our presented analysis could be extended to a higher number of sensors, upon a hardware update, and it could include testing on colored noise as well.