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 PDF

Info

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
attenuator
compressed sensing
radiofrequency signal
frequency
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.)
Active
Application number
CN201610435298.8A
Other languages
Chinese (zh)
Other versions
CN106130571A (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.)
Wenjing College Yantai University
Original Assignee
Wenjing College Yantai 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 Wenjing College Yantai University filed Critical Wenjing College Yantai University
Priority to CN201610435298.8A priority Critical patent/CN106130571B/en
Publication of CN106130571A publication Critical patent/CN106130571A/en
Application granted granted Critical
Publication of CN106130571B publication Critical patent/CN106130571B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/0003Software-defined radio [SDR] systems, i.e. systems wherein components typically implemented in hardware, e.g. filters or modulators/demodulators, are implented using software, e.g. by involving an AD or DA conversion stage such that at least part of the signal processing is performed in the digital domain
    • H04B1/0007Software-defined radio [SDR] systems, i.e. systems wherein components typically implemented in hardware, e.g. filters or modulators/demodulators, are implented using software, e.g. by involving an AD or DA conversion stage such that at least part of the signal processing is performed in the digital domain wherein the AD/DA conversion occurs at radiofrequency or intermediate frequency stage
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/06Receivers
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0837Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
    • H04B7/0842Weighted combining
    • H04B7/086Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a kind of signal sampling method of reseptance and device based on frequency domain and airspace compressed sensing, the present invention can reduce radar, the high-speed ADC signal sampling acceptance rate of the broadbands such as electronic reconnaissance and ultra-wideband communications application field, its principle is using the emitter Signals such as radar or ultra-wideband communications transmitter in frequency domain and airspace sparsity feature, devise tandem compressed sensing calculation matrix, it is this to design the sample rate for not only reducing time domain broadband radiofrequency signal ADC, system design cost and system power dissipation are saved, feasible scheme is provided for signal storage and further analysis, and receiver channel number is also greatly reduced, reduce radar, electronic reconnaissance and ultra-wideband communication system design complexities.Meanwhile signal reconstruction is instead of conventional channelization scheme.The receiver channel problem of inconsistency encountered in conventional design has obtained basic solution.

Description

A kind of signal sampling method of reseptance and device based on frequency domain and airspace compressed sensing
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 phase-interferometer or spatial spectrum angle-measuring 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 intermediate-freuqncy 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 frequency-measurement accuracy and constant false alarm threshold value The number of significant digit of analog-digital converter is both depended in setting to noise Estimation of Mean precision, resolving power is too low, and frequency-measurement accuracy is deteriorated, It is unfavorable for the DOA estimation based on frequency measurement system.False-alarm 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 Frequency-region 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 high-precision 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 fL-fH, receive system and use n antenna element, receive N road radiofrequency signal x (t).Assuming that the frequency-domain sparse degree of radar emitter signal is K, then frequency domain compressed sensing calculation matrix A ∈ RM ×NRow 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 ∈ Rm×nAfterwards, 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 time-domain sampling signal is transformed into frequency domain, frequency sweep is carried out to frequency-region signal, to determine which emitter Signals are located at In a sub-band.The receiver channel problem of inconsistency encountered in conventional electrical reconnaissance system has obtained basic solution.Because Traditional electronic reconnaissance system mostly uses greatly multi-passage design scheme to realize that a certain function, usual process are antenna radio-frequency signal warps It after crossing broadband receiver, is exported after radiofrequency signal is transformed into intermediate-freuqncy signal, at this time radiofrequency signal bandwidth drops to high-speed 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 all-digital 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.I-th tunnel radiofrequency signal xi(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 ", yi(t)=AxiIt (t) is the i-th 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 " multiply-add " operation, exports m Road analog radio-frequency 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)

1. a kind of signal sampling method of reseptance based on frequency domain and airspace compressed sensing, which comprises the steps of:
Step 1, respectively progress frequency domain compression sense received per radiofrequency signal all the way for n antenna element in reception system Know processing, the road the n radiofrequency signal after obtaining frequency compression;Wherein, it uses with row vector as M, column vector is the frequency domain pressure of N Contracting perception 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 calculation matrix Spatial domain compression is carried out, is obtained by the compressed road the m radiofrequency signal of spatial domain;Wherein, the spatial domain compressed sensing measurement The row vector of 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, and the sampling for being achieved in signal connects Receipts process.
2. a kind of dress for realizing the frequency domain compressed sensing calculation matrix in signal sampling method of reseptance described in claim 1 It sets, which is characterized in that including M × N number of attenuator and M+1 adder;Frequency where the input terminal of each attenuator receives The corresponding radiofrequency signal of rate domain compressed sensing calculation matrix;M × N number of attenuator is divided into M group, every group of N number of attenuator it is defeated An adder is terminated out, and M group attenuator corresponds to M adder, another adder connects the output in the M adder End.
3. a kind of dress for realizing the spatial domain compressed sensing calculation matrix in signal sampling method of reseptance described in claim 1 It sets, which is characterized in that including m × n attenuator and m adder;M × n attenuator is divided into m group, and every group has n Attenuator;Wherein, n attenuator in each group, which corresponds, receives the road the n radiofrequency signal after frequency compression, every group of decaying The output end of device connects an adder jointly.
CN201610435298.8A 2016-06-07 2016-06-07 A kind of signal sampling method of reseptance and device based on frequency domain and airspace compressed sensing Active CN106130571B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610435298.8A CN106130571B (en) 2016-06-07 2016-06-07 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) 2016-06-07 2016-06-07 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) 2016-11-16
CN106130571B true CN106130571B (en) 2018-12-21

Family

ID=57470866

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610435298.8A Active CN106130571B (en) 2016-06-07 2016-06-07 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)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10630511B2 (en) * 2017-12-13 2020-04-21 Hon Hai Precision Industry Co., Ltd. Methods and devices for channel estimation
CN111464184B (en) * 2020-04-22 2022-03-15 电子科技大学 Time-interleaved ADC based on compressed sensing
CN112731486B (en) * 2020-12-24 2023-07-04 深圳大学 Signal angle and signal frequency estimation method, device, equipment and storage medium
CN112929303A (en) * 2021-01-21 2021-06-08 哈尔滨工程大学 Broadband compressed sensing direction-finding method of double-chain quantum charged system search mechanism
CN113503954B (en) * 2021-05-27 2022-04-22 南京大学 phi-OTDR signal measurement and noise suppression method based on multi-domain compressed sensing
CN113784381B (en) * 2021-10-13 2024-02-27 厦门大学 Ultra-wideband communication system measurement matrix design method based on chaotic compressed sensing

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103117818A (en) * 2013-01-16 2013-05-22 南京邮电大学 Broadband spectrum sensing method based on space-frequency joint compressed sensing
US9015007B2 (en) * 2010-12-13 2015-04-21 Southwest Research Institute Sensor array processor with multichannel reconstruction from random array sampling
CN105158735A (en) * 2015-06-03 2015-12-16 西安电子科技大学 Space frequency two-dimensional spectrum estimation method based on compressed sampling array

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9015007B2 (en) * 2010-12-13 2015-04-21 Southwest Research Institute Sensor array processor with multichannel reconstruction from random array sampling
CN103117818A (en) * 2013-01-16 2013-05-22 南京邮电大学 Broadband spectrum sensing method based on space-frequency joint compressed sensing
CN105158735A (en) * 2015-06-03 2015-12-16 西安电子科技大学 Space frequency two-dimensional spectrum estimation method based on compressed sampling array

Also Published As

Publication number Publication date
CN106130571A (en) 2016-11-16

Similar Documents

Publication Publication Date Title
CN106130571B (en) A kind of signal sampling method of reseptance and device based on frequency domain and airspace compressed sensing
US8723723B2 (en) Dual mode ground penetrating radar (GPR)
CN110058193B (en) Digital multi-beam angle measurement method and system based on single receiving channel
CN105158735B (en) Null tone Two-Dimensional Spectral Estimation method based on compression sampling array
CN103596637B (en) The detection method of antenna failure and device
CN102647241A (en) Non-coherent detection system and method for short-wave broad-band channel
CN106597405A (en) Multi-carrier signal form-based ionosphere detection method and system
CN113126096B (en) Satellite-borne hyperspectral multifunctional microwave atmosphere detector
CN103338024B (en) The complementary Kalman filtering apparatus and method of time delay in antenna array
CN116400387B (en) Weak interference detection method and system based on double-array element cross correlation
CN108494508A (en) The high-efficiency frequency spectrum detection method restored based on MWC associated support collection
CN104122273A (en) Radiometer based on multichannel frequency band synthesis
CN103399030B (en) A kind of system and method for realizing answer signal phase difference detection under blank pipe 3/A pattern
CN109936855A (en) Method, apparatus, equipment and the medium being distributed using base station signal retrieving precipitation
CN104950282B (en) Sparse reconstruct is realized in continuous domain broadband signal super-resolution direction-finding method and device
CN110554440A (en) Satellite-borne microwave radiation measuring system and measuring method
CN113671448A (en) Universal broadband radar signal monitoring equipment
CN102130731B (en) High-resolution multi-path channel delay spectrum determination method
CN105099503A (en) Pulse compression method of CHIRP signals, and wireless signal transceiver thereof
CN103312639A (en) Least-square estimation device and least-square estimation method for time delay in antenna array
Christofilakis et al. Rainfall measurements due to radio frequency signal attenuation at 2 GHz
JP4294634B2 (en) Calibration table creation method and arrival direction estimation method
CN104967491B (en) Multichannel width tests system signal reception processing method mutually
CN104090984A (en) Design method for heterogeneous fingerprint database used for indoor positioning
CN103412303A (en) Communication satellite reflected signal remote sensing monitoring system and monitoring method thereof

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