A Fast Simulation Method and Apparatus for Sea Clutter Based on Backscattering Morchin Model

Abstract

This invention relates to a method and apparatus for rapid simulation of sea clutter based on the backscattering Morchin model. The method includes: determining a radar illumination area; dividing the radar illumination area into multiple scattering units according to the range ring resolution width and azimuth resolution; dividing each scattering unit into two triangular-like regions, and calculating the scattering cross-section of each triangular-like region using the backscattering Morchin model; obtaining the echo signal of the corresponding scattering unit based on the scattering cross-section of each triangular-like region; and superimposing the echo signals of scattering units on the same range ring to obtain the simulated sea clutter echo signal. This invention solves the problem of rapid sea clutter simulation. The generated simulated sea clutter echo signal exhibits strong temporal correlation, improving the effectiveness and realism of sea clutter simulation, and demonstrating high efficiency for rapid sea clutter simulation.

Description

Technical Field

[0001] This invention relates to the field of radar technology, and more specifically, to a method and apparatus for rapid simulation of sea clutter based on the backscattering Morchin model. Background Technology

[0002] With the development of radar technology and the increasing demand for marine surveillance and detection missions, sea clutter simulation technology has become an important application area in radar systems. For marine radars, the sea environment is complex and variable; waves, wind speed, and other factors can cause strong sea clutter interference, affecting the accurate detection and tracking of targets. While traditional radars can complete basic detection tasks in simple sea conditions, sea clutter simulation technology is needed to optimize radar detection performance in complex sea conditions to improve the accuracy and efficiency of target identification. Summary of the Invention

[0003] To address the aforementioned problems in the existing technology, this invention provides a method and apparatus for rapid simulation of sea clutter based on the backscattering Morchin model.

[0004] According to a first aspect of the present invention, a fast simulation method for sea clutter based on the backscattering Morchin model is provided, the method comprising:

[0005] Determine the radar illumination area;

[0006] The radar illumination area is divided into multiple scattering units based on the range ring resolution width and azimuth resolution.

[0007] Each scattering unit is divided into two triangular regions, and the scattering cross-section of each triangular region is calculated using the backscattering Morchin model.

[0008] The echo signal of the corresponding scattering unit is obtained based on the scattering cross-sectional area of ​​each triangular region.

[0009] By superimposing the echo signals of scattering units located on the same distance ring, a simulated sea clutter echo signal is obtained.

[0010] Optionally, the scattering cross-section of each triangular region is represented as follows:

[0011]

[0012] Where, σ m,k Let represent the scattering cross-section of the triangular region of the scattering unit at the k-th azimuth resolution of the m-th distance ring. dA represents the local backscattering coefficient of the triangular region corresponding to the scattering surface at the k-th azimuth resolution of the scattering unit in the m-th range ring. m,k ζ represents the area of ​​the triangular region of the scattering unit at the k-th azimuth resolution of the m-th range ring. m,k This represents the angle between the triangular region of the scattering unit at the k-th azimuth resolution of the m-th distance ring and the standard horizontal plane.

[0013] Alternatively, the echo signal of the scattering unit is represented as follows:

[0014] S(θ i,k ,R i,m )=S1(θ i,k ,R i,m )+S2(θ i,k ,R i,m );

[0015] Where S(θ) i,k ,R i,m R represents the echo signal received by the radar at time i on the scattering element. i,m θ represents the distance from the origin of the rectangular coordinate system to the m-th scattering unit on the distance ring at time i. i,k S1(θ) represents the angle between the scattering element at the k-th azimuth resolution at time i and the x-axis. i,k ,R i,m S2(θ) represents the echo signal received by the radar at time i on the first type of triangular region. i,k ,R i,m ) represents the echo signal received by the radar at time i on the second type of triangular region.

[0016] Optionally, the echo signal received by the radar at time i on the triangular region is represented as follows:

[0017]

[0018] Where, σ m,k This represents the scattering cross-section of the triangular region of the scattering unit at the k-th azimuth resolution of the m-th range ring, where j represents the imaginary unit. The triangular scattering unit echo phase represents the scattering unit at the k-th azimuth angle resolution of the m-th range ring.

[0019] Optionally, the simulated sea clutter echo signal is represented as follows:

[0020]

[0021] in, P represents the simulated sea clutter echo signal corresponding to the m-th range loop at time i. cTo simulate the average power required for sea clutter, P sc S represents the sea clutter power in the radar illumination area. i,m This represents the initial echo signal corresponding to the m-th distance loop at time i.

[0022] Optionally, the initial echo signal corresponding to the m-th distance loop at time i is represented as follows:

[0023]

[0024] Where, N θ For a single range ring containing azimuth resolution, S(θ) i,k ,R i,m R represents the echo signal received by the radar at time i on the scattering element. i,m θ represents the distance from the origin of the rectangular coordinate system to the m-th scattering unit on the distance ring at time i. i,k This represents the angle between the scattering element at the k-th azimuth resolution at time i and the x-axis.

[0025] Optionally, the sea clutter power of the radar illumination area is expressed as follows:

[0026]

[0027] Where, N p N represents the radar pulse number. R This indicates the number of range rings in the radar illumination area.

[0028] According to a second aspect of the present invention, a rapid simulation device for sea clutter based on a backscattering Morchin model is provided, the device comprising:

[0029] The area determination module is used to determine the radar illumination area;

[0030] The region division module is used to divide the radar illumination area into multiple scattering units according to the range ring resolution width and azimuth resolution.

[0031] The scattering cross-section calculation module is used to divide each scattering unit into two triangular regions and calculate the scattering cross-section of each triangular region using the backscattering Morchin model.

[0032] The echo signal acquisition module is used to obtain the echo signal of the corresponding scattering unit based on the scattering cross-sectional area of ​​each triangular region.

[0033] The simulation module is used to superimpose the echo signals of scattering units on the same distance ring to obtain simulated sea clutter echo signals.

[0034] The technical solution provided by this invention may include the following beneficial effects:

[0035] Based on the above technical solution, a fast simulation method for sea clutter based on the backscattering Morchin model is proposed. This invention solves the problem of fast simulation of sea clutter. The generated simulated sea clutter echo signal has strong temporal correlation, which improves the effectiveness and realism of sea clutter simulation. It has high efficiency and can perform fast simulation of sea clutter.

[0036] Other features and advantages of the present invention will be described in detail in the following detailed description section. Attached Figure Description

[0037] The accompanying drawings are provided to further illustrate the invention and form part of the specification. They are used together with the following detailed description to explain the invention, but do not constitute a limitation thereof. In the drawings:

[0038] Figure 1 This is a flowchart illustrating a rapid simulation method for sea clutter based on the backscattering Morchin model, according to an exemplary embodiment.

[0039] Figure 2 This is a schematic diagram illustrating a radar illumination area according to an exemplary embodiment.

[0040] Figure 3 This is a schematic diagram of a scattering unit according to an exemplary embodiment.

[0041] Figure 4 This is a schematic diagram of yet another scattering unit according to an exemplary embodiment.

[0042] Figure 5 This is a schematic diagram illustrating the simulation results of sea level with a wind speed of 6 m / s, according to an exemplary embodiment.

[0043] Figure 6 This is a schematic diagram illustrating the simulation results of a simulated sea clutter echo signal according to an exemplary embodiment.

[0044] Figure 7 This is an amplitude distribution histogram of a simulated sea clutter echo signal, illustrated according to an exemplary embodiment.

[0045] Figure 8a This is a schematic diagram illustrating the comparison between a simulated sea clutter echo signal and a classical statistical model, according to an exemplary embodiment.

[0046] Figure 8b This is a schematic diagram illustrating the comparison between a simulated sea clutter echo signal and a classical statistical model, according to an exemplary embodiment.

[0047] Figure 9a This is a schematic diagram of a sea clutter spectral profile according to an exemplary embodiment.

[0048] Figure 9b This is a schematic diagram illustrating an average spectral profile according to an exemplary embodiment.

[0049] Figure 10a This is a schematic diagram illustrating the time autocorrelation coefficient variation curve of a simulated sea clutter echo signal according to an exemplary embodiment.

[0050] Figure 10b This is a schematic diagram illustrating the range autocorrelation coefficient variation curve of a simulated sea clutter echo signal according to an exemplary embodiment.

[0051] Figure 11 This is a block diagram illustrating a rapid simulation device for sea clutter based on the backscattering Morchin model, according to an exemplary embodiment. Detailed Implementation

[0052] Figure 1 This is a flowchart illustrating a rapid simulation method for sea clutter based on the backscattering Morchin model, according to an exemplary embodiment. Figure 1 As shown, the method includes the following steps.

[0053] S101. Determine the radar illumination area.

[0054] Understandably, before determining the radar illumination area, a sea level model can be established. In a rectangular coordinate system, the coordinates of the sea level scatterer at a certain moment can be represented as:

[0055]

[0056] Where g represents the gravitational acceleration at the sea surface, the y-coordinate is fixed, and x0 is the coordinate of x at the initial moment; γ j They follow an independent distribution between 0 and 2π; It follows a Gaussian distribution with a standard deviation of 0.2 rad; α j Let ω be a Gaussian random variable reflecting the height of ocean waves. j These are the wave spectrum distribution parameters.

[0057] In one implementation, Figure 2 This is a schematic diagram illustrating a radar illumination area according to an exemplary embodiment, such as... Figure 2 As shown, assuming that the beam illumination forms an elliptical coverage area under down-looking radar mode, the area of ​​the radar antenna's coverage area is:

[0058] A c =ΔRθ 3dBRsec(ψ g );

[0059] Where ΔR is the distance resolution, ΔR = c / 2B, c is the speed of light, B is the bandwidth, and θ 3dB R is the half-power beamwidth of the radar antenna, R is the radar detection range, and ψ is the half-power beamwidth of the radar antenna. g For secting, sec() represents a secant.

[0060] S102. Divide the radar illumination area into multiple scattering units according to the range ring resolution width and azimuth resolution.

[0061] Understandable Figure 3 This is a schematic diagram of a scattering unit according to an exemplary embodiment, such as... Figure 3 As shown, a coordinate system is established with the sea surface projection point at the radar's initial location as the origin and the direction of wave motion as the x-axis. The range ring resolution width is set to ΔR, and the azimuth resolution size is Δθ. The radar's main beam illumination area is divided into several scattering units of ΔR × Δθ. The number of range rings is N. R A single range ring contains N azimuth resolution units. θ .

[0062] S103. Divide each scattering unit into two triangular regions and calculate the scattering cross-section of each triangular region using the backscattering Morchin model.

[0063] Understandable Figure 4 This is a schematic diagram illustrating yet another scattering unit according to an exemplary embodiment, such as... Figure 4 As shown, Figure 3 The scattering unit can be partially magnified to obtain Figure 4 Each scattering unit is divided into two triangular-like regions, and the scattering cross-section of each triangular-like region is calculated based on the backscattering Morchin model. Taking time i as an example, A m,k dA represents the spatial scattering unit at the k-th angular resolution of the m-th range ring. m,k For A m,k The projected area at the standard level; R i,m Let θ be the distance from the scattering unit on the m-th distance ring to the origin of the rectangular coordinate system. i,k Let be the angle between the k-th azimuth sector scattering unit and the x-axis. Each scattering unit is arranged according to... Figure 4 The form shown will dA m,k It is divided into two triangular-like forms. The areas of the two triangular-like scattering units satisfy:

[0064]

[0065] Let the direction vector of the radar beam at time i be... The radar coordinates are (x p ,y p ,z p ),Will Figure 4 The cylindrical coordinates of a, b, c, and d, when converted to rectangular coordinates, are (x, y, c, d) respectively. a ,y a ,z a ), (x b ,y b ,z b ), (x c ,y c ,z c ), (x d ,y d ,z d A m,k,1 The triangle formed by points a, b, and c has a normal vector. and

[0066] The local ground grazing angle of the radar main beam on this triangular scattering surface can be expressed as:

[0067]

[0068] The sea clutter backscattering coefficient σ in the Morchin model 0 It can be represented as:

[0069]

[0070] Where λ is the radar wavelength, ψ g For the angle of the brush, S represents the sea state class, ψ c This is the wiping angle factor.

[0071] Scrubbing angle factor ψ c The definition is expressed as follows:

[0072]

[0073] Among them, h e The roughness of the sea surface (m) is defined by the following formula:

[0074]

[0075] The scattering cross-section of each triangular region is represented as follows:

[0076]

[0077] Where, σ m,kLet represent the scattering cross-section of the triangular region of the scattering unit at the k-th azimuth resolution of the m-th distance ring. dA represents the local backscattering coefficient of the triangular region corresponding to the scattering surface at the k-th azimuth resolution of the scattering unit in the m-th range ring. m,k ζ represents the area of ​​the triangular region of the scattering unit at the k-th azimuth resolution of the m-th range ring. m,k This represents the angle between the triangular region of the scattering unit at the k-th azimuth resolution of the m-th distance ring and the standard horizontal plane.

[0078] Using the first type of triangular region dA m,k,1 For example, consider the local rubbing angle ψ. m,k,1 Substitute into the Morchin model to calculate the backscattering coefficient Let the normal vector of the standard horizontal plane be... The angle ζ between the scattering unit region and the standard horizontal plane can be calculated. m,k,1 It can be represented as follows:

[0079]

[0080] The obtained ζ m,k,1 Substituting into the formula for the scattering cross-section of a triangular region, we can obtain the scattering cross-section σ of the first type of triangular region. m,k,1 .

[0081] S104. Obtain the echo signal of the corresponding scattering unit based on the scattering cross-sectional area of ​​each triangular region.

[0082] Understandably, firstly, according to the formula for calculating the centroid of a triangle, A... m,k,1 Triangular region and A m,k,2 The centroids of the triangular regions are as follows:

[0083]

[0084] The distances from the radar to the two triangular-shaped scattering elements are as follows:

[0085]

[0086] The echo phase of each triangular scattering unit is:

[0087]

[0088] Based on the echo phase of each triangular scattering unit, the echo signal received by the radar at time i in the triangular region can be represented as follows:

[0089]

[0090] Where, σ m,k This represents the scattering cross-section of the triangular region of the scattering unit at the k-th azimuth resolution of the m-th range ring, where j represents the imaginary unit. The triangular scattering unit echo phase represents the scattering unit at the k-th azimuth angle resolution of the m-th range ring.

[0091] Furthermore, the echo signal of the scattering unit is represented as follows:

[0092] S(θ i,k ,R i,m )=S1(θ i,k ,R i,m )+S2(θ i,k ,R i,m );

[0093] Where S(θ) i,k ,R i,m R represents the echo signal received by the radar at time i on the scattering element. i,m θ represents the distance from the origin of the rectangular coordinate system to the m-th scattering unit on the distance ring at time i. i,k S1(θ) represents the angle between the scattering element at the k-th azimuth resolution at time i and the x-axis. i,k ,R i,m S2(θ) represents the echo signal received by the radar at time i on the first type of triangular region. i,k ,R i,m ) represents the echo signal received by the radar at time i on the second type of triangular region.

[0094] S105. Superimpose the echo signals of the scattering units on the same distance ring to obtain the simulated sea clutter echo signal.

[0095] It is understandable that by superimposing the vectors of the scattering elements within the same range ring of the radar, the initial echo signal corresponding to the m-th range ring at time i is represented as follows:

[0096]

[0097] Where, N θ For a single range ring containing azimuth resolution, S(θ) i,k ,R i,m R represents the echo signal received by the radar at time i on the scattering element. i,m θ represents the distance from the origin of the rectangular coordinate system to the m-th scattering unit on the distance ring at time i. i,k This represents the angle between the scattering element at the k-th azimuth resolution at time i and the x-axis.

[0098] Let the number of pulses be N. pThe power of sea clutter in the radar illumination area is expressed as follows:

[0099]

[0100] Where, N p N represents the radar pulse number. R This indicates the range ring number within the radar illumination area.

[0101] Based on the sea clutter power in the radar-illuminated area and the initial echo signal corresponding to the m-th range loop at time i, the simulated sea clutter echo signal can be represented as follows:

[0102]

[0103] in, P represents the simulated sea clutter echo signal corresponding to the m-th range loop at time i. c To simulate the average power required for sea clutter, P sc S represents the sea clutter power in the radar illumination area. i,m This represents the initial echo signal corresponding to the m-th distance loop at time i.

[0104] In one implementation, to verify the effectiveness and superiority of this fast sea clutter simulation method, the method is used to simulate sea clutter signals. Specific simulation parameters are shown in Table 1:

[0105] Table 1

[0106] Parameter name Parameter settings Radar height / m 100 Scrub corner / ° 5 Pulse width / μs 2 Pulse count 100 Bandwidth / MHz 1 Sampling frequency / MHz 10 Carrier frequency / GHz 2 Beamwidth / ° 3.3 Azimuth resolution / ° 0.05 polarization mode HH Pulse repetition period / μs 480 sea ​​feeling 3 Wind speed / m / s 6 Spurious-to-noise ratio / dB 40

[0107] Figure 5 This is a schematic diagram illustrating the sea level simulation results with a wind speed of 6 m / s, based on an exemplary embodiment. Figure 6 This is a schematic diagram illustrating the simulation results of a simulated sea clutter echo signal according to an exemplary embodiment. The sea clutter signal is located at distance cells 50 to 100. Figure 7 This is an amplitude distribution histogram of a simulated sea clutter echo signal, illustrated according to an exemplary embodiment.

[0108] Statistical characteristics of simulated sea clutter signals were analyzed, and the parameters of Rayleigh, Log-Normal, Weibull and K distributions were estimated using the maximum likelihood estimation method. The parameter estimation results are shown in Table 2.

[0109] Table 2

[0110]

[0111] Substituting the parameters into the statistical model yields a comparative diagram. Figure 8aThis is a schematic diagram illustrating a comparison between a simulated sea clutter echo signal and a classical statistical model, according to an exemplary embodiment. Figure 8b This is a schematic diagram illustrating a comparison between a simulated sea clutter echo signal and a classical statistical model, according to an exemplary embodiment. Figure 8a This is a schematic diagram comparing simulated sea clutter echo signals with Rayleigh and Log-Normal distributions. Figure 8b This is a schematic diagram comparing the simulated sea clutter echo signal with the Weibull and K distributions.

[0112] The MSD test was performed on the fitting results of the four statistical models, and the test results are shown in Table 3.

[0113] Table 3

[0114] Statistical distribution model MSD test value Rayleigh distribution 0.0304 Log-Normal distribution 0.0400 Weibull distribution 0.0270 K-distribution 0.0231

[0115] The MSD goodness-of-fit test method indicates that a smaller MSD value indicates a better fit. As shown in Table 3, it can be seen that the simulated sea clutter signal fits the K-distribution statistical model best.

[0116] Coherent accumulation was used to perform spectral analysis on the simulated sea clutter signal. The sea clutter signal had a total of 100 pulses, and the FFT points were set to 256. Figure 9a This is a schematic diagram illustrating a sea clutter spectrum according to an exemplary embodiment. Figure 9b This is a schematic diagram of an average spectrum according to an exemplary embodiment, showing the sea clutter spectrum and an average spectral profile of the sea clutter region (50-100 distance cells). Figure 9a and Figure 9b As shown in the figure. Figure 9a and Figure 9b It can be seen that the sea clutter Doppler distribution is around 29.41 Hz, and the clutter spectrum is broadened. This is because the wave motion causes the clutter spectrum to broaden and the clutter Doppler center is not zero, which is consistent with the actual situation.

[0117] The correlation between the time and range dimensions of simulated sea clutter signals was analyzed. Figure 10a This is a schematic diagram illustrating the time autocorrelation coefficient variation curve of a simulated sea clutter echo signal according to an exemplary embodiment. Figure 10b This is a schematic diagram illustrating the range autocorrelation coefficient variation curve of a simulated sea clutter echo signal according to an exemplary embodiment. It is generally assumed that if the correlation coefficient is greater than e... -1 The system assumes that signals are correlated, and vice versa. Among these, Figure 10aThe different colored curves represent the time autocorrelation coefficient variation curves at different distance units, and the thick black line represents the variation curve of the mean time autocorrelation coefficient over time; it can be seen that the decorrelation time of the simulated sea clutter echo signal is about 12ms. Figure 10b In the diagram, one color represents the distance autocorrelation coefficient at a pulse location as a function of distance cells, and the thick black line represents the mean distance autocorrelation coefficient as a function of distance cells. It can be seen that the mean distance autocorrelation coefficient of sea clutter decreases to e after one distance cell. -1 The value quickly stabilizes below 0.3, indicating that the simulated sea clutter echo signal is random in the distance dimension.

[0118] The simulation results above show that the simulated sea clutter echo signal generated by this invention conforms to the K-distribution statistical model; and the sea clutter has a relatively strong correlation in time, proving the effectiveness and superiority of the fast sea clutter simulation method.

[0119] Traditional sea clutter simulation methods, as described above, have several limitations in terms of efficiency and accuracy. For example, spectral models require complex marine environmental parameters, which are time-consuming to acquire and process, and struggle to capture subtle features under complex sea conditions. Statistical models rely on large amounts of historical data, have poor generalization ability, and are ill-suited for handling unseen extreme sea conditions. Numerical simulations are computationally intensive and difficult to apply in real-time. Neural network simulations are prone to overfitting when training data is insufficient, resulting in unstable performance in real-world environments. In contrast, this invention enables rapid simulation, and the generated sea clutter statistical features conform to the K-distribution model, exhibiting strong time correlation. This significantly improves the effectiveness and realism of the simulation, enhancing its speed and accuracy, and providing more reliable technical support for radar detection in marine environments.

[0120] Figure 11 This is a block diagram illustrating a rapid simulation device for sea clutter based on the backscattering Morchin model, according to an exemplary embodiment. See also... Figure 11 The device 1100 may include:

[0121] Area determination module 1101 is used to determine the radar illumination area;

[0122] The region division module 1102 is used to divide the radar illumination area into multiple scattering units according to the range ring resolution width and azimuth resolution.

[0123] The scattering cross-section calculation module 1103 is used to divide each scattering unit into two triangular regions and calculate the scattering cross-section of each triangular region using the backscattering Morchin model.

[0124] The echo signal acquisition module 1104 is used to obtain the echo signal of the corresponding scattering unit based on the scattering cross-sectional area of ​​each triangular region.

[0125] Simulation module 1105 is used to superimpose the echo signals of scattering units on the same distance ring to obtain simulated sea clutter echo signals.

[0126] Regarding the apparatus in the above embodiments, the specific manner in which each module performs its operation has been described in detail in the embodiments related to the method, and will not be elaborated upon here.

[0127] The preferred embodiments of the present invention have been described in detail above with reference to the accompanying drawings. However, the present invention is not limited to the specific details of the above embodiments. Within the scope of the technical concept of the present invention, various simple modifications can be made to the technical solution of the present invention, and these simple modifications all fall within the protection scope of the present invention.

[0128] It should also be noted that the various specific technical features described in the above specific embodiments can be combined in any suitable manner without contradiction. In order to avoid unnecessary repetition, the present invention will not describe the various possible combinations separately.

[0129] Furthermore, various different embodiments of the present invention can be combined in any way, as long as they do not violate the spirit of the present invention, they should also be regarded as the content disclosed by the present invention.

Claims

1. A fast simulation method for sea clutter based on the backscattering Morchin model, characterized in that, The method includes: Determine the radar illumination area; The radar illumination area is divided into multiple scattering units based on the range ring resolution width and azimuth resolution. Each scattering unit is divided into two triangular regions, and the scattering cross-section of each triangular region is calculated using the backscattering Morchin model. The echo signal of the corresponding scattering unit is obtained based on the scattering cross-sectional area of ​​each triangular region. By superimposing the echo signals of scattering units located on the same distance ring, a simulated sea clutter echo signal is obtained.

2. The fast simulation method for sea clutter based on the backscattering Morchin model according to claim 1, characterized in that, The scattering cross-sectional area of ​​each triangular region is represented as follows: Where, σ m,k Let represent the scattering cross-section of the triangular region of the scattering unit at the k-th azimuth resolution of the m-th distance ring. dA represents the local backscattering coefficient of the triangular region corresponding to the scattering surface at the k-th azimuth resolution of the scattering unit in the m-th range ring. m,k ζ represents the area of ​​the triangular region of the scattering unit at the k-th azimuth resolution of the m-th range ring. m,k This represents the angle between the triangular region of the scattering unit at the k-th azimuth resolution of the m-th distance ring and the standard horizontal plane.

3. The fast simulation method for sea clutter based on the backscattering Morchin model according to claim 1, characterized in that, The echo signal of the scattering unit is represented as follows: S(θ i,k ,R i,m )=S1(θ i,k ,R i,m )+S2(θ i,k ,R i,m ); Where S(θ) i,k ,R i,m R represents the echo signal received by the radar at time i on the scattering element. i,m θ represents the distance from the origin of the rectangular coordinate system to the m-th scattering unit on the distance ring at time i. i,k S1(θ) represents the angle between the scattering element at the k-th azimuth resolution at time i and the x-axis. i,k ,R i,m S2(θ) represents the echo signal received by the radar at time i on the first type of triangular region. i,k ,R i,m ) represents the echo signal received by the radar at time i on the second type of triangular region.

4. The fast simulation method for sea clutter based on the backscattering Morchin model according to claim 3, characterized in that, The echo signal received by the radar at time i on the triangular region is represented as follows: Where, σ m,k This represents the scattering cross-section of the triangular region of the scattering unit at the k-th azimuth resolution of the m-th range ring, where j represents the imaginary unit. The triangular scattering unit echo phase represents the scattering unit at the k-th azimuth angle resolution of the m-th range ring.

5. The fast simulation method for sea clutter based on the backscattering Morchin model according to claim 1, characterized in that, The simulated sea clutter echo signal is represented as follows: in, P represents the simulated sea clutter echo signal corresponding to the m-th range loop at time i. c To simulate the average power required for sea clutter, P sc S represents the sea clutter power in the radar illumination area. i,m This represents the initial echo signal corresponding to the m-th distance loop at time i.

6. The fast simulation method for sea clutter based on the backscattering Morchin model according to claim 5, characterized in that, The initial echo signal corresponding to the m-th distance loop at time i is represented as follows: Where, N θ For a single range ring containing azimuth resolution, S(θ) i,k ,R i,m R represents the echo signal received by the radar at time i on the scattering element. i,m θ represents the distance from the origin of the rectangular coordinate system to the m-th scattering unit on the distance ring at time i. i,k This represents the angle between the scattering element at the k-th azimuth resolution at time i and the x-axis.

7. The fast simulation method for sea clutter based on the backscattering Morchin model according to claim 6, characterized in that, The power of sea clutter in the radar-illuminated area is expressed as follows: Where, N p N represents the radar pulse number. R This indicates the number of range rings in the radar illumination area.

8. A rapid simulation device for sea clutter based on the backscattering Morchin model, characterized in that, The device includes: The area determination module is used to determine the radar illumination area; The region division module is used to divide the radar illumination area into multiple scattering units according to the range ring resolution width and azimuth resolution. The scattering cross-section calculation module is used to divide each scattering unit into two triangular regions and calculate the scattering cross-section of each triangular region using the backscattering Morchin model. The echo signal acquisition module is used to obtain the echo signal of the corresponding scattering unit based on the scattering cross-sectional area of ​​each triangular region. The simulation module is used to superimpose the echo signals of scattering units on the same distance ring to obtain simulated sea clutter echo signals.

Images