CN106130571B  A kind of signal sampling method of reseptance and device based on frequency domain and airspace compressed sensing  Google Patents
A kind of signal sampling method of reseptance and device based on frequency domain and airspace compressed sensing Download PDFInfo
 Publication number
 CN106130571B CN106130571B CN201610435298.8A CN201610435298A CN106130571B CN 106130571 B CN106130571 B CN 106130571B CN 201610435298 A CN201610435298 A CN 201610435298A CN 106130571 B CN106130571 B CN 106130571B
 Authority
 CN
 China
 Prior art keywords
 signal
 compressed sensing
 attenuator
 frequency
 radiofrequency signal
 Prior art date
Links
 238000005070 sampling Methods 0 abstract claims description title 24
 239000011159 matrix materials Substances 0 abstract claims description 29
 238000004364 calculation methods Methods 0 abstract claims description 23
 229920002574 CR39 Polymers 0 abstract description 16
 238000007906 compression Methods 0 claims description 15
 238000000034 methods Methods 0 claims description 5
 230000002829 reduced Effects 0 abstract description 3
 238000003860 storage Methods 0 abstract description 3
 238000004458 analytical methods Methods 0 abstract description 2
 239000003637 basic solution Substances 0 abstract description 2
 238000004891 communication Methods 0 abstract description 2
 230000001595 contractor Effects 0 claims description 2
 230000001702 transmitter Effects 0 abstract description 2
 230000000875 corresponding Effects 0 claims 1
 230000001603 reducing Effects 0 abstract 1
 241001442055 Vipera berus Species 0 description 8
 230000004304 visual acuity Effects 0 description 3
 238000007792 addition Methods 0 description 2
 230000018109 developmental process Effects 0 description 2
 238000005516 engineering processes Methods 0 description 2
 238000001914 filtration Methods 0 description 2
 238000001228 spectrum Methods 0 description 2
 206010008531 Chills Diseases 0 description 1
 241000018646 Pinus brutia Species 0 description 1
 241001282153 Scopelogadus mizolepis Species 0 description 1
 239000000969 carrier Substances 0 description 1
 239000004567 concrete Substances 0 description 1
 238000010276 construction Methods 0 description 1
 238000007796 conventional methods Methods 0 description 1
 230000000694 effects Effects 0 description 1
 238000002156 mixing Methods 0 description 1
 230000004048 modification Effects 0 description 1
 238000006011 modification Methods 0 description 1
 230000004224 protection Effects 0 description 1
Abstract
Description
Technical field
The present invention relates to radar, electronic reconnaissance and wideband communication signal process field, more particularly to it is a kind of based on frequency domain and The signal sampling method of reseptance and device of airspace compressed sensing.
Background technique
Conventional electrical reconnaissance system is passed through frequently with phaseinterferometer or spatial spectrum anglemeasuring method to aerial/Ground emitter signals Carry out direction of arrival (DOA) estimation and/or positioning.Common workflow: Radar emitter radiofrequency signal is defeated through aerial array Enter, by receiver front end clipping, low noise, filtering and mixing, (on) become intermediatefreuqncy signal output after down coversion, ADC is adopted Channelizing is done after sample, is become baseband signal and is further processed.
With the fast development of Radar Technology, multifunctional novel radar is continued to bring out, radar transmitter radiofrequency signal frequency Rate is higher and higher, and bandwidth is more and more wider, very high to electronic reconnaissance system structure serious challenge, such as ADC sample rate, And the resolving power of ADC can not too low (general number of significant digit be not less than 8bit), because of frequencymeasurement accuracy and constant false alarm threshold value The number of significant digit of analogdigital converter is both depended in setting to noise Estimation of Mean precision, resolving power is too low, and frequencymeasurement accuracy is deteriorated, It is unfavorable for the DOA estimation based on frequency measurement system.Falsealarm probability increases, and is unfavorable for radar emitter signal detection.But ADC technology Development can not infinitely cater to this demand, and cannot usually make progress simultaneously between sample rate and resolving power because This can bring the great expense incurred in cost, and this requires new method, new principle must be used to solve electronic reconnaissance system physics reality Existing feasibility and complexity problem.
As long as although bandpass sample theory show ADC sample rate meet be more than or equal to twice of signal bandwidth, without with respect to Signal carrier frequency, so that it may restore original signal without distortions, but for a wideband electronic reconnaissance system, this condition is put Pine does not tackle the problem at its root.
If directly sampled to radiofrequency signal, so that it may solve high band, big bandwidth, low noise wideband radar signal and detect Rigors, the compressive sensing theories such as receipts, sampling, storage provide theoretical foundation for direct RF sampling, its principle is benefit With radar emitter signal in the sparsity of frequency domain, being realized by construction frequency domain measurement matrix can be reconstructed with compared with low sampling rate Frequencyregion signal purpose, the realization for wideband electronic reconnaissance system provide feasible scheme.However, this when of receiver channel Complexity is still quite high, and conditional electronic reconnaissance system antenna receiving unit and receiver channel correspond, with antenna number Mesh increases significantly, and receiver channel number increases therewith, and system cost and complexity are got worse, for example spatial spectrum angle measurement is using big Scale linear array, angle measurement accuracy is quite high, but brings receiver design excessively challenge.
Summary of the invention
In view of this, the present invention provides a kind of signal sampling method of reseptance and dress based on frequency domain and airspace compressed sensing It sets, reduces digital receiver port number while ADC sample rate can be reduced, greatly reduce system design cost.
The invention has the following beneficial effects:
(1) present invention is right using frequency domain compressed sensing matrix and the concatenated form of spatial domain compressed sensing calculation matrix The compression of radiofrequency signal progress frequency domain and spatial domain, wherein after the compression of frequency domain, ADC sample rate becomes adopting originally M/N times of sample rate has saved system design cost and system power dissipation, provides feasibility side for signal storage and further analysis Case；After spatial domain is compressed, digital receiver port number becomes original m/n times, significantly reduces system complexity.
Detailed description of the invention
Fig. 1 is the flow chart of sampling method of reseptance of the invention.
Fig. 2 is the realization device schematic diagram of frequency domain compressed sensing calculation matrix of the invention.
Fig. 3 is the realization device schematic diagram of spatial domain compressed sensing calculation matrix of the invention.
Specific embodiment
The present invention will now be described in detail with reference to the accompanying drawings and examples.
Radar emitter signal is that compressed sensing building low complex degree receiver brings hope in the sparsity of spatial domain, Theoretical basis has been established for building low cost wide having electronic reconnaissance system, has realized that highprecision DOA estimates using less receiver channel Meter.
Signal is the foundation for constructing compressed sensing calculation matrix dimension size, the determination of degree of rarefication in the degree of rarefication of not same area It needs combined with concrete application background, determination method can be achieved by the prior art.We assume that radar emitter signal Frequency domain and airspace degree of rarefication it is known that the present invention proposes a kind of signal sampling recipient based on frequency domain and airspace compressed sensing Method, as shown in Figure 1, including the following steps:
Step 1 is directed to the received every radiofrequency signal all the way of n antenna element in reception system, carries out frequency domain pressure respectively Contracting perception processing, the road the n radiofrequency signal after obtaining frequency compression；Wherein, it uses with row vector as M, column vector is the frequency of N Domain compressed sensing calculation matrix carries out frequency compression, and M ≈ 3K, K are degree of rarefication of the target in frequency domain；N is number of samples；
The road n radiofrequency signal after step 2, the frequency compression obtained for step 1 is sent into spatial domain compressed sensing and measures square Spatial domain compression is carried out in battle array, is obtained by the compressed road the m radiofrequency signal of spatial domain；Wherein, the spatial domain compressed sensing The row vector of calculation matrix is m, and column vector n, m ≈ 3H, H are degree of rarefication of the target in spatial domain；
Step 3 is sampled to by the compressed road the m radiofrequency signal of spatial domain respectively using m Sampling Machine；
Step 4 respectively receives the road the m radiofrequency signal after sampling using m receiver, is achieved in adopting for signal Sample receive process.
In Fig. 1, electronic reconnaissance system covering frequence range is f_{L}f_{H}, receive system and use n antenna element, receive N road radiofrequency signal x (t).Assuming that the frequencydomain sparse degree of radar emitter signal is K, then frequency domain compressed sensing calculation matrix A ∈ R^{M} ^{×N}Row vector meet M ≈ 3K, due to K < < N, M < < N, N are to exist to the radiofrequency signal for not carrying out frequency domain compression The number of samples of ADC in unit time, M are the radiofrequency signal compressed to frequency domain within the same unit time after compression sampling ADC number of samples.The radiofrequency signal of the road n antenna receiving unit input is x (t), then passes through frequency domain compressed sensing calculation matrix A Afterwards, become n road radiofrequency signal y (t)=Ax (t), ADC sample rate is M/N times of original sample rate, and sample rate greatly reduces.
After spatial domain compressed sensing calculation matrix, n road radiofrequency signal y (t) is surveyed using spatial domain compressed sensing Moment matrix B ∈ R^{m×n}Afterwards, become the road m radiofrequency signal output z (t)=By (t)=BAx (t), wherein m ≈ 3H, H indicate that target exists The degree of rarefication of spatial domain.Due to H < < n, then m < < n, digital receiver port number becomes original m/n times at this time, subtracts significantly Receiver channel number is lacked, Receiver Complexity is greatly lowered.Signal reconstruction in diagram 1 realizes frequency domain sparse signal weight Structure, channelization scheme are realized by signal reconstruction, common channelization scheme in conventional electrical reconnaissance system are saved, by number Low noise, filtering, (on) eventually become baseband signal subsequent processing after down coversion.
Signal reconstruction is instead of channelization scheme in conventional electrical reconnaissance system；Channelizing side in conventional electrical reconnaissance system Case is substantially that timedomain sampling signal is transformed into frequency domain, frequency sweep is carried out to frequencyregion signal, to determine which emitter Signals are located at In a subband.The receiver channel problem of inconsistency encountered in conventional electrical reconnaissance system has obtained basic solution.Because Traditional electronic reconnaissance system mostly uses greatly multipassage design scheme to realize that a certain function, usual process are antenna radiofrequency signal warps It after crossing broadband receiver, is exported after radiofrequency signal is transformed into intermediatefreuqncy signal, at this time radiofrequency signal bandwidth drops to highspeed ADC When the acceptable range of sample rate, starting to execute ADC sampling, this is conventional method, actually in conventional scheme, Duo Gejie It receives channel and uses a large amount of analog devices, be not exclusively digital receiver, so bringing the amplitude phase unbalance problem of signal.And The direct sampling plan of new radio frequency then uses alldigital receiver, and multichannel sensor gain and phase uncertainties problem, other links are not present The error of introducing can be ignored, such as sample quantization error etc..
The present invention also provides the realization of frequency domain compressed sensing calculation matrix and spatial domain compressed sensing calculation matrix dresses It sets, as shown in Fig. 2, realizing the device of frequency domain compressed sensing calculation matrix, including M × N number of attenuator and M+1 addition Device；Indicate that attenuator, ∑ () indicate adder；The radiofrequency signal on road where the input terminal of each attenuator receives；M× N number of attenuator is divided into M group, and the output of every group of N number of attenuator terminates an adder, and M group attenuator corresponds to M addition Device, another adder connect the output end in the M adder.Ith tunnel radiofrequency signal x_{i}(t) be input to parallel M × N number of by In the structured LDPC binary measure matrix that attenuator array is constituted, " focusing " of output signal passes through first adder of diagram It realizes, second adder realizes frequency domain compressed sensing calculation matrix totally " focusing ", y_{i}(t)=Ax_{i}It (t) is the ith road radio frequency Signal output.
As shown in figure 3, realizing the device of spatial domain compressed sensing calculation matrix, including m × n attenuator and m add Musical instruments used in a Buddhist or Taoist mass；Indicate that attenuator, ∑ () indicate adder；M × n attenuator is divided into m group, and every group has n attenuator；Its In, n attenuator in any one group, which corresponds, receives the road n radiofrequency signal after frequency compression, every group of attenuator Output end connects an adder jointly.M is input to by the n road radiofrequency signal y (t) that frequency domain compressed sensing calculation matrix exports parallel In × n attenuator, the row vector of LDPC binary measure matrix B is made of n attenuator, after " multiplyadd " operation, exports m Road analog radiofrequency signal z (t).Analog radio frequency after frequency domain compressed sensing calculation matrix and airspace compressed sensing calculation matrix Signal z (t), after the road m ADC sampling, not only ADC sample rate is reduced, but also receiver channel number decreases, still DOA estimated accuracy is unaffected.
In conclusion the above is merely preferred embodiments of the present invention, being not intended to limit the scope of the present invention. All within the spirits and principles of the present invention, any modification, equivalent replacement, improvement and so on should be included in of the invention Within protection scope.
Claims (3)
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201610435298.8A CN106130571B (en)  20160607  20160607  A kind of signal sampling method of reseptance and device based on frequency domain and airspace compressed sensing 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201610435298.8A CN106130571B (en)  20160607  20160607  A kind of signal sampling method of reseptance and device based on frequency domain and airspace compressed sensing 
Publications (2)
Publication Number  Publication Date 

CN106130571A CN106130571A (en)  20161116 
CN106130571B true CN106130571B (en)  20181221 
Family
ID=57470866
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201610435298.8A CN106130571B (en)  20160607  20160607  A kind of signal sampling method of reseptance and device based on frequency domain and airspace compressed sensing 
Country Status (1)
Country  Link 

CN (1)  CN106130571B (en) 
Citations (3)
Publication number  Priority date  Publication date  Assignee  Title 

CN103117818A (en) *  20130116  20130522  南京邮电大学  Broadband spectrum sensing method based on spacefrequency joint compressed sensing 
US9015007B2 (en) *  20101213  20150421  Southwest Research Institute  Sensor array processor with multichannel reconstruction from random array sampling 
CN105158735A (en) *  20150603  20151216  西安电子科技大学  Space frequency twodimensional spectrum estimation method based on compressed sampling array 

2016
 20160607 CN CN201610435298.8A patent/CN106130571B/en active IP Right Grant
Patent Citations (3)
Publication number  Priority date  Publication date  Assignee  Title 

US9015007B2 (en) *  20101213  20150421  Southwest Research Institute  Sensor array processor with multichannel reconstruction from random array sampling 
CN103117818A (en) *  20130116  20130522  南京邮电大学  Broadband spectrum sensing method based on spacefrequency joint compressed sensing 
CN105158735A (en) *  20150603  20151216  西安电子科技大学  Space frequency twodimensional spectrum estimation method based on compressed sampling array 
Also Published As
Publication number  Publication date 

CN106130571A (en)  20161116 
Similar Documents
Publication  Publication Date  Title 

Brown  The Galactic Radio Spectrum Between 130 and 26OOkHz  
Shepherd  Radio wave loss deviation and shadow loss at 900 MHz  
Czechowsky et al.  The mobile SOUSY Doppler radar: Technical design and first results  
JP2007529000A (en)  Direction detection  
EP1626287A1 (en)  Technique for estimating signal arrival direction not by utilizing eigenvalue decomposition and reception beam shaper  
Perrenoud  The computercontrolled solar radio spectrometer ‘IKARUS’  
Ott et al.  Urban pathloss characteristics at 820 MHz  
Hoffman et al.  Propagation observations at 3.2 millimeters  
Reid et al.  First VHF radar measurements of mesopause summer echoes at mid‐latitudes  
Cottony et al.  A highresolution rapidscan antenna  
US8548032B2 (en)  Method of detecting a radar signal, radar detection module, and wireless transceiver including the same  
CN101149429A (en)  Array mutual coupling calibration and source direction estimation method suitable for uniform circular array  
CN101105525A (en)  Pure phase type broad frequency band microwave radiation source direction finding system and method  
Lo et al.  A new approach for estimating indoor radio propagation characteristics  
Nuteson et al.  Performance characterization of FPGA techniques for calibration and beamforming in smart antenna applications  
US6700537B2 (en)  Method for calibrating a wideband direction finding system  
Lundgren  A study of mutual coupling effects on the direction finding performance of ESPRIT with a linear microstrip patch array using the method of moments  
Zhang et al.  Design of the second generation microwave humidity sounder (MWHSII) for Chinese meteorological satellite FY3  
Hamidi et al.  Evaluation of signal to noise ratio (SNR) of log periodic dipole antenna (LPDA)  
EP2717481A2 (en)  System and Method for determining an interferer transmitting an interfering signal  
CN103517285A (en)  Method and system for coverage performance evaluation of antenna  
US8723723B2 (en)  Dual mode ground penetrating radar (GPR)  
CN102879781B (en)  Distributed synthetic aperture radiometer array imaging method and system  
Hamidi et al.  Modification and Performance of Log Periodic Dipole Antenna  
Grubb et al.  A new general purpose high performance HF Radar 
Legal Events
Date  Code  Title  Description 

DD01  Delivery of document by public notice 
Addressee: Wu Riheng Document name: Notification of Acceptance of Patent Application 

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