Sea clutter rapid simulation method and device based on back scattering Morchin model

Abstract

The invention relates to a sea clutter rapid simulation method and device based on a back scattering Morchin model, and the method is characterized in that the method comprises the steps: determining a radar irradiation region; dividing a radar irradiation area into a plurality of scattering units according to the distance ring resolution width and the azimuth angle resolution; dividing each scattering unit into two triangle-like regions, and calculating the scattering sectional area of each triangle-like region through a back scattering Morchin model; obtaining an echo signal of the corresponding scattering unit according to the scattering sectional area of each triangle-like area; and superposing the echo signals of the scattering units on the same distance ring to obtain a simulated sea clutter echo signal. According to the method, the problem of rapid simulation of the sea clutter is solved, the generated simulated sea clutter echo signal has high correlation in time, the effectiveness and authenticity of sea clutter simulation are improved, the efficiency is high, and the sea clutter can be rapidly simulated.

Description

Technical Field

[0001] The invention relates to the field of radar technology, and in particular to a method and a device for rapid simulation of sea clutter based on a backscattering Morchin model. Background Art

[0002] With the development of radar technology, due to the increasing demand for ocean surveillance and detection tasks, sea clutter simulation technology has become an important application direction in radar systems. For marine radars, the sea surface environment is complex and changeable. Factors such as waves and wind speed will cause strong sea clutter interference to radar signals, thus affecting the accurate detection and tracking of targets. Although traditional radars can complete basic detection tasks in simple sea conditions, in order to improve the accuracy and efficiency of target recognition in complex sea conditions, it is necessary to use sea clutter simulation technology to optimize the radar's detection performance. Summary of the invention

[0003] In order to solve the above problems existing in the prior art, the present invention provides a method and device for rapid simulation of sea clutter based on a backscattering Morchin model.

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

[0005] Determine the radar illumination area;

[0006] Dividing the radar illumination area into a plurality of scattering units according to the range ring resolution width and the azimuth resolution;

[0007] Dividing each of the scattering units into two triangular-like regions, and calculating the scattering cross-sectional area of ​​each triangular-like region by a backscattering Morchin model;

[0008] Obtaining an echo signal of a corresponding scattering unit according to the scattering cross-sectional area of ​​each triangular-like region;

[0009] The echo signals of scattering units on the same range ring are superimposed to obtain the simulated sea clutter echo signal.

[0010] Optionally, the scattering cross-sectional area of ​​each triangular-like region is expressed as follows:

[0011]

[0012] Among them, σ m,k represents the scattering cross-sectional area of ​​the triangular region of the scattering unit at the kth azimuth resolution of the mth range ring, represents the local backscattering coefficient of the scattering surface corresponding to the triangular area of ​​the scattering unit at the kth azimuth resolution of the mth range ring, dA m,k represents the area of ​​the triangular region of the scattering unit at the kth azimuth resolution of the mth range ring, ζ m,k It represents the angle between the triangular area of ​​the scattering unit at the kth azimuth resolution of the mth range ring and the standard horizontal plane.

[0013] Optionally, 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] Among them, S(θ i,k ,R i,m ) represents the echo signal received by the radar on the scattering unit at the i-th moment, R i,m represents the distance between the scattering unit on the mth range ring and the origin of the rectangular coordinate system at the i-th time, θ i,k represents the angle between the scattering unit at the kth azimuth resolution at the i-th moment and the x-axis, S1(θ i,k ,R i,m ) represents the echo signal of the first type of triangle area received by the radar at the i-th moment, S2(θ i,k ,R i,m ) represents the echo signal received by the radar on the second type triangular area at the i-th moment.

[0016] Optionally, the echo signal received by the radar on the triangular-like area at the i-th moment is expressed as follows:

[0017]

[0018] Among them, σ m,k represents the scattering cross-sectional area of ​​the triangular region of the scattering unit at the kth azimuth resolution of the mth range ring, j represents the imaginary unit, Represents the echo phase of the triangular scattering unit of the scattering unit at the kth azimuth resolution of the mth range ring.

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

[0020]

[0021] in, represents the simulated sea clutter echo signal corresponding to the mth range ring at the i-th time, P cis the average power required to simulate sea clutter, P sc is the sea clutter power in the radar illumination area, S i,m Represents the initial echo signal corresponding to the mth range ring at the i-th time.

[0022] Optionally, the initial echo signal corresponding to the m-th range ring at the i-th time is expressed as follows:

[0023]

[0024] Among them, N θ is the azimuth resolution number contained in a single range ring, S(θ i,k ,R i,m ) represents the echo signal received by the radar on the scattering unit at the i-th moment, R i,m represents the distance between the scattering unit on the mth range ring and the origin of the rectangular coordinate system at the i-th time, θ i,k It represents the angle between the scattering unit and the x-axis at the k-th azimuth resolution at the i-th moment.

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

[0026]

[0027] Among them, N p Indicates the number of radar pulses, N R Indicates the number of range rings in the radar illumination area.

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

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

[0030] An area division module, used for dividing the radar illumination area into a plurality of scattering units according to the range ring resolution width and the azimuth resolution;

[0031] A scattering cross-sectional area calculation module, used for dividing each of the scattering units into two triangular-like regions, and calculating the scattering cross-sectional area of ​​each triangular-like region by a backscattering Morchin model;

[0032] An echo signal acquisition module, used for obtaining the echo signal of the corresponding scattering unit according to the scattering cross-sectional area of ​​each triangular-like region;

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

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

[0035] Through the above technical scheme, a method for rapid simulation of sea clutter based on the backscattering Morchin model is proposed. The present invention solves the problem of rapid simulation of sea clutter. The generated simulated sea clutter echo signal meets the requirements of having a strong correlation in time, improves the effectiveness and authenticity of sea clutter simulation, has high efficiency, and can rapidly simulate sea clutter.

[0036] Other features and advantages of the present invention will be described in detail in the following detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] The accompanying drawings are used to provide a further understanding of the present invention and constitute a part of the specification. Together with the following specific embodiments, they are used to explain the present invention but do not constitute a limitation of the present invention. In the accompanying drawings:

[0038] Figure 1 The present invention is a flow chart showing a method for rapid simulation of sea clutter based on a backscattering Morchin model according to an exemplary embodiment.

[0039] Figure 2 The figure is a schematic diagram showing a radar illumination area according to an exemplary embodiment.

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

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

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

[0043] Figure 6 It is a schematic diagram of simulation results of simulating sea clutter echo signals according to an exemplary embodiment.

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

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

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

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

[0048] Figure 9b is a schematic diagram showing an average frequency spectrum profile according to an exemplary embodiment.

[0049] Fig.10a The figure is a schematic diagram of a time autocorrelation coefficient variation curve of a simulated sea clutter echo signal according to an exemplary embodiment.

[0050] Fig.10b The figure is a schematic diagram of a curve showing a variation of the range autocorrelation coefficient of a simulated sea clutter echo signal according to an exemplary embodiment.

[0051] Fig.11 It is a block diagram of a fast simulation device for sea clutter based on a backscattering Morchin model according to an exemplary embodiment. DETAILED DESCRIPTION

[0052] Figure 1 is a flow chart of a method for rapid simulation of sea clutter based on a 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] It is understandable that before determining the radar illumination area, a sea level model can be established first. In the rectangular coordinate system, the coordinates of the sea level scatterer at a certain moment can be expressed as:

[0055]

[0056] Where g represents the gravitational acceleration of the sea surface, the y coordinate is fixed, and x0 is the x coordinate at the start time; γ j Obey an independent distribution between 0 and 2π; Obeys a Gaussian distribution with a standard deviation of 0.2rad; α j is a Gaussian random variable reflecting the wave height, ω j is the wave spectrum distribution parameter.

[0057] In one embodiment, Figure 2 is a schematic diagram showing a radar illumination area according to an exemplary embodiment. Figure 2 As shown, assuming that in the downward-looking radar mode, the beam illumination forms an elliptical coverage area, then the coverage area of ​​the radar antenna 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, θ 3dB is the half-power beamwidth of the radar antenna, R is the radar detection range, ψ g is the ground-rubbing angle, and sec() means secant.

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

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

[0062] S103, dividing each scattering unit into two triangular-like regions, and calculating the scattering cross-sectional area of ​​each triangular-like region by using a backscattering Morchin model.

[0063] Understandably, Figure 4 is a schematic diagram of another scattering unit according to an exemplary embodiment. Figure 4 As shown, Figure 3 The scattering unit part can be enlarged to get Figure 4 , each scattering unit is divided into two triangular regions, and the scattering cross-sectional area of ​​each triangular region is calculated based on the backscattering Morchin model. Taking the i-th moment as an example, A m,k is the spatial scattering unit at the kth angle resolution of the mth range ring; dA m,k A m,k Projected area at standard level; R i,m is the distance between the scattering unit on the mth range ring and the origin of the rectangular coordinate system, θ i,k is the angle between the scattering unit on the kth azimuth sector and the x-axis. Figure 4 The form shown converts dA m,k Divided into two triangular forms. The areas of the two triangular scattering units satisfy:

[0064]

[0065] Assume that the direction vector of the radar beam at the i-th moment is The radar coordinates are (x p ,y p ,z p ),Will Figure 4 The cylindrical coordinates of a, b, c, and d are converted into rectangular coordinates respectively (x a ,y a ,z a )、(x b ,y b ,z b )、(x c ,y c ,z c )、(x d ,y d ,z d ). m,k,1 The normal vector of the triangle in this area is composed of three points: a, b, and c. ,and ,

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

[0067]

[0068] Backscatter coefficient σ of sea clutter in Morchin model 0 It can be expressed as:

[0069]

[0070] Where λ is the radar wavelength, ψ g is the ground-brushing angle, S represents the sea state level, ψ c is the floor rubbing angle factor.

[0071] Floor rubbing factor ψ c The definition is as follows:

[0072]

[0073] Among them, h e Represents the roughness of the sea surface (m), and its definition satisfies:

[0074]

[0075] The scattering cross-sectional area of ​​each triangular-like region is expressed as follows:

[0076]

[0077] Among them, σ m,krepresents the scattering cross-sectional area of ​​the triangular region of the scattering unit at the kth azimuth resolution of the mth range ring, represents the local backscattering coefficient of the scattering surface corresponding to the triangular area of ​​the scattering unit at the kth azimuth resolution of the mth range ring, dA m,k represents the area of ​​the triangular region of the scattering unit at the kth azimuth resolution of the mth range ring, ζ m,k It represents the angle between the triangular area of ​​the scattering unit at the kth azimuth resolution of the mth range ring and the standard horizontal plane.

[0078] Take the first kind of triangle area dA m,k,1 For example, 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 area and the standard horizontal plane can be obtained m,k,1 It can be expressed as follows:

[0079]

[0080] The obtained ζ m,k,1 Substituting the scattering cross-sectional area formula of the triangular region into the scattering cross-sectional area formula of the first type of triangular region, we can get the scattering cross-sectional area σ m,k,1 .

[0081] S104. Obtain the echo signal of the corresponding scattering unit according to the scattering cross-sectional area of ​​each quasi-triangle region.

[0082] It can be understood that, first, according to the triangle centroid calculation formula, A m,k,1 Quasi-triangle area and A m,k,2 The centroids of the triangular-like regions are:

[0083]

[0084] The distances of the two triangular scattering units from the radar are:

[0085]

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

[0087]

[0088] According to the echo phase of each triangular scattering unit, the echo signal received by the radar on the triangular area at the i-th moment can be expressed as follows:

[0089]

[0090] Among them, σ m,k represents the scattering cross-sectional area of ​​the triangular region of the scattering unit at the kth azimuth resolution of the mth range ring, j represents the imaginary unit, Represents the echo phase of the triangular scattering unit of the scattering unit at the kth azimuth resolution of the mth range ring.

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

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

[0093] Among them, S(θ i,k ,R i,m ) represents the echo signal received by the radar on the scattering unit at the i-th moment, R i,m represents the distance between the scattering unit on the mth range ring and the origin of the rectangular coordinate system at the i-th time, θ i,k represents the angle between the scattering unit at the kth azimuth resolution at the i-th moment and the x-axis, S1(θ i,k ,R i,m ) represents the echo signal of the first type of triangle area received by the radar at the i-th moment, S2(θ i,k ,R i,m ) represents the echo signal received by the radar on the second type triangular area at the i-th moment.

[0094] S105 , superimposing the echo signals of the scattering units on the same range ring to obtain a simulated sea clutter echo signal.

[0095] It can be understood that the vector superposition of the scattering units in the same range ring of the radar can be expressed as follows:

[0096]

[0097] Among them, N θ is the azimuth resolution number contained in a single range ring, S(θ i,k ,R i,m ) represents the echo signal received by the radar on the scattering unit at the i-th moment, R i,m represents the distance between the scattering unit on the mth range ring and the origin of the rectangular coordinate system at the i-th time, θ i,k It represents the angle between the scattering unit and the x-axis at the k-th azimuth resolution at the i-th moment.

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

[0099]

[0100] Among them, N p Indicates the number of radar pulses, N R Indicates the number of range rings in the radar illumination area.

[0101] According to the sea clutter power in the radar illumination area and the initial echo signal corresponding to the mth range ring at the i-th time, the simulated sea clutter echo signal can be expressed as follows:

[0102]

[0103] in, represents the simulated sea clutter echo signal corresponding to the mth range ring at the i-th time, P c is the average power required to simulate sea clutter, P sc is the sea clutter power in the radar illumination area, S i,m Represents the initial echo signal corresponding to the mth range ring at the i-th time.

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

[0105] Table 1

[0106] Parameter name Parameter settings Radar altitude / m 100 Wiping angle / ° 5 Pulse width / μs 2 Pulse number 100 Bandwidth / MHz 1 Sampling frequency / MHz 10 Carrier frequency / GHz 2 Beam width / ° 3.3 Azimuth resolution / ° 0.05 Polarization HH Pulse repetition period / μs 480 sea ​​feeling 3 Wind speed / m / s 6 C / N ratio / dB 40

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

[0108] The statistical characteristics of the sea clutter simulation signal are analyzed, and the parameters of Rayleigh distribution, Log-Normal distribution, Weibull distribution and K distribution are 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, we can get a comparative diagram. Figure 8ais a schematic diagram showing a comparison between a simulated sea clutter echo signal and a classical statistical model according to an exemplary embodiment. Figure 8b is a schematic diagram showing 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 of the comparison between the simulated sea clutter echo signal and the Rayleigh distribution and Log-Normal distribution. Figure 8b Schematic diagram of the comparison between the simulated sea clutter echo signal and the Weibull distribution and K distribution.

[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 models MSD test value Rayleigh distribution 0.0304 Log-Normal Distribution 0.0400 Weibull distribution 0.0270 K-distribution 0.0231

[0115] MSD goodness of fit test method, the smaller the MSD test value, the better the goodness of fit. As shown in Table 3, it can be found that the simulated sea clutter signal has the best fitting effect with the K distribution statistical model.

[0116] Coherent accumulation is used to perform spectrum analysis on the simulated sea clutter signal. The number of sea clutter signal pulses is 100 and the number of FFT points is 256. Figure 9a is a schematic diagram of a sea clutter spectrum according to an exemplary embodiment. Figure 9b is a schematic diagram of an average spectrum according to an exemplary embodiment, a sea clutter spectrum and an average spectrum profile of a sea clutter region (50 to 100 distance units) are shown as follows Figure 9a and Figure 9b As shown in 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 be broadened and the clutter Doppler center is not zero, which is consistent with the actual situation.

[0117] The correlation between the time dimension and distance dimension of the sea clutter simulation signal is analyzed. Fig.10a is a schematic diagram of a time autocorrelation coefficient variation curve of a simulated sea clutter echo signal according to an exemplary embodiment. Fig.10b FIG. 1 is a schematic diagram of a distance autocorrelation coefficient variation curve of a simulated sea clutter echo signal according to an exemplary embodiment. It is generally believed that if the correlation coefficient is greater than e -1 , it is believed that the signals are correlated, otherwise, it is believed that there is no correlation between the signals. Fig.10aThe different colored curves in the figure represent the time autocorrelation coefficient variation curves at different distance units, and the thick black line represents the variation curve of the mean value of the time autocorrelation coefficient over time; it can be seen that the decorrelation time of the simulated sea clutter echo signal is about 12ms. Fig.10b One color represents the curve of the distance autocorrelation coefficient at a pulse with the distance unit, and the thick black line represents the curve of the mean distance autocorrelation coefficient with the distance unit. It can be seen that the mean distance autocorrelation coefficient of sea clutter drops to e after 1 distance unit. -1 It is below 0.3 and quickly stabilizes below 0.3. It can be seen that the simulated sea clutter echo signal is random in the distance dimension.

[0118] It can be seen from the above simulation results that the simulated sea clutter echo signal generated by the present invention conforms to the K distribution statistical model; and the sea clutter has a relatively strong correlation in time, which proves the effectiveness and superiority of the sea clutter fast simulation method.

[0119] Through the above technical solutions, the traditional sea clutter simulation method has many limitations in efficiency and accuracy. For example, the spectral model requires complex ocean environment parameters, which is time-consuming to obtain and process, and it is difficult to capture subtle features under complex sea conditions; the statistical model relies on a large amount of historical data, has poor generalization ability, and is difficult to cope with unprecedented extreme sea conditions; the numerical simulation has a large amount of calculation and is difficult to apply in real time; and the neural network simulation is prone to overfitting when the training data is insufficient, and the performance is unstable in the real environment. In contrast, the present invention can achieve fast simulation, and the generated sea clutter statistical characteristics conform to the K distribution model, with strong time correlation, which greatly improves the effectiveness and authenticity of the simulation, can significantly improve the speed and accuracy of sea clutter simulation, and provides more reliable technical support for radar detection in marine environments.

[0120] Fig.11 is a block diagram of a fast sea clutter simulation device based on a backscattering Morchin model according to an exemplary embodiment, see Fig.11 , the apparatus 1100 may include:

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

[0122] The area division module 1102 is used to divide the radar illumination area into a plurality of scattering units according to the range ring resolution width and the azimuth resolution;

[0123] A scattering cross-sectional area calculation module 1103 is used to divide each scattering unit into two triangular-like regions, and calculate the scattering cross-sectional area of ​​each triangular-like region by using a backscattering Morchin model;

[0124] An echo signal acquisition module 1104 is used to obtain the echo signal of the corresponding scattering unit according to the scattering cross-sectional area of ​​each triangular-like region;

[0125] The simulation module 1105 is used to superimpose the echo signals of the scattering units on the same range ring to obtain a simulated sea clutter echo signal.

[0126] Regarding the device in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0127] The preferred embodiments of the present invention are described in detail above in conjunction with the accompanying drawings. However, the present invention is not limited to the specific details in the above embodiments. Within the technical concept of the present invention, a variety of simple modifications can be made to the technical solution of the present invention, and these simple modifications all belong to 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 further describe various possible combinations.

[0129] In addition, various embodiments of the present invention may be arbitrarily combined, and as long as they do not violate the concept of the present invention, they should also be regarded as the contents disclosed by the present invention.

Claims

1. A fast simulation method for sea clutter based on backscattering Morchin model, characterized in that: The method comprises: Determine the radar illumination area; Dividing the radar illumination area into a plurality of scattering units according to the range ring resolution width and the azimuth resolution; Dividing each of the scattering units into two triangular-like regions, and calculating the scattering cross-sectional area of ​​each triangular-like region by a backscattering Morchin model; Obtaining an echo signal of a corresponding scattering unit according to the scattering cross-sectional area of ​​each triangular-like region; The echo signals of scattering units on the same range ring are superimposed to obtain the simulated sea clutter echo signal.

2. The method for rapid simulation of 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 expressed as follows: Among them, σ m,k represents the scattering cross-sectional area of ​​the triangular region of the scattering unit at the kth azimuth resolution of the mth range ring, represents the local backscattering coefficient of the scattering surface corresponding to the triangular area of ​​the scattering unit at the kth azimuth resolution of the mth range ring, dA m,k represents the area of ​​the triangular region of the scattering unit at the kth azimuth resolution of the mth range ring, ζ m,k It represents the angle between the triangular area of ​​the scattering unit at the kth azimuth resolution of the mth range ring and the standard horizontal plane.

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

4. The method for rapid simulation of sea clutter based on the backscattering Morchin model according to claim 3, characterized in that: The echo signal received by the radar on the triangular area at the i-th moment is expressed as follows: Among them, σ m,k represents the scattering cross-sectional area of ​​the triangular region of the scattering unit at the kth azimuth resolution of the mth range ring, j represents the imaginary unit, Represents the echo phase of the triangular scattering unit of the scattering unit at the kth azimuth resolution of the mth range ring.

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

6. The method for rapid simulation of sea clutter based on backscattering Morchin model according to claim 5, characterized in that: The initial echo signal corresponding to the mth range ring at the i-th time is expressed as follows: Among them, N θ is the azimuth resolution number contained in a single range ring, S(θ i,k ,R i,m ) represents the echo signal received by the radar on the scattering unit at the i-th moment, R i,m represents the distance between the scattering unit on the mth range ring and the origin of the rectangular coordinate system at the i-th time, θ i,k It represents the angle between the scattering unit and the x-axis at the k-th azimuth resolution at the i-th moment.

7. The method for rapid simulation of sea clutter based on backscattering Morchin model according to claim 6, characterized in that: The sea clutter power in the radar illumination area is expressed as follows: Among them, N p Indicates the number of radar pulses, N R Indicates the number of range rings in the radar illumination area.

8. A fast simulation device for sea clutter based on backscattering Morchin model, characterized in that: The device comprises: An area determination module is used to determine the radar illumination area; An area division module, used for dividing the radar illumination area into a plurality of scattering units according to the range ring resolution width and the azimuth resolution; A scattering cross-sectional area calculation module, used for dividing each of the scattering units into two triangular-like regions, and calculating the scattering cross-sectional area of ​​each triangular-like region by a backscattering Morchin model; An echo signal acquisition module, used for obtaining the echo signal of the corresponding scattering unit according to the scattering cross-sectional area of ​​each triangular-like region; The simulation module is used to superimpose the echo signals of the scattering units on the same range ring to obtain the simulated sea clutter echo signal.

Images