Microwave detection performance simulation analysis method, equipment and medium

Abstract

The present invention relates to a microwave detection performance simulation analysis method, equipment and storage medium. The microwave detection performance simulation analysis method comprises: step S1, calculating the detection parameters of the microwave detection load and the coverage of the antenna beam according to the microwave detection load parameters; step S2, simulating and calculating the radar echo of the target according to the relative position and attitude relationship between the microwave detection load and the target; step S3, calculating the scattered echo of each clutter block to generate environmental clutter; step S4, performing spatiotemporal adaptive processing on the target, the environment and the interference echo, and calculating the microwave detection performance parameters according to the target fluctuation type and the spatiotemporal processing gain. The present invention can greatly improve the accuracy of the simulation and the consistency with the actual working performance of the equipment, improve the accuracy of the system simulation, avoid the large calculation requirements in the signal level simulation, can well ensure the real-time performance of the simulation analysis, and concisely and efficiently analyze the performance of the microwave detection load accurately.

Description

Technical Field

[0001] The present invention relates to the technical field of system modeling and simulation, and in particular to a microwave detection effectiveness simulation analysis method, equipment and medium. Background Art

[0002] In recent years, with the widespread application of model-based system engineering (MBSE) in complex engineering fields such as aviation, aerospace and space systems, the accuracy of subsystem models and their consistency with actual systems have become key factors for the successful design of complex system engineering. However, in the system simulation of microwave detection payloads in aircraft and other systems, the model accuracy is not enough, and it is difficult to accurately simulate and analyze its performance through system simulation technology. On the one hand, in order to adapt to different simulation requirements, each subsystem model needs to use parameters and simulation models to quickly give simulation results, which makes it difficult for the current simulation model to fully consider the impact of the environment and interference on the microwave detection payload. If the signal-level simulation process in the design of microwave detection payload is directly used, the simulation calculation time is too long, which seriously reduces the efficiency and real-time performance of the entire system simulation.

[0003] Therefore, in response to the current simulation needs of model-based system engineering, there is an urgent need for a performance simulation analysis method that takes both simulation speed and simulation accuracy into consideration. Summary of the invention

[0004] In order to solve the technical problems existing in the above-mentioned prior art, the purpose of the present invention is to provide a microwave detection performance simulation analysis method, equipment and storage medium, which can simulate the advantages of signal-level simulation, comprehensively consider the impact of the environment and interference on the performance in the entire workflow of the equipment, reduce the calculation amount of model simulation, and meet the real-time requirements in system simulation.

[0005] To achieve the above-mentioned object of the invention, the present invention provides a microwave detection efficiency simulation analysis method, which comprises the following steps:

[0006] Step S1, calculating the detection parameters of the microwave detection payload and the antenna beam coverage;

[0007] Step S2, simulating and calculating the radar echo of the target according to the relative position and attitude relationship between the microwave detection payload and the target;

[0008] Step S3, calculating the scattered echo of each clutter block to generate environmental clutter;

[0009] Step S4: Perform spatiotemporal adaptive processing on the target, environment and interference echoes, and calculate microwave detection effectiveness parameters according to the target fluctuation type and spatiotemporal processing gain.

[0010] According to a technical solution of the present invention, the step S1 specifically includes:

[0011] Step S11, according to the load platform parameters of the microwave detection payload, determine the current velocity vector and height parameters of the load platform, and determine the geometric coordinate relationship between the detection area and the load platform;

[0012] Step S12: Modeling the three-dimensional radiation pattern of the microwave detection payload according to the antenna parameters of the microwave detection payload to obtain the antenna radiation pattern of the microwave detection payload that satisfies the main lobe beam width;

[0013] Step S13: Calculate the beam coverage of the main beam according to the geometric coordinate relationship between the detection area and the payload platform and the antenna pattern.

[0014] According to a technical solution of the present invention, the step S2 specifically includes:

[0015] Step S21, calculating the radar cross-section data of the target according to the material parameters and shape parameters of the target, and determining the scattering type;

[0016] Step S22: Calculate the echo amplitude of the backscattered echo of the target according to the radar cross section data of the target and the radar echo equation.

[0017] According to a technical solution of the present invention, the step S22 specifically includes:

[0018] According to the attitude of the microwave detection payload and its relative position to the target, the antenna gain is obtained through the antenna pattern model. The RCS in the corresponding direction is obtained through the scattering cross-section data of the target, that is,

[0019] The echo amplitude of the target's backscattered echo is calculated according to the radar equation. The formula is:

[0020]

[0021] Among them, P av is the average transmitted power, is the antenna gain in the detection direction in the detection payload coordinate system, λ is the wavelength, is the RCS of the target in its local coordinate system in the direction of the microwave detection load, k is the Boltzmann constant, T0 is the equivalent temperature, B n is the Doppler bandwidth, F n is the receiver noise figure, K s is the system loss.

[0022] According to a technical solution of the present invention, step S3 specifically includes:

[0023] Step S31, establishing a calculation model for the surface clutter backscattering coefficient of the microwave detection payload irradiation area according to the ground type;

[0024] Step S32, calculating the area and number of clutter blocks according to the distance relationship between the target and the radar and the resolution of the microwave detection payload;

[0025] Step S33, calculating the Doppler shift of the clutter block, and calculating the reflected echo of each clutter block according to the area, azimuth and back reflection coefficient of the clutter block;

[0026] Step S34: Calculate the echo parameters of the environmental clutter according to the number of clutter blocks and the reflected echo of each clutter block.

[0027] According to a technical solution of the present invention, the step S32 specifically includes:

[0028] Step S321, completing distance grid division;

[0029] The distance resolution of the microwave detection payload on the ground is used to divide the distance units. The distance resolution is calculated using the following formula:

[0030]

[0031] Where ψ is the angle between the incident direction of the radar pulse and the ground; R Within the range, two different targets can be distinguished by the microwave detection payload. In the range direction, the main beam coverage area is divided into multiple equally spaced range rings. The distance from any range ring to the microwave detection payload is expressed as R. i express;

[0032] Step S322, completing azimuth grid division;

[0033] According to the Mitchell model, the azimuth direction is divided according to the angular resolution of the radar, f r represents the pulse repetition frequency, and the Doppler bandwidth of the microwave detection payload is represents the elevation angle of the microwave detection payload radar antenna, and the Doppler resolution is f r / M, then in the azimuth direction, the number of units is the number of clutter ring units, and we have:

[0034] K=Δf d M / f r

[0035] The azimuth division size is:

[0036]

[0037] The area of ​​the clutter ring is calculated as:

[0038]

[0039] Step S323, calculating the distance blur quantity;

[0040] The radar transmits a pulse train at a fixed frequency, and any pulse train forms a series of circular areas on the ground.

[0041] According to the period of the pulse signal, the maximum unambiguous distance R is calculated u :

[0042]

[0043] Among them, T p Represents the pulse period, which is the same as the pulse repetition frequency f r The two are reciprocal;

[0044] Assuming that the grazing angles in the radar main beam are approximately equal, the amount of range ambiguity in the main lobe is expressed as:

[0045]

[0046] According to a technical solution of the present invention, the step S4 specifically includes:

[0047] Step S41, setting the number, power and movement of interference sources, and calculating the power reaching the position of the microwave detection payload antenna according to the power attenuation equation;

[0048] Step S42, calculating the space-time steering vector and generating a space-time adaptive processing matrix;

[0049] Step S43: Calculate the target detection probability according to the target fluctuation type and the signal-to-noise ratio after spatiotemporal processing.

[0050] According to a technical solution of the present invention, the step S42 specifically includes:

[0051] Step S421: Calculate the spatial steering vector

[0052] In the microwave detection payload coordinate system, the microwave detection payload flies along the X-axis, the array antenna is distributed along the X-axis, and the position of the first antenna unit is taken as the coordinate origin. The coordinate position of the i-th channel antenna is expressed as:

[0053] coor i rad =((i-1)d, 0, 0)

[0054] For the position at a pitch angle of θ EL , the azimuth angle is θ AZThe target at the position is represented by coor=R(sinθ EL cosθ AZ , sinθ EL sinθ AZ , cosθ EL ), its vector distance to the ith channel antenna is expressed as

[0055] Dis i rad =(Rsinθ EL cosθ AZ -(i-1)d, Rsinθ EL sinθ AZ , Rcosθ EL );

[0056] The vector difference between the i-th channel antenna and the first vector antenna is expressed as:

[0057] Dis i rad -Dis0 rad =((i-1)d, 0, 0)

[0058] The projection of the i-th channel antenna and the first vector antenna in the target direction is expressed as:

[0059] (i-1)dsinθ EL cosθ AZ

[0060] The phase difference between the i-th channel antenna and the first vector antenna is expressed as:

[0061]

[0062] The normalized angle is defined as:

[0063]

[0064] Then the spatial steering vector is expressed as:

[0065]

[0066] Step S422, calculating the time-oriented vector;

[0067] Assuming the target is stationary, the Doppler frequency is mainly determined by the relative motion speed between the target and the microwave detection payload in the radial direction. The projection of the velocity vector on the distance vector between the target and the radar is the radial velocity.

[0068] The velocity vector is along the X-axis in the satellite coordinate system, so the velocity vector is:

[0069] V rad=(V, 0, 0)

[0070] The radial velocity is:

[0071] v r rad =V*sinθ EL cosθ AZ

[0072] The Doppler frequency is:

[0073]

[0074] Then the Doppler time steering vector is:

[0075]

[0076] Step S423: Generate a space-time covariance matrix R according to the target, environment and interference power at the microwave detection payload, and obtain a space-time adaptive processing matrix by using space-time optimization processing:

[0077]

[0078] According to one aspect of the present invention, an electronic device is provided, comprising: one or more processors, one or more memories, and one or more computer programs; wherein the processor is connected to the memory, and the one or more computer programs are stored in the memory, and when the electronic device is running, the processor executes the one or more computer programs stored in the memory, so that the electronic device performs a microwave detection effectiveness simulation analysis method as described in any one of the above technical solutions.

[0079] According to one aspect of the present invention, a computer-readable storage medium is provided for storing computer instructions. When the computer instructions are executed by a processor, a microwave detection effectiveness simulation analysis method as described in any one of the above technical solutions is implemented.

[0080] Compared with the prior art, the present invention has the following beneficial effects:

[0081] The present invention proposes a microwave detection effectiveness simulation and analysis method, equipment and storage medium. The microwave detection effectiveness simulation and analysis method models and simulates the main signal processing process of the microwave detection payload through core methods such as time-space adaptive processing, and finally gives the detection probability of the microwave detection payload under the influence of environment and interference according to the relationship between the target fluctuation type and the detection probability. It can simulate the data processing process of the entire microwave detection payload, and analyze and model the main work process of the microwave detection payload, which solves the accuracy problem of the model in the current system simulation. By modeling and analyzing the environment such as microwave detection payload signal transmission, environmental influence, interference influence and main signal processing algorithm, the accuracy of the simulation and the consistency with the actual working efficiency of the equipment can be greatly improved, thereby improving the accuracy of the system simulation.

[0082] Furthermore, the present invention performs statistical modeling through various environments for signal sequence processing without processing the actual signal sequence, thus avoiding the large computing requirements in signal-level simulation, and can well ensure the real-time performance of the simulation analysis, and can accurately analyze the effectiveness of the microwave detection payload in a concise and efficient manner. BRIEF DESCRIPTION OF THE DRAWINGS

[0083] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments are briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention, and for ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0084] Figure 1 A schematic diagram of a flow chart of a microwave detection effectiveness simulation analysis method provided in accordance with an embodiment of the present invention;

[0085] Figure 2 A specific flow chart of step S1 in the flow chart of the microwave detection effectiveness simulation analysis method according to one embodiment of the present invention is schematically shown;

[0086] Figure 3 A specific flow chart schematically showing step S2 in the flow chart of the microwave detection effectiveness simulation analysis method according to one embodiment of the present invention;

[0087] Figure 4 A specific flow chart schematically showing step S3 in the flow chart of the microwave detection effectiveness simulation analysis method according to one embodiment of the present invention;

[0088] Figure 5 A specific flow chart of step S4 in the flow chart of the microwave detection effectiveness simulation analysis method according to one embodiment of the present invention is schematically shown;

[0089] Figure 6Schematically represents the coordinate system in the embodiment of the present invention;

[0090] Figure 7 The diagram schematically shows the antenna beam coverage in the embodiment of the present invention. DETAILED DESCRIPTION

[0091] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0092] like Figures 1 to 5 As shown, a microwave detection performance simulation analysis method of the present invention comprises the following steps:

[0093] Step S1, calculating the detection parameters of the microwave detection payload and the antenna beam coverage according to the microwave detection payload parameters;

[0094] Step S2, simulating and calculating the radar echo of the target according to the relative position and attitude relationship between the microwave detection payload and the target;

[0095] Step S3, calculating the scattered echo of each clutter block to generate environmental clutter;

[0096] Step S4: Perform spatiotemporal adaptive processing on the target, environment and interference echoes, and calculate microwave detection effectiveness parameters according to the target fluctuation type and spatiotemporal processing gain.

[0097] like Figure 2 As shown, in one embodiment of the present invention, preferably, the step S1 specifically includes:

[0098] Step S11, according to the load platform parameters of the microwave detection payload, determine the current velocity vector and height parameters of the load platform, and determine the geometric coordinate relationship between the detection area and the load platform; specifically, Figure 6 As shown, the coordinate system in this embodiment involves 3 coordinate systems in total.

[0099] The first is a ground-fixed coordinate system. A global coordinate system is established at any fixed position on the ground, with the XY plane parallel to the ground and the Z axis pointing to the sky.

[0100] The second is the microwave detection payload coordinate system. A local coordinate system is established on the microwave detection payload, with the positive direction of the Z axis perpendicular to the ground, the direction of the X axis being the same as the velocity direction, and the direction of the Y axis being determined according to the right-hand relationship. In this embodiment, it is represented by the superscript rad.

[0101] The third is the target coordinate system. Its X-axis is the projection of the line connecting the target and the microwave detection payload on the ground, the Z-axis is parallel to the global coordinate system, and the Y-axis is determined according to the right-hand relationship. In this embodiment, it is represented by the superscript tar.

[0102] Assume that the detection load moves in the horizontal direction, the height from the ground is H, the speed is V, and the yaw angle between the carrier speed direction and the axial direction is φ v , in the global coordinate system, the carrier velocity vector At the same time, the carrier frequency of the microwave detection payload is determined to be f c , wavelength is λ, pulse repetition frequency is f r , M pulses are emitted in one CPI.

[0103] Step S12: Modeling the three-dimensional radiation pattern of the microwave detection payload according to the antenna parameters of the microwave detection payload to obtain the antenna radiation pattern of the microwave detection payload that satisfies the main lobe beam width;

[0104] There are generally two cases to obtain the antenna pattern of the microwave detection payload that meets the main lobe beam width:

[0105] One is to import the external antenna pattern, which is suitable for the gain data of each direction of the known measured antenna;

[0106] The other is the key parameter function simulation of the directional pattern, which is suitable for antenna gain simulation without measured data or design data, and the antenna direction of a single channel Figure 1 Generally, the form of Sinque function is adopted. Its mathematical expression is

[0107]

[0108] Among them, γ rad and are the azimuth and elevation angles of the microwave detection payload in its local coordinate system; and are the beam widths (first zero point widths) of the microwave detection payload antenna beam in azimuth and elevation directions respectively; and are the elevation angle and azimuth angle of the radar antenna respectively. θ and In order to fine-tune the antenna pattern, it is generally determined by calculating the value of the half-power beamwidth.

[0109] If the sum and difference beams are considered, a four-channel antenna pattern can be used to divide the antenna array into four channels: upper left, lower left, upper right, and lower right. Since the four channels only point in slightly different directions, one channel can be calculated and the other channels can be converted from it. In fact, the difference between the channels is the phase difference in the elevation and azimuth directions.

[0110] Step S13: Calculate the beam coverage of the main beam according to the geometric coordinate relationship between the detection area and the payload platform and the antenna pattern.

[0111] like Figure 7 As shown in the figure, taking the calculation of the main lobe beam coverage of the microwave detection payload pitch plane as an example, in the detection payload coordinate system, the angle between the beam pointing and the Z axis is The null beam width is Therefore, the pitch angles of the main lobe beam in the pitch direction corresponding to the nearest and farthest points of the ground coverage are:

[0112]

[0113] According to the height H, the distances from the main lobe coverage in the pitch direction to the projection point of the microwave detection on the ground can be calculated as follows:

[0114]

[0115] Since the echo of the microwave detection payload on the ground is composed of a series of annular areas, the azimuth direction is equivalent to direct coverage using the entire antenna pattern.

[0116] like Figure 3 As shown, in one embodiment of the present invention, preferably, the step S2 specifically includes:

[0117] Step S21, calculating the radar cross-section data of the target according to the material parameters and shape parameters of the target, and determining the scattering type;

[0118] Step S22: Calculate the echo amplitude of the backscattered echo of the target according to the radar cross section data of the target and the radar echo equation.

[0119] In one embodiment of the present invention, preferably, the step S22 specifically includes:

[0120] According to the attitude of the microwave detection payload and its relative position to the target, the antenna gain is obtained through the antenna pattern model. The RCS in the corresponding direction is obtained through the scattering cross-section data of the target, that is,

[0121] The echo amplitude of the target's backscattered echo is calculated according to the radar equation. The formula is:

[0122]

[0123] Among them, P av is the average transmit power, is the antenna gain in the detection direction in the detection payload coordinate system, λ is the wavelength, is the RCS of the target in its local coordinate system in the direction of the microwave detection load, k is the Boltzmann constant, T0 is the equivalent temperature, B n is the Doppler bandwidth, F n is the receiver noise figure, L s is the system loss.

[0124] like Figure 4 As shown, in one embodiment of the present invention, preferably, the step S3 specifically includes:

[0125] Step S31, establishing a calculation model for the surface clutter backscattering coefficient of the microwave detection payload irradiation area according to the ground type;

[0126] The irradiation area is mainly modeled by surface clutter, which includes ground clutter and sea clutter. The surface clutter reflectivity v is generally used. 0 A measure of the size of surface clutter. Surface clutter reflectivity is also called radar backscatter coefficient, which is the backscatter cross section per unit area. 0 With the clutter grazing angle θ cr The changing relationship is as follows.

[0127] σ 0 =γsinθ cr

[0128] θ cr It represents the clutter area-radar ray elevation angle. γ is the backscattering radar cross section per unit projected area. For a specific surface type, γ can be statistically approximated to be a constant.

[0129] For ground clutter, the statistical parameters of γ in the model for typical landforms are shown in Table 1. The surface reflectivity of ground clutter can be directly calculated based on the data in the table.

[0130]

[0131] Table 1

[0132] For sea clutter, the calculation of γ is slightly more complicated. γ is related to the radar wavelength λ and the sea state S (or Beaufort wind scale K B ) are as follows.

[0133] [γ] dB =6S-10logλ-58=6K B -10logλ-64.

[0134] Step S32, calculating the area and number of clutter blocks according to the distance relationship between the target and the radar and the resolution of the microwave detection payload, specifically includes:

[0135] Step S321, completing distance grid division;

[0136] The distance resolution of the microwave detection payload on the ground is used to divide the distance units. The distance resolution is calculated using the following formula:

[0137]

[0138] Where ψ is the angle between the incident direction of the radar pulse and the ground; R Within the range, two different targets can be distinguished by the microwave detection payload. In the range direction, the main beam coverage area is divided into multiple equally spaced range rings. The distance from any range ring to the microwave detection payload is expressed as R. i express;

[0139] Step S322, completing azimuth grid division;

[0140] According to the Mitchell model, the Doppler bandwidth of the microwave detection payload is divided according to the angular resolution of the radar in the azimuth direction. The Doppler resolution is f r / M, then in the azimuth direction, the number of units is the number of clutter ring units, and we have:

[0141] K=Δf d M / f r

[0142] The azimuth division size is:

[0143]

[0144] The area of ​​the clutter ring is calculated as:

[0145]

[0146] Step S323, calculating the distance blur quantity;

[0147] The radar transmits a pulse train at a certain frequency, and each pulse train forms a series of circular areas on the ground. The interval between each circular area is determined by the period of the pulse signal, also called the maximum unambiguous distance, and is represented by the symbol R u express.

[0148] According to the period of the pulse signal, the maximum unambiguous distance R is calculated u :

[0149]

[0150] Among them, T p Represents the pulse period, which is the same as the pulse repetition frequency f r The two are reciprocal;

[0151] Assuming that the grazing angles in the radar main beam are approximately equal, the amount of range ambiguity in the main lobe is expressed as:

[0152]

[0153] Step S33, calculating the Doppler shift of the clutter block, and calculating the reflected echo of each clutter block according to the area and orientation of each clutter block; the details are as follows:

[0154] For the kth clutter unit on the current i-th range ring, its Doppler frequency shift can be written as:

[0155]

[0156] Its backscattered signal can be written as:

[0157]

[0158] After being substituted into the radar equation, the signal s at the microwave detection payload is calculated i (k).

[0159] Step S34: According to the number of range rings and the number of clutter units, the number of clutter blocks in the entire calculation process can be obtained. According to the number of clutter blocks and the reflected echo, the echo parameter calculation process of the environmental clutter is calculated as follows:

[0160]

[0161] like Figure 5 As shown, in one embodiment of the present invention, preferably, the step S4 specifically includes:

[0162] Step S41, setting the number, power and movement of interference sources, and calculating the power reaching the position of the microwave detection payload antenna according to the power attenuation equation;

[0163] Step S42, calculating the space-time steering vector and generating a space-time adaptive processing matrix;

[0164] Step S43: Calculate the target detection probability according to the target fluctuation type and the signal-to-noise ratio after spatiotemporal processing.

[0165] In one embodiment of the present invention, preferably, the step S42 specifically includes:

[0166] Step S421: Calculate the spatial steering vector

[0167] In the microwave detection payload coordinate system, the microwave detection payload flies along the X-axis, the array antenna is distributed along the X-axis, and the position of the first antenna unit is taken as the coordinate origin. The coordinate position of the i-th channel antenna is expressed as:

[0168] coor i rad =((i-1)d, 0, 0)

[0169] For the position at a pitch angle of θ EL , the azimuth angle is θ AZ The target at the position is represented by coor=R(sinθ EL cosθ AZ , sinθ EL sinθ AZ , cosθ EL ), its vector distance to the ith channel antenna is expressed as

[0170] Dis i rad =(Rsinθ EL cosθ AZ -(i-1)d, Rsinθ EL sinθ AZ , Rcosθ EL );

[0171] The vector difference between the i-th channel antenna and the first vector antenna is expressed as:

[0172] Dis i rad -Dis0 rad =((i-1)d, 0, 0)

[0173] The projection of the i-th channel antenna and the first vector antenna in the target direction is expressed as:

[0174] (i-1)dsinθ EL cosθ AZ

[0175] The phase difference between the i-th channel antenna and the first vector antenna is expressed as:

[0176]

[0177] The normalized angle is defined as:

[0178]

[0179] Then the spatial steering vector is expressed as:

[0180]

[0181] Step S422, calculating the time-oriented vector;

[0182] Assuming the target is stationary, the Doppler frequency is mainly determined by the relative motion speed between the target and the microwave detection payload in the radial direction. The projection of the velocity vector on the distance vector between the target and the radar is the radial velocity.

[0183] The velocity vector is along the X-axis in the satellite coordinate system, so the velocity vector is:

[0184] V rad =(V, 0, 0)

[0185] The radial velocity is:

[0186] v r rad =V*sinθ EL cosθ AZ

[0187] The Doppler frequency is:

[0188]

[0189] Then the Doppler time steering vector is:

[0190]

[0191] Step S423: Generate a space-time covariance matrix R according to the target, environment and interference power at the microwave detection payload, and obtain a space-time adaptive processing matrix by using space-time optimization processing:

[0192]

[0193] According to one aspect of the present invention, there is provided an electronic device, comprising: one or more processors, one or more memories, and one or more computer programs; wherein the processor is connected to the memory, and the one or more computer programs are stored in the memory, and when the electronic device is running, the processor executes the one or more computer programs stored in the memory, so that the electronic device performs a microwave detection effectiveness simulation analysis method as any one of the above technical solutions.

[0194] According to one aspect of the present invention, a computer-readable storage medium is provided for storing computer instructions. When the computer instructions are executed by a processor, a microwave detection effectiveness simulation analysis method as described in any one of the above technical solutions is implemented.

[0195] Computer-readable storage media may include any medium capable of storing or transmitting information. Examples of computer-readable storage media include electronic circuits, semiconductor memory devices, ROM, flash memory, erasable ROM (EROM), floppy disks, CD-ROMs, optical disks, hard disks, optical fiber media, radio frequency (RF) links, etc. Code segments may be downloaded via computer networks such as the Internet, intranets, etc.

[0196] The present invention provides a microwave detection performance simulation analysis method, device and storage medium, the microwave detection performance simulation analysis method includes: step S1, calculating the detection parameters of the microwave detection load and the antenna beam coverage according to the microwave detection load parameters; step S2, simulating and calculating the radar echo of the target according to the relative position and attitude relationship between the microwave detection load and the target; step S3, calculating the scattered echo of each clutter block to generate environmental clutter; step S4, performing spatiotemporal adaptive processing on the target, the environment and the interference echo, and calculating the microwave detection performance parameters according to the target fluctuation type and spatiotemporal processing gain. Through the above simulation analysis steps, the data processing flow of the entire microwave detection load can be simulated, and the analysis modeling of the main work flow of the microwave detection load can be solved, which solves the accuracy problem of the model in the current system simulation. Through the modeling analysis of the environment such as microwave detection load signal transmission, environmental influence, interference influence and main signal processing algorithm, the accuracy of the simulation and the consistency with the actual working efficiency of the equipment can be greatly improved, and the accuracy of the system simulation can be improved.

[0197] Furthermore, the present invention performs statistical modeling through various environments for signal sequence processing without processing the actual signal sequence, thus avoiding the large computing requirements in signal-level simulation, and can well ensure the real-time performance of the simulation analysis, and can accurately analyze the effectiveness of the microwave detection payload in a concise and efficient manner.

[0198] In addition, it should be noted that the present invention can be provided as a method, an apparatus or a computer program product. Therefore, the embodiments of the present invention can take the form of a complete hardware embodiment, a complete software embodiment or an embodiment combining software and hardware. Moreover, the embodiments of the present invention can take the form of a computer program product implemented on one or more computer-usable storage media containing computer-usable program code.

[0199] The embodiments of the present invention are described with reference to the flowcharts and / or block diagrams of the methods, terminal devices (systems), and computer program products according to the embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of the processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, an embedded processor, or other programmable data processing terminal device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing terminal device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0200] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing terminal device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 These computer program instructions can also be loaded onto a computer or other programmable data processing terminal device, so that a series of operation steps are executed on the computer or other programmable terminal device to produce a computer-implemented process, so that the instructions executed on the computer or other programmable terminal device provide for implementing the process in the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0201] It should also be noted that, in this article, the terms "include", "comprises" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or terminal device including a series of elements includes not only those elements, but also includes other elements not explicitly listed, or also includes elements inherent to such process, method, article or terminal device. In the absence of further restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, article or terminal device including the elements.

[0202] Finally, it should be noted that the above is a preferred embodiment of the present invention. It should be pointed out that although the preferred embodiment of the present invention has been described, for those skilled in the art, once the basic creative concept of the present invention is known, several improvements and modifications can be made without departing from the principles of the present invention, and these improvements and modifications should also be regarded as the protection scope of the present invention. Therefore, the attached claims are intended to be interpreted as including the preferred embodiment and all changes and modifications that fall within the scope of the embodiments of the present invention.

Claims

1. A microwave detection performance simulation analysis method, characterized in that: The following steps are involved: Step S1, calculating the detection parameters of the microwave detection payload and the antenna beam coverage; Step S2, simulating and calculating the radar echo of the target according to the relative position and attitude relationship between the microwave detection payload and the target; Step S3, calculating the scattered echo of each clutter block to generate environmental clutter, specifically including: Step S31, establishing a calculation model for the surface clutter backscattering coefficient of the microwave detection payload irradiation area according to the ground type; Step S32, calculating the area and number of clutter blocks according to the distance relationship between the target and the radar and the resolution of the microwave detection payload, specifically includes: Step S321, completing distance grid division; The distance resolution of the microwave detection payload on the ground is used to divide the distance units. The distance resolution is calculated using the following formula: Where ψ is the angle between the incident direction of the radar pulse and the ground; R Within the range, two different targets can be distinguished by the microwave detection payload. In the range direction, the main beam coverage area is divided into multiple equally spaced range rings. The distance from any range ring to the microwave detection payload is expressed as R. i express; Step S322, completing azimuth grid division; According to the Mitchell model, the azimuth direction is divided according to the angular resolution of the radar, f r represents the pulse repetition frequency, and the Doppler bandwidth of the microwave detection payload is represents the elevation angle of the microwave detection payload radar antenna, and the Doppler resolution is f r / M, then in the azimuth direction, the number of units is the number of clutter ring units, and we have: K=Δf d M / f r The azimuth division size is: The area of ​​the clutter ring is calculated as: Step S323, calculating the distance blur quantity; The radar transmits a pulse train at a fixed frequency, and any pulse train forms a series of circular areas on the ground. According to the period of the pulse signal, the maximum unambiguous distance R is calculated u : Among them, T p Represents the pulse period, which is the same as the pulse repetition frequency f r The two are reciprocal; Assuming that the grazing angles in the radar main beam are approximately equal, the amount of range ambiguity in the main lobe is expressed as: Step S33: Calculate the Doppler shift of the clutter block, and calculate the reflected echo s of each clutter block according to the area, orientation and back reflection coefficient of the clutter block. i (k); Step S34: Calculate the echo parameters of the environmental clutter according to the number of clutter blocks and the reflected echo of each clutter block. The process of calculating the echo parameters of the environmental clutter is as follows: Step S4, performing spatiotemporal adaptive processing on the target, environment and interference echoes, and calculating microwave detection effectiveness parameters according to the target fluctuation type and spatiotemporal processing gain; The step S4 specifically includes: Step S41, setting the number, power and movement of interference sources, and calculating the power reaching the position of the microwave detection payload antenna according to the power attenuation equation; Step S42, calculating the space-time steering vector and generating a space-time adaptive processing matrix; Step S43, calculating the target detection probability according to the target fluctuation type and the signal-to-noise ratio after spatiotemporal processing; Specifically include: Step S421: Calculate the spatial steering vector In the microwave detection payload coordinate system, the microwave detection payload flies along the X-axis, the array antenna is distributed along the X-axis, and the position of the first antenna unit is taken as the coordinate origin. The coordinate position of the i-th channel antenna is expressed as: coor i rad =((i-1)d,0,0) For the position at a pitch angle of θ EL , the azimuth angle is θ AZ The target at the position is represented by coor=R(sinθ EL cosθ AZ , sinθ EL sinθ AZ , cosθ EL ), its vector distance to the ith channel antenna is expressed as Dis i rad =(Rsinθ EL cosθ AZ -(i-1)d,Rsinθ EL sinθ AZ ,Rcosθ EL ); The vector difference between the i-th channel antenna and the first vector antenna is expressed as: Dis i rad -Dis0 rad =((i-1)d,0,0) The projection of the i-th channel antenna and the first vector antenna in the target direction is expressed as: (i-1)dsinθ EL cosθ AZ The phase difference between the i-th channel antenna and the first vector antenna is expressed as: The normalized angle is defined as: Then the spatial steering vector is expressed as: Step S422, calculating the time-oriented vector; Assuming the target is stationary, the Doppler frequency is mainly determined by the relative motion speed between the target and the microwave detection payload in the radial direction. The projection of the velocity vector on the distance vector between the target and the radar is the radial velocity. The velocity vector is along the X-axis in the satellite coordinate system, so the velocity vector is: V rad =(V,0,0) The radial velocity is: v r rad =V*sinθ EL cosθ AZ The Doppler frequency is: Then the Doppler time steering vector is: Step S423: Generate a space-time covariance matrix according to the target, environment and interference power at the microwave detection payload, and obtain a space-time adaptive processing matrix by using space-time optimization processing.

2. The microwave detection effectiveness simulation analysis method according to claim 1, characterized in that: The step S1 specifically includes: Step S11, according to the load platform parameters of the microwave detection payload, determine the current velocity vector and height parameters of the load platform, and determine the geometric coordinate relationship between the detection area and the load platform; Step S12: Modeling the three-dimensional radiation pattern of the microwave detection payload according to the antenna parameters of the microwave detection payload to obtain the antenna radiation pattern of the microwave detection payload that satisfies the main lobe beam width; Step S13: Calculate the beam coverage of the main beam according to the geometric coordinate relationship between the detection area and the payload platform and the antenna pattern.

3. The microwave detection effectiveness simulation analysis method according to claim 1, characterized in that: The step S2 specifically includes: Step S21, calculating the radar cross-section data of the target according to the material parameters and shape parameters of the target, and determining the scattering type; Step S22: Calculate the echo amplitude of the backscattered echo of the target according to the radar cross section data of the target and the radar echo equation.

4. The microwave detection effectiveness simulation analysis method according to claim 3 is characterized in that: The step S22 specifically includes: According to the attitude of the microwave detection payload and its relative position to the target, the antenna gain is obtained through the antenna pattern model. The RCS in the corresponding direction is obtained through the scattering cross-section data of the target, that is, The echo amplitude of the target's backscattered echo is calculated according to the radar equation. The formula is: Among them, P av is the average transmitted power, is the antenna gain in the detection direction in the detection payload coordinate system, θ rad , represents the pitch angle and azimuth angle of the target in the detection payload coordinate system, λ is the wavelength, is the RCS of the target in its local coordinate system in the direction of the microwave detection load, k is the Boltzmann constant, T0 is the equivalent temperature, B n is the Doppler bandwidth, F n is the receiver noise figure, L s is the system loss.

5. An electronic device, characterized in that: include: One or more processors, one or more memories, and one or more computer programs; wherein the processor is connected to the memory, and the one or more computer programs are stored in the memory, and when the electronic device is running, the processor executes the one or more computer programs stored in the memory, so that the electronic device performs the microwave detection effectiveness simulation analysis method as described in any one of claims 1 to 4.

6. A computer-readable storage medium, characterized in that: Used to store computer instructions, which, when executed by a processor, implement the microwave detection effectiveness simulation analysis method as described in any one of claims 1 to 4.

Images