CN113219433A - Knowledge-aided SR-STAP method and storage medium - Google Patents
Knowledge-aided SR-STAP method and storage medium Download PDFInfo
- Publication number
- CN113219433A CN113219433A CN202110577866.9A CN202110577866A CN113219433A CN 113219433 A CN113219433 A CN 113219433A CN 202110577866 A CN202110577866 A CN 202110577866A CN 113219433 A CN113219433 A CN 113219433A
- Authority
- CN
- China
- Prior art keywords
- doppler frequency
- space
- frequency
- clutter
- knowledge
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses a knowledge-aided SR-STAP method and a storage medium, wherein the method comprises the following steps: acquiring sample data of a distance unit to be detected through an airborne radar; uniformly dispersing the spatial frequency of the clutter into a grid point set; acquiring a discrete Doppler frequency grid point set of the clutter; adjusting the Doppler frequency based on a preset condition to obtain a new Doppler frequency group; adjusting the spatial frequency to obtain a new discretization spatial frequency group; constructing a space-time guide vector dictionary according to the new Doppler frequency group and the space-time guide vector corresponding to the new discretization space frequency group; carrying out sparse decomposition on the sample data by utilizing the space-time guiding vector dictionary to obtain a sparse recovery vector; and performing clutter covariance matrix estimation by using the sparse recovery vector. The invention jointly adjusts the space frequency and the Doppler frequency, constructs an ultra-complete space-time guiding vector dictionary, effectively solves the problem of grid mismatch and improves the performance of the STAP method.
Description
Technical Field
The invention relates to the technical field of airborne radars, in particular to a knowledge-aided SR-STAP method and a storage medium.
Background
The airborne radar has excellent moving target detection capability and has great application value in the fields of military and civil use. However, the airborne radar is in a downward-looking working state, and the clutter intensity in the received signal is high. In addition, the movement of the carrier leads the distribution range of the noise waves to be widened and the intensity to be increased. Meanwhile, due to the influence of factors such as terrain variation, electromagnetic interference, system errors and the like, clutter presents obvious non-uniform distribution, and the target detection capability of the radar is seriously influenced. In order to solve the above problems, effectively suppress clutter and improve the target detection capability of an airborne radar, a scholars has proposed a Space-time Adaptive Processing (STAP) method. The method utilizes the coupling of the clutter in space and time to carry out self-adaptive filtering in a space-time two-dimensional domain, thereby accurately detecting the moving target. In order to ensure better clutter suppression and moving target detection effects, the conventional STAP algorithm generally needs to acquire enough training samples from a uniform clutter environment. However, in practical scenarios, especially in non-homogeneous clutter environments, it is difficult to meet this requirement. Thus, a drastic performance degradation of the conventional STAP method results.
With the rise of the compressive sensing theory, some scholars propose a sparse recovery method for space-time adaptive filtering, and gradually form a sparse power spectrum recovery-based STAP method (SR-STAP). In the method, a redundant base dictionary formed by space-time guide vectors is used for carrying out sparse representation on sample data, elements with larger coefficients are selected to reconstruct the sample data, high-resolution power spectrum estimation of clutter can be formed, and better clutter suppression performance is realized. The dictionary used by sparse representation is generally composed of space-time two-dimensional guide vectors corresponding to uniformly discretized spatial frequency and Doppler frequency, and the precondition is that guide vectors of clutter in sample data are accurately matched with part of guide vectors in the dictionary, namely, the clutter can accurately fall on discretized space-time grid points. In practical application, however, due to factors such as an array working mode and system parameters, a difference exists between a clutter guide vector and a dictionary guide vector, that is, a part of clutter components are no longer located on a discretization space-time plane grid point, so-called off-grid effect occurs, and thus sparse recovery errors are significantly increased and the performance of the STAP filter is seriously reduced. The conventional SR-STAP method uses a fixed dictionary formed by space-time plane uniform discretization, so that clutter sparse recovery has large errors, and clutter suppression and moving target detection are not facilitated.
In order to improve the performance of the SR-STAP method, researchers at home and abroad explore a method for solving the off-grid effect to increase the sparse recovery accuracy. The knowledge-assisted SR-STAP method becomes a research hotspot of scholars at home and abroad in recent years, and scientific research results show that the knowledge-assisted SR-STAP method can better solve the problems of insufficient training samples and performance reduction of the STAP caused by an off-network effect, and relieve the pressure caused by the operation amount and data storage. Researchers at the university of Wuhan Engineers in 2018 propose an algorithm (KASR-STAP) for alleviating the off-net effect in the SR-STAP method by adjusting a dictionary by using the prior knowledge of the clutter ridge. The algorithm utilizes the prior knowledge of radar system parameters, inertial navigation system information and the like to calculate the clutter ridge distribution, and disperses grids on an angle-Doppler plane, so that the off-grid effect can be well relieved under different array scenes. However, the KASR-STAP method has the problems that when the prior knowledge is used for discretizing grids on an angle-doppler plane, the grid position of doppler frequency is adjusted only according to clutter distribution, the grid position of spatial frequency is not correspondingly adjusted, improvement on off-grid effect (grid mismatch problem) is incomplete, and improvement of the SR-STAP method performance is not ideal. The existing dictionary correction method based on knowledge assistance does not fully consider the coupling of Doppler frequency and spatial frequency, the dictionary correction is not complete, and the performance improvement degree of STAP is not high enough.
Disclosure of Invention
In view of the above technical problems in the prior art, embodiments of the present invention provide a knowledge-aided SR-STAP method, which can solve the technical problems that the existing SR-STAP method is not favorable for clutter suppression and moving target detection, and the existing knowledge-aided dictionary correction method is not comprehensive in improving off-grid effect and the improvement degree of STAP performance is not high enough.
In order to solve the technical problem, the embodiment of the invention adopts the following technical scheme:
a knowledge-aided SR-STAP method comprises the following steps:
acquiring sample data of a distance unit to be detected through an airborne radar, wherein the airborne radar is in front side view work;
uniformly dispersing the spatial frequency of the clutter into a grid point set;
acquiring a discrete Doppler frequency grid point set of the clutter;
adjusting the Doppler frequency based on a preset condition to obtain a new Doppler frequency group;
adjusting the spatial frequency based on the new Doppler frequency group to obtain a new discretization spatial frequency group;
constructing a space-time guide vector dictionary according to the new Doppler frequency group and the space-time guide vector corresponding to the new discretization space frequency group;
carrying out sparse decomposition on the sample data by utilizing the space-time guiding vector dictionary to obtain a sparse recovery vector;
and performing clutter covariance matrix estimation by using the sparse recovery vector.
Further, uniformly discretizing the spatial frequency of the clutter into a set of grid points, comprising:
according to the formula fsUniformly dispersing the normalized spatial frequency into NsSet of grid points { fs,iIn which N iss＝ρsN，ρsAnd dividing coefficients for space frequency grids, wherein N is the number of array elements of the airborne radar.
Further, obtaining a set of discrete doppler frequency grid points for the clutter comprises:
according to the formula
Obtaining a corresponding discrete Doppler frequency grid point set { fd,iIn which θαIs the angle theta between the array axis and the direction of flight of the aircraftpDirecting the transmitted beam of the radar antennaThe angle of the aircraft wing-to-wing,
angle of elevation for directing the beam emitted by the radar antenna, beta 4V/lambda frIs a folding factor, V is the aircraft flight speed, lambda is the radar operating wavelength, frIs the pulse repetition frequency.
Further, adjusting the doppler frequency based on the preset condition to obtain a new doppler frequency group, including:
refining the Doppler frequency interval based on a first preset condition;
combining adjacent Doppler frequency intervals based on a second preset condition;
dispersing Doppler frequency meeting a third preset condition by using space frequency interval;
and combining the Doppler frequencies adjusted based on the first preset condition, the second preset condition and the third preset condition to obtain a new Doppler frequency group.
Further, the refining the doppler frequency interval based on the first preset condition includes:
according to the formula
Refining the Doppler frequency interval to make the original Doppler frequency interval deltad,iNot greater than spatial frequency interval ΔsWherein, iss＝|fs,i+1-fs|，Δd,i＝|fd,i+1-fd,i|,i＝1,...,Ns-1。
Further, the merging adjacent doppler frequency intervals based on a second preset condition includes:
according to the formula
Combining adjacent Doppler frequency intervals to ensure that the original Doppler frequency interval is greater than a deltasWherein a is an adjusting parameter, and a is more than 0 and less than 1.
Further, the discretizing, by using the spatial frequency interval, the doppler frequency meeting the third preset condition includes:
using spatial frequency separation ΔsDiscrete coincidence a Δs≤Δd,i≤ΔsThe doppler frequency of the condition.
Further, adjusting the spatial frequency based on the new doppler frequency set to obtain a new discretized spatial frequency set, including:
according to the formula
Adjusting the spatial frequency to obtain a new set of discretized spatial frequency
Further, the method further comprises:
from the estimated clutter covariance matrix
Using formulas
Calculating adaptive filtering weights, wherein StIs a space-time steering vector, mu is a non-zero constant;
and filtering the sample data of the distance unit to be detected according to the calculated adaptive filtering weight.
Embodiments of the present invention further provide a computer-readable storage medium, on which computer-executable instructions are stored, and when executed by a processor, the computer-executable instructions implement the steps of the knowledge-aided SR-STAP method described above.
Compared with the prior art, the knowledge-aided SR-STAP method provided by the embodiment of the invention combines knowledge assistance and the SR-STAP method, jointly adjusts the space frequency and the Doppler frequency interval by means of the priori knowledge of the airborne radar array, the space-time coupling of clutter distribution and other priori information, constructs an ultra-complete space-time guide vector dictionary, comprehensively corrects the dictionary, effectively solves the problem of grid mismatch, improves the sparse recovery precision of clutter, and improves the performance of the STAP method. In addition, the invention reduces the number of samples required by the algorithm, and is particularly suitable for moving target detection in a clutter non-uniform environment.
Drawings
FIG. 1 is a flow diagram of a knowledge-aided SR-STAP method according to an embodiment of the present invention;
FIG. 2(a) is a diagram of a clutter power spectrum simulation result obtained by using a knowledge-aided SR-STAP method according to an embodiment of the present invention;
FIG. 2(b) is a diagram of a clutter power spectrum simulation result obtained by using a conventional SR-STAP method;
FIG. 2(c) is a graph showing the simulation result of clutter power spectrum obtained by the KASR-STAP method;
FIG. 3 is a graph of the signal to noise ratio loss of the output of the method of the present invention, the conventional SR-STAP method and the KASR-STAP method;
FIG. 4 is a comparison graph of the detection capability of the method of the embodiment of the present invention, the conventional SR-STAP method and the KASR-STAP method for a moving object in a sample of distance units to be detected.
Detailed Description
In order to make the technical solutions of the present invention better understood, the present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.
It will be understood that various modifications may be made to the embodiments disclosed herein. Accordingly, the foregoing description should not be construed as limiting, but merely as exemplifications of embodiments. Those skilled in the art will envision other modifications within the scope and spirit of the application.
The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate embodiments of the application and, together with a general description of the application given above and the detailed description of the embodiments given below, serve to explain the principles of the application.
These and other characteristics of the present application will become apparent from the following description of preferred forms of embodiment, given as non-limiting examples, with reference to the attached drawings.
It should also be understood that, although the present application has been described with reference to some specific examples, a person of skill in the art shall certainly be able to achieve many other equivalent forms of application, having the characteristics as set forth in the claims and hence all coming within the field of protection defined thereby.
The above and other aspects, features and advantages of the present application will become more apparent in view of the following detailed description when taken in conjunction with the accompanying drawings.
Specific embodiments of the present application are described hereinafter with reference to the accompanying drawings; however, it is to be understood that the disclosed embodiments are merely examples of the application, which can be embodied in various forms. Well-known and/or repeated functions and constructions are not described in detail to avoid obscuring the application of unnecessary or unnecessary detail. Therefore, specific structural and functional details disclosed herein are not to be interpreted as limiting, but merely as a basis for the claims and as a representative basis for teaching one skilled in the art to variously employ the present application in virtually any appropriately detailed structure.
Fig. 1 is a flowchart of a knowledge-aided SR-STAP method according to an embodiment of the present invention.
As shown in fig. 1, the SR-STAP method based on knowledge assistance provided by the embodiment of the present invention includes:
step S1: and acquiring sample data of the distance unit to be detected through an airborne radar, wherein the airborne radar is in front side view work.
The method comprises the steps that an airborne radar transmits pulse signals to a target to be detected and receives radar echo signals, and the radar echo signals of each distance unit of the airborne radar are recorded as sample data of one distance unit to be detected.
The sample data X of an arbitrary distance unit can be expressed as: clutter + noise, according to the sparse representation model, X can be represented as follows:
X＝Φ×α+n
wherein phi is a space-time steering vector matrix, alpha is a sparse recovery vector, and n is noise. By utilizing the sparsity of the clutter power spectrum on a space-time plane, the clutter can be sparsely decomposed into phi multiplied by alpha.
In this embodiment, the airborne radar is an airborne uniform linear array radar, and the flight speed V of the airborne vehicle and the included angle θ between the array axis and the flight direction of the airborne vehicle are determinedαWhen theta isαWhen the angle is 0 degrees, the side view uniform linear array is formed, and the airborne radar works in the side view. At this time, the spatial cone angle ψ of the clutter and the normalized Doppler frequency f thereofdThere is a coupling relationship between them, and therefore, the clutter should be located on the clutter ridge.
Step S2: and uniformly dispersing the spatial frequency of the clutter into a grid point set.
According to the formula fsUniformly dispersing the normalized spatial frequency into NsSet of grid points { fs,iIn which N iss＝ρsN，ρsAnd dividing coefficients for space frequency grids, wherein N is the number of array elements of the airborne radar.
Step S3: a set of discrete doppler frequency grid points of clutter is acquired.
According to the formula
Obtaining a corresponding discrete Doppler frequency grid point set { fd,iIn which θαIs the angle theta between the array axis and the direction of flight of the aircraftpThe yaw angle at which the beam is transmitted by the radar antenna,
angle of elevation for directing the beam emitted by the radar antenna, beta 4V/lambda frIs a folding factor, V is the aircraft flight speed, lambda is the radar operating wavelength, frIs the pulse repetition frequency. Namely, the Doppler frequency is calculated in a discrete space-time plane by using a knowledge-aided method.
Step S4: and adjusting the Doppler frequency based on a preset condition to obtain a new Doppler frequency group.
Step S4 specifically includes:
(1) the doppler frequency interval is refined based on a first preset condition.
According to the formula
Refining the Doppler frequency interval to ensure the original Doppler frequency interval deltad,iNot greater than spatial frequency interval ΔsWherein, iss＝|fs,i+1-fs,i|，Δd,i＝|fd,i+1-fd,i|,i＝1,...,Ns-1。
Spatial frequency separation ΔsCan be expressed as:
by the above-mentioned spatial frequency spacing ΔsBy the expression (2), the spatial frequency interval DeltasIs constant, therefore, the raw Doppler interval Δ between the i +1 th grid point and the i grid point before the un-refinement is scaledd,iAnd the spatial frequency interval ΔsComparison at Δd,iGreater than ΔsThen, the interval between the ith grid point and the (i + 1) th grid point is divided by an interval deltasAnd (6) uniformly dividing.
(2) Adjacent doppler frequency intervals are combined based on a second preset condition.
According to the formula
Combining adjacent Doppler frequency intervals to ensure the original Doppler frequency interval deltad,iGreater than a ΔsWherein a is an adjusting parameter, and a is more than 0 and less than 1. Due to the original doppler frequency separation Δ obtained from step S3d,iIs not evenly divided, sod,iThe sizes are different, and the adjustment parameter a is valued in the interval of 0 < a < 1, so that deltad,iNot necessarily greater than a Δs. Therefore, in this step, when Δ isd,iLess than a ΔsThen, the ith and (i + 1) th Doppler intervals are adjusted according to the formula
Combining, i.e. may be expressed as combining adjacent doppler frequency intervals.
In this embodiment, by adaptively adjusting the doppler frequency interval of the parameter a and performing theoretical analysis, the optimal value interval of a is as follows:
wherein, Deltad1,mFor the Doppler frequency interval of the echo signal, taking four values of 0.1, 0.2, 0.3 and 0.4 in the optimal value range of a, and obtaining the loss of the output signal-to-noise ratio (SCNR) of the radar signal through experimental verification when the value of a is 0.2Loss) Therefore, in this embodiment, a is 0.2.
(3) And dispersing the Doppler frequency meeting the third preset condition by using space frequency interval.
In this step, when the original Doppler frequency interval is Δd,iSatisfies the condition a Deltas≤Δd,i≤ΔsThe original doppler frequency and separation are unchanged. Using spatial frequency separation ΔsDiscrete coincidence a Δs≤Δd,iThe Doppler frequency is less than or equal to delta. That is, the third preset condition is determined based on the first preset condition and the second preset condition.
(4) Combining the Doppler frequencies adjusted based on the first preset condition, the second preset condition and the third preset condition to obtain a new Doppler frequency group
The Doppler interval is refined and combined and adjusted through the steps (1) and (2) to obtain the condition deltad,i≥ΔsAnd Δd,i≤aΔsAnd (3) the new doppler frequency is compared with the condition a Δ satisfied in step (3)s≤Δd,i≤ΔsCombining the temporally unadjusted Doppler frequencies to obtain a new set of adjusted Doppler frequencies
Step S5: and adjusting the spatial frequency based on the new Doppler frequency group to obtain a new discretization spatial frequency group.
According to the formula
Adjusting the spatial frequency to obtain a new set of discretized spatial frequency
That is, each Doppler frequency is adjusted to obtain a new Doppler frequency in step S4
Then, the space frequency is adjusted by the prior knowledge to obtain each new space frequency to form a new discretization space frequency group
Step S6: according to new Doppler frequency set
And new set of discretized spatial frequencies
And constructing a space-time steering vector dictionary phi by the corresponding space-time steering vectors.
Step S7: and carrying out sparse decomposition on the sample data X by utilizing the space-time guiding vector dictionary phi to obtain a sparse recovery vector alpha.
Step S1 shows that the clutter component can be sparsely decomposed into Φ × α, and after the space-time steering vector dictionary Φ is determined, the sparse vector α can be obtained by matrix inversion.
Step S8: and performing clutter covariance matrix estimation by using the sparse recovery vector.
Specifically, the sample x is replaced with the sparse recovery vector α obtained in step S7lAccording to the formula
Estimating clutter covariance matrix
In some embodiments, after step S8, the method further comprises:
from the estimated clutter covariance matrix
Using formulas
Calculating adaptive filtering weights, wherein StIs a space-time steering vector, mu is a non-zero constant;
and filtering the sample data of the distance unit to be detected according to the calculated adaptive filtering weight.
And performing adaptive filtering on the sample data according to the adaptive filtering weight, so that clutter signals in the sample can be suppressed, and a moving target signal can be detected.
The SR-STAP method based on knowledge assistance provided by the embodiment of the invention combines the knowledge assistance and the SR-STAP method, and performs joint adjustment on the space frequency and the Doppler frequency interval by means of the priori knowledge of the airborne radar array, the space-time coupling of clutter distribution and other priori information, so as to construct an ultra-complete space-time guide vector dictionary, comprehensively correct the dictionary, effectively solve the problem of grid mismatch, improve the sparse recovery precision of clutter and improve the performance of the STAP method. In addition, the invention reduces the number of samples required by the algorithm, and is particularly suitable for moving target detection in a clutter non-uniform environment.
According to the embodiment of the invention, simulation experiments are carried out according to the steps, and experimental verification and analysis are carried out on the method of the embodiment of the invention. Table 1 shows simulation parameters when the airborne radar is in the front-side view operating mode and the simulation experiment is performed according to the above steps.
TABLE 1 simulation parameters
Fig. 2(a) to 2(c) show schematic diagrams of simulation results of clutter power spectrum, wherein fig. 2(a) shows a simulation result diagram of clutter power spectrum calculated by the method of the present invention, fig. 2(b) shows a simulation result diagram of clutter power spectrum obtained by the conventional SR-STAP method, and fig. 2(c) shows a simulation result diagram of clutter power spectrum obtained by the KASR-STAP method. As shown in fig. 2(a) to 2(c), clutter power spectrums estimated by the method of the present invention are all concentrated on clutter ridges, and the occurrence of off-grid effect is suppressed, so that the influence of grid mismatch SR-STAP method is suppressed.
Fig. 3 shows output signal-to-noise ratio loss-contrast diagrams of the method of the present invention, the conventional SR-STAP method and the KASR-STAP method. As shown in fig. 3, the curve of the present invention has small fluctuation, and only the main clutter region generates nulls and has narrow notches, the depth of the curve is deepened by about 4dB compared with the depth of the other two algorithms, and the curve outside the main clutter region is also very smooth, which indicates that the adaptive filtering weight calculated by the method of the present invention is not affected by the off-grid effect, the clutter suppression performance is good, and the detection performance for the moving target is improved.
FIG. 4 shows a comparison of the detection capabilities of the method of the present invention, the conventional SR-STAP method and the KASR-STAP method for a moving object in a sample of range cells to be detected. As shown in fig. 4, the residual clutter power of the method of the present invention is smaller than that of the other two methods. That is to say, the SCNR output by the method is higher than that of the other two methods, the clutter suppression capability is better, and the stronger moving target detection capability is embodied.
Therefore, the knowledge-assisted SR-STAP method can inhibit the influence of grid mismatch on the SR-STAP method, effectively improve the detection capability of the moving target and improve the performance of the SR-STAP method.
An embodiment of the present invention further provides a computer-readable storage medium, on which computer-executable instructions are stored, and when the computer-executable instructions are executed by a processor, the SR-STAP method based on knowledge assistance in the embodiment according to the present invention is implemented.
The processor executing the computer-executable instructions described above may be a processing device, such as a microprocessor, Central Processing Unit (CPU), Graphics Processing Unit (GPU), etc., including one or more general purpose processing devices. More specifically, the processor may be a Complex Instruction Set Computing (CISC) microprocessor, Reduced Instruction Set Computing (RISC) microprocessor, Very Long Instruction Word (VLIW) microprocessor, processor running other instruction sets, or processors running a combination of instruction sets. The processor may also be one or more special-purpose processing devices such as an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA), a Digital Signal Processor (DSP), a system on a chip (SoC), or the like.
The storage medium may be a memory such as Read Only Memory (ROM), Random Access Memory (RAM), phase change random access memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), Electrically Erasable Programmable Read Only Memory (EEPROM), other types of Random Access Memory (RAM), flash disk or other forms of flash memory, cache, registers, static memory, compact disk read only memory (CD-ROM), Digital Versatile Disks (DVD) or other optical storage, magnetic cassettes or other magnetic storage devices, or any other potentially non-transitory medium that may be used to store information or instructions that may be accessed by a computer device, and so forth.
The above embodiments are only exemplary embodiments of the present invention, and are not intended to limit the present invention, and the scope of the present invention is defined by the claims. Various modifications and equivalents may be made by those skilled in the art within the spirit and scope of the present invention, and such modifications and equivalents should also be considered as falling within the scope of the present invention.
Claims (10)
1. A knowledge-aided SR-STAP method is characterized by comprising the following steps:
acquiring sample data of a distance unit to be detected through an airborne radar, wherein the airborne radar is in front side view work;
uniformly dispersing the spatial frequency of the clutter into a grid point set;
acquiring a discrete Doppler frequency grid point set of the clutter;
adjusting the Doppler frequency based on a preset condition to obtain a new Doppler frequency group;
adjusting the spatial frequency based on the new Doppler frequency group to obtain a new discretization spatial frequency group;
constructing a space-time guide vector dictionary according to the new Doppler frequency group and the space-time guide vector corresponding to the new discretization space frequency group;
carrying out sparse decomposition on the sample data by utilizing the space-time guiding vector dictionary to obtain a sparse recovery vector;
and performing clutter covariance matrix estimation by using the sparse recovery vector.
2. The knowledge-aided SR-STAP method according to claim 1, wherein uniformly discretizing the spatial frequency of the clutter into a set of grid points comprises:
according to the formula fsUniformly dispersing the normalized spatial frequency into NsSet of grid points { fs,iIn which N iss＝ρsN，ρsAnd dividing coefficients for space frequency grids, wherein N is the number of array elements of the airborne radar.
3. The knowledge-aided SR-STAP method according to claim 2, wherein obtaining the set of discrete doppler frequency grid points of the clutter comprises:
according to the formula
Obtaining a corresponding discrete Doppler frequency grid point set { fd,iIn which θαIs the angle theta between the array axis and the direction of flight of the aircraftpThe yaw angle at which the beam is transmitted by the radar antenna,
angle of elevation for directing the beam emitted by the radar antenna, beta 4V/lambda frIs a folding factor, V is the aircraft flight speed, lambda is the radar operating wavelength, frIs the pulse repetition frequency.
4. The knowledge-aided SR-STAP method of claim 2, wherein adjusting the doppler frequency based on a preset condition to obtain a new doppler frequency group comprises:
refining the Doppler frequency interval based on a first preset condition;
combining adjacent Doppler frequency intervals based on a second preset condition;
dispersing Doppler frequency meeting a third preset condition by using space frequency interval;
and combining the Doppler frequencies adjusted based on the first preset condition, the second preset condition and the third preset condition to obtain a new Doppler frequency group.
5. The knowledge-aided SR-STAP method according to claim 4, wherein said refining Doppler frequency interval based on a first preset condition comprises:
according to the formula
Refining the Doppler frequency interval to make the original Doppler frequency interval deltad,iNot greater than spatial frequency interval ΔsWherein, iss＝|fs,i+1-fs|，Δd,i＝|fd,i+1-fd,i|,i＝1,...,Ns-1。
6. The knowledge-aided SR-STAP method according to claim 5, wherein the merging of adjacent doppler frequency intervals based on a second preset condition comprises:
according to the formula
Combining adjacent Doppler frequency intervals to ensure that the original Doppler frequency interval is greater than a deltasWherein a is an adjusting parameter, and a is more than 0 and less than 1.
7. The knowledge-aided SR-STAP method according to claim 6, wherein said discretizing Doppler frequencies meeting a third predetermined condition by spatial frequency separation comprises:
using spatial frequency separation ΔsDiscrete coincidence a Δs≤Δd,i≤ΔsThe doppler frequency of the condition.
8. The knowledge-aided SR-STAP method of claim 1, wherein adjusting spatial frequencies based on the new doppler frequency set to obtain a new discretized spatial frequency set comprises:
according to the formula
Adjusting the spatial frequency to obtain a new set of discretized spatial frequency
9. The knowledge-aided SR-STAP method according to claim 1, wherein the method further comprises:
from the estimated clutter covariance matrix
Using formulas
Calculating adaptive filtering weights, wherein StIs a space-time steering vector, mu is a non-zero constant;
and filtering the sample data of the distance unit to be detected according to the calculated adaptive filtering weight.
10. A computer-readable storage medium having stored thereon computer-executable instructions that, when executed by a processor, implement the knowledge-based assisted SR-STAP method according to any of claims 1-9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110577866.9A CN113219433B (en) | 2021-05-26 | 2021-05-26 | Knowledge-aided SR-STAP method and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110577866.9A CN113219433B (en) | 2021-05-26 | 2021-05-26 | Knowledge-aided SR-STAP method and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113219433A true CN113219433A (en) | 2021-08-06 |
| CN113219433B CN113219433B (en) | 2023-03-31 |
Family
ID=77098647
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110577866.9A Active CN113219433B (en) | 2021-05-26 | 2021-05-26 | Knowledge-aided SR-STAP method and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113219433B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113655458A (en) * | 2021-09-02 | 2021-11-16 | 内蒙古工业大学 | Dictionary correction-based space-time adaptive processing method and device and storage medium |
| CN113820679A (en) * | 2021-08-10 | 2021-12-21 | 西安电子科技大学 | Filtering method, device, equipment and storage medium of radar echo data |
Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7193558B1 (en) * | 2003-09-03 | 2007-03-20 | The United States Of America As Represented By The Secretary Of The Navy | Radar processor system and method |
| CN102313885A (en) * | 2010-06-30 | 2012-01-11 | 中国科学院电子学研究所 | Multi-dimensional microwave imaging system and method |
| US20140035776A1 (en) * | 2011-04-14 | 2014-02-06 | Nir Asher Sochen | Method and system for target detection |
| CN104105449A (en) * | 2011-12-01 | 2014-10-15 | 毛伊图像公司 | Motion detection using ping-based and multiple aperture Doppler ultrasound |
| CN104345300A (en) * | 2014-10-30 | 2015-02-11 | 河海大学 | Onboard non-positive side view array radar STAP (Space Time Adaptive Processing) method for clutter space-time spectrum linear compensation |
| CN105223560A (en) * | 2015-10-13 | 2016-01-06 | 中国人民解放军空军工程大学 | Based on the airborne radar object detection method of the sparse recovery of clutter pitching azimuth spectrum |
| CN106324596A (en) * | 2016-09-09 | 2017-01-11 | 深圳大学 | Random pulse Doppler radar angle-Doppler imaging method based on compressed sensing |
| CN107544061A (en) * | 2017-09-11 | 2018-01-05 | 深圳大学 | Knowledge based engineering is sparse to recover space-time adaptive processing method and system |
| CN108931766A (en) * | 2018-04-28 | 2018-12-04 | 河海大学 | A kind of non-homogeneous STAP jamming target filtering method based on sparse reconstruct |
| CN109116311A (en) * | 2018-09-19 | 2019-01-01 | 西安电子科技大学 | Knowledge based assists the clutter suppression method of sparse iteration covariance estimation |
| CN109375179A (en) * | 2018-10-29 | 2019-02-22 | 中国电子科技集团公司第十四研究所 | A kind of strong clutter suppression method of airborne radar short range based on rarefaction representation |
| CN109557539A (en) * | 2018-12-27 | 2019-04-02 | 中国人民解放军空军工程大学 | Airborne Passive radar Ground moving targets detection method, apparatus and electronic equipment |
| CN110703209A (en) * | 2019-09-30 | 2020-01-17 | 西安电子科技大学 | Method for suppressing ground distance fuzzy clutter of high repetition frequency airborne forward looking array radar |
| CN110764069A (en) * | 2019-11-14 | 2020-02-07 | 内蒙古工业大学 | Sparse recovery STAP color loading method based on knowledge assistance |
| CN110907910A (en) * | 2019-11-27 | 2020-03-24 | 中国船舶重工集团公司第七二四研究所 | Distributed coherent radar moving target echo coherent synthesis method |
| CN111308431A (en) * | 2020-02-27 | 2020-06-19 | 西北大学 | Two-dimensional two-pulse cancellation method based on estimation error |
| CN111580065A (en) * | 2020-07-06 | 2020-08-25 | 内蒙古工业大学 | Sparse STAP recovery method based on knowledge assistance |
| CN112415476A (en) * | 2020-11-13 | 2021-02-26 | 中国民航大学 | Dictionary mismatch clutter space-time spectrum estimation method based on sparse Bayesian learning |
| CN112415475A (en) * | 2020-11-13 | 2021-02-26 | 中国民航大学 | Non-grid sparse recovery non-front side array STAP method based on atomic norm |
-
2021
- 2021-05-26 CN CN202110577866.9A patent/CN113219433B/en active Active
Patent Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7193558B1 (en) * | 2003-09-03 | 2007-03-20 | The United States Of America As Represented By The Secretary Of The Navy | Radar processor system and method |
| CN102313885A (en) * | 2010-06-30 | 2012-01-11 | 中国科学院电子学研究所 | Multi-dimensional microwave imaging system and method |
| US20140035776A1 (en) * | 2011-04-14 | 2014-02-06 | Nir Asher Sochen | Method and system for target detection |
| CN104105449A (en) * | 2011-12-01 | 2014-10-15 | 毛伊图像公司 | Motion detection using ping-based and multiple aperture Doppler ultrasound |
| CN104345300A (en) * | 2014-10-30 | 2015-02-11 | 河海大学 | Onboard non-positive side view array radar STAP (Space Time Adaptive Processing) method for clutter space-time spectrum linear compensation |
| CN105223560A (en) * | 2015-10-13 | 2016-01-06 | 中国人民解放军空军工程大学 | Based on the airborne radar object detection method of the sparse recovery of clutter pitching azimuth spectrum |
| CN106324596A (en) * | 2016-09-09 | 2017-01-11 | 深圳大学 | Random pulse Doppler radar angle-Doppler imaging method based on compressed sensing |
| CN107544061A (en) * | 2017-09-11 | 2018-01-05 | 深圳大学 | Knowledge based engineering is sparse to recover space-time adaptive processing method and system |
| CN108931766A (en) * | 2018-04-28 | 2018-12-04 | 河海大学 | A kind of non-homogeneous STAP jamming target filtering method based on sparse reconstruct |
| CN109116311A (en) * | 2018-09-19 | 2019-01-01 | 西安电子科技大学 | Knowledge based assists the clutter suppression method of sparse iteration covariance estimation |
| CN109375179A (en) * | 2018-10-29 | 2019-02-22 | 中国电子科技集团公司第十四研究所 | A kind of strong clutter suppression method of airborne radar short range based on rarefaction representation |
| CN109557539A (en) * | 2018-12-27 | 2019-04-02 | 中国人民解放军空军工程大学 | Airborne Passive radar Ground moving targets detection method, apparatus and electronic equipment |
| CN110703209A (en) * | 2019-09-30 | 2020-01-17 | 西安电子科技大学 | Method for suppressing ground distance fuzzy clutter of high repetition frequency airborne forward looking array radar |
| CN110764069A (en) * | 2019-11-14 | 2020-02-07 | 内蒙古工业大学 | Sparse recovery STAP color loading method based on knowledge assistance |
| CN110907910A (en) * | 2019-11-27 | 2020-03-24 | 中国船舶重工集团公司第七二四研究所 | Distributed coherent radar moving target echo coherent synthesis method |
| CN111308431A (en) * | 2020-02-27 | 2020-06-19 | 西北大学 | Two-dimensional two-pulse cancellation method based on estimation error |
| CN111580065A (en) * | 2020-07-06 | 2020-08-25 | 内蒙古工业大学 | Sparse STAP recovery method based on knowledge assistance |
| CN112415476A (en) * | 2020-11-13 | 2021-02-26 | 中国民航大学 | Dictionary mismatch clutter space-time spectrum estimation method based on sparse Bayesian learning |
| CN112415475A (en) * | 2020-11-13 | 2021-02-26 | 中国民航大学 | Non-grid sparse recovery non-front side array STAP method based on atomic norm |
Non-Patent Citations (7)
| Title |
|---|
| GAO, ZQ等: "Robust STAP algorithm based on knowledge-aided SR for airborne radar", 《IET RADAR SONAR AND NAVIGATION 》 * |
| LI ZHIHUI: "Off-grid STAP algorithm based on local search orthogonal matching pursuit", 《SYSTEMS ENGINEERING AND ELECTRONICS 》 * |
| PENGYUAN HE,SHUN HE,ZHIWEI YANG,PENGHUI HUANG: "An Off-Grid STAP Algorithm Based On Local Mesh Splitting With Bistatic Radar System", 《IEEE》 * |
| QIANG WANG等: "Clutter Reduced-Dimension Sparse Recovery Method on Knowledge-Aided for Airborne Phased Array Radar", 《IEEE》 * |
| 张欢欢,高志奇,黄平平,徐伟: "基于知识辅助的网格失配下SR-STAP字典校正方法", 《信号处理》 * |
| 林雪芳: "知识辅助空时自适应处理技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 高志奇: "机载雷达的稳健稀疏空时自适应处理方法研究", 《中国博士学位论文全文数据库信息科技辑》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113820679A (en) * | 2021-08-10 | 2021-12-21 | 西安电子科技大学 | Filtering method, device, equipment and storage medium of radar echo data |
| CN113820679B (en) * | 2021-08-10 | 2023-12-26 | 西安电子科技大学 | Filtering method, device, equipment and storage medium for radar echo data |
| CN113655458A (en) * | 2021-09-02 | 2021-11-16 | 内蒙古工业大学 | Dictionary correction-based space-time adaptive processing method and device and storage medium |
| CN113655458B (en) * | 2021-09-02 | 2023-06-02 | 内蒙古工业大学 | Space-time adaptive processing method and device based on dictionary correction and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113219433B (en) | 2023-03-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113219433B (en) | Knowledge-aided SR-STAP method and storage medium | |
| CN110045323B (en) | Matrix filling-based co-prime matrix robust adaptive beamforming algorithm | |
| Degurse et al. | Reduced-rank STAP for target detection in heterogeneous environments | |
| CN111948619B (en) | Passive cooperative detection method and system for targets under irradiation of multiple types of external radiation sources | |
| CN105302936A (en) | Self-adaptive beam-forming method based on related calculation and clutter covariance matrix reconstruction | |
| CN108919220B (en) | Missile-borne SAR front-side-view imaging method based on embedded GPU | |
| CN103984676A (en) | Rectangular projection adaptive beamforming method based on covariance matrix reconstruction | |
| CN113219432B (en) | Moving object detection method based on knowledge assistance and sparse Bayesian learning | |
| CN105335336A (en) | Sensor array steady adaptive beamforming method | |
| CN113030843B (en) | Low-altitude target direction finding method based on improved IAA beam forming | |
| CN105044705A (en) | Multibeam depth measuring sonar multi-subarray wave beam sharpening method | |
| CN108931766A (en) | A kind of non-homogeneous STAP jamming target filtering method based on sparse reconstruct | |
| CN113422630A (en) | Self-adaptive focusing broadband beam forming method and system | |
| CN109212466B (en) | Quantum dragonfly evolution mechanism-based broadband direction finding method | |
| CN112763972B (en) | Sparse representation-based double parallel line array two-dimensional DOA estimation method and computing equipment | |
| Li et al. | Bayesian linear regression with cauchy prior and its application in sparse mimo radar | |
| Aboutanios et al. | Fast iterative interpolated beamforming for high fidelity single snapshot DOA estimation | |
| CN110850421A (en) | Underwater target detection method based on space-time adaptive processing of reverberation symmetric spectrum | |
| CN109490859B (en) | Robust detector for orientation and Doppler disturbances in a partially homogeneous environment | |
| Guan et al. | A CFAR detector for MIMO array radar based on adaptive pulse compression-Capon filter | |
| CN114563760B (en) | Second-order super-beam forming method, equipment and medium based on SCA array | |
| CN113702899B (en) | Propagation algorithm of covariance difference based on phase fraction low-order moment | |
| Wang et al. | Fast and adaptive method for SAR superresolution imaging based on point scattering model and optimal basis selection | |
| CN115575906A (en) | Fusion detection method for intelligent interference suppression under non-uniform background | |
| CN113655458B (en) | Space-time adaptive processing method and device based on dictionary correction and storage medium |
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 |