CN113987824B - Antenna isolation full-wave calculation method based on transmitting antenna and receiving antenna line ports - Google Patents

Antenna isolation full-wave calculation method based on transmitting antenna and receiving antenna line ports Download PDF

Info

Publication number
CN113987824B
CN113987824B CN202111316920.0A CN202111316920A CN113987824B CN 113987824 B CN113987824 B CN 113987824B CN 202111316920 A CN202111316920 A CN 202111316920A CN 113987824 B CN113987824 B CN 113987824B
Authority
CN
China
Prior art keywords
antenna
transmitting antenna
receiving
basis function
line
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
CN202111316920.0A
Other languages
Chinese (zh)
Other versions
CN113987824A (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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN202111316920.0A priority Critical patent/CN113987824B/en
Publication of CN113987824A publication Critical patent/CN113987824A/en
Application granted granted Critical
Publication of CN113987824B publication Critical patent/CN113987824B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

The invention discloses an antenna isolation full-wave calculation method based on transmitting antenna and receiving antenna line ports, which comprises the following steps: s1, loading an initialized CAD grid model; s2, mounting a receiving/transmitting antenna on the target grid model; s3, extracting the coordinates of the basis function list, and extracting the serial number of the basis function according to the relation between the coordinates of the center point of the basis function and the feeding coordinates of the receiving antenna and the feeding coordinates of the transmitting antenna; s4, calculating line port excitation and port current of the transmitting antenna and the receiving antenna; and S5, calculating the antenna isolation between the transmitting antenna and the receiving antenna. The feed source of the receiving and transmitting antenna only needs to obtain the point coordinates on the antenna port, and a split point is arranged at a certain position in the middle of a line segment from the starting point to the end point of the feed source for feeding energy, so that the establishment of the line feed port can be completed, the model is very convenient and fast to construct, and the isolation degree is very convenient to calculate.

Description

Antenna isolation full-wave calculation method based on transmitting antenna and receiving antenna line ports
Technical Field
The invention relates to antenna isolation calculation, in particular to an antenna isolation full-wave calculation method based on ports of a transmitting antenna and a receiving antenna wire.
Background
With the rapid development of electronic technology and communication technology, tens of hundreds of electronic systems are integrated on the same platform, and meanwhile, an electronic system may need several or even more than ten antennas operating in different bands to receive or transmit electronic signals, for example, an airplane or a naval vessel may integrate various electronic devices and antennas thereof for performing different functions. The near-far field coupling of different antennas in the same system is very strong, and even the normal operation of each transceiver station can be seriously interfered, so that the performance of an electronic system is degraded, and therefore, how to predict and avoid the interference is very important for the normal operation of communication equipment. In addition, when the antenna has a large transmitting power, the surrounding electronic devices are also strongly interfered and cannot work normally, so that the antenna isolation analysis is an important problem of electromagnetic compatibility. The antenna isolation provides an important index for system electromagnetic compatibility (EMC), is a description of the space radiation coupling of the transceiver device, is a factor to be considered in the overall electromagnetic compatibility design, and provides an important basis for the antenna.
Therefore, in the electromagnetic fast calculation scene, it is significant to be able to accurately calculate the isolation of the antenna by using a numerical method, but when analyzing the line port in the large-size electromagnetic fast calculation field, the traditional MOM cannot directly analyze the isolation of the antenna, and the problem of calculating the isolation between the antennas excited by the line port is difficult to be effectively solved.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides an antenna isolation full-wave calculation method based on the ports of a transmitting antenna and a receiving antenna, the feed source of the transmitting antenna and the receiving antenna only needs to obtain point coordinates on the antenna ports, and a split point is arranged at a certain position in the middle of a line segment from the starting point to the end point of the feed source and used for energy feed-in, so that the establishment of the line feed port can be completed, the model establishment is very convenient and fast, and the isolation calculation is very convenient.
The purpose of the invention is realized by the following technical scheme: the antenna isolation full-wave calculation method based on the ports of the transmitting antenna and the receiving antenna wire comprises the following steps:
s1, loading an initialized CAD grid model:
extracting a line feed transmitting antenna model and a line feed receiving antenna model of the CAD grid model from the nanostran format grid model file, and designating a feed source of a receiving/transmitting antenna; importing an electrically large-size target grid model from a lattice model file in a nastran format;
s2, mounting a receiving/transmitting antenna on the target grid model:
respectively setting working planes of the receiving/transmitting antennas according to the antenna installation positions and directions of the receiving/transmitting antennas on the target, aligning the receiving/transmitting antennas and corresponding feed sources to the respective working planes, and completing the installation of the receiving/transmitting antennas on the target;
s3, extracting the coordinates of the basis function list, and extracting the serial number of the basis function according to the relation between the coordinates of the center point of the basis function and the feeding coordinates of the receiving antenna and the feeding coordinates of the transmitting antenna;
s4, calculating line port excitation and port current of the transmitting antenna and the receiving antenna;
and S5, calculating the antenna isolation between the transmitting antenna and the receiving antenna.
Further, the step S3 includes:
the step S3 includes:
s301. receiving antenna from transmitting antenna,Reading point coordinates and point connections of all grids in a grid model of a feed grid and a target into a memory, and selecting a grid point coordinate list [ Pt ] from a grid point coordinate list according to an rwg basis function form]NnodeAnd Point connection List [ Elem]NelemExtracted basis function list fmM is 1,2, …, nbase, wherein nbase is the total number of basis functions; nnode is the number of grid points, and Nelem is the total number of point connection lists, namely the number of units;
s302, reading a line feed point of a transmitting antenna
Figure GDA0003543393540000021
And an end point
Figure GDA0003543393540000022
Adding feed point coordinate P of transmitting antenna at the end of grid point coordinate listt＝(xt,yt,zt) Wherein
Figure GDA0003543393540000023
Figure GDA0003543393540000024
At this time, the total number of nodes in the grid point coordinate list Pt becomes Nnode +1, and the last node is Pt；
S303. query [ Pt ]]NnodeIn the list
Figure GDA0003543393540000025
And
Figure GDA0003543393540000026
are respectively numbered
Figure GDA0003543393540000027
And
Figure GDA0003543393540000028
known as PtiP oftFor Nnode +1, two line segment cells are added at the end of the point join list Elem: each of the units Nelem +1
Figure GDA0003543393540000029
Line segment to Nnode +1 and Nnode +1 to Nlem +2 cells
Figure GDA00035433935400000210
A line segment of (a); at this time, the total number of the point connection list Elem becomes Nelem +2, and the last two units are two newly added line segments;
s304, adding a new basis function at the end of the rwg basis function list, wherein the basis function is composed of a Nelem +1 unit and a Nelem +2 unit, and the basis function list f is formed at the momentmThe total number of the antenna ports is changed into nbase +1, and the construction of the port basis function of the transmitting antenna is completed;
s305. reading line feed point of receiving antenna
Figure GDA00035433935400000211
And an end point
Figure GDA00035433935400000212
Adding feed point coordinate P of transmitting antenna at the end of grid point coordinate listr＝(xr,yr,zr) Wherein
Figure GDA00035433935400000213
Figure GDA00035433935400000214
At this time, the total number of nodes in the grid point coordinate list Pt becomes Nnode +2, and the last node is Pr；
S306, query [ Pt ]]NnodeIn the list
Figure GDA00035433935400000215
And
Figure GDA00035433935400000216
are respectively numbered
Figure GDA00035433935400000217
And
Figure GDA00035433935400000218
known as PriP ofrFor Nnode +2, two line segment cells are added at the end of the point join list Elem: each of the units Nelem +3
Figure GDA00035433935400000219
Line segment to Nnode +2 and Nnode +2 to Nlem +4 th cell
Figure GDA00035433935400000220
A line segment of (a); at this time, the total number of the point connection list Elem becomes Nelem +4, and the last two units are two newly added line segments;
s307, adding a new basis function at the end of the rwg basis function list, wherein the basis function is composed of a Nelem +3 unit and a Nelem +4 unit, and the basis function list f is formed at the momentmThe total number of the antenna ports is changed into nbase +2, and the construction of the port basis function of the receiving antenna is completed;
s308, reading feed point coordinate P of transmitting antennat＝(xt,yt,zt) And feed point coordinates P of the receiving antennar＝(xr,yr,zr)；
S309, traversing the base function list to find out the base function fmCenter point coordinate (x)m,ym,zm) The number of the basis function coinciding with the coordinates of the feed point of the transmitting antenna line is marked as m0And finding the basis function fmCenter point coordinate (x)m,ym,zm) The number of the basis function coinciding with the coordinate of the feed point of the receiving antenna line is marked as m1。
Further, the step S4 includes:
s401, a pair matrix [ A ]]nbase×nbaseFilling to obtain a system matrix A, a required by moment method calculationmnThe mth row and nth column element of A; the fill formula is as follows:
Figure GDA0003543393540000031
wherein f ismAnd fnIs the mth and nth basis functions, G is a three-dimensional Green's function;
s402, exciting item rhs]nbaseFilling to obtain a right vector required by moment method calculation;
Figure GDA0003543393540000032
s403, calculating and solving Ax ═ rhs to obtain x, and obtaining the current x [ ibase ] on the unknown quantity of each basis function ibase]Ibase is an integer subscript from 1 to nbase; obtaining the excitation of the port of the transmitting antenna line, wherein the current of the port of the transmitting antenna line is as follows:
Figure GDA0003543393540000033
the current at the port of the receiving antenna line is:
Figure GDA0003543393540000034
wherein the subscript t represents transmit antenna excitation, the superscript t or r represents transmit antenna or receive antenna port position, and I represents port current;
s404, re-pairing the matrix [ A']nbase×nbaseFilling to obtain a system matrix A ', a ' required by moment method calculation 'mnRow m and column n of A'; the filling formula is as follows:
Figure GDA0003543393540000035
wherein f ismAnd fnIs the mth and nth basis functions, G is a three-dimensional Green's function;
s405, re-pairing the excitation term [ rhs']nbaseFilling to obtain a right vector required by moment method calculation;
Figure GDA0003543393540000036
s406, calculating and solving the value A 'x-rhs' to obtain x ', and obtaining the current x' on the unknown quantity of each basis function ibase [ ibase]Ibase is 1An integer subscript to nbase; obtaining the excitation of the port of the receiving antenna wire, wherein the current of the port of the transmitting antenna wire is as follows:
Figure GDA0003543393540000037
the current at the port of the receiving antenna line is:
Figure GDA0003543393540000038
where the subscript r denotes receive antenna excitation, the superscript t or r denotes transmit antenna or receive antenna port position, and I denotes port current.
Further, the antenna isolation c (db) between the transmitting antenna and the receiving antenna is calculated in step S5 as follows:
Figure GDA0003543393540000039
wherein the content of the first and second substances,
Figure GDA00035433935400000310
Z0is the characteristic impedance of the antenna system.
The invention has the beneficial effects that: in the invention, the feed source of the receiving and transmitting antenna only needs to obtain the point coordinates on the antenna port, and the split point is arranged at a certain position in the middle of a line segment of the starting point and the end point of the feed source for feeding energy, so that the establishment of the line feed port can be completed, and the model construction is very convenient and fast; the antenna isolation C (dB) obtained by the method has high precision when the length of the wire section meets the requirement of electric small size, and meets the requirement of most engineering calculation; the right vector is filled only by configuring the right quantity corresponding to the basis function at the feeding point as 1, and the right quantities corresponding to other basis functions as 0, so that the filling calculation efficiency of the right vector is high.
Drawings
FIG. 1 is a flow chart of a method of the present invention;
FIG. 2 is a schematic diagram of a wire-fed transmit antenna model;
FIG. 3 is a schematic diagram of a wire-fed receive antenna model;
fig. 4 is a schematic view showing the mounting of the transmitting/receiving antenna on the target.
Detailed Description
The technical solutions of the present invention are further described in detail below with reference to the accompanying drawings, but the scope of the present invention is not limited to the following.
As shown in fig. 1, the method for calculating the full-wave antenna isolation based on the ports of the transmitting antenna and the receiving antenna comprises the following steps:
s1, loading an initialized CAD grid model:
extracting a line feed transmitting antenna model (an ultrashort wave element antenna shown in figure 2) and a line feed receiving antenna model (an ultrashort wave saber antenna shown in figure 3) of the CAD grid model from the nanostran format grid model file, and designating a feed source of a receiving/transmitting antenna; importing an electrically large-size target grid model from a lattice model file in a nanostran format;
s2, mounting a receiving/transmitting antenna on the target grid model:
respectively setting a receiving/transmitting antenna working plane according to the antenna mounting position and direction of the receiving/transmitting antenna on the target, aligning the receiving/transmitting antenna and the corresponding feed source to the respective working plane, and completing the mounting of the receiving/transmitting antenna on the target, as shown in fig. 4;
s3, extracting the coordinates of the basis function list, and extracting the serial number of the basis function according to the relation between the coordinates of the center point of the basis function and the feeding coordinates of the receiving antenna and the feeding coordinates of the transmitting antenna;
s4, calculating line port excitation and port current of the transmitting antenna and the receiving antenna;
and S5, calculating the antenna isolation between the transmitting antenna and the receiving antenna.
Further, the step S3 includes:
s301, reading point coordinates and point connections of all grids from grid models of transmitting antennas, receiving antennas, feed grids and targets into a memory, and according to the rwg basis function form, reading point coordinates and point connections from a grid point coordinate list [ Pt ]]NnodeAnd Point connection List [ Elem]NelemExtract the basis function list fmM is 1,2, …, nbase, wherein nbase is the total number of basis functions; nnode is the number of grid points, Nelem is the point connectionThe total number of the lists, namely the number of units;
s302, reading a line feed point of a transmitting antenna
Figure GDA0003543393540000041
And an end point
Figure GDA0003543393540000042
Adding feed point coordinate P of transmitting antenna at the end of grid point coordinate listt＝(xt,yt,zt) Wherein
Figure GDA0003543393540000051
Figure GDA0003543393540000052
At this time, the total number of nodes in the grid point coordinate list Pt becomes Nnode +1, and the last node is Pt；
S303. query [ Pt ]]NnodeIn the list
Figure GDA0003543393540000053
And
Figure GDA0003543393540000054
are respectively numbered
Figure GDA0003543393540000055
And
Figure GDA0003543393540000056
known as PtSerial number iP oftFor Nnode +1, two line segment cells are added at the end of the point join list Elem: each of the units Nelem +1
Figure GDA0003543393540000057
Line segment to Nnode +1 and Nnode +1 to Nlem +2 cells
Figure GDA0003543393540000058
A line segment of (a); the total number of point-to-point connection lists Elem now becomes Nelem +2, the last two elementsTwo newly added line segments;
s304, adding a new basis function at the end of the rwg basis function list, wherein the basis function is composed of a Nelem +1 unit and a Nelem +2 unit, and the basis function list f is formed at the momentmThe total number of the antenna ports is changed into nbase +1, and the construction of the port basis function of the transmitting antenna is completed;
s305. reading line feed point of receiving antenna
Figure GDA0003543393540000059
And an end point
Figure GDA00035433935400000510
Adding feed point coordinate P of transmitting antenna at the end of grid point coordinate listr＝(xr,yr,zr) Wherein
Figure GDA00035433935400000511
Figure GDA00035433935400000512
At this time, the total number of nodes in the grid point coordinate list Pt becomes Nnode +2, and the last node is Pr；
S306, query [ Pt ]]NnodeIn the list
Figure GDA00035433935400000513
And
Figure GDA00035433935400000514
are respectively numbered
Figure GDA00035433935400000515
And
Figure GDA00035433935400000516
known as PrSerial number iP ofrFor Nnode +2, two line segment cells are added at the end of the point join list Elem: each of the units Nelem +3
Figure GDA00035433935400000517
To Nnode +2Segment and Nnode +2 to Nelem +4 cells
Figure GDA00035433935400000518
A line segment of (a); at this time, the total number of the point connection list Elem becomes Nelem +4, and the last two units are two newly added line segments;
s307, adding a new basis function at the end of the rwg basis function list, wherein the basis function is composed of a Nelem +3 unit and a Nelem +4 unit, and the basis function list f is formed at the momentmThe total number of the antenna ports is changed into nbase +2, and the construction of the port basis function of the receiving antenna is completed;
s308, reading feed point coordinate P of transmitting antennat＝(xt,yt,zt) And feed point coordinates P of the receiving antennar＝(xr,yr,zr)；
S309, traversing the base function list to find out the base function fmCenter point coordinate (x)m,ym,zm) The number of the basis function coinciding with the coordinates of the feed point of the transmitting antenna line is marked as m0And finding the basis function fmCenter point coordinate (x)m,ym,zm) The number of the basis function coinciding with the coordinate of the feed point of the receiving antenna line is marked as m1。
Further, the step S4 includes:
s401. Pair matrix [ A ]]nbase×nbaseFilling to obtain a system matrix A, a required by moment method calculationmnThe mth row and nth column element of A; the fill formula is as follows:
Figure GDA00035433935400000519
wherein f ismAnd fnIs the mth and nth basis functions, G is a three-dimensional Green's function;
s402, exciting item rhs]nbaseFilling to obtain a right vector required by moment method calculation;
Figure GDA0003543393540000061
s403, calculating and solving Ax ═ rhs to obtain x, and obtaining the current x [ ibase ] on the unknown quantity of each basis function ibase]Ibase is an integer subscript from 1 to nbase; obtaining the excitation of the port of the transmitting antenna line, wherein the current of the port of the transmitting antenna line is as follows:
Figure GDA0003543393540000062
the current at the port of the receiving antenna line is:
Figure GDA0003543393540000063
wherein the subscript t represents transmit antenna excitation, the superscript t or r represents transmit antenna or receive antenna port position, and I represents port current;
s404, re-pairing the matrix [ A']nbase×nbaseFilling to obtain a system matrix A ', a ' required by moment method calculation 'mnRow m and column n of A'; the fill formula is as follows:
Figure GDA0003543393540000064
wherein f ismAnd fnIs the mth and nth basis functions, G is a three-dimensional Green's function;
s405, exciting item [ rhs 'is newly paired']nbaseFilling to obtain a right vector required by moment method calculation;
Figure GDA0003543393540000065
s406, calculating and solving the value A 'x-rhs' to obtain x ', and obtaining the current x' on the unknown quantity of each basis function ibase [ ibase]Ibase is an integer subscript from 1 to nbase; obtaining the excitation of the port of the receiving antenna wire, wherein the current of the port of the transmitting antenna wire is as follows:
Figure GDA0003543393540000066
the current at the port of the receiving antenna line is:
Figure GDA0003543393540000067
where the subscript r denotes receive antenna excitation, the superscript t or r denotes transmit antenna or receive antenna port position, and I denotes port current.
Further, the antenna isolation c (db) between the transmitting antenna and the receiving antenna is calculated in step S5 as follows:
Figure GDA0003543393540000068
wherein the content of the first and second substances,
Figure GDA0003543393540000069
Z0the characteristic impedance of the antenna system is typically taken to be 50 ohms.
The above examples are only intended to illustrate the technical solution of the present invention, but not to limit it; although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions of the embodiments of the present invention.

Claims (4)

1. The antenna isolation full-wave calculation method based on the ports of the transmitting antenna and the receiving antenna line is characterized in that: the method comprises the following steps:
s1, loading an initialized CAD grid model:
extracting a line feed transmitting antenna model and a line feed receiving antenna model of the CAD grid model from the nanostran format grid model file, and designating a feed source of a receiving/transmitting antenna; importing an electrically large-size target grid model from a lattice model file in a nastran format;
s2, mounting a receiving/transmitting antenna on the target grid model:
respectively setting receiving/transmitting antenna working planes according to the antenna installation position and direction of the receiving/transmitting antenna on the target, aligning the receiving/transmitting antenna and the corresponding feed source to the respective working planes, and completing the installation of the receiving/transmitting antenna on the target;
s3, extracting the coordinates of the basis function list, and extracting the serial number of the basis function according to the relation between the coordinates of the center point of the basis function and the feeding coordinates of the receiving antenna and the feeding coordinates of the transmitting antenna;
s4, calculating line port excitation and port current of the transmitting antenna and the receiving antenna;
and S5, calculating the antenna isolation between the transmitting antenna and the receiving antenna.
2. The method of claim 1, wherein the antenna isolation full wave calculation is based on ports of a transmitting antenna and a receiving antenna line: the step S3 includes:
s301, reading point coordinates and point connections of all grids from grid models of transmitting antennas, receiving antennas, feed grids and targets into a memory, and according to the rwg basis function form, reading point coordinates and point connections from a grid point coordinate list [ Pt ]]NnodeAnd Point connection List [ Elem]NelemExtract the basis function list fm1, 2.., nbase, wherein nbase is the total number of basis functions; nnode is the number of grid points, and Nelem is the total number of point connection lists, namely the number of units;
s302, reading a line feed point of a transmitting antenna
Figure FDA0003543393530000011
And an end point
Figure FDA0003543393530000012
Adding feed point coordinate P of transmitting antenna at the end of grid point coordinate listt＝(xt，yt，zt) Wherein
Figure FDA0003543393530000013
Figure FDA0003543393530000014
At this time, the total number of nodes in the grid point coordinate list Pt becomes Nnode +1, and the last nodeIs Pt；
S303, inquiring [ Pt ]]NnodeIn the list
Figure FDA0003543393530000015
And
Figure FDA0003543393530000016
are respectively numbered
Figure FDA0003543393530000017
And
Figure FDA0003543393530000018
known as PtiP oftFor Nnode +1, two line segment cells are added at the end of the point join list Elem: each of the units Nelem +1
Figure FDA0003543393530000019
Line segment to Nnode +1 and Nnode +1 to Nlem +2 cells
Figure FDA00035433935300000110
A line segment of (a); at this time, the total number of the point connection list Elem becomes Nelem +2, and the last two units are two newly added line segments;
s304, adding a new basis function at the end of the rwg basis function list, wherein the basis function is composed of a Nelem +1 unit and a Nelem +2 unit, and the basis function list f is formed at the momentmThe total number of the antenna ports is changed into nbase +1, and the construction of the port basis function of the transmitting antenna is completed;
s305. read line of receiving antenna and feed the starting point
Figure FDA00035433935300000111
And an end point
Figure FDA00035433935300000112
Adding feed point coordinate P of transmitting antenna at the end of grid point coordinate listr＝(xr，yr，zr) Wherein
Figure FDA00035433935300000113
Figure FDA0003543393530000021
At this time, the total number of nodes in the grid point coordinate list Pt becomes Nnode +2, and the last node is Pr；
S306, query [ Pt ]]NnodeIn the list
Figure FDA0003543393530000022
And
Figure FDA0003543393530000023
are respectively numbered
Figure FDA0003543393530000024
And
Figure FDA0003543393530000025
known as PriP ofrFor Nnode +2, two line segment cells are added at the end of the point join list Elem: each of the units Nelem +3
Figure FDA0003543393530000026
Line segment to Nnode +2 and Nnode +2 to Nlem +4 th cell
Figure FDA0003543393530000027
A line segment of (a); at this time, the total number of the point connection list Elem becomes Nelem +4, and the last two units are two newly added line segments;
s307, adding a new basis function at the end of the rwg basis function list, wherein the basis function is composed of a Nelem +3 unit and a Nelem +4 unit, and the basis function list f is formed at the momentmThe total number of the antenna ports is changed into nbase +2, and the construction of the port basis function of the receiving antenna is completed;
s308, reading feed point coordinate P of transmitting antennat＝(xt，yt，zt) And feed point coordinates P of the receiving antennar＝(xr，yr，zr)；
S309, traversing the base function list to find out the base function fmCenter point coordinate (x)m，ym，zm) The number of the basis function coinciding with the coordinates of the feed point of the transmitting antenna line is marked as m0And finding the basis function fmCenter point coordinate (x)m，ym，zm) The number of the basis function coinciding with the coordinate of the feed point of the receiving antenna line is marked as m1。
3. The method of claim 1, wherein the antenna isolation full wave calculation is based on ports of a transmitting antenna and a receiving antenna line: the step S4 includes:
s401, a pair matrix [ A ]]nbase×nbaseFilling to obtain a system matrix A, a required by moment method calculationmnThe mth row and nth column element of A; the fill formula is as follows:
Figure FDA0003543393530000028
wherein f ismAnd fnIs the mth and nth basis functions, G is a three-dimensional Green's function; epsilonlIs the free space dielectric constant; mu.slIs free space permeability;
Figure FDA0003543393530000029
in order to be a spatial linear differential operator,
Figure FDA00035433935300000210
omega is angular frequency;
s402, exciting item rhs]nbaseFilling to obtain a right vector required by moment method calculation;
Figure FDA00035433935300000211
s403, calculating and solving Ax ═ rhs to obtain x, and obtaining the current x [ ibase ] on the unknown quantity of each basis function ibase]Ibase is an integer subscript from 1 to nbase; obtaining the excitation of the port of the transmitting antenna line, wherein the current of the port of the transmitting antenna line is as follows:
Figure FDA00035433935300000212
the current at the port of the receiving antenna line is:
Figure FDA00035433935300000213
wherein the subscript t represents transmit antenna excitation, the superscript t or r represents transmit antenna or receive antenna port position, and I represents port current;
s404, re-pairing the matrix [ A']nbase×nbaseFilling to obtain a system matrix A ', a ' required by moment method calculation 'mnRow m and column n of A'; the fill formula is as follows:
Figure FDA0003543393530000031
wherein f ismAnd fnIs the mth and nth basis functions, G is a three-dimensional Green's function; epsilonlIs the free space dielectric constant; mu.slIs free space permeability;
Figure FDA0003543393530000032
in order to be a spatial linear differential operator,
Figure FDA0003543393530000033
omega is angular frequency;
s405, exciting item [ rhs 'is newly paired']nbaseFilling to obtain a right vector required by moment method calculation;
Figure FDA0003543393530000034
s406, calculating and solvingSolving A 'x ═ rhs' to obtain x ', and obtaining the current x' on the unknown quantity of each basis function ibase [ ibase]Ibase is an integer subscript from 1 to nbase; obtaining the excitation of the port of the receiving antenna wire, wherein the current of the port of the transmitting antenna wire is as follows:
Figure FDA0003543393530000035
the current at the port of the receiving antenna line is:
Figure FDA0003543393530000036
where the subscript r denotes receive antenna excitation, the superscript t or r denotes transmit antenna or receive antenna port position, and I denotes port current.
4. The method of claim 1, wherein the antenna isolation full wave calculation is based on ports of a transmitting antenna and a receiving antenna line: the manner of calculating the antenna isolation c (db) between the transmitting antenna and the receiving antenna in step S5 is as follows:
Figure FDA0003543393530000037
wherein the content of the first and second substances,
Figure FDA0003543393530000038
Z0is the characteristic impedance of the antenna system.
CN202111316920.0A 2021-11-09 2021-11-09 Antenna isolation full-wave calculation method based on transmitting antenna and receiving antenna line ports Active CN113987824B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111316920.0A CN113987824B (en) 2021-11-09 2021-11-09 Antenna isolation full-wave calculation method based on transmitting antenna and receiving antenna line ports

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111316920.0A CN113987824B (en) 2021-11-09 2021-11-09 Antenna isolation full-wave calculation method based on transmitting antenna and receiving antenna line ports

Publications (2)

Publication Number Publication Date
CN113987824A CN113987824A (en) 2022-01-28
CN113987824B true CN113987824B (en) 2022-05-03

Family

ID=79747220

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111316920.0A Active CN113987824B (en) 2021-11-09 2021-11-09 Antenna isolation full-wave calculation method based on transmitting antenna and receiving antenna line ports

Country Status (1)

Country Link
CN (1) CN113987824B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115661285B (en) * 2022-10-27 2023-09-22 东莘电磁科技(成都)有限公司 Regular curved surface-borne antenna time-harmonic near-field induction characteristic image generation method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101877022A (en) * 2010-05-21 2010-11-03 西安电子科技大学 Carrier antenna structure meshing method based on method of moment (MOM) for triangular vector basis function
CN106951607A (en) * 2017-03-06 2017-07-14 电子科技大学 A kind of electrically large platform modeling method for being used to calculate isolation between antenna
CN113111524A (en) * 2021-04-21 2021-07-13 北京航空航天大学 Antenna isolation degree prediction method based on receiving and transmitting antenna far field test data, storage medium and device

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160099493A1 (en) * 2014-10-02 2016-04-07 Richard Smith Antenna Device

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101877022A (en) * 2010-05-21 2010-11-03 西安电子科技大学 Carrier antenna structure meshing method based on method of moment (MOM) for triangular vector basis function
CN106951607A (en) * 2017-03-06 2017-07-14 电子科技大学 A kind of electrically large platform modeling method for being used to calculate isolation between antenna
CN113111524A (en) * 2021-04-21 2021-07-13 北京航空航天大学 Antenna isolation degree prediction method based on receiving and transmitting antenna far field test data, storage medium and device

Also Published As

Publication number Publication date
CN113987824A (en) 2022-01-28

Similar Documents

Publication Publication Date Title
CN113987824B (en) Antenna isolation full-wave calculation method based on transmitting antenna and receiving antenna line ports
CN101060199A (en) A design method for compact and lamination structure microstrip antenna integrated with the receiving and sending functions
CN106415316A (en) Method and device for indoor positioning
CN105553584A (en) 3DMIMO channel modeling method
CN102073754B (en) Comprehensive electromechanical analysis method of reflector antenna based on error factor
CN101930495B (en) Planar slot array antenna mechanical and electronic comprehensive analysis method based on grid conversion
CN103065007A (en) Antenna modeling method based on test data of directional diagram and input impedance
CN104573271A (en) Parameter-variable precision plasma antenna radiation performance simulation method
CN112099064A (en) Network type GNSS original data acquisition, transmission and positioning system
CN113987825B (en) Antenna isolation full wave calculation method based on surface element edge port
CN113987826B (en) Waveguide port-based antenna isolation full-wave calculation method
CN109344456A (en) A kind of phase compensation of high integration submatrix and verification method
CN106226664A (en) Transformer station partial discharge signals detecting system
CN111562547A (en) 3D visualization method and system for monitoring element
Darawade et al. A review on dielectric resonator antenna and its analysis setup
CN104462651A (en) Method for quickly determining mounting position of microwave device bolts based on sparse array
Nguyen et al. Genetic algorithm for optimization of L-shaped PIFA antennas
Shalini et al. The Design, Simulation, and Measurement of a 30 GHz Microstrip Patch Antenna for Millimeter-Wave Applications
CN108170883A (en) Manufacturer's method under a kind of spacecraft cable based on structured parameter model
Badulescu Simulation of Inverted F-type Antenna Arrays for Evaluation of Bluetooth Connectivity between a wBMU and wCMUs for Electric Vehicles
CN108768552A (en) A kind of spuious test method of 1800M RF consistencies transmitter and device
Karibayev et al. The Possibilities of Using Fractal Antennas in Modern Wireless Communication Technologies
Bai et al. Wide‐Band High‐Gain DGS Antenna System for Indoor Robot Positioning
Zheng et al. Design of High-Performance Antenna System for High-Speed Railways
CN117648527A (en) Cable crosstalk accurate prediction method based on improved sparrow search algorithm

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