CN111142097A - Target direction estimation method in amplitude modulation broadcast external radiation source radar - Google Patents

Target direction estimation method in amplitude modulation broadcast external radiation source radar Download PDF

Info

Publication number
CN111142097A
CN111142097A CN201910606032.9A CN201910606032A CN111142097A CN 111142097 A CN111142097 A CN 111142097A CN 201910606032 A CN201910606032 A CN 201910606032A CN 111142097 A CN111142097 A CN 111142097A
Authority
CN
China
Prior art keywords
signal
echo signal
echo
virtual array
direct wave
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.)
Pending
Application number
CN201910606032.9A
Other languages
Chinese (zh)
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.)
Wuhan Institute of Technology
Original Assignee
Wuhan Institute of Technology
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 Wuhan Institute of Technology filed Critical Wuhan Institute of Technology
Priority to CN201910606032.9A priority Critical patent/CN111142097A/en
Publication of CN111142097A publication Critical patent/CN111142097A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/418Theoretical aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a target direction estimation method in an amplitude modulation broadcast external radiation source radar, which comprises the following steps: determining a mutual ambiguity function of the direct wave signal and the received signal of each antenna; based on the fuzzy function, obtaining a virtual array signal corresponding to the echo signal according to the time delay and the corresponding Doppler frequency shift of the echo signal to be determined; and acquiring the arrival direction of the echo signal according to the virtual array signal. The method is based on the property of a fuzzy function, and obtains a virtual array signal corresponding to an original echo signal according to the delay of the echo signal to be determined and the corresponding Doppler frequency shift, so that the signal-to-noise ratio and the signal-to-interference ratio of the echo signal are enhanced, the virtual array signal which obviously improves the signal-to-noise ratio and the signal-to-interference ratio of an expected echo signal is constructed, and the direction of the echo signal with the extremely low signal-to-noise ratio is more accurately estimated.

Description

Target direction estimation method in amplitude modulation broadcast external radiation source radar
Technical Field
The invention belongs to the field of wireless signal positioning, and particularly relates to a target direction estimation method in an amplitude modulation broadcast external radiation source radar.
Background
An external radiation source radar system detects and tracks targets by processing signals from non-cooperative sources. It is difficult to measure the Direction-of-Arrival (DOA) of the echo signals because the Signal-to-Noise Ratio (SNR) of the echo signals is usually very low and there is strong direct Signal interference.
Most external radiation source radar systems use time differences and doppler frequency shifts to detect objects, which are difficult to obtain for radiation sources with limited bandwidth, such as amplitude modulated radio signals. Therefore, the orientation of the target becomes an important parameter for detecting, tracking, and positioning the target. Currently, phase interferometer technology, Adcock antenna arrays and range-doppler domain array signal processing methods have been used for some external radiation source radar systems. However, these methods have the defects of large system error, serious coupling effect, lack of theoretical analysis and the like, and are difficult to estimate the direction of the echo signal with extremely low signal-to-noise ratio.
Disclosure of Invention
The invention aims to overcome the defects of the prior art, and provides a target direction estimation method in an amplitude modulation broadcast external radiation source radar, which can solve the problems of enhancing the signal-to-noise ratio and the signal-to-interference ratio of an echo signal, and can construct a virtual array signal for remarkably improving the signal-to-noise ratio and the signal-to-interference ratio of an expected echo signal, so that the direction of the echo signal with extremely low signal-to-noise ratio is more accurately estimated.
The invention provides a target direction estimation method in an amplitude modulation broadcast external radiation source radar, which comprises the following steps: determining a mutual ambiguity function of the direct wave signal and the received signal of each antenna; based on the fuzzy function, obtaining a virtual array signal corresponding to the echo signal according to the time delay and the corresponding Doppler frequency shift of the echo signal to be determined; and acquiring the arrival direction of the echo signal according to the virtual array signal.
The invention has the following beneficial effects: by calculating the cross-ambiguity function of the direct wave and the array, a virtual array signal is constructed according to the delay of the echo signal to be determined and the corresponding Doppler frequency shift, so that the signal-to-noise ratio and the signal-to-interference ratio of the echo signal are enhanced, and the direction of the echo signal with the extremely low signal-to-noise ratio is more accurately estimated.
Drawings
The invention will be further described with reference to the accompanying drawings and embodiments, in which:
FIG. 1 is a flow chart of an embodiment of the present invention;
FIG. 2 is a two-dimensional time-Doppler plot after conventional processing;
FIG. 3 is a graph of azimuth-Doppler shift in an embodiment of the present invention;
fig. 4 is an elevation-doppler shift diagram in an embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1, an embodiment of the present invention provides a method for estimating a target direction in an am broadcast external radiation source radar, including:
and S10, determining the mutual ambiguity function of the direct wave signal and the received signal of each antenna.
In S10, it is assumed that there are 1 direct wave and L echo signals in the received signal, and the received signal x (k) is composed of the direct signal and its echo signals in the form of time delay and doppler shift, which can be expressed as:
Figure BDA0002120796550000021
wherein K is a sampling point, and the total number is K; s (k) is a time domain signal; a isi、ηi、τiV and viRespectively the steering vector, amplitude attenuation, time delay and Doppler frequency shift of the ith echo signal;
Figure BDA0002120796550000022
Fsis the sampling rate; n (k) represents noise; matrices x (k), n (k) and aiIs mx 1, let i-0 denote the direct signal, which corresponds to η0=1,τ 00 and ν0=0。
Let the number of antenna elements be M, and the signal received by each antenna be xm(k) And M is 1,2, …, M, performing a beam forming technology on the array signal, and acquiring a direct wave signal as a reference signal: y isR(k)=wHx (k) s (k), where w is the weight vector of the beamformer. Reference signal yR(k) And the reception signal x of the m-th antennam(k) Is a mutual fuzzy function Fm(τ, ν) is:
Figure BDA0002120796550000023
where τ represents time shift and ν represents frequency shift, indicates conjugation.
And S20, based on the fuzzy function, obtaining a virtual array signal corresponding to the original echo signal according to the delay of the echo signal to be determined and the corresponding Doppler frequency shift.
In S20, FmThe combination of (τ, v) is expressed as a matrix-form array mutual ambiguity function F (τ, v):
F(τ,ν)=[F1(τ,ν),F2(τ,ν),...,FM(τ,ν)]T
where T denotes transpose.
F (τ, ν) is an M × 1 dimensional matrix, which can be expressed as:
Figure BDA0002120796550000031
for an echo signal with i-L, if τL,ν=νLAnd known yR(k) S (k), then:
Figure BDA0002120796550000032
wherein tau isLV and vLTime shift and Doppler shift of the L-th echo signal, aLAnd ηLA steering vector and amplitude attenuation representing the lth echo signal.
From this, a virtual array signal corresponding to the original echo signal, i.e. an enhanced virtual echo signal, is obtained.
And S30, acquiring the arrival direction of the echo signal according to the virtual array signal.
According to the virtual array signal, the direction of arrival of the corresponding echo signal to be determined can be obtained by adopting a direction of arrival estimation method, such as a rotation invariant algorithm, a maximum likelihood estimation algorithm and the like.
According to the target direction estimation method in the amplitude modulation broadcast external radiation source radar, the virtual array signal corresponding to the original echo signal is obtained according to the delay of the echo signal to be determined and the corresponding Doppler frequency shift based on the fuzzy function, so that the signal-to-noise ratio and the signal-to-interference ratio of the echo signal are enhanced, the virtual array signal which obviously improves the signal-to-noise ratio and the signal-to-interference ratio of the expected echo signal is constructed, and the direction of the echo signal with the extremely low signal-to-noise ratio is estimated more accurately.
On the basis of the foregoing embodiment, as an alternative embodiment, before S10, the method further includes: and S00, applying a beam forming technology to the direct wave direction, and acquiring the direct wave signal from the received combined signal of the direct wave signal and the echo signal.
According to the above-mentioned combined signal, the direct wave is obtained by applying the beam forming technique to the direction of the direct wave, i.e. the direct signal y in the reference channelR(k) S (k) and echo signals in the monitoring channels
Figure BDA0002120796550000041
May be obtained by beam forming techniques.
On the basis of the above embodiment, as an alternative embodiment, S30 includes: and according to the virtual array signal, applying a single sampling point MUSIC method to the virtual array signal to obtain the direction of arrival of the echo signal. By the MUSIC method of a single sampling point, more accurate direction of arrival can be obtained.
On the basis of the above embodiment, as an optional embodiment, applying a single sampling point MUSIC method to the virtual array signal, and acquiring the direction of arrival of the echo signal includes: acquiring an array signal of a single sampling point of the virtual array signal, and performing cross-correlation calculation on the array signal of the single sampling point to further obtain a covariance matrix of a recovery rank; and applying an MUSIC method to the covariance matrix of the recovery rank to acquire the direction of arrival of the echo signal.
Specifically, in S30, the array signal of one sampling point is set as:
FL=[F1LL),F2LL),...,FMLL)]T
the following cross-correlation calculation is made:
Figure BDA0002120796550000042
wherein FL(m + i) denotes FLThe m + i th element in (b).
The covariance matrix of the recovery rank is:
Figure BDA0002120796550000043
wherein, KrIs a constant between the number of sources and the number of array elements.
The direction of arrival is estimated using the conventional MUSIC method:
Figure BDA0002120796550000044
wherein ENIs RFThe superscript H denotes the conjugate transpose, a (θ) is the steering vector:
Figure BDA0002120796550000045
where d is the array pitch, λ is the wavelength, and θ is the direction of arrival.
According to the target direction estimation method in the amplitude modulation broadcast external radiation source radar, the direction of arrival of the echo signal is obtained through the single sampling point MUSIC method, and the more accurate direction of arrival can be obtained.
On the basis of the foregoing embodiment, as an alternative embodiment, in S10, the applying a beamforming technique to the direct wave direction includes: a robust adaptive beamforming technique is applied to the direct wave direction. By using a robust adaptive beamforming technique, strong interference can be eliminated, making the finally obtained direction of arrival more accurate.
On the basis of the above embodiment, as an alternative, the dimension of the noise subspace is M-2.
Specifically, F (τ)LL) There are three signal components: a with guide vectorLWith a steered virtual undesired signal vector aiAnd n (k) virtual noise; the ratio DUR of the amplitude of the virtual desired echo signal to the amplitude of the undesired signal is:
Figure BDA0002120796550000051
wherein i is 0,1, …, L-1.
The ratio DNR of the amplitude of the virtual desired echo signal to the virtual noise amplitude is:
Figure BDA0002120796550000052
wherein, PSAnd PNThe power of the direct signal and the noise, respectively.
Since the direct signal is very strong, the direct wave cannot be completely cancelled even if the ratio of the amplitude of the virtual desired echo signal to the amplitude of the virtual undesired signal is very large. In the embodiment of the present invention, the dimension of the noise subspace is set to be M-2, that is, the signal subspace is composed of the desired signal and the residual direct wave.
In a specific application of the invention, a Uniform Circular Array (UCA) of 16 antennas is provided, with the radiation source being more than 1000 kilometers away, so that the transmitted signal reaches the receiving array by ionospheric reflection. The frequency of the source signal is 15.5MHz and the bandwidth <8 kHz.
Setting the receiving array signal as x (k), the processing flow is as follows:
step 1: forming a beam in the direct wave direction to obtain a reference signal: y isR(k)=wHx(k)≈s(k)。
Step 2: calculating a cross-ambiguity function of the reference signal and the array signal:
Figure BDA0002120796550000061
and 3, step 3: since the bandwidth of the AM broadcast signal is very narrow, a time difference τ is arbitrarily selectedLAt a Doppler shift of vLForming a single sample point array signal FL=[F1LL),F2LL),...,FMLL)]T
And 4, step 4: calculating a covariance matrix of the recovery rank:
Figure BDA0002120796550000062
wherein
Figure BDA0002120796550000063
And 5, step 5: to RFCarrying out eigenvalue decomposition to obtain a signal subspace EN
And 6, step 6: computing spatial spectra
Figure BDA0002120796550000064
Actual time of calculation vLCan be obtained by scanning, and can be combined to obtain a Doppler-DOA image.
The results of the conventional processing of the test data and the results of the inventive processing are shown in FIGS. 2-4. Figure 2 is a time-doppler plot of the combination of data extracted in the doppler dimension from a 2 minute conventional range-doppler plot, from which figure 2 a number of clear target tracks can be seen. And selecting the target at 85 seconds and 18.3Hz for target direction estimation. Since the uniform circular array has two directional dimensions, fixing one of them to scan the other, for example fixing a pitch angle of 65 °, can obtain the degree of the azimuth angle corresponding to the pitch angle. Fig. 3 and 4 are the results of the azimuth and elevation angle estimates of the target at 85 seconds, respectively, where fig. 3 is a doppler-azimuth plot for a fixed elevation angle of 65 ° and fig. 4 is a doppler-elevation plot for a fixed elevation angle of 213 °.

Claims (6)

1. A target direction estimation method in an amplitude modulation broadcast external radiation source radar is characterized by comprising the following steps:
s10, determining a mutual ambiguity function of the direct wave signal and the received signal of each antenna;
s20, based on the fuzzy function, obtaining a virtual array signal corresponding to the echo signal according to the time delay and the corresponding Doppler frequency shift of the echo signal to be determined;
and S30, acquiring the arrival direction of the echo signal according to the virtual array signal.
2. The method of claim 1, wherein before S10, the method further comprises:
and S00, applying a beam forming technology to the direct wave direction, and acquiring the direct wave signal from the received combined signal of the direct wave signal and the echo signal.
3. The method of claim 1, wherein the step S30 comprises:
and according to the virtual array signal, applying an MUSIC method of a single sampling point to the virtual array signal to obtain the direction of arrival of the echo signal.
4. The method of claim 3, wherein the applying the single-sample-point MUSIC method to the virtual array signal to obtain the direction of arrival of the echo signal comprises:
acquiring an array signal of a single sampling point of the virtual array signal, and performing cross-correlation calculation on the array signal of the single sampling point to obtain a cross-correlation function;
acquiring a covariance matrix of a recovery rank according to a cross-correlation function;
and acquiring the arrival direction of the echo signal based on the MUSIC method according to the covariance matrix and the corresponding noise subspace.
5. The method of claim 2, wherein the step of applying a beamforming technique to the direct wave direction in S10 comprises:
a robust adaptive beamforming technique is applied to the direct wave direction.
6. The method of claim 4, wherein the dimension of the noise subspace is M-2;
wherein M is the number of antenna elements.
CN201910606032.9A 2019-07-05 2019-07-05 Target direction estimation method in amplitude modulation broadcast external radiation source radar Pending CN111142097A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910606032.9A CN111142097A (en) 2019-07-05 2019-07-05 Target direction estimation method in amplitude modulation broadcast external radiation source radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910606032.9A CN111142097A (en) 2019-07-05 2019-07-05 Target direction estimation method in amplitude modulation broadcast external radiation source radar

Publications (1)

Publication Number Publication Date
CN111142097A true CN111142097A (en) 2020-05-12

Family

ID=70516732

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910606032.9A Pending CN111142097A (en) 2019-07-05 2019-07-05 Target direction estimation method in amplitude modulation broadcast external radiation source radar

Country Status (1)

Country Link
CN (1) CN111142097A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115508780A (en) * 2022-11-23 2022-12-23 杭州兆华电子股份有限公司 Synthetic aperture acoustic imaging method

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1253434A1 (en) * 2001-04-27 2002-10-30 Mitsubishi Electric Information Technology Centre Europe B.V. Method for estimating a direction of arrival
CN103353595A (en) * 2013-06-18 2013-10-16 西安电子科技大学 Meter wave radar height measurement method based on array interpolation compression perception
CN104515969A (en) * 2013-09-29 2015-04-15 电子科技大学 Hexagonal array-based coherent signal two-dimensional DOA (Direction of Arrival) estimation method
CN104698433A (en) * 2015-03-16 2015-06-10 电子科技大学 Single-snapshot data-based coherent signal DOA (direction of arrival) estimating method
CN106842115A (en) * 2017-01-25 2017-06-13 天津大学 The two-dimentional direction-finding method of principle is damaged based on ROOT MUSIC algorithms and order
CN108680910A (en) * 2018-05-15 2018-10-19 河海大学 Frequency modulation broadcasting external illuminators-based radar object detection method based on waveform cognition
CN108872946A (en) * 2018-04-20 2018-11-23 西安电子科技大学 The robust ada- ptive beamformer method of steering vector and covariance matrix Joint iteration
CN108957388A (en) * 2018-05-21 2018-12-07 南京信息工程大学 A kind of MIMO radar coherent DOA estimation method based on covariance matching SL0 algorithm

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1253434A1 (en) * 2001-04-27 2002-10-30 Mitsubishi Electric Information Technology Centre Europe B.V. Method for estimating a direction of arrival
CN103353595A (en) * 2013-06-18 2013-10-16 西安电子科技大学 Meter wave radar height measurement method based on array interpolation compression perception
CN104515969A (en) * 2013-09-29 2015-04-15 电子科技大学 Hexagonal array-based coherent signal two-dimensional DOA (Direction of Arrival) estimation method
CN104698433A (en) * 2015-03-16 2015-06-10 电子科技大学 Single-snapshot data-based coherent signal DOA (direction of arrival) estimating method
CN106842115A (en) * 2017-01-25 2017-06-13 天津大学 The two-dimentional direction-finding method of principle is damaged based on ROOT MUSIC algorithms and order
CN108872946A (en) * 2018-04-20 2018-11-23 西安电子科技大学 The robust ada- ptive beamformer method of steering vector and covariance matrix Joint iteration
CN108680910A (en) * 2018-05-15 2018-10-19 河海大学 Frequency modulation broadcasting external illuminators-based radar object detection method based on waveform cognition
CN108957388A (en) * 2018-05-21 2018-12-07 南京信息工程大学 A kind of MIMO radar coherent DOA estimation method based on covariance matching SL0 algorithm

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
REN ET AL: "EXTENDING MUSIC TO SINGLE SNAPSHOT AND ON LINE DIRECTION FINDING APPLICATIONS", 《RADAR》 *
YANG LI ET AL: "YANG LI,DOA estimation for echo signals and experimental results in the AM radio-based passive radar", 《IEEE ACCESS》 *
YANG LI: "DOA estimation for echo signals and experimental results in the AM radio-based passive radar" *
吴海洲 等: "数字电视辐射源雷达基于空域滤波的直达波获取", 兵工学报 *
李磊 等: "通道幅相误差条件下的相干信源波达方向估计", 系统工程与电子技术 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115508780A (en) * 2022-11-23 2022-12-23 杭州兆华电子股份有限公司 Synthetic aperture acoustic imaging method

Similar Documents

Publication Publication Date Title
Vasanelli et al. Calibration and direction-of-arrival estimation of millimeter-wave radars: A practical introduction
CN110488255B (en) Monopulse high-resolution angle measurement system and method for phased array radar
CN108693511B (en) Moving target angle calculation method of time division multiplexing MIMO radar
CN106226754A (en) The low elevation angle based on time reversal Wave arrival direction estimating method
CN107390197B (en) Radar self-adaption sum-difference beam angle measurement method based on feature space
CN110426670B (en) Super-resolution DOA estimation method for external radiation source radar based on TLS-CS
CN103926555B (en) A kind of method that utilization not rounded signal measuring antenna array receiver machine width is mutually responded
Bekkerman et al. Spatially coded signal model for active arrays
CN108828504B (en) MIMO radar target direction fast estimation method based on partial correlation waveform
CN114002673A (en) Satellite-borne passive SAR non-cooperative signal sensing system and multi-dimensional parameter estimation method
CN105158754B (en) A kind of method that target positioning is carried out using multiple input single output radio system
CN111693979B (en) Digital array monopulse angle measurement method based on Taylor expansion
US11269052B2 (en) Signal processing method
CN111142097A (en) Target direction estimation method in amplitude modulation broadcast external radiation source radar
Sakhnini et al. An experimental evaluation of robust near-field radar localization using a massive MIMO testbed
CN114488142A (en) Radar two-dimensional angle imaging method and system based on difference-sum beam
CN114325560A (en) Super-resolution target direction finding method for beam scanning radar
CN112444790A (en) Method for detecting target under strong interference condition
CN113126076A (en) Collaborative distance and angle parameter estimation method and device based on dual-mode radar
Huang et al. In-field calibration of passive array receiver using detected target
Huang et al. Using time reversal for direction of arrival estimation in multipath environment
Berry et al. Generalised phase monopulse for multi-target DoA estimation and extended target spatial imaging
CN109752688B (en) Method for calculating angle difference of adjacent information sources for sensor array system
Tang et al. Bistatic MIMO radar height estimation method based on adaptive beam-space RML data fusion
CN112986899B (en) DOA estimation method of airborne MIMO radar in multipath environment

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200512