CN114265018A - Short-range clutter suppression method based on multi-frequency split radar - Google Patents

Short-range clutter suppression method based on multi-frequency split radar Download PDF

Info

Publication number
CN114265018A
CN114265018A CN202210184291.9A CN202210184291A CN114265018A CN 114265018 A CN114265018 A CN 114265018A CN 202210184291 A CN202210184291 A CN 202210184291A CN 114265018 A CN114265018 A CN 114265018A
Authority
CN
China
Prior art keywords
clutter
data
fuzzy
vector
frequency
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202210184291.9A
Other languages
Chinese (zh)
Other versions
CN114265018B (en
Inventor
王华柯
刘成
全英汇
廖桂生
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xi'an Hesai Electronic Technology Co ltd
Original Assignee
Xidian University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN202210184291.9A priority Critical patent/CN114265018B/en
Publication of CN114265018A publication Critical patent/CN114265018A/en
Application granted granted Critical
Publication of CN114265018B publication Critical patent/CN114265018B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a short-range clutter suppression method based on a multi-frequency split radar, which comprises the following steps: obtaining clutter echo data based on a radar system of superposition step frequency; analyzing two-dimensional distribution according to the receiving space frequency and the Doppler frequency in the clutter echo data; finding out the relation between the emission space frequency and the fuzzy area; constructing a principal value distance compensation vector based on the relation between the emission space frequency and the fuzzy region to obtain compensated clutter echo data; constructing a blocking vector according to the compensated clutter echo data to obtain the blocked data of each fuzzy area; and (4) carrying out simultaneous equations on expressions of the blocked data of each fuzzy region to obtain clutter data of each fuzzy region. The invention has the advantages of only one transmitting array element, hardware cost saving, low calculation complexity, easy realization, low requirement on system parameters, no need of considering the condition of transmitting spatial frequency ambiguity, clutter separation and good inhibition effect.

Description

Short-range clutter suppression method based on multi-frequency split radar
Technical Field
The invention belongs to the technical field of radar signal processing, and particularly relates to a short-range clutter suppression method based on a multi-frequency split system radar.
Background
The airborne radar is usually in the downward-looking work, the clutter is strong, due to the high-speed motion of the platform, the clutter has certain relative speed, the relative speed in different directions is different, and weak low-speed targets are submerged in the strong clutter, so that the performance of the radar is greatly reduced. Ground clutter echoes show correlation in both space and time dimensions, while interference is mainly correlated in angle, and signal to noise ratio (SCNR) can be effectively maximized by using space-time freedom. Space-time adaptive processing (STAP) involves adaptively adjusting a two-dimensional space-time filter response to maximize the output signal-to-noise ratio. By introducing the STAP algorithm into the airborne radar, the radar detection performance is improved.
In an airborne radar, clutter Doppler spread is serious, in order to avoid Doppler ambiguity as far as possible, a signal with a higher pulse repetition frequency is generally transmitted, but as is well known, the maximum unambiguous distance and the pulse repetition frequency are in inverse proportion, the higher pulse repetition frequency can bring a serious distance ambiguity problem, clutter in different ambiguity intervals can be superposed together at the moment, echo data of the clutter can not meet the independent and same-distribution condition any more, the direct application of the STAP can also suppress the notch broadening of the clutter factor, and the suppression performance of the STAP on the clutter can be remarkably reduced. For the non-ortholateral array, the distribution of the clutter has dependence on the distance, the dependence of the short-range clutter on the distance is very serious, and the dependence of the long-range clutter on the distance is relatively light.
Although some existing clutter compensation methods can reduce the distance dependency by compensating clutter data, when distance ambiguity is involved, the enhancement effect of the compensation methods on the STAP of the non-positive side array is not obvious due to different distance dependencies of different fuzzy areas.
Disclosure of Invention
In order to solve the problems in the prior art, the invention provides a short-range clutter suppression method based on a multi-frequency split system radar. The technical problem to be solved by the invention is realized by the following technical scheme:
a short-range clutter suppression method based on a multi-frequency division radar, the short-range clutter suppression method comprising:
step 1, obtaining clutter echo data based on a radar system with superposed step frequency;
step 2, analyzing two-dimensional distribution of the clutter echo data on the receiving space frequency and the Doppler frequency according to the receiving space frequency and the Doppler frequency in the clutter echo data;
step 3, obtaining the relation between the transmitting space frequency and the distance according to the transmitting space frequency of the clutter echo data so as to find out the relation between the transmitting space frequency and the fuzzy area;
step 4, constructing a principal value distance compensation vector based on the relation between the emission space frequency and the fuzzy area to obtain compensated clutter echo data;
step 5, constructing a blocking vector according to the compensated clutter echo data to obtain the blocked data of each fuzzy area;
and 6, establishing a simultaneous equation set for the expressions of the blocked data of the fuzzy areas to obtain the clutter data of each fuzzy area.
In an embodiment of the present invention, the transmitting antennas are configured to transmit M orthogonal signals, each of the orthogonal signals transmits K pulses, there are N receiving antennas in total, and there are P fuzzy areas in total, then the clutter data received by the ith distance unit is:
Figure 497357DEST_PATH_IMAGE001
wherein,
Figure 352180DEST_PATH_IMAGE002
the transmitting space guide vector, the receiving space guide vector and the Doppler guide vector of the clutter echo data are respectively represented, and Nc represents the equivalent number of clutter sources in each range unit.
In one embodiment of the present invention, the two-dimensional distribution is:
Figure 928655DEST_PATH_IMAGE003
Figure 279125DEST_PATH_IMAGE004
wherein,
Figure 31181DEST_PATH_IMAGE005
representing the received spatial frequency and the normalized Doppler frequency, respectively
Figure 791326DEST_PATH_IMAGE006
The space-time two-dimensional steering vector of (a),Ra covariance matrix representing clutter echo data, L represents a total number of training samples, L represents an ith training sample,Hrepresenting a conjugate transpose, and the training samples represent clutter echo data at the ith range bin.
In one embodiment of the invention, the relationship between the transmit spatial frequency and the distance is:
Figure 996043DEST_PATH_IMAGE007
Figure 867047DEST_PATH_IMAGE008
wherein,
Figure 473609DEST_PATH_IMAGE009
which represents the spatial frequencies of the reception of the signal,
Figure 404656DEST_PATH_IMAGE010
the increment of the frequency is represented by,cthe speed of light is indicated and is,Rthe distance is represented as a function of time,
Figure 96668DEST_PATH_IMAGE011
indicating a dominant distance, P indicating the number of blurred regions,
Figure 771363DEST_PATH_IMAGE012
maximum unambiguous distance.
In one embodiment of the present invention, the principal value distance compensation vector is:
Figure 232431DEST_PATH_IMAGE013
where M is the number of transmit waveforms.
In one embodiment of the present invention, the compensated clutter echo data is:
Figure 599959DEST_PATH_IMAGE014
wherein,
Figure 779267DEST_PATH_IMAGE015
representing an identity matrix of dimension N,
Figure 992074DEST_PATH_IMAGE016
representing an identity matrix of dimension K.
In one embodiment of the present invention, the blur area includes a first blur area, a second blur area, a third blur area and a fourth blur area, wherein the step 4 includes:
step 4.1, compensating the dominant value distance of the first fuzzy area through the corresponding dominant value distance compensation vector to obtain clutter echo data after compensation in the first fuzzy area;
step 4.2, compensating the dominant value distance of the second fuzzy area through the corresponding dominant value distance compensation vector to obtain clutter echo data after compensation in the second fuzzy area;
4.3, compensating the dominant value distance of the third fuzzy area through the corresponding dominant value distance compensation vector to obtain clutter echo data after compensation in the third fuzzy area;
and 4.4, compensating the principal value distance of the fourth fuzzy area through the corresponding principal value distance compensation vector to obtain the clutter echo data after compensation in the fourth fuzzy area.
In one embodiment of the present invention, the step 5 comprises:
step 5.1, rearranging the compensated clutter echo data in the first fuzzy area, and constructing a first blocking vector to obtain data after the clutter echo data are blocked in the first fuzzy area;
step 5.2, rearranging the compensated clutter echo data in the second fuzzy area, and constructing a second blocking vector to obtain data after the clutter blocking in the second fuzzy area;
step 5.3, rearranging the compensated clutter echo data in the third fuzzy area, and constructing a third blocking vector to obtain data after the third fuzzy area blocks the clutter;
and 5.4, rearranging the compensated clutter echo data in the fourth fuzzy area, and constructing a fourth blocking vector to obtain data after the clutter blocking in the fourth fuzzy area.
In one embodiment of the present invention, the first occlusion vector is:
Figure 573228DEST_PATH_IMAGE017
the second occlusion vector is:
Figure 846078DEST_PATH_IMAGE018
the third occlusion vector is:
Figure 512682DEST_PATH_IMAGE019
the fourth occlusion vector is:
Figure 523320DEST_PATH_IMAGE020
the data after the clutter is blocked in the first fuzzy region is data;
Figure 958981DEST_PATH_IMAGE021
the data after the clutter blocking of the second fuzzy region is data;
Figure 668311DEST_PATH_IMAGE022
the data after the third fuzzy region blocks the clutter is;
Figure 556633DEST_PATH_IMAGE023
the data after the fourth fuzzy region blocks the clutter is;
Figure 111242DEST_PATH_IMAGE024
wherein,
Figure 666988DEST_PATH_IMAGE025
the dimensions of the matrix are represented by,C1C2C3C4the kronecker products of the received space steering vector and the doppler steering vector respectively corresponding to different fuzzy areas,
Figure 281640DEST_PATH_IMAGE026
m denotes the mth transmit waveform, and p denotes the number of blur areas.
In an embodiment of the present invention, after the step 6, the method further includes:
and processing clutter data of the first fuzzy area by adopting a Doppler frequency shift compensation method, and then processing the clutter data compensated by the Doppler frequency shift compensation method in the first fuzzy area and the clutter data of the second fuzzy area, the third fuzzy area and the fourth fuzzy area by adopting a STAP algorithm so as to eliminate the clutter data of each fuzzy area.
The invention has the beneficial effects that:
the invention has the advantages of only one transmitting array element, low hardware cost, low calculation complexity, easy realization, low requirement on system parameters, no need of considering the fuzzy condition of transmitting spatial frequency, good clutter separation and good inhibition effect.
The present invention will be described in further detail with reference to the accompanying drawings and examples.
Drawings
Fig. 1 is a schematic flowchart of a short-range clutter suppression method based on a multi-frequency diversity radar according to an embodiment of the present invention;
fig. 2 is a schematic flowchart of another short-range clutter suppression method based on a multi-frequency diversity radar according to an embodiment of the present invention;
FIG. 3 is a schematic diagram illustrating the relationship between the transmission spatial frequency and the distance after the first blurring region compensation is performed according to an embodiment of the present invention;
FIG. 4a is a schematic diagram of the power spectrum distribution of clutter in the first fuzzy region after mixed clutter data separation according to an embodiment of the present invention;
FIG. 4b is a schematic diagram of the power spectrum distribution of clutter in the second fuzzy region after mixed clutter data separation according to an embodiment of the present invention;
FIG. 4c is a schematic diagram of the power spectrum distribution of clutter in the third fuzzy region after mixed clutter data separation according to an embodiment of the present invention;
FIG. 4d is a schematic diagram of the power spectrum distribution of clutter in the fourth fuzzy region after mixed clutter data separation according to an embodiment of the present invention;
FIG. 5a is a diagram illustrating the relationship between the improvement factor and the normalized Doppler frequency after STAP is performed on the data of the first blurred region after separation according to an embodiment of the present invention;
FIG. 5b is a diagram illustrating the relationship between the improvement factor and the normalized Doppler frequency after STAP is performed on the data of the second blurred region after separation according to the embodiment of the present invention;
FIG. 5c is a diagram illustrating the relationship between the improvement factor and the normalized Doppler frequency after STAP is performed on the data of the third blurred region directly after separation according to an embodiment of the present invention;
FIG. 5d is a diagram illustrating the relationship between the improvement factor and the normalized Doppler frequency after STAP is directly performed on the data of the fourth blurred region after separation according to the embodiment of the present invention;
fig. 6 is a schematic diagram of a clutter power spectrum distribution of the first blurred region after doppler shift compensation is performed on the first blurred region according to an embodiment of the present invention;
fig. 7 is a comparison graph of the relationship between the normalized doppler frequency and the improvement factor of the STAP performed after the first blurred region is compensated, and the STAP performed directly in the past.
Detailed Description
The present invention will be described in further detail with reference to specific examples, but the embodiments of the present invention are not limited thereto.
Example one
Referring to fig. 1 and fig. 2, fig. 1 is a schematic flow chart of a short-range clutter suppression method based on a multi-frequency system radar according to an embodiment of the present invention, and fig. 2 is a schematic flow chart of another short-range clutter suppression method based on a multi-frequency system radar according to an embodiment of the present invention. The embodiment of the invention provides a short-range clutter suppression method based on a multi-frequency split radar, which comprises the following steps:
and step 1, obtaining clutter echo data based on a radar system of superposition step frequency.
Specifically, setting a transmitting antenna to transmit M orthogonal signals, each of which transmits K pulses, N receiving antennas in total and P fuzzy areas in total, and then receiving clutter data by the ith distance unit:
Figure 657258DEST_PATH_IMAGE027
wherein,
Figure 749979DEST_PATH_IMAGE028
the transmitting space guide vector, the receiving space guide vector and the Doppler guide vector of the clutter echo data are respectively represented, and Nc represents the equivalent number of clutter sources in each range unit.
And 2, analyzing the two-dimensional distribution of the clutter echo data on the receiving space frequency and the Doppler frequency according to the receiving space frequency and the Doppler frequency in the clutter echo data.
Specifically, the property of the clutter may be determined according to the two-dimensional distribution, which is used as a basis for the property change of the clutter before and after the back compensation, and the two-dimensional power spectrum distribution of the clutter echo data is:
Figure 160232DEST_PATH_IMAGE029
Figure 70419DEST_PATH_IMAGE030
wherein,
Figure 667753DEST_PATH_IMAGE031
representing the received spatial frequency and the normalized Doppler frequency, respectively
Figure 829744DEST_PATH_IMAGE032
The space-time two-dimensional steering vector of (a),Ra covariance matrix representing clutter echo data, L represents a total number of training samples, L represents an ith training sample,Hrepresenting a conjugate transpose, and the training samples represent clutter echo data at the ith range bin.
And 3, obtaining the relation between the transmitting space frequency and the distance according to the transmitting space frequency of the clutter echo data so as to find out the relation between the transmitting space frequency and the fuzzy area.
In this embodiment, the mth transmission signal of the transmission array element can be expressed as:
Figure 563345DEST_PATH_IMAGE033
wherein,
Figure 785379DEST_PATH_IMAGE034
which represents the envelope of the transmitted signal,Ewhich represents the total power of the transmitted signal,Twhich represents the duration of the transmitted pulse or pulses,ta time variable is represented by a time variable,f m representing the carrier frequency of the mth transmitted signal. Assuming that the envelope of the transmitted signal satisfies the orthogonality condition, there are:
Figure 870010DEST_PATH_IMAGE035
wherein,
Figure 835692DEST_PATH_IMAGE036
which represents an arbitrary time delay, is,
Figure 814012DEST_PATH_IMAGE037
it is indicated that the conjugate operation is performed,
Figure 669929DEST_PATH_IMAGE038
representing the envelope of the signal.
Considering a far-field point source, the time delay for the nth unit to receive a signal transmitted via the far-field point source can be expressed as:
Figure 507435DEST_PATH_IMAGE039
wherein,
Figure 745650DEST_PATH_IMAGE040
which represents a common time delay, is,Rrepresenting the distance between the receiving array element to the far field point,cthe speed of light is indicated and is,
Figure 985001DEST_PATH_IMAGE041
representing the time delay difference between the cells,
Figure 17679DEST_PATH_IMAGE042
which is indicative of the speed of the platform,
Figure 201536DEST_PATH_IMAGE043
representing the angle of the far field point.
The kth pulse received by the nth unit can be expressed as:
Figure 509021DEST_PATH_IMAGE044
wherein,
Figure 71720DEST_PATH_IMAGE045
representing the complex scattering coefficient of the point source. Taking into account the narrowband assumption, i.e.
Figure 806458DEST_PATH_IMAGE046
Then the received signal is at
Figure 87398DEST_PATH_IMAGE047
After matched filtering, it can be decomposed into:
Figure 198573DEST_PATH_IMAGE049
wherein,
Figure 146938DEST_PATH_IMAGE050
Figure 786997DEST_PATH_IMAGE051
which is indicative of the carrier frequency,dthe space between the array elements is shown,
Figure 820813DEST_PATH_IMAGE052
representing the doppler frequency of the point source,
Figure 329154DEST_PATH_IMAGE053
which represents the wavelength of the light emitted by the light source,
Figure 272971DEST_PATH_IMAGE054
representing the point source velocity. The receive snapshots in frequency diversity array MIMO can therefore be expressed in vector form as:
Figure 349511DEST_PATH_IMAGE055
wherein,
Figure 729677DEST_PATH_IMAGE056
the superscript T is the transpose operator,
Figure 911216DEST_PATH_IMAGE057
is the product of Kronecker.
Figure 568594DEST_PATH_IMAGE058
Figure 816035DEST_PATH_IMAGE059
And
Figure 558863DEST_PATH_IMAGE060
respectively representing transmit, receive steering vectors and time steering vectors, expressed as follows:
Figure 549953DEST_PATH_IMAGE061
Figure 61837DEST_PATH_IMAGE062
Figure 480180DEST_PATH_IMAGE063
note here that:
Figure 975884DEST_PATH_IMAGE064
in order to transmit the spatial frequencies,
Figure 770664DEST_PATH_IMAGE065
Figure 402634DEST_PATH_IMAGE066
the increment of the frequency is represented by,cthe speed of light is indicated and is,Rthe distance is represented as a function of time,
Figure 257457DEST_PATH_IMAGE067
indicating a dominant distance, P indicating the number of blurred regions,
Figure 506036DEST_PATH_IMAGE068
maximum unambiguous distance;
Figure 838928DEST_PATH_IMAGE069
to accept spatial frequencies;
Figure 466350DEST_PATH_IMAGE070
is the normalized doppler frequency.
And 4, constructing a principal value distance compensation vector based on the relation between the transmitting spatial frequency and the fuzzy region to obtain compensated clutter echo data, wherein fig. 3 is a schematic diagram of the relation between the transmitting spatial frequency and the distance after the first fuzzy region is compensated.
In this embodiment, the principal value distance compensation vector is:
Figure 757654DEST_PATH_IMAGE071
where M is the number of transmit waveforms.
Preferably, the blur area includes a first blur area, a second blur area, a third blur area, and a fourth blur area. Then:
the principal value distance compensation vector of the first fuzzy area is as follows:
Figure 979949DEST_PATH_IMAGE072
the principal value distance compensation vector of the second fuzzy area is as follows:
Figure 382111DEST_PATH_IMAGE073
the principal value distance compensation vector of the third fuzzy area is as follows:
Figure 988673DEST_PATH_IMAGE074
the principal value distance compensation vector of the fourth fuzzy area is as follows:
Figure 326244DEST_PATH_IMAGE075
in a specific embodiment, step 4 may specifically include:
and 4.1, compensating the dominant value distance of the first fuzzy area through the corresponding dominant value distance compensation vector to obtain the clutter echo data compensated in the first fuzzy area.
Specifically, the dominant value distance of the first fuzzy region is compensated through the corresponding dominant value distance compensation vector, clutter echo data irrelevant to the dominant value distance are obtained, the clutter echo data at the moment are only relevant to the number of fuzzy regions, and the transmitting space frequency of the clutter of the first fuzzy region is a constant irrelevant to the number of channels.
And 4.2, compensating the dominant value distance of the second fuzzy area through the corresponding dominant value distance compensation vector to obtain the clutter echo data compensated in the second fuzzy area.
Specifically, the dominant value distance of the second fuzzy area is compensated through the corresponding dominant value distance compensation vector, and the transmitting spatial frequency of the clutter of the second fuzzy area is a constant independent of the number of channels.
And 4.3, compensating the dominant value distance of the third fuzzy area through the corresponding dominant value distance compensation vector to obtain the clutter echo data compensated in the third fuzzy area.
And 4.4, compensating the dominant value distance of the fourth fuzzy area through the corresponding dominant value distance compensation vector to obtain the clutter echo data compensated in the fourth fuzzy area.
In one embodiment, the clutter echo data after the dominant distance compensation is:
Figure 18257DEST_PATH_IMAGE076
wherein,
Figure 958531DEST_PATH_IMAGE077
representing an identity matrix of dimension N,
Figure 419599DEST_PATH_IMAGE078
representing an identity matrix of dimension K.
And 5, constructing a blocking vector according to the compensated clutter echo data to obtain the blocked data of each fuzzy area.
In one embodiment, step 5 comprises:
and 5.1, rearranging the compensated clutter echo data in the first fuzzy area, and constructing a first blocking vector to obtain data after the clutter echo data are blocked in the first fuzzy area.
Specifically, one of the range units is selected, the clutter echo data is rearranged, a first blocking vector is constructed, and the clutter data of the first fuzzy area is eliminated from the clutter echo data.
After the principal distance compensation of the first fuzzy region, the clutter echo data can be represented again after being rearranged as:
Figure 52706DEST_PATH_IMAGE079
wherein,
Figure 232015DEST_PATH_IMAGE080
Figure 851346DEST_PATH_IMAGE081
Figure 432500DEST_PATH_IMAGE082
Figure 705350DEST_PATH_IMAGE083
and respectively representing the emission spatial frequency of each of the four fuzzy areas after the first fuzzy area compensation.
The first occlusion vector is:
Figure 637533DEST_PATH_IMAGE084
the data after the clutter is blocked in the first fuzzy region is as follows;
Figure 919610DEST_PATH_IMAGE085
wherein,
Figure 614991DEST_PATH_IMAGE086
representing data after the first blurred region blocks clutter.
And 5.2, rearranging the compensated clutter echo data in the second fuzzy area, and constructing a second blocking vector to obtain data after the clutter is blocked in the second fuzzy area.
Specifically, the same distance unit is selected, clutter echo data are rearranged, a second blocking vector is constructed again, and clutter data in a second fuzzy area are eliminated from the clutter echo data.
Wherein the second occlusion vector is:
Figure 58741DEST_PATH_IMAGE087
the data after the clutter is blocked in the second fuzzy region is as follows;
Figure 353588DEST_PATH_IMAGE088
wherein,
Figure 173776DEST_PATH_IMAGE089
representing data after the second blurred region blocks clutter.
And 5.3, rearranging the compensated clutter echo data in the third fuzzy area, and constructing a third blocking vector to obtain data after the third fuzzy area blocks the clutter.
Wherein the third occlusion vector is:
Figure 729522DEST_PATH_IMAGE090
the data after the third fuzzy region blocks the clutter is;
Figure 609754DEST_PATH_IMAGE091
wherein,
Figure 985371DEST_PATH_IMAGE092
representing data after the third blurred region blocks clutter.
And 5.4, rearranging the compensated clutter echo data in the fourth fuzzy area, and constructing a fourth blocking vector to obtain data after the clutter is blocked in the fourth fuzzy area.
Wherein the fourth occlusion vector is:
Figure 609251DEST_PATH_IMAGE093
the data after the fourth fuzzy region blocks the clutter is;
Figure 753924DEST_PATH_IMAGE094
wherein,
Figure 539478DEST_PATH_IMAGE095
representing data after the fourth blurred region blocks clutter.
Wherein,
Figure 402391DEST_PATH_IMAGE096
the dimensions of the matrix are represented by,C1C2C3C4the kronecker products of the received space steering vector and the doppler steering vector respectively corresponding to different fuzzy areas,
Figure 564383DEST_PATH_IMAGE097
m denotes the mth transmit waveform, and p denotes the number of blur areas.
And 6, establishing a simultaneous equation set for the expressions of the blocked data of each fuzzy area to obtain the clutter data of each fuzzy area.
And 7, traversing all the distance units to obtain clutter data of all the fuzzy areas, processing the clutter data of the first fuzzy area by adopting a Doppler frequency shift compensation method, and then processing the clutter data compensated by the Doppler frequency shift compensation method in the first fuzzy area and the clutter data of the second fuzzy area, the third fuzzy area and the fourth fuzzy area by adopting an STAP algorithm to eliminate the clutter data of each fuzzy area.
Next, the effect of the method for suppressing short-range clutter by using a multi-frequency division radar according to the present embodiment is further described through a simulation experiment.
Figure 94721DEST_PATH_IMAGE098
The power spectrum distribution of the clutter after the clutter is separated by the method of the invention is shown in fig. 4a-4d, as can be seen from fig. 4a-4d, the distance dependency of the short-range clutter is more serious, which causes the power spectrum distribution to spread, while the distance dependency of the clutter of the other fuzzy areas is smaller, and the clutter distributions of different distances are almost the same, which is also the reason that the notch of the improvement factor of the first fuzzy area is wider and the clutter suppression performance is poorer in fig. 5a-5d, while the notch of the improvement factor of the other fuzzy areas by directly using the STAP is narrower, so that the clutter of the areas can be directly suppressed by using the STAP. Fig. 6 shows the clutter power spectrum distribution after applying the doppler shift compensation algorithm to the clutter in the first blurred region, and it can be clearly seen that the distribution of the mainlobe clutter is more concentrated. Fig. 7 is an improved factor graph obtained by using the STAP algorithm after compensation, the notch is already obviously narrowed, and the suppression effect is improved.
The invention has the advantages of only one transmitting array element, low hardware cost, low calculation complexity, easy realization, low requirement on system parameters, no need of considering the fuzzy condition of transmitting spatial frequency, good clutter separation and good inhibition effect.
In the description of the present invention, the terms "first" and "second" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implying any number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of the present invention, "a plurality" means two or more unless specifically defined otherwise.
In the description herein, references to the description of the term "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., mean that a particular feature, structure, material, or characteristic data point described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, various embodiments or examples described in this specification can be combined and combined by those skilled in the art.
The foregoing is a more detailed description of the invention in connection with specific preferred embodiments and it is not intended that the invention be limited to these specific details. For those skilled in the art to which the invention pertains, several simple deductions or substitutions can be made without departing from the spirit of the invention, and all shall be considered as belonging to the protection scope of the invention.

Claims (10)

1. A short-range clutter suppression method based on a multi-frequency division radar is characterized by comprising the following steps:
step 1, obtaining clutter echo data based on a radar system with superposed step frequency;
step 2, analyzing two-dimensional distribution of the clutter echo data on the receiving space frequency and the Doppler frequency according to the receiving space frequency and the Doppler frequency in the clutter echo data;
step 3, obtaining the relation between the transmitting space frequency and the distance according to the transmitting space frequency of the clutter echo data so as to find out the relation between the transmitting space frequency and the fuzzy area;
step 4, constructing a principal value distance compensation vector based on the relation between the emission space frequency and the fuzzy area to obtain compensated clutter echo data;
step 5, constructing a blocking vector according to the compensated clutter echo data to obtain the blocked data of each fuzzy area;
and 6, establishing a simultaneous equation set for the expressions of the blocked data of the fuzzy areas to obtain the clutter data of each fuzzy area.
2. The short-range clutter suppression method based on multi-band radar of claim 1, wherein the transmit antenna is configured to transmit M orthogonal signals, each of the orthogonal signals transmits K pulses, there are N receive antennas in total, and there are P fuzzy areas in total, then the clutter data received by the ith distance unit is:
Figure 320919DEST_PATH_IMAGE001
wherein,
Figure 175743DEST_PATH_IMAGE002
the transmitting space guide vector, the receiving space guide vector and the Doppler guide vector of the clutter echo data are respectively represented, and Nc represents the equivalent number of clutter sources in each range unit.
3. The short-range clutter suppression method based on multi-frequency division radar according to claim 1, wherein the two-dimensional distribution is:
Figure 830846DEST_PATH_IMAGE003
Figure 163738DEST_PATH_IMAGE004
wherein,
Figure 853477DEST_PATH_IMAGE005
representing the received spatial frequency and the normalized Doppler frequency, respectively
Figure 613623DEST_PATH_IMAGE006
The space-time two-dimensional steering vector of (a),Ra covariance matrix representing clutter echo data, L represents a total number of training samples, L represents an ith training sample,Hrepresenting a conjugate transpose, and the training samples represent clutter echo data at the ith range bin.
4. The short-range clutter suppression method based on multi-band radar according to claim 1, wherein the relationship between the transmit spatial frequency and the distance is:
Figure 753092DEST_PATH_IMAGE007
Figure 889676DEST_PATH_IMAGE008
wherein,
Figure 496237DEST_PATH_IMAGE009
which represents the spatial frequencies of the reception of the signal,
Figure 427284DEST_PATH_IMAGE010
the increment of the frequency is represented by,cthe speed of light is indicated and is,Rthe distance is represented as a function of time,
Figure 56980DEST_PATH_IMAGE011
indicating a dominant distance, P indicating the number of blurred regions,
Figure 731675DEST_PATH_IMAGE012
maximum unambiguous distance.
5. The short-range clutter suppression method based on multi-frequency division radar according to claim 4, wherein the principal value distance compensation vector is:
Figure 189814DEST_PATH_IMAGE013
where M is the number of transmit waveforms.
6. The short-range clutter suppression method based on multi-frequency division radar according to claim 5, wherein the clutter echo data after compensation is:
Figure 495024DEST_PATH_IMAGE014
wherein,I N representing an identity matrix of dimension N,I K representing an identity matrix of dimension K.
7. The short-range clutter suppression method based on multi-frequency diversity radar according to claim 6, wherein the fuzzy areas comprise a first fuzzy area, a second fuzzy area, a third fuzzy area and a fourth fuzzy area, wherein said step 4 comprises:
step 4.1, compensating the dominant value distance of the first fuzzy area through the corresponding dominant value distance compensation vector to obtain clutter echo data after compensation in the first fuzzy area;
step 4.2, compensating the dominant value distance of the second fuzzy area through the corresponding dominant value distance compensation vector to obtain clutter echo data after compensation in the second fuzzy area;
4.3, compensating the dominant value distance of the third fuzzy area through the corresponding dominant value distance compensation vector to obtain clutter echo data after compensation in the third fuzzy area;
and 4.4, compensating the principal value distance of the fourth fuzzy area through the corresponding principal value distance compensation vector to obtain the clutter echo data after compensation in the fourth fuzzy area.
8. The short-range clutter suppression method based on multi-frequency division radar according to claim 6, wherein said step 5 comprises:
step 5.1, rearranging the compensated clutter echo data in the first fuzzy area, and constructing a first blocking vector to obtain data after the clutter echo data are blocked in the first fuzzy area;
step 5.2, rearranging the compensated clutter echo data in the second fuzzy area, and constructing a second blocking vector to obtain data after the clutter blocking in the second fuzzy area;
step 5.3, rearranging the compensated clutter echo data in the third fuzzy area, and constructing a third blocking vector to obtain data after the third fuzzy area blocks the clutter;
and 5.4, rearranging the compensated clutter echo data in the fourth fuzzy area, and constructing a fourth blocking vector to obtain data after the clutter blocking in the fourth fuzzy area.
9. The short-range clutter suppression method based on multi-band radar according to claim 8, wherein the first blocking vector is:
Figure 674333DEST_PATH_IMAGE015
the second occlusion vector is:
Figure 887139DEST_PATH_IMAGE016
the third occlusion vector is:
Figure 405976DEST_PATH_IMAGE017
the fourth occlusion vector is:
Figure 678826DEST_PATH_IMAGE018
the data after the clutter is blocked in the first fuzzy region is data;
Figure 342501DEST_PATH_IMAGE020
the data after the clutter blocking of the second fuzzy region is data;
Figure 358999DEST_PATH_IMAGE022
the data after the third fuzzy region blocks the clutter is;
Figure 732342DEST_PATH_IMAGE024
the data after the fourth fuzzy region blocks the clutter is;
Figure 441672DEST_PATH_IMAGE026
wherein,
Figure 329994DEST_PATH_IMAGE027
the dimensions of the matrix are represented by,C1C2C3C4kronecker product of received space steering vector and doppler steering vector corresponding to different fuzzy areas
Figure 87865DEST_PATH_IMAGE028
M denotes the mth transmit waveform, and p denotes the number of blur areas.
10. The short-range clutter suppression method based on multi-frequency division radar according to claim 7, further comprising, after said step 6:
and processing clutter data of the first fuzzy area by adopting a Doppler frequency shift compensation method, and then processing the clutter data compensated by the Doppler frequency shift compensation method in the first fuzzy area and the clutter data of the second fuzzy area, the third fuzzy area and the fourth fuzzy area by adopting a STAP algorithm so as to eliminate the clutter data of each fuzzy area.
CN202210184291.9A 2022-02-28 2022-02-28 Short-range clutter suppression method based on multi-frequency split radar Active CN114265018B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210184291.9A CN114265018B (en) 2022-02-28 2022-02-28 Short-range clutter suppression method based on multi-frequency split radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210184291.9A CN114265018B (en) 2022-02-28 2022-02-28 Short-range clutter suppression method based on multi-frequency split radar

Publications (2)

Publication Number Publication Date
CN114265018A true CN114265018A (en) 2022-04-01
CN114265018B CN114265018B (en) 2022-05-31

Family

ID=80833700

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210184291.9A Active CN114265018B (en) 2022-02-28 2022-02-28 Short-range clutter suppression method based on multi-frequency split radar

Country Status (1)

Country Link
CN (1) CN114265018B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104635214A (en) * 2015-02-12 2015-05-20 西安电子科技大学 Airborne forward-looking frequency diversity array radar distance fuzzy clutter suppression method
CN106569212A (en) * 2016-11-09 2017-04-19 西安空间无线电技术研究所 Multichannel SAR-GMTI range ambiguity clutter suppression method
CN108020817A (en) * 2017-09-28 2018-05-11 西安电子科技大学 Air-borne Forward-looking battle array radar clutter suppression method based on registration
CN109814070A (en) * 2019-01-31 2019-05-28 西安电子科技大学 Range ambiguity clutter suppression method based on false impulse
CN113253223A (en) * 2021-03-30 2021-08-13 北京理工大学 Target detection method for non-stationary clutter suppression based on step frequency signal
CN113376599A (en) * 2021-01-19 2021-09-10 西安电子科技大学 FDA distance fuzzy clutter suppression method based on mainlobe correction

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104635214A (en) * 2015-02-12 2015-05-20 西安电子科技大学 Airborne forward-looking frequency diversity array radar distance fuzzy clutter suppression method
CN106569212A (en) * 2016-11-09 2017-04-19 西安空间无线电技术研究所 Multichannel SAR-GMTI range ambiguity clutter suppression method
CN108020817A (en) * 2017-09-28 2018-05-11 西安电子科技大学 Air-borne Forward-looking battle array radar clutter suppression method based on registration
CN109814070A (en) * 2019-01-31 2019-05-28 西安电子科技大学 Range ambiguity clutter suppression method based on false impulse
CN113376599A (en) * 2021-01-19 2021-09-10 西安电子科技大学 FDA distance fuzzy clutter suppression method based on mainlobe correction
CN113253223A (en) * 2021-03-30 2021-08-13 北京理工大学 Target detection method for non-stationary clutter suppression based on step frequency signal

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
兰岚等: "波形分集阵列雷达抗干扰进展", 《***工程与电子技术》 *

Also Published As

Publication number Publication date
CN114265018B (en) 2022-05-31

Similar Documents

Publication Publication Date Title
EP3589970B1 (en) Method and system for obtaining an adaptive angle-doppler ambiguity function in mimo radars
CN110412568B (en) Distance fuzzy clutter suppression method based on extended azimuth phase coding
CN111880171B (en) Pulse segment coding method for eliminating radar target blind speed
CN104297734B (en) Deceiving interference suppressing method based on the MIMO radar of frequency diversity array
CA2901610C (en) Surface wave radar
EP0815470B1 (en) Adaptive filtering of matched-filter data
US20070247353A1 (en) Method and Apparatus for Performing Bistatic Radar Functions
US5748143A (en) Adaptive post-doppler sequential beam processor
Bergin et al. Radar waveform optimization for colored noise mitigation
CN113253223B (en) Target detection method for non-stationary clutter suppression based on step frequency signal
Blunt et al. Multistatic adaptive pulse compression
Rabideau Clutter and jammer multipath cancellation in airborne adaptive radar
CN103969641B (en) A kind of beam transmitting three-D imaging method
CN113376599B (en) FDA distance fuzzy clutter suppression method based on mainlobe correction
CN114895261A (en) Clutter suppression method based on multi-frequency sub-pulse coding array
EP2281325B1 (en) A process for minimising jammer noise in receiver systems
CN114265018B (en) Short-range clutter suppression method based on multi-frequency split radar
Vijaykumar Mahamuni Space-time adaptive processing techniques (STAP) for mitigation of jammer interference and clutter suppression in airborne radar systems: A MATLAB implementation-based study
CN114609595A (en) Frequency division orthogonal MIMO radar signal processing method
CN112834991B (en) MIMO radar slow target detection method based on time domain frequency diversity
CN111505600B (en) STPC-based FDA-MIMO radar signal processing method, device and medium
Gerlach et al. Combined multistatic adaptive pulse compression and adaptive beamforming for shared-spectrum radar
CN113253222A (en) Airborne FDA-MIMO bistatic radar distance fuzzy clutter suppression and dimension reduction search method
Heng et al. Modified JDL with Doppler compensation for airborne bistatic radar
CN114966571B (en) Noise convolution interference suppression method based on frequency diversity MIMO radar

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230818

Address after: 710065 Room 71203, Unit 7, Building 1, Tianyuan Apartment, East Section of Art Street, New Industrial Park, High tech Zone, Xi'an City, Shaanxi Province

Patentee after: Xi'an Hesai Electronic Technology Co.,Ltd.

Address before: 710071 Taibai South Road, Yanta District, Xi'an, Shaanxi Province, No. 2

Patentee before: XIDIAN University