CN101951271B - Compressive sampling based ultra wideband (IR-UWB) signal detection method - Google Patents

Abstract

The invention relates to a compressive sampling based ultra wideband (IR-UWB) signal detection method, which comprises an IR-UWB signal compressive sampling receiving system, and the IR-UWB signal compressive sampling receiving system comprises a multi-channel parallel sampling unit for carrying out sampling on a plurality of channels of the IR-UWB signals, a measuring waveform generator for respectively sending measuring waveforms to the channels of the multi-channel sampling unit, and a digital back-end processing component for receiving the measuring values of the channels sampled by the multi-channel sampling unit; and the channels carry out linear projection on the IR-UWB signals respectively according to the measuring waveforms generated by the measuring waveform generator. The method comprises the following steps: obtaining a compressive sampling sequence of pilot-symbol receiving signals, obtaining a compressive sampling sequence of data-symbol receiving signals, and obtaining a symbol decision. The compressive sampling based IR-UWB (ultra wideband) signal detection method can realize the low-sampling-rate and low-hardware-cost on the detection of the IR-UWB signals without high sampling rate, analog delay line and compressive sampling receiving mode with accurate channel estimation.

Description

Impulse ultra-wideband signal detection method based on compression sampling

Technical field

The present invention relates to a kind of impulse ultra-wideband signal detection method, relate in particular to a kind of impulse ultra-wideband signal detection method based on compression sampling.

Background technology

In the short-distance wireless communication that with IR-UWB (Impulse Radio-Ultra WideBand, pulse ultra-broad band are called for short " IR-UWB ") is physical-layer techniques was used, " pilot tone is assisted+the BPM modulation " this signaling method used the most generally.In traditional IR-UWB receiving side signal case, full resolution digital receiver, simulation correlation receiver, Rake receiver etc. all can be used as the symbol detection of this kind signaling method, are summarized as follows respectively:

The full resolution digital correlation receives

Carry out to received signal the sampling of Nyquist or higher rate, the correlate template of sample sequence acquisition that receives signal according to pilot pulse (can use MLSE or GML CLEAN scheduling algorithm to obtain clean template, also can use the time average of a plurality of reception of impulse sequences as noisy template), then the sample sequence with data pulse reception signal carries out related operation, and relevant output is as the symbol judgement statistic.

Simulation correlation reception

Directly use the reception waveform of pilot pulse as the simulation correlate template of subsequent data pulses reception signal, the sampling of relevant output is used for symbol judgement, required sampling rate is frame rate, Nyqusit speed greatly reduces relatively, but template signal is noisy, so the relative ideal template has loss on the performance; If use the time average waveform of a plurality of pilot pulses reception waveforms as template, approach when the time average of use sample sequence is as template in the performance that then obtains and the full resolution digital received.

Rake receives

Use the multipath composition in a plurality of associated branch coupling reception signals, the signal energy that each branch road is collected and as the symbol judgement amount, required sampling rate also is frame rate, but performance depends critically upon the levels of precision of channel estimating, when branch road quantity is abundant, approaching with the performance of simulation correlation reception.

All there is actual realization hard problem in above-mentioned several signal detecting method.The full resolution digital correlation receives needs the Nyquist sampling rate, and this can not " low-power consumption, low cost " realize with the existing hardware level basically for the wide IR-UWB signal of very bandwidth; Simulation correlation reception needs long broadband analog delay line, and hardware is realized also very difficult; Rake receives required accurate channel estimating take larger pilot-frequency expense as cost, and needs that Rake branch road quantity is abundant, each branch road timing controlled is enough accurate, and hardware complexity is very high.

Summary of the invention

The technical problem that the present invention solves is: a kind of impulse ultra-wideband signal detection method based on compression sampling is provided, overcomes in the prior art, hardware complexity, the technical problem that requires high sampling rate and accurately need accurate channel estimating.

Technical scheme of the present invention is: a kind of impulse ultra-wideband signal detection method based on compression sampling is provided, comprise impulse ultra-wideband signal compression sampling receiving system, described impulse ultra-wideband signal compression sampling receiving system comprises the multi-channel parallel sampling unit that divides a plurality of passages to sample to described impulse ultra-wideband signal, send respectively the measured waveform generator of measured waveform to each passage of described multi-channel sampling unit, receive the digital back-end reception ﹠ disposal assembly through the measured value of described multi-channel sampling unit sampling, each passage carries out linear projection according to the measured waveform that described measured waveform generator produces to described impulse ultra-wideband signal respectively, comprises the steps:

Obtain the compression sampling sequence that frequency pilot sign receives signal: measure by the compression that continuous multiframe is received signal, obtain the sample sequence that frequency pilot sign receives signal, as the template sequence of follow-up coherent detection;

Obtain the compression sampling sequence that data symbol receives signal: measure by the compression that continuous multiframe is received signal, obtain the compression sampling sequence that data symbol receives signal, as the correlated series for the treatment of of follow-up coherent detection;

Obtain symbol judgement: by coherent detection, obtain symbol judgement, i.e. testing result;

Further technical scheme of the present invention is: establishing used number of pilot symbols is N _p, the pulse that each symbol uses repeats transmission times and is N _f, M required measured value of symbol judgement shared on the measurement to continuous D frame signal, makes N _D=N _f/ D, then a total N _pN _DCriticize compression and measure sequence, it is y that n criticizes sequence _p[n]=Ф r _Prj+ Ф w _p[n], wherein: Ф is that matrix, r are measured in corresponding compression _PrjThe reception signal that is a pulse drops on the virtual sample sequence that compresses in the drop shadow spread, w[n] be the virtual sample sequence that receives the noise in the signal.

Further technical scheme of the present invention is: also comprise feedback loop, described feedback loop with the result partial feedback of described digital back-end processing components to described measured waveform generator, the measured waveform generator produces new measured waveform according to feedback information, the template sequence of follow-up coherent detection obtain and treat obtaining of correlated series, all use the new measured waveform that produces to obtain.。

Further technical scheme of the present invention is: for the situation that feedback loop is arranged, in the compression sampling sequence step of obtaining frequency pilot sign reception signal, frequency pilot sign receives signal and is divided into first and second portion, comprises the steps:

Obtain the compression sampling sequence that first's frequency pilot sign receives signal: transmitting terminal sends

Individual frequency pilot sign compresses measurement with the measurement matrix Ф of completely random to corresponding reception signal, and is total

Criticize the compression measured value, n criticizes sequence and is designated as Wherein

It is the virtual sample sequence that this part frequency pilot sign receives noise in the signal.

Obtain subspace estimation: the compression sampling sequence that adopts first's frequency pilot sign to receive signal is estimated signal subspace H.

Obtain the subspace compression sampling sequence that the second portion frequency pilot sign receives signal: the digital back-end processing components is after the estimation that obtains signal subspace H, it is fed back to the measured waveform generator, the measured waveform generator produces new measured waveform according to the signal subspace estimated information that obtains, namely

Wherein G is random matrix, to remaining in the frequency pilot sign

Individual reception signal compresses measurement, can obtain altogether

Criticize compression and measure sequence, be designated as

Wherein

The virtual sample sequence that this part frequency pilot sign receives noise in the signal,

Mean sequence Will be as the template sequence of coherent detection afterwards.

Further technical scheme of the present invention is: in obtaining the subspace estimation step, the compression sampling sequence that adopts first's frequency pilot sign to receive signal is carried out the signal subspace estimation, comprises the steps:

Compression sampling sequence initialization with the subspace: dictionary matrix V=Ф Ψ (v _iThe i row of representing matrix V), measure residual error

Estimated result L=[] and iterations t=1;

Obtain the sequence number with the matched atoms of residual error: from the dictionary matrix, seek the sequence number with the matched atoms of residual error,

Obtain the coefficient of current atom and upgrade residual error: upgrade residual error

Judge whether to stop: establishing the iterations thresholding is K ', if t＞K ', then iteration stops, and skips to for the 4th step; Otherwise upgrade t=t+1, with current search to the atom sequence number be increased in the estimated result i.e. L=[L, l _t], and the rebound first step;

Obtain the estimation of subspace: the estimated result that obtains signal subspace according to L is

Technique effect of the present invention is: a kind of impulse ultra-wideband signal detection method based on compression sampling is provided, by obtain compression sampling sequence that frequency pilot sign receives signal as correlate template, obtain data symbol and receive the compression sampling sequence of signal as treating correlated series, calculating impulse ultra-wideband signal is finished in relevant output as the symbol judgement parameter detection.The impulse ultra-wideband signal detection method that the present invention is based on compression sampling need not high sampling rate, need not analog delay line, need not accurate channel estimating, so that the reception of IR-UWB signal can be hanged down sampling rate, low hardware cost is realized.

Description of drawings

Fig. 1 is structural representation of the present invention.

Fig. 2 is flow chart of the present invention.

Fig. 3 is that the present invention has the flow chart that obtains sample sequence in the feedback loop situation.

Fig. 4 is the flow chart that the present invention obtains subspace estimation.

Embodiment

Below in conjunction with specific embodiment, technical solution of the present invention is further specified.

Such as Fig. 1, Fig. 2, Fig. 3, shown in Figure 4, the specific embodiment of the present invention is: the impulse ultra-wideband signal detection method that the present invention is based on compression sampling, comprise impulse ultra-wideband signal compression sampling receiving system, described impulse ultra-wideband signal compression sampling receiving system comprises the multi-channel parallel sampling unit 2 that divides a plurality of passages to sample to described impulse ultra-wideband signal, send respectively the measured waveform generator 1 of measured waveform to each passage of described multi-channel sampling unit, receive the digital back-end processing components 3 through the measured value of described multi-channel sampling unit sampling, each passage carries out linear projection according to the measured waveform that described measured waveform generator produces to described impulse ultra-wideband signal respectively.

Among the present invention, pilot tone and data symbol signal arrive being described below of receiving terminal from the transmitting terminal channel:

Transmitting is s (t),

$s\left(t\right)=d\left(t\right)*p\left(t\right)$

$=(\underset{i=0}{\overset{N_p-1}{\mathrm{\Σ}}}\underset{j=0}{\overset{N_f-1}{\mathrm{\Σ}}}\mathrm{\δ}(t-{\mathrm{iN}}_f{T}_f-{\mathrm{jT}}_f)+\underset{i=0}{\overset{N_s-1}{\mathrm{\Σ}}}\underset{j=0}{\overset{N_f-1}{\mathrm{\Σ}}}{b}_i\mathrm{\δ}(t-{\mathrm{iN}}_f{T}_f-{\mathrm{jT}}_f-N_p{N}_f{T}_f))*p\left(t\right)$

$=\underset{i=0}{\overset{N_p-1}{\mathrm{\Σ}}}\underset{j=0}{\overset{N_f-1}{\mathrm{\Σ}}}p(t-i{N}_f{T}_f-{\mathrm{jT}}_f)+\underset{i=0}{\overset{N_s-1}{\mathrm{\Σ}}}\underset{j=0}{\overset{N_f-1}{\mathrm{\Σ}}}{b}_{i}p(t-{\mathrm{iN}}_f{T}_f-{\mathrm{jT}}_f-N_p{N}_f{T}_f)---\left(1\right)$

N _pExpression frequency pilot sign number;

N _sExpression data symbol number

N _fRepresent that the pulse that each symbol uses repeats transmission times;

P (t) represents the transmitted waveform, usually uses each rank derived function waveform of Gaussian pulse, and pulse duration is designated as T _p

T _fIndicating impulse sends the interval.

Channel is h (t)

$h\left(t\right)=\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{\mathrm{\α}}_l\mathrm{\δ}(t-{\mathrm{\τ}}_l)---\left(2\right)$

L represents the multipath sum,

α _lExpression single footpath gain,

τ _lThe expression list is the time of advent directly.

The reception signal is r (t):

$r\left(t\right)=d\left(t\right)*p\left(t\right)*a_1\left(t\right)*h\left(t\right)*a_2\left(t\right)$

$=s\left(t\right)*h\left(t\right)*a_1\left(t\right)*a_2\left(t\right)$

$=(\underset{i=0}{\overset{N_p-1}{\mathrm{\Σ}}}\underset{j=0}{\overset{N_f-1}{\mathrm{\Σ}}}p(t-{\mathrm{iN}}_f{T}_f{\mathrm{jT}}_f)+\underset{i=0}{\overset{N_s-1}{\mathrm{\Σ}}}\underset{j=0}{\overset{N_f-1}{\mathrm{\Σ}}}{b}_{i}p(t-{\mathrm{iN}}_f{T}_f-{\mathrm{jT}}_f-N_p{N}_f{T}_f))*\left(\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{\mathrm{\α}}_l\mathrm{\δ}(t-{\mathrm{\τ}}_l)\right)*a_1\left(t\right)*a_2\left(t\right)---\left(3\right)$

$=(\underset{i=0}{\overset{N_p-1}{\mathrm{\Σ}}}\underset{j=0}{\overset{N_f-1}{\mathrm{\Σ}}}\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{\mathrm{\α}}_{l}p(t-{\mathrm{\τ}}_l-{\mathrm{iN}}_f{T}_f-{\mathrm{jT}}_f)+\underset{i=0}{\overset{N_s-1}{\mathrm{\Σ}}}\underset{j=0}{\overset{N_f-1}{\mathrm{\Σ}}}{b}_i\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{\mathrm{\α}}_{l}p(t-{\mathrm{\τ}}_l-{\mathrm{iN}}_f{T}_f-{\mathrm{iT}}_f-N_p{N}_f{T}_f))*a_1\left(t\right)*a_2\left(t\right)$

A wherein ₁(t), a ₂(t) be respectively transmitting antenna, reception antenna response, do not consider the wave distortion problem, formula (3) is reduced to:

$r\left(t\right)=\underset{i=0}{\overset{N_p-1}{\mathrm{\Σ}}}\underset{j=0}{\overset{N_f-1}{\mathrm{\Σ}}}\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{\mathrm{\α}}_{l}p(t-{\mathrm{\τ}}_l-{\mathrm{iN}}_f{T}_f-{\mathrm{jT}}_f)+\underset{i=0}{\overset{N_s-1}{\mathrm{\Σ}}}\underset{j=0}{\overset{N_f-1}{\mathrm{\Σ}}}{b}_i\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{\mathrm{\α}}_{l}p(t-{\mathrm{\τ}}_l-{\mathrm{iN}}_f{T}_f-{\mathrm{jT}}_f-N_p{N}_f{T}_s)---\left(4\right)$

Through band pass filter g (t) filtering out-of-band noise:

$x\left(t\right)=(r\left(t\right)+n\left(t\right))*g\left(t\right)=r\left(t\right)+n\left(t\right)*g\left(t\right)$

$=\left(\underset{i=0}{\overset{N_p-1}{\mathrm{\Σ}}}\underset{j=0}{\overset{N_f-1}{\mathrm{\Σ}}}\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{\mathrm{\α}}_{l}p(t-{\mathrm{\τ}}_l-{\mathrm{iN}}_f{T}_f-{\mathrm{jT}}_f)\right)+\underset{i=0}{\overset{N_s-1}{\mathrm{\Σ}}}\underset{j=0}{\overset{N_f-1}{\mathrm{\Σ}}}{b}_i\underset{l=1}{\overset{L}{\mathrm{\Σ}}}{\mathrm{\α}}_{l}p(t-{\mathrm{\τ}}_l-{\mathrm{iN}}_f{T}_f-{\mathrm{jT}}_f-N_p{N}_f{T}_f)+w\left(t\right)---\left(5\right)$

N (t) expression AWGN noise, variance is σ ², bilateral power spectral density is designated as N ₀/ 2,

G (t) represents band pass filter, and passband is [f _c-B/2, f _c+ B/2], play the filtering out-of-band noise, can not cause distortion to receiving waveform,

The band limit output of w (t) expression n (t), auto-correlation function is R _w(τ)=BN ₀Sinc (B τ) cos (2 π f _cτ).

The impulse ultra-wideband signal detection method that the present invention is based on compression sampling comprises the steps:

Step 100: obtain the compression sampling sequence that frequency pilot sign receives signal, that is, measure by the compression that continuous multiframe is received signal, obtain the sample sequence that frequency pilot sign receives signal, as the template sequence of follow-up coherent detection;

In specific implementation process, the present invention is divided into without feedback loop and two kinds of working methods of feedback loop is arranged, feedback loop is with the result partial feedback of the described digital back-end processing components loop to described measured waveform generator, the measured waveform generator produces new measured waveform according to feedback information, the template sequence of follow-up coherent detection obtain and treat obtaining of correlated series, all use the new measured waveform that produces to obtain.Wherein:

During without feedback loop: establish the required M of a symbol judgement measured value and share on the measurement to continuous D frame signal, make N _D=N _f/ D, then a total N _pN _DCriticize compression and measure sequence, be designated as

y _p[n]＝Фr _prj+Фw _p[n]，1≤n≤N _pN _D (6)

Wherein Ф is that matrix, w are measured in compression _p[n] is that variance is σ ²The virtual sample sequence of WGN, r _PrjThe reception signal that is a pulse drops on the virtual sample sequence that compresses in the drop shadow spread.

When feedback loop was arranged, the process that obtains the compression sampling sequence was as follows:

Step 110: obtain the compression sampling sequence that first's frequency pilot sign receives signal, that is, this stage is used N _pBefore in the individual frequency pilot sign

Individual, compress to received signal measurement with the measurement matrix Ф of completely random.Total Criticize compression and measure sequence, be designated as

$y_{p_1}\left[n\right]={\mathrm{\Φr}}_{\mathrm{prj}}+{\mathrm{\Φw}}_{p_1}\left[n\right],1\≤n\≤N_{p_1}{N}_D---\left(7\right)$

Wherein

That variance is σ ²The virtual sample sequence of WGN.

Step 120: obtain subspace estimation, that is, use Carry out the estimation to signal subspace H, its flow process is complete to be described below:

Step 121: initialization.

Make dictionary matrix V=Ф Ψ, v _iThe i row of representing matrix V, i.e. i atom in the dictionary; Measure residual error

Estimated result L=[], iterations t=1.

Step 122: obtain the sequence number of the atom that mates most with residual error, that is: from dictionary, seek the sequence number of the atom that mates most with residual error.

$l_t=\underset{i=\mathrm{1,2},...,N}{\arg \max }\frac{|<e_{t-1},v_i>|}{\left|\right|v_i\left|\right|}---\left(8\right)$

Step 123: upgrade residual error.

$e_t=e_{t-1}-\frac{|<e_{t-1},v_{l_t}>|}{{\left|\right|v_{l_t}\left|\right|}^2}{v}_{l_t}---\left(9\right)$

Step 124: judge whether to stop.

If t＞K ', then iteration stops; Otherwise upgrade t=t+1, with current search to the atom sequence number be increased in the estimated result i.e. L=[L, l _t], and return step 121 initialization.

Step 125: obtain the estimation of subspace, the estimated result that obtains signal subspace according to L is

Step 130: obtain the subspace compression sampling sequence that the second portion frequency pilot sign receives signal.

Digital back-end reception ﹠ disposal assembly 3 feeds back to measured waveform generator 2 with it after the estimation that obtains H, measured waveform generator 2 produces new measured waveform according to the signal subspace estimated information that obtains, namely Wherein G is random matrix.Use new measured waveform to remaining in the frequency pilot sign

Individual reception signal compresses measurement, can obtain altogether

Criticize compression and measure sequence, be designated as

$y_{p_2}\left[n\right]=\stackrel{\&OverBar;}{\mathrm{\&Phi;}}{r}_{\mathrm{prj}}+\stackrel{\&OverBar;}{\mathrm{\&Phi;}}{w}_{p_2}\left[n\right],1\&le;n\&le;N_{p_2}{N}_D---\left(10\right)$

Wherein

That variance is σ ²The virtual sample sequence of WGN.The mean sequence of formula (10)

To in subspace compression domain coherent detection subsequently, be used as correlate template.

Step 200: obtain the compression sampling sequence that data symbol receives signal, that is, measure the compression sampling sequence that matrix obtains data symbol reception signal according to the compression that the measured waveform generator produces;

During without feedback loop, j (1≤j≤N _s) the total N of individual data symbol _DCriticize compression and measure sequence, be designated as

y _s|j[n]＝b _jФr _prj+Фw _s[n]，1≤n≤N _D，1≤j≤N _s (11)

W wherein _sIt is σ that [n] data symbol receives the variance that adds in the signal ²The virtual sample sequence of WGN.When feedback loop is arranged, the result partial feedback of described digital back-end processing components during to described measured waveform generator, is measured matrix with the subspace compression

To N _sThe compression of the j in the individual data symbol reception signal is measured sequence and (is total to N _DBatch) be

$y_{s|j}\left[n\right]=b_j\stackrel{\&OverBar;}{\mathrm{\&Phi;}}{r}_{\mathrm{prj}}+\stackrel{\&OverBar;}{\mathrm{\&Phi;}}{w}_{s|j}\left[n\right],1\&le;n{\&le;N}_D,1\&le;j\&le;N_s---\left(12\right)$

W wherein _S|j[n] is that variance is σ ²The virtual sample sequence of WGN.Mean sequence shown in the formula (12)

Will be as reception burst to be correlated with in subspace compression domain correlation detector subsequently.

Step 300: obtain symbol judgement: by Hypothesis Testing Problem derivation coherent detection judgement amount, obtain symbol judgement according to detecting judgement amount.

Specific as follows:

During without feedback loop, to b _jJudgement be following Hypothesis Testing Problem:

$\left\{\begin{array}{cc}{H}_0:y_s\left[n\right]=-{\mathrm{\&Phi;r}}_{\mathrm{prj}}+\mathrm{\&Phi;}{w}_s\left[n\right],& (b_j=-1)\\ H_1:y_s\left[n\right]={\mathrm{\&Phi;r}}_{\mathrm{prj}}+\mathrm{\&Phi;}{w}_s\left[n\right],& (b_j=1)\end{array}\right.---\left(13\right)$

Based on generalized likelihood-ratio test, can derive following coherent detection judgement amount

$T\left(y_s\right)={(\mathrm{\&Phi;}{r}_{\mathrm{prj}}+\mathrm{\&Phi;}{\stackrel{\&OverBar;}{w}}_p)}^T(b_j\mathrm{\&Phi;}{r}_{\mathrm{prj}}+\mathrm{\&Phi;}{\stackrel{\&OverBar;}{w}}_s)\left\{\begin{array}{c}>0\&RightArrow;{\mathrm{<figure-callout\; id="1"\; label="H"\; filenames="HSA00000240840900011.png"\; state="\{\{state\}\}">H</figure-callout>}}_1\\ <0\&RightArrow;H_0\end{array}\right.---\left(14\right)$

In the following formula, correlate template is the mean sequence of formula (6), treats that correlated series is the mean sequence of formula (11).

Feedback loop is arranged, can derive the coherent detection judgement amount and be

$T\left(y_s\right)={(\stackrel{\&OverBar;}{\mathrm{\&Phi;}}{r}_{\mathrm{prj}}+\stackrel{\&OverBar;}{\mathrm{\&Phi;}}{\stackrel{\&OverBar;}{w}}_{p_2})}^T(\stackrel{\&OverBar;}{\mathrm{\&Phi;}}{r}_{\mathrm{prj}}+\stackrel{\&OverBar;}{\mathrm{\&Phi;}}{\stackrel{\&OverBar;}{w}}_s)\left\{\begin{array}{c}>0\&RightArrow;{\mathrm{<figure-callout\; id="1"\; label="H"\; filenames="HSA00000240840900011.png"\; state="\{\{state\}\}">H</figure-callout>}}_1\\ <0\&RightArrow;H_0\end{array}\right.---\left(15\right)$

Correlate template is the mean sequence of formula (10), treats that correlated series is the mean sequence of formula (12).

Above content is the further description of the present invention being done in conjunction with concrete preferred implementation, can not assert that implementation of the present invention is confined to these explanations.For the general technical staff of the technical field of the invention, without departing from the inventive concept of the premise, can also make some simple deduction or replace, all should be considered as belonging to protection scope of the present invention.

Claims (5)

1. impulse ultra-wideband signal detection method based on compression sampling, it is characterized in that, impulse ultra-wideband signal compression sampling receiving system comprises the multi-channel parallel sampling unit that divides a plurality of passages to sample to described impulse ultra-wideband signal, send respectively the measured waveform generator of measured waveform to each passage of described multi-channel parallel sampling unit, receive the digital back-end processing components through the measured value of described multi-channel parallel sampling unit sampling, each passage carries out linear projection according to the measured waveform that described measured waveform generator produces to described impulse ultra-wideband signal respectively, comprises the steps:

Obtain the compression sampling sequence that frequency pilot sign receives signal: measure by the compression that continuous multiframe is received signal, obtain the sample sequence that frequency pilot sign receives signal, as the template sequence of follow-up coherent detection;

Obtain the compression sampling sequence that data symbol receives signal: measure by the compression that continuous multiframe is received signal, obtain the compression sampling sequence that data symbol receives signal, as the correlated series for the treatment of of follow-up coherent detection;

Obtain symbol judgement: by coherent detection, obtain symbol judgement.

2. the impulse ultra-wideband signal detection method based on compression sampling according to claim 1 is characterized in that, establishing used number of pilot symbols is N _p, the pulse that each symbol uses repeats transmission times and is N _f, M required measured value of symbol judgement shared on the measurement to continuous D frame signal, makes N _D=N _fD, then a total N _pN _DCriticize compression and measure sequence, it is y that n criticizes sequence _p[n]=Φ r _Prj+ Φ w _p[n], wherein: Φ is that matrix, r are measured in corresponding compression _PrjThe reception signal that is a pulse drops on the virtual sample sequence that compresses in the drop shadow spread, w _p[n] is that variance is σ ²The virtual sample sequence of WGN.

3. the impulse ultra-wideband signal detection method based on compression sampling according to claim 1, it is characterized in that, described impulse ultra-wideband signal compression sampling receiving system also comprises feedback loop, described feedback loop with the result partial feedback of described digital back-end processing components to described measured waveform generator, the measured waveform generator produces new measured waveform according to feedback information, the template sequence of follow-up coherent detection obtain and treat obtaining of correlated series, all use the new measured waveform that produces to obtain.

4. the impulse ultra-wideband signal detection method based on compression sampling according to claim 3, it is characterized in that, for the situation that feedback loop is arranged, in the compression sampling sequence step of obtaining frequency pilot sign reception signal, frequency pilot sign receives signal and is divided into first and second portion, comprises the steps:

Obtain the compression sampling sequence that first's frequency pilot sign receives signal: transmitting terminal sends

Individual frequency pilot sign compresses measurement with the measurement matrix Φ of completely random to corresponding reception signal, and establishing used number of pilot symbols is N _p, the pulse that each symbol uses repeats transmission times and is N _f, M required measured value of symbol judgement shared on the measurement to continuous D frame signal, makes N _D=N _fD, total Criticize the compression measured value, n criticizes sequence and is designated as Wherein

Be the virtual sample sequence that this part frequency pilot sign receives noise in the signal, Φ is that matrix is measured in corresponding compression;

Obtain subspace estimation: the compression sampling sequence that adopts first's frequency pilot sign to receive signal is estimated signal subspace H;

Obtain the subspace compression sampling sequence that the second portion frequency pilot sign receives signal: the digital back-end processing components is after the estimation that obtains signal subspace H, it is fed back to the measured waveform generator, the measured waveform generator produces new measured waveform according to the signal subspace estimated information that obtains, namely

Wherein G is random matrix, to remaining in the frequency pilot sign

Individual reception signal compresses measurement, can obtain altogether

Criticize compression and measure sequence, be designated as

Wherein The virtual sample sequence that this part frequency pilot sign receives noise in the signal,

Mean sequence

Will be as the template sequence of coherent detection afterwards, r _PrjThe reception signal that is a pulse drops on the virtual sample sequence that compresses in the drop shadow spread.

5. described impulse ultra-wideband signal detection method based on compression sampling according to claim 4, it is characterized in that, in obtaining the subspace estimation step, the compression sampling sequence that adopts first's frequency pilot sign to receive signal is carried out signal subspace and is estimated that the first's frequency pilot sign that adopts has p ₁Individual, it is total that it receives compression measurement sequence corresponding to signal

Batch, n criticizes sequence and is

Wherein

That variance is σ ²The virtual sample sequence of WGN, the signal subspace algorithm for estimating that adopts comprises the steps:

Step 121, with the compression sampling sequence initialization of subspace: dictionary matrix V=Φ Ψ, measure residual error

Estimated result L=[] and iterations t=1, v _iThe i row of representing matrix V;

Step 122, obtain the sequence number with the matched atoms of residual error: from the dictionary matrix, seek the sequence number with the matched atoms of residual error, $l_t=\underset{i=\mathrm{1,2},...,N}{\arg \max }\frac{|<e_{t-1},v_i>|}{\left|\right|v_i\left|\right|};$

The coefficient of step 123, the current atom of acquisition also upgrades residual error: upgrade residual error

Step 124, judge whether to stop: establishing maximum iteration time is K', if t〉K', then iteration stops, and skips to step 125; Otherwise upgrade t=t+1, with current search to the atom sequence number be increased in the estimated result i.e. L=[L, l _t], and rebound step 122;

The estimation of step 125, acquisition subspace: the estimated result that obtains signal subspace according to L is $H=[{{\mathrm{\&psi;}}_l}_1,{{\mathrm{\&psi;}}_l}_2,...{{\mathrm{\&psi;}}_l}_{K\&prime;}],$

Wherein, Ψ represents sparse expression matrix, ψ _lThe l row of expression sparse expression matrix Ψ,

Expression is owned

Criticize the mean sequence that receives the compression measurement sequence of signal for first's frequency pilot sign, |＜e _T-1, v _i| the residual error e when representing the t-1 time iteration _T-1I row v with matrix V _iCarry out inner product operation, then take absolute value, || v _i|| the i row v of representing matrix V _iTwo norm value,

Residual error e when representing the t-1 time iteration _T-1L with matrix V _tRow

Carry out inner product operation, then take absolute value.

Images