CN105578480A - Undersampling frequency spectrum perception pre-decision method orienting broadband modulation converter - Google Patents

Undersampling frequency spectrum perception pre-decision method orienting broadband modulation converter Download PDF

Info

Publication number
CN105578480A
CN105578480A CN201510939335.4A CN201510939335A CN105578480A CN 105578480 A CN105578480 A CN 105578480A CN 201510939335 A CN201510939335 A CN 201510939335A CN 105578480 A CN105578480 A CN 105578480A
Authority
CN
China
Prior art keywords
branch road
variance
test statistics
frequency spectrum
average
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510939335.4A
Other languages
Chinese (zh)
Other versions
CN105578480B (en
Inventor
齐佩汉
耿雨晴
李赞
高锐
司江勃
关磊
熊天意
王盛云
王思勉
申鹏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201510939335.4A priority Critical patent/CN105578480B/en
Publication of CN105578480A publication Critical patent/CN105578480A/en
Application granted granted Critical
Publication of CN105578480B publication Critical patent/CN105578480B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W16/00Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
    • H04W16/14Spectrum sharing arrangements between different networks
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/382Monitoring; Testing of propagation channels for resource allocation, admission control or handover
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/0006Assessment of spectral gaps suitable for allocating digitally modulated signals, e.g. for carrier allocation in cognitive radio
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses an undersampling frequency spectrum perception pre-decision method orienting a broadband modulation converter, and mainly solves the problems in the prior art that false alarm probability is high and computation cost is high. The technical scheme is that 1. a receiving end performs compressive sampling on received signals so that m paths of samples are obtained and time domain energy of each branch of acquired samples is computed; 2. test statistics r<j,k> is computed according to time domain energy; 3. the decision threshold gamma<j,k> of test statistics is computed; and 4. test statistics is compared with the decision threshold, and existence of signals is decided: the signals exist through decision if r<j,k> is greater than gamma<j,k>, and the signals do not exist through decision if r<j,k> is less than or equal to gamma<j,k>. The undersampling frequency spectrum perception pre-decision method has advantages that influence of noise power on perception performance is low, computation complexity is low and influence of broadband Gaussian white noise non-sparsity on undersampling frequency spectrum perception can be effectively reduced, and the undersampling frequency spectrum perception pre-decision method can be used for broadband frequency spectrum perception of analog signal compressive sampling.

Description

Towards the pre-decision method of lack sampling frequency spectrum perception of wide-band modulation converter
Technical field
The invention belongs to communication technical field, relate to frequency spectrum perception technology, further relate to a kind of pre-decision method of lack sampling frequency spectrum perception towards wide-band modulation converter, can be used for the frequency spectrum perception in cognitive radio.
Background technology
Frequency spectrum perception is as the founder key technology of cognitive radio system, can be the feedback information that Dynamic Closed Loop system provides necessary, for cognitive radio system/network provides environmental information and decision-making foundation, for in order dynamically spectrum allocation may and reasonably spectrum regulatory provide safeguard, be the raising availability of frequency spectrum, improve wireless channel transmission conditions, intelligent network, the basis of initiatively hiding various interference, effectively carrying out dynamic spectrum resource management.Current frequency spectrum sensing method is mostly absorbed in narrower frequency range, is found frequency spectrum access chance, from Shannon theorem, available bandwidth directly determines theoretical maximum bit rate, the access frequency spectrum that narrow band spectrum perception provides obviously can not carry the traffic rate of user's request, and cognitive user needs to adopt broader frequency spectrum perception to find out in wider frequency range simultaneously and more accesses chance.Classical Nyquist sampling thheorem is pointed out, in order to undistorted reconstruct analog signal, the sample frequency of analog signal at least should be more than or equal to the twice of signal spectrum bandwidth, when cognitive user carries out broader frequency spectrum perception, need with higher resolution and extremely low power consumption collecting sample from very wide frequency range, this digital signal processor system of giving based on Shannon-nyquist sampling theorem brings and gathers front end conversion speed, height handles up large capacity cache space and the 'bottleneck' restrictions of back-end logic device operation highest frequency.
Compression sampling CS theory only needs the non-self-adapting linear measurement sample point of only a few, just can recover analog signal sparse in time domain or other transform domains by the method for convex optimization with great probability, directly can realize the conversion of signal from analog to information according to the method, this is reduce Analog-digital Converter speed and alleviate Digital Signal Processing pressure to provide new theoretical foundation.Given this, compression sampling is combined with broader frequency spectrum perception, make full use of the sparse characteristic treating perceived spectral, to complete the information gathering of broader frequency spectrum far below the analog-to-digital conversion rate of Nyquist sample rate, and completing broader frequency spectrum perception in real time with extremely low Digital Signal Processing expense, is solve one of broadband analog signal collection and high-speed digital signal transmission, storage and the effective way processing bottleneck.
At present, the broader frequency spectrum cognitive method that can be used for analog signal compression sampling is broadly divided into two classes:
The first kind is based on analog information transducer AIC, AIC lack sampling framework with finite discrete multi-tone signal add and for analog input signal model, by analog input signal is multiplied with random mark sequence, then the sampling of integration cumulative sum low pass is carried out, after obtaining compression collecting sample, reconstructing method can be utilized to draw primary signal or its statistical property, degree judgement is completed again by time domain or frequency domain energy detection method, this collection framework realizes simple, but when recovering continuous spectra signal in broadband, the predicament that input signal model mismatch produces the higher-dimension complex matrix computing that great reconstructed error or input signal model accurately bring can be faced.
Equations of The Second Kind is based on wide-band modulation converter MWC, MWC lack sampling framework is using the limited union with translation invariant subspace as analog input signal model, by input signal is multiplied with cycle random mark sequence on multiple branch road, realize the different weight factor down-conversions of input signal Subspace Decomposition, thus reduce required sampling rate, by converting unlimited measurement vector corresponding for continuous analog signals collecting to multiple measurement vector system to lack sampling sample framework establishment, then orthogonal matching pursuit algorithm is utilized to obtain support corresponding to occupied frequency range, it is reasonable that MWC lack sampling framework has signal model, support determines the advantages such as real-time and the realization of useful commercial device.But owing to there is white noise in the frequency range of broadband, broadband white noise is not sparse on time domain, frequency domain or other transform domains, adopts the lack sampling sample of wide-band modulation converter, directly carry out broader frequency spectrum perception, easily cause following problem:
(1) although there is not any primary user's signal in the broadband frequency range treating perception, only there is white Gaussian noise, the restructing algorithm of compressed sensing still can take situation according to openness providing in frequency band, and this will cause serious false alarm probability;
(2) process utilizing lack sampling sample to carry out signal reconstruction and support to determine, bring great calculation cost, and these computings is insignificant;
(3) strategy of the frequency spectrum perception Influence on test result communication of mistake, brings secondary user's communication and frequently interrupts, reduce the utilance of frequency band.
Summary of the invention
The object of the invention is to the above-mentioned deficiency of the frequency spectrum perception technology overcome based on wide-band modulation converter, a kind of pre-decision method of lack sampling frequency spectrum perception towards wide-band modulation converter is proposed, to effectively reduce the non-openness impact on lack sampling frequency spectrum perception of broadband Gaussian white noise, promote frequency spectrum perception accuracy, reduce calculation cost, improve band efficiency.
In order to complete above-mentioned purpose, the pre-decision method of lack sampling frequency spectrum perception towards wide-band modulation converter that the present invention proposes, comprises the steps:
(1) receiving terminal utilizes wide-band modulation converter to carry out compression sampling to the received signal, obtains m road sample y i(n), i=1,2 ..., m, n=0,1 ..., N-1, N are sample points, and calculate the time domain energy of each branch road collecting sample respectively
(2) test statistics is calculated according to time domain energy
Wherein, with be respectively the time domain energy of jth branch road and kth branch road, j=2,3 ..., m, k=1,2 ..., j-1, obtains m (m-1)/2 test statistics;
(3) test statistics r is calculated j,kdecision threshold γ j,k:
(3a) by the time domain energy of i-th branch road be transformed into frequency domain energy by test statistics r j,kfrequency domain is transformed into by time domain,
(3b) test statistics r is built j,kcumulative distribution function P f:
P f = 1 - &Phi; ( N &gamma; j , k - N 1 - 2 &gamma; j , k &rho; j , k + &gamma; j , k 2 )
Wherein ρ j,kbe with coefficient correlation, value is Φ () function expression is &Phi; ( x ) = &Integral; - &infin; x e - t 2 d t ;
(3c) according to the cumulative distribution function P of test statistics f, utilize CFAR criterion, calculate decision threshold γ j,k:
&gamma; j , k = &lsqb; &Phi; - 1 ( 1 - P f ) &rsqb; 2 &rho; j , k - N + ( &rho; j , k 2 - 1 ) ( &lsqb; &Phi; - 1 ( 1 - P f ) &rsqb; 4 - 2 N &lsqb; &Phi; - 1 ( 1 - P f ) &rsqb; 2 ) &lsqb; &Phi; - 1 ( 1 - P f ) &rsqb; 2 - N ,
Wherein, P ffor the default false alarm probability of each branch decision, Φ -1() is the inverse function of Φ () function;
(4) by test statistics r that step (2) obtains j,kthe threshold value γ obtained with step (3) j,krelatively, signal is determined whether, if there is statistic r j,kbe greater than thresholding γ j,k, then signal has been judged to; If there is not any statistic r j,kbe greater than thresholding γ j,k, then no signal is judged to.
The present invention has the following advantages:
1, the present invention is owing to utilizing a small amount of compression sampling sample, can effectively reduce the non-openness impact on lack sampling frequency spectrum perception of broadband Gaussian white noise;
2, the present invention has nothing to do due to its threshold value and noise variance, so perceptual performance is less by the impact of noise power, and effectively can resist the impact on noise uncertainty;
3, the present invention carries out simple computing module-square owing to only needing a small amount of sample, therefore computation complexity is low, can meet the real-time of frequency spectrum perception.
Accompanying drawing explanation
Fig. 1 is realization flow figure of the present invention;
Fig. 2 is the present invention's correct detection probability analogous diagram under different signal to noise ratio;
Fig. 3 is the present invention's false alarm probability analogous diagram under different signal to noise ratio;
Fig. 4 is in the present invention when there is noise uncertainty, detection probability simulation comparison figure under different signal to noise ratio;
Fig. 5 is the present invention's receiver performance curve analogous diagram under different false alarm probability.
Embodiment
The present invention is used for the broader frequency spectrum perception of analog signal compression sampling, and perception end is at each received over subchannels signal, and sampling to received signal processes.
With reference to Fig. 1, performing step of the present invention is as follows:
Step 1, perception end calculates the time domain energy T of each branch road collecting sample i td.
Perception end utilizes wide-band modulation converter to carry out compression sampling to the received signal, obtains m road sample y i(n), i=1,2 ..., m, n=0,1 ..., N-1;
Utilize branch road sample y in () calculates time domain energy
T i t d = &Sigma; n = 0 N - 1 | y i ( n ) | 2 ,
Wherein, m is the collection circuitry number of wide-band modulation converter, and N is sample points.
Step 2, calculates test statistics r according to time domain energy j,k
r j , k = T j t d T k t d ,
Wherein, with be respectively the time domain energy of jth and k branch road, j=2,3 ..., m, k=1,2 ..., j-1,
M (m-1)/2 test statistics is calculated according to the different values of j and k.
Step 3, according to the false alarm probability preset, when primary user's signal does not exist, calculates decision threshold γ j,k.
When primary user's signal does not exist, Received signal strength only comprises noise, i.e. input signal x (t)=ω (t), suppose ω (t) for average be 0, variance is white Gaussian noise, its decision threshold γ j,kcalculation procedure as follows:
(3a) by the time domain energy of i-th branch road be transformed into frequency domain energy by test statistics r j,kfrequency domain is transformed into by time domain,
(3b) test statistics r is built j,kcumulative distribution function P f:
(3b1) Y is calculated ithe statistical property of (k):
(3b11) structure is gathered, to branch road sample y according to wide-band modulation converter MWC compression in () carries out discrete Fourier transform, obtain branch road sample spectra Y i(k):
Y i ( k ) = &Sigma; l = - L 0 L 0 C i l X ( k - ln ) , k &Element; &lsqb; 0 , N ) ,
Wherein, c ilfor cycle pseudorandom ± 1 sequence p ithe cycle Fourier series expansion coefficient of (t), L is the subspace number required for complete expression input signal x (t) Fourier transform X (j Ω), and L can calculate according to following formula,
Wherein, f nYQthe equivalent Nyquist sampling rate of MWC, f pthe frequency of random mark sequence, f sit is low pass sampling rate;
(3b12) in the non-existent situation of primary user, input signal spectrum X (k) is average is 0, variance is discrete random variable, utilize average and the variance of X (k), calculate branch road sample spectra Y ithe statistical property of (k), it comprises:
Calculate Y iaverage E [the Y of (k) i(k)] and variance D [Y i(k)] be:
E[Y i(k)]=0, D &lsqb; Y i ( k ) &rsqb; = &sigma; &omega; 2 ,
Calculate Y ithe average of (k) real part and variance for:
Calculate Y ithe average of (k) imaginary part and variance for:
Calculate Y i(k) real part and Y ithe coefficient correlation of (k) imaginary part: (3b2) branch road sample spectra Y is calculated ithe average E of (k) mould square [|| Y i(k) || 2] and variance D [|| Y i(k) || 2]:
E &lsqb; | | Y i ( k ) | | 2 &rsqb; = &sigma; &omega; 2 , D &lsqb; | | Y i ( k ) | | 2 &rsqb; = &sigma; &omega; 4 ;
Wherein, it is the variance of input signal spectrum X (k);
(3b3) branch road sample spectra Y is calculated ithe coefficient correlation of (k) mould square:
(3b31) coefficient correlation on different frequent points is:
cov[||Y p(a)|| 2,||Y q(b)|| 2]=0,
(3b32) coefficient correlation on identical frequency is:
cov &lsqb; | | Y p ( k ) | | 2 , | | Y q ( k ) | | 2 &rsqb; = &sigma; w 4 ( &rho; 0 2 + &rho; 1 2 + &rho; 2 2 + &rho; 3 2 ) 2 ,
Wherein p, q=1,2 ..., m, p ≠ q, k=0,1 ..., N-1,
C in formula qmfor the cycle Fourier series expansion coefficient of cycle pseudorandom ± 1 sequence, represent and get real part, represent and get imaginary part;
(3b4) the average E [T of branch road frequency domain energy is calculated i fd], variance D [T i fd] and correlation coefficient ρ j,k:
E &lsqb; T i f d &rsqb; = E &lsqb; 1 N &Sigma; k = 0 N - 1 | Y i ( k ) | 2 &rsqb; = &sigma; &omega; 2
D &lsqb; T i f d &rsqb; = D &lsqb; 1 N &Sigma; k = 0 N - 1 | Y i ( k ) | 2 &rsqb; = &sigma; &omega; 4 N
&rho; j , k = C o v &lsqb; T j f d , T k f d &rsqb; D &lsqb; T j f d &rsqb; D &lsqb; T k f d &rsqb; = 1 N 2 &lsqb; &Sigma; k = 0 N - 1 ( C o v &lsqb; | Y j ( k ) | 2 , | Y k ( k ) | 2 &rsqb; ) &rsqb; D &lsqb; T j f d &rsqb; D &lsqb; T k f d &rsqb; = ( &rho; 0 2 + &rho; 1 2 + &rho; 2 2 + &rho; 3 2 ) 2
(3b5) according to the cumulative distribution function F of any two correlated Gaussian stochastic variable G and H ratio rr (), calculates test statistics r j,kcumulative distribution function P f:
F rr the expression formula of () is:
F R ( r ) = Pr ( G / H < r ) = &Phi; ( r&mu; H - &mu; G &sigma; G 2 - 2 r&rho;&sigma; G &sigma; H + r 2 &sigma; H 2 )
Wherein, μ gfor the average of Gaussian random variable G, for the variance of Gaussian random variable G, μ hfor the average of Gaussian random variable H, for the variance of Gaussian random variable H, ρ is the coefficient correlation of G and H.
As N>20, T i fdapproximate Gaussian distributed, its average and variance are respectively with with coefficient correlation be ρ j,k, test statistics r can be obtained thus j,kcumulative distribution function be:
P f = P ( T j f d T k f d &GreaterEqual; &gamma; ) = 1 - P ( T j f d T k f d < &gamma; ) = 1 - &Phi; ( &gamma; E &lsqb; T k f d &rsqb; - E &lsqb; T j f d &rsqb; ( D &lsqb; T j f d &rsqb; ) 2 - 2 &gamma;&rho; j , k D &lsqb; T j f d &rsqb; D &lsqb; T k f d &rsqb; + &gamma; 2 ( D &lsqb; T k f d &rsqb; ) ) = 1 - &Phi; ( N &gamma; - N 1 - 2 &gamma;&rho; j , k + &gamma; 2 )
Wherein ρ j,kit is jth branch road frequency domain energy with kth branch road frequency domain energy coefficient correlation, γ be inspection system
Metering r j,kdecision threshold, E [] represents average, and D [] represents variance;
(3c) according to the cumulative distribution function P of test statistics f, utilize CFAR criterion, calculate decision threshold γ j,k:
&gamma; j , k = &lsqb; &Phi; - 1 ( 1 - P f ) &rsqb; 2 &rho; j , k - N + ( &rho; j , k 2 - 1 ) ( &lsqb; &Phi; - 1 ( 1 - P f ) &rsqb; 4 - 2 N &lsqb; &Phi; - 1 ( 1 - P f ) &rsqb; 2 ) &lsqb; &Phi; - 1 ( 1 - P f ) &rsqb; 2 - N ,
Wherein, P ffor the default false alarm probability of each branch decision, Φ -1() is the inverse function of Φ () function.
Step 4: the test statistics r that step (2) is obtained j,kthe judging threshold γ obtained with step (3) j,kcompare, determine whether signal:
If there is statistic r j,kbe greater than thresholding γ j,k, then signal has been judged to;
If there is not any statistic r j,kbe greater than thresholding γ j,k, then no signal is judged to.
Effect of the present invention can be further illustrated by following emulation:
A, simulated conditions
The equivalent sampling speed of the wide-band modulation converter that analogue system adopts is f nYQ=6GHz, collection circuitry number is m=20, and the cycle of random mark ± 1 sequence is T p=7.5ns, frequency are f p=1/T p, in one-period, the random chip number of equivalence is L=45, and each passage low-pass filtering cut-off frequency used is f s/ 2, single channel sampling rate is f s=f p; At frequency range (0, f nYQ/ 2) in, co-exist in N=2 signal, the chip rate of each signal is s r=1.024MBaud, the carrier frequency of each signal produces at random, and the power of N number of signal is identical, and the ratio of the gross power and noise power that define N number of signal is signal to noise ratio, and under each signal to noise ratio, carry out 1000 emulation, the false alarm probability preset is P fa, utilize MWC converter to carry out compression collection, sample points K=100 used.
B, emulation content and result
Emulation 1: be-20dB ~-10dB and default false alarm probability P in signal to noise ratio faunder the condition of=0.01, emulate correct detection probability of the present invention, emulation as shown in Figure 2.
As seen from Figure 2, utilize K=100 lack sampling sample points, when signal to noise ratio is more than or equal to-8dB, correct detection probability of the present invention can reach and be greater than 98%, and the present invention can obtain very superior pre-judgement performance in wider SNR ranges.
Emulation 2: be-20dB ~-10dB and default false alarm probability P in signal to noise ratio faunder the condition of=0.01, emulate false alarm probability of the present invention, emulation as shown in Figure 3.
As seen from Figure 3, utilize K=100 lack sampling sample points, in the scope of signal to noise ratio-20dB ~-10dB, false alarm probability of the present invention suppresses substantially near default false alarm probability value, and the threshold value that the present invention calculates is reasonable, effectively can overcome the generation of false-alarm.
Emulation 3: signal to noise ratio be-20dB ~-10dB, false alarm probability be 0.01 and noise uncertainty ρ=1.4 condition under, correct detection probability of the present invention is emulated, and itself and the simulation curve that there is not noise uncertainty are contrasted, simulation result is as shown in Figure 4.
As seen from Figure 4, when there is noise uncertainty, correct detection probability of the present invention not too large change, performance curve of the present invention overlaps substantially when there is noise uncertainty and there is not noise uncertainty, illustrates that the present invention can effectively to antinoise uncertainty.
Emulation 4: be under the condition of-16dB in signal to noise ratio, the present invention is emulated towards the receiver performance curve of the pre-decision method of lack sampling frequency spectrum perception of wide-band modulation converter, and itself and the simulation curve that there are noise uncertainty ρ=1.4 are contrasted, simulation result is as shown in Figure 5.
As seen from Figure 5, the pre-decision algorithm that the present invention provides can under less false alarm probability and lower signal to noise ratio condition, obtain comparatively superior perceptual performance, when there is noise uncertainty, receiver performance curve of the present invention is also substantially constant, illustrates that the present invention has robustness to noise uncertainty.
Comprehensive above-mentioned analysis of simulation result, the present invention under less false alarm probability and lower signal to noise ratio condition, can obtain comparatively superior perceptual performance, when there is noise uncertainty, perceptual performance is substantially constant, illustrates that the present invention can effectively to antinoise uncertainty.

Claims (4)

1., towards the pre-decision method of lack sampling frequency spectrum perception of wide-band modulation converter, comprise the steps:
(1) receiving terminal utilizes wide-band modulation converter to carry out compression sampling to the received signal, obtains m road sample y i(n), i=1,2 ..., m, n=0,1 ..., N-1, N are sample points, and calculate the time domain energy of each branch road collecting sample respectively
(2) test statistics is calculated according to time domain energy
Wherein, with be respectively the time domain energy of jth branch road and kth branch road, j=2,3 ..., m, k=1,2 ..., j-1, obtains m (m-1)/2 test statistics;
(3) test statistics r is calculated j,kdecision threshold γ j,k:
(3a) by the time domain energy of i-th branch road be transformed into frequency domain energy T i fd, by test statistics r j,kfrequency domain is transformed into by time domain,
(3b) test statistics r is built j,kcumulative distribution function P f:
Wherein ρ j,kbe with coefficient correlation, value is Φ () function expression is
(3c) according to the cumulative distribution function P of test statistics f, utilize CFAR criterion, calculate decision threshold γ j,k:
Wherein, P ffor the default false alarm probability of each branch decision, Φ -1() is the inverse function of Φ () function;
(4) by test statistics r that step (2) obtains j,kthe threshold value γ obtained with step (3) j,kcompare, determine whether signal, if there is statistic r j,kbe greater than thresholding γ j,k, then signal has been judged to; If there is not any statistic r j,kbe greater than thresholding γ j,k, then no signal is judged to.
2. the pre-decision method of lack sampling frequency spectrum perception towards wide-band modulation converter according to claim 1, wherein calculates each branch road collecting sample y in step (1) ithe time domain energy of (n) calculate according to following formula:
3. the pre-decision method of lack sampling frequency spectrum perception towards wide-band modulation converter according to claim 1, wherein in step (3a) by the time domain energy of i-th branch road be transformed into frequency domain energy T i fd, be utilize handkerchief Savall law of conservation of energy, calculate according to following formula:
Wherein, Y ik () is y ithe discrete Fourier transform of (n), i=1,2 ..., m, n=0,1 ..., N-1.
4. the pre-decision method of lack sampling frequency spectrum perception towards wide-band modulation converter according to claim 1, wherein builds test statistics r in step (3b) j,kcumulative distribution function P f, carry out as follows:
(3b1) Y is calculated ithe statistical property of (k):
(3b11) structure is gathered, to branch road sample y according to wide-band modulation converter MWC compression in () carries out discrete Fourier transform, obtain branch road sample spectra Y i(k):
Wherein, c ilfor cycle pseudorandom ± 1 sequence p ithe cycle Fourier series expansion coefficient of (t), L is the subspace number required for complete expression input signal x (t) Fourier transform X (j Ω), and L can calculate according to following formula,
Wherein, f nYQthe equivalent Nyquist sampling rate of MWC, f pthe frequency of random mark sequence, f sit is low pass sampling rate;
(3b12) in the non-existent situation of primary user, input signal spectrum X (k) is average is 0, variance is discrete random variable, utilize average and the variance of X (k), calculate branch road sample spectra Y ithe statistical property of (k):
Calculate Y iaverage E [the Y of (k) i(k)] and variance D [Y i(k)] be:
E[Y i(k)]=0,
Calculate Y ithe average of (k) real part and variance for:
Calculate Y ithe average of (k) imaginary part and variance for:
Calculate Y i(k) real part and Y ithe coefficient correlation of (k) imaginary part:
(3b2) branch road sample spectra Y is calculated iaverage E [the ‖ Y of (k) mould square i(k) ‖ 2] and variance D [‖ Y i(k) ‖ 2]:
Wherein, it is the variance of input signal spectrum X (k);
(3b3) branch road sample spectra Y is calculated ithe coefficient correlation of (k) mould square:
(3b31) coefficient correlation on different frequent points is:
cov[‖Y p(a)‖ 2,‖Y q(b)‖ 2]=0,
(3b32) coefficient correlation on identical frequency is:
Wherein p, q=1,2 ..., m, p ≠ q, k=0,1 ..., N-1,
c pm, c qmfor the cycle Fourier series expansion coefficient of cycle pseudorandom ± 1 sequence, represent and get real part, represent and get imaginary part;
(3b4) the average E [T of branch road frequency domain energy is calculated i fd], variance D [T i fd] and correlation coefficient ρ j,k:
(3b5) according to the cumulative distribution function F of any two correlated Gaussian stochastic variable G and H ratio rr (), calculates test statistics r j,kcumulative distribution function P f:
F rr the expression formula of () is:
Wherein, μ gfor the average of Gaussian random variable G, for the variance of Gaussian random variable G, μ hfor the average of Gaussian random variable H, for the variance of Gaussian random variable H, ρ is the coefficient correlation of G and H.
As N>20, T i fdapproximate Gaussian distributed, its average and variance are respectively with with coefficient correlation be ρ j,k, test statistics r can be obtained thus j,kcumulative distribution function be:
Wherein ρ j,kit is jth branch road frequency domain energy with kth branch road frequency domain energy coefficient correlation, γ is test statistics r j,kdecision threshold, E [] represents average, and D [] represents variance.
CN201510939335.4A 2015-12-15 2015-12-15 The pre- decision method of lack sampling frequency spectrum perception towards wide-band modulation converter Active CN105578480B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510939335.4A CN105578480B (en) 2015-12-15 2015-12-15 The pre- decision method of lack sampling frequency spectrum perception towards wide-band modulation converter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510939335.4A CN105578480B (en) 2015-12-15 2015-12-15 The pre- decision method of lack sampling frequency spectrum perception towards wide-band modulation converter

Publications (2)

Publication Number Publication Date
CN105578480A true CN105578480A (en) 2016-05-11
CN105578480B CN105578480B (en) 2019-04-19

Family

ID=55888034

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510939335.4A Active CN105578480B (en) 2015-12-15 2015-12-15 The pre- decision method of lack sampling frequency spectrum perception towards wide-band modulation converter

Country Status (1)

Country Link
CN (1) CN105578480B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106255125A (en) * 2016-09-30 2016-12-21 广州粤讯信息科技有限公司 A kind of spectrum regulatory system
CN107181548A (en) * 2017-05-15 2017-09-19 西安电子科技大学 One kind compression frequency spectrum perception performance improvement method
CN107332566A (en) * 2017-06-19 2017-11-07 四川大学 Supported collection quick recovery method based on MWC
CN107426737A (en) * 2017-09-22 2017-12-01 哈尔滨工业大学 Broader frequency spectrum cognitive method based on single-pass configuration modulation wide-band transducer
CN109495198A (en) * 2019-01-25 2019-03-19 西安电子科技大学 Lack sampling frequency spectrum sensing method based on the detection of residual error correlation matrix
CN113765607A (en) * 2021-03-15 2021-12-07 全球能源互联网研究院有限公司 Frequency spectrum detection method and device
CN115276799A (en) * 2022-07-27 2022-11-01 西安理工大学 Decision threshold self-adapting method for undersampling modulation and demodulation in optical imaging communication

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102946288A (en) * 2012-11-23 2013-02-27 西安电子科技大学 Compressed spectrum sensing method based on autocorrelation matrix reconstitution
CN104270234A (en) * 2014-09-12 2015-01-07 中国科学院上海微系统与信息技术研究所 Broadband signal detection and identification method based on Nyquist under-sampling
WO2015077371A1 (en) * 2013-11-19 2015-05-28 Massachusetts Institute Of Technology METHODS AND APPARATUSES FOR MONITORING OCCUPANCY OF WIDEBAND GHz SPECTRUM, AND SENSING RESPECTIVE FREQUENCY COMPONENTS

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102946288A (en) * 2012-11-23 2013-02-27 西安电子科技大学 Compressed spectrum sensing method based on autocorrelation matrix reconstitution
WO2015077371A1 (en) * 2013-11-19 2015-05-28 Massachusetts Institute Of Technology METHODS AND APPARATUSES FOR MONITORING OCCUPANCY OF WIDEBAND GHz SPECTRUM, AND SENSING RESPECTIVE FREQUENCY COMPONENTS
CN104270234A (en) * 2014-09-12 2015-01-07 中国科学院上海微系统与信息技术研究所 Broadband signal detection and identification method based on Nyquist under-sampling

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
OR YAIR,ET AL.: "Uniform Linear Array Based Spectrum Sensing from Sub-Nyquist Samples", 《IEEE XPLORE DIGITAL LIBRARY》 *
陈坤,等: "MWC欠采样方法中伪随机序列的改进", 《信息技术》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106255125A (en) * 2016-09-30 2016-12-21 广州粤讯信息科技有限公司 A kind of spectrum regulatory system
CN107181548A (en) * 2017-05-15 2017-09-19 西安电子科技大学 One kind compression frequency spectrum perception performance improvement method
CN107332566A (en) * 2017-06-19 2017-11-07 四川大学 Supported collection quick recovery method based on MWC
CN107332566B (en) * 2017-06-19 2020-09-08 四川大学 MWC-based support set rapid recovery method
CN107426737A (en) * 2017-09-22 2017-12-01 哈尔滨工业大学 Broader frequency spectrum cognitive method based on single-pass configuration modulation wide-band transducer
CN107426737B (en) * 2017-09-22 2020-10-02 哈尔滨工业大学 Broadband spectrum sensing method based on single-channel structure modulation broadband converter
CN109495198A (en) * 2019-01-25 2019-03-19 西安电子科技大学 Lack sampling frequency spectrum sensing method based on the detection of residual error correlation matrix
CN113765607A (en) * 2021-03-15 2021-12-07 全球能源互联网研究院有限公司 Frequency spectrum detection method and device
CN113765607B (en) * 2021-03-15 2022-08-30 全球能源互联网研究院有限公司 Frequency spectrum detection method and device
CN115276799A (en) * 2022-07-27 2022-11-01 西安理工大学 Decision threshold self-adapting method for undersampling modulation and demodulation in optical imaging communication

Also Published As

Publication number Publication date
CN105578480B (en) 2019-04-19

Similar Documents

Publication Publication Date Title
CN105578480A (en) Undersampling frequency spectrum perception pre-decision method orienting broadband modulation converter
CN103873170B (en) A kind of frequency spectrum detecting method of compressed sensing under blind sparse condition
CN101630982A (en) Broadband cognitive radio frequency spectrum detection method based on parallelly compressed sensing
CN103346798B (en) Method for acquiring signal at sampling frequency lower than Nyquist frequency
CN109495198A (en) Lack sampling frequency spectrum sensing method based on the detection of residual error correlation matrix
CN103117818B (en) Broadband spectrum sensing method based on space-frequency joint compressed sensing
CN103532645B (en) The compression frequency spectrum sensing method that a kind of observing matrix is optimized
CN107634768A (en) A kind of MWC compression sampling wideband digital receiver PDW forming methods based on FPGA
CN104703216A (en) Energy-efficient-based detecting method of multi-task Bayes compressed sensing broadband spectrum
CN112287796B (en) Radiation source identification method based on VMD-Teager energy operator
CN104780008A (en) Broadband spectrum sensing method based on self-adaptive compressed sensing
CN101988935A (en) Digital down converter-Hilbert-Huang transform based instantaneous frequency measuring method
CN106301631B (en) A kind of coprime lack sampling frequency spectrum sensing method and its device based on Subspace Decomposition
CN106027179A (en) Wideband frequency spectrum sensing method based on comprehensive co-prime analysis and device thereof
CN106170139A (en) A kind of frequency spectrum detecting method and system
CN102421105A (en) Method for modeling ultra wide-band (UWB) compressive sensing wireless channel
CN103684634B (en) Based on the compression frequency spectrum sensing method of locating information in heterogeneous wireless sensor net
CN103138859B (en) Cognition wireless broadband frequency spectrum compressed sensing method based on backtracking and centralized type cooperation
CN105429918A (en) Rapid acquisition method used for deep space measurement and control signals
CN103929256B (en) A kind of multiframe compressed sensing signal spectrum detection method
CN109379745A (en) A kind of multi-user Cooperation frequency spectrum perception strategy study method towards 5G
CN106452626A (en) Broadband spectrum compression sensing based on multi-group relatively-prime sampling
CN103220016A (en) Generation system and method of pulse ultra wideband system orthogonal sparse dictionary
CN110944336B (en) Time spectrum sensing method based on limited new information rate
Gong et al. Property investigation on the additive white Gaussian noise after sub-Nyquist sampling

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant