CN113625062A - Antenna housing electrical property estimation method based on Taylor expansion method - Google Patents
Antenna housing electrical property estimation method based on Taylor expansion method Download PDFInfo
- Publication number
- CN113625062A CN113625062A CN202110865829.8A CN202110865829A CN113625062A CN 113625062 A CN113625062 A CN 113625062A CN 202110865829 A CN202110865829 A CN 202110865829A CN 113625062 A CN113625062 A CN 113625062A
- Authority
- CN
- China
- Prior art keywords
- antenna
- antenna housing
- dielectric constant
- radome
- electrical property
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
- G01R29/08—Measuring electromagnetic field characteristics
- G01R29/10—Radiation diagrams of antennas
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
The invention discloses an antenna housing electrical property estimation method based on a Taylor expansion method, which comprises the following steps: generating a grid of an antenna housing model, generating an antenna housing thickness matrix and a dielectric constant matrix, establishing an antenna housing electrical property analysis model, estimating the electrical property after the antenna housing thickness and the dielectric constant are changed, judging whether the antenna electrical property after the antenna housing thickness and the dielectric constant are changed meets the antenna electrical property design requirement, if so, outputting a far field directional diagram and an electrical property index, otherwise, generating the antenna housing thickness matrix and the dielectric constant matrix again and repeating the process until the electrical property design requirement is met. According to the antenna housing electrical property analysis method, the antenna housing thickness matrix and the dielectric constant are respectively substituted into the antenna housing electrical property analysis model, the calculated amount is greatly reduced, and the real-time estimation of the antenna housing electrical property when the antenna housing thickness and the dielectric constant change can be realized.
Description
Technical Field
The invention belongs to the technical field of communication, and further relates to an antenna cover electrical property estimation method based on a Taylor expansion method in the technical field of antennas. The method can be used for estimating the influence of the physical property parameters of the antenna housing on the electrical property of the antenna after the physical property parameters of the antenna housing are changed.
Background
The radome is an important component of a radar system, and the radome provides an all-weather working environment for a radar antenna. For foundation and ship-borne radars, the antenna housing enables the radar to work with high precision under various severe weather conditions, the reliability and the service life of the radar antenna are greatly improved, and the maintenance and repair cost is reduced. When the antenna housing is preliminarily designed, firstly, the margin range of physical parameters is determined, the physical parameters of the antenna housing are subjected to fluctuation type change inevitably under the influence of subsequent processing and manufacturing errors and severe environment in the service process, the electrical performance of an antenna far field is further influenced, a far field directional diagram is distorted, and the antenna housing is subjected to high-precision antenna-housing integrated system due to the high processing and manufacturing cost.
Qin Yujian, a published article "Analysis of tangent based on Matlab" (published as IEEE, published date 2006.03.20), discloses a method for estimating electrical properties after calculating the change in dielectric parameters of a lidding material using planar spectral integration. According to the method, an electromagnetic field of the inner surface of the antenna housing is calculated by performing quadruple integration on a plane wave spectrum, then a tangential electromagnetic field of the outer surface of the housing wall is calculated, and finally a curved surface integral is performed on the tangential field of the outer surface of the antenna housing along the outer surface of the antenna housing. The method has the disadvantages that the quadruple integration is carried out on the plane wave spectrum, wherein the quadruple integration causes a complex calculation process to greatly reduce the calculation efficiency, so that when the electrical property of the antenna housing is predicted, the method is complex, and the same calculation step is repeatedly carried out for many times, so that the procedure is complex when the electrical property of the antenna housing is estimated.
The patent document "radome electrical property prediction method based on far field" (application number 201210009491.7, application date 2012.01.03, and publication number CN 102590656B) applied by the university of west ampere electronics science and technology discloses a radome electrical property prediction method based on far field. The method comprises the steps of firstly establishing a geometric model of an antenna and a feed source according to structural parameters of the parabolic antenna, dividing the antenna model into grids according to one fifth of wavelength, solving a far-field value of the antenna through a physical optical method, then taking the obtained far-field value of the antenna as excitation, solving a two-dimensional far-field value of the antenna after covering the antenna through the geometric optical method, and finally comparing a two-dimensional far-field directional diagram and electrical performance indexes of the antenna to obtain an electrical performance change value of the antenna after covering the antenna. The method has the disadvantages that when the method is used for predicting the electrical property of the antenna after the antenna is covered under the severe environment, the physical property parameter of the antenna cover can be repeatedly changed, and the integral equation of a far field needs to be recalculated once every time the physical property parameter is changed, so that the method has the problems of large calculation time and calculation amount, and is not suitable for estimating the situation that the physical property parameter of the antenna cover is changed for many times.
Disclosure of Invention
The invention aims to provide an antenna housing electrical property estimation method based on a Taylor expansion method aiming at the defects of the prior art, and the method is used for solving the problem that the existing method for calculating the electrical property of the antenna housing occupies a large amount of calculation time and calculation amount when the thickness and the dielectric constant of the antenna housing are analyzed in real time.
The specific idea for achieving the purpose of the invention is that the antenna housing thickness matrix and the dielectric constant matrix are respectively formed by all variable quantities of the antenna housing thickness and the dielectric constant change caused by the real-time measurement of the erosion of the natural environment, the antenna housing thickness matrix and the dielectric constant matrix are directly substituted into the antenna housing electrical property analysis model, the electrical property of the antenna with the antenna housing is obtained, the calculation time and the calculation amount of the prior art for analyzing the electrical property of the antenna housing when the thickness and the dielectric constant change in real time are reduced, the complex program for estimating the electrical property of the antenna housing is simplified, and the antenna housing thickness matrix and the dielectric constant change are suitable for the situation when the physical property parameters of the antenna housing are repeatedly changed. Aiming at the condition that the process for analyzing the change of the thickness and the dielectric constant of the antenna housing is complex and is not suitable for estimating the repeated change of the physical property parameters of the antenna housing in the prior art, the invention analyzes the electrical property when the thickness and the dielectric constant of the antenna housing are changed by combining the antenna housing electrical property analysis model and the real-time monitoring result, and respectively obtains the field intensity and the electrical property index of an antenna far field when the thickness and the dielectric constant of the antenna housing are changed.
In order to achieve the purpose, the technical scheme adopted by the invention comprises the following steps:
with the wavelength length of the electromagnetic wave of 1/5 as the size of each grid cell, randomly dividing the radome model into a plurality of triangular cell grids;
(2a) respectively calculating transmission coefficients in the vertical polarization direction and the horizontal direction of the electromagnetic waves on each unit grid;
(2b) calculating the transmission coefficient of each unit grid along the main polarization direction of the electromagnetic waves;
(2c) calculating far field intensities of the antenna with the antenna housing in spatial observation directions with different pitch angles and azimuth angles, drawing the far field intensities in all the spatial observation directions into a far field directional diagram, and extracting electrical performance indexes of the antenna with the antenna housing;
respectively forming an antenna housing thickness matrix and a dielectric constant matrix by using all variable quantities of antenna housing thickness and dielectric constant change caused by natural environment erosion measured in real time;
step 4, establishing an antenna housing electrical property analysis model:
(4a) the antenna housing electrical property analysis model when the thickness of the antenna housing changes is as follows:
wherein, F'θ,φThe field intensity of the antenna with the radome after the nth column element in the thickness matrix is added in the far field in the space observation direction of the pitch angle theta azimuth angle phi,
it is shown that the operation of derivation is performed,
representing the field intensity of the antenna in the far field in the space observation direction of the pitch angle theta and the azimuth angle phi before the change of the antenna housing, h representing the dielectric constant of the antenna housing [. cndot.)]TDenotes a transpose operation, Δ Xh,nRepresenting the variation of the thickness of the antenna housing at the nth time;
(4b) the antenna housing electrical property analysis model when the antenna housing dielectric constant changes is as follows:
wherein, F ″)θ,φRepresenting the field intensity of the antenna with the radome after the addition of the nth row element in the dielectric constant matrix in the far field in the space observation direction of the pitch angle theta azimuth angle phi, wherein epsilon represents the dielectric constant of the radome, and delta Xε,nRepresenting the variation of the dielectric constant of the antenna housing at the nth time;
step 5, estimating the electrical property after the thickness and the dielectric constant of the antenna housing are changed:
substituting the radome thickness matrix and the dielectric constant matrix into a radome electrical property analysis model, respectively calculating the antenna far field intensity after the radome thickness and the dielectric constant are changed, drawing a far field directional diagram after the radome thickness and the dielectric constant are changed, and extracting electrical property indexes after the radome thickness and the dielectric constant are changed;
step 6, comparing the electrical performance indexes before and after the thickness and the dielectric constant of the antenna housing are changed;
step 7, judging whether the antenna electrical property after the thickness and the dielectric constant of the antenna housing are changed meets the design requirement of the antenna electrical property, if so, executing step 8, otherwise, executing step 3;
and 8, outputting far field directional diagrams and electrical performance indexes before and after the thickness and the dielectric constant of the antenna housing change.
Compared with the prior art, the invention has the following advantages:
1, because the antenna housing electrical property analysis model is established, the far field intensity of the antenna is calculated only after the antenna housing is added once, and the thickness variation and the dielectric constant variation of the antenna housing are respectively substituted into the antenna housing electrical property analysis model, so that the electrical property after the antenna housing is changed is obtained, the step of analyzing the electrical property after the antenna housing is changed is simplified, the problem of the complex procedure of analyzing the electrical property when the thickness variation and the dielectric constant variation of the antenna housing are analyzed in the prior art is solved, and the antenna housing electrical property analysis model has the advantages of reducing the occupied calculation amount and the calculation time of analyzing the electrical property when the thickness variation and the dielectric constant variation of the antenna housing are analyzed.
2, because all the variable quantities of the thickness and the dielectric constant of the antenna housing caused by the erosion of the natural environment measured in real time respectively form an antenna housing thickness matrix and a dielectric constant matrix, the antenna housing thickness matrix and the dielectric constant matrix are substituted into an antenna housing electrical property analysis model, whether the telecommunication energy can meet the design requirements of the antenna electrical property when the thickness and the dielectric constant of the antenna housing change is judged in real time, the problem that the electrical property of the antenna housing cannot be estimated in time when the thickness and the dielectric constant of the antenna housing change are monitored in real time in the prior art is solved, and the method has the advantage that the electrical property estimation can be carried out when the thickness and the dielectric constant of the antenna housing change are monitored in real time.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a comparison graph of antenna far field patterns corresponding to thickness changes of an antenna housing in simulation experiment 1 of the present invention;
fig. 3 is a comparison graph of antenna far field patterns corresponding to the dielectric constant change of the radome in simulation experiment 2 of the present invention.
Detailed Description
The invention is described in further detail below with reference to the figures and the specific embodiments.
The implementation steps of the present invention are described in further detail with reference to fig. 1.
The radome model was randomly divided into a plurality of triangular cell meshes using a commercially available mesh division software with an electromagnetic wave wavelength length of 1/5 as the size of each mesh cell.
And 2, generating a far field intensity directional diagram of the antenna with the antenna housing.
And respectively calculating the transmission coefficients of each unit grid along the vertical polarization direction and the horizontal direction of the electromagnetic wave according to the following formula:
wherein, TVmRepresents a transmission coefficient in the horizontal polarization direction of the electromagnetic wave, T, on the m-th unit cellHmDenotes a transmission coefficient in a direction perpendicular to the polarization direction of the electromagnetic wave on the m-th cell grid, ch (·) denotes a hyperbolic cosine operation, j denotes an imaginary unit, V denotes an intermediate variable, and V ═ 2 pi (epsilon-sin)2θm)1/2Table of/lambda, piCircumference ratio, ε represents the dielectric constant of the radome, sin represents the sinusoidal operation, θmThe incident angle of the electromagnetic wave on the m-th cell grid, λ represents the wavelength of the frequency at which the antenna operates,
c denotes the speed of light, f denotes the frequency at which the antenna operates, Z1Representing the wave impedance in free space.
And calculating the transmission coefficient of each unit grid along the main polarization direction of the electromagnetic wave according to the following formula:
wherein, TmDenotes the transmission coefficient in the main polarization direction of the electromagnetic wave on the m-th grid, cos denotes the cosine operation, betamRepresenting the polarization angle of the incident electromagnetic wave on the m-th grid, exp representing the exponential operation with e as the base, ηHmRepresents the transmission coefficient T in the vertical polarization direction of the electromagnetic wave on the m-th gridHmThe phase of (a) is determined,
as indicated by the intermediate variable(s),
δ′m＝ηHm-ηVm，ηVmrepresents the transmission coefficient T in the horizontal polarization directionVmThe phase of (c).
According to the following formula, calculating the far field intensity of the antenna after the antenna housing in the space observation directions with different pitch angles and azimuth angles, drawing the far field intensity in all the space observation directions into a far field directional diagram, and extracting the electrical performance index of the antenna after the antenna housing is added:
wherein S' represents the area of the aperture field of the antenna, E0Representing the field strength at the aperture field of the antenna, e(·)Denotes an exponential operation with a natural constant e as base, k0Which represents the propagation constant of the free space,
xmand ymRespectively representing the position coordinates of the m-th grid projected on the antenna aperture field.
And the electrical performance indexes of the antenna after the antenna housing is added comprise antenna gain, a left pair lobe level and a right pair lobe level before the antenna housing is changed.
And 3, generating an antenna housing thickness matrix and a dielectric constant matrix.
And respectively forming a radome thickness matrix and a dielectric constant matrix by using all the variable quantities of the radome thickness and the dielectric constant change caused by the real-time measurement of the erosion of the natural environment.
And 4, establishing an antenna housing electrical property analysis model.
The antenna housing electrical property analysis model when the thickness of the antenna housing changes is as follows:
wherein, F'θ,φThe field intensity of the antenna with the radome after the nth column element in the thickness matrix is added in the far field in the space observation direction of the pitch angle theta azimuth angle phi,
it is shown that the operation of derivation is performed,
representing the field intensity of the antenna in the far field in the space observation direction of the pitch angle theta and the azimuth angle phi before the change of the antenna housing, h representing the dielectric constant of the antenna housing [. cndot.)]TDenotes a transpose operation, Δ Xh,nThe variation of the thickness of the radome at the nth time is shown.
The antenna housing electrical property analysis model when the antenna housing dielectric constant changes is as follows:
wherein, F ″)θ,φRepresenting the field intensity of the antenna with the radome after the addition of the nth row element in the dielectric constant matrix in the far field in the space observation direction of the pitch angle theta azimuth angle phi, wherein epsilon represents the dielectric constant of the radome, and delta Xε,nThe change amount of the dielectric constant of the radome at the nth time is shown.
And 5, estimating the electrical property of the antenna housing after the thickness and the dielectric constant of the antenna housing are changed.
And substituting the radome thickness matrix and the dielectric constant matrix into a radome electrical property analysis model, respectively calculating the antenna far field intensity after the radome thickness and the dielectric constant are changed, drawing a far field directional diagram after the radome thickness and the dielectric constant are changed, and extracting electrical property indexes after the radome thickness and the dielectric constant are changed.
And extracting the electrical performance indexes after the thickness and the dielectric constant of the antenna housing change, wherein the electrical performance indexes comprise the gain, the left-side lobe level and the right-side lobe level of the antenna after the thickness and the dielectric constant of the antenna housing change, and the gain, the left-side lobe level and the right-side lobe level of the antenna after the dielectric constant of the antenna housing change.
And 6, comparing the electrical performance indexes before and after the thickness and the dielectric constant of the antenna housing change according to the following formula:
wherein, Δ Gh,nRepresents the gain loss value, Δ G, of the antenna after the nth radome thickness changeε,nShows the gain loss value G of the antenna after the dielectric constant of the antenna housing changes for the nth time1Denotes the gain, Δ L, of the antenna before the change of the radomeh,nThe left side lobe level rise value L of the antenna after the nth antenna housing thickness change is represented1Indicating the left side lobe level, L, of the antenna before the radome changesh,nRepresents the left side lobe level, Δ L, of the antenna after the nth radome thickness changeε,nThe level rise value L of the left side lobe of the antenna after the dielectric constant of the antenna housing is changed for the nth timeε,nRepresents the left side lobe level, Δ R, of the antenna after the nth radome dielectric constant changeh,nRepresents the right minor lobe level rise value, R, of the antenna after the nth radome thickness change1Indicating the right side lobe level, R, of the antenna before the radome changeh,nRepresents the right minor lobe level, Δ R, of the antenna after the nth radome thickness changeε,nRepresents the right minor lobe level rise value, R, of the antenna after the dielectric constant of the antenna housing changes for the nth timeε,nThe right-lobe level of the antenna after the nth radome permittivity change is shown.
And 7, judging whether the antenna electrical property after the thickness and the dielectric constant of the antenna cover are changed meets the design requirement of the antenna electrical property, if so, executing the step 8, and otherwise, executing the step 3.
And 8, outputting far field directional diagrams and electrical performance indexes before and after the thickness and the dielectric constant of the antenna housing change.
The effect of the present invention is further explained by combining the simulation experiment as follows:
1. simulation experiment conditions are as follows:
the hardware platform of the simulation experiment of the invention is as follows: the processor is an Intel i 58500 k CPU, the main frequency is 3.0GHz, and the memory is 16 GB.
The software platform of the simulation experiment of the invention is as follows: windows 10 operating system and Matlab 2020 b.
The embodiment used in the simulation experiment of the invention is a reflector antenna with the working center frequency of 16GHz, an additional antenna housing of the antenna adopts a single-layer medium structure, the aperture of the antenna and the diameter of the antenna housing are both 1.2m, the dielectric constant and the loss tangent value of the antenna housing are respectively 4 and 0.01568, and the thickness of the antenna housing is 1.2 mm.
2. Simulation content and result analysis thereof:
the simulation experiment of the invention is to respectively analyze the electrical property of the antenna with the antenna housing after the antenna housing is added by adopting the method and the prior art (geometric optics-ray tracing method) to obtain the directional diagram and the electrical property index of the far field of the antenna with the antenna housing.
The geometrical optics-ray tracing method in the prior art in the simulation experiment refers to the electrical property analysis method of the antenna with the antenna housing, which is called GO-RT calculation method for short, proposed by TRICOLES G et al in the fields of Radiation patterns of a microwave antenna enclosed by a hollow dielectric wedge [ J ]. Journal of the optical facility of America,1963,53(5):545-557.
The effect of the present invention will be further described below with reference to the simulation diagrams of fig. 2 and 3.
Fig. 2 is a comparison diagram of antenna far-field patterns corresponding to antenna covers with thicknesses of 1.2mm and 2.2mm respectively in simulation experiment 1 of the present invention. The diagram is obtained by uniformly selecting 200 sampling points in the range of a far field observation angle of [ -3 degrees, 3 degrees ], calculating the far field intensity power of the antenna on each sampling point after the antenna housing is added, and normalizing the maximum value of the far field intensity power per se. The abscissa represents the observation angle in degrees in the far-field direction of the antenna, and the ordinate represents the normalized far-field strength power in dB. The curve marked by a solid line in fig. 2 represents a curve of all normalized far-field intensity powers calculated by the prior art when the radome thickness is 1.2mm within the observation angle range of [ -3 °, 3 ° ], the curve marked by a dot-dash line represents a curve of all normalized far-field intensity powers calculated by the prior art when the radome thickness is 2.2mm within the observation angle range of [ -3 °, 3 ° ], and the curve marked by a dotted line represents a curve of all normalized far-field intensity powers calculated by the method of the present invention when the radome thickness is 2.2mm within the observation angle range of [ -3 °, 3 ° ].
Fig. 3 is a comparison graph of far-field patterns of the antenna corresponding to the antenna cover with dielectric constants of 4 and 5 respectively in simulation experiment 2 of the present invention. The diagram is obtained by uniformly selecting 200 sampling points in the range of a far field observation angle of [ -3 degrees, 3 degrees ], calculating the far field intensity power of the antenna on each sampling point after the antenna housing is added, and normalizing the maximum value of the far field intensity power per se. The abscissa represents the observation angle in degrees in the far-field direction of the antenna, and the ordinate represents the normalized far-field strength power in dB. The curve marked by a solid line in fig. 3 represents the curve of all normalized far-field intensity powers when the radome dielectric constant is 4 in the range of the calculated observation angle of the prior art of [ -3 °, 3 ° ], the curve marked by a dot-dash line represents the curve of all normalized far-field intensity powers when the radome dielectric constant is 5 in the range of the calculated observation angle of the prior art of [ -3 °, 3 ° ], and the curve marked by a dotted line represents the curve of all normalized far-field intensity powers when the radome dielectric constant is 5 in the range of the calculated observation angle of [ -3 °, 3 ° ].
It can be seen from fig. 2 and 3 that the curves respectively drawn by the result of the present invention and the result of the prior art are basically coincident, and the graphs show good consistency, which proves that the method of the present invention has certain precision for estimating the electrical property of the antenna after covering, and the estimation result is more accurate.
In order to further prove the effect of the simulation experiment of the invention, the calculation results of the invention and the prior art are respectively evaluated by utilizing the change conditions (gain change, left-side lobe level change and right-side lobe level change) of three electrical performance indexes.
According to the calculation formula of the step 6 of the invention, the results of the electrical property index changes before and after the thickness and dielectric constant of the radome change are obtained, and all the calculation results are drawn as table 1:
TABLE 1 Electrical Performance analysis of the present invention and Prior Art estimation results in simulation experiments
It can be seen from table 1 that the error of gain change caused by the thickness change of the radome is estimated to be 0.06dB, the error of estimation of the level change of the left side lobe is 0.28dB, and the error of estimation of the level change of the right side lobe is 0.28dB, meanwhile, the error of gain change caused by the dielectric constant change of the radome is estimated to be 0.05dB, the error of estimation of the level change of the left side lobe is 0.18dB, and the error of estimation of the level change of the right side lobe is 0.17 dB.
The above simulation experiments show that: the method comprises the steps of respectively forming an antenna housing thickness matrix and a dielectric constant matrix by all variable quantities of antenna housing thickness and dielectric constant change caused by real-time measurement of natural environment erosion, substituting the variable quantities into an antenna housing electrical property analysis model to obtain electrical property after the antenna housing is changed, simplifying the step of analyzing the electrical property complexity after the antenna housing is changed, and judging whether the telecommunication energy meets the design requirements of the antenna electrical property when the thickness of the antenna housing and the dielectric constant are changed or not in real time, so that the problems that in the prior art, the physical property parameter real-time change of the antenna housing can only be analyzed by repeating the same calculation step for multiple times, a large amount of calculated quantity and calculation time are caused, and the electrical property when the thickness of the antenna housing and the dielectric constant are changed cannot be estimated in real time are solved, and the method is very practical and efficient.
Claims (7)
1. An antenna housing electrical property estimation method based on a Taylor expansion method is characterized in that: establishing an antenna housing electrical property analysis model based on a Taylor expansion method, and combining the antenna housing electrical property analysis model and a real-time monitoring result to analyze the electrical property of the antenna housing when the thickness and the dielectric constant change; the steps of the estimation method include the following:
step 1, generating a grid of an antenna housing model:
with the wavelength length of the electromagnetic wave of 1/5 as the size of each grid cell, randomly dividing the radome model into a plurality of triangular cell grids;
step 2, generating a far field intensity directional diagram of the antenna with the antenna housing:
(2a) respectively calculating transmission coefficients in the vertical polarization direction and the horizontal direction of the electromagnetic waves on each unit grid;
(2b) calculating the transmission coefficient of each unit grid along the main polarization direction of the electromagnetic waves;
(2c) calculating far field intensities of the antenna with the antenna housing in spatial observation directions with different pitch angles and azimuth angles, drawing the far field intensities in all the spatial observation directions into a far field directional diagram, and extracting electrical performance indexes of the antenna with the antenna housing;
step 3, generating an antenna housing thickness matrix and a dielectric constant matrix:
respectively forming an antenna housing thickness matrix and a dielectric constant matrix by using all variable quantities of antenna housing thickness and dielectric constant change caused by natural environment erosion measured in real time;
step 4, establishing an antenna housing electrical property analysis model:
(4a) the antenna housing electrical property analysis model when the thickness of the antenna housing changes is as follows:
wherein, Fθ′,φThe field intensity of the antenna with the radome after the nth column element in the thickness matrix is added in the far field in the space observation direction of the pitch angle theta azimuth angle phi,
it is shown that the operation of derivation is performed,
representing the field intensity of the antenna in the far field in the space observation direction of the pitch angle theta and the azimuth angle phi before the change of the antenna housing, h representing the dielectric constant of the antenna housing [. cndot.)]TDenotes a transpose operation, Δ Xh,nRepresenting the variation of the thickness of the antenna housing at the nth time;
(4b) the antenna housing electrical property analysis model when the antenna housing dielectric constant changes is as follows:
wherein, F ″)θ,φRepresenting the field intensity of the antenna with the radome after the addition of the nth row element in the dielectric constant matrix in the far field in the space observation direction of the pitch angle theta azimuth angle phi, wherein epsilon represents the dielectric constant of the radome, and delta Xε,nRepresenting the variation of the dielectric constant of the antenna housing at the nth time;
step 5, estimating the electrical property after the thickness and the dielectric constant of the antenna housing are changed:
substituting the radome thickness matrix and the dielectric constant matrix into a radome electrical property analysis model, respectively calculating the antenna far field intensity after the radome thickness and the dielectric constant are changed, drawing a far field directional diagram after the radome thickness and the dielectric constant are changed, and extracting electrical property indexes after the radome thickness and the dielectric constant are changed;
step 6, comparing the electrical performance indexes before and after the thickness and the dielectric constant of the antenna housing are changed;
step 7, judging whether the antenna electrical property after the thickness and the dielectric constant of the antenna housing are changed meets the design requirement of the antenna electrical property, if so, executing step 8, otherwise, executing step 3;
and 8, outputting far field directional diagrams and electrical performance indexes before and after the thickness and the dielectric constant of the antenna housing change.
2. The radome electrical property estimation method based on the taylor expansion method as claimed in claim 1, wherein the transmission coefficients in the vertical polarization direction and the horizontal direction of the electromagnetic wave on each unit cell in the step (2a) are calculated by the following formula:
wherein, TVmRepresents a transmission coefficient in the horizontal polarization direction of the electromagnetic wave, T, on the m-th unit cellHmDenotes a transmission coefficient in a direction perpendicular to the polarization direction of the electromagnetic wave on the m-th cell grid, ch (·) denotes a hyperbolic cosine operation, j denotes an imaginary unit, V denotes an intermediate variable, and V ═ 2 pi (epsilon-sin)2θm)1/2[ lambda ], [ pi ] denotesCircumference ratio, ε represents the dielectric constant of the radome, sin represents the sinusoidal operation, θmThe incident angle of the electromagnetic wave on the m-th cell grid, λ represents the wavelength of the frequency at which the antenna operates,
c denotes the speed of light, f denotes the frequency at which the antenna operates, Z1Representing the wave impedance in free space.
3. The radome electrical property estimation method based on the taylor expansion method as claimed in claim 1, wherein the transmission coefficient in the main polarization direction of the electromagnetic wave on each unit cell in the step (2b) is calculated by the following formula:
wherein, TmDenotes the transmission coefficient in the main polarization direction of the electromagnetic wave on the m-th grid, cos denotes the cosine operation, betamRepresenting the polarization angle of the incident electromagnetic wave on the m-th grid, exp representing the exponential operation with e as the base, ηHmRepresents the transmission coefficient T in the vertical polarization direction of the electromagnetic wave on the m-th gridHmThe phase of (a) is determined,
as indicated by the intermediate variable(s),
δ′m＝ηHm-ηVm，ηVmrepresents the transmission coefficient T in the horizontal polarization directionVmThe phase of (c).
4. The method for estimating the electrical property of the radome based on the taylor expansion method as claimed in claim 1, wherein the far field strength of the radome antenna after the radome is added in the step (2c) in the space observation directions with different pitch angles and azimuth angles is calculated by the following formula:
wherein S' represents the area of the aperture field of the antenna, E0Representing the field strength at the aperture field of the antenna, e(·)Denotes an exponential operation with a natural constant e as base, k0Which represents the propagation constant of the free space,
xmand ymRespectively representing the position coordinates of the m-th grid projected on the antenna aperture field.
5. The radome electrical property estimation method based on the taylor expansion method as claimed in claim 1, wherein the radome electrical property indicator in step (2c) refers to: the antenna gain, the left pair lobe level and the right pair lobe level are changed before the antenna housing.
6. The radome electrical property estimation method based on the taylor expansion method as claimed in claim 1, wherein the electrical property index after the radome thickness and the dielectric constant are changed in the step (5) refers to: the antenna comprises the gain, the left-side lobe level and the right-side lobe level of the antenna after the thickness of the antenna housing is changed, and the gain, the left-side lobe level and the right-side lobe level of the antenna after the dielectric constant of the antenna housing is changed.
7. The radome electrical property estimation method based on the taylor expansion method as claimed in claim 1, wherein the electrical property index before and after the radome thickness and dielectric constant change is compared in step (6) is obtained by the following formula:
wherein, Δ Gh,nRepresents the gain loss value, Δ G, of the antenna after the nth radome thickness changeε,nShows the gain loss value G of the antenna after the dielectric constant of the antenna housing changes for the nth time1Denotes the gain, Δ L, of the antenna before the change of the radomeh,nThe left side lobe level rise value L of the antenna after the nth antenna housing thickness change is represented1Indicating the left side lobe level, L, of the antenna before the radome changesh,nRepresents the left side lobe level, Δ L, of the antenna after the nth radome thickness changeε,nThe level rise value L of the left side lobe of the antenna after the dielectric constant of the antenna housing is changed for the nth timeε,nRepresents the left side lobe level, Δ R, of the antenna after the nth radome dielectric constant changeh,nRepresents the right minor lobe level rise value, R, of the antenna after the nth radome thickness change1Indicating the right side lobe level, R, of the antenna before the radome changeh,nRepresents the right minor lobe level, Δ R, of the antenna after the nth radome thickness changeε,nRepresents the right minor lobe level rise value, R, of the antenna after the dielectric constant of the antenna housing changes for the nth timeε,nThe right-lobe level of the antenna after the nth radome permittivity change is shown.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110865829.8A CN113625062B (en) | 2021-07-29 | 2021-07-29 | Antenna housing electrical property estimation method based on Taylor expansion method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110865829.8A CN113625062B (en) | 2021-07-29 | 2021-07-29 | Antenna housing electrical property estimation method based on Taylor expansion method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113625062A true CN113625062A (en) | 2021-11-09 |
| CN113625062B CN113625062B (en) | 2022-05-17 |
Family
ID=78381670
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110865829.8A Active CN113625062B (en) | 2021-07-29 | 2021-07-29 | Antenna housing electrical property estimation method based on Taylor expansion method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113625062B (en) |
Citations (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5359338A (en) * | 1989-09-20 | 1994-10-25 | The Boeing Company | Linear conformal antenna array for scanning near end-fire in one direction |
| US20030038745A1 (en) * | 2001-08-24 | 2003-02-27 | Farzin Lalezari | Antenna apparatus including compound curve antenna structure and feed array |
| JP2010177983A (en) * | 2009-01-29 | 2010-08-12 | Japan Radio Co Ltd | Patch array antenna |
| US20110192643A1 (en) * | 2008-09-12 | 2011-08-11 | Pilar Marin Palacios | Electromagnetic radiation attenuator and method for controlling the spectrum thereof |
| CN102590656A (en) * | 2012-01-03 | 2012-07-18 | 西安电子科技大学 | Antenna cover electric property forecasting method based on distant field |
| FR2970569A1 (en) * | 2011-01-17 | 2012-07-20 | Eads Europ Aeronautic Defence | Method for designing or repairing radome for antenna, involves repeating steps of determination of radioelectric performances and modification of structure of radome until modified structure is consistent with terms of use of radome |
| CN103401070A (en) * | 2013-07-13 | 2013-11-20 | 西安电子科技大学 | Far field-based wall thickness grinding method for thin shell-type radome |
| CN105977631A (en) * | 2016-06-03 | 2016-09-28 | 西安电子科技大学 | Composite material radome far field directional diagram interval analysis |
| CN106469850A (en) * | 2016-09-14 | 2017-03-01 | 西安电子科技大学 | A kind of Thickness Design Method of antenna house |
| CN107219410A (en) * | 2017-06-21 | 2017-09-29 | 西安空间无线电技术研究所 | A kind of Planar Near-Field Measurement modification method based on probe frequency sweep shift offset |
| CN107391880A (en) * | 2017-08-18 | 2017-11-24 | 上海无线电设备研究所 | A kind of accurate efficient radome electrical property energy analysis method |
| CN108984880A (en) * | 2018-07-04 | 2018-12-11 | 西安电子科技大学 | Anamorphic array antenna electric performance fast-compensation method under dynamic load based on least square method |
| CN109613346A (en) * | 2018-12-28 | 2019-04-12 | 中国科学院新疆天文台 | The reconstructing method and system of reflector antenna far-field pattern |
| CN110470916A (en) * | 2019-08-16 | 2019-11-19 | 西安电子科技大学 | Reflector antenna electrical property rapid evaluation and compensation method under solar temperature field action |
| CN110581368A (en) * | 2019-09-18 | 2019-12-17 | 湖南大学 | Hydrological monitoring radar flat microstrip array antenna and design method thereof |
| CN110658390A (en) * | 2018-06-28 | 2020-01-07 | Aptiv技术有限公司 | Method for determining alignment errors of an antenna and a vehicle using an antenna and a detection device |
| US20200049751A1 (en) * | 2017-04-21 | 2020-02-13 | The Board Of Regents Of The University Of Oklahoma | Apparatus and method(s) for wet radome characterization and radar calibration |
| CN111144048A (en) * | 2019-12-23 | 2020-05-12 | 北京机电工程研究所 | Design optimization method of wave-transparent, invisible and heat-proof radome |
| CN111177890A (en) * | 2019-12-11 | 2020-05-19 | 西安电子科技大学 | Radome electrical performance uncertainty analysis method based on polygonal interval |
| CN111276822A (en) * | 2020-02-24 | 2020-06-12 | 西安电子科技大学 | Antenna pattern main lobe controllable antenna array design method |
| CN211700553U (en) * | 2020-03-13 | 2020-10-16 | 上海几何伙伴智能驾驶有限公司 | Low sidelobe radar antenna and vehicle-mounted radar antenna |
| CN112840225A (en) * | 2018-10-12 | 2021-05-25 | 佩里森股份有限公司 | Antenna housing measuring system and method |
-
2021
- 2021-07-29 CN CN202110865829.8A patent/CN113625062B/en active Active
Patent Citations (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5359338A (en) * | 1989-09-20 | 1994-10-25 | The Boeing Company | Linear conformal antenna array for scanning near end-fire in one direction |
| US20030038745A1 (en) * | 2001-08-24 | 2003-02-27 | Farzin Lalezari | Antenna apparatus including compound curve antenna structure and feed array |
| US20110192643A1 (en) * | 2008-09-12 | 2011-08-11 | Pilar Marin Palacios | Electromagnetic radiation attenuator and method for controlling the spectrum thereof |
| JP2010177983A (en) * | 2009-01-29 | 2010-08-12 | Japan Radio Co Ltd | Patch array antenna |
| FR2970569A1 (en) * | 2011-01-17 | 2012-07-20 | Eads Europ Aeronautic Defence | Method for designing or repairing radome for antenna, involves repeating steps of determination of radioelectric performances and modification of structure of radome until modified structure is consistent with terms of use of radome |
| CN102590656A (en) * | 2012-01-03 | 2012-07-18 | 西安电子科技大学 | Antenna cover electric property forecasting method based on distant field |
| CN103401070A (en) * | 2013-07-13 | 2013-11-20 | 西安电子科技大学 | Far field-based wall thickness grinding method for thin shell-type radome |
| CN105977631A (en) * | 2016-06-03 | 2016-09-28 | 西安电子科技大学 | Composite material radome far field directional diagram interval analysis |
| CN106469850A (en) * | 2016-09-14 | 2017-03-01 | 西安电子科技大学 | A kind of Thickness Design Method of antenna house |
| US20200049751A1 (en) * | 2017-04-21 | 2020-02-13 | The Board Of Regents Of The University Of Oklahoma | Apparatus and method(s) for wet radome characterization and radar calibration |
| CN107219410A (en) * | 2017-06-21 | 2017-09-29 | 西安空间无线电技术研究所 | A kind of Planar Near-Field Measurement modification method based on probe frequency sweep shift offset |
| CN107391880A (en) * | 2017-08-18 | 2017-11-24 | 上海无线电设备研究所 | A kind of accurate efficient radome electrical property energy analysis method |
| CN110658390A (en) * | 2018-06-28 | 2020-01-07 | Aptiv技术有限公司 | Method for determining alignment errors of an antenna and a vehicle using an antenna and a detection device |
| CN108984880A (en) * | 2018-07-04 | 2018-12-11 | 西安电子科技大学 | Anamorphic array antenna electric performance fast-compensation method under dynamic load based on least square method |
| CN112840225A (en) * | 2018-10-12 | 2021-05-25 | 佩里森股份有限公司 | Antenna housing measuring system and method |
| CN109613346A (en) * | 2018-12-28 | 2019-04-12 | 中国科学院新疆天文台 | The reconstructing method and system of reflector antenna far-field pattern |
| CN110470916A (en) * | 2019-08-16 | 2019-11-19 | 西安电子科技大学 | Reflector antenna electrical property rapid evaluation and compensation method under solar temperature field action |
| CN110581368A (en) * | 2019-09-18 | 2019-12-17 | 湖南大学 | Hydrological monitoring radar flat microstrip array antenna and design method thereof |
| CN111177890A (en) * | 2019-12-11 | 2020-05-19 | 西安电子科技大学 | Radome electrical performance uncertainty analysis method based on polygonal interval |
| CN111144048A (en) * | 2019-12-23 | 2020-05-12 | 北京机电工程研究所 | Design optimization method of wave-transparent, invisible and heat-proof radome |
| CN111276822A (en) * | 2020-02-24 | 2020-06-12 | 西安电子科技大学 | Antenna pattern main lobe controllable antenna array design method |
| CN211700553U (en) * | 2020-03-13 | 2020-10-16 | 上海几何伙伴智能驾驶有限公司 | Low sidelobe radar antenna and vehicle-mounted radar antenna |
Non-Patent Citations (3)
| Title |
|---|
| WANYE XU等: "A New Efficient Thickness Profile Design Method for Streamlined Airborne Radomes", 《 IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》 * |
| 许万业等: "基于远场的薄壳式天线罩壁厚修磨量反求方法", 《系统工程与电子技术》 * |
| 陶海志等: "多层天线罩的电性能分析", 《现代电子技术》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113625062B (en) | 2022-05-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109472066B (en) | Reflecting surface antenna random error analysis method based on unit central point displacement | |
| CN113625062B (en) | Antenna housing electrical property estimation method based on Taylor expansion method | |
| CN113964549A (en) | Design method and device of space sampling antenna based on interference cancellation | |
| Le Coq et al. | IETR millimeter-wave Compact Antenna Test Range implementation and validation | |
| Park et al. | Analysis of Curved Frequency selective surface for radome using characteristic basis function method | |
| Mrdakovic et al. | A method for full wave analysis of electrically large transparent radomes | |
| Parini et al. | Computational electromagnetic modelling of compact antenna test range quiet zone probing: A comparison of simulation techniques | |
| Fuchs et al. | Compressive sensing approach for fast antenna far field characterization | |
| Granet et al. | Reflectarray compact antenna test range concept | |
| CN109301500A (en) | The design method of Chebyshev's micro-strip array antenna | |
| Rotgerink et al. | Modelling of effects of nose radomes on radar antenna performance | |
| Koshkid'ko et al. | Design and investigation of a linear equidistant slotted waveguide antenna | |
| Ma et al. | Low-cost lens antenna design for microwave moisture detection | |
| He et al. | Design of microstrip array antenna by using active element pattern technique combining with taylor synthesis method | |
| Huang et al. | An array of two periodic leaky-wave antennas with sum and difference beam scanning for application in target detection and tracking | |
| Tolin et al. | Uwb millimeter-wave 1d mimo array for non-destructive testing | |
| Boriskina et al. | Effect of the imperfect flat Earth on the vertically polarized radiation of a cylindrical reflector antenna | |
| Addamo et al. | Bloch wave analysis of long leaky coaxial cables | |
| Alves et al. | Radar cross section of simple and complex targets in the C-band: A comparison between anechoic chamber measurements and simulations | |
| Tian et al. | Examination of spherical antenna far-field scattering suppression through electromagnetic simulation | |
| Fuscaldo et al. | Radiating properties of 1-D bidirectional leaky-wave antennas | |
| Qiang et al. | Application of AI-SI technique to large airborne radome analysis | |
| La et al. | 3.5 kW wind turbine for cellular base station: Radar cross section modelling and measurement | |
| Shao et al. | Field Strength Prediction of 230 MHz Urban Macro Cellular Electric Wireless Private Network Based on Uniform Geometrical Theory of Diffraction | |
| Oğuzer et al. | Analysis of the nonconcentric reflector antenna-in-radome system by the iterative reflector antenna and radome interaction |
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 |