CN116125423A - Scattered field characterization method of electromagnetic target - Google Patents

Scattered field characterization method of electromagnetic target Download PDF

Info

Publication number
CN116125423A
CN116125423A CN202310079398.1A CN202310079398A CN116125423A CN 116125423 A CN116125423 A CN 116125423A CN 202310079398 A CN202310079398 A CN 202310079398A CN 116125423 A CN116125423 A CN 116125423A
Authority
CN
China
Prior art keywords
unit
triangle
vector
scattering
field
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202310079398.1A
Other languages
Chinese (zh)
Other versions
CN116125423B (en
Inventor
李尧尧
蔡少雄
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dongshen Electromagnetic Technology Chengdu Co ltd
Original Assignee
Dongshen Electromagnetic Technology Chengdu Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Dongshen Electromagnetic Technology Chengdu Co ltd filed Critical Dongshen Electromagnetic Technology Chengdu Co ltd
Priority to CN202310079398.1A priority Critical patent/CN116125423B/en
Publication of CN116125423A publication Critical patent/CN116125423A/en
Application granted granted Critical
Publication of CN116125423B publication Critical patent/CN116125423B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Aerials With Secondary Devices (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a scattered field characterization method of an electromagnetic target, which comprises the following steps: s1, giving a triangular mesh file of a target scatterer, wherein the triangular mesh file comprises a point list { P } of the target scatterer i Sum triangle cell list { T } j -a }; s2, applying frequency domain uniform plane electromagnetic waves to the target scatterer, and for { T } j J-th triangle unit T in } j Calculating the scattering power parameter Dv under the condition of vertical polarization j And scattering power parameter Dh under horizontal polarization conditions j And calculates j triangle units T j Is a scattering field of (2); s3, traversing a triangle unit list { T } j Each triangle element in the sequence and determines the entire target fringe field. The invention effectively obtains the scattering field distribution of the target scattering body, has simple calculation and does not need to consume a large amount of calculation resources.

Description

Scattered field characterization method of electromagnetic target
Technical Field
The invention relates to a scattered field, in particular to a scattered field characterization method of an electromagnetic target.
Background
Fringe field computation is a key element in the research of radar target characteristics, the radar scattering cross section of a radar target can be calculated by using the fringe field, and the parameter describes the physical characteristics of the radar target and is an important characteristic of electromagnetic scattering of the target in a high-frequency region. Along with the continuous improvement of the performance of a radar system and the gradual deepening of the understanding of researchers on electromagnetic scattering mechanisms, the stealth and recognition of radar targets are mutually independent, but the development of the technology of the radar system and the recognition of the radar targets are mutually promoted, the scattering characteristics of the targets can be controlled by recognizing the constraint rules of the scattering characteristics of the targets, and the stealth design and the target recognition of the radar targets are further serviced.
The conventional method for calculating the scattered field is to obtain a target within the range of the working frequency bandwidth through calculation electromagnetism, excite the target by an incident plane wave with a certain radar beam angle and receive the target with a certain radar beam angle, however, the method relates to the transmission and the reception of electromagnetic waves with a wide frequency band and a wide angle range, a large number of frequency sampling points and incidence angle sampling points are needed, and the electric size of the target to be calculated is very large, so that huge calculation resources are needed to be consumed for calculating the scattered field by adopting the conventional calculation electromagnetism.
Disclosure of Invention
The invention aims to overcome the defects of the prior art, and provides a scattered field characterization method of an electromagnetic target, which effectively obtains the scattered field distribution of a target scatterer, is simple to calculate and does not need to consume a large amount of calculation resources.
The aim of the invention is realized by the following technical scheme: a method of fringe field characterization of an electromagnetic target, comprising the steps of:
s1, giving a triangular mesh file of a target scatterer, wherein the triangular mesh file comprises a point list { P } of the target scatterer i Sum triangle cell list { T } j }；
S2, applying frequency domain uniform plane electromagnetic waves to the target scatterer, and for { T } j J-th triangle unit T in } j Calculating the scattering power parameter Dv under the condition of vertical polarization j And scattering power parameter Dh under horizontal polarization conditions j And calculates j triangle units T j Is a scattering field of (2);
s3, traversing a triangle unit list { T } j Each triangular cell in the sequence, the fields of the individual triangular cells are obtained in step S2 and the fields of the whole target scatterer are determined.
The beneficial effects of the invention are as follows: the invention considers the reason that the scattering field is generated on the surface of a complex target, and the surface induction field is displayed in the form of surface scalar visualization, thus, by aiming at the target scattering bodyThe planar electromagnetic wave is uniformly distributed in the frequency domain, and then the scattering capacity parameter Dv under the condition of vertical polarization is calculated j And scattering power parameter Dh under horizontal polarization conditions j And then determining the scattered field of the whole target scatterer according to the scattered fields of the triangular units, so that the scattered field distribution of the target scatterer is effectively obtained, the calculation is simple, and a large amount of calculation resources are not required.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a schematic diagram of a scatterer model corresponding to a triangle mesh file of a target scatterer;
FIG. 3 is a schematic diagram of the scalar field distribution of the surface Dv obtained in the example;
fig. 4 is a schematic diagram of the scalar field distribution of the surface Dh obtained in the example.
Detailed Description
The technical solution of the present invention will be described in further detail with reference to the accompanying drawings, but the scope of the present invention is not limited to the following description.
As shown in fig. 1, a method for characterizing a scattered field of an electromagnetic target includes the following steps:
s1, giving a triangular mesh file of a target scatterer, wherein the triangular mesh file comprises a point list { P } of the target scatterer i Sum triangle cell list { T } j }；
Point list { P } of the target scatterer i Comprises N in } p The point coordinates, i-th point P i Is marked P by the seat i :(P i .x,P i .y,P i .z)；
Triangle cell list { T } of the target scatterer j Comprises N in } e The three point numbers of the j-th triangle unit defined according to the right-hand rotation direction are denoted as (T) j .i 1 ,T j .i 2 ,T j .i 3 )；
When the coordinates of the kth point in the jth triangle unit are queried, the coordinates are first selected from the triangle unit list { T } j Find j-th in }Triangle unit data T j Then from T j The sequence number ip=t of the kth point obtained in (a) j .i k Finally from the point list { P i The coordinate of the kth point in the jth triangle unit is found to be (P) ip .x,P ip .y,P ip .z)。
S2, applying frequency domain uniform plane electromagnetic waves to the target scatterer, and for { T } j J-th triangle unit T in } j Calculating the scattering power parameter Dv under the condition of vertical polarization j And scattering power parameter Dh under horizontal polarization conditions j And calculates j triangle units T j Is a scattering field of (2);
the step S2 includes:
s201 for { T ] j J-th triangle unit T in } j Acquiring coordinates of a central point of a triangle unit and a normal vector of the triangle unit, calculating the area of the triangle unit and the vector of the area of the triangle unit, and constructing a j-th triangle unit T j Surface local coordinate system C of (2) j ；
S202, calculating the surface current distribution at the center position of the jth unit
Figure BDA0004066998790000021
S203, calculating the radiation density of a scattered field formed by the jth triangle unit, and calculating a scattering capacity parameter Dv under the condition of vertical polarization j And scattering power parameter Dh under horizontal polarization conditions j 。
Wherein, the step S201 includes:
a1, calculating coordinates of a central point of the triangle unit:
from the triangle element list { T j The sequence numbers of the three points in the j-th triangle element are i1 = T, respectively j .i 1 、i2＝T j .i 2 and i3＝Tj .i 3 From the list of points { P i The coordinates of the three points obtained in the process are respectively P i1 、P i2 and Pi3 Then there is a cell center point coordinate
Figure BDA0004066998790000031
Is that;
Figure BDA0004066998790000032
a2, calculating a normal vector of the triangle unit:
obtaining triangle unit list { T } j Coordinates of three points in the j-th triangle unit are respectively: p (P) i1 、P i2 and Pi3 Calculating unit normal vector of jth triangle unit
Figure BDA0004066998790000033
The method comprises the following steps:
Figure BDA00040669987900000314
wherein norm is a vector normalization function: norm (a) =a/len (a), len being a vector length function; x is the vector cross divisor;
a3, calculating the area of the triangle unit:
obtaining triangle unit list { T } j Coordinates of three points in the j-th triangle unit are respectively: p (P) i1 、P i2 and Pi3 Calculating the cell area S of the jth triangle cell j The method comprises the following steps:
S j ＝0.5len[(P i1 -P i2 )×(P i3 -P i2 )] (3)
a4, calculating a triangle area vector:
obtaining triangle unit list { T } j Coordinates of three points in the j-th triangle unit are respectively: p (P) i1 、P i2 and Pi3 Calculating the cell area vector S of the jth triangle cell j The method comprises the following steps:
S j ＝0.5[(P i1 -P i2 )×(P i3 -P i2 )] (4)
a5, establishing a surface local coordinate system:
surface localization of the jth triangle unitCoordinate system C j Comprising four elements: center point coordinates
Figure BDA0004066998790000034
Horizontal axis unit vector->
Figure BDA0004066998790000035
Vertical axis unit vector->
Figure BDA0004066998790000036
Height axis unit vector +.>
Figure BDA0004066998790000037
The surface part of the triangle unit can be simplified by the coordinate and vector transformation of the four elements of the surface part coordinate system, which C j The four quantities in (a) are calculated as follows:
Figure BDA0004066998790000038
Figure BDA0004066998790000039
Figure BDA00040669987900000310
Figure BDA00040669987900000311
wherein ,
Figure BDA00040669987900000312
is the center point coordinate of the triangle unit, +.>
Figure BDA00040669987900000313
Is from vertex P of triangle unit number i1 i1 Point to vertex number i 2P i2 The unit vector of the edge of (2); />
Figure BDA0004066998790000041
The normal vector is the outside unit normal vector of the plane where the triangle unit is located; />
Figure BDA0004066998790000042
Is->
Figure BDA0004066998790000043
and />
Figure BDA0004066998790000044
Forming right-hand orthogonal coordinate system, the right-hand sequence is +.>
Figure BDA0004066998790000045
Wherein, the step S202 includes:
let the electromagnetic field of the frequency domain planar electromagnetic wave be expressed as:
Figure BDA0004066998790000046
Figure BDA0004066998790000047
Figure BDA0004066998790000048
Figure BDA0004066998790000049
Figure BDA00040669987900000410
wherein ,k0 ＝2πf/c 0 Is free space wave number; c 0 Is the electromagnetic wave speed in free space; f (f)Is the working frequency; e (E) θ Is the vertically polarized component of the planar electromagnetic wave electric field,
Figure BDA00040669987900000411
for the horizontal polarization components of the planar electromagnetic field, θ and +.>
Figure BDA00040669987900000412
Is the incident angle of plane electromagnetic wave, +.>
Figure BDA00040669987900000413
Is the incident direction unit vector of the plane electromagnetic wave, +.>
Figure BDA00040669987900000414
Vertical polarization unit vector for incident electromagnetic wave, +.>
Figure BDA00040669987900000415
Horizontal polarization unit vector for incident electromagnetic wave, +.>
Figure BDA00040669987900000416
For complex number, η=120pi is free space electromagnetic wave impedance, E is electric field, H is magnetic field, r is field point coordinate vector; />
Figure BDA00040669987900000417
and />
Figure BDA00040669987900000418
A unit coordinate vector is a global coordinate system;
in the embodiment of the application, the frequency domain planar electromagnetic wave is used as an excitation source for solving a target scattered field, induced electromagnetic field distribution is formed on the surface of the target scattered field through excitation of the frequency domain planar electromagnetic wave, then the scattered field in a specified scattered direction is obtained through radiation integral calculation of the induced electromagnetic field distribution, and visual characterization is that integral terms of radiation integral calculation are characterized, wherein the integral terms describe the contribution strength of the scattered field in the specified scattered direction.
Article of clothingOptical approximation, induced current J at any point r on the area of the target scatterer illuminated by a uniform planar electromagnetic wave PO (r) is:
Figure BDA00040669987900000419
wherein ,
Figure BDA00040669987900000420
is the unit normal vector of the surface of the scattering body at the r position; x is a vector cross symbol;
coordinates of the center point of the jth triangle unit
Figure BDA00040669987900000421
Calculating the surface current distribution at the center of the jth cell as r +.>
Figure BDA00040669987900000422
Wherein, the step S203 includes:
the jth triangle element to the diffuser surface is at θ s And
Figure BDA00040669987900000423
analyzing the scattered field in the direction to obtain the far field contribution of the unit to the scattered field, namely the core characterization quantity Ra of electromagnetic scattering phenomenon j Expressed as:
Figure BDA00040669987900000424
Figure BDA00040669987900000425
Figure BDA0004066998790000051
Figure BDA0004066998790000052
wherein ω=2pi f, μ is the free space permeability,
Figure BDA0004066998790000053
is the scattering direction unit vector of the plane electromagnetic wave, +.>
Figure BDA0004066998790000054
and />
Figure BDA0004066998790000055
Is the vertical and horizontal polarization unit vector in the scattering direction of electromagnetic wave, R s Monitoring the radius of the sphere where the scattered field is located; ra (Ra) j Radiation factor for forming a scattering field for the j-th triangle element,>
Figure BDA0004066998790000056
is a green function;
the radiation density of the scattered field formed by the jth triangle unit is Da j ＝Ra j /S j Due to Da j Is an omnidirectional vector, and is required to be subjected to vector decomposition to obtain a scattering capability parameter Dv under the condition of measuring vertical polarization j And measuring the scattering power parameter Dh under the condition of horizontal polarization j ：
Figure BDA0004066998790000057
Figure BDA0004066998790000058
Wherein, is a vector dot product symbol, dv j and Dhj As scalar quantity, the scattering body is measured at theta s And
Figure BDA0004066998790000059
scattering power parameters in the vertical and polarization directions.
And collecting the four types of parameters required to be provided for visualizing the scalar quantities such as Dv and Dh along with the distribution of the scatterer model, and completing the visualization of the scalar quantities such as Dv and Dh.
Triangle cell list { T } of the target scatterer j Co-containing N e Triangle units each T j The scattering power parameters of the vertical polarization and the horizontal polarization of (C) are respectively Dv j and Dhj . Use of TECPLOT360 application to Dv j and Dhj Performing post-processing display, connecting triangle unit points according to a data format given in a TECPLOT360 application program using manual, outputting point coordinates and the value of each triangle unit into a DAT file, and loading the DAT file in the TECPLOT360 to finish Dv j and Dhj (rendering and displaying each triangle cell grid and according to the Dv of each triangle cell) j and Dhj Carrying out gray value discrimination on the triangle by taking the value to obtain a target scatterer post-processing result after gray value discrimination
The j-th unit is at theta s And
Figure BDA00040669987900000510
scattering field ∈in the direction>
Figure BDA00040669987900000511
Expressed as:
Figure BDA00040669987900000512
in the embodiment of the present application, a scatterer model corresponding to a triangle mesh file of a target scatterer is set as shown in fig. 2, and the working frequency is as follows: f=4 GHz, incidence angle: θ=35.0 and,
Figure BDA00040669987900000513
scattering angle: θ s ＝65.6,/>
Figure BDA00040669987900000514
Amplitude of incident electromagnetic wave: e (E) θ ＝1.0,/>
Figure BDA00040669987900000515
After the treatment according to the scheme of the present application, the obtained scalar field distribution of the surface Dv is shown in fig. 3, the obtained scalar field distribution of the surface Dh is shown in fig. 4, and in the embodiment of the present application, the scattering center which can be displayed is shown in the circle position as above, that is, the present application does not need to perform time-consuming integral operation, but only needs to collect the scattering capability parameter Dv under the condition of measuring vertical polarization through the post-treatment technology j And measuring the scattering power parameter Dh under the condition of horizontal polarization j And the file format supported by the post-processing application software TECPLOT360 is formed together with the grid data, so that the gray value of the scattering capability parameter on the target scatterer is judged, the obtained scattering capability parameter distribution on the scatterer characterizes the scattering capability of the target under specific incidence and scattering angles, and the position of the scattering center (such as a fluctuation point in a circle of fig. 3 and 4) can be directly and conveniently seen in the figure.
S3, traversing a triangle unit list { T } j Each triangular cell in the sequence, the fields of the individual triangular cells are obtained in step S2 and the fields of the whole target scatterer are determined.
The step S3 includes:
s301, traversing a triangle unit list { T } j Each triangular unit in the array, obtaining a scattering field of each triangular unit according to the step S2;
s302, calculating the theta of the whole target s And
Figure BDA0004066998790000061
scattering field ∈in the direction>
Figure BDA0004066998790000062
Expressed as:
Figure BDA0004066998790000063
according to the radiation density Dv of each cell under vertical and polarized conditions j and Dhj Is able to acquire a scattered field from (21) and (22)
Figure BDA0004066998790000064
Finally, it should be noted that: the above description is only a preferred embodiment of the present invention, and the present invention is not limited thereto, but it is to be understood that the present invention is described in detail with reference to the above embodiments, and modifications, for example, variations in the names of the methods, etc. may be made to the methods described in the above embodiments by those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (7)

1. A method for characterizing a scattered field of an electromagnetic target, comprising: the method comprises the following steps:
s1, giving a triangular mesh file of a target scatterer, wherein the triangular mesh file comprises a point list { P } of the target scatterer i Sum triangle cell list { T } j }；
S2, applying frequency domain uniform plane electromagnetic waves to the target scatterer, and for { T } j J-th triangle unit T in } j Calculating the scattering power parameter Dv under the condition of vertical polarization j And scattering power parameter Dh under horizontal polarization conditions j And calculates j triangle units T j Is a scattering field of (2);
s3, traversing a triangle unit list { T } j Each triangular cell in the sequence, the fields of the individual triangular cells are obtained in step S2 and the fields of the whole target scatterer are determined.
2. A method of fringe field characterization of an electromagnetic object as recited in claim 1 wherein: point list { P } of the target scatterer i Comprises N in } p The point coordinates, i-th point P i Is marked P by the seat i :(P i .x,P i .y,P i .z)；
Triangle cell list { T } of the target scatterer j Comprises N in } e The three point numbers of the j-th triangle unit defined according to the right-hand rotation direction are denoted as (T) j .i 1 ,T j .i 2 ,T j .i 3 )；
When the coordinates of the kth point in the jth triangle unit are queried, the coordinates are first selected from the triangle unit list { T } j Find the j-th triangle unit data T j Then from T j The sequence number ip=t of the kth point obtained in (a) j .i k Finally from the point list { P i The coordinate of the kth point in the jth triangle unit is found to be (P) ip .x,P ip .y,P ip .z)。
3. A method of fringe field characterization of an electromagnetic object as recited in claim 2, wherein: the step S2 includes:
s201 for { T ] j J-th triangle unit T in } j Acquiring coordinates of a central point of a triangle unit and a normal vector of the triangle unit, calculating the area of the triangle unit and the vector of the area of the triangle unit, and constructing a j-th triangle unit T j Surface local coordinate system C of (2) j ；
S202, calculating the surface current distribution J at the center position of the jth unit PO (T j c )；
S203, calculating the radiation density of a scattered field formed by the jth triangle unit, and calculating a scattering capacity parameter Dv under the condition of vertical polarization j And scattering power parameter Dh under horizontal polarization conditions j 。
4. A method of fringe field characterization of an electromagnetic object as recited in claim 2, wherein: the step S201 includes:
a1, calculating coordinates of a central point of the triangle unit:
from the triangle element list { T j The sequence numbers of the three points in the j-th triangle element are i1 = T, respectively j .i 1 、i2＝T j .i 2 and i3＝Tj .i 3 From the list of points { P i The coordinates of the three points obtained in the process are respectively P i1 、P i2 and Pi3 Then there is a cell center point coordinate T j c Is that;
T j c ＝(P i1 +P i2 +P i3 )/3 (1)
a2, calculating a normal vector of the triangle unit:
obtaining triangle unit list { T } j Coordinates of three points in the j-th triangle unit are respectively: p (P) i1 、P i2 and Pi3 Calculating unit normal vector of jth triangle unit
Figure FDA0004066998670000021
The method comprises the following steps:
Figure FDA0004066998670000022
wherein norm is a vector normalization function: norm (a) =a/len (a), len being a vector length function; x is the vector cross divisor;
a3, calculating the area of the triangle unit:
obtaining triangle unit list { T } j Coordinates of three points in the j-th triangle unit are respectively: p (P) i1 、P i2 and Pi3 Calculating the cell area S of the jth triangle cell j The method comprises the following steps:
S j ＝0.5len[(P i1 -P i2 )×(P i3 -P i2 )] (3)
a4, calculating a triangle area vector:
obtaining triangle unit list { T } j Coordinates of three points in the j-th triangle unit are respectively: p (P) i1 、P i2 and Pi3 Calculate j' thCell area vector S of triangle cell j The method comprises the following steps:
S j ＝0.5[(P i1 -P i2 )×(P i3 -P i2 )] (4)
a5, establishing a surface local coordinate system:
surface local coordinate system C of jth triangle unit j Comprising four elements: center point coordinates
Figure FDA0004066998670000023
Horizontal axis unit vector
Figure FDA0004066998670000024
Vertical axis unit vector->
Figure FDA0004066998670000025
Height axis unit vector +.>
Figure FDA0004066998670000026
The surface part of the triangle unit can be simplified by the coordinate and vector transformation of the four elements of the surface part coordinate system, which C j The four quantities in (a) are calculated as follows:
Figure FDA0004066998670000027
Figure FDA0004066998670000028
Figure FDA0004066998670000029
Figure FDA00040669986700000210
wherein ,
Figure FDA00040669986700000211
is the center point coordinate of the triangle unit, +.>
Figure FDA00040669986700000212
Is from vertex P of triangle unit number i1 i1 Point to vertex number i 2P i2 The unit vector of the edge of (2); />
Figure FDA00040669986700000213
The normal vector is the outside unit normal vector of the plane where the triangle unit is located; />
Figure FDA00040669986700000214
Is->
Figure FDA00040669986700000215
and />
Figure FDA00040669986700000216
Forming right-hand orthogonal coordinate system, the right-hand sequence is +.>
Figure FDA00040669986700000217
5. A method of fringe field characterization of an electromagnetic object as recited in claim 2, wherein: the step S202 includes:
let the electromagnetic field of the frequency domain planar electromagnetic wave be expressed as:
Figure FDA0004066998670000031
Figure FDA0004066998670000032
Figure FDA0004066998670000033
Figure FDA0004066998670000034
Figure FDA0004066998670000035
wherein ,k0 ＝2πf/c 0 Is free space wave number; c 0 Is the electromagnetic wave speed in free space; f is the working frequency; e (E) θ Is the vertically polarized component of the planar electromagnetic wave electric field,
Figure FDA0004066998670000036
for the horizontal polarization components of the planar electromagnetic field, θ and +.>
Figure FDA0004066998670000037
Is the incident angle of plane electromagnetic wave, +.>
Figure FDA0004066998670000038
Is the incident direction unit vector of the plane electromagnetic wave, +.>
Figure FDA0004066998670000039
Vertical polarization unit vector for incident electromagnetic wave, +.>
Figure FDA00040669986700000310
Horizontal polarization unit vector for incident electromagnetic wave, +.>
Figure FDA00040669986700000311
For complex number, η=120pi is free space electromagnetic wave impedance, E is electric field, H is magnetic field, r is field point coordinate vector；/>
Figure FDA00040669986700000312
and />
Figure FDA00040669986700000313
A unit coordinate vector is a global coordinate system;
induced current J at any point r on the area of the target scatterer illuminated by the uniform planar electromagnetic wave PO (r) is:
Figure FDA00040669986700000314
wherein ,
Figure FDA00040669986700000315
is the unit normal vector of the surface of the scattering body at the r position; x is a vector cross symbol;
coordinates of the center point of the jth triangle unit
Figure FDA00040669986700000325
Calculating the surface current distribution at the center of the jth cell as r +.>
Figure FDA00040669986700000326
6. A method of fringe field characterization of an electromagnetic object as recited in claim 2, wherein: the step S203 includes:
the jth triangle element to the diffuser surface is at θ s And
Figure FDA00040669986700000316
analyzing the scattered field in the direction to obtain the far field contribution of the unit to the scattered field, namely the core characterization quantity Ra of electromagnetic scattering phenomenon j Expressed as:
Figure FDA00040669986700000317
Figure FDA00040669986700000318
Figure FDA00040669986700000319
Figure FDA00040669986700000320
wherein ω=2pi f, μ is the free space permeability,
Figure FDA00040669986700000321
is the scattering direction unit vector of the plane electromagnetic wave, +.>
Figure FDA00040669986700000322
and />
Figure FDA00040669986700000323
Is the vertical and horizontal polarization unit vector in the scattering direction of electromagnetic wave, R s Monitoring the radius of the sphere where the scattered field is located; ra (Ra) j Radiation factor for forming a scattering field for the j-th triangle element,>
Figure FDA00040669986700000324
is a green function;
the radiation density of the scattered field formed by the jth triangle unit is Da j ＝Ra j /S j Due to Da j Is an omnidirectional vector, and is required to be subjected to vector decomposition to obtain a scattering capability parameter Dv under the condition of measuring vertical polarization j And measuring scattering power under horizontal polarization conditionNumber Dh j ：
Figure FDA0004066998670000041
Figure FDA0004066998670000042
Wherein, is a vector dot product symbol, dv j and Dhj As scalar quantity, the scattering body is measured at theta s And
Figure FDA0004066998670000043
scattering capability parameters in vertical and polarization in the direction;
the j-th unit is at theta s And
Figure FDA0004066998670000044
scattering field ∈in the direction>
Figure FDA0004066998670000045
Expressed as:
Figure FDA0004066998670000046
7. a method of fringe field characterization of an electromagnetic object as recited in claim 5, wherein: the step S3 includes:
s301, traversing a triangle unit list { T } j Each triangular unit in the array, obtaining a scattering field of each triangular unit according to the step S2;
s302, calculating the theta of the whole target s And
Figure FDA0004066998670000047
scattering field ∈in the direction>
Figure FDA0004066998670000048
Expressed as:
Figure FDA0004066998670000049
according to the radiation density Dv of each cell under vertical and polarized conditions j and Dhj Is able to acquire a scattered field from (21) and (22)
Figure FDA00040669986700000410
/>
CN202310079398.1A 2023-01-13 2023-01-13 Scattered field characterization method of electromagnetic target Active CN116125423B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310079398.1A CN116125423B (en) 2023-01-13 2023-01-13 Scattered field characterization method of electromagnetic target

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310079398.1A CN116125423B (en) 2023-01-13 2023-01-13 Scattered field characterization method of electromagnetic target

Publications (2)

Publication Number Publication Date
CN116125423A true CN116125423A (en) 2023-05-16
CN116125423B CN116125423B (en) 2023-09-01

Family

ID=86302521

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310079398.1A Active CN116125423B (en) 2023-01-13 2023-01-13 Scattered field characterization method of electromagnetic target

Country Status (1)

Country Link
CN (1) CN116125423B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5588032A (en) * 1992-10-14 1996-12-24 Johnson; Steven A. Apparatus and method for imaging with wavefields using inverse scattering techniques
CN104112044A (en) * 2014-07-04 2014-10-22 同济大学 Efficient analyzing method for superfine line structure object electromagnetic property
CN108388732A (en) * 2018-02-27 2018-08-10 中国人民解放军空军工程大学 Plunder extra large Target multipath scattering properties emulated computation method and system
CA2964548A1 (en) * 2017-04-18 2018-10-18 University Of Manitoba Method for determining electric field and related applications
KR20190059384A (en) * 2017-11-23 2019-05-31 인하대학교 산학협력단 Apparatus and method for measuring of radar cross section based on physical theory of diffraction
CN113376597A (en) * 2021-06-18 2021-09-10 西安电子科技大学 Complex terrain electromagnetic scattering rapid simulation method based on digital elevation map and GPU
WO2021260405A1 (en) * 2020-06-26 2021-12-30 Surf Technology As Methods and instrumentation for estimation of wave propagation and scattering parameters
CN114488133A (en) * 2022-03-08 2022-05-13 北京卫星信息工程研究所 Method for extracting and classifying multidimensional scattering characteristics of satellite-borne GNSS-S radar ship
CN115267720A (en) * 2022-08-18 2022-11-01 中国人民解放军91977部队 Method for calculating composite electromagnetic scattering RCS (radar cross section) of marine ship target

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5588032A (en) * 1992-10-14 1996-12-24 Johnson; Steven A. Apparatus and method for imaging with wavefields using inverse scattering techniques
CN104112044A (en) * 2014-07-04 2014-10-22 同济大学 Efficient analyzing method for superfine line structure object electromagnetic property
CA2964548A1 (en) * 2017-04-18 2018-10-18 University Of Manitoba Method for determining electric field and related applications
KR20190059384A (en) * 2017-11-23 2019-05-31 인하대학교 산학협력단 Apparatus and method for measuring of radar cross section based on physical theory of diffraction
CN108388732A (en) * 2018-02-27 2018-08-10 中国人民解放军空军工程大学 Plunder extra large Target multipath scattering properties emulated computation method and system
WO2021260405A1 (en) * 2020-06-26 2021-12-30 Surf Technology As Methods and instrumentation for estimation of wave propagation and scattering parameters
CN113376597A (en) * 2021-06-18 2021-09-10 西安电子科技大学 Complex terrain electromagnetic scattering rapid simulation method based on digital elevation map and GPU
CN114488133A (en) * 2022-03-08 2022-05-13 北京卫星信息工程研究所 Method for extracting and classifying multidimensional scattering characteristics of satellite-borne GNSS-S radar ship
CN115267720A (en) * 2022-08-18 2022-11-01 中国人民解放军91977部队 Method for calculating composite electromagnetic scattering RCS (radar cross section) of marine ship target

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
宋东安;邢芳;温定娥;张崎;: "舰船RCS整形设计中的电磁分析", 中国舰船研究, no. 03, pages 52 - 55 *
张兆;沈孟育;张涵信;: "基于非结构Cartesian网格的电磁散射场计算", 清华大学学报(自然科学版), no. 11, pages 2068 - 2071 *
李晓峰;谢拥军;樊君;王元源;: "考虑棱边散射的半空间复杂导体目标高频分析方法", 物理学报, no. 02, pages 908 - 913 *

Also Published As

Publication number Publication date
CN116125423B (en) 2023-09-01

Similar Documents

Publication Publication Date Title
CN102508220B (en) Method for obtaining radar cross section (RCS) of homogeneous bi-isotropic medium object
CN109100692B (en) Coarse surface and multiple target composite scattering simulation method based on iterative physical optics
US5812434A (en) Electromagnetic field strength calculator having function of displaying currents to be analyzed
CN102385655B (en) Simulation method for radiation coupling electromagnetic susceptibility of electronic equipment
CN112446152B (en) Antenna far-field directional pattern analysis method based on infinitesimal small dipole model deformation array
CN112733364A (en) Foil strip cloud scattering rapid calculation method based on impedance matrix blocking
CN115546015A (en) Method for generating characteristic image of harmonic induction field when grid curved surface is subjected to plane electromagnetic wave
He et al. Closed-form solutions for analysis of cylindrically conformal microstrip antennas with arbitrary radii
CN113221370B (en) FSS radome modeling method based on conical surface projection
CN116125423B (en) Scattered field characterization method of electromagnetic target
CN105302992A (en) Quick optimization and simulation design method for stirrer in reverberation chamber
CN106777536A (en) Electro-magnetic far-field two, three-dimensional visual processing method based on fine Electromagnetic Simulation
CN109783829B (en) Electromagnetic field prediction method for three-dimensional FEM (field emission modeling) and two-dimensional FMM (frequency modulation modeling)
CN114896868A (en) Estimation method for radiation field of ultra-wide spectrum electromagnetic pulse radiation Vivaldi antenna array
CN112285435B (en) Equivalent simulation method of high-power magnetic field radiation source
CN115983053A (en) Electromagnetic simulation method of antenna with dielectric substrate based on moment method
CN115659607A (en) Method for determining multilayer wave-transparent structure
CN111931353A (en) Scattered field solving method applied to simulated FSS structure
Wang et al. Edge diffraction in NURBS-UTD method
CN111735996A (en) Multipath interference suppression method and device for mathematical wave-absorbing darkroom construction
CN116484724B (en) Antenna detection method based on electromagnetic modeling parameters of high-order moment method
Hassan et al. Time-Domain Sensitivity Analysis for Conductivity Distribution in Maxwell's Equations
CN114076924B (en) Method for realizing large-scale complex target classification recognition based on multi-frequency echo data
CN116203505B (en) Orthogonal matching pursuit sound source identification method and device based on block sparse Bayes
CN111832157B (en) Large-scale quasi-periodic structure electromagnetic scattering characteristic analysis method based on sub-global basis function method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant