CN107783086A - Method for the distorted position of diagnostic antenna battle array bore width phase field - Google Patents
Method for the distorted position of diagnostic antenna battle array bore width phase field Download PDFInfo
- Publication number
- CN107783086A CN107783086A CN201710966768.8A CN201710966768A CN107783086A CN 107783086 A CN107783086 A CN 107783086A CN 201710966768 A CN201710966768 A CN 201710966768A CN 107783086 A CN107783086 A CN 107783086A
- Authority
- CN
- China
- Prior art keywords
- formula
- field
- spectrum
- antenna
- plane waves
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4004—Means for monitoring or calibrating of parts of a radar system
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/021—Auxiliary means for detecting or identifying radar signals or the like, e.g. radar jamming signals
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
A kind of method that the present invention provides distorted position for diagnostic antenna battle array bore width phase field, it is to carry out filed-close plane wave spectrum to Near-field Data to convert to obtain the spectrum of plane waves component of k-space first, then the directional diagram for implementing probe and unit in k-space spectrum of plane waves corrects, and finally carries out the inverse transformation of spectrum of plane waves aperture field and obtains Aperture field distribution.Finally, in order to obtain more accurate data, it is necessary to carry out the reconstruct of aperture field, the width distributed mutually of actual physical location on mouth face is obtained, so as to judge position and corresponding radiating element that aperture field or exciting current be distorted, reach the purpose that " diagnosis " is carried out to antenna.The present invention can be accurately positioned antenna failure cell position, investigate failure in time, save a large amount of testing times, and improve direction plot quality by way of inverting compensation.In addition, the present invention makees efficiently to calculate using FFT (FFT), thus there is stronger engineering practicability.
Description
Technical field
The present invention relates to the distortion of antenna near-field technical field of measurement and test, more particularly, to diagnostic antenna battle array bore width phase field
The method of position.
Background technology
In recent years, array antenna especially phased array antenna technology has obtained faster development, various dual-use
It is widely used in radio system.At present, the development task of domestic array antenna especially phased array antenna is very
Heavy, but due to reasons such as processing and manufacturings, antenna has a large amount of inevitable errors in actual design, array shakes in addition
The aging of coupling, active device between dynamic, antenna element and its change of temperature characterisitic etc. can all introduce unpredictable to antenna
Error or cause disabling unit, this just needs that the width of radar antenna is mutually diagnosed and debugged, to determine that antenna respectively radiates
The width distributed mutually situation of electric current on unit.
At present, mainly there are two methods the country when diagnosing, debugging large-scale array antenna：The first is close with probe
Antenna array, the region of the unit to be failed in battle array or unusual change is measured roughly.The method can not accurately measure Antenna aperture
Width distributed mutually, because stronger multipath effect between probe and AUT be present, this effect can be as probe be on scanning plane
Sample motion and change, so as to result in very big Error Diagnostics, and such a method needs constantly artificial mobile probe
With manually recorded data, substantial amounts of manpower and time cost can be expended.Second is embedding using the business software of NSI companies of the U.S.
The bore inverting module of set carries out the diagnosis of aperture field width phase, but external using bore inverting module as business secret, undisclosed
The principle of its inverting, mode and approach are had no way of understanding by its principle and source code, the external world, no calligraphy learning, references or even basic at it
Upper enhancement and optimization；Moreover, bore inverting module is mainly used on radar, it is related to national defense industry, completely using other country's software
The risk that radar parameter is divulged a secret be present；Finally, the bore inverting module is as software part encapsulation outlet, expensive, state
It is interior independently to apply.No matter from independent research, independent intellectual property right, the angle for ensuring radar test parameter secrecy, or from
Product secondary development, the angle for reducing production cost, shortening the manufacturing cycle, are required for proposing a set of suitable China's radar (antenna
Battle array) be actually needed, technical scheme that there is perfect knowledge property right, also the only source code write on this basis just have and pacify
Quan Xing, prevent back door.
In view of the above-mentioned problems, domestic there are a large amount of scientific research institutions to carry out beneficial research and probe and trial, mainly have at present
The method of equivalent magnetic current and Plane wave spectrum.But in practice, our the invention above methods still suffer from the limitation of use environment,
Specially：
The method of equivalent magnetic current, which is mainly disadvantageous in that, does not carry out probe orientation figure correction, so that this method moves towards real
With being restricted.In addition, in the method, due to calculator memory can not when being limited in and determining equivalent magnetic current by Near-field Data
Consider the contribution of equivalent magnetic current in whole plane, the tribute of equivalent magnetic current in antenna opening diametric plane to be measured or slightly larger region can only be considered
Offer, plus the influence of Finite scanning plane in Planar Near-Field Measurement so that the precision of the Aperture field distribution obtained using this method by
Limit.In addition, this method needs solution matrix equation when determining equivalent magnetic current by Near-field Data, therefore, when antenna aperture to be measured
Face and scanning plane are very big, and during so as to cause matrix element a lot, computational efficiency is not high, adds time cost.
And existing spectrum of plane waves rule is schemed with vertical direction in only the horizontal direction when being odd number point, can just make
With.And in actual antennas test, due to frequency, sampling face and probe to the factors such as antenna distance, sampling number is probably even
Several is also likely to be odd number, and this results in domestic traditional method using being above very restricted.
The content of the invention
In order to solve above-mentioned technological deficiency present in prior art, the present invention provides one kind and is used for diagnostic antenna battle array bore
The method of the distorted position of width phase field.
The present invention is achieved by the following technical solutions：
For the method for the distorted position of diagnostic antenna battle array bore width phase field, the antenna array includes the radiation of more than 4
Unit, the radiating element are generated and transmitted by exciting current respectively, and bore width phase field is collectively formed by the exciting current；The party
Method comprises the following steps：
Initial step：
Obtain the planar near-field test data of the antenna array；
Aligning step：
Near field-spectrum of plane waves conversion of probe orientation figure correction is carried out to the planar near-field test data, is carried out
The spectrum of plane waves component of k-space after probe correction, the spectrum of plane waves directional diagram of the k-space after being corrected；
Inversion step：
Spectrum of plane waves-aperture field inverse transformation is carried out to the spectrum of plane waves directional diagram of the k-space after the correction, obtains inverting
Bore width phase field pattern afterwards；
Compare step：
The actual physical location of each radiating element of the antenna array is measured, by the actual physics position of each radiating element
Put and be compared with the position in the bore width phase field pattern after the inverting, so as to draw the width of each radiating element of antenna array
Distributed mutually, and then judge the position that is distorted and corresponding radiating element.
The present invention is relative to the beneficial effect of prior art：
(1) present invention is the distribution by measuring antenna near field, and backstepping goes out distribution or each spoke of antenna of antenna aperture field
The exciting current of unit is penetrated, so as to judge position and corresponding radiation list that aperture field or exciting current be distorted
Member, reach the purpose that " diagnosis " is carried out to antenna, while take this as a foundation and each passage of antenna is corrected, make bore face width phase
Distribution meets design requirement, so that antenna performance hit the target requires.
(2) present invention can be accurately positioned antenna failure cell position, investigate failure in time, save a large amount of testing times,
And improve direction plot quality by way of inverting compensation.In addition, the present invention is made efficiently using FFT (FFT)
Calculate, thus there is stronger engineering practicability.
(3) present invention thoroughly overcomes China when diagnosing, debugging large-scale array antenna, mainly manually close with probe
Antenna array, the method for measuring the region of the unit to be failed in battle array or unusual change roughly.Overcoming can not be accurate using this method
The problem of measuring the width distributed mutually of Antenna aperture, is greatly reduced Error Diagnostics.
(4) present invention firstly provides the distorted position that diagnostic antenna battle array bore width phase field is carried out by the way of computer inverting
Method, break the external commercial monopoly to the high-accuracy technology in China.Using the technical program fidelity height, avoid being artificially introduced
Error, inverting efficiency high, use range it is wide, the software based on the inventive method, there is complete intellectual property, after being not present
The risk of door.
(5) it the software based on the inventive method, independently can on computers run, reduce purchase cost, and have and open again
Hair, again improved upgrading advantage.
(6) present invention realizes the inverting and diagnosis of antenna array width phase in planar near-field test, realizes the country near
Intellectual property autonomy-oriented in terms of the diagnosis of field, invention can be nested near field test system as standalone module, be that China is near
Building for test system of field lays a solid foundation.
(7) present invention realizes the inverting and diagnosis of antenna array width phase in planar near-field test, realizes the country near
Intellectual property autonomy-oriented in terms of the diagnosis of field, invention can be nested near field test system as standalone module, be that China is near
Building for test system of field lays a solid foundation.
(8) the present invention overcomes the deficiency of the method for equivalent magnetic current, existing Plane wave spectrum, it is contemplated that probe orientation figure
Influence, quick FFT is used in calculating process, and any kind of sampling number is applicable, thus with precision it is high,
Calculating speed is fast and the advantages that versatility.In other words, can solve array antenna aperture field width simultaneously by the present invention mutually to examine
The problems such as disconnected precision is not high, calculating speed is slow and algorithm is poor for applicability.
Brief description of the drawings
Fig. 1 is the general flow chart of the method for the distorted position that the present invention is used for diagnostic antenna battle array bore width phase field.
Fig. 2 is probe when being horizontal linear polarization, the corresponding relation figure of probe co-ordinate system and antenna coordinate system.
Fig. 3 isSchematic diagram of the vector in space coordinates.
When Fig. 4 positions probe is vertical linear polarization, the corresponding relation figure of probe co-ordinate system and antenna coordinate system.
Fig. 5 isSchematic diagram of the vector in space coordinates.
Fig. 6 is the 3 D stereo amplitude pattern using embodiments of the invention 1.
Fig. 7 is the aperture field horizontal amplitude distribution curve using embodiments of the invention 1.
Fig. 8 is the aperture field horizontal phase distribution curve using embodiments of the invention 1.
Fig. 9 is to navigate to each element amplitude distribution curve using embodiments of the invention 1.
Figure 10 is the 3 D stereo amplitude pattern using embodiments of the invention 2.
Figure 11 is using the aperture field horizontal amplitude distribution curve of embodiments of the invention 2.
Figure 12 is using the aperture field horizontal phase distribution curve of embodiments of the invention 2.
Figure 13 uses embodiments of the invention 2 to navigate to each element amplitude distribution curve.
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, it is right below in conjunction with drawings and Examples
The present invention is further elaborated.It should be appreciated that embodiment described herein is only to explain the present invention,
It is not intended to limit the present invention.
As shown in figure 1, the method for the distorted position for diagnostic antenna battle array bore width phase field, the antenna array includes 4
Radiating element above, the radiating element are generated and transmitted by exciting current, bore are collectively formed by the exciting current respectively
Width phase field；Characterized in that, this method comprises the following steps：
Step S1, initial step：
Obtain the planar near-field test data of the antenna array
Step S2, aligning step：
Near field-spectrum of plane waves conversion of probe orientation figure correction is carried out to the planar near-field test data, is carried out
The spectrum of plane waves component of k-space after probe correction, the spectrum of plane waves directional diagram of the k-space after being corrected.
Specifically, the aligning step specifically includes：
If main electric field polarization of popping one's head in is orientated, scanned in z=d planes in the x-direction, obtain when probe is horizontal polarization and adopt
Near-field Data and antenna wave spectrum between relation function：
Wherein, the implication of parameter is respectively that k is wave number；BxIt is for probe
The signal measured during horizontal polarization in z=d planes；F[Bx] it is BxTWO-DIMENSIONAL FOURIER inverse transformation；For the plane wave of antenna
Spectrum, and have For probe be horizontal polarization when spectrum of plane waves, have
Above-mentioned function is normalized, is converted to a11Ax+a12Ay=Ix, the formula is on spectrum of plane waves component
An equation；
If main electric field polarization of popping one's head in is orientated, scanned in z=d planes in the y-direction, obtain when probe is vertical polarization and adopt
Near-field Data and antenna wave spectrum between relation function：
Wherein, the implication of parameter is respectively ByThe signal that is measured in z=d planes when for probe being vertical polarization；F[By]
It is ByTWO-DIMENSIONAL FOURIER inverse transformation；For probe be vertical polarization when spectrum of plane waves, have
Above-mentioned function is normalized, is converted to a21Ax+a22Ay=Iy, the formula be on spectrum of plane waves component another
Equation；
2 equations of simultaneous spectrum of plane waves component, obtained equation group are as follows：
A can be solved by above formulaxAnd Ay, and then try to achieve Az, so as to obtain the spectrum of plane waves after progress probe orientation figure correctionTangential componentFurther obtainValue.
In other words, in step sl, main electric field of popping one's head in is only lower in one direction to be scanned, and is that cannot get spectrum of plane waves component,
Only two kinds of situations join together to form a two dimensional equation group, can just obtain spectrum of plane waves.
Furtherly, the aligning step further comprises：
Referring to Fig. 2, if main electric field polarization of popping one's head in is orientated in the x-direction, when being scanned in z=d planes, the signal that measures is Bx
(x, y, d), then it can be obtained by the coupling formula of planar near-field antenna measurement
In formula, k is wave number；kzFor components of the k in z directions；For the vector in k directions；For the spectrum of plane waves of antenna；For probe be horizontal polarization when spectrum of plane waves；F[Bx(x, y, d)] it is BxThe TWO-DIMENSIONAL FOURIER contravariant of (x, y, d)
Change, that is, have：
In formula, π=3.14159；kxAnd kyRespectively k x and y-component；
The relation of probe co-ordinate system and antenna coordinate system is：
Referring to Fig. 3,Vector is in the relation of space coordinates：
Here,ForUnit Vector；
And
It can then be obtained by formula 303, formula 304 and formula 305,
As available from the above equation：
Have：
Then in (formula 301) formulaFollowing form can be turned to：
In formula, A 'θ′For The component in direction,For The component in direction；Formula 309 is substituted into
Formula 301, and according to formula 303, can obtain
Because during test, antenna is located at passive region, then has：Again because electric field and ripple
The relational expression of spectrum is：
In formula, a0Represent the complex amplitude of incidence wave；
Formula 312 is substituted into formula 311 to obtain
I.e.
I.e.
Then it can obtain
WillSubstitution formula 310, then have：
A ' can be obtained by formula 309x′, A 'y′, A 'z′With A 'θ′,Relation be：
A′z′=-A 'θ′Sin θ ' (formula 320)
Because the far-field pattern of probe can be approached by following formula：
And it is known that the relation between normalization far-field pattern function and spectrum of plane waves is：
So：
Above in two formulas,For far-field patternIn θ ' component；For far-field pattern Component, fE(θ ') and fH(θ ') is respectively the E faces and H faces normalization far-field pattern function popped one's head in；Then have
Two formulas formula 318, formula 319, formula 320 will be substituted into above, then had
Order
Then formula 317 can turn to following form：
a11Ax+a12Ay=Ix(formula 333)
Above formula is an equation on spectrum of plane waves component, and equation has two unknown quantity AxAnd Ay, it is necessary to pop one's head in
Rotation, which once measures, to be resettled an equation simultaneous solution and can obtain AxWith Ay。
Furtherly, the aligning step further comprises：
The signal that if main electric field polarization of popping one's head in is orientated in the y-direction, is measured when being scanned in z=d planes is By(x, y, d),
It can then be obtained by the coupling formula of planar near-field antenna measurement：
In formula,For probe be vertical polarization when spectrum of plane waves；It is ByThe two dimension of (x, y, d)
Fourier inverse transformations, that is, have
Referring to Fig. 4, in the situation, i.e., formula 401 it is corresponding in the case of, the relation of probe co-ordinate system and antenna coordinate system
For：
Referring to Fig. 5,Vector is under space coordinates：
Then
And under current coordinate system,
It can then be obtained by formula 403, formula 404 and formula 405,
As available from the above equation
Have
Then in formula 401Following form can be turned to：
In formula, A "θ″For The component in direction,For The component in direction；By the generation of formula 409
Enter formula 401, and according to formula 403, can obtain
WillAbove formula is substituted into, then is had：
A " can be obtained by formula 409x″, A "y″, A "z″With A "θ″,Relation be
A″z″=-A "θ″Sin θ " (formula 414)
Again because the far-field pattern of probe can be approached by following formula：
And it is known that the relation between normalization far-field pattern function and spectrum of plane waves is：
So
Above in two formulas,For far-field patternIn θ " component；For far field side
Xiang Tu Component, fE(θ ") and fH(θ ") is respectively the E faces and H faces normalization far-field pattern letter popped one's head in
Number；So have
Two formulas formula 412, formula 413,414 will be substituted into above, then had
Order
Then formula 411 can turn to following form：a21Ax+a22Ay=Iy(formula 427).
Simultaneous formula 333 and formula 427, obtained equation group are as follows：
A can be solved by above formulaxAnd Ay, and then try to achieve Az, so as to obtain the spectrum of plane waves after progress probe orientation figure correctionTangential componentFurther obtainValue.
Step S3, inversion step：
Spectrum of plane waves-aperture field inverse transformation is carried out to the spectrum of plane waves directional diagram of the k-space after the correction, obtains inverting
Bore width phase field pattern afterwards.
Specifically, the inversion step specifically includes：
Obtain carrying out the spectrum of plane waves after probe orientation figure correctionTangential componentCalculation Plane is near afterwards
The tangential component of field, and the value of its each unit opening position on array antenna bore face is extracted, specific implementation is：
The plane wave expansion formula of field is produced by antenna, the tangential component that can obtain electric field in z=h planes is
By parameterBy Fourier transform, the tangential component of electric field is converted intoFurtherly,
By parameterIt is converted into by Fourier transformConcretely comprise the following steps：
In formula, Δ x is x directions sampling interval, and Δ y is y directions sampling interval, and N is x directions sampling number, and M is y directions
Sampling number,
Therefore the tangential component of electric field is at each uniform lattice point in z=h planes
In formula, p and q are integer, are numeral corresponding to each sampled point；
If taking h=0, that is, the tangential electric field on array antenna bore face is can obtain, that is, obtain the bore width phase after inverting
Field distribution
Step S4, compare step.
The actual physical location of each radiating element of the antenna array is measured, by the actual physics position of each radiating element
Put and be compared with the position in the bore width phase field pattern after the inverting, so as to draw the width of each radiating element of antenna array
Distributed mutually, and then judge the position that is distorted and corresponding radiating element.
Obtained after inverting be 2 groups on position and the three-dimensional data of amplitude-phase value, one of which is range value, another
Group is phase value, and by interpolation method, the position that radiating element is inserted in above-mentioned three-dimensional data can obtain each radiating element
Amplitude-phase value.
Specifically included specifically, comparing step：
The actual physical location of each radiating element of antenna array is measured, to what is obtained in the inversion stepEnter row interpolation, you can draw the width distributed mutually of each radiating element of antenna array, and then judge what is be distorted
Position and corresponding radiating element.
Embodiment 1
For the accuracy that diagnostic method provided by the invention diagnoses to front disabling unit, certain L-band is hung down in darkroom
Linear polarization antenna is tested.Tested front horizontal cell number is 42, and horizontal cell spacing is 0.138m, vertical direction list
First number is 1, test frequency 990Mhz.Disabling unit is introduced simultaneously at the 36th and the 40th cell position, test and sharp
Inverting diagnosis is carried out with programmed algorithm, obtained result is with being analyzed as follows：
Fig. 6 is 3 D stereo amplitude pattern, and Fig. 7 and Fig. 8 are aperture field horizontal amplitude scatter chart and bore respectively
Field horizontal phase scatter chart.Fig. 9 is to navigate to each element amplitude scatter chart.
As can be seen from Figures 7 and 8, after introducing disabling unit, corresponding width distributed mutually is changed, in relevant position
Range value is decreased obviously.And the position of disabling unit can be gone out using Accurate Diagnosis from auditory localization cues Fig. 9 as the 36th and the 40th list
Member, so as to demonstrate the correctness that algorithm diagnoses to disabling unit.
Embodiment 2
Disabling unit diagnostic test, test frequency 3140MHz are carried out to certain S-band array antenna, horizontal cell number is
128, unit spacing is 0.046m.Disabling unit is introduced directly at the 5th and the 100th cell position, carries out the active shape near field
State is tested.Inverting diagnostic result is carried out with being analyzed as follows to test result using the present invention：
Figure 10, Figure 11 and Figure 12 are aperture field amplitude distribution figure, aperture field horizontal amplitude scatter chart, aperture field respectively
Horizontal phase scatter chart.Figure 13 is to navigate to each element amplitude scatter chart.
Figure 10 to Figure 12 has reacted the Aperture field distribution after introducing disabling unit, can determine the position of disabling unit substantially
Put.The position that auditory localization cues in Figure 13 are then accurately diagnosed to be disabling unit is the 5th and the 100th unit, so as to further
Demonstrate the accuracy that the algorithm of the present invention diagnoses to disabling unit.
Relative and prior art, compared with prior art, beneficial effect is as shown in table 1 by the present invention：
The beneficial effect of the invention relative to prior art of table 1
As it will be easily appreciated by one skilled in the art that these are only presently preferred embodiments of the present invention, not limiting
The present invention, all any modification, equivalent and improvement made within the spirit and principles of the invention etc., it should be included in this
Within the protection domain of invention.
Claims (6)
1. the method for the distorted position for diagnostic antenna battle array bore width phase field, the antenna array include the radiation list of more than 4
Member, the radiating element are generated and transmitted by exciting current respectively, and bore width phase field is collectively formed by the exciting current；Its feature
It is, this method comprises the following steps：
Initial step：
Obtain the planar near-field test data of the antenna array；
Aligning step：
Near field-spectrum of plane waves conversion of probe orientation figure correction is carried out to the planar near-field test data, obtains being popped one's head in
The spectrum of plane waves component of k-space after correction, the spectrum of plane waves directional diagram of the k-space after being corrected；
Inversion step：
Spectrum of plane waves-aperture field inverse transformation is carried out to the spectrum of plane waves directional diagram of the k-space after the correction, after obtaining inverting
Bore width phase field pattern；
Compare step：
Measure the actual physical location of each radiating element of the antenna array, by the actual physical location of each radiating element with
The position in bore width phase field pattern after the inverting is compared, so as to draw the width phase point of each radiating element of antenna array
Cloth, and then judge the position that is distorted and corresponding radiating element.
2. the method for the distorted position according to claim 1 for diagnostic antenna battle array bore width phase field, it is characterised in that
The aligning step specifically includes：
If main electric field polarization of popping one's head in is orientated, scanned in z=d planes in the x-direction, acquisition probe is adopted near when being horizontal polarization
The function of relation between field data and antenna wave spectrum：
Wherein, the implication of parameter is respectively that k is wave number；BxIt is horizontal for probe
The signal measured during polarization in z=d planes；F[Bx] it is BxTWO-DIMENSIONAL FOURIER inverse transformation；For the spectrum of plane waves of antenna,
And have For probe be horizontal polarization when spectrum of plane waves, have
Above-mentioned function is normalized, is converted to a11Ax+a12Ay=Ix, the formula is one on spectrum of plane waves component
Individual equation；
If main electric field polarization of popping one's head in is orientated, scanned in z=d planes in the y-direction, acquisition probe is adopted near when being vertical polarization
The function of relation between field data and antenna wave spectrum：
Wherein, the implication of parameter is respectively ByThe signal that is measured in z=d planes when for probe being vertical polarization；F[By] it is By
TWO-DIMENSIONAL FOURIER inverse transformation；For probe be vertical polarization when spectrum of plane waves, have
Above-mentioned function is normalized, is converted to a21Ax+a22Ay=Iy, the formula be on spectrum of plane waves component another
Equation；
2 equations of simultaneous spectrum of plane waves component, obtained equation group are as follows：
A can be solved by above formulaxAnd Ay, and then try to achieve Az, so as to obtain the spectrum of plane waves after progress probe orientation figure correction
Tangential componentFurther obtainValue.
3. the method for the distorted position according to claim 1 or 2 for diagnostic antenna battle array bore width phase field, its feature exist
In the aligning step further comprises：
The signal that if main electric field polarization of popping one's head in is orientated in the x-direction, is measured when being scanned in z=d planes is Bx(x, y, d), then by putting down
The coupling formula of face near field antenna measurements can obtain
In formula, k is wave number；kzFor components of the k in z directions；For the vector in k directions；For the spectrum of plane waves of antenna；For probe be horizontal polarization when spectrum of plane waves；F[Bx(x, y, d)] it is BxThe TWO-DIMENSIONAL FOURIER contravariant of (x, y, d)
Change, that is, have：
In formula, π=3.14159；kxAnd kyRespectively k x and y-component；
The relation of probe co-ordinate system and antenna coordinate system is：
Vector is in the relation of space coordinates：
Here,ForUnit Vector；
And
It can then be obtained by formula 303, formula 304 and formula 305,
As available from the above equation：
Have：
Then in (formula 301) formulaFollowing form can be turned to：
In formula, A 'θ′For The component in direction,For The component in direction；
Formula 309 is substituted into formula 301, and according to formula 303, can be obtained
Because during test, antenna is located at passive region, then has：
Again because the relational expression of electric field and wave spectrum is：
In formula, a0Represent the complex amplitude of incidence wave；
Formula 312 is substituted into formula 311 to obtain
I.e.
I.e.
Then it can obtain
WillSubstitution formula 310, then have：
A ' can be obtained by formula 309x′, A 'y′, A 'z′With A 'θ′,Relation be：
A′z′=-A 'θ′Sin θ ' (formula 320)
Because the far-field pattern of probe can be approached by following formula：
And it is known that the relation between normalization far-field pattern function and spectrum of plane waves is：
So：
Above in two formulas,For far-field patternIn θ ' component；For far-field pattern Component, fE(θ ') and fH(θ ') is respectively the E faces and H faces normalization far-field pattern function popped one's head in；
Then have
Two formulas formula 318, formula 319, formula 320 will be substituted into above, then had
Order
Then formula 317 can turn to following form：
a11Ax+a12Ay=Ix(formula 333)
Above formula is an equation on spectrum of plane waves component.
4. the method for the distorted position according to claim 1 or 2 for diagnostic antenna battle array bore width phase field, its feature exist
In the aligning step further comprises：
The signal that if main electric field polarization of popping one's head in is orientated in the y-direction, is measured when being scanned in z=d planes is By(x, y, d), then by putting down
The coupling formula of face near field antenna measurements can obtain：
In formula,For probe be vertical polarization when spectrum of plane waves；It is ByThe two dimension of (x, y, d)
Fourier inverse transformations, that is, have
In the situation, i.e., formula 401 it is corresponding in the case of, the relation of probe co-ordinate system and antenna coordinate system is：
Vector is under space coordinates：
Then
And under current coordinate system,
It can then be obtained by formula 403, formula 404 and formula 405,
As available from the above equation
Have
Then in formula 401Following form can be turned to：
In formula, A "θ″For The component in direction,For The component in direction；
Formula 409 is substituted into formula 401, and according to formula 403, can be obtained
WillAbove formula is substituted into, then is had：
A " can be obtained by formula 409x″, A "y″, A "z″With A "θ″,Relation be
A″z″=-A "θ″Sin θ " (formula 414)
Again because the far-field pattern of probe can be approached by following formula：
And it is known that the relation between normalization far-field pattern function and spectrum of plane waves is：
So
Above in two formulas,For far-field patternIn θ " component；For far-field pattern Component, fE(θ ") and fH(θ ") is respectively the E faces and H faces normalization far-field pattern function popped one's head in；Institute
To have
Two formulas formula 412, formula 413,414 will be substituted into above, then had
Order
Then formula 411 can turn to following form：a21Ax+a22Ay=Iy(formula 427).
5. the method for the distorted position according to claim 1 or 2 for diagnostic antenna battle array bore width phase field, its feature exist
In the inversion step specifically includes：
Obtain carrying out the spectrum of plane waves after probe orientation figure correctionTangential componentCalculation Plane near field afterwards
Tangential component, and the value of its each unit opening position on array antenna bore face is extracted, specific implementation is：
The plane wave expansion formula of field is produced by antenna, the tangential component that can obtain electric field in z=h planes is
By parameterBy Fourier transform, the tangential component of electric field is converted intoFurtherly, will join
NumberIt is converted into by Fourier transformConcretely comprise the following steps：
In formula, Δ x is x directions sampling interval, and Δ y is y directions sampling interval, and N is x directions sampling number, and M samples for y directions
Points,
Therefore the tangential component of electric field is at each uniform lattice point in z=h planes
In formula, p and q are integer, are numeral corresponding to each sampled point；
If taking h=0, that is, the tangential electric field on array antenna bore face is can obtain, that is, obtain the bore width phase field after inverting point
Cloth
6. being used for the method for the distorted position of diagnostic antenna battle array bore width phase field according to claim 1 or 5, its feature exists
In the comparison step specifically includes：
The actual physical location of each radiating element of antenna array is measured, to what is obtained in the inversion step
Enter row interpolation, you can draw the width distributed mutually of each radiating element of antenna array, and then judge that the position that is distorted and institute are right
The radiating element answered.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710966768.8A CN107783086B (en) | 2017-10-17 | 2017-10-17 | Method for diagnosing distorted position of antenna array aperture amplitude phase field |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710966768.8A CN107783086B (en) | 2017-10-17 | 2017-10-17 | Method for diagnosing distorted position of antenna array aperture amplitude phase field |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107783086A true CN107783086A (en) | 2018-03-09 |
| CN107783086B CN107783086B (en) | 2021-01-05 |
Family
ID=61433766
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710966768.8A Active CN107783086B (en) | 2017-10-17 | 2017-10-17 | Method for diagnosing distorted position of antenna array aperture amplitude phase field |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107783086B (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109521374A (en) * | 2018-10-18 | 2019-03-26 | 北京航空航天大学 | It is embedded in integrated balun and loads the tangential electric field probe of resonance of modified dipole |
| CN111044802A (en) * | 2019-12-04 | 2020-04-21 | 上海霍莱沃电子系统技术股份有限公司 | Phased array antenna unit amplitude and phase balancing method and device based on aperture field |
| CN111812418A (en) * | 2020-07-17 | 2020-10-23 | 中国汽车工程研究院股份有限公司 | Tire pressure monitoring antenna performance testing system and method |
| CN112385156A (en) * | 2018-07-16 | 2021-02-19 | 莱特普茵特公司 | System and method for over-the-air (OTA) testing to detect faulty elements in an active array antenna of an Extremely High Frequency (EHF) wireless communication device |
| CN113589051A (en) * | 2021-08-04 | 2021-11-02 | 西安电子科技大学 | Clutter suppression field electromagnetic measurement device and method |
| US20220404462A1 (en) * | 2021-06-22 | 2022-12-22 | Src, Inc. | Method for calibrating a phased array |
| US11982761B2 (en) * | 2021-06-22 | 2024-05-14 | Src, Inc. | Method for calibrating a phased array |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5343204A (en) * | 1993-07-29 | 1994-08-30 | Unisys Corporation | Auto-focusing correction for rotational acceleration effects on inverse synthetic aperture radar images |
| US5726656A (en) * | 1996-12-19 | 1998-03-10 | Hughes Electronics | Atmospheric correction method for interferometric synthetic array radar systems operating at long range |
| GB2318941A (en) * | 1996-10-07 | 1998-05-06 | Advantest Corp | Radio hologram image reconstruction using fast fourier transforms |
| US5809063A (en) * | 1996-10-25 | 1998-09-15 | General Electric Company | Coherent detection architecture for remote calibration of coherent systems using direct sequence spread spectrum transmission of reference and calibration signals |
| CN1664611A (en) * | 2005-03-24 | 2005-09-07 | 武汉大学 | Method for correcting passive channels based on non-linear antenna array |
| CN103411528A (en) * | 2013-08-26 | 2013-11-27 | 中国科学院空间科学与应用研究中心 | Method for calculating electric field probe rotation offset through circular polarization antenna axial ratio directional diagram |
| CN103424066A (en) * | 2013-08-26 | 2013-12-04 | 中国科学院空间科学与应用研究中心 | Method for using circularly polarized antenna far field phase difference for calculating rotation offset of electric field probe |
| CN103698753A (en) * | 2013-12-19 | 2014-04-02 | 武汉大学 | Passive passage correcting method of small-size array |
| CN104063544A (en) * | 2014-06-25 | 2014-09-24 | 电子科技大学 | Near-field effect error analysis method |
| CN104765022A (en) * | 2015-03-18 | 2015-07-08 | 中船重工鹏力（南京）大气海洋信息系统有限公司 | Methods for characteristic value probability statistics model establishment and antenna self-check based on echo spectra |
| CN105717493A (en) * | 2016-03-17 | 2016-06-29 | 哈尔滨工程大学 | Passive radar radiation source polarization identification and direction-finding system based on synthesis polarization method |
| CN106169654A (en) * | 2016-06-08 | 2016-11-30 | 中国电子科技集团公司第三十八研究所 | A kind of broadband active multibeam antenna system and bearing calibration thereof |
| CN106526562A (en) * | 2016-10-17 | 2017-03-22 | 西北工业大学 | Antenna directional pattern-based method for correcting RCS of calculation target |
| CN106546828A (en) * | 2016-10-31 | 2017-03-29 | 西北工业大学 | A kind of sphere Multi probe antenna measurement data processing method based on microwave network analysis |
-
2017
- 2017-10-17 CN CN201710966768.8A patent/CN107783086B/en active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5343204A (en) * | 1993-07-29 | 1994-08-30 | Unisys Corporation | Auto-focusing correction for rotational acceleration effects on inverse synthetic aperture radar images |
| GB2318941A (en) * | 1996-10-07 | 1998-05-06 | Advantest Corp | Radio hologram image reconstruction using fast fourier transforms |
| US5809063A (en) * | 1996-10-25 | 1998-09-15 | General Electric Company | Coherent detection architecture for remote calibration of coherent systems using direct sequence spread spectrum transmission of reference and calibration signals |
| US5726656A (en) * | 1996-12-19 | 1998-03-10 | Hughes Electronics | Atmospheric correction method for interferometric synthetic array radar systems operating at long range |
| CN1664611A (en) * | 2005-03-24 | 2005-09-07 | 武汉大学 | Method for correcting passive channels based on non-linear antenna array |
| CN103424066A (en) * | 2013-08-26 | 2013-12-04 | 中国科学院空间科学与应用研究中心 | Method for using circularly polarized antenna far field phase difference for calculating rotation offset of electric field probe |
| CN103411528A (en) * | 2013-08-26 | 2013-11-27 | 中国科学院空间科学与应用研究中心 | Method for calculating electric field probe rotation offset through circular polarization antenna axial ratio directional diagram |
| CN103698753A (en) * | 2013-12-19 | 2014-04-02 | 武汉大学 | Passive passage correcting method of small-size array |
| CN104063544A (en) * | 2014-06-25 | 2014-09-24 | 电子科技大学 | Near-field effect error analysis method |
| CN104765022A (en) * | 2015-03-18 | 2015-07-08 | 中船重工鹏力（南京）大气海洋信息系统有限公司 | Methods for characteristic value probability statistics model establishment and antenna self-check based on echo spectra |
| CN105717493A (en) * | 2016-03-17 | 2016-06-29 | 哈尔滨工程大学 | Passive radar radiation source polarization identification and direction-finding system based on synthesis polarization method |
| CN106169654A (en) * | 2016-06-08 | 2016-11-30 | 中国电子科技集团公司第三十八研究所 | A kind of broadband active multibeam antenna system and bearing calibration thereof |
| CN106526562A (en) * | 2016-10-17 | 2017-03-22 | 西北工业大学 | Antenna directional pattern-based method for correcting RCS of calculation target |
| CN106546828A (en) * | 2016-10-31 | 2017-03-29 | 西北工业大学 | A kind of sphere Multi probe antenna measurement data processing method based on microwave network analysis |
Non-Patent Citations (7)
| Title |
|---|
| J.J. LEE ; E.M. FERREN ; D.P. WOOLLEN ; K.M. LEE: "《Near-field probe used as a diagnostic tool to locate defective elements in an array antenna》", 《 IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》 * |
| R.G. YACCARINO ; Y. RAHMAT-SAMII: "《Phaseless bi-polar planar near-field measurements and diagnostics of array antennas》", 《 IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》 * |
| 潘宇虎: "《平面近场测量和诊断技术的研究》", 《中国优秀博硕士学位论文全文数据库 (硕士) 信息科技辑》 * |
| 焦禹: "《相控阵天线平面近场误差分析与校准方法研究》", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 王彬虎: "《平面近场测量技术及误差分析》", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 陈玉林、李刚: "《带探头补偿近场-口径场的反演与诊断》", 《雷达科学与技术》 * |
| 陈玉林; 于丁: "《电子科技》", 《平面近场测量中截断误差对口径场的不确定性分析》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112385156A (en) * | 2018-07-16 | 2021-02-19 | 莱特普茵特公司 | System and method for over-the-air (OTA) testing to detect faulty elements in an active array antenna of an Extremely High Frequency (EHF) wireless communication device |
| CN112385156B (en) * | 2018-07-16 | 2023-05-16 | 莱特普茵特公司 | System and method for detecting over-the-air (OTA) testing of failed elements in an active array antenna of an Extremely High Frequency (EHF) wireless communication device |
| CN109521374A (en) * | 2018-10-18 | 2019-03-26 | 北京航空航天大学 | It is embedded in integrated balun and loads the tangential electric field probe of resonance of modified dipole |
| CN109521374B (en) * | 2018-10-18 | 2020-10-27 | 北京航空航天大学 | Resonance tangential electric field probe embedded with integrated balun and loaded improved dipole |
| CN111044802A (en) * | 2019-12-04 | 2020-04-21 | 上海霍莱沃电子系统技术股份有限公司 | Phased array antenna unit amplitude and phase balancing method and device based on aperture field |
| CN111812418A (en) * | 2020-07-17 | 2020-10-23 | 中国汽车工程研究院股份有限公司 | Tire pressure monitoring antenna performance testing system and method |
| CN111812418B (en) * | 2020-07-17 | 2021-09-28 | 中国汽车工程研究院股份有限公司 | Tire pressure monitoring antenna performance testing system and method |
| US20220404462A1 (en) * | 2021-06-22 | 2022-12-22 | Src, Inc. | Method for calibrating a phased array |
| US11982761B2 (en) * | 2021-06-22 | 2024-05-14 | Src, Inc. | Method for calibrating a phased array |
| CN113589051A (en) * | 2021-08-04 | 2021-11-02 | 西安电子科技大学 | Clutter suppression field electromagnetic measurement device and method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107783086B (en) | 2021-01-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107783086A (en) | Method for the distorted position of diagnostic antenna battle array bore width phase field | |
| US11422176B2 (en) | Method and device for measuring radiation pattern of antenna array, and computer readable storage medium | |
| CN109581078B (en) | Directional diagram measuring system and method suitable for antenna in half-space environment | |
| CN104063544B (en) | Near-field effect error analysis method | |
| CN107783092A (en) | To rcs measurement system and method behind near field based on Chain relation | |
| Zhang et al. | An efficient probe calibration based near-field-to-near-field transformation for EMI diagnosis | |
| CN111366793A (en) | Planar antenna near field measurement method and system for reducing truncation error | |
| CN113804985A (en) | Anti-interference antenna directional pattern measuring method based on hybrid shielding chamber | |
| CN116520035A (en) | Two-dimensional pattern testing method adopting plane near-field rapid inversion method | |
| CN109286069A (en) | Array antenna field source reconstructing method and device based on dipole | |
| Varela et al. | Fast spherical near-field to far-field transformation for offset-mounted antenna measurements | |
| CN111505394A (en) | Antenna spherical surface near field measurement method based on probe antenna position error correction | |
| Khashimov et al. | Dynamic patterns technique for circular phased array diagnostics | |
| CN111487474B (en) | Numerical twin electromagnetic measuring system | |
| US11901636B2 (en) | Compact antenna test range (CATR) alignment verification | |
| CN113821907B (en) | Amplitude and phase automatic calibration method for large planar antenna array system | |
| CN112255469B (en) | Antenna arbitrary curved surface near field test phase acquisition device and method without reference signal | |
| CN111965602B (en) | Phased array radar amplitude-phase consistency detection method and system | |
| Varela et al. | Single-cut phaseless near-field measurements using specialized probes | |
| Cornelius et al. | Spherical near-field far-field transformation with infinite ground plane at arbitrary height z | |
| Manohar et al. | Mimicking antenna near-field measurements using full wave solvers for error characterization | |
| Wang et al. | Correction for probe-positioning errors in digital holographic reconstruction phase-less near-field measurement | |
| Zhengguang et al. | An efficient lineal sampling method for RCS prediction | |
| CN115833972B (en) | System and method for calibrating target angular position output by array feed system | |
| Kaslis et al. | Simulating $\theta $-zero and Axes Intersection Errors in Spherical Near-Field Antenna Measurements |
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 |