Adaptive and Blind Wideband Spectrum Sensing Scheme Using Singular Value Decomposition

The Modulated Wideband Converter (MWC) can provide a sub-Nyquist sampling for continuous analog signal and reconstruct the spectral support. However, the existing reconstruction algorithms need a priori information of sparsity order, are not self-adaptive for SNR, and are not fault tolerant enough. These problems affect the reconstruction performance in practical sensing scenarios. In this paper, an Adaptive and Blind Reduced MMV (Multiple Measurement Vectors) Boost (ABRMB) scheme based on singular value decomposition (SVD) for wideband spectrum sensing is proposed. Firstly, the characteristics of singular values of signals are used to estimate the noise intensity and sparsity order, and an adaptive decision threshold can be determined. Secondly, optimal neighborhood selection strategy is employed to improve the fault tolerance in the solver of ABRMB. The experimental results demonstrate that, compared with ReMBo (Reduce MMV and Boost) and RPMB (Randomly Projecting MMV and Boost), ABRMB can significantly improve the success rate of reconstruction without the need to know noise intensity and sparsity order and can achieve high probability of reconstruction with fewer sampling channels, lower minimum sampling rate, and lower approximation error of the potential of spectral support.


Introduction
Spectrum resource has become increasingly scarce with emerging wireless services.Nevertheless, assigned radio spectrum to authorized users is mostly underutilized.Cognitive radio (CR) technology can improve frequency spectrum utilization by detecting and accessing the frequency range that has not been occupied by authorized users (primary users).Therefore, one of the crucial tasks in the cognitive radio system is to constantly monitor the potential spectrum bands and detect the activities of primary users [1,2].The methods of signal detection on narrow band mainly include energy detection [3], coherent detection [4], and feature detection [5].Since the wideband spectrum sensing can provide more spectrum access opportunities for cognitive users (secondary users), multiband based wideband spectrum sensing has gained much more research attention in recent years [6,7].
At the receiving side, the traditional method of acquiring RF (Radio Frequency) information is to use demodulation under the condition of known carrier.However, in practical scenarios, it is often required to directly perform blind sensing for high frequency wideband analog signals, which brings great challenges to spectrum detection.Firstly, extremely high sampling rate has exceeded the limit of physical ADC (Analog-to-Digital Converter) capability.Meanwhile, the storage and transmission for the sampled data can bring huge overhead.In addition, the carrier frequency of the received multiband RF signal is usually unknown.To address these challenges, compressed sensing theory is applied to wideband spectrum sensing based on compressed sampling and signal reconstruction [8,9].In order to perform sub-Nyquist sampling on continuous spectrum signal, Mishali et al. proposed the MWC (Modulated Wideband Converter) scheme which uses multichannel compressed sampling and reconstruction with parallel structure based on compressive sensing and time-frequency transform theory [10,11].In MWC system, the multiband analog sparse signal can be collected with sub-Nyquist sampling, and the spectral support of signal can be reconstructed by a CTF (continuous-to-finite) block.

Wireless Communications and Mobile Computing
From [10], it is known that the time domain reconstruction model of MWC can be attributed to the Multiple Measurement Vectors (MMV) problem.MWC has been considered to be an effective approach for compressed sampling and signal reconstruction with its applications on radar, broadband communication, and cognitive radio spectrum sensing [12][13][14].
The accurate reconstruction of spectral support of signal is the key to realize the spectrum sensing.At present, there is still much room for improvement in the success rate of sensing, the required minimum number of channels, and the maximum number of subbands that can be reconstructed.The solvability of MMV problems is important for MWC reconstruction ability.To enhance the solvability of MMV problem, Reduce MMV and Boost (ReMBo) algorithm, which can transform MMV problem into single measurement vectors (SMV) problem is proposed in [15].In ReMBo, the sampling value matrix  is projected on a random vector  which is independently and identically distributed in the interval [−1, 1] and obtains a column vector V, V = .Then, it tries to obtain the estimated spectral support of signal by iteratively computing SMV problem (V = Φ = Φ) for many times.Although ReMBo is better than most of the existing reconstruction algorithms in the success rate of reconstruction and computational overhead, there is still a big gap between the maximum number of reconstructed frequency subbands and the theoretical upper limit [16].The total sampling rate and the minimum number of required channels of ReMBo algorithm are still much higher than the theoretical lower bound [10,17].To improve the performance of MWC system, a dimension adjustable Randomly Projecting MMV and Boost (RPMB) framework algorithm, which transforms the initial MMV problem into a series of low dimensional MMV problems with the same sparsity order, was proposed in [18].In the solving process, the solution is tested several times until a satisfactory solution is obtained.Compared with ReMBo, RPMB can reduce the number of required hardware channels and improve the maximum number of reconstructed frequency subbands in the condition of precise reconstruction.In spite of this, the RPMB framework is still faced with difficult problems to be solved: (1) The reconstruction performance of RPMB is greatly reduced under low SNR because the influence of the noise intensity is not considered in the RPMB solver; (2) the RPMB is not flexible enough since the decision threshold of RPMB is a fixed value, which has no self-adaptive ability for the uncertainty of the noise; (3) RPMB needs to know the sparsity order of the signal in advance, which is difficult in some practical applications; (4) although the RPMB solver improves the algorithm's fault tolerance by finding more potential support bands, the effect of fault tolerance remains to be improved.
To address the above-mentioned problems, this paper proposes a Self-Adaptive and Blind Reduced MMV and Boost (ABRMB) scheme for MWC based on SVD [19].Our main contributions are summarized as follows.
(1) ABRMB can first reveal the intrinsic properties of signal using SVD and perform the linear fitting on the noise singular value by using the tail singular values as noise data.
Then, it can estimate the noise intensity, and the threshold of deciding support elements can be set adaptively.
(2) Our proposed scheme calculates the contribution of the estimated noise singular on the singular values of signal and figures out the sparsity order by using gradient and difference operation.Hence the number of signal subbands can be estimated.
(3) In the solver of ABRMB, the fault tolerance can be improved by the optimal neighborhood selection method.Compared with ReMBo and RPMB, when noise intensity and sparsity order are unknown, our proposed scheme can use fewer channels to achieve the high probability reconstruction of spectral support in low SNR range and increase the maximum subband number of reconstructed signal.
The rest of this paper is organized as follows.In Section 2, the related work is introduced.In Section 3, the multiband sparse signal model is introduced, and sampling principles of MWC are described as well as some of the present major issues.Section 4 introduces the key ideas of the noise intensity estimation, sparsity order estimation, and optimal neighborhood selection.In Section 5, the ABRMB scheme includes its solver which are specifically described, the convergence proof of which is also given.Section 6 performs validation and analysis on ABRMB from four important performance evaluation indexes.Conclusions and discussion are in Section 7.
The main notations used in this paper are listed in Notations to make them clearer and easier to read.

Related Work
Current research on the performance of spectrum sensing can be mainly classified into the following two categories.The first category, especially in the design of CRN, focuses on investigating the optimal trade-off between energy consumption and throughput for secondary users (SUs).The second category, especially for wideband spectrum sensing, mainly focuses on studying how to achieve a good balance between noise interference and sensing accuracy.
As for the first category, some earlier researches focus on the trade-off between the interference and throughput.Specifically, for the unslotted CRN, Yang et al. [20] designed an interference-constrained novel simultaneous sensing and transmission scheme.This scheme exploited the statistic information of the activities of primary users (PUs), and the transmission duration is adaptively adjusted to avoid the interference.In the same period, energy efficient techniques, in cooperative spectrum sensing (CSS), have also gained a lot of attention [6].In [21], considering the energy consumption problems of each sensing node in cooperative CR sensor network, a user selection scheme was proposed to minimize the energy consumed on each CR node based on the binary knapsack problem and its dynamic solution.After that, Ejaz et al. [22] indicated that an optimal sensing, reporting, and transmission duration can be found while providing a given throughput and reliability constraint.The solution can obtain the best trade-off between energy consumption and throughput for SUs.As for the second category, in order to overcome the limitation of hardware technology and achieve wideband spectrum sensing, the research work mainly concentrated on how to get a good trade-off between the noise interference and sensing accuracy [4].Firstly, based on the compressed sensing theory, a parallel MWC architecture for wideband spectrum sensing was designed by Mishali et al. [10,11], which can achieve sub-Nyquist sampling for the wideband signal.Then, a MWC system with SwSOMP algorithm was proposed in [23].Given the noise interference, the algorithm can get higher reconstruction precision of spectral support.In [24], in order to reduce the high hardware complexity of MWC, a compressive wideband spectrum sensing scheme based on a single channel was designed.
The work of this paper belongs to the second category.All of the aforementioned literatures only consider ideal noise interference model (NIM) and assume that the signal sparsity is known.However, in reality, NIM is always nonideal, where the intensity of noise varies with the change of communication environment.Moreover, the signal sparsity is always unknown in reality.In addition, we also hope to adopt a better reconstruction strategy and use fewer parallel channels to obtain higher reconstruction accuracy.Therefore, these problems need to be further investigated.

Problem Statement of MWC Model and Our
Proposed Problem Solution 3.1.MWC System Model.The reconstruction principles of the spectral support of MWC system is shown in Figure 1. Figure 1(a) shows the sparse multiband analog signal model.Sparse multiband analog signal is common in cognitive radio [25].Suppose the received signal () is a sparse bandpass analog signal.The spectrum is distributed in frequency interval [− nyq /2,  nyq /2], and  nyq is Nyquist rate, which usually can reach GSPS order of magnitude.Suppose the spectrum of () only contains  subbands whose bandwidth is  푖 ≤  ( ≥  > 0), and there is no overlap between subbands. is the maximal bandwidth in subbands, and the center carrier frequency of each subband is unknown. 푁 and  can be defined as In ( 1),  푁 is union of all subbands whose amplitude is not zero, which represents the effective frequency component of ().
As shown in Figure 1(a), the whole band is equally divided into  consecutive narrow band channels, every bandwidth of which is no less than .Then the spectrum of () has at most 2 parts which get energy in the whole frequency band.If the channel is marked [1, . . ., ], the set of channel numbers corresponding to each subband  푖 () is the spectral support of (), which is defined as Λ = supp(()).|Λ| is the number of nonzero elements in Λ, that is, the potential of Λ.The frequency bands corresponding to these numbers are called support bands.Since 2 is much smaller than , () can be considered as sparse multiband signal.Figure 1(b) shows the reconstruction process of the support of ().Assume that the number of bands is 4,  푠 =  푝 , and  푝 ≥ .We divide the wideband spectrum into  spectrum slices, where  = 2 0 + 1.In order to ensure that the discrete Fourier transformation result of sampling sequence contains all the components of the original signal spectrum (),  0 must satisfy  0 = ⌊( nyq +  푠 )/2 푝 ⌋ − 1.After mixing and lowpass filtering, the spectrum information of the original signal appears in the sampling interval [− 푠 /2,  푠 /2], and the mixing Wireless Communications and Mobile Computing coefficient of each spectrum slice is  푖푙 , where  is the index of spectrum slices.According to the theory of compressed sensing, we can obtain the spectrum support of the multiband signal.
In conclusion, the support bands of () must meet two conditions: (1) it is distributed in an extremely wide frequency range; (2) the support bands of signal only exist in a few discrete frequency spectrum.

Sampling in MWC System
. MWC contains multiple parallel sampling channels, and each channel has same structure with mixer, low-pass filter, and ADC [10].The received signal () is input to  parallel channels at the same time, and each channel is multiplied by periodic mixing signal  푖 () with different mode, which makes spectrum of () move to baseband.Each of the channels'  푖 () is not related.The period of  푖 () is  푝 = 1/ 푝 , and  is used to show the number of random ±1 switches in a cycle. 푝 is defined as the switching frequency of mixed signals.The mixed signals pass through the low-pass filter whose cut-off frequency is 1/2 푠 .At last, it passes through ADC whose sampling rate is  푠 = 1/ 푠 and obtains  groups low-speed digital sampling sequences  푖 [].
On analysis of the th channel, the Fourier series expansion of the random mixing function is where When  = 0,  0 = 1/, and when  ̸ = 0,  푙 = (1 −  −푗(2휋푙/퐿) )/2.In (2),  푖 () denotes a pseudorandom sequence of ±1, which is used as a mixing signal of the th sampling channel;  is index of spectrum slice;  푖푙 is the coefficient of Fourier series expansion.
Next, after passing through the low-pass filter with frequency characteristic of () = { 1 |푓|≤푓  /2 0 |푓|>푓  /2 , the relationship between DTFT (Discrete Time Fourier Transform) of  푖 [] and the Fourier transform () of () is (3) and  0 is the smallest integer that can satisfy  = 2 0 + 1 ≥  nyq / 푝 .Formula (3) shows that the spectrum of the output sequence  푖 [] is equivalent to the shift weighted sum of the original signal spectrum () with  푝 as its step, which is intercepted into spectral fragments with  푠 width by low-pass filter.If  푖 ( 푗2휋푓푇  ) is considered as the th component of m-dimensional column vector () and ( −  푝 ) is considered as the th component of 2 0 + 1dimensional column vector (), then (3) can be expressed as In ( 4), Φ is a  ×  measurement matrix, For any frequency  ∈ [− 푠 /2,  푠 /2], (5) is a typical compressed sensing problem.When sampled values matrix  and an observation matrix Φ are known, a sparse vector  can be recovered.Since MWC is a MMV problem, the one-dimensional measured value vector turns into a twodimensional matrix.As a result, paper [15] gets a onedimensional vector after projecting  ingeniously and then reconstructs the signal support bands by using compressed sensing technique.Furthermore, the reconstruction by building a CTF (continuous-to-finite) module is described in [10].

Problem Statement.
The principle of reconstructing spectral support for MWC is shown in Figure 1 where Figure 1(b) gives a specific reconstruction process.As can be seen in Figure 1, if the number of support bands of () is , the maximum reconstructed subbands are 2 after passing through the MWC system; that is, the maximum possible sparsity order  is 2.In the MOMPMMV (Modified Orthogonal Matching Pursuit for MMV), the solver mentioned in paper [18], the initial solution Z is obtained via the OMPMMV, which needs to know  beforehand.In addition, in the condition judgment of RPMB, the algorithm uses | Λ| ≤  to determine whether the potential of spectral support meets the requirements.However, in cognitive radio applications, it is difficult to know the sparsity order of signal.Therefore, aiming at the practical application, an effective method to estimate the sparsity order of signal is needed.
Secondly, in [18], the condition of deciding spectral support is ‖ Z푖 → ‖ 2 ≥ ; and the condition of determining whether the spectral support meets the requirements in the iterative process of RPMB is ‖ Ŷ − Φ Ẑ‖ 퐹 ≤ .As mentioned in Section 1, in the case where there is no noise or high SNR,  can be figured by taking a fixed small value.If the SNR is low or fluctuates greatly, the value of  will have to perform self-adaptive change according to the intensity of the noise.In paper [18],  is a fixed value without considering SNR.Hence, RPMB algorithm is neither suitable nor flexible under low or fluctuated SNR.Thus an effective method is needed to estimate the intensity of the noise and adaptively change the decision threshold .
Finally, in MOMPMMV, to improve the fault tolerance ability, the method is Λ = max(‖ Z 푖→ ‖ 2 , ), where  is the number of sampling channels, that is, taking  maximum norm 2 as the initial spectral support.This method is effective, but it is not the best way.For example, as shown in Figure 1(a), the /2th subband  푁/2 () of () is located at the 14th and 15th narrow channels, and the spectrum energy of the 15th narrowband channel is quite small, as a result of which if the fault tolerance method in [18] is used, under low SNR, number 15 will not be able to merge into the initial support Λ.So it is necessary to design a better fault tolerance strategy.

Brief Description of Our Proposed Problem Solution.
Our proposed problem solution is shown in Figure 2 where multiband analog signal () passes through the MWC system to obtain compressed sampling data .Then, in order to reduce complexity, we reduce the dimension of  by means of random projection.Subsequently, the signal sparsity and noise intensity are estimated in the preprocessing, namely, K and  Noise .Meanwhile, to improve the fault tolerance ability, we adopt the optimal neighborhood selection strategy, which is described in detail in Section 4.3.Finally, by using the prior knowledge and the compression sensing reconstruction algorithm, we can obtain the estimated spectral support Λ.It is worth noting that, in the reconstruction algorithm, we can adaptively adjust the decision threshold by using the estimated noise intensity  Noise .

Preprocessing of ABRMB Scheme
4.1.Estimation of Noise Intensity.Singular value decomposition is one of the most basic and important tools in modern numerical analysis.The SVD is performed on the measured values  with noise using (6). ∈ R 푚×푟 ,  ≥  > 0,  ≥  > 0, and  < , where  is sampling length.
In (6), unitary matrix  ∈ R 푚×푚 is the left singular vector of , and unitary matrix  ∈ R 푟×푟 is right singular vector.Σ ∈ R 푚×푟 is a diagonal matrix, and the main diagonal element is the singular value with a descending order.Since rank() = , (6) can be simplified as Let Σ 퐾 = diag( 1 ,  2 , . . .,  퐾 ) and Σ 푚−퐾 = diag( 퐾+1 ,  퐾+3 , . . .,  푚 ). 푖 is th singular value of , and  ≥  ≥ 1.After SVD, an important characteristic is that most of the energy of signal is concentrated in the first  singular values, while the energy of the noise is distributed in all of the singular values; but noise can be reflected in the tail singular value.Assume that the singular value of  can be decomposed into the singular value Σ 푠 of the original signal and the singular value Σ 푛 of the noise, and Σ 푠 = diag( 푠1 ,  푠2 , . . .,  푠푚 ), Σ 푛 = diag( 푛1 ,  푛2 , . . .,  푛푚 ). Figure 3 describes the contribution of Σ 푠 and Σ 푛 to Σ under different SNR, which are represented by  푠 and  푛 respectively.The definitions are as follows: The contribution of the original signal to the tail of the singular values is small.The energy of the signal is concentrated in the first  singular values, and the tail of the singular values is mainly determined by noise.Figure 3 shows that the best data source for noise intensity estimation is bottom 50% of the singular values.Because the available singular values are also reduced with the decrease of the sampling channels, in order to ensure the accuracy of estimation, the bottom 30% of the singular values is adopted to estimate the noise intensity in this paper.
Figure 4 illustrates the variation of the singular values of  with different noise intensity and different sampling channel number.As shown in Figure 4, the singular value is larger when the SNR is lower or the sampling channels number is more.It should be noted that since the initial part of the singular value is mainly determined by the energy of the signal, the initial part of the singular value has little change versus different SNR and channel numbers.Yet the tail of the singular values changes a lot under different SNR because it is influenced significantly by noise intensity as the contribution of real signal is much less for those values compared with noise.This fact is utilized in this paper to estimate the noise intensity.
Next, the distribution of the noise singular values is investigated.Since the noise is the uniformly distributed white Gaussian noise (WGN), it can be seen in Figure 5, the singular values of noise are largely distributed on a straight line under the circumstances with different SNR and different sampling channel number.Using this fact, the singular values of tail 30% of  are used as singular values of noise to perform linear fitting.The fitting line is used to estimate the noise singular values, and the noise intensity is estimated by (9).The fitting result is shown in Figure 6.Since there is little The value of The value of energy of real signal in the tail of the singular values, there is a certain deviation between the fitted singular values and original singular values of noise.However, under the low SNR, the deviation will be diluted due to the relatively large noise intensity.In addition, in order to reduce the influence of the deviation, this paper uses the optimal neighborhood selection strategy in Section 3.3 to carry out fault tolerance processing.Suppose Nosie ∈ R 푚×푟 is a uniformly distributed WGN, then the intensity  Nosie of the noise can be calculated according to the formula where  푛 = { 푛 (1),  푛 (2), . . .,  푛 (), . . . 푛 ()};  푛 () denotes the th noise singular value.

Estimation of Sparsity Order.
According to the previous analysis, we can get the estimated noise singular values Σ푛 = diag(σ 푛1 , σ푛2 , . . ., σ푛푚 ).The contribution of noise singular values on the singular values of  is calculated by If the noise intensity is strong,  푅  is first obtained by performing gradient operation on R푛 , and then  푅  is obtained by performing difference operation on  푅  .The results of the operation are in an ascending order.The position of the minimum value in  푅  plus 1 is the estimated sparsity order K.If the sampling channel number is close to the theoretical lower limit, K need to add one adjustable parameter , and empirical value  is 1.The calculation methods of  푅  and  푅  are shown in (11).Figure 7 is a sketch of estimation of sparsity order under different SNR.
When noise is low, the signal energy is dominant.The sparsity order can be estimated directly by using the singular values of .First, all singular values are shifted to the left one time, and the last empty position is filled by  푚 , represented as Σ 푎 = diag( 2 ,  3 , . . .,  푚 ,  푚 ).Then,  is calculated by (12) and is listed in a descending order.The index of maximum value in  is the estimated sparsity order K.Although K has some deviations, the deviation has little impact on the success rate of reconstruction due to the joint sparse reconstruction and fault tolerance mechanism.

Optimal Neighborhood Selection Strategy.
From the signal model showed in Figure 1(a), we know that the energy of subband can only locate on two adjacent narrow channels.If support Λ 푖1 is obtained, support Λ 푖2 can only be its neighborhood channel; that is, Λ 푖2 ∈ nei(Λ 푖1 ).However, there are usually two neighborhood channels of Λ 푖1 .Obviously, the channel with the largest norm 2 value is chosen to be the other support channel, as given in where the function of pos() is to obtain the optimal neighborhood support of Λ 푖1 .‖ nei(Λ 1 )→ ‖ 2 denotes the signal energy on two neighborhood channels of Λ 푖1 .Obviously, this strategy is very helpful to obtain the approximate spectral support.

Proposed ABRMB Optimization Scheme
5.1.Structure of ABRMB Scheme.It was pointed out in [18,26] that random projection can be used for joint reconstruction to improve reconstruction performance.In particular, it can be known from [18] that the sensing performance tends to be stable if the number of preserved measured vectors is equal to  after the operation of projection.The ABRMB also inherits this idea, and it is combined with the pretreatment methods mentioned in Section 3 to improve the integrated system performance.The structure of ABRMB is shown in Figure 8.

Pseudocode of ABRMB Scheme
. Pseudocode 1 provides the pseudocode to describe the ABRMB scheme. represents a uniformly distributed random matrix whose values are continuously selected from [−1, 1].
As can be seen in Pseudocode 1, the ABRMB allows a certain degree approximation deviation for spectral support with the aid of estimated noise intensity and sparsity order.In addition, if the loop iteration is ended, the optimal spectral support is still not found; the ABRMB can select the spectral support, with minimum potential, as the optimal support set from all of the stored spectral support.If there are more than one minimum potential, then choose the one with smallest residual as the best spectral support.Because of the guarantee of fault tolerance mechanism, the experiment shows that this method is very effective.
The solver of ABRMB scheme is described in Pseudocode 2.

The Convergence of ABRMB Scheme
Theorem 1. Assume that () is multiband signal described as in Figure 1(a); we use the MWC structure for sample signal as shown in Figure 8.If the following conditions are established, for any  ∈ [− 푠 /2,  푠 /2], () is the only N-sparse solution of (4).

𝑚 ≥ 2𝑁.
3. The number  of ±1 symbols in a periodic sequence 4. Any 2 column of Φ is linearly independent.
The proof of Theorem 1 is shown in [10].
Next, we prove that the spectral support of the original MMV problem Λ = supp() can be obtained by solving the new MMV problem  = Φ, where  is the result after projection.
Lemma 2. Suppose  ∈  퐿 is a known nonzero vector, and  ∈  퐿 is a random set of vectors obtained from a continuous probability distribution.Then, the probability that the event " on  is not 0" is 1.
In addition, MWC can achieve sub-Nyquist sampling.Without affecting the success rate of reconstruction, the subsequent processing of the system is more favorable with higher degree of compression.The total sampling rate of MWC is  Σ =  푠 .The theoretical minimum of sampling rate for multiband signal, that is, the Landau rate [15], is defined as where ( 푁 ) is Landau rate; that is, it is the sum of all subband frequency widths,  푁 denotes the union of  effective frequency components in multiband signal, and ( 푖 −  푖 ) is frequency width of th subband.
Since the number of the channels  and total sampling rate  Σ are directly relevant, obviously, the cost of the system and the corresponding sampling rate are lower for smaller .

The Maximum Number of Subbands That Can Be
Reconstructed in Signal.From [16], the upper bound of the reconstruction capability of the MMV problem is In (15), () is the joint sparsity order of , (Φ) is the Kruskal rank of Φ, and Rank() is the rank of .As can be learned from (15), the smaller the rank , the lower the sparsity order of the signal that can be reconstructed, which is also the reason that the upper bound of the sparsity order in [15] is particularly small.The number of subbands is  = ()/2.Therefore, the performance of the system is better when the number of subbands that can be reconstructed is larger.

The
In (16), | Λ| is the potential of estimated support, (| Λ| − |Λ|) represents the difference between the estimated support potential and the actual support potential, and "s.t." is the abbreviation of "subject to".In (17), |Λ| max is the potential of maximum spectral support for multiband sparse signals. is the number of narrow bands, and  is the number of support bands of the original signal.Apparently, smaller approximate error is better.Since  ≫ |Λ|, in the application of wideband spectrum sensing in cognitive radio, as long as  Λ ≤  upper , the impact on the secondary users can be negligible.

The Successful Probability of the Reconstruction.
In the analysis of recovered spectral support, we refer to the successful recovery criteria in [10]; that is, when the estimated support Λ and the actual support Λ satisfy (18), where Λ ⊇ Λ, and Φ ↓ Λ is with full column rank, it is considered as a successful reconstruction.
In (19),  푖 ,  푖 ,  푖 , and  푖 represent the energy coefficient, bandwidth, carrier frequency, and time offset of the th band, respectively. is the number of subbands in the signal.() is white Gaussian noise.The following procedure is repeated 500 times to calculate the probability of success.4. Estimate the spectral support using ReMBo, RPMB, and ABRMB, respectively, and determine whether it is successfully recovered.
Firstly, under the same conditions, the changes of the success rate of reconstruction with channel number using ABRMB, RPMB, and ReMBo are studied under different SNR.As can be seen in Figure 9, when SNR equals 10, 20, and 30, respectively, the success rate of the reconstruction of spectral support using ABRMB is better than using ReMBo or RPMB.Particularly, when  = 20 and SNR = {10, 20, 30} dB, the improvements of the success probability of ABRMB compared with RPMB and ReMBo are shown in Table 1.Therefore, ABRMB scheme is more effective and has better adaptability in the case when SNR is low or fluctuated due to the reasonable estimation of the noise intensity and the sparsity order and the optimal neighborhood selection strategy.
Secondly, in the same circumstances, the needed minimum channel number and the minimum sampling rate of high probability reconstruction are studied.As shown in Figure 9 and Table 2, under different SNR, the number of hardware channels and the sampling rate of ABRMB for the   the system can be saved as the number of hardware channels needed by the scheme which can be reduced.
Next, under the same conditions, the changes of success rate of reconstruction with the number of the bands are studied using these three algorithms.The number of the subbands  is chosen in the range [2,12] with 2 as its advancing step.The parameters are SNR = 20 dB,  = 20,  푖 ∈ {1, 2, 3, 4, 5, 6}, and  푖 ∈ {0.4,0.7, 0.2, 0.9, 1.2, 1.5} s.As shown in Figure 10, since the number of the bands is directly related to the sparsity order of the signal, with the increase of , the signal is not sparse enough and the success rate of reconstruction for three algorithms is reduced.Nevertheless, the reconstruction performance of ABRMB is obviously better than that of RPMB and ReMBo for  = 4, 6, and 8.
At last, we compare the average approximation error on the potential of estimated spectral support.We set  = 6, SNR = 20 dB, and  is in the range [15,40] with 1 as its step.As can be seen in Figure 11, since the potential of the support reconstructed by RPMB is related to the number of channels, when  = 24, the average approximation error of the spectral support of RPMB has already reached the upper limit.From Figure 11, the average approximation error of ABRMB is the smallest in these three algorithms, and the potential of the spectral support reconstructed by ABRMB is the closest to the actual number of frequency channels.In this way, the  ABRMB can provide more spectrum access opportunities for the secondary users in cognitive radio networks.
From the analysis of above four performance metrics, we can see that ABRMB scheme can achieve low-speed sampling by utilizing the sparse characteristics of wideband signals in frequency domain.The sampling rate can be reduced to 14.9% of the Nyquist sampling rate when SNR = 10 dB.Meanwhile, for noise uncertainty, the detection threshold can be adaptively changed as the noise power changes in ABRMB scheme.Based on the neighborhood selection strategy and sparsity estimation, ABRMB can obtain higher success reconstruction probability.Then, the projection reduction operation can reduce the computational complexity, and it does not affect the convergence of the scheme.Finally, ABRMB scheme can find better trade-off between noise interference and sensing accuracy.

Conclusions
By using SVD theory, this paper proposes a self-adaptive and blind wideband spectrum sensing MWC scheme which leads to a flexible and high performance solver.The SVD is performed before signal reconstruction in the preprocessing block.The noise intensity and sparsity order of the signal are estimated in this block, and then the subbands number of the signal can be obtained.ABRMB scheme can use the estimated noise intensity and sparsity order to process the multiband

15 cFigure 1 :
Figure 1: The reconstruction principles of the spectral support of MWC system.

Figure 2 :
Figure 2: The flowchart of our proposed problem solution.

iFigure 3 :
Figure 3: The contribution of the singular values of the signal and noise to total singular values.

Figure 4 :Figure 5 :
Figure 4: The comparison of singular value of  under different SNR.

Figure 6 :i 12 − 20 Figure 7 :
Figure 6: The linear fitting of the noise singular values.

Figure 8 :
Figure 8: The support reconstruction scheme based on ABRMB.

Figure 9 :
Figure 9: The comparison of recovery success rate of support under different SNR.

Figure 10 :
Figure 10: The effects of subbands number on the support recovery.

Figure 11 :
Figure 11: The comparison of the approximation error of the support potential.
is performed on both ends of (4), we can get the corresponding relationship between the sequence [] = [ 1 [],  2 [], . . .,  퐿 []] 푇 and the sampling data 1 ≤  ≤ , and  < .If IDTFT (Inverse Discrete Time Fourier Transform) Approximate Error of the Support Potential.For the original signal, the potential of spectral support is |Λ|, that is, the length of spectral support.Obviously, |Λ| = ‖Λ‖ 0 .If | Λ| = |Λ|, then the approximate error is 0. If | Λ| > |Λ|, then the approximate error existed.In this paper, the approximate error  Λ and the upper bound of the error  upper are defined as
3. Generate new sinc signal according to  푖 .

Table 1 :
Comparison of maximum promoting rate of the reconstruction success probability.
high probability reconstruction are smaller than those of the ReMBo and RPMB.Thus, our proposed scheme can use fewer hardware channels and lower sampling rate to achieve the high success rate of reconstruction.Obviously, the cost of

Table 2 :
The comparison of the minimum channel and minimum sampling rate needed.