Missile-borne phased array multichannel radar clutter data refined simulation method

Abstract

The invention particularly relates to a missile-borne phased array multichannel radar clutter data refinement simulation method, which comprises the following steps of: obtaining coordinates of a radar transmitter and receiving channels at each pulse moment according to a trajectory of a missile body platform within coherent accumulation time; dividing a detection scene into super-resolution ground points, and combining the super-resolution ground points according to the slant-range resolution and the Doppler resolution in sequence to obtain each distance ring and a corresponding scattering unit; then, calculating a directional diagram gain and clutter energy of each scattering unit; and finally, performing four-stage element-by-element simulation by adopting a channel-pulse-distance ring-scattering unit mode based on a strict slant range process, and performing echo superposition on the scattering units of the same distance ring to obtain clutter data of each pulse and each distance ring under each receiving channel. Through the method, high-fidelity missile-borne phased array multichannel radar clutter echo data under a given detection scene and a radar working condition can be obtained, and simulation data support is provided for further research on a missile-borne STAP technology.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of signal processing, in particular to a missile-borne phased array multi-channel radar clutter data refinement simulation method, which can be used for missile-borne phased array multi-channel radar. BACKGROUND

[0002] Space-time adaptive processing (STAP) is an effective technology widely used in airborne / spaceborne early warning radar moving target detection. Based on array multi-channel theory, the technology realizes adaptive suppression of clutter by combining two-dimensional data in space and time. In missile-borne radar, due to the influence of objective factors such as platform size limitation and strong motion maneuverability, STAP technology has not been widely used. Traditional missile-borne radar usually adopts non-adaptive Doppler velocity clutter suppression processing methods such as moving target detection (MTD) and pulse Doppler (PD), which are difficult to meet the increasingly urgent demand for weak slow moving target detection in strong clutter background. Therefore, applying multi-channel theory and STAP technology to missile-borne platforms is an important direction for future development.

[0003] The performance of STAP technology is closely related to the complexity of clutter multi-channel characteristics. At present, domestic and foreign researches have been carried out on the clutter characteristics of missile-borne phased array multi-channel radar in non-normal side array and different motion states, and a series of improved STAP methods suitable for missile-borne platforms have been proposed, such as registration compensation method based on radar parameters, robust difference STAP method, etc.

[0004] However, the clutter signal model used in the current research is a space-time steering vector model, which does not consider the influence of ground scattering element resolution and high-speed motion of the missile body, and the obtained clutter data characteristics may be different from the real echo. Obtaining missile-borne multi-channel radar clutter measured data requires a lot of manpower and material resources. In order to quickly evaluate the performance of STAP technology under different working conditions in the early stage of engineering demonstration, it is an efficient means to obtain high-fidelity clutter data through modeling simulation, so it is urgent to establish a complete and accurate missile-borne phased array multi-channel radar clutter data simulation method.

[0005] It should be noted that the information disclosed in the above background section is only used to strengthen the understanding of the background of the present application, and therefore can include information that does not constitute prior art known to those of ordinary skill in the art. SUMMARY

[0006] In order to solve the above problems existing in the prior art, the application provides a fine simulation method for clutter data of a missile-borne phased array multi-channel radar, which combines the characteristics of the missile-borne radar, is based on a strict slant range history, simultaneously considers the influence of an actual resolution, and strictly simulates the radar signal transmission-reflection-reception process through a way of sequentially simulating echo data of each channel, each pulse, each distance ring, and each scattering unit, so that high-fidelity clutter echo data under a given detection scene and radar working condition can be obtained.

[0007] Other characteristics and advantages of the application will become apparent from the following detailed description, or will be learned by practice of the application.

[0008] According to a first aspect of the application, a fine simulation method for clutter data of a missile-borne phased array multi-channel radar is provided, and the method comprises: obtaining coordinates of a radar transmitter and each receiving channel at each pulse moment according to a trajectory of a missile platform within a coherent accumulation time; dividing the detection scene into super-resolution ground points, sequentially merging according to a slant range resolution and a Doppler resolution, obtaining each distance ring and a corresponding scattering unit, and calculating a center coordinate and a scattering cross-sectional area of the scattering unit; calculating a transmission pattern gain, a receiving pattern gain, and clutter energy of each scattering unit; generating a transmission signal in a frequency domain, calculating a transmission-reception two-way slant range of the scattering unit based on the coordinates of the radar transmitter and each receiving channel at each pulse moment, simulating frequency domain clutter echo data of each receiving channel, each pulse, each distance ring, and each scattering unit based on the transmission-reception two-way slant range, superimposing echo data of all scattering units on the same distance ring to obtain clutter data of each receiving channel, each pulse, and each distance ring, and performing inverse Fourier transform to obtain time domain clutter data.

[0009] In some exemplary embodiments, the method for obtaining the coordinates of the radar transmitter and each receiving channel at each pulse moment according to the trajectory of the missile platform within the coherent accumulation time specifically comprises: the trajectory of the missile platform comprises a motion velocity and an acceleration , the coordinates of the missile-borne phased array multi-channel radar transmitter at each pulse moment and the coordinates of each receiving channel are calculated based on the motion velocity and the acceleration , and the method specifically comprises: calculating the coordinates of the radar transmitter at each pulse moment and the coordinates of each receiving channel .

[0010]

[0011] wherein,​ for The transmitter coordinates at the nth pulse moment, the nth pulse moment For the first The pulse moment The coordinates of each receiving channel.

[0012] In some exemplary embodiments, the step of dividing the detection scene into super-resolution ground points and merging them sequentially according to slant range resolution and Doppler resolution to obtain each range ring and its corresponding scattering unit specifically involves: For detection scenarios in both azimuth and elevation directions, super-resolution ground points are defined according to the highest beam resolution of the receiving antenna; Super-resolution ground points belonging to the same slant range resolution cell are merged, resulting in a total of [number missing]. The first distance ring is used to obtain the super-resolution ground points contained in each distance ring. The set of super-resolution ground points contained in a distance ring is denoted as ; For each distance ring Within a single Doppler resolution unit, super-resolution ground points are further merged, resulting in a total of [number missing]. The scattering unit is used to obtain the super-resolution ground points contained in each scattering unit within each range ring. On the distance ring, the first The set of super-resolution ground points contained in a scattering unit is denoted as ; Calculate the center coordinates and scattering cross-section of each scattering unit on each range ring.

[0013] In some exemplary embodiments, the calculation of the emission pattern gain of each scattering unit is performed using the following formula:

[0014] in, Indicates the first On the distance ring, the first The emission pattern of each scattering unit and These represent the number of array elements in the azimuth and elevation directions of the antenna, respectively. and These represent the weighting coefficients of the row submatrix and the column submatrix, respectively. Indicates the spacing between array elements. Indicates the signal wavelength. and They represent the first On the distance ring, the first The azimuth and elevation angles of each scattering unit. and These represent the azimuth and elevation angles of the beam center, respectively.

[0015] In some exemplary embodiments, the calculation of the receive pattern gain for each scattering element is performed using the following formula:

[0016] in, Indicates the first On the distance ring, the first The receiving pattern of each scattering unit.

[0017] In some exemplary embodiments, the calculation of the clutter energy of each scattering unit is performed using the following formula:

[0018] in, Indicates the first On the distance ring, the first clutter energy of each scattering unit Indicates average transmit power. and These represent the transmit and receive antenna gains, respectively. Indicates the pulse repetition period. Indicates system loss. Indicates the signal pulse width.

[0019] In some exemplary embodiments, the simulated frequency domain clutter echo data for each receiving channel, each pulse, each range loop, and each scattering unit is derived using the following formula:

[0020] in, and They represent the first Inter-channel amplitude and phase errors of each channel Indicates signal frequency. This represents the transmit / receive slant range of the scattering unit. Indicates the first The receiving channel, the first The pulse, the first The distance ring, the first Frequency domain echo signal of each scattering unit.

[0021] According to a second aspect of the present invention, a storage medium is provided having a computer program stored thereon, which, when executed by a processor, implements the refined simulation method for missile-borne phased array multi-channel radar clutter data described in the first aspect.

[0022] According to a third aspect of the present application, a computer program product is provided, and the computer program product has a computer program stored thereon, and the computer program is executed by a processor to implement the missile-borne phased array multi-channel radar clutter data refinement simulation method of the first aspect.

[0023] According to a fourth aspect of the present application, an electronic device is provided, and the electronic device comprises: a processor; and a memory configured to store executable instructions of the processor; wherein the processor is configured to implement the missile-borne phased array multi-channel radar clutter data refinement simulation method of the first aspect by executing the executable instructions.

[0024] The missile-borne phased array multi-channel radar clutter data refinement simulation method provided by the embodiments of the present application has the following advantages compared with the prior art: 1. The missile-borne phased array multi-channel radar clutter data refinement simulation method provided by the embodiments of the present application is combined with the characteristics of a missile-borne radar, is based on a strict slant range history, considers the influence of an actual resolution, and simulates echo data in a way of channel by channel, pulse by pulse, distance ring by distance ring, and scattering unit by scattering unit, so that the missile-borne phased array multi-channel radar clutter data refinement simulation method can strictly simulate a radar signal transmission-reflection-reception process and make up for the limitations of the prior art that is based on a space-time steering vector and does not consider two-dimensional resolution.

[0025] 2. The echo data clutter energy, clutter Doppler spectrum width, clutter multi-channel synthesis effect, and clutter space-time characteristics obtained by simulation are consistent with theoretical results, so that the missile-borne phased array multi-channel radar clutter data refinement simulation method can accurately and effectively obtain high-fidelity clutter echo data of a given scene and provide simulation data support for further development of missile-borne phased array multi-channel radar detection capability evaluation and STAP technology research.

[0026] It should be understood that the foregoing general description and the following detailed description are only exemplary and explanatory, and cannot limit the present application. BRIEF DESCRIPTION OF DRAWINGS

[0027] The accompanying drawings, which are incorporated into and form part of the specification, illustrate embodiments consistent with the present application and, together with the specification, serve to explain the principles of the present application. Obviously, the drawings in the following description are only some embodiments of the present application, and other drawings can be obtained by those skilled in the art without creative labor on the basis of these drawings.

[0028] Figure 1 is a flowchart of a missile-borne phased array multi-channel radar clutter data refinement simulation method of the present application; Figure 2 is a schematic diagram of the relationship between a missile platform, an antenna array surface and a detection scene in the embodiments of the present application; Figure 3This is the clutter distance-Doppler plot of the first receiving channel in this embodiment of the invention; Figure 4 This is a range-Doppler distribution diagram of each scattering unit in an embodiment of the present invention; Figure 5 This is the azimuth-elevation two-dimensional normalized transmit / receive direction map of the detection scene in this embodiment of the invention; Figure 6 This is a clutter normalized energy distribution diagram of each scattering unit in the embodiment of the present invention; Figure 7 This is the spatial characteristic spectrum of clutter after channel coherent synthesis processing in this embodiment of the invention; Figure 8 This is a comparison diagram of the target energy before and after channel coherent synthesis processing in an embodiment of the present invention; Figure 9 This is a spatiotemporal spectrum image of clutter in an embodiment of the present invention; Figure 10 This is a schematic diagram of the theoretical spatial frequency-Doppler curve of clutter in an embodiment of the present invention. Detailed Implementation

[0029] Exemplary embodiments will now be described more fully with reference to the accompanying drawings. However, these exemplary embodiments can be implemented in many forms and should not be construed as limited to the examples set forth herein; rather, they are provided so that the invention will be more comprehensive and complete, and will fully convey the concept of the exemplary embodiments to those skilled in the art. The described features, structures, or characteristics may be combined in any suitable manner in one or more embodiments.

[0030] Furthermore, the accompanying drawings are merely illustrative of the invention and are not necessarily drawn to scale. The same reference numerals in the drawings denote the same or similar parts, and therefore repeated descriptions of them will be omitted. Some block diagrams shown in the drawings are functional entities and do not necessarily correspond to physically or logically independent entities. These functional entities can be implemented in software, in one or more hardware modules or integrated circuits, or in different network and / or processor devices and / or microcontroller devices.

[0031] To address the shortcomings and deficiencies of existing technologies, this example embodiment provides a refined simulation method for clutter data from a missile-borne phased array multi-channel radar. This simulation method includes: obtaining the coordinates of the radar transmitter and each receiving channel at each pulse moment based on the trajectory of the missile platform during the coherent accumulation time; dividing the detection scene into super-resolution ground points, and merging them sequentially according to slant range resolution and Doppler resolution to obtain each range loop and its corresponding scattering unit; then calculating the pattern gain and clutter energy of each scattering unit; finally, based on a rigorous slant range history, performing a four-level element-wise simulation using a "channel-pulse-range loop-scattering unit" approach, and superimposing the echoes of scattering units within the same range loop to obtain clutter data for each pulse and each range loop under each receiving channel. This invention can obtain high-fidelity clutter echo data from a missile-borne phased array multi-channel radar under a given detection scene and radar operating conditions, providing simulation data support for further research on missile-borne STAP technology.

[0032] Please see Figure 1 , Figure 1 This is a flowchart of a method for refined simulation of clutter data from a missile-borne phased array multi-channel radar, provided by an embodiment of the present invention. The method includes the following steps: Step 1: Based on the actual trajectory of the missile platform within one coherent accumulation time, calculate the coordinates of the onboard phased array multi-channel radar transmitter and each receiving channel at each pulse moment.

[0033] Furthermore, based on the velocity of the projectile platform during the coherent accumulation time... and acceleration ,calculate The coordinates of the missile-borne phased array multi-channel radar transmitter at each pulse moment and Coordinates of each receiving channel ; The specific explanation is as follows: Reference Figure 2 This is a schematic diagram illustrating the relationship between the projectile platform, antenna array, and detection scenario in an embodiment of the present invention. With the center of mass of the projectile platform as the origin, the axial direction of the projectile is... The axis, perpendicular to the projectile's axis and within the projectile's plane, is... axis, The axes are determined by the right-hand screw rule to establish a spatial rectangular coordinate system. Therefore, the initial coordinates of the radar transmitter can be expressed as:

[0034] Assuming the antenna array of the missile-borne phased array multi-channel radar is composed of It consists of several array elements, with the spacing between each element being [missing information]. ,in Indicates the signal wavelength. and These represent the number of array elements in the azimuth and elevation directions of the antenna, respectively. In signal processing, the antenna arrays are typically first synthesized to form a matrix... A linear array consisting of n equivalent channels, at the initial time n The coordinates of each receiving channel can be represented as:

[0035] in, Indicates the first The distance between each receiving channel and the center of the antenna array.

[0036] The velocity and acceleration vectors of the projectile platform can be represented in this coordinate system as:

[0037]

[0038] in, , and Indicates the projectile in the coordinate system axis, shaft and Magnitude of velocity in the axial direction; , and Indicates the projectile in the coordinate system axis, shaft and The magnitude of acceleration in the axial direction.

[0039] Assuming coherent accumulation time Inside, the missile-borne phased array radar launched a total of A pulse, then in The coordinates of the radar receiving channel at each moment within a given time period can be represented as follows:

[0040]

[0041] The coordinates of the radar transmitter at each moment can be represented as:

[0042]

[0043] Step 2: Divide the detection scene into super-resolution ground points, and merge them according to slant range resolution and Doppler resolution to obtain each range ring and the scattering units it contains, and calculate the center coordinates and scattering cross-section of the scattering units.

[0044] Furthermore, the ranges for the azimuth and pitch directions are respectively and The detection scenario is based on the highest beam resolution of the receiving antenna. Super-resolution ground points are defined; then, points belonging to a single slant-range resolution cell are identified. The super-resolution ground points were merged, resulting in a total of [number] points. The first distance ring is used to obtain the super-resolution ground points contained in each distance ring. The set of super-resolution ground points contained in a distance ring is denoted as For each distance ring Within, it belongs to a Doppler resolution unit. The super-resolution ground points are then merged to form a total of [number]. The scattering unit is used to obtain the super-resolution ground points contained in each scattering unit within each range ring. On the distance ring, the first The set of super-resolution ground points contained in a scattering unit is denoted as Finally, the center coordinates of each scattering unit on each range ring are calculated. and scattering cross-section .

[0045] The specific explanation is as follows: Step 2.1: Within the detection scene range, use the highest beam resolution of the receiving antenna. The initial detection scene is split to obtain super-resolution ground points.

[0046] The angular ranges corresponding to the azimuth and range swaths of the detection scene relative to the antenna array are respectively and The azimuth and elevation angles corresponding to each super-resolution ground point are respectively

[0047]

[0048] Then the first The coordinates of a super-resolution ground point can be represented as:

[0049] in, Indicates the first The slant range of each scattering unit, and They represent the first The elevation and azimuth angles of each scattering unit.

[0050] Step 2.2: For cells belonging to a single slant range resolution unit The super-resolution ground points are merged into one. Each distance loop contains super-resolution ground points.

[0051] No. The set of super-resolution ground points contained in the distance ring It can be represented as:

[0052] in, Indicates the minimum slant range of the detection scene; Indicates slant range resolution; Indicates the first The slant distance of each super-resolution ground point.

[0053] Step 2.3: Within each range ring, super-resolution ground points belonging to a single Doppler resolution cell are merged into a single super-resolution ground point. ( ( ) scattering units. Obtain the super-resolution ground points contained in each scattering unit within each range ring.

[0054] No. On the distance ring, the first The subscripts of the super-resolution ground points contained in a scattering unit can be represented as:

[0055] in, Indicates Doppler resolution. Indicates the first The minimum Doppler frequency of the distance loop, Indicates the first On the distance ring, the first Doppler frequencies at super-resolution ground points.

[0056] Step 2.4: Calculate the center coordinates and scattering cross-section of each scattering unit on each range ring.

[0057] No. On the distance ring, the first The center coordinates of a scattering unit can be represented as:

[0058] in, Indicates the first On the distance ring, the first The number of super-resolution ground points contained in a scattering unit Indicates the first On the distance ring, the first The scattering unit contains the first The coordinates of a super-resolution ground point.

[0059] No. On the distance ring, the first The target scattering cross-section of a scattering unit can be expressed as:

[0060] in, Indicates the first On the distance ring, the first The scattering coefficient of each scattering unit Indicates the first On the distance ring, the first The area of ​​each scattering unit.

[0061] Step 3: Calculate the pattern gain and clutter energy of each scattering unit on each range ring.

[0062] Furthermore, the emission pattern gain of each scattering unit is calculated. Receiver pattern gain and clutter energy .

[0063] The specific explanation is as follows: No. On the distance ring, the first The emission pattern of a scattering unit can be represented as:

[0064] in, and These represent the number of array elements in the azimuth and elevation directions of the antenna, respectively. and These represent the weighting coefficients of the row submatrix and the column submatrix, respectively. Indicates the spacing between array elements. Indicates the signal wavelength. and They represent the first On the distance ring, the first The azimuth and elevation angles of each scattering unit. and These represent the azimuth and elevation angles of the beam center, respectively.

[0065] No. On the distance ring, the first The receiving pattern of a scattering unit can be represented as:

[0066] Therefore, the first On the distance ring, the first The clutter energy of a scattering unit can be expressed as:

[0067] in, Indicates average transmit power. and These represent the transmit and receive antenna gains, respectively. Indicates the pulse repetition period. Indicates system loss. Indicates the signal pulse width.

[0068] Step 4: Generate frequency domain transmission signals and simulate frequency domain clutter echo signals channel by channel, pulse by pulse, range loop by range loop, and scattering unit by scattering unit. Superimpose the frequency domain clutter echo signals of all scattering units on the same range loop and perform inverse Fourier transform to obtain the time domain clutter echo signals of each channel, pulse, and range loop.

[0069] Furthermore, a transmitted signal is generated in the frequency domain. Simulated frequency domain clutter echo data for each receiving channel, each pulse, each range loop, and each scattering unit. By superimposing the echo data from all scattering units on the same range loop, clutter data for each pulse and each range loop under each receiving channel can be obtained. Performing an inverse Fourier transform yields the time-domain clutter data. .

[0070] The specific explanation is as follows: Step 4.1: Generate a spurious frequency domain transmission signal .

[0071] Step 4.2: Initialize the simulation channels .

[0072] Step 4.3: Simulation to obtain the first Time-domain clutter echo data for all pulses and all distance loops in each channel.

[0073] Furthermore, step 4.3 specifically includes: Step 4.3.1: Initialize the simulation pulse .

[0074] Step 4.3.2: Calculate the first... At the pulse moment, the first The first channel to the first On the distance ring, the first The transmit / receive slant range of each scattering unit:

[0075] Step 4.3.3: Simulation of the first step The receiving channel, the first The pulse, the first The distance ring, the first The frequency domain echo signal of each scattering unit is:

[0076] in, and They represent the first Inter-channel amplitude and phase errors of each channel Indicates signal frequency. This represents the slant distance for both transmission and reception of the scattering unit.

[0077] Step 4.3.4: [The text appears to be incomplete and contains several grammatical errors. A more accurate translation The echo signals from all scattering units on the range ring are superimposed to obtain the first... At the pulse moment, the first The first channel One distance-loop frequency domain echo signal:

[0078] Step 4.3.5: [The text appears to be incomplete and contains several grammatical errors. A more accurate translation would require the full context.] At the pulse moment, the first The first channel Transformation of distance loop echo data to the time domain:

[0079] Among them, symbols This indicates the inverse Fourier transform operation.

[0080] Step 4.3.6: Update the simulated pulse sequence number Repeat steps 4.3.3 to 4.3.5 until... . (To obtain the first) Time-domain clutter echo data for all pulses and all distance loops in each channel.

[0081] Step 4.3.7: Additive noise can be added to the clutter echo signal. ,in Represents the Boltzmann constant. Indicates temperature. Represents the noise figure. Indicates by A random complex number, corresponding to pulse and There are three distance rings, both real and imaginary, with a mean of 0 and a variance of 1.

[0082] Step 4.4: Update the channel numbers in the simulation. Repeat step 4.3 until... .

[0083] The effectiveness of the present invention's embodiment based on a refined simulation method for clutter data from a missile-borne phased array multi-channel radar can be illustrated by the following simulation experiments: (1) Experimental conditions: The equivalent receiving channels of the missile-borne phased array multi-channel radar are 4, and they are evenly arranged; the azimuth and elevation dimensions of the antenna are both 0.25m; the azimuth and elevation angles of the beam center are 6° and 80°, respectively; the yaw, elevation and roll angles of the antenna installation are 90°, 0° and 0°, respectively; the height of the missile platform is 1km, moving horizontally at a speed of Mach 2 and an acceleration of 0; the detection scenario is the area contained by the main lobe of the antenna in the elevation direction and the first side lobe in the azimuth direction, with the elevation and azimuth angles of the beam center being 6° and 10°, respectively, and the terrain being sea state 5; the signal form is LFM signal, the signal pulse width is 30μs, the carrier frequency is 16GHz, the average transmit power is 4.2KW, the system loss is 3dB, the pulse repetition frequency is 10000Hz, the coherent accumulation time is 30ms, and there are no amplitude errors or phase errors between channels.

[0084] (2) Experimental content and results analysis: Experiment 1: In this experiment, the clutter echo data was generated using a refined simulation method for multi-channel phased array radar clutter data proposed in this invention. After pulse compression, the range-Doppler plot of the first receiving channel was plotted. Please refer to [link to relevant documentation]. Figure 3 This is the clutter range-Doppler plot for the first receiving channel. It can be seen that... Figure 3 Due to pattern modulation, the clutter distance-Doppler graph shows that the clutter energy in the near-range and far-range sidelobe regions is reduced to noise levels.

[0085] To verify the accuracy of the range-Doppler image of this data, a range-Doppler distribution diagram of each scattering unit is plotted. Please refer to [link to relevant documentation]. Figure 4 This is the range-Doppler distribution diagram for each scattering unit. (And...) Figure 3 In contrast, the clutter distribution in the range-Doppler plane is correct.

[0086] To draw a normalized two-dimensional transmit / receive pattern of azimuth and elevation for the probe scene, please refer to [link / reference]. Figure 5 This represents the azimuth-elevation two-dimensional normalized transmit / receive pattern of the detection scenario. Since the simulated detection scenario encompasses the region including the antenna's main lobe in the elevation direction and the first side lobe in the azimuth direction, the theoretical beamwidth can be calculated as follows: The actual simulated beamwidth is The theoretical values ​​for the beam center in the elevation and azimuth directions are 6° and 80°, respectively. Figure 5 The simulation results are in agreement. Please refer to the diagram for the normalized energy distribution of clutter in each scattering element. Figure 6 This is the normalized energy distribution diagram of clutter in each scattering unit. Since the theoretical CNR of a single channel at the beam center is 16.3dB, Figure 6 The scattering element with normalized energy lower than the beam center CNR in the middle clutter wave Figure 5This is reflected in the noise level. Therefore, the simulated range-Doppler characteristics of clutter are in good agreement with the theory.

[0087] Experiment 2: In this experiment, the accuracy of the energy characteristics of the simulation results is further verified based on the clutter data obtained from Experiment 1. The generated clutter echo data is further processed by channel coherent synthesis to compensate for the phase difference caused by the time delays of different channels. After beam center spectrum shifting, the spatial characteristic spectrum of the clutter after channel coherent synthesis is plotted. Please refer to [link to relevant documentation]. Figure 7 This represents the spatial characteristic spectrum of clutter after coherent channel synthesis. Since the theoretical CNR after coherent processing of the beam center channel is... ,pass Figure 7 The CNR obtained from the statistical analysis of each feature value is as follows:

[0088]

[0089] in, express The number of large eigenvalues ​​among the eigenvalues ​​is obviously in Figure 7 middle Therefore, the calculated CNR after beam center coherent processing is 22.5 dB, which is within 0.2 dB of the theoretical CNR.

[0090] To further verify whether the channel coherent processing effect has achieved the theoretical result. At the beam center, a single moving target point is simulated using the same method as in this invention to separate it from the clutter region. An energy comparison diagram of the moving target point before and after channel coherent synthesis is plotted. Please refer to [link to relevant documentation]. Figure 8 This is a comparison of the target energy before and after channel coherent synthesis. The target energy is increased to... This matches the theoretical value.

[0091] Experiment 3: In this experiment, based on the clutter data obtained from Experiment 1, the accuracy of the spatiotemporal characteristics of the simulation results is further verified. The clutter data after pulse compression is arranged according to the following formula. Dimensional Data:

[0092] The clutter plus noise covariance matrix is ​​estimated using the following formula:

[0093] Among them, symbols This represents the conjugate transpose of a matrix.

[0094] Using the estimated clutter-noise covariance matrix, a spatiotemporal spectrum image of the clutter is plotted. Please refer to [link to relevant documentation]. Figure 9 This is the spatiotemporal spectrum image of the clutter. Calculate the Doppler frequency and spatial frequency of each scattering element, and plot the spatial frequency-Doppler curves of each scattering element on each range ring. Please refer to [link to relevant documentation]. Figure 10 The diagram shows the theoretical spatial frequency-Doppler curve of the clutter. It can be seen that the distribution of the spatiotemporal spectrum image is consistent with the distribution of the theoretical spatial frequency-Doppler curve of the clutter, proving that the spatiotemporal characteristics of the clutter are in good agreement with the theory.

[0095] In summary, the present invention provides a refined simulation method for clutter data of a missile-borne phased array multi-channel radar. Based on the characteristics of the missile-borne platform and a strict slant range history, it can simulate the radar signal transmission-reflection-reception process through a channel-by-channel, pulse-by-pulse, range-by-range loop, and scattering unit-by-scattering method. The generated clutter echo data matches the theoretical clutter characteristics, which facilitates the generation of clutter data under various radar operating conditions and detection scenarios, and provides simulation data support for further research on missile-borne STAP technology.

[0096] It should be noted that, as another aspect, this application also provides a storage medium, which may be included in an electronic device or may exist independently without being assembled into the electronic device. The aforementioned storage medium carries one or more programs, which, when executed by an electronic device, cause the electronic device to perform the methods described in the following embodiments. For example, the electronic device may perform... Figure 1 The steps of the method shown.

[0097] In one embodiment, this application provides a computer program product including a computer program that, when executed by a processor, implements the steps in the above-described method embodiments.

[0098] Furthermore, the above figures are merely illustrative of the processes included in the method according to exemplary embodiments of the present invention, and are not intended to be limiting. It is readily understood that the processes shown in the above figures do not indicate or limit the temporal order of these processes. Additionally, it is readily understood that these processes may be executed synchronously or asynchronously, for example, in multiple modules.

[0099] Other embodiments of the invention will readily occur to those skilled in the art upon consideration of the specification and practice of the invention herein. This application is intended to cover any variations, uses, or adaptations of the invention that follow the general principles of the invention and include common knowledge or customary techniques in the art not disclosed herein. The specification and embodiments are to be considered exemplary only, and the true scope and spirit of the invention are indicated by the claims.

[0100] It should be understood that the present invention is not limited to the precise structure described above and shown in the accompanying drawings, and various modifications and changes can be made without departing from its scope. The scope of the invention is defined only by the appended claims.

Claims

1. A method for fine simulation of clutter data of a missile-borne phased array multichannel radar, characterized in that, The method comprises: According to the trajectory of the projectile platform in the coherent accumulation time, the coordinates of the radar transmitter and each receiving channel at each pulse moment are obtained; The detection scene is divided into super-resolution ground points, which are successively merged according to the slant range resolution and the Doppler resolution to obtain each distance ring and the corresponding scattering unit, and the center coordinates and the scattering cross-sectional area of the scattering unit are calculated; The transmit pattern gain, the receive pattern gain and the clutter energy of each scattering unit are calculated; The transmission signal is generated in the frequency domain, the transmit-receive two-way slant range of the scattering unit is calculated based on the coordinates of the radar transmitter and each receiving channel at each pulse moment, and the frequency domain clutter echo data are simulated for each receiving channel, each pulse, each distance ring and each scattering unit based on the transmit-receive two-way slant range; the echo data of all scattering units on the same distance ring are superimposed to obtain the clutter data of each receiving channel, each pulse and each distance ring; and the time domain clutter data are obtained through inverse Fourier transform.

2. The method of claim 1, wherein, The coordinates of the radar transmitter and each receiving channel at each pulse moment are obtained according to the trajectory of the projectile platform in the coherent accumulation time, and specifically: The trajectory of the projectile platform includes its velocity. and acceleration Based on motion speed and acceleration calculate The coordinates of the missile-borne phased array multi-channel radar transmitter at each pulse moment and Coordinates of each receiving channel Specifically: in, for The transmitter coordinates at the nth pulse moment, the nth pulse moment For the first At the nth pulse moment The coordinates of each receiving channel.

3. The method of claim 1, wherein, The detection scene is divided into super-resolution ground points, which are successively merged according to the slant range resolution and the Doppler resolution to obtain each distance ring and the corresponding scattering unit, and specifically: The detection scene in the azimuth and pitch directions is divided into super-resolution ground points according to the highest beam resolution of the receiving antenna; The super-resolution ground points belonging to one slant range resolution unit are merged, and are merged into distance rings, and super-resolution ground points contained in each distance ring are obtained, and the super-resolution ground points contained in the first distance ring are denoted as ;​ For each distance ring Within a single Doppler resolution unit, super-resolution ground points are further merged, resulting in a total of [number missing]. The scattering unit is used to obtain the super-resolution ground points contained in each scattering unit within each range ring. On the distance ring, the first The set of super-resolution ground points contained in a scattering unit is denoted as ; The center coordinates and the scattering cross-sectional area of each scattering unit on each distance ring are calculated.

4. The method of claim 3, wherein, The transmit pattern gain of each scattering unit is calculated by using the following formula: wherein, denotes the transmit pattern of the m-th scatterer on the n-th distance ring, denotes the transmit pattern of the m-th scatterer on the n-th distance ring, denotes the transmit pattern of the m-th scatterer on the n-th distance ring, and denote the number of elements in the azimuth and elevation direction, respectively, and denote the row and column subarray weighting coefficients, respectively, denotes the element spacing, denotes the signal wavelength, and denote the azimuth and elevation angle of the m-th scatterer on the n-th distance ring, denote the azimuth and elevation angle of the m-th scatterer on the n-th distance ring, denote the azimuth and elevation angle of the m-th scatterer on the n-th distance ring, and denote the azimuth and elevation angle of the beam center.

5. The method of claim 4, wherein, The receive pattern gain of each scattering unit is calculated by using the following formula: wherein denotes the th receive pattern of the th scattering element on the th distance ring.

6. The method of claim 5, wherein, The clutter energy of each scattering unit is calculated by using the following formula: where represents the clutter energy of the jth scatterer on the ith range bin, represents the average transmitted power, and represents the transmit and receive antenna gains, respectively, represents the pulse repetition period, represents the system loss, represents the signal pulse width.​​ 7. The method of claim 1, wherein, The frequency domain clutter echo data are simulated for each receiving channel, each pulse, each distance ring and each scattering unit by using the following formula: wherein, and denote the inter-channel amplitude error and the phase error of the channel, denotes the signal frequency, denotes the transceiver two-way slant range of the scattering element, denotes the th received channel, the th pulse, the th range bin, the th frequency domain echo signal of the scattering element.

8. A storage medium having stored thereon a computer program, characterized in that The computer program is executed by the processor to realize the missile-borne phased array multi-channel radar clutter data fine simulation method according to any one of claims 1 to 7.

9. A computer program product comprising a computer program, characterized in that, The computer program is executed by the processor to realize the missile-borne phased array multi-channel radar clutter data fine simulation method according to any one of claims 1 to 7.

10. An electronic device, comprising: Comprise: A processor; And A memory for storing executable instructions of the processor; Wherein the processor is configured to execute the missile-borne phased array multi-channel radar clutter data fine simulation method according to any one of claims 1 to 7 by executing the executable instructions.

Images