CN104635214B - Air-borne Forward-looking frequency diversity array radar range ambiguity clutter suppression method - Google Patents
Air-borne Forward-looking frequency diversity array radar range ambiguity clutter suppression method Download PDFInfo
- Publication number
- CN104635214B CN104635214B CN201510076007.6A CN201510076007A CN104635214B CN 104635214 B CN104635214 B CN 104635214B CN 201510076007 A CN201510076007 A CN 201510076007A CN 104635214 B CN104635214 B CN 104635214B
- Authority
- CN
- China
- Prior art keywords
- distance
- frequency diversity
- frequency
- diversity array
- representing
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000001629 suppression Effects 0.000 title claims abstract description 16
- 239000013598 vector Substances 0.000 claims abstract description 67
- 230000001419 dependent effect Effects 0.000 claims description 19
- 238000012545 processing Methods 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 abstract description 5
- 230000008569 process Effects 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 24
- 238000001228 spectrum Methods 0.000 description 24
- 238000004088 simulation Methods 0.000 description 19
- 230000008878 coupling Effects 0.000 description 14
- 238000010168 coupling process Methods 0.000 description 14
- 238000005859 coupling reaction Methods 0.000 description 14
- 230000001427 coherent effect Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000000926 separation method Methods 0.000 description 3
- 230000003595 spectral effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
Landscapes
- Engineering & Computer Science (AREA)
- Radar Systems Or Details Thereof (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
Abstract
The invention belongs to radar clutter suppression technology field, particularly to Air-borne Forward-looking frequency diversity array radar range ambiguity clutter suppression method, it concretely comprises the following steps: fast beat of data x when obtaining the echo sky of Air-borne Forward-looking frequency diversity array radarc;Determine that the secondary range in unambiguous distance district relies on the secondary range in compensation vector and fuzzy distance district and relies on the secondary range that compensation vector draws in Doppler's Combined Frequency district, spatial domain and rely on compensation vector and show that fast beat of data carries out space-time adaptive process to unambiguous distance district fast beat of data and fuzzy distance district fast beat of data when secondary range relies on the echo sky after compensating when secondary range relies on the echo sky after compensating respectively when secondary range relies on the echo sky after compensating.
Description
Technical Field
The invention belongs to the technical field of radar clutter suppression, and particularly relates to a method for suppressing range ambiguity clutter of an airborne forward-looking frequency diversity array radar, which is mainly used for distinguishing the clutter in a non-ambiguity range region from the clutter in a spatial frequency domain.
Background
Space-time adaptive processing (STAP) has important application value in an onboard early warning radar, and can detect a target from clutter and interference by combining time and space information. In the airborne positive side-looking array radar, the coupling relation between the spatial frequency and the Doppler frequency of clutter is not changed along with the distance. However, in some special airborne radar applications (such as forward looking array airborne radar, circular array airborne radar, etc.), the distance invariance of the clutter distribution is no longer true, and the clutter distribution has a severe distance dependency. The range dependence of clutter will lead to a severe degradation of STAP method performance. When range ambiguity occurs, the slow target needs to compete not only with the unambiguous range clutter but also with the ambiguous range clutter, and the low signal-to-noise ratio signal will be completely submerged in the clutter. Fig. 1 shows a range ambiguity diagram of an airborne radar in the prior art. Referring to fig. 1, the clutter of the unambiguous range region and the ambiguity range region are mixed together in the time domain, and it should be noted that the clutter characteristics of the two regions are very different. Therefore, the conventional STAP method is severely degraded in such severe cases.
Disclosure of Invention
The invention aims to provide a method for suppressing distance fuzzy clutter of an airborne forward-looking frequency diversity array radar, explores the controllable degree of freedom of a Frequency Diversity Array (FDA) in a distance dimension, constructs a signal model of FDA-STAP (frequency diversity array-space-time adaptive processing), and provides a method for suppressing the distance fuzzy clutter. The frequency diversity array is adopted, so that the array guide vector is a function of the distance, and therefore, the distance information provided by the array guide vector can be utilized to realize the separation of the clutter of the unambiguous distance area and the ambiguity distance area. Considering that the guide vector can generate quadratic distance dependency, the invention provides a quadratic distance dependency compensation method. The clutter of the unambiguous distance region and the clutter of the ambiguous distance region are distinguished in a space frequency domain, secondary distance dependence compensation is combined with a traditional clutter spectrum compensation method, and the problems of distance dependence and distance ambiguity are solved.
In order to achieve the technical purpose, the invention is realized by adopting the following technical scheme.
The method for suppressing the range ambiguity clutter of the airborne forward-looking frequency diversity array radar comprises the following steps:
step 1, transmitting signals by using each array element of a frequency diversity array, and receiving corresponding echo signals by using the frequency diversity array; obtaining echo space-time snapshot data x of airborne forward-looking frequency diversity array radarc;
Step 2, determining the quadratic distance dependent compensation vector of the unambiguous distance area asAnd the quadratic distance-dependent compensation vector of the blurred distance region is represented asDeriving quadratic distance-dependent compensation vectors in the spatial domain-Doppler joint frequency region Wherein,representing the Kronecker product, α being 1 or 2; 1KColumn vectors of N rows, 1KAll elements in the array are 1, and N represents the array element number in the frequency diversity array;
step 3, obtaining the secondary distance of the non-fuzzy distance areaRelying on compensated echo space-time snapshot dataEcho space-time snapshot data of fuzzy distance area after secondary distance dependence compensation ⊙ represents Hadamard product;
step 4, respectively carrying out echo space-time snapshot data after secondary distance dependence compensation on unambiguous distance areasEcho space-time snapshot data of fuzzy distance area after secondary distance dependence compensationAnd performing space-time self-adaptive processing.
The invention is characterized by further improvement:
in step 1, echo space-time snapshot data x of airborne forward-looking frequency diversity array radarcComprises the following steps:
wherein, M represents the distance unit number of the airborne forward-looking frequency diversity array radar, and M is 1,2.. M; n is a radical ofaDenotes the number of scattering points in each distance unit, i 1,2a;γ{i,m}Represents the complex scattering coefficient of the ith scattering point in the mth range unit,representing the Kronecker product, ⊙ representing the Hadamard product;represents the normalized doppler frequency of the ith scattering point in the mth range bin,representing a time domain guide vector corresponding to the ith scattering point in the mth distance unit; f. ofR {m}Representing the normalized distance frequency, s, of each scattering point in the m-th range binR(fR {m}) Representing a distance guide vector corresponding to each scattering point in the mth distance unit; f. ofa {i,m}Representing the normalized directional frequency, s, of the ith scattering point in the mth range bina(fa {i,m}) Representing a direction guide vector corresponding to the ith scattering point in the mth distance unit; f. ofs {i,m}Representing the normalized spatial frequency, s, of the ith scattering point in the mth range bins(fs {i,m}) And representing the full-space steering vector corresponding to the ith scattering point in the mth distance unit.
In the step 2, the process is carried out,comprises the following steps:
wherein the superscript T represents a matrix or vector means,Δ f denotes the frequency increment of the frequency diversity array, c denotes the speed of light, and R is the speed of light when α is equal to 1c=R1,R1For setting the distance value, when α is equal to 2, Rc=R2=R1+Ru,Ru=c/2fr,frRepresenting the pulse repetition frequency of the airborne forward looking frequency diversity array radar transmitted signal.
The array element spacing d of the frequency diversity array satisfies
Wherein λ is0Which represents the carrier wavelength of the reference array element transmission signal in the frequency diversity array.
The invention has the beneficial effects that: 1) the distance dependence characteristic of an airspace and the quadratic dependence relation of space-time two-dimensional clutter of the forward-looking frequency diversity array radar are analyzed and indicated; 2) a new quadratic distance dependence compensation method is provided; 3) the quadratic distance dependence compensation and the traditional clutter spectrum compensation method are combined, so that distance fuzzy clutter is effectively suppressed. 4) The forward looking frequency diversity array radar can separate the clutter in the unambiguous distance region from the clutter in the ambiguous distance region in the airspace frequency domain, so that the non-uniform clutter suppression performance can be greatly improved.
Drawings
FIG. 1 is a prior art airborne radar range ambiguity diagram;
FIG. 2 is a schematic diagram of a geometric model and an echo receiving process of the airborne forward-looking frequency diversity array radar of the present invention;
FIG. 3 is a schematic diagram of the relationship between the spatial frequency of the scattering point of the airborne forward-looking frequency diversity array radar and the distance and angle;
FIG. 4 is a schematic diagram of a clutter space Doppler coupling relationship of an airborne forward-looking frequency diversity array radar of the present invention in the presence of range ambiguity;
FIG. 5a is a schematic diagram of a clutter space Doppler coupling relationship of a space Doppler frequency coupling relationship obtained by applying quadratic dependence compensation proposed by the present invention to an unambiguous distance region;
FIG. 5b is a schematic diagram of a clutter space Doppler coupling relationship of a space Doppler frequency coupling relationship obtained by applying quadratic dependence compensation proposed by the present invention to a fuzzy distance region;
FIG. 6 is a schematic flow chart of a method for suppressing range ambiguity clutter of an airborne forward-looking frequency diversity array radar according to the present invention;
FIG. 7a is a schematic diagram of the original clutter spectrum distribution when detecting a moving target in a non-blurred distance region in a simulation experiment;
FIG. 7b is a schematic diagram of clutter spectrum distribution after secondary distance-dependent compensation obtained by the present invention when detecting a moving target in a non-ambiguous distance region in a simulation experiment;
FIG. 7c is a schematic diagram of a clutter spectrum compensated by a conventional clutter spectrum when detecting a moving target in a non-ambiguity distance region in a simulation experiment;
FIG. 8a is a schematic diagram of the original clutter spectrum distribution when detecting a moving target in a fuzzy distance region in a simulation experiment;
FIG. 8b is a schematic diagram of clutter spectrum distribution after secondary distance-dependent compensation obtained by the present invention when detecting a moving target in a fuzzy distance region in a simulation experiment;
fig. 8c is a schematic diagram of a clutter spectrum after performing conventional clutter spectrum compensation when detecting a moving target in a fuzzy distance region in a simulation experiment.
Detailed Description
The invention will be further described with reference to the accompanying drawings in which:
referring to fig. 2, a schematic diagram of a geometric model and an echo receiving process of the airborne forward looking frequency diversity array radar of the present invention is shown. The height (platform height) of the airborne frequency diversity array radar is H, the speed of the airborne frequency diversity array radar is V, the number of coherent pulses of the received signals of the airborne forward-looking frequency diversity array radar is K, and the interval of the coherent pulses of the transmitted signals of the airborne forward-looking frequency diversity array radar is TrPulse repetition frequency f of radar transmitted signals of airborne forward-looking frequency diversity arrayr=1/Tr。
In the airborne front video frequency diversity array radar, a frequency diversity array is a uniform linear array consisting of N omnidirectional array elements, and the array element interval in the frequency diversity array is d. Carrier frequency f of nth array element transmitting signal in frequency diversity arraynComprises the following steps:
fn=f0+(n-1)Δf,n=1,2,…,N (1)
wherein f is0Δ f is the frequency increment of the frequency diversity array for the reference frequency of the frequency diversity array. Under the premise, the method for suppressing the range ambiguity clutter of the airborne forward-looking frequency diversity array radar of the inventionThe method comprises the following steps:
step 1, transmitting signals by using each array element of a frequency diversity array, and receiving corresponding echo signals by using the frequency diversity array; obtaining echo space-time snapshot data x of airborne forward-looking frequency diversity array radarc。
The method comprises the following specific steps:
transmitting signals by using each array element of the frequency diversity array, and receiving corresponding echo signals by using the frequency diversity array; in the embodiment of the invention, the frequency diversity array is a uniform linear array with N array elements, and the reference frequency (the carrier frequency of the 1 st array element transmitting signal of the frequency diversity array) of the frequency diversity array is f0The frequency increment of the frequency diversity array is Δ f. The waveforms transmitted by any two array elements of the frequency diversity array are orthogonal to each other, and the signals transmitted by each array element of the frequency diversity array can be effectively separated in the received echo data. To achieve sufficient signal gain, receive beamforming is performed on each transmit waveform in the echo data as the frequency diversity array receives the signal. This step can also be described as being transmitted by a single array element and received by a full aperture array.
For any ground scattering point, the slope from the ground scattering point to the nth array element of the frequency diversity array is expressed asWhere N is 1,2, …, N, theta denotes the azimuth angle of the corresponding ground scatter point with respect to the frequency diversity array,the pitch angle of the corresponding ground scattering point relative to the frequency diversity array is shown, R represents a given reference pitch (e.g., the pitch from the ground scattering point to the nth array element), and d represents the array element spacing in the frequency diversity array.
The phase psi of the echo signal received by the frequency diversity array (full array aperture) after the transmitted signal of the nth array element of the frequency diversity array is scatterednComprises the following steps:
wherein c represents the speed of light, fnThe carrier frequency of the signal transmitted by the nth array element in the frequency diversity array is shown. Taking the 1 st array element in the frequency diversity array as the reference array element, the phase difference delta psi between the nth array element and the reference array element in the frequency diversity arraynComprises the following steps:
wherein f is0Representing a reference frequency of said frequency diversity array, af representing a frequency increment of said frequency diversity array, λ0Representing the carrier wavelength, λ, of a reference array element transmission signal in a frequency diversity array0=c/f0. As can be seen from the content on the right side of the last equal sign of equation (3), the first term is a function of distance and frequency increment, the second term is the same as the conventional phased array, and the third term is a quadratic modulation term. In practice, the third term can be ignored since the carrier frequency is negligible compared to the frequency increment. It can be seen that Δ ψ in equation (3) is different from that of the conventional phased arraynBoth in relation to angle and distance.
For convenience of consideration, the invention only considers the horizontal speed of the motion platform, and then the normalized Doppler frequency f of the signal reflected by the corresponding ground scatterer from the transmission signal of the nth array element of the frequency diversity arraytComprises the following steps:
wherein V represents the speed of the carrier, frRepresenting the pulse repetition frequency of the airborne forward looking frequency diversity array radar transmitted signal. From equation (4), it can be seen that the frequency diversity array is notIn the same array element, the normalized doppler frequencies of the signals reflected back by the corresponding ground scatterers are slightly different, however, in general, the difference can be ignored due to the small frequency increment of the frequency diversity array. In the embodiment of the present invention, the transmitted signal of the nth array element of the frequency diversity array is scattered, and then the k pulse echo signal x is received by the frequency diversity array (full array aperture)nkComprises the following steps:
wherein, γpFor complex coefficients, gamma, independent of the carrier frequencypIs a set complex coefficient. ZetakIs composed of
The echo signal x of the corresponding scattering point received by the airborne forward-looking frequency diversity array radar is:
wherein,representing the Kronecker product, ⊙ representing the Hadamard product, ftIndicating the normalized Doppler frequency, f, of the corresponding scattering pointsRepresenting the normalized spatial frequency, f, of the corresponding scattering pointRRepresenting the normalized distance frequency, f, of the corresponding scattering pointaRepresenting the normalized directional frequency of the corresponding scattering point. f. oft、fs、fRAnd faRespectively as follows:
fs=fR+fa;
strepresenting a time-domain steering vector and,stthe column vector of the K rows represents the coherent pulse number of the received signal of the airborne forward-looking frequency diversity array radar; sRA distance-oriented vector is represented which is,sRis a column vector of N rows, N representing the number of array elements of the frequency diversity array. saWhich represents a directional steering vector, is,saa column vector of N rows. ssIs a full-airspace guide vector synthesized by a distance guide vector and a direction guide vector,ssa column vector of N rows. st(ft)、sR(fR)、sa(fa) And ss(fs) Respectively as follows:
where the superscript T represents a matrix or vector arrangement. It can be seen that st(ft) The same time domain steering vector as that of the traditional phased array radar, sa(fa) The same as the directional steering vector of a conventional phased array radar. But sa(fa) Different from the full airspace guide vector of the traditional phased array radar, sR(fR) Depending on the ramp distance R and the frequency increment af of the frequency diversity array.
In the embodiment of the invention, the echo signal of the airborne forward-looking frequency diversity array radar is formed by scattering and superposing clutter in an equidistant ring, so that the echo space-time snapshot data x of the airborne forward-looking frequency diversity array radarcComprises the following steps:
where M denotes the number of range units (fuzzy range number), M being 1,2.. M; n is a radical ofaDenotes the number of independent scattering points in each distance unit, i 1,2a。γ{i,m}Represents the complex scattering coefficient (as a known quantity) of the ith scattering point in the mth range bin,representing the Kronecker product, ⊙ representing the Hadamard product, ft {i,m}Representing the normalized Doppler frequency, s, of the ith scattering point in the mth range bint(ft {i,m}) Representing the time domain steering vector corresponding to the ith scattering point in the mth distance unit,st(ft {i,m}) And the column vector of K rows represents the coherent pulse number of the signals received by the airborne forward-looking frequency diversity array radar. f. ofR {m}Representing the normalized distance frequency, s, of each scattering point in the m-th range binR(fR {m}) Representing the range steering vector corresponding to each scattering point in the mth range bin,sR(fR {m}) Is a column vector of N rows, N representing the number of array elements of the frequency diversity array. f. ofa {i,m}Representing the normalized directional frequency, s, of the ith scattering point in the mth range bina(fa {i,m}) Representing the corresponding square of the ith scattering point in the mth range unitA direction vector is guided to the direction vector,sa(fa {i,m}) A column vector of N rows. f. ofs {i,m}Representing the normalized spatial frequency, s, of the ith scattering point in the mth range bins(fs {i,m}) Representing the full-space steering vector corresponding to the ith scattering point in the mth distance unit,ss(fs {i,m}) A column vector of N rows.
In the examples of the present invention, ft {i,m}、fR {m}、fa {i,m}And fs {i,m}Are respectively:
wherein V represents the speed of the carrier, frRepresenting the pulse repetition frequency, λ, of radar transmitted signals of an airborne forward-looking frequency diversity array0Representing the carrier wavelength, λ, of a reference array element transmission signal in a frequency diversity array0=c/f0,f0A reference frequency representing the frequency diversity array; Δ f represents the frequency increment of the frequency diversity array, c represents the speed of light, and d represents the array element spacing in the frequency diversity array. Theta{i,m}Represents the azimuth angle of the ith scattering point in the mth range bin relative to the frequency diversity array,represents the elevation angle of the ith scattering point in the mth range bin relative to the frequency diversity array,R{m}a given reference slope corresponding to each scattering point in the mth range bin (e.g., the slope from any scattering point in the mth range bin to the nth array element) is indicated.
In the examples of the present invention, st(ft {i,m})、sR(fR {m})、sa(fa {i,m}) And ss(fs {i,m}) Are respectively:
where the superscript T represents a matrix or vector arrangement.
Step 2, determining the quadratic distance dependent compensation vector of the unambiguous distance area asAnd the quadratic distance-dependent compensation vector of the blurred distance region is represented asDeriving quadratic distance-dependent compensation vectors in the spatial domain-Doppler joint frequency region
Wherein,denotes the Kronecker product, α ═ 1, 2; 1KColumn vectors of N rows, 1KAll elements in (1).
The method comprises the following specific steps:
for ground scatter points, the full-space steering vector consists of two parts: a direction guide vector and a distance guide vector. For an airborne forward-looking frequency diversity array radar, the direction steering vector of the ground scattering point is the same as that of a traditional phased array radar. It should be noted that the distance guide vector is a function of distance and frequency increment, and will correspond to the normalized spatial frequency f of the ground scattering pointsRewritten as:
from equation (9), the spatial frequency depends on the distance and angle.
Referring to fig. 3, it is a schematic diagram of the relationship between the spatial frequency of the scattering point of the airborne forward-looking frequency diversity array radar and the distance and angle. As shown in FIG. 3, the solid line represents a distance R1With the dashed line representing a distance R2Wave front of, R1And R2Are given values; r1-R2Δ R. For a given direction, the spatial frequency of the ground scattering point of the airborne forward-looking frequency diversity array radar varies with distance.
In the embodiment of the invention, the airborne forward-looking frequency diversity array radar relates to R1Normalization of ground scattering pointsChange the spatial frequency fs(R1) Is shown as
Airborne forward looking frequency diversity array radar with respect to R2Normalized spatial frequency f of ground scattering pointss(R2) Is shown as
Then there are:
equation (12) gives the spatial frequency difference in dependence on the ramp distance difference Δ R at a given direction and frequency increment. Therefore, a more problematic distance dependence problem, i.e., a secondary distance dependence problem, is encountered in the airborne forward looking frequency diversity array radar of the present invention.
When the ambiguity distance occurs, the spatial frequency difference of the unambiguous distance and the ambiguity distance can be expressed as follows:
wherein R is2=R1+Ru,Ru=c/2fr,RuRepresenting the maximum unambiguous distance. R1Is a given value. Although clutter in the unambiguous distance region and the ambiguous distance region cannot be separated in the time domain, the airborne forward-looking frequency diversity array radar can still utilize the spatial frequency difference DeltafsThey are separated. This also verifies the superiority of the frequency diversity array over conventional phased arrays.
In the airborne forward-looking frequency diversity array radar, the angle Doppler frequency relation of the clutter corresponding to the ground scattering point is expressed as
Wherein,h is the height of the carrier. Equation (14) shows the elliptical coupling relationship of the clutter distribution. Referring to fig. 4, a schematic diagram of a clutter space doppler coupling relationship of the airborne forward looking frequency diversity array radar in the presence of range ambiguity is shown. In fig. 4, the horizontal axis represents normalized doppler frequency (normalized doppler frequency) and the vertical axis represents normalized spatial frequency (normalized spatial frequency), different lines represent different distance elements, respectively, with the unit m, A, B, C representing three different scattering points, and R represents the slant range. It can be seen that when the slant distance changes, the curve of the elliptical coupling relationship (clutter space Doppler coupling relationship) also follows the space frequency fRAnd changes accordingly. In addition, the ellipse increases with increasing slope distance.
Therefore, unlike the traditional phased array radar, the clutter distribution of the airborne forward looking frequency diversity array radar of the invention changes along with the change of the slant range in the airspace frequency domain. Therefore, the airborne forward looking frequency diversity array radar can realize the separation of the unambiguous range clutter and the ambiguous range clutter in the airspace frequency domain. It should be noted that although the Frequency Diversity Array (FDA) provides additional information in the spatial and frequency domains, it also causes quadratic distance dependence problems.
The airborne forward-looking frequency diversity array radar can provide additional information in a space frequency domain, and an effective distance fuzzy clutter suppression method is discussed in the invention. Firstly, in order to overcome the problem of secondary distance dependence, a secondary distance dependence compensation method is provided; two groups of compensation vectors are respectively constructed according to the fuzzy distance area and the fuzzy distance area, and then distance fuzzy clutter suppression is realized by combining a traditional clutter spectrum compensation method.
Specifically, in the embodiment of the present invention, for a given distance (the distance units are the same), the quadratic distance-dependent compensation vector of the unambiguous distance region is expressed asThe quadratic distance-dependent compensation vector of the blurred distance region is represented asLet α be 1,2, then there are:
wherein the superscript T represents a matrix or vector means, a column vector of N rows. f. ofcWhich is indicative of the frequency of the compensation,Δ f represents the frequency increment of the frequency diversity array, c represents the speed of light, RcIs corresponding slope distance, when α is equal to 1, Rc=R1When α is equal to 2, RcFor fuzzy distance, Rc=R2=R1+Ru,RuDenotes the maximum unambiguous distance, Ru=c/2fr,frRepresenting the pulse repetition frequency of the airborne forward looking frequency diversity array radar transmitted signal.
In the embodiment of the invention, the compensation vectorThe quadratic distance dependence of clutter can be compensated. After the secondary distance dependence compensation processing, the clutter spectrum meets the independent same distribution characteristic in an airspace.
In the embodiment of the invention, quadratic distance dependence compensation vector in airspace-Doppler joint frequency regionComprises the following steps:
wherein,denotes the Kronecker product, α ═ 1, 2; 1KColumn vectors of N rows, 1KAll elements in (1).
Step 3, obtaining echo space-time snapshot data of the unambiguous distance area after secondary distance dependence compensationEcho space-time snapshot data of fuzzy distance area after secondary distance dependence compensation Wherein ⊙ represents the Hadamard product.
The method comprises the following specific steps:
in an on-board front video frequency diversity array radar, diversityRespectively carrying out secondary distance dependence compensation on the echo space-time snapshot data in the unambiguous distance region and the ambiguous distance region to obtain the echo space-time snapshot data in the unambiguous distance region after the secondary distance dependence compensationEcho space-time snapshot data of fuzzy distance area after secondary distance dependence compensation
Let α be 1,2, then:
wherein, the first and second guide rollers are arranged in a row, a column vector of N rows. f. ofcWhich is indicative of the frequency of the compensation,Δ f represents the frequency increment of the frequency diversity array, c represents the speed of light, RcIs corresponding slope distance, when α is equal to 1, Rc=R1When α is equal to 2, RcFor fuzzy distance, Rc=R2=R1+Ru,RuDenotes the maximum unambiguous distance, Ru=c/2fr,frRepresenting the pulse repetition frequency of the airborne forward looking frequency diversity array radar transmitted signal.Is a compensated spatial domain steering vector and,
for the corresponding compensated spatial frequencies:
as is clear from the above, the angular doppler frequency relation of the clutter in equation (14) can be rewritten as follows:
when α is 1, there are:
when α is 2, there are:
referring to fig. 5a, a clutter space doppler coupling relationship schematic diagram of a space doppler frequency coupling relationship obtained after the secondary dependence compensation proposed by the present invention is applied to a non-ambiguity distance region; referring to fig. 5b, a clutter space doppler coupling relationship diagram of a space doppler frequency coupling relationship obtained after the secondary dependence compensation proposed by the present invention is applied to the fuzzy distance region; in fig. 5a and 5b, the horizontal axis represents the normalized doppler frequency (normalized doppler frequency) and the vertical axis represents the normalized spatial frequency (normalized spatial frequency), and different lines represent different distance elements, respectively. As can be seen from fig. 5a and 5b, clutter without and in the ambiguity range can be effectively resolved in the spatial frequency domain. As shown in fig. 5a, from the spatial frequency domain, the clutter distribution of the unambiguous range region of the airborne forward looking frequency diversity array radar will be aligned to the middle part, and the ambiguous range region will be aligned to both sides. In contrast, as shown in fig. 5(b), the clutter of the unambiguous range region is located on both sides, and the clutter of the ambiguous range region is located in the middle of the spatial frequency domain.
In addition, in the embodiment of the invention, in order to ensure the effective classification of the unambiguous distance clutter and the ambiguous distance clutter, the array element spacing d of the frequency diversity array of the invention should satisfy
Wherein λ is0Representing the carrier wavelength, λ, of a reference array element transmission signal in a frequency diversity array0=c/f0,f0Representing a reference frequency of the frequency diversity array. Without loss of generality, the array element space d is usually a quarter of a wavelength, i.e., d is 0.25 λ0。
Step 4, respectively carrying out echo space-time snapshot data after secondary distance dependence compensation on unambiguous distance areasEcho space-time snapshot data of fuzzy distance area after secondary distance dependence compensationPerforming space-time self-adaptive processing to obtain echo space-time snapshot data of the fuzzy distance area after the corresponding clutter suppression after the secondary distance dependence compensation and echo space-time snapshot data of the fuzzy distance area after the corresponding clutter suppression after the secondary distance dependence compensationAnd (4) data.
In step 4, the procedure of the space-time adaptive processing is well known to those skilled in the art and will not be described in detail herein.
In summary, the space-time adaptive processing (STAP) performance of conventional phased array radars may be severely degraded when range ambiguities are present. Firstly, due to the fact that clutter non-stationarity is serious due to distance ambiguity, samples required by clutter covariance matrix estimation are difficult to obtain; secondly, in a non-front side-view geometrical structure, the compensation of the unambiguous range clutter and the fuzzy range clutter is mutually restricted, so that the traditional clutter spectrum compensation method is invalid; finally, clutter freedom degree is greatly increased due to distance ambiguity, and burden of system clutter suppression is increased. In the embodiment of the invention, after secondary distance dependent compensation (clutter spectrum compensation), the airborne forward-looking frequency diversity array radar realizes the separation of distance fuzzy clutter in an airspace frequency domain, the clutter spectrum compensation of a non-fuzzy distance area and a fuzzy distance area is not mutually restricted any more, and the clutter compensation and the clutter suppression can be realized by utilizing a traditional clutter spectrum compensation method. Referring to fig. 6, a schematic flow chart of the method for suppressing the airborne forward looking frequency diversity array radar range ambiguity clutter according to the present invention is shown.
The effect of the present invention can be further illustrated by the following simulation results:
1) simulation conditions
TABLE 1 simulation parameters
In the simulation experiment, the test distance of the non-fuzzy distance area is set to be 8000m, and the test distance of the fuzzy distance area is set to be 18000 m.
2) Emulated content
Simulation experiment 1: the target is set in the unambiguous distance zone. Referring to fig. 7a, it is a schematic diagram of the original impurity spectrum distribution when detecting a moving target in a non-fuzzy distance region in a simulation experiment; fig. 7b is a schematic diagram of clutter spectrum distribution after secondary distance-dependent compensation obtained by the method when a moving target is detected in a non-fuzzy distance area in a simulation experiment; referring to fig. 7c, a clutter spectrum schematic diagram after conventional clutter spectrum compensation is performed when a moving target is detected in a non-fuzzy distance region in a simulation experiment; in fig. 7a to 7c, the horizontal axis represents the normalized doppler frequency and the vertical axis represents the normalized spatial frequency, and different gray values represent different clutter spectral amplitudes. As can be seen from fig. 7a to 7c, in the on-board front video frequency diversity array radar, the unambiguous range clutter and the ambiguous range clutter in the spatial frequency region can be discriminated. Due to quadratic distance dependence, clutter is heavily spread in the spatial frequency domain. After the quadratic distance dependent compensation and clutter compensation are carried out by using the invention, the distribution of clutter in an unambiguous region in a spatial frequency domain is [ -1/4,1/4], and the distribution of clutter in an ambiguous distance region is [ -1/2, -1/4 ]. U [1/4,1/2 ]. At this time, the clutter in the fuzzy distance region does not have much influence on the target detection in the non-fuzzy distance region.
Simulation experiment 2: the target is set to be in the fuzzy distance zone. In fact, the distance dependence of clutter varies little in ambiguous distance regions, but varies significantly in unambiguous distance regions. Referring to fig. 8a, a schematic diagram of original clutter spectrum distribution when detecting a moving target in a fuzzy distance region in a simulation experiment; referring to fig. 8b, a schematic diagram of clutter spectrum distribution after secondary distance-dependent compensation obtained by the present invention when detecting a moving target in a fuzzy distance region in a simulation experiment is shown; referring to fig. 8c, a clutter spectrum schematic diagram after performing conventional clutter spectrum compensation when detecting a moving target in a fuzzy distance region in a simulation experiment is shown; in fig. 8a to 8c, the horizontal axis represents the normalized doppler frequency and the vertical axis represents the normalized spatial frequency, and different gray values represent different clutter spectral amplitudes. As can be seen from fig. 8a to 8c, after the fuzzy distance region clutter is compensated by the present invention, the distribution of the fuzzy distance region clutter is [ -1/4,1/4], the distribution of the non-fuzzy distance region clutter is [ -1/2, -1/4], [ u ], [1/4,1/2], and the distance fuzzy clutter can also be resolved in the spatial frequency domain. It follows that the unambiguous distance clutter and the ambiguous distance clutter after quadratic distance dependence compensation are separable. And by combining the traditional clutter compensation method, the clutter suppression performance can be obviously improved.
In conclusion, the simulation experiment verifies the correctness, reliability and effectiveness of the method.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include such modifications and variations.
Claims (4)
1. The method for suppressing the range ambiguity clutter of the airborne forward-looking frequency diversity array radar is characterized by comprising the following steps of:
step 1, transmitting signals by using each array element of a frequency diversity array, and receiving corresponding echo signals by using the frequency diversity array; obtaining echo space-time snapshot data x of airborne forward-looking frequency diversity array radarc;
Step 2, determining the quadratic distance dependent compensation vector of the unambiguous distance area asAnd the quadratic distance-dependent compensation vector of the blurred distance region is represented asDeriving quadratic distance-dependent compensation vectors in the spatial domain-Doppler joint frequency region Wherein,representing the Kronecker product, α being 1 or 2; 1KColumn vectors of N rows, 1KAll elements in (1), N represents the number of array elements in the frequency diversity array, fcRepresents the compensation frequency;
step 3, obtaining echo space-time snapshot data of the unambiguous distance area after secondary distance dependence compensationEcho space-time snapshot data of fuzzy distance area after secondary distance dependence compensation Wherein ⊙ represents the Hadamard product,representing the echo space-time snapshot data after quadratic distance-dependent compensation, α is 1 or 2, when α equals 1,is a non-fuzzy distance zone twiceThe distance-dependent compensated echo space-time snapshot data is updated, when α is 2,echo space-time snapshot data of the fuzzy distance area after secondary distance dependence compensation;
step 4, respectively carrying out echo space-time snapshot data after secondary distance dependence compensation on unambiguous distance areasEcho space-time snapshot data of fuzzy distance area after secondary distance dependence compensationAnd performing space-time self-adaptive processing.
2. The method for range ambiguity suppression for an airborne forward looking frequency diversity array radar as claimed in claim 1 wherein in step 1, the echo space time snapshot data x for the airborne forward looking frequency diversity array radar is obtainedcComprises the following steps:
wherein, M represents the distance unit number of the airborne forward-looking frequency diversity array radar, and M is 1,2.. M; n is a radical ofaDenotes the number of scattering points in each distance unit, i 1,2a;γ{i,m}Represents the complex scattering coefficient of the ith scattering point in the mth range unit,representing the Kronecker product, ⊙ representing the Hadamard product, ft {i,m}Representing the normalized Doppler frequency, s, of the ith scattering point in the mth range bint(ft {i,m}) Representing a time domain guide vector corresponding to the ith scattering point in the mth distance unit; f. ofR {m}Representing the normalized distance frequency, s, of each scattering point in the m-th range binR(fR {m}) Representing a distance guide vector corresponding to each scattering point in the mth distance unit; f. ofa {i,m}Representing the normalized directional frequency, s, of the ith scattering point in the mth range bina(fa {i,m}) Representing a direction guide vector corresponding to the ith scattering point in the mth distance unit; f. ofs {i,m}Representing the normalized spatial frequency, s, of the ith scattering point in the mth range bin (fs {i,m}) And representing the full-space steering vector corresponding to the ith scattering point in the mth distance unit.
3. The method for range ambiguity suppression for an airborne forward looking frequency diversity array radar according to claim 1, wherein in step 2,comprises the following steps:
wherein the superscript T represents a matrix or vector means,Δ f denotes the frequency increment of the frequency diversity array, c denotes the speed of light, and R is the speed of light when α is equal to 1c=R1,R1For setting the distance value, when α is equal to 2, Rc=R2=R1+Ru,Ru=c/2 fr,frRepresenting the pulse repetition frequency of the airborne forward looking frequency diversity array radar transmitted signal.
4. The method for range ambiguity suppression for an airborne forward looking frequency diversity array radar according to any of claims 1 to 3 wherein the array element spacing d of the frequency diversity array is such that
Wherein λ is0Which represents the carrier wavelength of the reference array element transmission signal in the frequency diversity array.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510076007.6A CN104635214B (en) | 2015-02-12 | 2015-02-12 | Air-borne Forward-looking frequency diversity array radar range ambiguity clutter suppression method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510076007.6A CN104635214B (en) | 2015-02-12 | 2015-02-12 | Air-borne Forward-looking frequency diversity array radar range ambiguity clutter suppression method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104635214A CN104635214A (en) | 2015-05-20 |
CN104635214B true CN104635214B (en) | 2017-01-04 |
Family
ID=53214170
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510076007.6A Expired - Fee Related CN104635214B (en) | 2015-02-12 | 2015-02-12 | Air-borne Forward-looking frequency diversity array radar range ambiguity clutter suppression method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104635214B (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105116385B (en) * | 2015-09-14 | 2017-06-20 | 电子科技大学 | A kind of method of generation FDA radar emission signals, device and FDA radars |
CN106529440B (en) * | 2016-11-02 | 2019-08-13 | 西安电子科技大学 | Coincidence frequency diversity battle array radar segmented matched filter method |
CN106569212B (en) * | 2016-11-09 | 2019-06-18 | 西安空间无线电技术研究所 | A kind of multichannel SAR-GMTI range ambiguity clutter suppression method |
CN108196239B (en) * | 2018-01-18 | 2022-03-18 | 西安电子科技大学 | Unambiguous parameter estimation method for frequency diversity MIMO radar |
CN109471100A (en) * | 2018-10-16 | 2019-03-15 | 湖北航天技术研究院总体设计所 | A kind of SAR doppler frequency rate estimation method and system |
CN109765536B (en) * | 2018-10-22 | 2023-04-21 | 西北大学 | FDA-MIMO dimension-reduction space-time adaptive clutter suppression method and device based on auxiliary channel |
CN114265018B (en) * | 2022-02-28 | 2022-05-31 | 西安电子科技大学 | Short-range clutter suppression method based on multi-frequency split radar |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3781878A (en) * | 1971-07-29 | 1973-12-25 | G Kirkpatrick | Means for suppressing ground clutter in airborne radar |
US7460059B1 (en) * | 2006-10-25 | 2008-12-02 | Sandia Corporation | Removing interfering clutter associated with radar pulses that an airborne radar receives from a radar transponder |
CN102721947A (en) * | 2012-06-14 | 2012-10-10 | 河海大学 | Efficient three-dimensional space-time adaptive clutter suppression method for airborne digital array radar |
CN103605114A (en) * | 2013-12-03 | 2014-02-26 | 西安电子科技大学 | Non-broadside array airborne radar short range clutter suppression method based on multiple frequencies |
CN104076337A (en) * | 2014-07-04 | 2014-10-01 | 西安电子科技大学 | Airborne radar clutter suppression method based on array element amplitude and phase error correction |
-
2015
- 2015-02-12 CN CN201510076007.6A patent/CN104635214B/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3781878A (en) * | 1971-07-29 | 1973-12-25 | G Kirkpatrick | Means for suppressing ground clutter in airborne radar |
US7460059B1 (en) * | 2006-10-25 | 2008-12-02 | Sandia Corporation | Removing interfering clutter associated with radar pulses that an airborne radar receives from a radar transponder |
CN102721947A (en) * | 2012-06-14 | 2012-10-10 | 河海大学 | Efficient three-dimensional space-time adaptive clutter suppression method for airborne digital array radar |
CN103605114A (en) * | 2013-12-03 | 2014-02-26 | 西安电子科技大学 | Non-broadside array airborne radar short range clutter suppression method based on multiple frequencies |
CN104076337A (en) * | 2014-07-04 | 2014-10-01 | 西安电子科技大学 | Airborne radar clutter suppression method based on array element amplitude and phase error correction |
Non-Patent Citations (2)
Title |
---|
一种稳健的机载前视阵雷达杂波谱补偿方法;李明等;《西安电子科技大学学报(自然科学版)》;20090831;第36卷(第4期);第633-638页 * |
前视阵频率分集雷达空时杂波特性研究;胡柏林等;《电子与信息学报》;20131130;第35卷(第11期);第2693-2699页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104635214A (en) | 2015-05-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104635214B (en) | Air-borne Forward-looking frequency diversity array radar range ambiguity clutter suppression method | |
CN106093870B (en) | The SAR-GMTI clutter suppression methods of hypersonic aircraft descending branch | |
CN109116311B (en) | Clutter suppression method based on knowledge-aided sparse iteration covariance estimation | |
CN104360325B (en) | Space-time adaptive processing method for airborne forward-looking array radar | |
CN112444811A (en) | Target detection and imaging method integrating MIMO radar and ISAR | |
US6720910B2 (en) | Pri-staggered post-doppler adaptive monopulse processing for detection and location of a moving target in ground clutter | |
Belfiori et al. | 2D-MUSIC technique applied to a coherent FMCW MIMO radar | |
CN109212489B (en) | FDA-MIMO radar fuzzy clutter suppression method based on auxiliary pulse | |
CN108490443B (en) | Multi-subarray synthetic aperture sonar omega-k imaging algorithm based on analytic solution and NUFFT | |
CN110632587A (en) | Weak moving object monitoring method based on rapid FMCW radar | |
CN108387884B (en) | Airborne radar clutter suppression method based on knowledge-assisted sparse progressive minimum variance | |
CN111007503A (en) | Moving target focusing and positioning method and system based on frequency spectrum accurate positioning | |
Jansen et al. | Practical multichannel SAR imaging in the maritime environment | |
CN107153191B (en) | Double-base ISAR imaging detection method for invisible airplane | |
CN103760540B (en) | Based on moving target detect and the method for parameter estimation of reconstruction signal and 1-norm | |
CN106526544B (en) | MIMOSAR clutter suppression method based on hypersonic platform | |
US11125871B2 (en) | Azimuth estimation device and method | |
CN111630410A (en) | FMCW radar sensor | |
He et al. | An adaptive coding-angle-Doppler clutter suppression approach with extended azimuth phase coding array | |
CN106802408B (en) | Airborne non-positive side array short range clutter distance fuzzy suppression method based on sparse recovery | |
Gao et al. | Static background removal in vehicular radar: Filtering in azimuth-elevation-doppler domain | |
CN110456342B (en) | Far-field multi-moving-object detection method of single-transmitting-antenna radar | |
CN109581366B (en) | Discrete sidelobe clutter identification method based on target steering vector mismatch | |
CN114428228B (en) | Clutter suppression method for high-repetition-frequency sum-difference antenna radar seeker | |
CN115856810A (en) | Method for reducing minimum detectable speed of intensive offset scanning and electronic equipment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170104 Termination date: 20170212 |