A Fast Power Spectrum Sensing Solution for Generalized Coprime Sampling

The growing scarcity of spectrum resources, wideband spectrum sensing is required to process a prohibitive volume of data at a high sampling rate. For some applications, spectrum estimation only requires second-order statistics. In this case, a fast power spectrum sensing solution is proposed based on the generalized coprime sampling. By exploring the sensing vector inherent structure, the autocorrelation sequence of inputs can be reconstructed from sub-Nyquist samples by only utilizing the parallel Fourier transform and simple multiplication operations. Thus, it takes less time than the state-of-the-art methods while maintaining the same performance, and it achieves higher performance than the existing methods within the same execution time, without the need for pre-estimating the number of inputs. Furthermore, the influence of the model mismatch has only a minor impact on the estimation performance, which allows for more efficient use of the spectrum resource in a distributed swarm scenario. Simulation results demonstrate the low complexity in sampling and computation, making it a more practical solution for real-time and distributed wideband spectrum sensing applications.


I. INTRODUCTION
T HE demand for spectrum resources is increasing due to the rapid development of low-orbit satellite constellation systems (e.g.SpaceX, OneWeb), 5G/6G networks, and the Internet of Things (IoT), etc [1]- [2].These applications are driving a growing demand for wideband spectrum sensing.Correspondingly, direct sampling requires a high-speed analog-to-digital converter (ADC) based on the Shannon-Nyquist sampling theorem, leading to a prohibitive volume of data and high energy cost.
Currently, the most widely used methods are sweep-tune sampling and filter band sampling [3], both of which fall under the category of low-speed sampling.Nevertheless, the scanning scheme has a detection latency and may miss shortlived signals.Additionally, the filter band scheme has a complicated structure and is prone to serious channel crosstalk.
Manuscript received November 23, 2023; revised month date, year; accepted month date, year.Date of publication month date, year; date of current version month date, year.This work was supported in part by the National Natural Science Foundation of China under Grant 62301119, and in part by the Key Project of the National Defense Science and Technology Foundation Strengthening Plan 173 under Grand 2022-JCJQ-ZD-010-12.The associate editor coordinating the review of this manuscript and approving it for publication was xxx.(Corresponding author: Kaili Jiang).
Digital Object Identifier As a result, there has been a trend towards using wideband spectrum sensing as a guide.
The recent compressive sensing (CS) theory provides a sub-sampling scheme with low-speed and large instantaneous bandwidth by utilizing the sparsity in the frequency domain [4]- [5].The typical CS schemes include analog-to-information converter (AIC), multi-coset sampling (MCS), and multi-rate sampling (MRS).However, the sparse signal recovery algorithms have high computational complexity and are extremely sensitive to noise.On the other hand, the engineering implementation of these schemes also presents a major difficulty, restricting their application.
To overcome these difficulties, the approach to spectrum estimation has been shifted processing the original signal to analyzing its second-order statistics.In this context, the compressive covariance sensing (CCS) theory provides a wideband spectrum sensing scheme that operates at low-speed and has a large instantaneous bandwidth.It is reliable even in environments with a low signal-to-noise ratio (SNR) and nonsparse conditions.The typical CCS methods include entropy functions minimization [6], matrix norm minimization [7], and Toeplize matrix completion [8], among others.However, all these methods are based on the reconstruction of the covariance matrix and the use of mutiple signal classification (MUSIC) algorithm, which have high computational complexity and require pre-estimation the number of signals.
Furthermore, a computationally efficient method is developed based on the relationship between the autocorrelation sequence and sub-sampling samples of the MCS scheme [9].Building on this, a fast solution for generalized coprime sampling is introduced, which utilizes only parallel FFT and multiplication operations.As a result, it achieves reduced time and low estimation error, presenting a tradeoff between system performance and the number of degrees of freedom (DOFs).Moreover, model mismatch has little effect on performance, making a more practical solution for real-time and distributed wideband spectrum sensing applications.The generalized coprime sampling architecture consists of two uniform sub-Nyquist sampling channels, whose sampling periods are coprime multiples of the Nyquist sampling period.The introduction of two other operations, the multiple coprime unit factor p ∈ N + and the non-overlapping factor q ∈ N + , increases the number of DOFs and improves the estimation accuracy.Therewith, the coprime sampling scheme is presented with sampling intervals r 0 T s and r 1 T s , as shown in Figure 1.Without loss of generality, it is assumed that r 0 < r 1 with r 0 , r 1 ∈ N + , where r 0 and r 1 are coprime.And the sampling interval T s corresponds to the Nyquist sampling rate f s .
For a wide-sense stationarity or cyclostationary process X(t), t ∈ R, it consists of a number of frequencies that are confined wihtin the bandwidth B s ≤ f s /2.The outputs of two uniform sub-Nyquist sampling channels can then be expressed as where x[n], n ∈ N denotes the Nyquist sampling samples, and the highest sampling rate of the coprime sampling system is given by 1/(r 0 T s ) = f s /r 0 .Accordingly, the elements of the sensing vector corresponding to the two coprime samplers can be donoted as and where l 0 = 0, 1, . . ., r 1 − 1, l 1 = 0, 1, . . ., r 0 − 1, and k = 0, 1, . . ., p − 1.
From a data acquisition perspective, the output samples obtained from the generalized coprime sampling scheme are a subset of the Nyquist samples, positioned at Based on the sensing vectors, the connection between the elements of the generalized coprime sampling vector and the Nyquist sampling vector can be expressed as and the elements of the connection sensing vector are defined as As a result, there is where T are all vectors of size N × 1.

III. PROPOSED FAST SOLUTION
Considering the widely-used unbiased estimation of the autocorrelation sequence for the output of generalized coprime sampling, the elements r y [m] can be given as Substituting equation ( 7) into equation ( 8) yields where |m| ≤ N −1.Thus, the power spectrum can be obtained by performing FFT on the autocorrelation sequence {r x [m]}, which is derived from the autocorrelation sequence {r y [m]} and {r a [m]}.Then, a practical solution that computationally efficient is employed to obtain the estimation of the autocorrelation sequence.The steps are as follows: Step 1: Pad vectors a[n] and y[n] with an additional Nlength zeros.
Step 2: Calculate the autocorrelation sequence based on the power spectrum estimation of vector a 2N T by involving FFT and IFFT.
Step 5: Obtain the power spectrum estimation by taking the FFT of vector rx The block diagram illustrating the fast power spectrum sensing solution for generalized coprime sampling is depicted in Figure 2. The proposed solution is efficient, as it only involves FFT/IFFT operations and some basic multiplication operations.Moreover, the autocorrelation sequence {r a [m]} of the sensing vector can be pre-calculated offline, as shown in the red section of Figure 2.This calculation is solely dependent on the generalized coprime sampling scheme.
Consequently, the computational complexity of the proposed solution involves the FFT of a (2N )-point sequence twice, resulting in (2N )log(2N ) floating-point operations according to (13).In addition, the FFT of a (2M − 1)-point sequence is executed once, leading to (2M − 1)log(2M − 1) floating-point operations according to (17).Upon incorporating 2M −1 multiplication calculations, the total computational complexity necessitates (4N )log(2N ) + (2M − 1)log(2M − 1) + (2M − 1) floating-point operations.This results in a lower computational complexity compared to the state-of-theart methods.Meanwhile, it is feasible to compute the FFT operations in parallel effectively, which is a more suitable practical solution to meet the real-time wideband power spectrum sensing applications.

IV. SIMULATION RESULTS
In the experiments, it is assume that there are I inputs with identical powers, which are distributed in the frequency band [2,18]GHz.Subsequently, the coprime integers r 0 = 3, r 1 = 4 and the Nyquist sampling rate f s = 32GHz are set.Furthermore, the relative root mean square error (RMSE) is adopted to evaluate the performance of the proposed fast power spectrum sensing method, which is defined as follows where fi (j) is the estimation of f i from the jth Monte Carlo trial, and five hundred Monte Carlo trials are conducted.
Herein, the estimated power spectrum results are first displayed in Figure 3  As shown in Figure 4, the RMSE results are compared as a function of the input SNR, where I = 18 MP signals are used and the frequency is randomly selected.It can be seen that the RMSE tends to be stabilized when SNR is greater than -2dB for p = 300.Under the same condition, the original method has the better performance at a low SNR by using the Toeplize matrix completion.Meanwhile, the estimation performance is improved as p increases, due to that the DOF increases with p and the resolution also improves.As expected, the Figure 5 has the same result.However, the selection of coprime sampling rate makes less difference to performance when p is greater than 1000.That is because the system redundancy under simulation is enough.Furthermore, the multiple coprime unit factor p not only affects the resolution, but also determines the execution time of algorithms.Thus, the execution time results are compared as a function of p as displayed in Figure 6, where I = 10 MP signals are used and the frequency is randomly selected with the 0dB SNR.It can be seen that the execution time of the proposed method has obvious advantages over the matrix completion method.Similar to the Figure 7, the proposed method outperforms the matrix completion method under the same execution time.Furthermore, the proposed method takes less time than the matrix completion method under the same performance.However, both methods rely on the larger number of samples.Therefore, there is a tradeoff between execution time and the system performance.
At last, the influence of the model matching degree between the sensing vector and measurements on the performance is discussed in Figure 8 with the 5dB SNR.Here, I = 18 frequencies are randomly selected for MP signals, and different time delay are used which is unknown for the sensing vector.As a result, the influence of the model mismatch has little fluctuation of RMSE within 200µs.That is interesting, which means that the fact potentially enables the spectrum resource in a distributed swarm scenario to be more efficiently utilized.

V. CONCLUSIONS
A fast power spectrum sensing solution for generalized coprime sampling is proposed, that only uses the parallel FFT and simple multiplication operations.It has obvious advantages over existing methods in terms of spectrum estimation performance and execution time, and there is no need to preestimate the number of inputs.Moreover, the influence of the model mismatch has little impact on the spectrum estimation performance, making it more suitable for further discussion on its application in the distributed swarm scenario.
Notations: The bold characters denote vectors.The notations R, N and N + represent the set of real numbers, nonnegative integers, and positive integers, respectively.The superscripts (•) T and (•) H indicate the transpose and conjugate transpose of a vector or a matrix, respectively.The operator • signifies the Hadamard product, | • | 2 signifies the element-wise square modulus of a vector, and ceil(•) signifies round up to an integer.The symbols F a and F −1 a mean the a-point fast Fourier transofrm (FFT) and the inverse fast Fourier transofrm (IFFT), respectively.arXiv:2311.13787v1[eess.SP] 23 Nov 2023 II.SIGNAL MODEL

Fig. 2 .
Fig. 2. Block diagram of the proposed fast power spectrum sensing solution.
, with p = 3000 and an input SNR of 15dB.There are I = 50 mono-frequency pulse (MP) signals randomly distributed in (a), I = 20 binary phase shift keying (BPSK) signals for 1M symbols per second with random frequency and code depicted in (b), I = 2 linear frequency modulation (LFM) signals with 10GHz bandwidth under ±6GHz initial frequencies given in (c) and I = 21 mixture signals of three types for (d).As can be observed, all frequencies are estimated accurately with the proposed method.