A Novel Approach to Wideband Spectrum Compressive Sensing Based on DST for Frequency Availability in LEO Mobile Satellite Systems

In LEOmobile satellite network, the L/S frequency availability is an essential task for global communication but entails severalmajor technical challenges: high sampling rate required forwideband sensing, limited power and computing resources for processing load, and frequency-selective wireless fading. This paper investigates the issue of frequency availability in LEO mobile satellite system, and a novel wideband spectrum compressed signal detection approach is proposed to obtain active primary users (PUs) subbands and their locations that should be avoided during frequency allocation.We define the novel wideband spectrum compressed sensing method based on discrete sine transform (DST-WSCS), which significantly improves the performance of spectrum detection and recovery accuracy compared with conventional discrete Fourier transform based wideband spectrum compressed sensing scheme (DFT-WSCS). Additionally, with the help of intersatellite links (ISL), the scheme ofmultiple satellites cooperative sensing according to OR and MAJ decision fusion rules is presented to achieve spatial diversity against wireless fading. Finally, in-depth numerical simulations are performed to demonstrate the performance of the proposed scheme in aspect of signal detection probability, reconstruction precision, processing time, and so forth.


Introduction
Currently, there is a renewed R&D interest and market opportunities for mobile satellite systems (MSS), which have been utilized for providing communication services to fixed and mobile users for a vast range of sectors (e.g., land-mobile, aeronautical, maritime, transports, rescue and disaster relief, and military) where terrestrial infrastructures are congested or infeasible to be deployed [1].According to the orbit types, the operational MSS can be categorized into GEO based MSS and non-GEO based MSS.In those systems, the satellite antenna has multiple spot beams that irradiate cells on the earth, thus creating a cellular-like coverage.Each orbit has some desirable properties as well as a set of drawbacks for efficient service delivery.However, it is impossible to be equipped with large antennas and high power supply for mobile terminals that need to communicate with each other around the world, such as handsets; the idea of employing lower altitude satellites is put into practice.Details of the satellite orbit, user link frequency band, physical layer, multiple access, satellite features, and intersatellite links (ISL) can be found in [2][3][4][5].
Frequency for MSS is assigned at the World Radiocommunication Conferences (WRC) periodically organized by the International Telecommunication Union-Radiocommunication (ITU-R) sector.While fixed services use high C and K frequency band, mobile services are better suited for L and S frequency band allocated at the World Administrative Radio Conference (WARC) 92 [6].The limitation of L/S band resources and the rapid growth demand for broadband data services are a push toward the use of Ku/Ka band for MSS.ITU-R has assigned portion of Ka frequency band to MSS and fixed satellite systems (FSS) on a primary basis in all regions (29.9-30GHz for earth-to-space link and 20.1-21.3GHz for space-to-earth link) and Ku frequency band to MSS on a secondary basis in all regions (14-14.5 GHz for earth-tospace link and 10-12 GHz for space-to-earth link).In fact, the systems adopted Ku/Ka are mainly intended for fixed users with big antennas for higher signal attenuation and impact of atmospheric effects and cannot provide a good coverage over seas, because antenna spot-beams footprints are focused on landmass [7].L/S frequency band has been exploited by MSS for a long time due to permission for small on-board antennas and personal communication anywhere with anyone at anytime.In other words, it is obvious that almost all the operational MSS regard L/S band as the preferential frequency resource to support mobile users.
According to the final documentation of WRC held in Geneva, Switzerland, 2012 (WRC-12), the allotment of L/S band for mobile satellite services promulgated by ITU Radio Regulation (ITU-RR) is listed in Table 1.In order to reuse frequency and alleviate interference, ITU-RR divides the earth into three regions: Region I (e.g., Europe, Africa, Soviet Russia situated in Asia, and Asia countries lying in the west of Mongolia and Iran's western borderline), Region II (e.g., Latin America, North America, Greenland, and Hawaii), and Region III (e.g., Oceania and the remaining district of Asia).
The statistical results of L/S band used by operational systems obviously show that it is impossible to find a completely pure L/S subband for another planned global mobile satellite system, which has gained much attention in those countries, for example, China, Russia, and European Union, that have not owned such an independently controllable system.The scarcity of L/S frequency resources becomes one of the crucial reasons that some MSS have to be postponed, such as LEO-based Courier panned by Germany, European Union, and Russia and the 48 polar-orbiting satellite constellation system proposed by China.Inevitably, whatever strategy of frequency availability for any mobile satellite system should consider the frequency allotment and region division for mobile satellite service in a global context and avoid interference with the operational MSS.
In general, the space segment satellites from different countries adopt different signal frequency when providing satellite mobile service to ground segment terminal users.The analysis of each L/S subband assigned to the operational MSS implies that it is not absolutely impossible to actualize user link frequency for those forthcoming systems on the basis of discovering the available vacancy spectrum.Concretely, it can be switched to subband 1668-1675 MHz/1518-1525 MHz when the satellite moves out of the coverage area of AlphaSat I-XL.Contrarily, in the area of Europe, Africa, and Middle East, it can be changed to subband 1980-2010 MHz/2170-2200 MHz.For high latitude area (i.e., outside ±70 ∘ latitude and polar region), it is manipulable to choose subband 1610-1626.5/2483.5-2500MHz because Globalstar provides communication services in mid-low latitude zone.Obviously, performing signal detecting for identification of active spectrum becomes the first and the key task to obtain frequency autonomously.In this aspect, it is similar to cognitive radio (CR) which has been found as a competitive solution to the issue of underutilization of the allocated spectrum [8].The purpose of spectrum sensing technology in CR is to discover the vacant spectrum to actualize utilization of the unoccupied spectrum by cognitive radio users, and the objective for application of spectrum sensing in MSS is to seek available frequency to provide mobile communication services.
In cognitive radio networks (CRN), secondary users (SUs) need sense the wideband spectrum fleetly and accurately in order to utilize the unused spectrum, as well as avoiding interference to PUs.The methodology significantly improves the use efficiency of the limited radio frequency resource as a result of a low percentage of the whole wideband spectrum taken by PUs.Recently, wideband spectrum sensing represents an effective solution to the issue of simultaneous sensing of the spectrum occupancy of multiple signals spreading over a wide frequency range [9].However, in wireless network, including satellite network, current analogto-digital hardware technology faces tremendous challenge posed by conventional spectral estimation method which operates in a way equal to or greater than the Nyquist sampling rate.Besides, time-varying spectrum sensing only collects a small quantity of measurements for detecting, which makes it extremely difficult to fulfill high-resolution signal reconstruction reliably.
With the recent advances in the theory of compressed sensing (CS), signal sparsity has been applied to reduce the number of measurements while maintaining high recovery accuracy [10,11].Compressed sampling (i.e., combined sampling with compression operations) as the efficient method to achieve real-time spectrum estimation at a sub-Nyquist sampling rate has attracted substantial research efforts, which is motivated by the fact that there are plenty of temporarily unemployed spectrum holes in the wideband scenario.The compressed sampling theory points out that any sparse signal can be compressed and sampled simultaneously.Preliminary work has been done to exploit CS for wideband spectrum compressed sensing based upon signal sparsity in frequency domain.[12] firstly carried out the work to take advantage of CS theory to acquire wideband signals using sub-Nyquist sampling rate.Recently, [13] presented a structured compressed sensing.And [14] proposed a power spectrum blind sampling (PSBS) algorithm trying to reconstruct the power spectrum.However, all the aforementioned methods belonged to the category of discrete Fourier transform (DFT) based wide spectrum compressed sensing with complexexponential functions set as orthogonal sparse base, which is not the only orthogonal basis used to reconstruct wideband communication signals.Furthermore, a wideband spectrum estimating scheme based on cyclic spectrum, considering the 2D sparse signal of interest in the cyclic spectrum domain, was proposed in [15].With the purpose of optimizing the compressive sensing recovery of image and audio signal, [16] put forward an adaptive tree structured dictionary of orthogonal bases.For cognitive power line communication (PLC) systems, [17] applied several time-frequency transforms included DFT, discrete cosine transform (DCT), and discrete sine transform (DST) to spectrum sensing.Additionally, the cooperative sensing technique has attracted considerable research attention in wideband cognitive radio networks as it is capable of significantly improving the performance of wireless communication systems [18][19][20][21][22].To the best of our knowledge, these works are mostly concerning CR in terrestrial wireless network.However, there have been few corresponding works about LEO-based mobile satellite communication network in the literature up to now.The time-variant environment and limited processing capability in MSS make it markedly different from terrestrial wireless network.The scheme of wideband spectrum compressed sensing (WSCS) not only implies lower computation complexity and shorter processing time but also means lower sampling rate and more accurate reconstruction.Although an enormous amount of research has focused on wideband spectrum sensing in CRN, there still exit the following challenges to provide reliable spectrum detecting in MSS.Firstly, the whole spectrum of L/S band for MSS is extremely wide, which may cause sensing problems of immense computational costs and long processing delay; secondly, the acquirable spectrum in MSS is time-varying because of high speed motion of the satellite; finally, the channel in MSS is typical fading and shadowing channel due to variable space circumstances.Herein, we propose a brand new approach to resolve these problems: the novel wideband spectrum compressed sensing method based on DST, briefly referred to as DST-WSCS, with a set of sinusoidal functions as an orthogonal sparse basis, to provide more accurate recovery and lower processing complexity.Essentially, the superiority of the novel algorithm is attributed to the following two aspects.On the one hand, the representation of signal in DST domain is much sparser than that of DFT domain.On the other hand, the DST uses only real arithmetic, as opposed to the complex-valued DFT.Finally, we also present collaborative compressed sensing architecture to deal with channel fading in the scenario where satellite may be incapable of detecting active signals because there are no means to distinguish white space and fading effect apart.
The rest of this paper is organized as follows.Section 2 describes the system model.Then, Section 3 proposes novel wide spectrum compressed sensing scheme based on DST.Immediately following, a collaborative signal sensing scheme is developed in Section 4. Extensive simulation results are given in Section 5. Finally, conclusions are drawn in Section 6. composed of space segment LEO satellite and terrestrial segment PU, in which LEO satellite is taken as secondary system, while the primary system consisted of PU licensed in the incumbent mobile satellite services.There are broadly three steps to seek available frequency, first in which individual LEO satellite implements compressed sensing to make local decision on the state of PU via cognitive satellite link based on energy detection, second in which LEO choose vacant frequency to communicate with terminal according to sensing decision, and third in which it is selective to actualize cooperative sensing with the help of ISL in the overlapped coverage.

System Model Descriptions
Suppose that the received signal () at SU is spanning over a wide spectrum within (  ,   ) depicted in Figure 2, as well as the periodic spectrum sensing time framework, where 0 ∼  0 is used for implementing compressed spectrum sensing and  0 ∼  is reserved for transmitting data.
In particular, the basic assumptions are given as follows: (i) A slotted frequency segmentation structure is utilized, in which the entire wideband spectrum is equivalently divided into  orthogonal subbands (also known as slots) centered at {  }  =1 , and the subband index is denoted by  ∈ [1,], each of which has bandwidth   = (/) (∀ ∈ [1, 𝐿]).The spectrum set { PU1 ,  PU2 , . . .,  PU , . . .,  PU } occupied by potentially active PUs lied within (  ,   ) consecutively, whose frequency locations and power spectral density (PSD) levels are to be detected and identified.Emphatically, only the unoccupied subbands which accord with the frequency allocated to MSS listed in Table 1 are available.
(ii) It is assumed that signal PSD over each spectrum band  PU is smooth and exhibited independently from its two neighboring bands  PU(−1) and  PU(+1) .However, the PSD levels of PUs are dynamically varying, depending on whether they are active or not in a particular time and geographical area.Those temporarily idle subbands are termed spectral holes and available for opportunistic spectrum access by SUs.
(iii) Neither the number of occupied subbands  ( ≤ ) nor their locations are known to SUs.For simplicity, let Ω  denote the set of subband indices overlapped by { PU } ≤ =1 in the entire bandwidth, and Ω  = { |  represents the subband indices of  PU projection in   }.Obviously, the active PUs may present at any subband in Ω  .
(iv) In the spectrum sensing period, all SUs keep quiet as enforced by protocols, for example, at the media access control layer, so that only PUs are emitting power.Meanwhile, the locations and the number of active subbands remain constant but may vary among different sensing period.
Thus, the continuous signal received at the RF frontend of SUs () is composed of only active PUs signals and background noise.Suppose that there are  active PUs during the sensing interval, whose transmitted signals are expressed by   (),  = 1, . . ., .We define the channel impulse response between the th PU transmitter and the th SU receiver as ℎ  ().Thus, the signal   () is reached at the th satellite detector in the form of the convolution operator ℎ  () *   ().We assume that the channel gain can be varying slowly and flat similar to the time-invariant slotted narrowband subchannel during each detection time window.The received signal at the th detector is given as follows: where the environmental noise around the detector   () is white Gaussian with zero mean and variance .For the time window as  ∈ [0, 0 ],   () have discrete form as an length signal    by using sampling rate   over the observation time.According to the Nyquist sampling theory, the sampling rate is required at or above 2(  −   ) samples per second.However, it leads to excessive memory requirement and prohibitive energy cost for wideband spectrum sensing.This dilemma motivates us to employ CS to reduce sampling rate.To represent the discrete signal response on the  channels, we collect the frequency-domain samples into  × 1 vector, and writing every term in discrete version, we have where ,  ()  , and  ()  are the complex-valued frequency-domain discrete versions of   (), ℎ  (),   (), and   (), respectively.Given (2), the spectrum sensing aims at estimating the transmitted spectrum ∑  =1  ()  from PUs, which depends on the a priori knowledge of channel state information (CSI).Identifying spectrum holes is a more important task for SUs, rather than the signal spectrum itself.It turns out to be a spectrum detection problem, aiming at determining the binary spectrum occupancy state   (∀ ∈ [1, ]) of active PUs, one-one correspondence with { PU } ≤ =1 .The spectrum detection problem can be solved whether CSI has been acquired or not, because ∑  =1  ()  shares the same nonzero support as ()  , provided that the channel does not experience deep fading and ∑  =1 ℎ () does not have zero elements.In this case, (2) can be concisely written as As a result of the aforementioned property, we can use the same sampling and reconstruction method to recover either the transmitted ∑  =1  ()  for spectrum estimation or the received  ()  for detection of the spectrum support.Besides, we develop collaborative compressive sensing scheme to overcome the drawbacks arising from severe fading and shadowing corrupted in signal with noise in Section 4 later.

Compressive Sampling.
Due to low percentage of spectral occupancy by active PUs, it is believed that the received signal    at th SU is sparse in the Fourier domain.As the motivation of this paper, the received signal can be sparsely represented in the frequency domain as follows: where the  × 1 vector is the sparse representation of    in frequency domain and has just  ≪  nonzero elements.
is an  ×  DFT matrix.
The CS theory states that    can be accurately reconstructed by using a sub-Nyquist sampling rate; that is,   < 2(  −   ).Mathematically, the compressed samples can be written as where where w() is the AWGN noise vector.As long as Φ −1 obeys the rule of restricted isometry property (RIP) or the sensing matrix Φ −1 is incoherent with  −1 , the occupancy of wideband spectrum could be converted to the reconstruction of In order to deal with the noisy signals, as (1) shows, LASSO algorithm has utilized some variants by minimizing the usual sum of squared errors with a bound on the sum of the absolute values of the coefficients.Thus, LASSO has additional advantage of denoising performance.To get the spectrum reconstruction more accurately, we formulate the model based on LASSO as follows: where  bounds the amount of noise in signals.Great many convex optimization software packages have been developed to tackle the LASSO problem, for example, [23].

Decision-Making Criterion. After 𝑟 (𝑗)
reconstruction, it is indispensable to estimate the average amplitude of r()  within each identified subband for spectrum hole detection.Spectrum detection in local satellite detector can be obtained by a binary hypothesis test problem.We make decision in the th subband by The summary statistic of estimated signal energy over the interval of each subband is computed as the following form: and r()  is the estimated power spectrum, the decision criterion of   = 1 or   = 0 is given by where  is an invariable threshold factor chosen bigger than one.For the whole band, we could get a decision vector d() = ( ) over  subbands made by the th satellite detector that is played by SU depicted in Figure 1.Suppose that there are  detectors simultaneously sensing network wide.Then, we obtain the decision over the th spectrum bandwidth  PU of  ()  = 1 or  ()  = 0 following the principle We could make decision vector d() = ( 2 , . . .,  ()  , . . .,  ()  ) over  spectrum bandwidth in the th detector.Particularly, in the case of deep fading or shadowing channel, inaccurate detection could be induced when PUs are actually present, which is weighed with detection probability   .And in another case of strong noise, false alarm probability   may be produced.Thus,   should be kept as high as possible, whereas the smaller   will be the better.

A Novel Scheme of Compressed Sensing Based on DST
For dynamic spectrum access in mobile satellite network, we propose novel algorithm of DST-based compressed spectrum sensing to fulfill better performance in detection probability, reconstruction accuracy, and processing time.

Mathematical Problems in Engineering
The DST sequence of -length discrete signal () is defined by Suppose () is a 2-point subtractive symmetry continuation of () represented by Obviously, the 2-point DFT of () is given by Substituting ( 14) into (15) yields where, by definition,    = ( −/ )  .If we change the second index of summation using  = 2 − 1 −  and recall that    = 1 for integer , then, () could be expressed as ) ] ,  = 0, 1, . . ., 2 − 1. ( Finally, we obtain the relationship between () and (): Im[ ] implies the imaginary part of the term enclosed.DST of () can be computed by taking the 2-point DFT of () and multiplying the result by (/2) /2  , as shown in (18).We note that () is real, whereas () is complex.Theoretical derivation indicates that the DST basis has better spectral compaction and energy concentration than DFT, which means signal representation is much sparser in DST domain than that of DFT as shown in Figure 3.
On the basis of preceding analysis, the mathematical description of the scheme of DST-WSCS is similar to DFT-WSCS scheme with Fourier matrix being replaced by the DST matrix, which can be solved via following three steps: (i) Estimating sparse sequence r()  through compressed measurements    at individual th satellite detector: from (4)-( 6), the compressed samples can be represented as where (iii) Getting frequency response r()  from r  via Fast Fourier Transform (FFT) to identify the state of spectrum usage: We get the frequency estimation given by r() Finally, we can get the state of subbands occupied in the whole sensing spectrum according to decision-making rules (10)- (12).

Collaborative Compressed Sensing Scheme
The wireless channel in MSS is typical fading channel.To obtain reliable detection of spectrum utilization, it is urgent to overcome detrimental effects of the fading.Taking Iridium and 48 polar-orbit constellation for example, they are abundant in satellite resources and intersatellite link (ISL).The probabilities of multiple-satellite coverage is very high; that is, global single-satellite coverage (SSC) rate is up to 100%, double-satellite coverage (DSC) rate is more than 58%, and triple-satellite coverage (TSC) rate is more than 9%, which improve proportionally as the latitude increases, illustrated in Figure 4.All those provide the key factor for the application of cooperative sensing via the user link (UL) and ISL.
The scenario depicted in Figure 5 demonstrates that the spectrum occupancy in the area of DSC and TSC can be cooperatively detected by several satellites located in the same  or adjacent orbit planes, where   , means the coverage area of satellite detector  , and the black solid line and red solid line stand for intrasatellite link and intersatellite link, respectively.Plenty of work has shown that collaborative compressive spectrum sensing which exploits the built-in spatial diversity stands out as an effective approach to alleviate the problem of corrupted detection and the sensing burden, as well as cost of individual SU in fading channels.
In distributed mobile satellite network, detector  , implements compressive sampling individually to collect reducedrate measurements and possibly make local estimation and decision on the spectrum occupancy of the network.When it comes to data fusion, there are two primary implementation options: centralized and decentralized fusion.However, centralized cooperation scheme requires a fusion center that collects compressed samples or their local decisions from all satellite detectors and makes global sensing decision, which incurs high communication cost and renders the entire network vulnerable to satellite node failure.In this paper, we develop a decentralized cooperative sensing approach to effectively address the abovementioned issue, which relies on the local measurement {   }, the local sensing estimation { ()  }, and detection decision { ()   }.The abovementioned collaborative spectrum sensing fusion approach can be, respectively, viewed as measurement fusion, estimation fusion, and decision fusion.A short assessment of their differences would be that measurement fusion needs higher communication overhead compared with the other two, while estimation fusion cannot make  , recover the transmitted spectrum ∑  =1  ()  individually for the consent on their estimates in the absence of CSI.After weighing the pros and cons, we turn to decision fusion to directly recover the spectrum occupancy vector  ()   of PUs, which is a binary version with the same nonzero support of ∑  =1  ()  .To further alleviate communication overhead in terms of processing time and power consumption, we only use local communication among onehop neighbors for cooperation.Let (, ) be undirected connected graph depicting the connectivity of the multihop mobile satellite networks, where satellite nodes form the set of vertices  = {V 1 , . . ., V  , . . ., V || }, where || represents the number of satellites in space segment, and each edge (V  , V  ) ∈ , ∀V  , V  ∈  and V  ̸ = V  , connects an unordered pair of distinct nodes within one-hop neighborhood, which consisted of intra-ISL and inter-ISL.Let  (,) denotes the set of one-hop neighbors of satellite detector  , .We define the channel of UL between active PUs in the ground segment and satellite detector  , as detecting channel, as well as reporting channel from the individual detector  , to the satellite fusion node that is chosen from  (,) sharing the overlapping area with  , .
The procedure of collaborative compressed sensing follows the steps outlined as follows: (i) In the processing of detecting channel, each satellite detector  , performs local reconstruction to attain r()  from compressed measurement    via (8) independently and then makes a binary decision vector  ()   from the nonzero support of r()  by ( 12) to decide whether the spectrum bandwidth { PU } ≤ =1 is employed or not.
(ii) In terms of reporting channel, each satellite detector  , only communicates with its one-hop neighbors  (,) to reduce the transmission power consumed during sensing, and the exchanged decision { ()  } information percolates over the mobile satellite networks through multiple rounds to make global decision on the spectrum occupancy.
(iii) Assume there are  satellite detectors simultaneously sensing in the mobile satellite network for a particular area, such as the region of SSC, DSC, and TSC.For the utilization state of the whole band, the chosen Mathematical Problems in Engineering center satellite fuses the local decisions and makes final decision vector d by taking OR logic rule, as (22) shows, d = OR ( (1)   ,  (2)   , . . .,  ()  ) , or taking majority vote from the local binary decisions on  ()  , termed as MAJ rule, as (23) shows, where the symbol ⌊/2⌋ denotes the operation of taking the minimum integer that is greater than or equal to /2.Finally, we could make a decision vector d = ( d1 , d2 , . . ., d , . . ., d ) over  spectrum bandwidth for the whole band detection.

Simulations
This section provides extensive computer simulations to validate the performance of the proposed wideband spectrum compressed sensing scheme and evaluate the effects of key design parameters and network operating condition.Firstly, we describe the simulation setup and relevant performance metrics.Then, for the scenario without cooperation, we compare steady-state performance of proposed novel scheme DST-WSCS with conventional DFT-WSCS and investigate the running time, as well as the impact of system parameters such as SNR, the number of active PUs, and the compression ratio.Finally, for collaborative compressed sensing scenario, we also research the effects of number of cooperators and decision fusion rules on the sensing performance under different values of SNR.  of the proposed DST-WSCS scheme with that of conventional DFT-WSCS scheme in terms of detection probability, false alarm probability, reconstruction accuracy, and running time.The Monte-Carlo simulation method of 1000 times is used to evaluate the performance of single detector.rate increasing.More measurements means higher sampling burden and energy consumption, and a higher SNR can afford less compression ratio during data acquisition.The proposed DST-WSCS provides higher successful detection than that of the conventional method DFT-WSCS, although it heightens the false alarm probability to a certain extent in low SNR.On the other hand, we could reduce the false alarm probability through heightening threshold factor, as depicted in Figure 10.
Figure 7 shows the relationship between the detection performance and the number of active PUs based on DST method when / = 30%.It is equivalent to comparing performances for different degree of sparsity of the wideband spectrum.It can be concluded that detection performance improves as degree of sparsity of the wideband spectrum increases.This is because the theory of compressed sensing is based on the sparsity of spectrum.Additionally, the number of active primary signals has little impact on the false alarm probability.
(ii) Sparse Spectrum Recovery Performance of NMSE versus SNR and Compression Ratio. Figure 8 depicts the NMSE of the estimated spectrum for various values of compression ratio and SNR, when there are four active PUs.From Figure 8, we can see DST-WSCS outperforms DFT-WSCS.And higher SNR or compression rate means higher reconstruction accuracy.
(iii) Running Time for Individual Satellite Detector versus SNR and Compression Ratio. Figure 9 gives the execution time of the proposed detection algorithm for different compression ratio /, with four active PUs.It shows that the detection time is closely related to the compressive ratio, while it has little connection with SNR.The average detection time of DST-based compressive sensing is smaller than that of the DFT-based by at least 20%, due to the real value of DST operation.Summarily, the DST-based method costs less processing time, as well as lower processing complexity than DFT-based method, which is of great significance in MSS for the limitation of processing capability and time-varying spectrum.
(iv) The Effect of Threshold Factor on Detecting and ROC Performance.Figure 10 focuses on the effect of threshold factor on detection performance and the receiver operating characteristic (ROC) performance, when / is 50% and SNR is 0 dB.It indicates that either detection or false alarm detection degrades as threshold factor increases.Obviously, the proposed scheme of DST-WSCS performs better than DFT-WSCS in detection probability for the same false alarm rate.It should be pointed out that we take 1.4 as threshold factor through the full paper for the sake of the balance between detection and false alarm probability, which also explains that the false alarm probability of DST-WSCS is a little higher than that of DFT-WSCS, as illustrated in Figure 6.

Cooperative Compressed Sensing Performance
Based on DFT and DST (i) The Effect of Different Number of Cooperators on   and   .Figure 11 compares the cooperative compressed sensing performance based on DFT and DST for different cooperators with OR fusion rule, when there are four subbands randomly occupied and compression rate / is 50%.It is clear that the probability of detection is enhanced as the number of cooperators increases.This is because collaborative scheme brings spatial diversity gain.However, the false alarm probability becomes worse, especially for low SNR; this is according to the theoretical analysis mentioned above.Fortunately, the absolute value of   is extremely small and would not bring serious interference to PUs.In order to reduce   with accepted   that provides sufficient protection to PUs, we can alternatively heighten threshold factor as shown in Figure 10 or employ MAJ fusion rule as demonstrated in Figure 12.
(ii) The Effect of Decision Fusion Rule on   and   .Figure 12 focuses on the effect of decision fusion rule on cooperative compressed sensing performance based on DFT and DST with three cooperators, when there are four subbands randomly occupied and compression rate / is 50%.It is obvious that the utilization of MAJ fusion rule can reduce both the probability of detection and the false alarm probability greatly.Nevertheless, the DST method performs better for the robustness to noise in the respect of detection probability.We suggest that the OR fusion rule and MAJ fusion rule be employed jointly in practice, which can provide higher probability of detection and lower false alarm probability, respectively.

Conclusion
In this paper, we proposed a novel wideband spectrum compressed sensing scheme based on DST to resolve the frequency resource scarcity problem in mobile satellite systems.Furthermore, we also presented collaborative compressed sensing scheme to mitigate the adverse effect of fading on sensing performance.Additionally, given consideration to balance between detection probability and false alarm probability, we suggested the joint fusion rule based on the OR and MAJ fusion rule.Simulation results show that notable improvement in the detection probability, processing time, and reconstruction accuracy can be obtained using DSTbased sensing algorithm compared with conventional DFTbased compressive sensing, and spatial diversity gain can be achieved via cooperative sensing at the expense of little loss in false alarm probability, which can be compromised by either heightening threshold factor or the jointed decision fusion rule.More applicable decision fusion rule for collaboration and universal reconstruction algorithm will be our future direction of research.

Figure 3 :
Figure 3: Comparison of sparsity between DST and DFT.

Figure 5 :
Figure 5: The scenario of collaborative spectrum detection in mobile satellite network.

Figure 6 :
Figure 6: Detection performance versus SNR for different compression rate.
(i) P  and P  versus SNR, Compression Ratio, and the Number of Active PUs.

Figure 9 :
Figure 9: Running time of individual satellite detector for different compressive sensing schemes.

Figure 10 :
Figure 10: The performances of detecting and ROC versus threshold factor.

Figure 11 :
Figure 11: Cooperative compressed sensing performance versus the number of cooperators.

Figure 12 :
Figure 12: Detection performance comparison at different decision fusion rules.

Table 1 :
L/S band assignment for MSS.
(5)× compressive measurement matrix, columns of Φ can be viewed as a set of basis functions or matched filters, and different Φ corresponds to different kinds of sampling.When Φ is the  ×  identity matrix, model(5)means Nyquist uniform sampling.The  × 1 measurement vector is the projection of    onto Φ.And the compression rate is defined /.Noting that