CN107037413A - The direct wave suppressing method offseted based on adaptive Doppler domain wave beam - Google Patents
The direct wave suppressing method offseted based on adaptive Doppler domain wave beam Download PDFInfo
- Publication number
- CN107037413A CN107037413A CN201710098338.9A CN201710098338A CN107037413A CN 107037413 A CN107037413 A CN 107037413A CN 201710098338 A CN201710098338 A CN 201710098338A CN 107037413 A CN107037413 A CN 107037413A
- Authority
- CN
- China
- Prior art keywords
- direct
- doppler
- domain
- signal
- array
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4052—Means for monitoring or calibrating by simulation of echoes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4052—Means for monitoring or calibrating by simulation of echoes
- G01S7/4082—Means for monitoring or calibrating by simulation of echoes using externally generated reference signals, e.g. via remote reflector or transponder
- G01S7/4095—Means for monitoring or calibrating by simulation of echoes using externally generated reference signals, e.g. via remote reflector or transponder the external reference signals being modulated, e.g. rotating a dihedral reflector or modulating a transponder for simulation of a Doppler echo
Abstract
The present invention relates to a kind of direct wave suppressing method offseted based on adaptive Doppler domain wave beam, comprising：Using the constant modulus property of direct-path signal, direct-path signal is extracted by the permanent mould method of space-time, and the two-dimensional array data that array time domain data is converted to range-Dopler domain are matched by being concerned with；Two-dimensional array data are subjected to piecemeal, and wave beam is carried out by adaptive diagonal loading Beamforming Method in Doppler domain offseting, suppress direct-path signal and extract target weak echo signal.The present invention carries out adaptive beam synthesis in " range Doppler " 2-D data domain, and dimension of adjusting the distance carries out piecemeal, and avoid the region that Doppler is zero, greatly improve when carrying out direct wave suppression based on conventional beam synthesizing method target weak echo signal cancellation and the algorithm free degree is mostly used in the problem of direct wave suppresses, lift target weak echo signal detection performance；Direct-path signal is extracted in the presence of multipath, interference, performance is stable, reliable, and efficiently.
Description
Technical field
The invention belongs to external sort algorithm direct wave suppression technology field, there is constant modulus property especially for direct-path signal
When there is provided a kind of direct wave suppressing method offseted based on adaptive Doppler domain wave beam.
Background technology
External illuminators-based radar is by the use of noncooperative radiation source as transmiting signal source to target acquisition, tracking.With biography
System active target Detection Techniques are compared, the energy with good disguised and very strong survival ability, good active jamming rejection
The advantages such as power, the ability for resisting stealth target.Therefore, the technology causes the extensive concern of countries in the world, as various countries' national defence neck
Domain to it is empty, over the ground, the powerful of opposite bank early warning.But need to solve in actual applications there are still some key issues,
The main waveform analysis that may be summarized to be external irradiation source signal and optimum choice, the purification of direct-path signal, direct wave and many
The suppression of the clutters such as footpath, the detection of target and identification and the positioning of target and tracking etc., wherein, direct wave purification is clutter recognition
Reference signal is provided with target detection, and the suppression of the clutter such as direct wave and multipath is the follow-up elder generation for carrying out target positioning and tracking
Certainly condition.
Being presently used for the algorithm of direct wave purification mainly has the balancing technique of time domain, on cepstrum and spatial domain on frequency domain
Adaptive beamformer.Although time domain equalization algorithm can track the characteristic of channel, the estimation to direct-path signal is completed,
Required exponent number is higher in the case of ensureing certain error；Cepstrum Technique can reach well when handling radar pulse signal
Refining effect, but for some signals, the resolution ratio in cepstrum domain is relatively low, therefore be no longer applicable；Adaptive beamformer
Method can effectively clutter reduction and interference, but be due to that the influence of signal correlation makes performance drastically decline.2008
F.Colone et al., which proposes space-time joint constant modulus algorithm, is used to removing the clutters such as multipath in reference channel, the algorithm it is through
Ripple purification performance is in most cases substantially better than simple spatial domain constant modulus algorithm and time domain constant modulus algorithm.On direct wave with
The suppression of multipath clutter, the method for mainly having spatial domain and time domain at present.Wherein, the method in spatial domain is mainly adaptive beam shape
Into this kind of method forms null to reach the purpose of clutter reduction at direct wave multipath clutter.But, under normal circumstances it
Need to know the prior information of target echo, weak mesh is easily offset such as coming to for target, and when suppressing strong direct wave
Mark echo-signal.The method of time domain mainly has two classes, and a class is the method for adaptive interference cancellation, such as BLMS algorithms, bulk
NLMS algorithms, EDS algorithms etc., this kind of method need to set the exponent number of wave filter, therefore it is whole to be more than filter order for time delay
The clutter recognition performance in several times sampling period drastically declines.Another kind of is extension phase elimination (Extensive Cancellation
Algorithm, ECA).This kind of algorithm clutter recognition better performances, but reference signal hydraulic performance decline and algorithm meter when including multipath
Calculation amount is big.Australian scholar G.Fabrizio proposes that one kind is based on adaptive Doppler domain Beamforming Method, and this method is led to
Coherent computing lifting signal to noise ratio is crossed, and selects the sample data for calculating weights effectively to realize the suppression of direct wave and weak mesh
The extraction of echo-signal is marked, but this method extracts direct-path signal by conventional beamformer method, can cause target in distance
That ties up is fuzzy.
The content of the invention
The deficiency with suppression technology is extracted for existing direct-path signal, the present invention provides a kind of based on adaptively how general
The direct wave suppressing method that domain wave beam is offseted is strangled, effectively realizes that direct wave suppresses and weak signal target echo-signal is extracted, improves target
Detection probability.
It is a kind of to be suppressed based on the direct wave that adaptive Doppler domain wave beam is offseted according to design provided by the present invention
Method, is comprised the following steps：
Step 1, the constant modulus property using direct-path signal, extract direct-path signal, and pass through by the permanent mould method of space-time
Relevant matching is converted to array time domain data the two-dimensional array data in distance-Doppler domain；
Step 2, two-dimensional array data are carried out to piecemeal, and pass through adaptive diagonal loading Wave beam forming side in Doppler domain
Method carries out wave beam and offseted, and suppresses direct-path signal and extracts target weak echo signal.
Above-mentioned, described step 1 includes following content：
Step 11, using L * channel array antenna and according to the radio of Nyquist sampling thheorem synchronous acquisition target emanations
Signal, obtains the array time domain data of signal；
Step 12, the DOA information for estimating according to the array time domain data of signal direct wave, pass through what is constrained based on direction
Space-time constant modulus algorithm extracts direct-path signal；
Step 13, using the direct-path signal extracted as reference signal, enter with every passage time domain data of aerial array
The relevant matching operation of row, array time domain data is converted to the two-dimensional array data in distance-Doppler domain.
It is preferred that, in described step 11, signal array model is expressed as：X (t)=AS (t)+N (t), wherein, N (t) is
The noise data vector that L × 1 is tieed up,
, A is L × (Nm+ 1) array manifold matrix is tieed up, L is element number of array, NmFor in reference channel multipath clutter disturb it is total
Number, θ0For the arrival bearing of direct-path signal, θiFor the arrival bearing of i-th multipath clutter；S (t) is (Nm+ 1) letter × 1 tieed up
Number vector；sd(t)=d (t) andRespectively direct-path signal and multipath clutter signal, τm_iFor
The time delay of noise jamming, fdmFor Doppler shift.
It is preferred that, described step 12 includes following content：The DOA information of direct wave is estimated, using permanent mould CMA algorithms,
Constrained by weighting matrix and introducing direct-path signal direction, extract direct-path signal.
It is preferred that, constrained, be embodied as by weighting matrix and introducing direct-path signal direction：
, wherein, E represents to expect, y (n) is the output at n-th of moment, and vec is the vectorization of matrix, and has：
wijPostpone the weight coefficient at moment at j-th for i-th of antenna, the weighting matrix being made from it is W=<wi,j
>L×P；X (n) is the receiving matrix at n-th of moment of reference edge, and P is the number of delay tap coefficient；Assuming that A (θ0)=a (θ0)·
C, c=[0,0 ..., R2,0,...,0]1×P,And R2Positioned at (P-1)/2 element, θ0For direct wave
Direction of arrival, vec (W)Hvec(A(θ0))=1 it is that the direction of direct wave is constrained, the more new formula of weighing vector is：
, wherein, XN=[vec (X (1)), vec (X (2)) ..., vec (X (N))], w=vec (W),
Pass throughObtain：
It is preferred that, array time domain data is converted to the two-dimensional array data in distance-Doppler domain in described step 13,
Conversion formula is：
, wherein, f (t) is window function, and x (t) is array channel output data, and y (t) is the direct-path signal extracted, M
For number of samples, l is range cell, and k represents doppler cells；Array format is expressed as
, wherein, c (l, k) is the two-dimensional array data in " distance-Doppler domain ".
It is preferred that, described step 2 is as follows comprising content：
Step 21, the two-dimensional array data progress Q deciles by every passage, determine that two do not include the protections that Doppler is zero
Region, and calculate the array covariance matrix of corresponding region after two-dimensional array data are divided；
Step 22, based on the Beamforming Method adaptively diagonally loaded, the not same period is solved by array covariance matrix
Hope the corresponding weight w of steering vector；
Step 23, according to direct wave the characteristics of Doppler domain is symmetrical, Doppler domain carry out spatial domain offset, gone directly
Range-Doppler spectrum after ripple suppression.
It is preferred that, described step 21 includes following content：Distance-Doppler two-dimensional array data to each passage are entered
The following operation of row：Range cell M is averagely divided into NbPart, every part includes Q range cell, i.e. Nb=M/Q；In each distance
On unit, the 0Hz protections region of Doppler domain is set, and is designated as grAnd gl, then the interval of positive doppler cells is k ∈ G+≡{gr,gr+ 1 ..., K }, it is k ∈ G to bear Doppler's interval-≡{-gl,-gl-1,...,-K}；Obtain each range cell
The corresponding covariance matrix of positive and negative Doppler be：
, wherein, R+(b)、R-(b) it is positive Doppler and the corresponding covariance matrix of negative Doppler respectively, b ∈ 1,2 ...,
Nb。
It is above-mentioned, the corresponding weight w of different expectation steering vectors is solved in described step 22, following content is included：Root
According to covariance matrix, the corresponding weights of each range cell, desired orientation are asked for according to calculation formula, calculation formula is：
Wherein, s (θs,φs) be desired signal steering vector, θs,φsFor the azimuth and the elevation angle of desired signal.
Above-mentioned, described step 23 includes following content：According to direct-path signal Doppler centered on 0Hz it is symmetrical special
Property, the positive weights in step 21 are used to bear the Beam synthesis of Doppler's correspondence array data, negative weights for positive Doppler couple
The Beam synthesis of array data is answered, the direct wave sidelobe cancellation of Doppler domain is realized and retains target echo signal, it is specific to represent
It is as follows：
Beneficial effects of the present invention：
The present invention carries out adaptive beam synthesis in " distance-Doppler " 2-D data domain, and dimension of adjusting the distance carries out piecemeal,
And the region that Doppler is zero is avoided, greatly improve weak time of target during based on conventional beam synthesizing method progress direct wave suppression
Ripple signal cancellation and the algorithm free degree are mostly used in the problem of direct wave suppresses, and lift target weak echo signal detection performance；
In addition, utilizing the space-time constant modulus algorithm constrained based on direction in the present invention, direct-path signal is extracted in the presence of multipath, interference,
Performance is stable, reliable, and efficiently, with stronger actual application value.
Brief description of the drawings：
Fig. 1 is schematic flow sheet of the invention；
Fig. 2 changes the two-dimensional array data flow diagram in distance-Doppler domain for array time domain data in the present invention；
Fig. 3 is carries out the schematic flow sheet that wave beam is offseted in the present invention by adaptive diagonal loading Beamforming Method；
Fig. 4 is theory diagram of the invention；
Fig. 5 is experimental principle schematic diagram；
Fig. 6 is that direct wave purifies result schematic diagram；
Fig. 7 is that direct wave suppresses comparing result schematic diagram.
Embodiment：
The present invention is further detailed explanation with technical scheme below in conjunction with the accompanying drawings, and detailed by preferred embodiment
Describe bright embodiments of the present invention in detail, but embodiments of the present invention are not limited to this.
Embodiment one, a kind of shown in Figure 1, direct wave suppressing method offseted based on adaptive Doppler domain wave beam,
Comprise the following steps：
Step 1, the constant modulus property using direct-path signal, extract direct-path signal, and pass through by the permanent mould method of space-time
Relevant matching is converted to array time domain data the two-dimensional array data in distance-Doppler domain；
Step 2, two-dimensional array data are carried out to piecemeal, and pass through adaptive diagonal loading Wave beam forming side in Doppler domain
Method carries out wave beam and offseted, and suppresses direct-path signal and extracts target weak echo signal.
The present invention carries out adaptive beam synthesis in " distance-Doppler " 2-D data domain, and dimension of adjusting the distance carries out piecemeal,
And the region that Doppler is zero is avoided, greatly improve weak time of target during based on conventional beam synthesizing method progress direct wave suppression
Ripple signal cancellation and the algorithm free degree are mostly used in the problem of direct wave suppresses, and lift target weak echo signal detection performance,
Reliability is high, with stronger actual application value.
Embodiment two, referring to shown in Fig. 1~3, a kind of direct wave suppression side offseted based on adaptive Doppler domain wave beam
Method, includes following content：
One) using the constant modulus property of direct-path signal, direct-path signal is extracted by the permanent mould method of space-time, and by relevant
Matching is converted to array time domain data the two-dimensional array data in distance-Doppler domain, and particular content is as follows：
1st, using L * channel array antenna and according to the radio signal of Nyquist sampling thheorem synchronous acquisition target emanations,
The array time domain data of signal is obtained, signal array model is expressed as：X (t)=AS (t)+N (t),
Wherein, N (t) is the noise data vector that L × 1 is tieed up,
, A is L × (Nm+ 1) array manifold matrix is tieed up, L is element number of array, NmFor in reference channel multipath clutter disturb it is total
Number, θ0For the arrival bearing of direct-path signal, θiFor the arrival bearing of i-th multipath clutter；S (t) is (Nm+ 1) letter × 1 tieed up
Number vector；sd(t)=d (t) andRespectively direct-path signal and multipath clutter signal, τm_iFor
The time delay of noise jamming, fdmFor Doppler shift.
2nd, the DOA information of direct wave is estimated according to the array time domain data of signal, passes through the space-time constrained based on direction
Constant modulus algorithm extracts direct-path signal
By estimating the DOA information of direct wave, using permanent mould CMA algorithms, pass through weighting matrix and introducing direct-path signal
Direction is constrained, and is embodied as：
, wherein, E represents to expect, y (n) is the output at n-th of moment, and vec is the vectorization of matrix, and has：
wijPostpone the weight coefficient at moment at j-th for i-th of antenna, the weighting matrix being made from it is W=<wi,j
>L×P；X (n) is the receiving matrix at n-th of moment of reference edge, and P is the number of delay tap coefficient；Assuming that A (θ0)=a (θ0)·
C, c=[0,0 ..., R2,0,...,0]1×P,And R2Positioned at (P-1)/2 element, θ0For direct wave
Direction of arrival, vec (W)Hvec(A(θ0))=1 it is that the direction of direct wave is constrained, the more new formula of weighing vector is：
, wherein, XN=[vec (X (1)), vec (X (2)) ..., vec (X (N))], w=vec (W),
Pass throughObtain：
3rd, the direct-path signal extracted is concerned with as reference signal with every passage time domain data of aerial array
Matching operation, array time domain data is converted to the two-dimensional array data in distance-Doppler domain, and conversion formula is：
, wherein, f (t) is window function, and x (t) is array channel output data, and y (t) is the direct-path signal extracted, M
For number of samples, l is range cell, and k represents doppler cells；Array format is expressed as
, wherein, c (l, k) is the two-dimensional array data in " distance-Doppler domain ".
Two) two-dimensional array data are subjected to piecemeal, and pass through adaptive diagonal loading Beamforming Method in Doppler domain
Carry out wave beam to offset, suppress direct-path signal and extract target weak echo signal, particular content is as follows：
1st, the two-dimensional array data of every passage are subjected to Q deciles, determine that two do not include the protection zone that Doppler is zero
Domain, and calculate the array covariance matrix of corresponding region after two-dimensional array data are divided.
Distance-Doppler two-dimensional array data to each passage are proceeded as follows：To prevent target echo signal
Offset, range cell M is averagely divided into NbPart, every part includes Q range cell, i.e. Nb=M/Q；To prevent excessive freedom
The suppression for direct wave is spent, in each range cell, the 0Hz protections region of Doppler domain is set, and be designated as grAnd gl,
Then the interval of positive doppler cells is k ∈ G+≡{gr,gr+ 1 ..., K }, it is k ∈ G to bear Doppler's interval-≡{-
gl,-gl-1,...,-K}；The corresponding covariance matrix of positive and negative Doppler for obtaining each range cell is：
, wherein, R+(b)、R-(b) it is positive Doppler and the corresponding covariance matrix of negative Doppler respectively, b ∈ 1,2 ...,
Nb。
2nd, based on the Beamforming Method adaptively diagonally loaded, different expectations is solved by array covariance matrix and are oriented to
The corresponding weight w of vector.
According to covariance matrix, the corresponding weights of each range cell, desired orientation are asked for according to calculation formula, calculate public
Formula is：
Wherein, s (θs,φs) be desired signal steering vector, θs,φsFor the azimuth and the elevation angle of desired signal.
3rd, according to direct wave the characteristics of Doppler domain is symmetrical, carry out spatial domain in Doppler domain and offset, obtain direct wave suppression
Range-Doppler spectrum after system.
According to direct-path signal Doppler centered on 0Hz symmetry characteristic, the positive weights in step 21 are used for negative how general
Strangling the Beam synthesis of correspondence array data, negative weights is used for the Beam synthesis of positive Doppler correspondence array data, realizes Doppler
The direct wave sidelobe cancellation and reservation target echo signal in domain, are specifically expressed as follows：
The direct wave that the present invention is directed in external illuminators-based radar technology suppresses problem and is based on adaptive Doppler there is provided one kind
The direct wave suppressing method that domain wave beam is offseted, using " distance-Doppler " data of array antenna realize strong direct wave suppress and
Target weak echo signal extraction.First with " constant modulus property " of direct-path signal, pass through the permanent mould Methods For Purification direct wave of space-time
Signal, is converted to array time domain data by the matching treatment that is concerned with based on this 2-D data in " distance-Doppler " domain, then
Two-dimensional array data are subjected to piecemeal, and in Doppler domain, wave beam pair is carried out using adaptive diagonal loading Beamforming Method
Disappear, to suppress direct-path signal and extract target weak echo signal.Suppress with direct wave of the tradition based on Adaptive beamformer
Algorithm is compared, and it is not enough by " counteracting " and the Beam synthesis algorithm free degree that the method that the present invention is provided solves target weak echo signal
The problems such as, significantly improve direct wave inhibition.
It is shown in Figure 5 to verify effectiveness of the invention, it is the measured test that direct wave suppresses, experiment uses sighting distance
Outer over-the-horizon radar signal is detected and positioned to airbound target in sighting distance as non-cooperation radiation source.Experiment uses straight
Footpath is 100 meters of 20 passage uniform circular arrays as receiving array, and utilizes Flightradar24 software monitors receiving array sighting distances
In the range of aircarrier aircraft be used as cooperative target；It is 225 degree, faces upward by measures the azimuth for obtaining non-cooperation radiation source with analysis
Angle is 16 degree, and centre frequency is 19.493MHz, and signal modulation mode is linear frequency modulation continuous wave (LFMCW).Receiving array is simultaneously
Direct wave and target reflection echo are received, the sample frequency of baseband signal is 62.5kHz.Below by method disclosed by the invention with
Direct wave restrainable algorithms based on conventional beamformer algorithm are compared, and direct-path signal purification is using based on spatial domain constraint
Space-time constant modulus algorithm, direct wave suppresses to offset algorithm using adaptive Doppler domain wave beam.In Fig. 6, Fig. 6-1 gives original straight
Up to ripple schematic diagram：Figure a is original direct-path signal schematic diagram, and figure b is the corresponding sonagrams of figure a；6-2 is given after direct wave purification
Schematic diagram：It is the signal time domain schematic diagram after direct wave is purified to scheme a, and figure b is the sonagram corresponding to figure a；Pass through Fig. 6-1 and 6-
Shown in 2, direct-path signal can preferably be extracted by demonstrating the present invention.In Fig. 7, b represents the doppler spectral of target position
Scheme, a and c are represented by method disclosed by the invention and the corresponding direct wave histamine result of conventional beamformer method respectively, from
Contrast is this it appears that algorithm proposed by the present invention has preferable direct wave rejection
The foregoing description of the disclosed embodiments, enables professional and technical personnel in the field to realize or use the application.
A variety of modifications to these embodiments will be apparent for those skilled in the art, as defined herein
General Principle can in other embodiments be realized in the case where not departing from spirit herein or scope.Therefore, the application
The embodiments shown herein is not intended to be limited to, and is to fit to and principles disclosed herein and features of novelty phase one
The most wide scope caused.
Claims (10)
1. a kind of direct wave suppressing method offseted based on adaptive Doppler domain wave beam, it is characterised in that include following content：
Step 1, the constant modulus property using direct-path signal, direct-path signal is extracted by the permanent mould method of space-time, and by relevant
Matching is converted to array time domain data the two-dimensional array data in distance-Doppler domain；
Step 2, two-dimensional array data are carried out to piecemeal, and entered in Doppler domain by adaptive diagonal loading Beamforming Method
Traveling wave beam is offseted, and is suppressed direct-path signal and is extracted target weak echo signal.
2. the direct wave suppressing method according to claim 1 offseted based on adaptive Doppler domain wave beam, its feature is existed
In described step 1 includes following content：
Step 11, using L * channel array antenna and according to the aerogram of Nyquist sampling thheorem synchronous acquisition target emanations
Number, obtain the array time domain data of signal；
Step 12, the DOA information for estimating according to the array time domain data of signal direct wave, pass through the space-time constrained based on direction
Constant modulus algorithm extracts direct-path signal；
Step 13, using the direct-path signal extracted as reference signal, carry out phase with every passage time domain data of aerial array
Dry matching operation, array time domain data is converted to the two-dimensional array data in distance-Doppler domain.
3. the direct wave suppressing method according to claim 2 offseted based on adaptive Doppler domain wave beam, its feature is existed
In in described step 11, signal array model is expressed as：X (t)=AS (t)+N (t),
Wherein, N (t) is the noise data vector that L × 1 is tieed up,
,
A is L × (Nm+ 1) array manifold matrix is tieed up, L is element number of array, NmThe sum disturbed for multipath clutter in reference channel, θ0
For the arrival bearing of direct-path signal, θiFor the arrival bearing of i-th multipath clutter；S (t) is (Nm+ 1) × 1 the signal of dimension is sweared
Amount；sd(t)=d (t) andRespectively direct-path signal and multipath clutter signal, τm_iFor clutter
The time delay of interference, fdmFor Doppler shift.
4. the direct wave suppressing method according to claim 2 offseted based on adaptive Doppler domain wave beam, its feature is existed
In described step 12 includes following content：The DOA information of direct wave is estimated, using permanent mould CMA algorithms, passes through weighting matrix
And the constraint of direct-path signal direction is introduced, extract direct-path signal.
5. the direct wave suppressing method according to claim 4 offseted based on adaptive Doppler domain wave beam, its feature is existed
In by weighting matrix and introducing direct-path signal direction and constrain, be embodied as：
,
Wherein, E represents to expect, y (n) is the output at n-th of moment, and vec is the vectorization of matrix, and has：
1
wijPostpone the weight coefficient at moment at j-th for i-th of antenna, the weighting matrix being made from it is W=<wi,j>L×P；X
(n) it is the receiving matrix at n-th of moment of reference edge, P is the number of delay tap coefficient；Assuming that A (θ0)=a (θ0) c, c=
[0,0,...,R2,0,...,0]1×P,And R2Positioned at (P-1)/2 element, θ0Reached for the ripple of direct wave
Direction, vec (W)Hvec(A(θ0))=1 it is that the direction of direct wave is constrained, the more new formula of weighing vector is：
,
Wherein, XN=[vec (X (1)), vec (X (2)) ..., vec (X (N))], w=vec (W),Pass throughObtain：
6. the direct wave suppressing method according to claim 2 offseted based on adaptive Doppler domain wave beam, its feature is existed
In array time domain data being converted to the two-dimensional array data in distance-Doppler domain, conversion formula is in described step 13：
,
Wherein, f (t) is window function, and x (t) is array channel output data, and y (t) is the direct-path signal extracted, and M is sampling point
Number, l is range cell, and k represents doppler cells；Array format is expressed as
,
Wherein, c (l, k) is the two-dimensional array data in " distance-Doppler domain ".
7. the direct wave suppressing method according to claim 2 offseted based on adaptive Doppler domain wave beam, its feature is existed
In described step 2 is as follows comprising content：
Step 21, the two-dimensional array data progress Q deciles by every passage, determine that two do not include the protection zones that Doppler is zero
Domain, and calculate the array covariance matrix of corresponding region after two-dimensional array data are divided；
Step 22, based on the Beamforming Method adaptively diagonally loaded, different expectations are solved by array covariance matrix and led
To the corresponding weight w of vector；
Step 23, according to direct wave the characteristics of Doppler domain is symmetrical, Doppler domain carry out spatial domain offset, obtain direct wave suppression
Range-Doppler spectrum after system.
8. the direct wave suppressing method according to claim 7 offseted based on adaptive Doppler domain wave beam, its feature is existed
In described step 21 includes following content：Distance-Doppler two-dimensional array data to each passage are proceeded as follows：
Range cell M is averagely divided into NbPart, every part includes Q range cell, i.e. Nb=M/Q；In each range cell, set
The 0Hz protections region of Doppler domain, and it is designated as grAnd gl, then the interval of positive doppler cells is k ∈ G+≡{gr,gr+
1 ..., K }, it is k ∈ G to bear Doppler's interval-≡{-gl,-gl-1,...,-K}；Obtain the positive and negative many of each range cell
It is general to strangle corresponding covariance matrix and be：
,
Wherein, R+(b)、R-(b) it is respectively positive Doppler and the corresponding covariance matrix of negative Doppler, b ∈ 1,2 ..., Nb。
9. the direct wave suppressing method according to claim 8 offseted based on adaptive Doppler domain wave beam, its feature is existed
In solving the corresponding weight w of different expectation steering vectors in described step 22, include following content：According to covariance matrix,
The corresponding weights of each range cell, desired orientation are asked for according to calculation formula, calculation formula is：
Wherein, s (θs,φs) be desired signal steering vector, θs,φsFor the azimuth and the elevation angle of desired signal.
10. the direct wave suppressing method according to claim 9 offseted based on adaptive Doppler domain wave beam, its feature is existed
In described step 23 includes following content：According to direct-path signal Doppler centered on 0Hz symmetry characteristic, by step 21
In positive weights be used to bear Beam synthesis, the negative weights of Doppler's correspondence array data and be used for positive Doppler correspondence array data
Beam synthesis, realizes the direct wave sidelobe cancellation of Doppler domain and retains target echo signal, be specifically expressed as follows：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710098338.9A CN107037413B (en) | 2017-02-22 | 2017-02-22 | The direct wave suppressing method offseted based on adaptive Doppler domain wave beam |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710098338.9A CN107037413B (en) | 2017-02-22 | 2017-02-22 | The direct wave suppressing method offseted based on adaptive Doppler domain wave beam |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107037413A true CN107037413A (en) | 2017-08-11 |
| CN107037413B CN107037413B (en) | 2019-05-03 |
Family
ID=59534096
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710098338.9A Active CN107037413B (en) | 2017-02-22 | 2017-02-22 | The direct wave suppressing method offseted based on adaptive Doppler domain wave beam |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107037413B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109901132A (en) * | 2019-03-24 | 2019-06-18 | 西安电子科技大学 | A kind of adaptive coherent Beamforming Method that multipath utilizes |
| CN112731305A (en) * | 2020-12-17 | 2021-04-30 | 国网四川省电力公司信息通信公司 | Direct wave suppression method and system based on adaptive Doppler domain beam cancellation |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4920347A (en) * | 1988-05-07 | 1990-04-24 | Mitsubishi Denki Kabushiki Kaisha | Pulse doppler radar system |
| CN104076343A (en) * | 2014-06-25 | 2014-10-01 | 西安电子科技大学 | Satellite-borne three-channel SAR-GMTI self-adaptive clutter suppression method |
| CN104777466A (en) * | 2015-03-23 | 2015-07-15 | 河海大学 | Passive bistatic radar signal processing method based on dimension reduction constant modulus blind equalization |
-
2017
- 2017-02-22 CN CN201710098338.9A patent/CN107037413B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4920347A (en) * | 1988-05-07 | 1990-04-24 | Mitsubishi Denki Kabushiki Kaisha | Pulse doppler radar system |
| CN104076343A (en) * | 2014-06-25 | 2014-10-01 | 西安电子科技大学 | Satellite-borne three-channel SAR-GMTI self-adaptive clutter suppression method |
| CN104076343B (en) * | 2014-06-25 | 2017-02-15 | 西安电子科技大学 | Satellite-borne three-channel SAR-GMTI self-adaptive clutter suppression method |
| CN104777466A (en) * | 2015-03-23 | 2015-07-15 | 河海大学 | Passive bistatic radar signal processing method based on dimension reduction constant modulus blind equalization |
| CN104777466B (en) * | 2015-03-23 | 2017-05-03 | 河海大学 | Passive bistatic radar signal processing method based on dimension reduction constant modulus blind equalization |
Non-Patent Citations (3)
| Title |
|---|
| FABRIZIO.G.A等: ""Robust adaptive beamforming for HF surface wave over-the-horizon radar"", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 * |
| ROBERTA CARDINALI等: ""Multipath cancellation on reference antenna for passive radar which exploits FM transmission"", 《2007 IET INTERNATIONAL CONFERENCE ON RADAR SYSTEMS》 * |
| 祝俊等: ""基于干扰抑制的恒模波束跟踪新算法"", 《电讯技术》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109901132A (en) * | 2019-03-24 | 2019-06-18 | 西安电子科技大学 | A kind of adaptive coherent Beamforming Method that multipath utilizes |
| CN109901132B (en) * | 2019-03-24 | 2022-12-02 | 西安电子科技大学 | Self-adaptive coherent beam forming method for multi-path utilization |
| CN112731305A (en) * | 2020-12-17 | 2021-04-30 | 国网四川省电力公司信息通信公司 | Direct wave suppression method and system based on adaptive Doppler domain beam cancellation |
| CN112731305B (en) * | 2020-12-17 | 2024-05-03 | 国网四川省电力公司信息通信公司 | Direct wave inhibition method and system based on self-adaptive Doppler domain beam cancellation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107037413B (en) | 2019-05-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103353591B (en) | Bistatic radar localization dimension reduction clutter suppression method based on MIMO | |
| CN104237883B (en) | Airborne radar space time self-adaptation processing method with sparse representation | |
| CN103353592B (en) | Bistatic radar multichannel combination dimension reduction clutter suppression method based on MIMO | |
| CN104360325B (en) | Space-time adaptive processing method for airborne forward-looking array radar | |
| CN106249219B (en) | SAR moving target detection methods based on adaptive matched filter response | |
| CN105807267B (en) | A kind of MIMO radar extends mesh object detection method | |
| CN105223560B (en) | Airborne radar object detection method based on the sparse recovery of clutter pitching azimuth spectrum | |
| CN103760529B (en) | Efficient cascading space-time adaptive processing method based on passive detection | |
| CN101881822B (en) | Method for inhibiting same frequency interference of shared-spectrum radars | |
| CN103176168B (en) | A kind of airborne non-working side battle array radar short range clutter cancellation method | |
| CN103605108B (en) | High-precision remote direction estimation method of acoustic vector array | |
| CN105510887B (en) | To the active suppressing formula disturbance restraining method of airborne radar under a kind of clutter background | |
| CN109375179A (en) | A kind of strong clutter suppression method of airborne radar short range based on rarefaction representation | |
| CN102156279A (en) | Method for detecting moving target on ground by utilizing bistatic radar based on MIMO (Multiple Input Multiple Output) | |
| CN106802402B (en) | DOA estimation method based on dual-layer Parallel circular array antenna | |
| CN109765562A (en) | A kind of three-dimensional looking forward sound sonar system and method | |
| CN107561512A (en) | A kind of polarization of pulse Doppler radar resistance to compression standard towing interference offsets method | |
| CN109633591A (en) | External illuminators-based radar is biradical away from localization method under a kind of observation station location error | |
| CN112596033A (en) | Array element-pulse coding MIMO radar deception jamming suppression method | |
| CN105319545A (en) | MIMO radar waveform design method for improving STAP detection performance | |
| CN106707247A (en) | High-frequency ocean radar target detection method based on compact antenna array | |
| CN109725285A (en) | A kind of DOA estimation method based on the adaptive phase angle conversion of MVDR covariance matrix element | |
| CN103731189B (en) | The dynamic Antenna Subarray Division of conformal array antenna and Wave arrival direction estimating method | |
| CN104020459A (en) | Waveform optimization method for improving MIMO-STAP detection performance | |
| CN108631851A (en) | The Adaptive beamformer method deepened based on uniform linear array null |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |