CN108445303A - A kind of Electromagnetic scatter simulated behavior method - Google Patents
A kind of Electromagnetic scatter simulated behavior method Download PDFInfo
- Publication number
- CN108445303A CN108445303A CN201810188886.5A CN201810188886A CN108445303A CN 108445303 A CN108445303 A CN 108445303A CN 201810188886 A CN201810188886 A CN 201810188886A CN 108445303 A CN108445303 A CN 108445303A
- Authority
- CN
- China
- Prior art keywords
- face element
- incident
- target
- face
- near 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention relates to a kind of Electromagnetic scatter simulated behavior methods, are related to electromagnetic scattering technical field.Wherein, this method includes:Subdivision is carried out to target, to obtain multiple face elements;Matrix equation of each face element under Near Field is generated according to Multilevel fast multipole algorithm, then, the electric current on face element is obtained by the matrix equation;The corresponding polarization reception electric field of face element is determined according to the electric current on the face element;Vector superposed processing is carried out to the corresponding polarization reception electric field of all face elements, to obtain the characterization parameter of Electromagnetic Scattering Characteristics of the target under Near Field.By above step, the precision and universality of Target near field Electromagnetic Scattering Characteristics analog result can be improved, is widely portable to Electromagnetic scatter characteristic research of the target under all kinds of antenna irradiations and missile-target encounter condition.
Description
Technical field
The present invention relates to electromagnetic scattering technical field more particularly to a kind of Electromagnetic scatter simulated behavior methods.
Background technology
The Electromagnetic scatter characteristic research of target is a key areas of target scattering characteristics research.With maturation
Far field electromagnetic scattering theory is different, and extended target effect is presented in target under Near Field, and radar return characteristic is by spherical wave
The multifactor influence such as effect, antenna radiation pattern, local irradiation.
Currently, mostly use greatly high-frequency approximation method (such as physical optical method) to the Electromagnetic scatter characteristic of target into
Row calculates.Existing method can not calculate the effect of electromagnetic coupling generation between each position of complex target, and not consider antenna side mostly
Influence to the factors such as figure and local irradiation.
With the increase of target local environment complexity and to the requirements for high precision of result of calculation, for Target near field
Electromagnetic Scattering Characteristics Research Requirements are further urgent.Therefore, there is an urgent need for study a kind of near field electromagnetic for accurately simulating complex target
The method of scattering properties, to improve the precision and universality of the calculating of Electromagnetic scatter characteristic.
Invention content
The technical problem to be solved in the present invention is, for the defects in the prior art, it is scattered to provide a kind of near field electromagnetic
Penetrate the analogy method of characteristic.
In order to solve the above technical problem, the present invention provides a kind of analogy methods of Electromagnetic scatter characteristic.
The analogy method of Electromagnetic scatter characteristic provided by the invention includes:Subdivision is carried out to target, it is multiple to obtain
Face element;Matrix equation of each face element under Near Field is generated according to Multilevel fast multipole algorithm, then, by described
Matrix equation obtains the electric current on face element;The corresponding polarization reception electric field of face element is determined according to the electric current on the face element;To institute
There is the corresponding polarization reception electric field of face element to carry out vector superposed processing, it is special to obtain electromagnetic scattering of the target under Near Field
The characterization parameter of property.
Optionally, described that matrix equation of each face element under Near Field is generated according to Multilevel fast multipole algorithm
The step of include:Each corresponding near field incident excitation item of face element and each face element are calculated according to Multilevel fast multipole algorithm
Corresponding impedance matrix, to generate matrix equation of each face element under Near Field.
Optionally, described that the corresponding near field incident excitation item of each face element is calculated according to Multilevel fast multipole algorithm
Step includes:The incident electric fields and incident magnetic on face element surface are obtained according to the polarization incident electric fields model of face element;According to described
The corresponding near field incident excitation item of incident electric fields and incident magnetic structure face element on face element surface.
Optionally, the polarization incident electric fields model includes:Linear polarization incident electric fields model, elliptic polarization incident electric fields
Model or circular polarisation incident electric fields model.
Optionally, the polarization incident electric fields model according to face element obtains the incident electric fields on face element surface and incident magnetic
The step of include:Following affecting parameters are calculated, including:The relative phase of the incident electromagnetic wave on face element surface increases, face element with
The power gain of the distance between transmitting antenna and the corresponding incident electromagnetic wave of face element;By the calculating knot of the affecting parameters
Fruit substitutes into the polarization incident electric fields model of face element, to obtain the incident electric fields on face element surface;According to the right hand of electric field and magnetic field
Coiled relationship obtains the incident magnetic of face element.
Optionally, the method further includes according to the power gain for calculating the corresponding incident electromagnetic wave of face element such as under type:
In the case of the directional diagram of known transmitting antenna, the corresponding incident electromagnetic wave of face element is calculated according to the angle information of incidence wave
Power gain;In the directional diagram but known transmitting antenna power gain on multiple discrete angulars of unknown transmitting antenna
In the case of, the power gain of the corresponding incident electromagnetic wave of face element is obtained according to Lagrangian one-dimensional interpolation.
Optionally, the method further includes according to the power gain for calculating the corresponding incident electromagnetic wave of face element such as under type:
In the case of the beam angle of known transmitting antenna, antenna is obtained the two of incidence wave according to sinc function graft features
Power direction function on a assigned direction, then according to the power direction function calculate face element it is corresponding enter radio
The power gain of magnetic wave.
Optionally, the step of electric current obtained by the matrix equation on face element includes:According to broad sense Minimum Residual
Poor method iteratively solves the matrix equation, to obtain the electric current on face element.
Optionally, the characterization parameter of the Electromagnetic Scattering Characteristics includes:The near field reflections power ratio of target, target it is close
Field RCS.
Optionally, the near field reflections power ratio of the target meets:
In formula, DiIndicate the direction coefficient of transmitting antenna greatest irradiation direction, DsIndicate reception antenna greatest irradiation direction
Direction coefficient, λ be incident electromagnetic wave wavelength, m indicate target m-th of face element, FimIndicate that m-th of face element of target exists
Directivity function on the incident electromagnetic wave direction of transmitting antenna transmitting, FsmIndicate that m-th of face element of target is received in reception antenna
Scattering electromagnetic wave direction on directivity function, RimIndicate the distance between m-th of face element and the transmitting antenna of target, Rsm
Indicate the distance between m-th of face element and the reception antenna of target, EimIndicate the polarization incident electric fields of the m face element, EsmTable
Show the polarization reception electric field of m-th of face element.
Optionally, the near field RCS of the target meets:
In formula, RiIndicate the distance between site reference point and the transmitting antenna of target, RsIndicate the site reference of target
The distance between point and reception antenna, Fi 2Indicate the normalized power directivity function reference value of transmitting antenna main lobe, Fs 2It indicates
The normalized power directivity function reference value of reception antenna main lobe, m indicate m-th of face element of target, FimIndicate the m of target
Directivity function of a face element on the incident electromagnetic wave direction that transmitting antenna emits, FsmIndicate that m-th of face element of target is connecing
Receive the directivity function on the scattering electromagnetic wave direction that antenna receives, RimBetween m-th of the face element and transmitting antenna that indicate target
Distance, RsmIndicate the distance between m-th of face element and the reception antenna of target, EimIndicate that the polarization of m-th of face element enters radio
, EsmIndicate the polarization reception electric field of m-th of face element.
Implement the present invention, has the advantages that:
The present invention generates each face element close by carrying out face element subdivision to target, according to Multilevel fast multipole algorithm
Then matrix equation under field condition obtains the electric current on face element, further according on the face element by the matrix equation
Electric current determines the corresponding polarization reception electric field of face element, and vector superposed processing is carried out to the corresponding polarization reception electric field of all face elements
And etc., the characterization parameter of Electromagnetic Scattering Characteristics of the target under Near Field can be obtained.Compared with prior art, this hair
The bright precision and universality that can improve Target near field Electromagnetic Scattering Characteristics analog result, is widely portable to target all kinds of
Near-field scattering characteristics research under antenna irradiation and missile-target encounter condition.
Description of the drawings
Fig. 1 is the key step schematic diagram of the Electromagnetic scatter simulated behavior method of the embodiment of the present invention one;
Fig. 2 is the key step schematic diagram of the Electromagnetic scatter simulated behavior method of the embodiment of the present invention two.
Fig. 3 is target-based coordinate system and antenna coordinate system schematic diagram;
Fig. 4 is incident schematic diagram of the electromagnetic wave that sends out of transmitting antenna at target difference face element;
Fig. 5 is the antenna beam schematic diagram in antenna coordinate system.
Specific implementation mode
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, ordinary skill
The every other embodiment that personnel are obtained without making creative work belongs to the model that the present invention protects
It encloses.
Embodiment one
Fig. 1 is the key step schematic diagram of the Electromagnetic scatter simulated behavior method of the embodiment of the present invention one.Such as Fig. 1
Shown, Electromagnetic scatter simulated behavior method provided in an embodiment of the present invention includes:
Step S101, subdivision is carried out to target, to obtain multiple face elements.
Illustratively, target can be split into multiple Triangular patch.
Step S102, matrix equation of each face element under Near Field is generated according to Multilevel fast multipole algorithm.
Multilevel fast multipole method is derived based on Maxwell equations and stringent electromagnetic theory, can cope with various structures
And scattering mechanism, computational accuracy is high, is a kind of i.e. efficient and accurate numerical computation method.
Specifically, when generating the matrix equation under Near Field according to multilevel fast multipole method, it is contemplated that spherical surface
The influence of the factors such as wave effect, antenna radiation pattern, local irradiation and electromagnetic coupling effect is avoided and is deposited in high-frequency approximation method
Due to that can not calculate detailed structure, electromagnetic coupling effect can not be calculated, be difficult to handle complex electromagnetic mechanism (such as row
Wave, peristaltic wave etc.) the problems such as the phenomenon that caused computational accuracy is low, calculated curve and measured value are misfitted, and then significantly
Improve the precision of complex target near-field scattering characteristics calculating.
Step S103, the electric current on face element is obtained by the matrix equation.
Illustratively, the matrix equation is iteratively solved according to broad sense minimum residual method (GMRES), to obtain on face element
Electric current.
Step S104, the corresponding polarization reception electric field of face element is determined according to the electric current on the face element.
Wherein, the corresponding polarization reception electric field of the face element can be regarded as:What reception antenna received comes from target face element
Scattering electric field.
Step S105, vector superposed processing is carried out to the corresponding polarization reception electric field of all face elements, to obtain target close
The characterization parameter of Electromagnetic Scattering Characteristics under field condition.
Wherein, the characterization parameter of the Electromagnetic Scattering Characteristics includes:The near field reflections power ratio of target, the near field of target
RCS (radar cross section).
In embodiments of the present invention, the accurate solution of the Electromagnetic scatter characteristic of target is established by above step
Model, the solving model can be used for calculating near-field scattering characteristics of the target under all kinds of antenna irradiations and missile-target encounter condition.
Embodiment two
Fig. 2 is the key step schematic diagram of the Electromagnetic scatter simulated behavior method of the embodiment of the present invention two.Such as Fig. 2
Shown, the Electromagnetic scatter simulated behavior method of the embodiment of the present invention includes:
Step S201, target-based coordinate system and antenna coordinate system are established.
Fig. 3 is target-based coordinate system and antenna coordinate system schematic diagram.As shown in figure 3, target-based coordinate system can use xtytztIt indicates,
Antenna coordinate system can use xayazaIt indicates.Wherein, coordinate of the antenna in target-based coordinate system is (x0,y0,z0).Meanwhile it can be used
Three Eulerian angles (α, β, γ) describe the relative rotational orientation of two coordinate systems.Wherein α is angle of precession, and β is nutational angle, γ
For angle of rotation, the spin matrix on each gyrobearing is represented by:
In turn, it can be obtained by the spin matrix of target-based coordinate system to antenna coordinate system and be:
R (α, β, γ)=Z (γ) N (β) Z (α)
It in turn, can be by position translation and the rotation transformation based on three Eulerian angles makes target-based coordinate system and antenna coordinate
System overlaps, and describes target and the relative position and posture of antenna with this, while can be used for each variable (such as wave vector vector) two
It is mutually converted in a coordinate system.
Step S202, subdivision is carried out to target, to obtain multiple face elements.
In this step, target can be split into multiple Triangular patch.Assuming that m-th of Triangular patch of target is in mesh
It is (x to mark the centre coordinate in coordinate systemm,ym,zm), m=1,2, Λ N0.Wherein, N0For the total of the Triangular object model included by target
Number.Further, it is assumed that coordinate of the transmitting antenna in target-based coordinate system is (x0,y0,z0).In turn, it is convenient in the next step
Calculate every affecting parameters corresponding to m-th of face element of target.
Step S203, the corresponding affecting parameters of each face element are calculated, including:The distance between face element and transmitting antenna, face
The relative phase increase of the incident electromagnetic wave on first surface, the power gain of the corresponding incident electromagnetic wave of face element.
As shown in figure 4, when the distance between antenna and target are unsatisfactory for far field condition, spherical wave is affected, main
Showing as the relative attitude between parameters and distance, target and antenna such as wave vector vector, relative phase increase of target has
It closes.In addition, due to target different parts between antenna at a distance from different (being unable to approximately equal), so, target difference face element with
The range attenuation of incidence and scattering field between antenna is also different.Therefore, it is necessary to be based on face element to calculate the affecting parameters,
It specifically includes:
A) distance between m-th of face element and transmitting antenna is calculated.
Illustratively, distance between m-th of face element and transmitting antenna can be calculated according to following formula:
In formula,The distance vector between m-th of face element and transmitting antenna, RinFor m-th face element and transmitting antenna it
Between distance value.
B) the unit wave vector vector k of the incident electromagnetic wave at m-th of face element center is calculatedin。
Illustratively, the unit wave vector vector k of the incident electromagnetic wave at m-th of face element center can be calculated according to following formulain
For:
The relative phase of the incident electromagnetic wave at each face element center is increased, then can be phase zero by aerial position
Point is calculated.In turn, the relative phase of the incident electromagnetic wave at m-th of face element center increases Δ φimFor:
Wherein,For wave number.
C) angle of the incident wave vector of m-th of face element in antenna coordinate system is calculated.
As shown in figure 5, the unit wave vector vector at m-th of face element center can be expressed as k in antenna coordinate systemim-a, m-th
Angle of the incident wave vector of face element in antenna coordinate system can be expressed as θim、
Illustratively, angle of the incident wave vector of m-th of face element in antenna coordinate system can be calculated according to following formula:
θim=arcos (z (kim-a))
kim-a=R (α, β, γ) kim
Wherein, θimFor the incident wave vector and z of m-th of face elementaThe angle of axis,Incident wave vector for m-th of face element exists
xaoyaThe projection of plane and xaThe angle of axis.
D) power gain of the incident electromagnetic wave on m-th of face element surface is calculated.
Illustratively, the power that can be divided into the incident electromagnetic wave that following two situations calculate separately m-th of face element surface increases
Benefit.
The first situation:The directional diagram of known transmitting antenna or the power gain on discrete angular.
1.1) in the directional diagram of known transmitting antenna (such as the power direction function of antennaOr antenna
Field strength directivity function) in the case of, the corresponding incoming electromagnetic of face element can be calculated according to the angle information of incidence wave
The power gain of wave, calculation formula are as follows:
In formula, GimFor the power gain of the incident electromagnetic wave on m-th of face element surface.θim、For the incidence of the m face elements
Angle (referred to as " incident angle of the m face element ") of the wave vector in antenna coordinate system.
1.2) in the directional diagram of unknown transmitting antenna but known transmitting antenna is in multiple discrete angular θj(j=1,2 ...
M)、Power gain in the case of, the corresponding incidence of face element can be obtained according to Lagrangian one-dimensional interpolation
The power gain of electromagnetic wave, calculation formula are as follows:
In formula, Pexample1(θj) it is in discrete angular θjOn power gain,For in discrete angularOn
Power gain, P (θim)、For the power gain on the incident angle of m-th of face element,It is m-th
The power gain of the corresponding incident electromagnetic wave of face element,For m-th of face element it is corresponding enter radio
The field strength gain of magnetic wave.
The second situation:The case where only limiting the beam angle of transmitting antenna.
In goal in research pole when Electromagnetic scatter characteristic, since the local irradiation effect under the Near Field of pole is especially bright
Aobvious, electromagnetic scattering result is extremely sensitive for antenna lobe width, therefore antenna radiation pattern is only needed to meet the lobe width set
It is required that.
Specifically, in the latter case, antenna can be obtained according to sinc function graft features at two of incidence wave
Then power direction function on assigned direction calculates the corresponding incoming electromagnetic of face element according to the power direction function
The power gain of wave.Wherein, the power direction function on described two assigned directions is as follows:
According to the power direction function on described two assigned directions, the incident electromagnetic wave of m-th of face element can be calculated
Power gainAnd the field strength gain of the incident electromagnetic wave of m-th of face element
Step S204, polarization incident electric fields model that the face element corresponding affecting parameters are substituted into face element, to obtain face
The incident electric fields on first surface.
Wherein, the polarization incident electric fields model of the face element includes:The linear polarization incident electric fields model of face element, face element
The circular polarisation incident electric fields model of elliptic polarization incident electric fields model or face element.Wherein, the linear polarization incident electric fields model of face element
The horizontal polarization model, vertical polarization model or general linear polarization model of face element can be divided into again.
Illustratively, horizontal polarization incident electric fields model of m-th of face element in antenna coordinate system is represented by:
EH-a(z)=0
EimH-a=EH-acos(ωt+φH)
EimH=R-1(α,β,γ)gEimH-a
In formula,For angle of the incident wave vector in antenna coordinate system of m-th of face element;ω is the angular frequency of electromagnetic wave;
T is the time;φHFor arbitrary constant;R-1(α, beta, gamma) is taken by R (α, beta, gamma) against acquiring, for antenna coordinate system to target-based coordinate system
Spin matrix;EH-a(x)、EH-a(y)、 EH-a(z) horizontal polarization unit of m-th of face element in antenna coordinate system is indicated respectively
Range weight of the incident electric fields in each reference axis;EimH-aThe horizontal polarization unit for being m-th of face element in antenna coordinate system enters
Radio field vector;EimHThe horizontal polarization per incident electric field intensity for being m-th of face element in target-based coordinate system.
In addition, vertical polarization incident electric fields model of m-th of face element in antenna coordinate system is represented by:
EV-a(z)=- sin θim
EimV-a=EV-acos(ωt+φV)
EimV=R-1(α,β,γ)gEimV-a
In formula, θim、For angle of the incident wave vector in antenna coordinate system of m-th of face element;ω is the angle of electromagnetic wave
Frequency;T is the time;φVFor arbitrary constant;R-1(α, beta, gamma) is taken by R (α, beta, gamma) against acquiring, for antenna coordinate system to target
The spin matrix of coordinate system;EV-a(x)、EV-a(y)、 EV-a(z) vertical pole of m-th of face element in antenna coordinate system is indicated respectively
Range weight of the change per incident electric field in each reference axis;EimV-aFor vertical polarization of m-th of face element in antenna coordinate system
Per incident electric field intensity;EimVThe vertical polarization per incident electric field intensity for being m-th of face element in target-based coordinate system.
In addition, general linear polarization incident electric fields model of m-th of face element in antenna coordinate system is represented by:
EimL=R-1(α,β,γ)gEimL-a
In formula, C0For arbitrary constant;ω is the angular frequency of electromagnetic wave;T is the time;φ is arbitrary constant;R-1(α,β,γ)
It is taken by R (α, beta, gamma) against acquiring, is spin matrix of the antenna coordinate system to target-based coordinate system;EimL-aIt is m-th of face element in antenna
General polarization per incident electric field intensity in coordinate system;EimLThe general polarization for being m-th of face element in target-based coordinate system is single
Position incident electric field vector.
It, can be by giving certain about elliptic polarization of m-th of face element in antenna coordinate system, circular polarisation incident electric fields model
Vertical polarization incident electric field vector and the linear combination of horizontal polarization incident electric field vector after one initial phase obtain, herein
No longer provide one by one.
Step S205, according to the right-handed helix relationship of electric field and magnetic field, the incident magnetic on face element surface is obtained.
Step S206, the corresponding near field incident excitation item of face element and face element pair are built according to Multilevel fast multipole algorithm
The impedance matrix answered, to generate matrix equation of the face element under Near Field.
Wherein, the corresponding near field incident excitation item of m-th of face element is represented by Cm, specific as follows:
In formula, gmFor the basic function of selection,For the incident electric fields on m-th of face element surface,
For the incident magnetic on m-th of face element surface, a1For combined field coefficient, η is wave impedance, and dS' is the face element differential of m-th of face element,
∫S'DS' indicates that integral and calculating is carried out in m-th of face element.
Wherein, corresponding impedance matrix is represented by [Z between m-th of face element and n-th of face elementmn], which reflects mesh
Mark the electromagnetic coupling attribute between itself different seamed edge electric current.
In turn, matrix equation of m-th of face element under Near Field is:
In formula, bnFor face element seamed edge current coefficient, N1To there is the sum of other face elements of coupling effect to m-th of face element.
Step S207, the matrix equation is iteratively solved according to broad sense minimum residual method, to obtain the electric current on face element.
Specifically, the electric current J on m-th of face elementm(r) it is represented by:
In formula, niIndicate the common edge corresponding to i-th side of m-th of face element, gniFor the basic function of selection, bniFor phase
The expansion coefficient for the basic function answered.
Step S208, the corresponding polarization reception electric field of face element is determined according to the electric current on the face element.
Specifically, Stratton-Chu integral formulas can be used and obtain the scattering electric field E' of m-th of face elementsm(r).Then,
Consider different antennae directional diagram and polarization reception, it can be according to E'sm(r) the polarization reception electric field E on the m face element is obtainedsm
(r)。
Step S209, vector superposed processing is carried out to the corresponding polarization reception electric field of all face elements, to obtain target close
The characterization parameter of Electromagnetic Scattering Characteristics under field condition.
Wherein, the characterization parameter of the Electromagnetic Scattering Characteristics may include:The near field reflections power ratio of target, target it is close
Field RCS (radar cross section).Further, the near field reflections power ratio of the target meets:
In formula, DiIndicate the direction coefficient of transmitting antenna greatest irradiation direction, DsIndicate reception antenna greatest irradiation direction
Direction coefficient, λ be incident electromagnetic wave wavelength, m indicate target m-th of face element, FimIndicate that m-th of face element of target exists
Directivity function on the incident electromagnetic wave direction of transmitting antenna transmitting, FsmIndicate that m-th of face element of target is received in reception antenna
Scattering electromagnetic wave direction on directivity function, RimIndicate the distance between m-th of face element and the transmitting antenna of target, Rsm
Indicate the distance between m-th of face element and the reception antenna of target, EimIndicate the polarization incident electric fields of the m face element, EsmTable
Show the polarization reception electric field of m-th of face element.
Further, the near field RCS of the target meets:
In formula, RiIndicate the distance between site reference point and the transmitting antenna of target, RsIndicate the site reference of target
The distance between point and reception antenna, Fi 2Indicate the normalized power directivity function reference value of transmitting antenna main lobe, Fs 2It indicates
The normalized power directivity function reference value of reception antenna main lobe, m indicate m-th of face element of target, FimIndicate the m of target
Directivity function of a face element on the incident electromagnetic wave direction that transmitting antenna emits, FsmIndicate that m-th of face element of target is connecing
Receive the directivity function on the scattering electromagnetic wave direction that antenna receives, RimBetween m-th of the face element and transmitting antenna that indicate target
Distance, RsmIndicate the distance between m-th of face element and the reception antenna of target, EimIndicate that the polarization of m-th of face element enters radio
, EsmIndicate the polarization reception electric field of m-th of face element.
The method of the embodiment of the present invention can at least realize following one or more technique effects:
1, establish complex target near-field scattering characteristics simulation model, can be used for calculating target in all kinds of antenna irradiations and
Near-field scattering characteristics under missile-target encounter condition.When it is implemented, being simulated according to the Target near field scattering properties simulation model
Certain guided missile model is compared in the near-field scattering characteristics under encounter conditions, and with actual measured results, comparing result
Fully demonstrate the accuracy of the complex target near-field scattering characteristics simulation model.
2, by using multilevel fast multipole algorithm, while considering spherical surface wave effect, antenna radiation pattern, local irradiation
And the influence factors such as electromagnetic coupling effect, it avoids present in high-frequency approximation method since detailed structure can not be calculated, can not
The caused calculating of the problems such as calculating electromagnetic coupling effect, being difficult to handle complex electromagnetic mechanism (such as traveling wave, peristaltic wave)
The phenomenon that precision is low, and calculated curve is misfitted with measured value greatly improves the precision of complex target near-field scattering characteristics calculating.
3. obtained near field RCS and near field reflections power ratio, it can be used for precise quantification and describe dissipating for target under Near Field
Penetrate characteristic.
Finally it should be noted that:The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations;To the greatest extent
Present invention has been described in detail with reference to the aforementioned embodiments for pipe, it will be understood by those of ordinary skill in the art that:It is still
Can be with technical scheme described in the above embodiments is modified, or which part technical characteristic is equally replaced
It changes;And these modifications or replacements, various embodiments of the present invention technical solution that it does not separate the essence of the corresponding technical solution
Spirit and scope.
Claims (10)
1. a kind of Electromagnetic scatter simulated behavior method, which is characterized in that the method includes:
Subdivision is carried out to target, to obtain multiple face elements;
Matrix equation of each face element under Near Field is generated according to Multilevel fast multipole algorithm, then, passes through the square
Battle array equation obtains the electric current on face element;
The corresponding polarization reception electric field of face element is determined according to the electric current on the face element;
Vector superposed processing is carried out to the corresponding polarization reception electric field of all face elements, to obtain electromagnetism of the target under Near Field
The characterization parameter of scattering properties.
2. according to the method described in claim 1, it is characterized in that, described generate each face according to Multilevel fast multipole algorithm
The step of matrix equation of the member under Near Field includes:
Each corresponding near field incident excitation item of face element and the corresponding resistance of each face element are calculated according to Multilevel fast multipole algorithm
Anti- matrix, to generate matrix equation of each face element under Near Field.
3. according to the method described in claim 2, it is characterized in that, described calculate each face according to Multilevel fast multipole algorithm
The step of first corresponding near field incident excitation item includes:
The incident electric fields and incident magnetic on face element surface are obtained according to the polarization incident electric fields model of face element;According to the face element table
The corresponding near field incident excitation item of incident electric fields and incident magnetic structure face element in face.
4. according to the method described in claim 3, it is characterized in that, the polarization incident electric fields model includes:Linear polarization is incident
Electric field model, elliptic polarization incident electric fields model or circular polarisation incident electric fields model.
5. according to the method described in claim 3, it is characterized in that, described obtain face according to the polarization incident electric fields model of face element
The step of incident electric fields and incident magnetic on first surface includes:
Following affecting parameters are calculated, including:The relative phase of the incident electromagnetic wave on face element surface increases, face element and transmitting antenna it
Between the corresponding incident electromagnetic wave of distance and face element power gain;The result of calculation of the affecting parameters is substituted into face element
Polarize incident electric fields model, to obtain the incident electric fields on face element surface;According to the right-handed helix relationship of electric field and magnetic field, face is obtained
The incident magnetic of member.
6. according to the method described in claim 5, it is characterized in that, the method further includes calculating face element pair according to such as under type
The power gain for the incident electromagnetic wave answered:
In the case of the directional diagram of known transmitting antenna, the corresponding incoming electromagnetic of face element is calculated according to the angle information of incidence wave
The power gain of wave;In the power gain on multiple discrete angulars of directional diagram but known transmitting antenna of unknown transmitting antenna
In the case of, the power gain of the corresponding incident electromagnetic wave of face element is obtained according to Lagrangian one-dimensional interpolation.
7. according to the method described in claim 5, it is characterized in that, the method further includes calculating face element pair according to such as under type
The power gain for the incident electromagnetic wave answered:
In the case of the beam angle of known transmitting antenna, antenna is obtained in incidence wave according to sinc function graft features
Power direction function on two assigned directions, then according to the power direction function calculate face element it is corresponding enter radio
The power gain of magnetic wave.
8. according to the method described in claim 1, it is characterized in that, the electric current obtained by the matrix equation on face element
The step of include:
The matrix equation is iteratively solved according to broad sense minimum residual method, to obtain the electric current on face element.
9. according to the method described in claim 1, it is characterized in that, the characterization parameter of the Electromagnetic Scattering Characteristics includes:Target
Near field reflections power ratio, target near field RCS.
10. according to the method described in claim 9, it is characterized in that, the near field reflections power ratio of the target meets:
In formula, DiIndicate the direction coefficient of transmitting antenna greatest irradiation direction, DsIndicate the direction of reception antenna greatest irradiation direction
Coefficient, λ are the wavelength of incident electromagnetic wave, and m indicates m-th of face element of target, FimIndicate m-th of face element of target in transmitting antenna
Directivity function on the incident electromagnetic wave direction of transmitting, FsmIndicate the scattering electricity that m-th of face element of target is received in reception antenna
Directivity function on magnetic wave direction, RimIndicate the distance between m-th of face element and the transmitting antenna of target, RsmIndicate target
The distance between m-th of face element and reception antenna, EimIndicate the polarization incident electric fields of m-th of face element, EsmIndicate m-th of face element
Polarization reception electric field.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810188886.5A CN108445303B (en) | 2018-03-08 | 2018-03-08 | Near-field electromagnetic scattering characteristic simulation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810188886.5A CN108445303B (en) | 2018-03-08 | 2018-03-08 | Near-field electromagnetic scattering characteristic simulation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108445303A true CN108445303A (en) | 2018-08-24 |
CN108445303B CN108445303B (en) | 2020-06-26 |
Family
ID=63193692
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810188886.5A Active CN108445303B (en) | 2018-03-08 | 2018-03-08 | Near-field electromagnetic scattering characteristic simulation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108445303B (en) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110826279A (en) * | 2019-11-07 | 2020-02-21 | 北京环境特性研究所 | Method and device for calculating optical characteristics of solid target |
CN111273281A (en) * | 2020-02-07 | 2020-06-12 | 北京环境特性研究所 | Short-distance electromagnetic measurement method and system |
CN111339650A (en) * | 2020-02-21 | 2020-06-26 | 北京环境特性研究所 | Electromagnetic basis function construction method and device with reduced orthogonal dimension |
CN112069713A (en) * | 2020-09-10 | 2020-12-11 | 上海无线电设备研究所 | Near-field scattering characteristic modeling method, electronic device and storage medium |
CN113092878A (en) * | 2021-03-31 | 2021-07-09 | 北京环境特性研究所 | Test method and detection device for electromagnetic radiation of W-band environment |
CN113408095A (en) * | 2020-03-16 | 2021-09-17 | 启碁科技股份有限公司 | Electromagnetic characteristic analysis method and electronic device |
CN113567941A (en) * | 2021-08-04 | 2021-10-29 | 上海无线电设备研究所 | Target near-field Doppler echo rapid simulation method |
CN113671269A (en) * | 2021-07-05 | 2021-11-19 | 南京航空航天大学 | Spherical near-far field method based on scattering matrix principle |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5225668A (en) * | 1991-06-06 | 1993-07-06 | The United States Of America As Represented By The Secretary Of The Navy | Photonic electromagnetic field sensor apparatus |
CN102054094A (en) * | 2010-11-24 | 2011-05-11 | 南京理工大学 | Fast directional multilevel simulation method for planar microstrip circuit |
CN103198227A (en) * | 2013-04-18 | 2013-07-10 | 南京理工大学 | Electromagnetic scattering analyzing method for superspeed flight targets |
CN103226644A (en) * | 2013-05-17 | 2013-07-31 | 南京理工大学 | Electromagnetic scattering property simulating method based on cylindrical surface equivalent source domain decomposition |
CN103675781A (en) * | 2012-09-25 | 2014-03-26 | 中国航天科工集团第二研究院二〇七所 | A method for accurately acquiring near field echo |
CN104199026A (en) * | 2014-08-28 | 2014-12-10 | 中国科学院电子学研究所 | Backscattering cross section measuring method based on trace scanning two-dimensional near field imaging |
CN104730503A (en) * | 2015-03-18 | 2015-06-24 | 中国科学院电子学研究所 | Method for determining influence on scaling by high-resolution SAR reference target RCS and compensation method |
CN105630740A (en) * | 2014-11-07 | 2016-06-01 | 南京理工大学 | Electromagnetic analysis method on the basis of matrix Taylor series expansion |
CN106093875A (en) * | 2016-06-23 | 2016-11-09 | 中国科学院电子学研究所 | Near field WB-RCS data acquisition and the compression method of denoising is followed the trail of based on base |
CN106156394A (en) * | 2015-04-22 | 2016-11-23 | 南京理工大学 | Electromagnetic property extracting method based on explicit difference scheme |
CN106526547A (en) * | 2016-12-09 | 2017-03-22 | 北京环境特性研究所 | InSAR technology-based linear scanning near-field RCS test clutter suppression method |
CN106646402A (en) * | 2016-11-11 | 2017-05-10 | 上海无线电设备研究所 | Electromagnetic pulse radar echo simulation method in offshore area |
-
2018
- 2018-03-08 CN CN201810188886.5A patent/CN108445303B/en active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5225668A (en) * | 1991-06-06 | 1993-07-06 | The United States Of America As Represented By The Secretary Of The Navy | Photonic electromagnetic field sensor apparatus |
CN102054094A (en) * | 2010-11-24 | 2011-05-11 | 南京理工大学 | Fast directional multilevel simulation method for planar microstrip circuit |
CN103675781A (en) * | 2012-09-25 | 2014-03-26 | 中国航天科工集团第二研究院二〇七所 | A method for accurately acquiring near field echo |
CN103198227A (en) * | 2013-04-18 | 2013-07-10 | 南京理工大学 | Electromagnetic scattering analyzing method for superspeed flight targets |
CN103226644A (en) * | 2013-05-17 | 2013-07-31 | 南京理工大学 | Electromagnetic scattering property simulating method based on cylindrical surface equivalent source domain decomposition |
CN104199026A (en) * | 2014-08-28 | 2014-12-10 | 中国科学院电子学研究所 | Backscattering cross section measuring method based on trace scanning two-dimensional near field imaging |
CN105630740A (en) * | 2014-11-07 | 2016-06-01 | 南京理工大学 | Electromagnetic analysis method on the basis of matrix Taylor series expansion |
CN104730503A (en) * | 2015-03-18 | 2015-06-24 | 中国科学院电子学研究所 | Method for determining influence on scaling by high-resolution SAR reference target RCS and compensation method |
CN106156394A (en) * | 2015-04-22 | 2016-11-23 | 南京理工大学 | Electromagnetic property extracting method based on explicit difference scheme |
CN106093875A (en) * | 2016-06-23 | 2016-11-09 | 中国科学院电子学研究所 | Near field WB-RCS data acquisition and the compression method of denoising is followed the trail of based on base |
CN106646402A (en) * | 2016-11-11 | 2017-05-10 | 上海无线电设备研究所 | Electromagnetic pulse radar echo simulation method in offshore area |
CN106526547A (en) * | 2016-12-09 | 2017-03-22 | 北京环境特性研究所 | InSAR technology-based linear scanning near-field RCS test clutter suppression method |
Non-Patent Citations (3)
Title |
---|
ARIF U¨C¸U¨NCU: "Multilayer Iterative Solutions of Large-Scale Electromagnetic Problems Using MLFMA", 《IV INTERNATIONAL ELECTROMAGNETIC COMPATIBILITY CONFERENCE》 * |
王 湃: "多层快速多极子分析三维复杂目标的谐振区电磁散射特性", 《北京大学学报》 * |
赵延文: "各向异性介质体目标电磁散射分析的快速算法", 《现代防御技术》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110826279A (en) * | 2019-11-07 | 2020-02-21 | 北京环境特性研究所 | Method and device for calculating optical characteristics of solid target |
CN111273281B (en) * | 2020-02-07 | 2021-12-07 | 北京环境特性研究所 | Short-distance electromagnetic measurement method and system |
CN111273281A (en) * | 2020-02-07 | 2020-06-12 | 北京环境特性研究所 | Short-distance electromagnetic measurement method and system |
CN111339650A (en) * | 2020-02-21 | 2020-06-26 | 北京环境特性研究所 | Electromagnetic basis function construction method and device with reduced orthogonal dimension |
CN113408095B (en) * | 2020-03-16 | 2023-08-29 | 启碁科技股份有限公司 | Electromagnetic characteristic analysis method and electronic device |
CN113408095A (en) * | 2020-03-16 | 2021-09-17 | 启碁科技股份有限公司 | Electromagnetic characteristic analysis method and electronic device |
CN112069713B (en) * | 2020-09-10 | 2022-12-23 | 上海无线电设备研究所 | Near-field scattering characteristic modeling method, electronic device and storage medium |
CN112069713A (en) * | 2020-09-10 | 2020-12-11 | 上海无线电设备研究所 | Near-field scattering characteristic modeling method, electronic device and storage medium |
CN113092878A (en) * | 2021-03-31 | 2021-07-09 | 北京环境特性研究所 | Test method and detection device for electromagnetic radiation of W-band environment |
CN113092878B (en) * | 2021-03-31 | 2023-01-20 | 北京环境特性研究所 | Test method and detection device for electromagnetic radiation of W-band environment |
CN113671269A (en) * | 2021-07-05 | 2021-11-19 | 南京航空航天大学 | Spherical near-far field method based on scattering matrix principle |
CN113567941A (en) * | 2021-08-04 | 2021-10-29 | 上海无线电设备研究所 | Target near-field Doppler echo rapid simulation method |
CN113567941B (en) * | 2021-08-04 | 2023-09-19 | 上海无线电设备研究所 | Rapid simulation method for near-field Doppler echo of target |
Also Published As
Publication number | Publication date |
---|---|
CN108445303B (en) | 2020-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108445303A (en) | A kind of Electromagnetic scatter simulated behavior method | |
Fang | Antenna theory and microstrip antennas | |
CN106291129B (en) | Phased array antenna far-field pattern measurement method | |
WO2016124898A1 (en) | Direction finding using signal power | |
Jeng | Near-field scattering by physical theory of diffraction and shooting and bouncing rays | |
CN104615854B (en) | A kind of beam-broadening and side lobe suppression method based on sparse constraint | |
CN103323845B (en) | Image inversion method of non-uniform sampling comprehensive bore diameter radiometer | |
CN106872970B (en) | A kind of multiple target change data transfer rate tracking device and its method based on differential evolution | |
Wan et al. | High Accuracy 2D‐DOA Estimation for Conformal Array Using PARAFAC | |
CN109597048A (en) | Metre wave radar DOA estimation method based on two-dimensional convolution neural network | |
CN106443571A (en) | Conformal mapping estimation of signal parameters via rotational invariance technique DOA estimation method | |
CN109085556A (en) | A kind of high-frequency ground wave radar wave field forming method based on a second order p-ratio | |
CN109212466A (en) | A kind of broadband direction-finding method based on quantum dragonfly mechanism of Evolution | |
Ayestaran et al. | Neural modeling of mutual coupling for antenna array synthesis | |
CN107515841B (en) | Method for calculating electromagnetic wave arrival angle | |
Dayal et al. | Pattern synthesis using accelerated particle swarm optimization | |
CN109818688A (en) | A kind of calibration method of array antenna, equipment, system and computer readable storage medium | |
Henault et al. | Null-steering sensitivity to ground properties in over-the-horizon radar with planar monopole array on finite radial ground plane | |
Jafri et al. | Scattering of electromagnetic plane wave by a circular disk with surface impedance | |
Zhang et al. | A Fast and Efficient Beamforming Algorithm Imitating Plant Growth Gene for Phased Array Antenna | |
CN113486550B (en) | Triplet feed coefficient determination method in semi-physical radio frequency simulation | |
Lin et al. | Simulation of suppression technology of angular glint based on MIMO radar | |
Ahmed et al. | Sampling schemes for accurate reconstruction and computation of performance parameters of antenna radiation pattern | |
Zhang et al. | A wideband digital polarization synthesis method | |
Pometcu et al. | Direction of Arrival Estimation for Radio-Frequency Seekers |
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 |