CN110765412B - Method for solving electromagnetic scattering of electric large object by centroid segmentation wavelet moment method - Google Patents

Method for solving electromagnetic scattering of electric large object by centroid segmentation wavelet moment method Download PDF

Info

Publication number
CN110765412B
CN110765412B CN201910917358.3A CN201910917358A CN110765412B CN 110765412 B CN110765412 B CN 110765412B CN 201910917358 A CN201910917358 A CN 201910917358A CN 110765412 B CN110765412 B CN 110765412B
Authority
CN
China
Prior art keywords
matrix
electric field
segmentation
wavelet
centroid
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.)
Active
Application number
CN201910917358.3A
Other languages
Chinese (zh)
Other versions
CN110765412A (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201910917358.3A priority Critical patent/CN110765412B/en
Publication of CN110765412A publication Critical patent/CN110765412A/en
Application granted granted Critical
Publication of CN110765412B publication Critical patent/CN110765412B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Measuring Magnetic Variables (AREA)

Abstract

The invention discloses a method for solving electromagnetic scattering of an electric large object by a centroid segmentation wavelet moment method, and belongs to the field of calculation electromagnetism in electromagnetic engineering. Aiming at the conditions that the formed matrix is complex, the calculated amount is large and errors are easy to occur when the electromagnetic scattering of the electric large-size conductor target is calculated, firstly, triangular dispersion is carried out on the target surface by utilizing RWG basis functions, then, each original triangle is divided into 9 identical sub-triangles, impedance matrix elements are efficiently filled by adopting a mass center segmentation method, and further, the impedance matrix without singular problem in calculation is obtained, secondly, the sparse matrix generated by discrete wavelet transformation is used for acting on the efficiently filled impedance matrix, and the dense impedance matrix is thinned by utilizing the multi-resolution and vanishing moment characteristics of the wavelet basis, so that the calculation time is reduced. Simulation proves that the method provided by the invention can be used for rapidly calculating the electromagnetic scattering characteristics of the large electric object under the condition of ensuring the calculation accuracy.

Description

Method for solving electromagnetic scattering of electric large object by centroid segmentation wavelet moment method
Technical Field
The invention belongs to the field of calculation electromagnetism in electromagnetic engineering, and particularly relates to a method for solving electromagnetic scattering of an electric large object by a centroid segmentation wavelet moment method.
Background
As an algorithm for solving the integral equation, the moment method directly carries out integral solving on the source region, does not need to set boundary conditions and does not generate dispersion errors, so that various electromagnetic field problems can be accurately solved, but the algorithm is limited by green functions, and a complex matrix equation needs to be solved, particularly when electromagnetic scattering of an electric large-size conductor target is calculated, a full-rank dense matrix can be formed, the inversion is complex, the calculated amount is larger, and errors are easy to occur. Based on this, there are several documents currently in the study of the rapid solution of electromagnetic scattering properties of electrically large objects.
J.M.Song, C.C.Lu proposes a fast multistage sub-method (Fast Multipole Method, FMM) and a multilayer fast multistage sub-method (Multilevel Fast Multipole Algorithm, MLFMA) for the problem of conductor scattering. These two algorithms are very effective in solving the free space target scattering problem, and can effectively reduce the time for calculating the impedance matrix, but such methods often require knowledge of the green's function of the actual problem when solving. When the green's function expression is not simple and clear enough, the complexity of formula derivation and programming is increased, and the algorithm has the defect of poor portability, and the formula is often required to be derived and programmed again for different problems. Ding, Z.Li, and R.S Chen propose filling of the matrix with higher order basis functions, which although a good description of the physical properties is possible, inevitably lead to multiple integration operations, and Anhui university Zhang Aikui proposes equivalent dipole moment method (EDM) which corresponds to the effect between two intervals as a function of the current basis and the near field generated by the electric dipole moment. The representation of the impedance matrix elements becomes very simple and the elements can be calculated directly, which avoids double integration, but is limited by the size of the triangular faces of the target surface segmentation
In summary, the calculation problem of the electromagnetic scattering characteristics of the electric large object in the existing literature is still to be further explored, and based on the method, the invention provides a rapid method for solving the electromagnetic scattering characteristics of the electric large object based on a wavelet moment method of centroid segmentation.
Disclosure of Invention
The invention aims to provide a high-precision and high-efficiency method for solving electromagnetic scattering of an electric large object by a wavelet moment method of centroid segmentation.
The aim of the invention is achieved by the following technical scheme:
the invention relates to a quick solving problem of electromagnetic scattering characteristics of an ideal electric large object, which is characterized in that an impedance matrix filled by the ideal electric large object is subjected to discrete processing by adopting a wavelet moment method based on a centroid segmentation method, so that the scattering characteristics of the electric large object can be quickly solved, impedance matrix elements are efficiently filled by utilizing the centroid segmentation method, and then an impedance matrix without singular problem in calculation is obtained, then a sparse matrix generated by discrete wavelet transformation is used for acting on the high-efficiency filled impedance matrix, and dense impedance matrix is thinned by utilizing the multi-resolution and vanishing moment characteristics of a wavelet base, so that the calculation time is reduced.
The method for solving electromagnetic scattering of the electric large object by using a centroid segmentation wavelet moment method specifically comprises the following steps:
step (1) obtaining the relation between a scattering electric field and an induced current according to a bit function and a potential function theory;
step (2), the tangential component of the electric field intensity of the conductor surface is zero according to the ideal conductor boundary condition;
step (3) adopts a gamma method, and uses RWG basis functions as weight functions to obtain a discrete electric field integral equation;
performing triangular division on the surface of the conductor by using RWG basis functions to obtain integral of a function g on an original triangle;
step (5) obtaining a matrix equation;
and (6) transforming the matrix equation.
The step (1) specifically comprises:
the relationship between the scattered electric field and the induced current obtained according to the bit function and the potential function theory is as follows:
Figure SMS_1
wherein A (r) represents the magnetic sagittal position,
Figure SMS_2
and (3) representing an electric sign bit, wherein the expression of A (r) is as follows:
Figure SMS_3
Figure SMS_4
the expression of (2) is as follows:
Figure SMS_5
j (R) represents the equivalent current of the target surface, G (R) represents the free-space Green's function, where the Green's function is
Figure SMS_6
η is the wave impedance in free space, k is the propagation constant in free space, R represents the distance from the field point to the origin, and the integration region is the surface on which the source current is located, i.e. the surface of the diffuser.
The step (2) specifically comprises:
substituting the formula (2) and the formula (3) into the formula (1) to obtain
Figure SMS_7
Wherein L represents an operator for calculating electric field radiation from current, and the tangential component of the electric field intensity of the conductor surface is zero as known from ideal conductor boundary conditions, namely
Figure SMS_8
wherein Ei Represents the incident electric field and η represents the wave impedance.
The step (3) specifically comprises: adopts a gamma method, uses RWG basis function as weight function, and can obtain discrete electric field integral equation
Figure SMS_9
The step (4) specifically comprises: when using centroid distribution, the conductor surface is triangulated using RWG basis functions, and the integral of the function g over the original triangle can be expressed as
Figure SMS_10
Figure SMS_11
Representing the center of each sub-angle, A m Representing the area of the original triangle, the magnetic and electrical vector bits can be rewritten into
Figure SMS_12
Figure SMS_13
The impedance matrix elements can thus be rewritten into
Figure SMS_14
The step (5) specifically comprises: the matrix equation can be obtained as follows
ZI=V (11)
Z, I and V respectively represent impedance matrixes, unknown vectors and excitation vectors, and U is a wavelet transformation matrix on the assumption that U is a non-singular matrix with a size.
The step (6) specifically comprises: the wavelet transformation matrix U is introduced, and a matrix equation can be transformed:
Z'I'=V' (12)
Z'=UZU T I'=(U T ) -1 I V'=UV (13)
t represents the transpose of the matrix.
The invention has the beneficial effects that:
the impedance matrix based on the centroid segmentation method and filled with the impedance matrix is subjected to discrete processing by adopting a wavelet moment method, so that the scattering characteristic of an electric main body can be rapidly solved, and the dense impedance matrix is thinned and thinned by utilizing the multi-resolution and vanishing moment characteristics of the wavelet matrix, so that the calculation time is reduced.
Drawings
FIG. 1 is a schematic diagram of a triangle centroid cut;
FIG. 2 is a rectangular flat plate after being split;
FIG. 3 is an ideal conductor sphere after dissection;
FIG. 4 is a double station RCS of rectangular flat panel;
FIG. 5 is a double station RCS for an ideal conductor sphere;
FIG. 6 is an impedance matrix fill time and calculated total time for a rectangular flat panel;
fig. 7 is an impedance matrix fill time and calculated total time for an ideal conductor sphere.
Detailed Description
The following is a further description of embodiments of the invention, taken in conjunction with the accompanying drawings:
first embodiment:
the electromagnetic scattering of an electric large object is calculated, and the relation between a scattering electric field and an induced current is firstly obtained according to a bit function theory as follows
Figure SMS_15
/>
Wherein A (r) represents the magnetic sagittal position,
Figure SMS_16
representing the electrical index, the expression of A (r)
Figure SMS_17
Figure SMS_18
The expression of (C) is as follows
Figure SMS_19
J (R) represents the equivalent current of the target surface, G (R) represents the free-space Green's function, where the Green's function is
Figure SMS_20
Eta isThe wave impedance in free space, k is the propagation constant in free space, R represents the distance from the field point to the origin, and the integration region is the surface on which the source current is located, i.e. the surface of the scatterer. Substituting the formulas (15) and (16) into (14) can obtain:
Figure SMS_21
l represents an operator for calculating electric field radiation by current, and tangential component of electric field intensity of the surface of the conductor is zero according to ideal conductor boundary conditions, namely:
Figure SMS_22
E i represents the incident electric field and η represents the wave impedance. By adopting a gamma method and using RWG basis functions as weight functions, a discrete electric field integral equation can be obtained:
Figure SMS_23
with centroid distribution, the conductor surface is triangulated using RWG basis functions, each of which is subdivided into nine identical small triangles, as shown in fig. 1. The black dot is located at the centroid of the triangle. The large black points represent the positions of the field points, the small black points represent the positions of the source points, and the situation that the field points and the far points are coincident at the moment can be clearly seen, namely the problem of singularity in the process of calculating impedance elements is avoided, and at the moment, the integration can be approximated by the weighted summation of the central values on each sub-triangle. Thus, the integral of the function g over the original triangle can be expressed as:
Figure SMS_24
Figure SMS_25
representing the center of each sub-angle shape,A m representing the area of the original triangle, the magnetic and electrical vector bits can be rewritten into
Figure SMS_26
Figure SMS_27
The impedance matrix elements can thus be rewritten into
Figure SMS_28
The mass center segmentation method can realize the efficient filling of the impedance matrix while solving the problem of singularity, but has no influence on the calculation of an impedance matrix equation. The root cause is that the moment method introduces a green's function as an integral kernel in calculating the scattering problem. Although it can accurately describe the propagation process of electromagnetic fields, the matrix it produces is a full rank dense matrix, which is computationally complex. Combining wavelet transforms with a moment method can effectively solve this problem. Since the complexity of solving the linear equation is related to the number of non-zero elements in the coefficient matrix, the more zero elements, the more sparse the coefficient matrix, the lower the computational complexity. The electric field integral equation can be solved by a moment method, and the unknown current is also expanded into a superposition form of RWG basis functions. The matrix equation can be obtained as follows
ZI=V (24)
Z, I and V respectively represent impedance matrixes, unknown vectors and excitation vectors, and U is a wavelet transformation matrix on the assumption that U is a non-singular matrix with a size. For any vector V, the wavelet basis has multi-resolution and vanishing moment characteristics. Therefore, u·v corresponds to a decomposition vector. The wavelet transformation matrix U is introduced into a conventional MoM, and the matrix equation can be transformed.
Z'I'=V' (25)
Z'=UZU T I'=(U T ) -1 I V'=UV (26)
T represents the transpose of the matrix, the effect of U is a sparse matrix, and for matrix Z, the matrix after wavelet transformation is UZU T Due to (Z -1 )'=UZ -1 U T As a sparse matrix, therefore, solve for i=u T (Z') -1 V' is greater than solving for i=z -1 V saves computation time, reducing the total computation time to compute an unknown quantity.
The method of the invention is used to calculate the two-dimensional rectangular plate and the two-dimensional RCS of the three-dimensional conductor sphere respectively, and the structure diagram after adopting Matlab function delaunay triangulation is shown in figures 2 and 3, the plate size is 6λ×6λ, the radius of the conductor sphere is 1λ, the incident plane wave is 300MHz, and the incident angle is
Figure SMS_29
The observation angle of the double-station RCS is
Figure SMS_30
Fig. 4 and fig. 5 show simulation results, and it can be seen that the method of the present invention is well matched with the conventional method, and fig. 6 and fig. 7 show comparison between the filling time and the total calculation time of the impedance matrix of the method of the present invention and the conventional method, and it can be seen that the filling time and the total calculation time of the impedance matrix elements of the method proposed by the present invention can be significantly reduced.
Specific embodiment II:
aiming at the problems that an impedance matrix formed when solving the electromagnetic scattering characteristic of an electric large object is complex and difficult to solve directly, the invention provides a calculation method for efficient filling and sparsification of the impedance matrix.
Establishing an electric field integral equation based on an ideal conductor:
(1) By potential function theory, we can get the relationship between the scattering electric field and the induced current as follows:
Figure SMS_31
wherein, by adopting a mass center segmentation method, the magnetic vector position and the electric mark position can be written into
Figure SMS_32
Figure SMS_33
Substituting equations (28) and (29) into (27) to obtain
Figure SMS_34
(2) The ideal conductor surface boundary condition requires zero tangential component of the electric field at the conductor surface, i.e
Figure SMS_35
(3) Adopting a gamma method, taking RWG basis function as a weight function, and dispersing an electric field integral equation to obtain an impedance matrix equation
Z·I=V (32)
Wherein the impedance matrix
Figure SMS_36
Excitation vector
Figure SMS_37
The sparse matrix generated by discrete wavelet transformation acts on the impedance matrix Z in (7) mn Thinning it, i.e
Z'=T·Z·T' I'=T·I V'=T·V (35)
Then
Z'·I'=V' (36)
And (5) performing simulation analysis. Fig. 1 and 2 are schematic diagrams after triangulation of two-dimensional and three-dimensional electric large object models, fig. 3 and 4 are schematic diagrams of radar scattering cross sections, and matrix filling time and total calculation time under two conditions are obtained through simulation, so that accuracy and high efficiency of the method are proved.
The above description is only of the preferred embodiments of the present invention and is not intended to limit the present invention, but various modifications and variations can be made to the present invention 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 (1)

1. The method for solving electromagnetic scattering of the electric large object by using a centroid segmentation wavelet moment method is characterized by comprising the following steps of:
(1) Establishing an electric field integral equation based on an ideal conductor;
through potential function theory, the relation between the scattered electric field and the induced current is:
Figure FDA0004145964960000011
wherein A (r) represents the magnetic vector position,
Figure FDA0004145964960000012
Figure FDA0004145964960000013
the electrical index is represented by the number of bits,
Figure FDA0004145964960000014
j (r) represents an equivalent current of the target surface; g (R) represents a green's function in free space,
Figure FDA0004145964960000015
k is a propagation constant in free space; r represents the distance from the field point to the origin; a (r), ->
Figure FDA0004145964960000016
Is introduced into the expression of (2), resulting in:
Figure FDA0004145964960000017
where η is the wave impedance in free space;
from the ideal conductor boundary conditions, the tangential component of the conductor surface electric field strength is zero, i.e
Figure FDA0004145964960000018
Adopting a gamma gold method, using RWG basis functions as weight functions, and an electric field integral equation based on an ideal conductor is as follows:
Figure FDA0004145964960000019
(2) Triangulating the conductor surface using RWG basis functions using a centroid segmentation method, each divided triangle being subdivided into nine identical small triangles, the integration being approximated by a weighted summation of the central values over each sub-triangle, the integral of the function g over the original triangle being expressed as:
Figure FDA00041459649600000110
wherein ,
Figure FDA00041459649600000111
representing the center of each sub-angle shape; a is that m Representing the area of the original triangle;
the magnetic vector and electrical label bits are rewritten as:
Figure FDA00041459649600000112
Figure FDA00041459649600000113
the impedance matrix element is therefore rewritten as:
Figure FDA0004145964960000021
(3) Solving an electric field integral equation by a moment method, and expanding the unknown current into a superposition form of RWG basis functions;
matrix equation:
ZI=V
wherein Z, I and V respectively represent an impedance matrix, an unknown vector and an excitation vector; introducing a wavelet transformation matrix U, and transforming a matrix equation:
Z′I′=V′
Z′=UZU T I′=(U T ) -1 I V′=UV
the wavelet transformation matrix U is a sparse matrix, and for the matrix Z, the matrix after wavelet transformation is UZU T The method comprises the steps of carrying out a first treatment on the surface of the Due to (Z -1 )'=UZ -1 U T As a sparse matrix, therefore solve for i=u T (Z') -1 V' is greater than solving for i=z -1 And V saves the calculation time and reduces the total calculation time for calculating the unknown quantity.
CN201910917358.3A 2019-09-26 2019-09-26 Method for solving electromagnetic scattering of electric large object by centroid segmentation wavelet moment method Active CN110765412B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910917358.3A CN110765412B (en) 2019-09-26 2019-09-26 Method for solving electromagnetic scattering of electric large object by centroid segmentation wavelet moment method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910917358.3A CN110765412B (en) 2019-09-26 2019-09-26 Method for solving electromagnetic scattering of electric large object by centroid segmentation wavelet moment method

Publications (2)

Publication Number Publication Date
CN110765412A CN110765412A (en) 2020-02-07
CN110765412B true CN110765412B (en) 2023-05-30

Family

ID=69330211

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910917358.3A Active CN110765412B (en) 2019-09-26 2019-09-26 Method for solving electromagnetic scattering of electric large object by centroid segmentation wavelet moment method

Country Status (1)

Country Link
CN (1) CN110765412B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101794356A (en) * 2010-03-26 2010-08-04 中国人民解放军空军工程大学 Computer extranuclear computing method for electromagnetic radiation and scattering of electric large object
CN101877022A (en) * 2010-05-21 2010-11-03 西安电子科技大学 Carrier antenna structure meshing method based on method of moment (MOM) for triangular vector basis function
CN103995935A (en) * 2014-05-26 2014-08-20 南京航空航天大学 Sparsification multi-layer adaptive crossover approximation method for analyzing conductor electromagnetic scattering
CN104914425A (en) * 2015-06-15 2015-09-16 哈尔滨工程大学 Super electrically large size strong electromagnetic pulse environment time-frequency-space multi-dimensional analysis model
CN106529082A (en) * 2016-12-02 2017-03-22 上海无线电设备研究所 Method for rapidly calculating electromagnetic scattering characteristics of electrically large targets

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2478591B1 (en) * 2009-09-16 2020-05-06 Agence Spatiale Européenne Aperiodic and non-planar array of electromagnetic scatterers and reflectarray antenna comprising the same

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101794356A (en) * 2010-03-26 2010-08-04 中国人民解放军空军工程大学 Computer extranuclear computing method for electromagnetic radiation and scattering of electric large object
CN101877022A (en) * 2010-05-21 2010-11-03 西安电子科技大学 Carrier antenna structure meshing method based on method of moment (MOM) for triangular vector basis function
CN103995935A (en) * 2014-05-26 2014-08-20 南京航空航天大学 Sparsification multi-layer adaptive crossover approximation method for analyzing conductor electromagnetic scattering
CN104914425A (en) * 2015-06-15 2015-09-16 哈尔滨工程大学 Super electrically large size strong electromagnetic pulse environment time-frequency-space multi-dimensional analysis model
CN106529082A (en) * 2016-12-02 2017-03-22 上海无线电设备研究所 Method for rapidly calculating electromagnetic scattering characteristics of electrically large targets

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈增帅 ; 赖奔 ; 黄春明 ; 梁昌洪 ; .电大散射问题的复合场积分方程快速傅里叶变换算法.西安电子科技大学学报.(第03期),第463-467页. *

Also Published As

Publication number Publication date
CN110765412A (en) 2020-02-07

Similar Documents

Publication Publication Date Title
CN106295120B (en) A kind of electromagnetic characteristic of scattering quick calculation method
Djordjevic et al. Higher order hybrid method of moments-physical optics modeling technique for radiation and scattering from large perfectly conducting surfaces
Liu et al. An efficient iterative MoM-PO hybrid method for analysis of an onboard wire antenna array on a large-scale platform above an infinite ground
Taygur et al. A bidirectional ray-tracing method for antenna coupling evaluation based on the reciprocity theorem
Zhang et al. MLFMA-PO hybrid technique for efficient analysis of electrically large structures
CN110765412B (en) Method for solving electromagnetic scattering of electric large object by centroid segmentation wavelet moment method
CN112446152B (en) Antenna far-field directional pattern analysis method based on infinitesimal small dipole model deformation array
Huang et al. Solid-angle error in the magnetic-field integral equation for perfectly electric conducting objects
CN105630740A (en) Electromagnetic analysis method on the basis of matrix Taylor series expansion
Sewraj et al. Directional method to compute reduced matrix system in MBF solvers
CN105184033B (en) Broadband electromagnetic property rapid Estimation method based on order stepping metal target
CN106503349B (en) One type periodic structure electromagnetic characteristic of scattering quick calculation method
Zhao et al. Modeling with NURBS surfaces used for the calculation of RCS
CN111368398B (en) Electromagnetic scattering characteristic analysis method and device for electrically large target with uncertain structure
Ahmad et al. Spectral domain fast multipole method for solving integral equations of electromagnetic wave scattering
Yin et al. Study on calculation and verification of radiowave propagation using parabolic equation for the antenna near the ground
Sun et al. A novel method for fast analysis of electromagnetic scattering characteristics of pec objects
Kolundzija et al. Ultra High Order Basis Functions in Analysis of Scattering from Large Metallic Structures
Guo et al. Analysis of antenna array with arbitrarily shaped radomes using fast algorithm based on VSIE
Xiang et al. Application of barycentric subdivision method for singularity integration in method of moments
Garcıa et al. Analysis of the parameters of an approach that combines the characteristic basis function method and the multilevel fast multipole
Huang et al. Numerical analysis of scattering by dielectric random rough surfaces using modified SMCG scheme and curvilinear RWG basis functions
Liu et al. A hybrid method of FEBI and PO for scattering analysis of inhomogeneous structures with large-scale platform
Delgado et al. Fast monostatic RCS computation using the near-field sparse approximate inverse and the multilevel fast multipole algorithm
Hu et al. H-Matrix Accelerated Direct Matrix Solver using Chebyshev-based Nyström Boundary Integral Equation 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