Missile-borne phased array multi-channel radar clutter data refinement simulation method

Abstract

The application particularly relates to a missile-borne phased array multichannel radar clutter data refinement simulation method, which comprises the following steps: obtaining the coordinates of a radar transmitter and each receiving channel at each pulse moment according to the trajectory of a missile platform within a coherent accumulation time; dividing a detection scene into super-resolution ground points, and sequentially performing merging according to slant range resolution and Doppler resolution to obtain each distance ring and its corresponding scattering unit; then, calculating the directional diagram gain and clutter energy of each scattering unit; finally, based on a strict slant range history, performing four-level element-by-element simulation in a "channel-pulse-distance ring-scattering unit" mode, and superimposing the echoes of the scattering units of the same distance ring to obtain the 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 radar working condition can be obtained, and simulation data support is provided for further research on missile-borne STAP technology.

Description

Technical Field

[0001] This invention relates to the field of signal processing technology, specifically to a refined simulation method for clutter data of missile-borne phased array multi-channel radar, which can be used in missile-borne phased array multi-channel radar. Background Technology

[0002] Space-time adaptive processing (STAP) is an effective technique widely used in airborne / spaceborne early warning radar for moving target detection. Based on array multi-channel theory, this technique adaptively suppresses clutter by combining spatial and temporal data. However, STAP technology has not been widely adopted in missile-borne radars due to limitations imposed by platform size and high maneuverability. Traditional missile-borne radars typically employ non-adaptive Doppler clutter suppression methods, such as Moving Target Detection (MTD) and Pulse Doppler (PD), which are insufficient to meet the increasingly urgent need for detecting weak, slow-moving targets in strong clutter environments. 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 researchers have conducted research on the clutter characteristics of airborne phased array multi-channel radars in non-frontal arrays and different motion states, and proposed a series of improved STAP methods suitable for airborne platforms, such as registration compensation methods based on radar parameters and robust sum-difference STAP methods.

[0004] However, the clutter signal model used in current research is a space-time steering vector model, which does not consider the effects of ground scattering unit resolution and high-speed projectile motion, and the characteristics of the obtained clutter data may differ from the actual echo. Obtaining measured clutter data from missile-borne multi-channel radar requires significant manpower and material resources. In order to quickly evaluate the performance of STAP technology under different operating conditions during the early engineering demonstration stage, obtaining high-fidelity clutter data through modeling and simulation is an efficient method. Therefore, it is urgent to establish a complete and accurate simulation method for clutter data from missile-borne phased array multi-channel radar.

[0005] It should be noted that the information disclosed in the background section above is only used to enhance the understanding of the background of the present invention, and therefore may include information that does not constitute prior art known to those skilled in the art. Summary of the Invention

[0006] To address the aforementioned problems in existing technologies, this invention provides a refined simulation method for clutter data from a missile-borne phased array multi-channel radar. This simulation method combines the characteristics of missile-borne radar, is based on a rigorous slant range history, and considers the impact of actual resolution. It simulates echo data sequentially by channel, pulse, range loop, and scattering unit, rigorously simulating the radar signal transmission-reflection-reception process, and can obtain high-fidelity clutter echo data under given detection scenarios and radar operating conditions.

[0007] Other features and advantages of the invention will become apparent from the following detailed description, or may be learned in part by practice of the invention.

[0008] According to a first aspect of the present invention, a method for refined simulation of clutter data from a missile-borne phased array multi-channel radar is provided, the method comprising: Based on the trajectory of the projectile platform during 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 then merged sequentially according to slant range resolution and Doppler resolution to obtain each range ring and its corresponding scattering unit, and the center coordinates and scattering cross-section of the scattering unit are calculated. Calculate the transmit pattern gain, receive pattern gain, and clutter energy for each scattering element; The transmitted signal is generated in the frequency domain. The slant range of the scattering unit is calculated based on the coordinates of the radar transmitter and the receiving channel at each pulse moment. Based on the slant range, frequency domain clutter echo data is simulated for each receiving channel, each pulse, each range loop, and each scattering unit. The echo data of all scattering units on the same range loop are superimposed to obtain the clutter data for each pulse and each range loop under each receiving channel. The inverse Fourier transform is then performed to obtain the time domain clutter data.

[0009] In some exemplary embodiments, obtaining the coordinates of the radar transmitter and each receiving channel at each pulse moment based on the trajectory of the projectile platform during the coherent accumulation time specifically involves: 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:

[0010]

[0011] 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.

[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 and 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 invention, a computer program product is provided, on which a computer program is stored, wherein when the computer program is executed by a processor, the fine simulation method for airborne phased array multi-channel radar clutter data described in the first aspect is implemented.

[0023] According to a fourth aspect of the present invention, an electronic device is provided, comprising: Processor; and Memory for storing the executable instructions of the processor; The processor is configured to implement the above-described method for fine simulation of missile-borne phased array multi-channel radar clutter data when executing the executable instructions.

[0024] The refined simulation method for airborne phased array multi-channel radar clutter data provided in the embodiments of the present invention has the following advantages compared with the prior art: 1. This invention combines the characteristics of missile-borne radar, is based on a strict slant range history, and takes into account the impact of actual resolution. By simulating echo data channel by channel, pulse by pulse, range loop by range, and scattering unit by scattering unit, it can strictly simulate the radar signal transmission-reflection-reception process, making up for the limitations of existing simulations that are based on space-time steering vectors and do not consider two-dimensional resolution.

[0025] 2. The clutter energy, clutter Doppler spectral width, clutter multi-channel synthesis effect, and clutter space-time characteristics obtained by simulation are in good agreement with the theoretical results. It can accurately and effectively obtain high-fidelity clutter echo data for a given scenario, providing simulation data support for further evaluation of missile-borne phased array multi-channel radar detection capabilities and STAP technology research.

[0026] It should be understood that the above general description and the following detailed description are exemplary and explanatory only, and are not intended to limit the invention. Attached Figure Description

[0027] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the invention and, together with the description, serve to explain the principles of the invention. It is obvious that the drawings described below are merely some embodiments of the invention, and those skilled in the art can obtain other drawings based on these drawings without any inventive effort.

[0028] Figure 1 This is a flowchart of a method for refining the simulation of clutter data from a missile-borne phased array multi-channel radar according to the present invention; 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; 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 composed of equivalent channels, at the initial time... 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. Divide the distance loops to obtain the super-resolution ground points contained in each distance loop.

[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 super-resolution ground point indices 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 would require the full context.] 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 Boltzmann's 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 refined simulation of clutter data from a missile-borne phased array multi-channel radar, characterized in that, The method includes: Based on the trajectory of the projectile platform within the coherent accumulation time, the coordinates of the radar transmitter and each receiving channel at each pulse moment are obtained; 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 The pulse moment The coordinates of each receiving channel; The detection scene is divided into super-resolution ground points, which are then merged sequentially according to slant range resolution and Doppler resolution to obtain each range ring and its corresponding scattering unit, and the center coordinates and scattering cross-section of the scattering unit are calculated. Calculate the transmit pattern gain, receive pattern gain, and clutter energy for each scattering element; The transmitted signal is generated in the frequency domain. The transmit and receive slant range of the scattering unit is calculated based on the coordinates of the radar transmitter and each receiving channel at each pulse moment. Based on the transmit and receive slant range, frequency domain clutter echo data is simulated for each receiving channel, each pulse, each range loop, and each scattering unit. The echo data of all scattering units on the same range loop are superimposed to obtain the clutter data for each pulse and each range loop under each receiving channel. The inverse Fourier transform is then performed to obtain the time domain clutter data. 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: in, and They represent the first Inter-channel amplitude and phase errors of each channel Indicates signal frequency. This represents the transmit and 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.

2. The method according to claim 1, characterized in that, The process 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 is as follows: 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.

3. The method according to claim 2, characterized in that, The emission pattern gain of each scattering unit is calculated using the following formula: 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.

4. The method according to claim 3, characterized in that, The receiver pattern gain of each scattering element is calculated using the following formula: in, Indicates the first On the distance ring, the first The receiving pattern of each scattering unit.

5. The method according to claim 4, characterized in that, The clutter energy of each scattering unit is calculated using the following formula: 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.

6. A storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the refined simulation method for missile-borne phased array multi-channel radar clutter data as described in any one of claims 1 to 5.

7. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by the processor, it implements the refined simulation method for multi-channel radar clutter data of missile-borne phased array as described in any one of claims 1 to 5.

8. An electronic device, characterized in that, include: processor; as well as Memory for storing the executable instructions of the processor; The processor is configured to execute the airborne phased array multi-channel radar clutter data refinement simulation method according to any one of claims 1 to 5 by executing the executable instructions.

Images