A three-dimensional real-time characterization method for dynamic target electromagnetic scattering characteristics of a high-speed aircraft

Abstract

The present application relates to a kind of high-speed aircraft dynamic target electromagnetic scattering characteristic three-dimensional real-time characterization method, belong to the field of digital processing, method includes: aircraft trajectory is discretized and posture acquisition, and the target phase position is calculated and coordinate conversion;Establish spherical scattering grid, respectively calculate the scattering field in the vertical polarization and horizontal polarization two polarization incidence conditions for each grid direction;Get the incident electromagnetic field on the target position and calculate the incident electromagnetic field on the target position by the scattering electromagnetic field of target in the local coordinate system direction of aircraft;After the vertical polarization component and horizontal polarization of scattering electromagnetic field are respectively valued, generate three-dimensional scattering direction diagram.The present application can adapt to various high dynamic flight trajectory and attitude change, provide high-fidelity electromagnetic environment support for hypersonic target identification, stealth design and radar echo simulation.

Description

Technical Field

[0001] This invention relates to the field of electronic digital processing, and in particular to a three-dimensional real-time characterization method for the electromagnetic scattering characteristics of dynamic targets of high-speed aircraft. Background Technology

[0002] Current characterization of the electromagnetic scattering characteristics of high-dynamic aircraft targets mainly relies on three types of technologies: static RCS (Radar Cross Section) databases, dynamic coordinate rotation correction, and simplified motion models. Static database technology (such as the AN / APG-81 radar system) presets scattering parameters under fixed coordinates. However, the local coordinate system of high-speed aircraft (Mach number > 5) changes by more than 10° within milliseconds, leading to severe inaccuracies in the scattering model. For example, in a hypersonic vehicle dive test, the measured RCS deviated by 15dB from the static database, triggering target misjudgment. Dynamic correction methods (such as the Meteor missile seeker algorithm) transform static database data using an Euler rotation matrix, resulting in computational complexity as high as... These existing technologies cannot meet the 10ms real-time response requirement and ignore changes in electromagnetic wave propagation paths (such as distance propulsion phase delay). In maneuvering states, the azimuth error is >5°, and the root mean square error of RCS prediction is ≥8dB. Simplified models (such as the FEKO motion module) sacrifice polarization integrity, causing the identification system to be unable to distinguish between real targets and decoys. The root causes of these defects are: the inherent conflict between the static nature of the coordinate system and high-speed motion; the incompatibility between real-time dynamic correction and accuracy; and the oversimplification of polarization characteristics. Therefore, addressing the problems of existing technologies is a pressing issue. Summary of the Invention

[0003] The purpose of this invention is to overcome the shortcomings of the prior art and provide a three-dimensional real-time characterization method for the electromagnetic scattering characteristics of dynamic targets of high-speed aircraft, thus solving the deficiencies of the prior art.

[0004] The objective of this invention is achieved through the following technical solution: a three-dimensional real-time characterization method for the electromagnetic scattering characteristics of dynamic targets of high-speed aircraft, the method comprising:

[0005] Step 1: Discretize the aircraft trajectory and acquire its attitude, and calculate and transform the target phase position;

[0006] Step 2: Establish a spherical scattering grid and calculate the scattering field for each grid direction under both vertical and horizontal polarization incident conditions.

[0007] Step 3: Obtain the incident electromagnetic field at the target location and calculate the scattered electromagnetic field of the incident electromagnetic wave at the target location via the target in the local coordinate system of the aircraft.

[0008] Step 4: After assigning values ​​to the vertical polarization component and the horizontal polarization of the scattered electromagnetic field, a three-dimensional scattering pattern is generated.

[0009] The discretization of the aircraft trajectory and attitude acquisition include:

[0010] Discretized flight trajectory: total duration By interval Divided into k time points, denoted as Obtain the time of flight trajectory of the aircraft Location = Vector information;

[0011] Attitude acquisition: Attitude information includes velocity direction Scrolling direction and pitch direction Three unit vectors.

[0012] The calculation and coordinate transformation of the target phase position includes:

[0013] Traverse the time series to obtain The location of the spacecraft at that moment is The projection point is At this time, the local coordinate system of the aircraft is ;

[0014] The target position is calculated based on the Cartesian coordinates in the aircraft's local coordinate system. ,in, Let be the coordinate components of the target position point relative to the aircraft's position vector on the velocity direction axis of the aircraft's local coordinate system. Let be the coordinate components of the target position point relative to the aircraft's position vector on the roll direction axis of the aircraft's local coordinate system. The coordinate components of the target position point relative to the aircraft's position vector on the pitch axis of the aircraft's local coordinate system;

[0015] Cartesian coordinates Convert to spherical coordinates And calculate the direction vector. , This represents the actual straight-line distance from the target to the aircraft. This indicates the angle between the target and the aircraft's pitch direction. This represents the angle between the target's projection point on the horizontal plane and the direction of the aircraft's velocity. Represents the radial vector. Represents the pitch direction vector. Represents the azimuth direction vector;

[0016] Generate M The angle is evenly distributed in and N The angle is evenly distributed in The spherical angle grid, its first indivual The j-th angle The direction angle is ,in, The pitch angle, It is the azimuth angle;

[0017] Generation dimension is 3D mesh The polarization channels used to record the four scattered field components on each spherical angle grid are respectively ,in, The polarization channel represents the vertically polarized incident and vertically polarized received scattered field components. The polarization channel represents the scattered field components incident vertically and received horizontally. The polarization channel represents the scattered field components incident horizontally and received vertically. The polarization channel represents the horizontally polarized incident and horizontally polarized received scattered field components.

[0018] Step two specifically includes:

[0019] Traverse the spherical angle grid to obtain the first... indivual The j-th angle The direction angle is Based on the local coordinate system of the aircraft From the origin in this direction Emitting electromagnetic waves ,in, This represents the vertically polarized wave component of the incident electromagnetic wave. Represents the horizontally polarized wave component of the incident electromagnetic wave, where, Represents the reference for the vertical polarization direction. Represents the reference for the horizontal polarization direction;

[0020] set up , Calculate according to the content of step three. Position and direction Vertical polarization component of the scattered field and horizontal polarization components .

[0021] Step three specifically includes:

[0022] A1. Plane electromagnetic waves For the direction of propagation Distance to advance to target position The incident electromagnetic field at the target location is obtained. ,in, The direction of propagation of the incident electromagnetic wave, For incident electromagnetic waves from the aircraft position Along direction propagate to target location The propagation distance;

[0023] A2. Calculate the incident electromagnetic field at the target position T where the incident electromagnetic wave arrives. via the target in the aircraft's local coordinate system of Scattered electromagnetic field in the direction ,in, Indicates the direction of propagation of the scattered wave.

[0024] The steps in A2 include:

[0025] A201. Calculate the local coordinate system of Spherical coordinates in the global coordinate system of orientation ,in, Indicates the same scattering direction The spherical coordinate angles obtained after transformation to a fixed global coordinate system;

[0026] A202, Extracting spherical coordinates The three nearest coordinate values ​​form the Lagrange coordinate interpolation grid on the sphere;

[0027] A203. Obtain spherical coordinates by interpolating the node values ​​on the interpolation grid. The four scattered field components are VV, VH, HV, and HH. ;

[0028] A204. Calculate the incident electromagnetic wave in the local coordinate system. of Incident direction in the global coordinate system of the direction ,in, Represents the local coordinate system of the aircraft The following describes the propagation direction of the incident electromagnetic wave. Indicates the local incident direction The spherical coordinate angle obtained after transformation to a fixed global coordinate system;

[0029] A205, Solving for the incident electromagnetic field Vertical polarization components in the world coordinate system and horizontal polarization components ;

[0030] A206, Calculation of incident electromagnetic field exist Scattered electromagnetic field in the direction Vertical polarization component and horizontal polarization components ;

[0031] A207. Synthesizing vertical polarization components in the world coordinate system. and horizontal polarization components The scattered electromagnetic field at the target location is obtained. ;

[0032] A208. When the scattered electromagnetic field at the target location Advance to At position point, the scattered electromagnetic field is scatter electromagnetic field Transform from world coordinate system to local coordinate system exist Vertical polarization component of the scattered field in the direction and horizontal polarization components .

[0033] Step four includes:

[0034] The vertical polarization component of the scattered electromagnetic field and horizontal polarization components Assign values ​​to respectively ,in, in, Representing a 3D mesh China direction ( , The vector consisting of the four polarization scattering field components corresponding to the given )

[0035] set up , Recalculate according to the content of step three. Position and direction Vertical polarization component of the scattered field and horizontal polarization components ;

[0036] The vertical polarization component of the scattered field and horizontal polarization components Assign values ​​to respectively ;

[0037] Collect M Angles and N on the corner Where i=0,1,…,M-1,j=0,1,…,N-1,l=0,1,2,3, four sets of data are obtained as follows: , , , Then we get , , and The scattering pattern distribution, where, Representing a 3D mesh China direction ( , The upper polarization channel number is The scattered field components, =0, 1, 2, 3 respectively correspond to , , , Four polarization channels, = Indicates direction ( , )superior The scattered field components of the polarization channel, = Indicates direction ( , )superior The scattered field components of the polarization channel, = Indicates direction ( , )superior The scattered field components of the polarization channel, = Indicates direction ( , )superior Scattered field components of the polarization channel.

[0038] The present invention has the following advantages:

[0039] 1. High accuracy and adaptability of dynamic scattering characterization: This invention accurately describes the changes in the relative geometric relationship between the target and radar during high-speed flight by establishing a time-varying local coordinate system synchronized with the aircraft's motion state in real time. Utilizing spherical coordinate mapping and vector decomposition, it effectively overcomes the failure problem of traditional static scattering models in high-speed scenarios, significantly improving the accuracy of dynamic RCS simulation. This method can adapt to various highly dynamic flight trajectories and attitude changes, providing high-fidelity electromagnetic environment support for hypersonic target identification, stealth design, and radar echo simulation.

[0040] 2. Excellent computational efficiency and ease of engineering implementation: By discretizing the time series and using a spherical angle grid, the continuous dynamic problem is transformed into a sequential local static calculation, significantly reducing computational complexity. Combined with the interpolation reuse of the target's inherent scattering parameters (S-matrix), repeated full-wave simulations are avoided, significantly improving computational efficiency while maintaining accuracy. This invention has a clear structure and strong modularity, making it easy to integrate into existing radar simulation systems. It supports large-scale parameter scanning and real-time simulation, making it suitable for rapid evaluation and optimization design in engineering practice.

[0041] 3. Supports full polarization scattering analysis with clear physical mechanisms: Employing dual-polarization incident and received modes, it fully acquires the four polarization scattering components (VV, VH, HV, and HH), comprehensively characterizing the polarization scattering characteristics of the target. Through rigorous path phase correction and coordinate system transformation, it clearly demonstrates the physical process of electromagnetic wave propagation, combining the robustness of numerical calculations with the transparency of physical meaning, providing a reliable tool for studying the electromagnetic properties of complex targets.

[0042] 4. High scalability and compatibility with existing databases and multiple algorithms: The framework of this invention is highly open and can be easily integrated with existing target RCS databases or various electromagnetic calculation engines (such as PO, MLFMM, etc.), exhibiting good compatibility and scalability. Its established dynamic mapping and interpolation mechanism also lays the foundation for the subsequent introduction of multiple scattering models, non-uniform targets, environmental interference, and other factors, making it suitable for a wider range of electromagnetic scattering and adversarial simulation scenarios. Attached Figure Description

[0043] Figure 1 This is a schematic diagram of the discretization of the aircraft's flight trajectory.

[0044] Figure 2 A schematic diagram defining a local coordinate system.

[0045] Figure 3 This is a schematic diagram of the spherical angle grid distribution.

[0046] Figure 4 This is a schematic diagram of the scattering field component matrix structure. Detailed Implementation

[0047] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of this application, and not all of the embodiments. The components of the embodiments of this application described and shown in the accompanying drawings can generally be arranged and designed in various different configurations. Therefore, the detailed description of the embodiments of this application provided below with reference to the accompanying drawings is not intended to limit the scope of protection of the claimed application, but merely represents selected embodiments of this application. All other embodiments obtained by those skilled in the art based on the embodiments of this application without inventive effort are within the scope of protection of this application. The present invention will be further described below with reference to the accompanying drawings.

[0048] This invention specifically relates to a three-dimensional real-time characterization method for the electromagnetic scattering characteristics of dynamic targets of high-speed aircraft. It addresses the mismatch between static models and high-speed motion by reconstructing a time-varying local coordinate system, fully compensates for propagation effects through accurate path modeling, reduces the query complexity of the global dispersion parameter library by utilizing spherical Lagrange interpolation, and independently calculates and stores four sets of dispersion components through a fully polarized independent solution mechanism. The specific method includes the following steps:

[0049] S1, Discretized flight trajectory: such as Figure 1 As shown, each moment in the flight trajectory of the aircraft is obtained. Each location point = Vector information:

[0050] ,

[0051] in , This represents a discrete time series, where k is the series index and K is the length of the time series. The time sampling interval is... Total duration.

[0052] for The position of the spacecraft on its trajectory at any given time is denoted as... , for The projection point on the XOY two-dimensional plane is denoted as . Let T be the altitude of the aircraft, and T be the target location point, i.e., the location of the target, represented as... ;

[0053] The tangential vector of the aircraft trajectory is in the direction of the flight velocity, denoted as: , The direction of the aircraft's roll angle is represented as , The pitch direction of the aircraft is represented as ,in, The aircraft are in At any given moment, its velocity direction is a unit vector. Cosines (components) of the three directions in the global coordinate system; The aircraft are in At any given moment, the unit vector of the rolling direction Component coordinates in the global coordinate system; The aircraft are in At any given moment, the unit vector in the pitch direction Component coordinates in the global coordinate system.

[0054] S2, such as Figure 2 As shown, by traversing the time series, we obtain... The location of the spacecraft at that moment is The projection point is At this time, the local coordinate system of the aircraft is The target position is calculated based on spherical coordinates in the local coordinate system of the aircraft. and the corresponding direction vector .

[0055] S201. Calculate the target position based on the rectangular coordinates in the local coordinate system of the aircraft. :

[0056] ,

[0057] ,

[0058] ;

[0059] in, The coordinate components of the target position point relative to the aircraft's position vector on the velocity direction axis (nose direction axis) of the aircraft's local coordinate system. The coordinate components of the target position point relative to the aircraft's position vector on the roll direction axis (lateral axis) of the aircraft's local coordinate system. The coordinate components of the target position point relative to the aircraft's position vector on the pitch axis (back axis) of the aircraft's local coordinate system; the directions of the velocity axis, roll axis, and pitch axis are respectively related to the unit velocity vector. , Rolling direction unit vector Pitch direction unit vector They are in the same direction.

[0060] S202, convert rectangular coordinates Convert to spherical coordinates :

[0061] ,

[0062] ,

[0063] ,

[0064] ;

[0065] in, This represents the actual straight-line distance from the target to the aircraft (the origin of the local coordinate system); Indicates the pitch direction of the target relative to the aircraft. The angle between the axis (i.e., the direction of the fuselage back); This represents the projection point of the target onto the horizontal plane, relative to the direction of the aircraft's velocity. The angle between the axis (i.e., the direction of the machine head); This indicates that the target's position is on the horizontal plane of the aircraft (i.e., in the direction of velocity). and scrolling direction The projected distance on the plane formed.

[0066] S203, Calculate the direction vector :

[0067] ,

[0068] ,

[0069] ;

[0070] in, This represents the radial vector, pointing in the direction of the target. It represents the pitch direction vector, indicating changes in the vertical direction; It represents the azimuth direction vector, indicating the change in the horizontal direction.

[0071] S204, Generate M items The angle is evenly distributed in and N The angle is evenly distributed in The spherical angle grid; its first indivual The j-th angle The direction angle is ;in, The pitch angle represents the angle between the grid sampling direction and the pitch direction axis (i.e., the dorsal direction axis) of the aircraft's local coordinate system, and its value ranges from [0,π]. The azimuth angle represents the angle between the projection of the grid sampling direction onto the horizontal plane of the aircraft (the plane formed by the velocity direction axis and the roll direction axis) and the velocity direction axis of the aircraft (i.e. the nose direction axis), with a value range of [0, 2π]. For the i-th pitch angle sample value, For the j-th azimuth angle sample value, and Together, determine the sampling direction of the (i,j)th spherical grid.

[0072] S205, Generation dimension is 3D mesh The polarization channels used to record the four scattered field components on each spherical angle grid are respectively ,in, The polarization channel represents the vertically polarized incident and vertically polarized received scattered field components. The polarization channel represents the scattered field components incident vertically and received horizontally. The polarization channel represents the scattered field components incident horizontally and received vertically. The polarization channel represents the horizontally polarized incident and horizontally polarized received scattered field components. It is a three-dimensional mesh with dimensions of That is, including × There are scattering directions, and each scattering direction records 4 polarization scattering field components. The total number of grid elements is . ; Representing a 3D mesh China direction ( , The vector formed by the four polarization scattering field components corresponding to ) has a direction ( , The four scattered field components recorded on the screen.

[0073] S3, such as Figure 3 As shown, by traversing the spherical angle grid, the first... indivual The j-th angle The direction angle is Based on the local coordinate system of the aircraft From the origin in this direction Emitting electromagnetic waves , Represented as:

[0074] ,

[0075] in, , ;

[0076] This represents the vertically polarized wave component of the incident electromagnetic wave. This represents the horizontally polarized wave component of the incident electromagnetic wave. For free space wavenumber, ; and This represents the unit vector in the local spherical coordinate system calculated in real time based on the aircraft's current attitude and the target orientation. Represents the reference for the vertical polarization direction. It represents the reference for the horizontal polarization direction.

[0077] S4, Settings , Calculate according to the steps S5-S6 below. Position and direction Vertical polarization component of the scattered field and horizontal polarization components ;

[0078] S5, this plane electromagnetic wave uses For the direction of propagation Distance to advance to target position The incident electromagnetic field at the target location is obtained. Represented as:

[0079] ,

[0080] in, The direction of propagation of the incident electromagnetic wave, For incident electromagnetic waves from the aircraft position Along direction propagate to target location propagation distance

[0081] S6, such as Figure 4 As shown, calculate the incident electromagnetic field at the target position T where the incident electromagnetic wave arrives. via the target in the aircraft's local coordinate system of Scattered electromagnetic field in the direction .

[0082] S601, Calculate the local coordinate system of Spherical coordinates in the global coordinate system of orientation Then we have:

[0083] ,

[0084] in, The local coordinate system of the aircraft describes the direction of the scattered wave propagation. Indicates the same scattering direction The spherical coordinate angle is obtained after transformation to a fixed global coordinate system.

[0085] Then we have:

[0086] ,

[0087] ,

[0088] in, , , ;

[0089] S602, Extracting spherical coordinates The three nearest coordinate values ​​form the Lagrange coordinate interpolation grid on the sphere;

[0090] S603. Obtain spherical coordinates by interpolating the node values ​​on the interpolation grid. The four scattered field components are VV, VH, HV, and HH. ;

[0091] S604. Calculate the incident electromagnetic wave in the local coordinate system. of Incident direction in the global coordinate system of the direction Then we have:

[0092] ,

[0093] in, Represents the local coordinate system of the aircraft The propagation direction of the incident electromagnetic wave described below; This indicates the aforementioned local incident direction. The spherical coordinate angle obtained after transformation to a fixed global coordinate system.

[0094] Then we have:

[0095] ,

[0096] ,

[0097] in, , , ;

[0098] S605, Solving for the incident electromagnetic field Vertical polarization components in the world coordinate system and horizontal polarization components :

[0099] ,

[0100] ,

[0101] in, ; ;

[0102] S606, Calculation of incident electromagnetic field exist Scattered electromagnetic field in the direction Vertical polarization component Horizontal polarization component ;

[0103] S607. Synthesizing vertical polarization components in the world coordinate system. and horizontal polarization components The scattered electromagnetic field at the target location is obtained. Represented as:

[0104] ,

[0105] in, , , ;

[0106] S608, Scattered electromagnetic field at the target location Advance to Location point, scattered electromagnetic field Represented as:

[0107] ,

[0108] in, .

[0109] S609, Scattering electromagnetic field Transform from world coordinate system to local coordinate system exist Vertical polarization component of the scattered field in the direction and horizontal polarization components :

[0110] ,

[0111] ,

[0112] in, , .

[0113] S7. The vertical polarization component of the scattered electromagnetic field and horizontal polarization components Assign values ​​to respectively :

[0114] ,

[0115] ;

[0116] S8, Settings , Recalculate according to steps S5-S6. Position and direction Vertical polarization component of the scattered field and horizontal polarization components ;

[0117] S9. The vertical polarization component of the scattered field and horizontal polarization components Assign values ​​to respectively :

[0118] ,

[0119] .

[0120] in, Represents a grid China direction ( , The vector formed by the four polarization scattering field components corresponding to ) , These are the 0th and 1st elements of the vector, respectively, representing the direction ( , )superior polarization channels and The scattered field component of the polarization channel, 20log represents the logarithmic gain operation with base 10, converting the scattered field amplitude into a value in dB.

[0121] S10, Collect M items Angles and N on the corner Where i=0,1,…,M-1,j=0,1,…,N-1,l=0,1,2,3, four sets of data are obtained as follows: , , , Then you can get , , and The distribution of the scattering pattern.

[0122] in, Representing a 3D mesh China direction ( , The upper polarization channel number is The scattered field components, These are the pitch angle numbers. This is the azimuth number. This is the polarization channel number. =0, 1, 2, 3 respectively correspond to , , , Four polarization channels; = Indicates direction ( , )superior The scattered field components of the polarization channel, = Indicates direction ( , )superior The scattered field components of the polarization channel, = Indicates direction ( , )superior The scattered field components of the polarization channel, = Indicates direction ( , )superior The scattered field components of the polarization channel, the four sets of data are all × The set of scattered field components of the four polarization channels in each scattering direction is used to generate the scattering pattern of each polarization channel.

[0123] The above description is merely a preferred embodiment of the present invention. It should be understood that the present invention is not limited to the forms disclosed herein and should not be construed as excluding other embodiments. It can be used in various other combinations, modifications, and improvements, and can be altered within the scope of the concept described herein through the above teachings or related technologies or knowledge. Modifications and variations made by those skilled in the art that do not depart from the spirit and scope of the present invention should be within the protection scope of the appended claims.

Claims

1. A three-dimensional real-time characterization method for the electromagnetic scattering characteristics of dynamic targets of high-speed aircraft, characterized in that: The method includes: Step 1: Discretize the aircraft trajectory and acquire its attitude, and calculate and transform the target phase position; Step 2: Establish a spherical scattering grid and calculate the scattering field for each grid direction under both vertical and horizontal polarization incident conditions. Step 3: Obtain the incident electromagnetic field at the target location, and calculate the scattered electromagnetic field of the incident electromagnetic wave at the target location via the target in the local coordinate system of the aircraft. Step 4: After assigning values ​​to the vertical polarization component and the horizontal polarization of the scattered electromagnetic field, a three-dimensional scattering pattern is generated.

2. The three-dimensional real-time characterization method for the electromagnetic scattering characteristics of dynamic targets of high-speed aircraft according to claim 1, characterized in that: The discretization of the aircraft trajectory and attitude acquisition include: Discretized flight trajectory: total duration By interval Divided into k time points, denoted as Obtain the time of flight trajectory of the aircraft Location = Vector information, and These represent the x-axis and y-axis coordinates, respectively. Indicates the altitude of the aircraft; Attitude acquisition: Attitude information includes velocity direction Scrolling direction and pitch direction Three unit vectors.

3. The three-dimensional real-time characterization method for the electromagnetic scattering characteristics of dynamic targets of high-speed aircraft according to claim 2, characterized in that: The calculation and coordinate transformation of the target phase position includes: Traverse the time series to obtain The location of the spacecraft at that moment is The projection point is At this time, the local coordinate system of the aircraft is ; The target position is calculated based on the Cartesian coordinates in the aircraft's local coordinate system. ,in, Let be the coordinate components of the target position point relative to the aircraft's position vector on the velocity direction axis of the aircraft's local coordinate system. Let be the coordinate components of the target position point relative to the aircraft's position vector on the roll direction axis of the aircraft's local coordinate system. The coordinate components of the target position point relative to the aircraft's position vector on the pitch axis of the aircraft's local coordinate system; Cartesian coordinates Convert to spherical coordinates And calculate the direction vector. , This represents the actual straight-line distance from the target to the aircraft. This indicates the angle between the target and the aircraft's pitch direction. This represents the angle between the target's projection point on the horizontal plane and the direction of the aircraft's velocity. Represents the radial vector. Represents the pitch direction vector. Represents the azimuth direction vector; Generate M The angle is evenly distributed in and N The angle is evenly distributed in The spherical angle grid, its first indivual The j-th angle The direction angle is ,in, The pitch angle, It is the azimuth angle; Generation dimension is 3D mesh The polarization channels used to record the four scattered field components on each spherical angle grid are respectively ,in, The polarization channel represents the vertically polarized incident and vertically polarized received scattered field components. The polarization channel represents the scattered field components incident vertically and received horizontally. The polarization channel represents the scattered field components incident horizontally and received vertically. The polarization channel represents the horizontally polarized incident and horizontally polarized received scattered field components.

4. The three-dimensional real-time characterization method for the electromagnetic scattering characteristics of dynamic targets of high-speed aircraft according to claim 3, characterized in that: Step two specifically includes: Traverse the spherical angle grid to obtain the first... indivual The j-th angle The direction angle is Based on the local coordinate system of the aircraft From the origin in this direction Emitting electromagnetic waves ,in, This represents the vertically polarized wave component of the incident electromagnetic wave. Represents the horizontally polarized wave component of the incident electromagnetic wave, where, Represents the reference for the vertical polarization direction. Represents the reference for the horizontal polarization direction; set up , Calculate according to the content of step three. Position and direction Vertical polarization component of the scattered field and horizontal polarization components .

5. The three-dimensional real-time characterization method for the electromagnetic scattering characteristics of dynamic targets of high-speed aircraft according to claim 4, characterized in that: Step three specifically includes: A1. Plane electromagnetic waves For the direction of propagation Distance to advance to target position The incident electromagnetic field at the target location is obtained. ,in, The direction of propagation of the incident electromagnetic wave, For incident electromagnetic waves from the aircraft position Along direction propagate to target location The propagation distance; A2. Calculate the incident electromagnetic field at the target position T where the incident electromagnetic wave arrives. via the target in the aircraft's local coordinate system of Scattered electromagnetic field in the direction ,in, Indicates the direction of propagation of the scattered wave.

6. The three-dimensional real-time characterization method for the electromagnetic scattering characteristics of dynamic targets of high-speed aircraft according to claim 5, characterized in that: The steps in A2 include: A201. Calculate the local coordinate system of Spherical coordinates in the global coordinate system of orientation ,in, Indicates the same scattering direction The spherical coordinate angles obtained after transformation to a fixed global coordinate system; A202, Extracting spherical coordinates The three nearest coordinate values ​​form the Lagrange coordinate interpolation grid on the sphere; A203. Obtaining spherical coordinates by interpolating the node values ​​on the interpolation grid. The four scattered field components corresponding to VV, VH, HV, and HH are: ,in, The interpolation results are respectively , , , The numerical values ​​of the scattered field components of the four polarization channels; A204. Calculate the incident electromagnetic wave in the local coordinate system. of Incident direction in the global coordinate system of the direction ,in, Represents the local coordinate system of the aircraft The following describes the propagation direction of the incident electromagnetic wave. Indicates the local incident direction The spherical coordinate angle obtained after transformation to a fixed global coordinate system; A205, Solving for the incident electromagnetic field Vertical polarization components in the world coordinate system and horizontal polarization components ; A206, Calculation of incident electromagnetic field exist Scattered electromagnetic field in the direction Vertical polarization component and horizontal polarization components ; A207. Synthesizing vertical polarization components in the world coordinate system. and horizontal polarization components The scattered electromagnetic field at the target location is obtained. ; A208. When the scattered electromagnetic field at the target location Advance to At position point, the scattered electromagnetic field is scatter electromagnetic field Transform from world coordinate system to local coordinate system exist Vertical polarization component of the scattered field in the direction and horizontal polarization components .

7. The three-dimensional real-time characterization method for the electromagnetic scattering characteristics of dynamic targets of high-speed aircraft according to claim 6, characterized in that: Step four includes: The vertical polarization component of the scattered electromagnetic field and horizontal polarization components Assign values ​​to respectively ,in, Representing a 3D mesh China direction ( , The vector consisting of the four polarization scattering field components corresponding to the given ) set up , Recalculate according to the content of step three. Position and direction Vertical polarization component of the scattered field and horizontal polarization components ; The vertical polarization component of the scattered field and horizontal polarization components Assign values ​​to respectively ; Collect M Angles and N on the corner Where i=0,1,…,M-1,j=0,1,…,N-1,l=0,1,2,3, four sets of data are obtained as follows: , , , Then we get , , and The scattering pattern distribution, where, Representing a 3D mesh China direction ( , The upper polarization channel number is The scattered field components, =0, 1, 2, 3 respectively correspond to , , , Four polarization channels, = Indicates direction ( , )superior The scattered field components of the polarization channel, = Indicates direction ( , )superior The scattered field components of the polarization channel, = Indicates direction ( , )superior The scattered field components of the polarization channel, = Indicates direction ( , )superior Scattered field components of the polarization channel.

Images