CN104699876B - The Performance Prediction method of the more sizes of antenna satellite platform, multiband mutual coupling antenna - Google Patents
The Performance Prediction method of the more sizes of antenna satellite platform, multiband mutual coupling antenna Download PDFInfo
- Publication number
- CN104699876B CN104699876B CN201310661990.9A CN201310661990A CN104699876B CN 104699876 B CN104699876 B CN 104699876B CN 201310661990 A CN201310661990 A CN 201310661990A CN 104699876 B CN104699876 B CN 104699876B
- Authority
- CN
- China
- Prior art keywords
- group
- layer
- matrix
- fft
- antenna
- 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
Links
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
The invention discloses a kind of more sizes of antenna satellite platform, the Performance Prediction methods of multiband mutual coupling antenna.Step is:Dummy satellite grid file is obtained, marks interested antenna port;Octree grouping, index sublayer group to father's layer group and father's layer group to the relation of sublayer group are carried out according to the number on discrete side in most sub-layers group to dummy satellite;Improved Electric Field Integral Equation and the right corresponding with antenna feed excitation vector are established, wherein the nearly effect partial of discrete integral equation uses moment method, low frequency contribution part uses multilayer FFT methods, and high frequency effect partial uses Fast multipole;More resolution sparse approximate inverse fore conditions are constructed using near-field region, using diagonal fore condition, broad sense RWG basic functions layer uses sparse approximate inverse fore condition more for basic function layers of differentiating;The mutual coupling matrix between institute's tag antenna port is solved by the current coefficient of BiCGStab solution by iterative method equations, and according to current coefficient.The present invention has the advantages that EMS memory occupation is small, it is parallel to be easy to, efficient.
Description
Technical Field
The invention relates to the technical field of electromagnetic simulation, in particular to a performance estimation method of a multi-size and multi-band mutual coupling antenna of an antenna satellite platform.
Background
The antennas on the satellite are responsible for the functions of ground communication, broadcasting, remote sensing and the like, but the coupling between the antennas can cause mutual interference, thereby influencing the respective normal work. The moment method based on the electric field integral equation has high calculation precision, and the small discrete unknown quantity is the preferential selection of full-wave analysis satellite platform antenna coupling (S.M.Rao, D.R.Wilton, and A.W.Glisson, "Electromagnetic scattering by surfaces of orbit shape," IEEE Trans. antennas propag, vol.30, no.3, pp.409-418, may 1982) but the impedance matrix discretely generated by the moment method is a dense matrix, and the complexity of storing and directly solving the dense matrix equation is O (N), respectively 2 ) And O (N) 3 ) For the problem of loading the multiband antenna on the satellite platform, the solution of the discrete matrix equation far exceeds the computing capability of a common computer. The other difficulty in analyzing the coupling of the multiband antenna loaded by the satellite platform through a moment method comes from matrix ill-condition caused by uneven grid dispersion, and convergence is slow or even not.
The matrix solves the problem. Common diagonal preconditions and algebraic method-based preconditions such as ILU preconditions are inefficient because algebraic-based preconditions are effective for normal grid discretization problems, very large construction preconditions for near-field regions of multi-scale problems consume computational resources, and the preconditions are sometimes not significant in effect. Based on the preconditions technology (f.vipiana, m.a.francvillea, and g.vecchi, "EFIE modeling of high definition multi-scale structures," IEEE trans.antennas processing, vol.57, pp.2362-2374, jul.2010.) of multi-resolution ILU (Incomplete LU, ILU), it has been proven that the matrix equation generated discretely by the multi-scale problem can be solved efficiently. However, in the method, the iterative solution process of the matrix is accelerated by using the traditional FFT method, and the computational complexity of solving the three-dimensional problem by using the FFT method is O (N) 1.5 logN), the computer resource consumption is large for the problem of large unknown quantity of multiband antenna coupling of the antenna platform. And ILU preconditions are not conducive to parallel computing.
The computational complexity of the multi-layer fast multipole method for analyzing the electromagnetic properties of a three-dimensional target is O (NlogN) (J.Song, C.Lu, and W.Chew, "Multi level fast multipole algorithm for electromagnetic calibration by large complex objects," IEEE trans. Antennas Propag., vol.45, no.10, pp.1488-1493, oct.1997.), but the direct use of the fast multipole method for multi-scale problems is less computationally efficient because the fast multipole method requires the finest layer group size to be no less than 0.2 wavelength, and the near field portion below 0.2 wavelength for dense mesh scattering for multi-scale problems is still large. The problem of too large low frequency near field matrix can be solved by low frequency fast multipole (L.J.Jiang and W.C.Chew, "A mixed-form fast multipole algorithm," IEEE trans.antennas Propag., vol.53, no.12, pp.4145-4156, dec.2005.) but the low frequency fast multipole method is very difficult to implement.
Disclosure of Invention
The invention aims to provide a performance estimation method of a multi-size and multi-band mutual coupling antenna of an antenna satellite platform, which has the advantages of high estimation precision and low resource consumption, so that coupling between the antennas can be efficiently obtained, and technical support is provided for the installation position of the satellite platform antenna.
The technical solution for realizing the purpose of the invention is as follows: a performance estimation method for a multi-size and multi-band cross-coupling antenna of an antenna satellite platform comprises the following steps:
step 1, acquiring a satellite model grid file and marking an interested antenna port;
step 2, according to the grid file information of the satellite model, carrying out octree grouping on the satellite model according to the number of the discrete edges in the finest layer group, counting the number of groups containing the discrete edges in each layer, and indexing the relationship from the sub-layer group to the parent layer group and from the parent layer group to the sub-layer group;
step 3, establishing an electric field integral equation and a right side excitation vector corresponding to an antenna feed source, wherein a near acting part of the discrete integral equation adopts a moment method, a low-frequency acting part adopts a multilayer FFT method, and a high-frequency acting part adopts a rapid multipole method;
step 4, constructing a multi-resolution sparse approximate inverse preconditions by using a near field region, wherein a diagonal preconditions is adopted by a multi-resolution basis function layer, and a sparse approximate inverse preconditions is adopted by a generalized RWG basis function layer;
and 5, solving current coefficients of the equation by a BiCGStab iterative method, and solving a mutual coupling matrix between the antenna ports marked in the step 1 according to the current coefficients.
Compared with the prior art, the invention has the following remarkable advantages: (1) The method is characterized in that a mixed multilayer FFT and multilayer fast multipole method is designed to solve the problem of high and low frequencies generated in the mutual coupling estimation process of the multi-size antenna of the satellite platform, and the multilayer FFT method has the advantages of simplicity and low complexity when used for solving a low-frequency region; (2) The near field region construction multi-resolution sparse approximate inverse pre-condition is designed, and the method has the advantages of low resource consumption, easiness in parallel processing and high convergence speed of the multi-size problem discrete equation.
The present invention is described in further detail below with reference to the attached drawing figures.
Drawings
FIG. 1 is a schematic diagram of a three-dimensional octree structure in the method of the present invention.
FIG. 2 is a schematic diagram of the low frequency interaction of the multi-layer FFT approximation in the method of the present invention.
FIG. 3 is a schematic diagram of a multi-size, multi-band antenna configuration of a satellite platform as analyzed by the present invention.
Fig. 4 is a schematic diagram of a port of a coupled antenna array according to the present invention.
Detailed Description
The invention is described in further detail below with reference to the figures and specific embodiments.
The invention relates to a performance estimation method of a multi-size and multi-band mutual coupling antenna of an antenna satellite platform, which is a simulation platform for estimating the multi-size and multi-band antenna coupling of the satellite platform, and is used for accelerating the multi-size and multi-band antenna dispersion generation of the satellite platform based on electric field integral equation dispersion, a multi-layer FFT method and a multi-layer fast multipole method, and has an equation with high and low frequency characteristics; a multi-resolution-Sparse Approximation Inverse (MR-SAI) precondition is used to accelerate an iterative process of discretely generated ill-conditioned equations, and the specific steps are as follows:
step 1, obtaining a satellite model grid file and marking interested antenna ports.
Step 2, according to the grid file information of the satellite model, carrying out octree grouping on the satellite model according to the number of the discrete edges in the finest layer group, counting the number of groups containing the discrete edges in each layer, and indexing the relationship from the sub-layer group to the parent layer group and from the parent layer group to the sub-layer group; with reference to fig. 1, octree grouping is performed on the satellite model according to the number of discrete edges in the finest layer group, which is specifically as follows:
(2.1) firstly, surrounding the satellite model by a cube to define a layer 0, then equally dividing the cube into eight small cubes to define a layer 1, and further equally dividing each small cube into eight small cubes until reaching the L-th layer, so that the number of average discrete edges in each group is not more than 50;
(2.2) the number of groups contained in the first layer is 8 l All groups in each layer are numbered 1 to 8 in sequence according to the position of the group center l In which 1 is<l<L;
(2.3) the l-1 st layer is defined as the parent group of the l-1 st layer, while the l-1 st layer is defined as the child group of the l-1 st layer, its parent group i is indexed by group i p The method comprises converting the number i into a binary sequence, removing the three digits on the right side of the binary sequence, and converting into a decimal system to obtain the number i p 。
And 3, establishing an electric field integral equation and a right excitation vector corresponding to an antenna feed source, wherein a near action part of the discrete integral equation adopts a moment method, a low-frequency action part adopts a multilayer FFT method, and a high-frequency action part adopts a rapid multipole method, so that a matrix vector in an iterative process is multiplied byWhereinA moment method part is adopted for a near action part,for the low-frequency multi-layer FFT portion,part of the method for high frequency multi-layer fast multipole, wherein the distance between the interaction groups is less than 0.3 wavelength is part of the low frequency interaction, and the interaction groups are part of the methodThe high-frequency action part with the wavelength longer than 0.3 adopts a multilayer FFT method in combination with the low-frequency action part shown in FIG. 2, which is as follows:
(3.1) projecting the basis functions on the irregular discrete grid in the group j onto regular cartesian grid points in the group j;
the base functions and the divergence of the base functions on the 1-n irregular discrete grids in the group j are interpolated by an interpolation factor pi of the scalar bit of the group j jA Interpolation factor pi of sum group j vector bit jD Project to group j 1 (M + 1) 3 On regular Cartesian grid points
In the formula (f) 1 (r)…f n (r) n basis functions included in the ith group, M is the number of interpolation points in the corresponding one-dimensional direction,for solving the divergence, s is the area of a discrete grid, and beta is a Lagrange interpolation polynomial;
lagrange polynomial in three-dimensional form at interpolation point uExpressed as the product of x-axis, y-axis and z-axis interpolating polynomials:
in the formula, x, y and z are coordinates of an interpolation point rectangular coordinate system;
(3.2) passing the basis functions on the regular Cartesian grid points in set j throughParent layer group j projected to group j p On the cartesian grid points of the medium rule:
in the formula (I), the compound is shown in the specification,to interpolate from group j to parent layer group j p The interpolation factor of (2);
(3.3) group j p Passing of electric potential on the points of the medium regular cartesian gridParent layer group i transitioning to group j p The regular cartesian grid points of (a) are,is a father layer group i p And parent layer group j p A green function matrix formed by the regular grid points of (a);
(3.4) group i p Passing of electric potential on the points of the medium regular cartesian gridInterpolating to regular cartesian grid points in group i:
in the formula (I), the compound is shown in the specification,for interpolating from group i to parent layer group i p The interpolation factor of (2);
(3.5) interpolating potentials on regular cartesian grid points in group i onto irregular discrete grids in group i:
1 to (M + 1) in group i 3 Regular cartesianInterpolation factor pi by group i scalar bit for potential interpolation at grid points iA Interpolation factor pi of sum group i vector bit iD A basis function interpolated on 1-n irregular discrete grids in the group i, and an impedance matrix Z formed between the group i and the group j ij FFF by monolayer is expressed as:
in the formula, k 0 A Green function matrix G formed by regular Cartesian grid points in group j and regular Cartesian grid points in group i for free space wavenumbers i,j Comprises the following steps:
in the formulaA green's function of free space formed between interpolation points in the group i and the group j;
(3.6) Green's function matrix G between group i and group j regular grid points in equation (7) i,j For the Toeplitz matrix, matrix vector multiplication is accelerated by using FFT, and matrix vector multiplication in the whole matrix equation iteration process is accelerated by single-layer FFT to be:
wherein Z is the whole moment method discrete impedance matrix, I is the current coefficient to be solved,a moment method part is adopted for the near acting part,for the low frequency multi-layer FFT section, FFT () is the FFT forward transform, IFFT () is the FFT inverse transform, Π A Interpolation factor, pi, being scalar bits D Is an interpolation factor of vector bits, G is a Green function matrix, k 0 Is the free space wavenumber;
the impedance matrix formed between group i and group j is approximated by a two-layer FFT as:
the matrix vector multiplication of the low frequency contribution portion is FFT accelerated to:
in the formula II i Interpolation factor of group i, Π j An interpolation factor for group j;
the high-frequency interaction part is accelerated by adopting a multilayer fast multipole method, and the matrix vector accelerated by the final fast integration method is multiplied by:
in the formula (I), the compound is shown in the specification,is part of a high frequency multi-layer fast multipole method.
Step 4, constructing a multi-resolution sparse approximation inverse preconditions by using a near field region, wherein a Diagonal (D) preconditions is adopted by a multi-resolution basis function (MR) layer, and a sparse approximation inverse preconditions is adopted by a generalized RWG (generalized Rao-Wilton-Glisson) basis function layer; the method comprises the following specific steps:
the multi-resolution basis function is a linear combination of RWG basis functions, and the non-zero elements of each row of the transformation matrix T represent linear combination coefficients of the RWG basis functions:
in the formula, T MR Linear combination coefficient matrix from RWG basis function to multi-resolution basis function, T gRWG For a linear combination coefficient matrix of RWG basis functions to generalized RWG basis functions, the basis functions transform a matrix that constructs preconditions into:
Z p =T(Z near )T T (12) Preconditions matrix Z p Written in the form of
In the formula, Z near | MR,MR Is an impedance matrix formed between the MR basis functions, Z near | MR,gRWG And Z near | gRWG,MR Impedance matrices formed for RWG and gRGG basis functions, Z near | gRWG,gRWG An impedance matrix formed between the basis functions of the gRGGs;
using diagonal preconditions D _ non-calculation in correspondence with portions of multi-resolution basis functions MR,MR :
Corresponding to the generalized RWG basis function part, adopting sparse approximate inverse preconditions Z p | gRWG,gRWG Taking the minimum value of:
wherein M is a preconditions matrix, Z e Is an identity matrix, e j Is a unit column vector, m j Is the column vector of M, | is a two-norm;
thereby obtaining a preconditioned matrix equation Z MR-SAI Comprises the following steps:
and 5, solving current coefficients of the equation by using a BiCGStab (coordinated differentiation Stabilized, biCGStab) iteration method, and solving a mutual coupling matrix between the antenna ports marked in the step 1 according to the current coefficients.
The advantages of the multi-layer FFT method over conventional FFT-like methods such as the Adaptive Integration Method (AIM) are: (1) The pre-calibration process is not needed, and compared with the traditional FFT method, the pre-treatment construction time is saved; (2) The method only establishes the regular grid points in the non-empty group, avoids the defect of computing resource waste caused by establishing the regular grid in the whole space by the traditional FFT method, and is suitable for multi-scale uneven target analysis. The advantage of the low-frequency action analysis part of the multilayer FFT method is that the operation is simple, and a complex mathematical processing process similar to the low-frequency rapid multipole method is not needed; compared with a low-rank compression method, the method has the advantages of high efficiency and low complexity; the multi-layer FFT method analyzes the complexity of the low-frequency action part to be O (N).
Example 1
A satellite platform including an antenna array and a communications dipole antenna is electromagnetically simulated according to the method of the present invention as shown in fig. 3. The coupling between the central element (port 1) and the upper right antenna element (port 2) of the array antenna as shown in fig. 4 was analyzed. The simulation frequency is 500MHz, and the electrical size of the corresponding antenna is 21 wavelengths; the subdivision size is changed from 1.7E-3 to 7.2E-2m, and the corresponding electrical size is 0.003 to 0.12 wavelength; the total unknowns of the discrete grid is 212537, but 111794 unknowns are concentrated on the antenna array and have a size of 1 square wavelength, so the subdivision is very non-uniform. The iterative solving step number of the equation by using the method is 657, and the consumption of the memory and the calculation time is 11.0GB and 14 hours. When diagonal preconditions are used, the matrix equation is not converged within 2000 steps of iteration, and if ILU preconditions are directly used, the memory consumption exceeds the maximum memory (96 GB) of a computer, the coupling size of the array center antenna unit and the array corner antenna unit simulated by the method is-70 dB.
Claims (1)
1. A performance estimation method for a multi-size and multi-band mutual coupling antenna of an antenna satellite platform is characterized by comprising the following steps:
step 1, acquiring a satellite model grid file and marking an interested antenna port;
step 2, according to the grid file information of the satellite model, carrying out octree grouping on the satellite model according to the number of the discrete edges in the finest layer group, counting the number of groups containing the discrete edges in each layer, and indexing the relationship from the sub-layer group to the parent layer group and from the parent layer group to the sub-layer group;
step 3, establishing an electric field integral equation and a right side excitation vector corresponding to an antenna feed source, wherein a near acting part of the discrete integral equation adopts a moment method, a low-frequency acting part adopts a multilayer FFT method, and a high-frequency acting part adopts a rapid multipole method;
step 4, constructing a multi-resolution sparse approximate inverse preconditions by using a near field region, wherein a diagonal preconditions is adopted by a multi-resolution basis function layer, and a sparse approximate inverse preconditions is adopted by a generalized RWG basis function layer;
step 5, solving current coefficients of the equation by a BiCGStab iterative method, and solving a cross-coupling matrix between the antenna ports marked in the step 1 according to the current coefficients;
and 2, according to the grid file information of the satellite model, grouping the satellite model into octree groups according to the number of the discrete edges in the finest layer group, specifically as follows:
(2.1) firstly, surrounding the satellite model by a cube to be defined as a layer 0, then equally dividing the cube into eight small cubes to be defined as a layer 1, and further equally dividing each small cube into eight small cubes until reaching the L-th layer, so that the number of average discrete edges in each group is not more than 50;
(2.2) the number of groups contained in the first layer is 8 l All groups in each layer are numbered 1 to 8 in sequence according to the position of the group center l In which 1 is<l<L;
(2.3) layer l-1 is defined as the parent layer group of layer l, and layer l is oppositeA sub-layer group defined as layer l-1, its parent group i being indexed by group i p The method comprises converting the serial number i into a binary sequence, removing the three bits on the right side of the binary sequence, and converting into a decimal number i p ;
And 3, establishing an electric field integral equation and a right side excitation vector corresponding to an antenna feed source, wherein a near acting part of the discrete integral equation adopts a moment method, a low-frequency acting part adopts a multilayer FFT method, and a high-frequency acting part adopts a rapid multipole method, so that a matrix vector in an iterative process is multiplied byWhereinA moment method part is adopted for the near acting part,for the low-frequency multi-layer FFT portion,the method part is a high-frequency multilayer rapid multipole method, and comprises the following steps:
(3.1) projecting the basis functions on the irregular discrete grid in the group j onto regular cartesian grid points in the group j;
the scatter of the basis function and the basis function on the 1-n irregular discrete grids in the group j is obtained through an interpolation factor pi of scalar bits of the group j jA Interpolation factor pi of sum group j vector bit jD Project to group j 1 (M + 1) 3 On regular cartesian grid points
In the formula (f) 1 (r)…f n (r) n basis functions included in the ith group, M is the number of interpolation points in the corresponding one-dimensional direction,for solving the divergence, s is the area of a discrete grid, and beta is a Lagrange interpolation polynomial;
lagrange polynomial in three-dimensional form at interpolation point uExpressed as the product of x-axis, y-axis and z-axis interpolation polynomials:
in the formula, x, y and z are coordinates of an interpolation point rectangular coordinate system;
(3.2) passing the basis functions on the regular Cartesian grid points in group j throughParent layer group j projected to group j p On the cartesian grid points of medium rule:
in the formula (I), the compound is shown in the specification,to interpolate from group j to parent layer group j p The interpolation factor of (2);
(3.3) group j p Passing of electric potential on the points of the medium regular cartesian gridParent layer group i transitioning to group j p On regular cartesian grid points in (a),is a parent layer group i p And parent layer group j p A green function matrix formed by the regular grid points of (a);
(3.4) group i p Passing of electric potential on the points of the medium regular cartesian gridInterpolating to regular cartesian grid points in group i:
in the formula (I), the compound is shown in the specification,for interpolating from group i to parent layer group i p The interpolation factor of (2);
(3.5) interpolating potentials on regular cartesian grid points in group i onto irregular discrete grids in group i:
1 to (M + 1) in group i 3 Interpolation of electric potential on regular cartesian grid points by interpolation factor pi of group i scalar bit iA Interpolation factor pi of sum group i vector bit iD A basis function interpolated on 1-n irregular discrete grids in the group i, and an impedance matrix Z formed between the group i and the group j ij FFF by monolayer is expressed as:
in the formula, k 0 A Green function matrix G formed by regular Cartesian grid points in group j and regular Cartesian grid points in group i for free space wavenumbers i,j Comprises the following steps:
in the formulaA green function of free space formed between interpolation points in the group i and the group j;
(3.6) Green's function matrix G between group i and group j regular grid points in equation (7) i,j For the Toeplitz matrix, matrix vector multiplication is accelerated by using FFT, and matrix vector multiplication in the whole matrix equation iteration process is accelerated by single-layer FFT to be:
wherein Z is the whole moment method discrete impedance matrix, I is the current coefficient to be solved,a moment method part is adopted for the near acting part,for the low frequency multi-layer FFT section, FFT () is the FFT forward transform, IFFT () is the FFT inverse transform, Π A Interpolation factor, pi, being scalar bits D Is an interpolation factor of vector bits, G is a Green function matrix, k 0 Is the free space wavenumber;
the impedance matrix formed between group i and group j is approximated by a two-layer FFT as:
the matrix vector multiplication of the low frequency contribution portion is FFT accelerated to:
in the formula II i Interpolation factor of group i, Π j Into a groupAn interpolation factor of j;
the high-frequency interaction part is accelerated by adopting a multilayer fast multipole method, and the matrix vector accelerated by the final fast integration method is multiplied by:
in the formula (I), the compound is shown in the specification,a high frequency multi-layer fast multipole method part;
in the step 4, a multi-resolution sparse approximate inverse preconditions is constructed by using the near-field action region, a diagonal preconditions is used for the multi-resolution basis function layer, and a sparse approximate inverse preconditions is used for the generalized RWG basis function layer, which are specifically as follows:
the multi-resolution basis function is a linear combination of RWG basis functions, and the non-zero elements of each row of the transformation matrix T represent linear combination coefficients of the RWG basis functions:
in the formula, T MR Linear combination coefficient matrix from RWG basis function to multi-resolution basis function, T gRWG For a linear combination coefficient matrix of RWG basis functions to generalized RWG basis functions, the basis functions transform a matrix that constructs preconditions into:
Z p =T(Z near )T T (12)
preconditions matrix Z p Written in the form of
In the formula, Z near | MR,MR Is an impedance matrix formed between the MR basis functions, Z near | MR,gRWG And Z near | gRWG,MR Is RWG basis function andimpedance matrix formed by the basis functions of gRGWG, Z near | gRWG,gRWG An impedance matrix formed between the basis functions of the gRGGs;
using diagonal preconditions D _ non-calculation in correspondence with portions of multi-resolution basis functions MR,MR :
Adopting sparse approximate inverse preconditions Z corresponding to generalized RWG basis function part p | gRWG,gRWG Taking the minimum value of:
wherein M is a preconditions matrix, Z e Is an identity matrix, e j Is a unit column vector, m j Is the column vector of M, | |. | | is a two-norm;
thereby obtaining a preconditioned matrix equation Z MR-SAI Comprises the following steps:
the low-frequency action part adopts a multilayer FFT method, the high-frequency action part adopts a rapid multipole method, the low-frequency action part is formed by the distance between interaction groups being less than 0.3 wavelength, and the high-frequency action part is formed by the distance between the interaction groups being more than 0.3 wavelength.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310661990.9A CN104699876B (en) | 2013-12-06 | 2013-12-06 | The Performance Prediction method of the more sizes of antenna satellite platform, multiband mutual coupling antenna |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310661990.9A CN104699876B (en) | 2013-12-06 | 2013-12-06 | The Performance Prediction method of the more sizes of antenna satellite platform, multiband mutual coupling antenna |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104699876A CN104699876A (en) | 2015-06-10 |
CN104699876B true CN104699876B (en) | 2018-05-29 |
Family
ID=53346996
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310661990.9A Active CN104699876B (en) | 2013-12-06 | 2013-12-06 | The Performance Prediction method of the more sizes of antenna satellite platform, multiband mutual coupling antenna |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104699876B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109858159B (en) * | 2019-02-01 | 2020-10-27 | 北京理工大学 | Electromagnetic analysis method based on ternary parallel multilayer rapid multipole |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102938492A (en) * | 2012-11-13 | 2013-02-20 | 北京遥测技术研究所 | Three-dimensional embattled miniaturization self-adaptation satellite navigation antenna |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005538623A (en) * | 2002-09-10 | 2005-12-15 | フラクトゥス・ソシエダッド・アノニマ | Combined multiband antenna |
-
2013
- 2013-12-06 CN CN201310661990.9A patent/CN104699876B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102938492A (en) * | 2012-11-13 | 2013-02-20 | 北京遥测技术研究所 | Three-dimensional embattled miniaturization self-adaptation satellite navigation antenna |
Non-Patent Citations (4)
Title |
---|
Cross-Band Mutual Coupling of Monopole Antennas on a Finite-Sized Ground Plane;Zhi Ning Chen, et al;《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》;20130831;第61卷(第8期);第4372-4375页 * |
利用混合模型法计算微带天线间互耦;张军等;《微波学报》;20040930;第20卷(第3期);第51-54页 * |
多尺度问题电磁特性叠层矩量法分析;陈如山等;《微波学报》;20131031;第29卷(第5和6期);第65-68页 * |
机载天线辐射特性并行计算与绘制技术;于荣欢等;《计算机工程》;20110228;第37卷(第3期);第1-3页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104699876A (en) | 2015-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhao et al. | A domain decomposition method for electromagnetic radiation and scattering analysis of multi-target problems | |
Wei et al. | A fast direct matrix solver for surface integral equation methods for electromagnetic wave scattering from non-penetrable targets | |
CN102081690B (en) | MDA (Matrix Decomposition Algorithm)-combined novel SVD (Singular Value Decomposition) method for complex circuit | |
CN106294894B (en) | Finite element boundary integration method for rapidly analyzing electromagnetic scattering characteristics of non-uniform target | |
Wan et al. | Hierarchical matrix techniques based on matrix decomposition algorithm for the fast analysis of planar layered structures | |
Wu et al. | A multiple-grid adaptive integral method for multi-region problems | |
CN106991222B (en) | Low-frequency electromagnetic characteristic simulation method based on laminated matrix decomposition | |
Tsang et al. | Wave scattering with UV multilevel partitioning method: 2. Three-dimensional problem of nonpenetrable surface scattering | |
CN102708229A (en) | Matrix decomposition and novel singular value decomposition combined method for complex layered medium structures | |
Burkholder et al. | Fast dual-MGS block-factorization algorithm for dense MoM matrices | |
Gedney et al. | The locally corrected Nyström method for electromagnetics | |
CN104699876B (en) | The Performance Prediction method of the more sizes of antenna satellite platform, multiband mutual coupling antenna | |
CN111651911B (en) | Lumped element impedance sensitivity rapid calculation method and optimization method | |
CN108038313A (en) | A kind of analysis method of the non-uniform electromagnetic characteristic of scattering of subdivision | |
Delgado et al. | An overview of the evolution of method of moments techniques in modern EM simulators | |
CN112668213A (en) | Rapid simulation analysis method for multi-scale antenna array | |
CN103150415B (en) | Microstrip circuit full-wave method | |
Zhao et al. | An improved full-wave multilevel Green’s function interpolation method with RBF-QR technique for fast field evaluation | |
Ding et al. | Multiresolution preconditioner combined with matrix decomposition–singular value decomposition algorithm for fast analysis of electromagnetic scatterers with dense discretisations | |
Liu et al. | Parallel fast time-domain integral-equation methods for transient electromagnetic analysis | |
CN106156479A (en) | Quickly analyze the nested multilamellar complex point source method of metal target Electromagnetic Scattering Characteristics | |
Wu et al. | Efficient incorporation of a PEC/PMC plane in the multiple-grid adaptive integral method | |
CN104699874B (en) | A kind of Multi-layer matrix compression method for analyzing micro-strip encapsulation interconnection line mutual coupling | |
Bekasiewicz et al. | Low-cost EM-simulation-driven multi-objective optimization of antennas | |
Mittra et al. | Some recent advances in the development of numerically efficient computational electromagnetic techniques |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |