CN102708257B - Method for fast determining tolerance of active phased array antenna structure - Google Patents

Abstract

The invention belongs to the technical field of radar antennas, and specifically relates to a method for fast determining tolerance of an active phased array antenna structure. The method comprises the steps of: 1, determining information of an antenna array surface radiation unit; 2, obtaining a height error of the array surface radiation unit; 3, obtaining a position mounting error of the array surface radiation unit; 4, calculating a phase error of an antenna port; 5, calculating electric field distribution of a far field of the antenna; 6, calculating an electric performance parameter related to the antenna; 7, judging and calculating whether the electric performance of the antenna meets the requirements under the current structure tolerance condition according to the requirements of tactical and technical targets, if not, carrying out the step 8 and the step 9, and if so, directly turning to the step 10, wherein the tolerance of the position and the height of the current radiation unit is the fast determined tolerance of the antenna structure. The problem of difficulty in fast determining the tolerance on the distribution antenna structure is effectively solved.

Description

A kind of fast determination method of active phase array antenna structure tolerance

Technical field

The invention belongs to Radar Antenna System field, be specifically related to a kind of fast determination method of active phase array antenna structure tolerance, can be used for instructing determining fast and the evaluation of organization plan of active phase array antenna structure tolerance.

Background technology

Active phase array antenna is in recent years just in the new technology of accelerated development, can meet antenna high-performance, failure-survival capability requirement simultaneously, is also the gordian technique that meets modern radar growth requirement.China develops and formally equipped in 2007 alert-2000 early warning planes of sky of China air force voluntarily as the No.1 military priority project of country second to none of air force, and the breakthrough of its most critical is exactly to have equipped Active Phased Array Radar.Along with radar is day by day urgent to increasing the demands such as operating distance, raising Stealth Fighter, active phase array antenna electrical property is just suggested higher requirement.The state-of-the-art fourth generation fighter that China is developing at present---destroy 20 for meeting the tactical and technical norms such as large detection range, high stealth, also will equip high performance Active Phased Array Radar.

And active phase array antenna electrical property not only depends on the amplitude-phase error of feed system, be more vulnerable to the site error impact of structure radiating element, and this is limited by structural design.As GBR ground radar in American National missile defense systems (NMD) has adopted large-scale flat surface active phased array antenna physique just; its radar weighs 1000 tons; and U.S. SBX sea base radar is as the current unique platform that can bring into play multiple important detecting function in ballistic missile defense system; its radar weighs 2000 tons; for supporting so heavy radar equipment; and meet radar antenna electrical performance indexes, and guarantee antenna array precision, Antenna Construction Design has been proposed to extremely strict requirement.

Meanwhile, because the athletic meeting of antenna mounting platform brings vibrational excitation, airborne radar for example, thereby cause, active phase array antenna is lower in working order there is malformation, causes structure radiating element error.Known, the front malformation that the manufacture of active phase array antenna structure, alignment error and environmental load cause, all will cause the position of distance between antenna array radiating element and each radiating element that relative variation occurs, thereby cause that radiating element exciting current and aperture field PHASE DISTRIBUTION etc. change, so that cause that antenna gain declines, minor level raises and beam position inaccurate etc., finally can reduce the electrical property of active phase array antenna.For this reason, how according to antenna electric performance index request, determine structure tolerance, carry out structural design, and evaluation structure scheme, be a difficult problem that will inevitably run in the development of high-performance active phase array antenna.

For solving this difficult problem, common way is the installation the way of restraint according to active phase array antenna front, suppose that the known distortion of certain rule occurs antenna array structure, typical form of distortion has flexural deformation, bowl-shape distortion, saddle facial disfigurement etc., and then according to the situation of change of antenna electric performance, provide the maximum deformation quantity of front, and then during definite structure tolerance, concrete visible Analysis of performance of active phased array antennas with distorted plane error.International Journal of Electronics, 2009, 96 volumes, 5 phases, 549-559..This method exists main problem to be, constrained when front require to be installed or in front center line, or front surrounding, or four jiaos of fronts, otherwise front deformation rule is difficult to adopt clear and definite mathematical formulae to be described.This causes its range of application to be very limited, in engineering, is difficult to promote, because antenna deformation form is unknown in reality, be not only three kinds of deformation patterns above.Other way, make exactly active phase array antenna exemplar, according to actual measurement electrical property result, through too much wheel scheme, repeatedly carry out to determine structure tolerance, this brings maximum problem is exactly that the whole solution formulation cycle is long, design cost is high, can not meet the demand of current China radar rapid development.

Summary of the invention

The object of the invention is to propose a kind of fast determination method of active phase array antenna structure tolerance, to effectively solve, when active phase array antenna Scheme design, be difficult to determine fast and the problem of distributing antenna structure tolerance, can instruct formulation and the distribution of active phase array antenna structure processing technology and installation accuracy, and the evaluation of antenna structure scheme.

The technical solution that realizes the object of the invention is that a kind of fast determination method of active phase array antenna structure tolerance, is characterized in that: at least comprise the steps:

Step 1, determines antenna array radiating element information;

Step 2, obtains the height error of front radiating element;

Step 3, obtains the position alignment error of front radiating element;

Step 4, calculates Antenna aperture phase error;

Step 5, calculates antenna far field Electric Field Distribution;

Step 6, calculates the relevant unit for electrical property parameters of antenna;

Step 7, according to tactical and technical norms requirement, whether the antenna electric performance that judgement calculates under current structure tolerance conditions meets the demands, if do not meet the demands, carries out step 8 and step 9; If meet the demands, directly forward step 10 to;

Step 8, revises front flatness;

Step 9, revises front installation accuracy, and re-starts step 2 to step 7;

Step 10, the position of current radiating element and height tolerance be exactly definite antenna structure tolerance fast.

Described step 1 determines that antenna array radiating element information comprises the steps:

1.1. obtain antenna length L _x, width L _y, frequency of operation f and wavelength X information, obtain line number M, columns N and the spacing d of radiating element in x, y direction of front radiating element _x, d _yinformation; Being numbered of radiating element (i, j) in front;

1.2. the information of above-mentioned antenna array radiating element is formed to the data file of set form according to the order of row, column numbering.

Described step 2 is obtained the height error of front radiating element;

2.1. making the front flatness average that distributing is 0, and standard deviation is σ _znormal distribution stochastic error, choose prompt radiation cell height error, i.e. σ _zit is 1 times of wavelength X;

2.2. making the design coordinate of (i, j) individual radiating element is (id _x, jd _y, 0), when front flatness z is when existing stochastic error Δ z, radiating element coordinate becomes:

(i·d _x,j·d _y,Δz _ij) （1）。

Described step 3 is obtained the position alignment error of front radiating element;

3.1. making the front installation accuracy average that distributing is 0, and standard deviation is σ _x=σ _ynormal distribution stochastic error, choose prompt radiation cell position alignment error, i.e. σ _x, σ _yit is 1 times of wavelength X;

3.2. when radiating element (i, j) position alignment error exists after stochastic error Δ x and Δ y in x, y direction, the corresponding 0xy plane of radiating element internal coordinate becomes

(i·d _x+Δx _ij,j·d _y+Δy _ij) （2）；

When the site error amount of radiating element is (Δ x _ij, Δ y _ij, Δ z _ij), the new coordinate of radiating element is:

(i·d _x+Δx _ij,j·d _y+Δy _ij,Δz _ij) （3）。

Described step 4 is the results according to step 2 and step 3, calculates Antenna aperture phase error:

4.1. target setting is expressed as (cos α with respect to the direction (θ, φ) at coordinate system O-xyz place with direction cosine _x, cos α _y, cos α _z), target with respect to the angle of coordinate axis and the pass of direction cosine is

$\left\{\begin{array}{c}\cos {\mathrm{\α}}_x=\sin \mathrm{\θ}\cos \mathrm{\φ}\\ {\cos \mathrm{\α}}_y=\sin \mathrm{\θ}\sin \mathrm{\φ}\\ \cos {\mathrm{\α}}_z=\cos \mathrm{\θ}\end{array}\right.---\left(4\right);$

4.2. between adjacent two radiating elements of antenna, at target place (θ, φ), along the space quadrature of x-axis, y-axis and z-axis, be respectively:

Wherein, wave constant k=2 π/λ

Therefore, (m, n) individual radiating element (0≤m≤M-1,0≤n≤N-1) with respect to the phase differential of (0,0) individual radiating element is:

Wherein, (Δ x ₀₀, Δ y ₀₀, Δ z ₀₀) be the displacement of (0,0) individual radiating element, β _mnbe phase differential in (m, n) individual radiating element battle array, by phase shifter, control and realize beam scanning;

4.3. the phase differential with respect to reference radiation unit by each radiating element, becomes the form of matrix according to the sequential storage of radiating element Position Number, this matrix represents the phase differential of active phase array antenna actinal surface.

Described step 5 is the result of calculation according to step 4, calculates antenna far field Electric Field Distribution:

Draw the electric field value of Antenna Far Field region point (θ, φ); Change (θ, φ) numerical value, double counting process, can draw in certain concrete scope of far-field region electric field value a little, field value is taken the logarithm, and is plotted in coordinate system, can obtain the directional diagram of far field regional extent.

Described step 6 is the result of calculation according to step 5, calculates the relevant unit for electrical property parameters of antenna and comprises the steps:

6.1. from directional diagram, obtain the electrical properties such as minor level SLL, beam position BS and beam angle HPBW;

6.2. according to Electric Field Distribution E (θ, φ), calculate antenna gain

$G=10{\mathrm{<figure-callout\; id="10"\; label="log"\; filenames="HDA00001676178500051.png"\; state="\{\{state\}\}">log</figure-callout>}}_{10}\frac{4\mathrm{\&pi;}}{{\&Integral;}_0^{2\mathrm{\&pi;}}{\&Integral;}_0^{\mathrm{\&pi;}}{\left[E(\mathrm{\&theta;},\mathrm{\&phi;})\right]}^2\sin \mathrm{\&theta;d\&theta;d\&phi;}}---\left(8\right);$

6.3. based on Antenna Design index, calculate the reducing amount of the electrical properties such as antenna gain G, minor level SLL, beam position BS.

Described step 8: revise front flatness, make standard deviation sigma _zmultiple for wavelength X is distinguished value successively in the following order

$[\frac{1}{2},\frac{1}{4},\frac{1}{5},\frac{1}{8},\frac{1}{10},\frac{1}{15},\frac{1}{20},\frac{1}{25},\frac{1}{30},\frac{1}{40},\frac{1}{60},\frac{1}{80},\frac{1}{100}]---\left(9\right).$

Described step 9: revise front installation accuracy, make standard deviation 0 _x, σ _y(σ _x=σ _y) for the multiple of wavelength X, distinguish in the following order successively value

$[\frac{1}{2},\frac{1}{4},\frac{1}{5},\frac{1}{8},\frac{1}{10},\frac{1}{15},\frac{1}{20},\frac{1}{25},\frac{1}{30},\frac{1}{40},\frac{1}{60},\frac{1}{80},\frac{1}{100}]---\left(10\right).$

The present invention is due to the site error of front radiating element is incorporated in the phase error of antenna opening diametric plane by phase differential, antenna structure displacement field and electromagnetic field are closely connected, thereby avoided only by structure tolerance, judging the deficiency of antenna performance, thereby realized the electrical and mechanical comprehensive analysis of antenna; Antenna electric performance computing method owing to having built that front radiating element exists arbitrary height error, position alignment error, thereby can realize by random generation radiating element site error the evaluation of electrical property, avoided prior art only can analyze the problem that is difficult to engineering application that three kinds of front constraint types (flexural deformation, bowl-shape distortion, saddle facial disfigurement) cause, precision and speed that antenna electric performance calculates have been improved simultaneously, and the active phase array antenna structure that can analyze different frequency range, has good applicability.

Simulation result shows, structure tolerance when method of the present invention can realize the active phase array antenna solution formulation of different frequency range is determined fast, instruct formulation and the distribution of active phase array antenna structure processing technology and installation accuracy, and the evaluation of antenna structure scheme, make determining of antenna array structure tolerance more have explanation of force and feasibility.

Accompanying drawing explanation

Referring to accompanying drawing, the present invention is described in further detail:

Fig. 1 is that active phase array antenna structure tolerance of the present invention is determined process flow diagram;

Fig. 2 is that the radiating element of planar rectangular active phase array antenna is arranged schematic diagram;

Fig. 3 is radiating element height error schematic diagram;

Fig. 4 is radiating element position alignment error schematic diagram;

Fig. 5 is the space geometry graph of a relation of target;

Fig. 6 is flat surface active phased array antenna front schematic diagram;

Fig. 7 is the antenna radiation pattern that Shi Butong front flatness is corresponding;

Fig. 8 is

the antenna radiation pattern that Shi Butong front flatness is corresponding;

Fig. 9 is

the antenna radiation pattern that Shi Butong front installation accuracy is corresponding;

Figure 10 is

the antenna radiation pattern that Shi Butong front installation accuracy is corresponding.

In figure, sequence number represents respectively: 1, antenna array radiating element; 2, antenna length; 3, day line width; 4, line number; 5, columns; 6, the spacing of radiating element in x direction; 7, the spacing of radiating element in y direction.

Embodiment

(with reference to Fig. 1) is as follows for concrete steps of the present invention:

Step 1, determines antenna array radiating element information.

1.1. obtain antenna length L _x, width L _y, frequency of operation f and wavelength X information, obtain line number M, columns N and the spacing d of radiating element in x, y direction of front radiating element _x, d _yinformation (see figure 2); Being numbered of radiating element (i, j) in front; Wherein i, j are the radiating element numbering in x, y direction respectively, and front lower right-hand corner is Base Serial Number, and the radiating element of front lower right-hand corner is numbered (0,0), and this is also the true origin that is positioned at the coordinate system Oxy of front simultaneously; Front normal direction is exactly the z axle of coordinate system Oxyz;

1.2. the information of above-mentioned antenna array radiating element is formed to the data file of set form according to the order of row, column numbering.

Step 2, obtains the height error of front radiating element.

2.1. active phase array antenna structure tolerance shows as front flatness and front installation accuracy at front, and front flatness is definite to height error by the z of radiating element;

2.2. making the front flatness average that distributing is 0, and standard deviation is σ _znormal distribution stochastic error, choose prompt radiation cell height error, i.e. σ _zit is 1 times of wavelength X;

2.3. making the design coordinate of (i, j) individual radiating element is (id _x, jd _y, 0), when front flatness z (see figure 3) when there is stochastic error Δ z, radiating element coordinate becomes:

(i·d _x,j·d _y,Δz _ij) （11）。

Step 3, obtains the position alignment error of front radiating element.

3.1. front installation accuracy is determined by radiating element position alignment error in x, y direction, and the mounting means of radiating element is identical in x, y direction, therefore its installation accuracy is identical;

3.2. making the front installation accuracy average that distributing is 0, and standard deviation is σ _x=σ _ynormal distribution stochastic error, choose prompt radiation cell position alignment error, i.e. σ _x, σ _yit is 1 times of wavelength X;

3.3. (see figure 4) after radiating element (i, j) position alignment error exists stochastic error Δ x and Δ y in x, y direction, the corresponding 0xy plane of radiating element internal coordinate becomes

(i·d _x+Δx _ij,j·d _y+Δy _ij) （12）；

When the site error amount of radiating element is (Δ x _ij, Δ y _ij, Δ x _ij), the new coordinate of radiating element is:

(i·d _x+Δx _ij,j·d _y+Δy _ij,Δz _ij) （13）。

Step 4, according to the result of step 2 and step 3, calculates Antenna aperture phase error.

4.1. target setting is expressed as (cos α with respect to the direction (θ, φ) at coordinate system O-xyz place with direction cosine _x, cos α _y, cos α _z), according to the known target of Fig. 5, with respect to the angle of coordinate axis and the pass of direction cosine, be

$\left\{\begin{array}{c}\cos {\mathrm{\&alpha;}}_x=\sin \mathrm{\&theta;}\cos \mathrm{\&phi;}\\ {\cos \mathrm{\&alpha;}}_y=\sin \mathrm{\&theta;}\sin \mathrm{\&phi;}\\ \cos {\mathrm{\&alpha;}}_z=\cos \mathrm{\&theta;}\end{array}\right.---\left(14\right);$

4.2. between adjacent two radiating elements of antenna, at target place (θ, φ), along the space quadrature of x-axis, y-axis and z-axis, be respectively:

Wherein, wave constant k=2 π/λ

Therefore, (m, n) individual radiating element (0≤m≤M-1,0≤n≤N-1) with respect to the phase differential of (0,0) individual radiating element is:

Wherein, (Δ x ₀₀, Δ y ₀₀, Δ z ₀₀) be the displacement of (0,0) individual radiating element, β _mnbe phase differential in (m, n) individual radiating element battle array, by phase shifter, control and realize beam scanning;

4.3. the phase differential with respect to reference radiation unit by each radiating element, becomes the form of matrix according to the sequential storage of radiating element Position Number, this matrix represents the phase differential of active phase array antenna actinal surface.

Step 5, according to the result of calculation of step 4, calculates antenna far field Electric Field Distribution.

5.1. because changing, radiating element spacing will make radiating element the mutual coupling coefficient S _{mn, pq}change into S _{mn, pq}(Δ x, Δ y), therefore (m, n) individual radiating element directional diagram is

$f_{\mathrm{mn}}(\mathrm{\&theta;},\mathrm{\&phi;})=f_e(\mathrm{\&theta;},\mathrm{\&phi;})\&CenterDot;\underset{p=0}{\overset{M-1}{\mathrm{\&Sigma;}}}\underset{q=0}{\overset{N-1}{\mathrm{\&Sigma;}}}{S}_{\mathrm{mn},\mathrm{pq}}^{\&prime;}(\mathrm{\&Delta;x},\mathrm{\&Delta;y})\&CenterDot;\mathrm{expjk}\left\{\right[(p-m)\mathrm{dx}+{\mathrm{\&Delta;x}}_{\mathrm{pq}}-\mathrm{\&Delta;}{x}_{\mathrm{mn}}]\&CenterDot;\cos {\mathrm{\&alpha;}}_x+---\left(17\right);$

$[(q-n)\mathrm{dy}+{\mathrm{\&Delta;y}}_{\mathrm{pq}}-{\mathrm{\&Delta;y}}_{\mathrm{mn}}]\&CenterDot;{\cos \mathrm{\&alpha;}}_y\}$

Wherein, f _e(θ, φ) is radiating element ideal orientation figure (all radiating elements are identical here), S ' _{mn, pq}(Δ x, Δ y) is the mutual coupling coefficient of (p, q) individual radiating element to (m, n) individual radiating element, and its value is

Here radiating element the mutual coupling coefficient can utilize the Electromagnetic Simulation softwares such as HFSS to calculate.

The exciting current amplitude of 5.2. applying to each radiating element in feeding network is phase place is β _mn, therefore, consider that the radiating element exciting current amplitude after mutual coupling is

$I_{\mathrm{mn}}=I_{\mathrm{mn}}^e+\underset{p=0}{\overset{M-1}{\mathrm{\&Sigma;}}}\underset{q=0}{\overset{N-1}{\mathrm{\&Sigma;}}}{S}_{\mathrm{mn},\mathrm{pq}}(\mathrm{\&Delta;x},\mathrm{\&Delta;y})\&CenterDot;I_{\mathrm{pq}}---\left(19\right);$

5.3. therefore according to formula (17), (18), (19), calculate radiating element and exist the far field Electric Field Distribution of the active phase array antenna of height error and position alignment error to be

5.4. according to calculating antenna far field Electric Field Distribution, draw the far field direction of an electric field figure of antenna.Draw the electric field value of Antenna Far Field region point (θ, φ); Change (θ, φ) numerical value, double counting process, can draw in certain concrete scope of far-field region electric field value a little, field value is taken the logarithm, and is plotted in coordinate system, can obtain the directional diagram of far field regional extent.

Step 6, according to the result of calculation of step 5, calculates the relevant unit for electrical property parameters of antenna.

6.1. from directional diagram, obtain the electrical properties such as minor level SLL, beam position BS and beam angle HPBW;

From directional diagram, can directly read out the unit for electrical property parameters such as controlling antenna wave beam to point BS and beam angle HPBW, to former directional diagram normalization, can obtain normalized radiation pattern, therefrom can read minor level SLL.As for other electrical property, as efficiency, on the basis of existing electric field value, by corresponding calculating, can obtain.

6.2. according to Electric Field Distribution E (θ, φ), calculate antenna gain G

$G=10{\mathrm{<figure-callout\; id="10"\; label="log"\; filenames="HDA00001676178500051.png"\; state="\{\{state\}\}">log</figure-callout>}}_{10}\frac{4\mathrm{\&pi;}}{{\&Integral;}_0^{2\mathrm{\&pi;}}{\&Integral;}_0^{\mathrm{\&pi;}}{\left[E(\mathrm{\&theta;},\mathrm{\&phi;})\right]}^2\sin \mathrm{\&theta;d\&theta;d\&phi;}}---\left(21\right);$

6.3. based on Antenna Design index, calculate reducing amount or the deterioration degree of the electrical properties such as antenna gain G, minor level SLL, beam position BS.

Step 8, revises front flatness, makes standard deviation sigma _zmultiple for wavelength X is distinguished value successively in the following order

$[\frac{1}{2},\frac{1}{4},\frac{1}{5},\frac{1}{8},\frac{1}{10},\frac{1}{15},\frac{1}{20},\frac{1}{25},\frac{1}{30},\frac{1}{40},\frac{1}{60},\frac{1}{80},\frac{1}{100}]---\left(22\right).$

Step 9, revises front installation accuracy, makes standard deviation sigma _x, σ _y(σ _x=σ _y) for the multiple of wavelength X, distinguish in the following order successively value

$[\frac{1}{2},\frac{1}{4},\frac{1}{5},\frac{1}{8},\frac{1}{10},\frac{1}{15},\frac{1}{20},\frac{1}{25},\frac{1}{30},\frac{1}{40},\frac{1}{60},\frac{1}{80},\frac{1}{100}]---\left(23\right).$

Advantage of the present invention can further illustrate by following emulation experiment:

1. simulated conditions:

Adopt the X frequency range planar rectangular active phase array antenna schematic diagram providing as Fig. 6, the front radiating element 1 that distributing on the antenna of length 2 and width 3, its line number 4 is that M, columns 5 are N, the spacing 6 of radiating element in x direction is d _x; The spacing 7 of radiating element in y direction is d _y.

Whole antenna package is containing 1024 radiating elements, and design parameter is as shown in table 1, and supposes the even weighting of the exciting current employing constant amplitude homophase of Antenna aperture.

Table 1 antenna parameter

There is the antenna electric performance in height error and two kinds of structure tolerance situations of position alignment error in emulation front radiating element, a kind of is z impact on electrical property to different front flatnesses, and another kind is x, the y impact on electrical property to different front installation accuracies.

2. simulation result:

Because antenna is operated in high-frequency band, so according to actual processing and manufacturing ability, the standard deviation Initial value choice of stochastic error is 3.2mm,

(1) in front z direction, add standard deviation sigma _zbe respectively

normal distribution stochastic error Δ z after, emulation obtains

with

as shown in Figure 7 and Figure 8, corresponding unit for electrical property parameters is as shown in table 2 for the antenna radiation pattern of plane.

(2) in front x, y direction, add identical standard poor

normal distribution stochastic error Δ x, after Δ y, emulation obtains with as shown in Figure 9 and Figure 10, corresponding unit for electrical property parameters is as shown in table 3 for the antenna radiation pattern of plane.

Antenna electric performance parameter under table 2 Different Plane degree tolerance

Antenna electric performance parameter under table 3 diverse location location tolerance

For this 32x32 active phase array antenna that works in 9.375GHz: from Fig. 7, Fig. 8 and table 2, flatness tolerance is larger, antenna gain loss is larger, and when flatness tolerance is λ/15(2.13mm) time, antenna gain loss is 0.454dB, minor level is raised as 0.447dB, meets the requirement that is less than 0.5dB in engineering for gain loss; From Fig. 9, Figure 10 and table 3, position location tolerance is larger, and antenna gain loss is larger, and when position location tolerance is λ/20(1.6mm) time, antenna gain loss is 0.471dB, and minor level is raised as 0.114dB, meets the requirement that is less than 0.5dB in engineering for gain loss.Data in table, the minor level deterioration degree under known Different Plane degree tolerance will be apparently higher than corresponding position location tolerance, and flatness and the installation accuracy of therefore tackling active phase array antenna structure tolerance propose respectively rational requirement.

Above-mentioned simulation numerical evidence, adopts the present invention can effectively determine fast the structure tolerance of active phase array antenna.

Claims (8)

1. a fast determination method for active phase array antenna structure tolerance, is characterized in that: at least comprise the steps:

Step 1, determines antenna array radiating element information;

Step 2, obtains the height error of front radiating element;

Step 3, obtains the position alignment error of front radiating element;

Step 4, calculates Antenna aperture phase error;

Step 5, calculates antenna far field Electric Field Distribution;

Step 6, calculates the relevant unit for electrical property parameters of antenna;

Step 7, according to tactical and technical norms requirement, whether the antenna electric performance that judgement calculates under current structure tolerance conditions meets the demands, if do not meet the demands, carries out step 8 and step 9; If meet the demands, directly forward step 10 to;

Step 8, revises front flatness;

Step 9, revises front installation accuracy, and re-starts step 2 to step 7;

Step 10, the position of current radiating element and height tolerance be exactly definite antenna structure tolerance fast;

In described step 1, determine that antenna array radiating element information comprises the steps:

1.1. obtain antenna length L _x, width L _y, frequency of operation f and wavelength X information, obtain line number M, columns N and the spacing d of radiating element in x, y direction of front radiating element _x, d _yinformation; Being numbered of radiating element (i, j) in front;

1.2. the information of above-mentioned antenna array radiating element is formed to the data file of set form according to the order of row, column numbering.

2. the fast determination method of a kind of active phase array antenna structure tolerance according to claim 1, is characterized in that: described step 2 is obtained the height error of front radiating element;

2.1. making the front flatness average that distributing is 0, and standard deviation is σ _znormal distribution stochastic error, choose prompt radiation cell height error, i.e. σ _zit is 1 times of wavelength X;

2.2. making the design coordinate of (i, j) individual radiating element is (id _x, jd _y, 0), when front flatness z is when existing stochastic error Δ z, radiating element coordinate becomes:

(i·d _x,j·d _y,Δz _ij) （1）。

3. the fast determination method of a kind of active phase array antenna structure tolerance according to claim 1, is characterized in that: described step 3 is obtained the position alignment error of front radiating element;

3.1. making the front installation accuracy average that distributing is 0, and standard deviation is σ _x=σ _ynormal distribution stochastic error, choose prompt radiation cell position alignment error, i.e. σ _x, σ _yit is 1 times of wavelength X;

3.2. when radiating element (i, j) position alignment error exists after stochastic error Δ x and Δ y in x, y direction, the corresponding Oxy plane of radiating element internal coordinate becomes

(i·d _x+Δx _ij,j·d _y+Δy _ij) （2）；

When the site error amount of radiating element is (Δ x _ij, Δ y _ij, Δ z _ij), the new coordinate of radiating element is:

(i·d _x+Δx _ij,j·d _y+Δy _ij,Δz _ij) （3）。

4. the fast determination method of a kind of active phase array antenna structure tolerance according to claim 1, is characterized in that: described step 4 is the results according to step 2 and step 3, calculates Antenna aperture phase error:

4.1. target setting is expressed as (cos α with respect to the direction (θ, φ) at coordinate system O-xyz place with direction cosine _x, cos α _y, cos α _z), target with respect to the angle of coordinate axis and the pass of direction cosine is

$\left\{\begin{array}{c}\cos {\mathrm{\&alpha;}}_x=\sin \mathrm{\&theta;}\cos \mathrm{\&phi;}\\ \cos {\mathrm{\&alpha;}}_y=\sin \mathrm{\&theta;}\sin \mathrm{\&phi;}\\ \cos {\mathrm{\&alpha;}}_z=\cos \mathrm{\&theta;}\end{array}\right.---\left(4\right);$

4.2. between adjacent two radiating elements of antenna, at target place (θ, φ), along the space quadrature of x-axis, y-axis and z-axis, be respectively:

Wherein, wave constant k=2 π λ

Therefore, (m, n) individual radiating element (0≤m≤M-1,0≤n≤N-1) with respect to the phase differential of (0,0) individual radiating element is:

Wherein, (Δ x ₀₀, Δ y ₀₀, Δ z ₀₀) be the displacement of (0,0) individual radiating element, β _mnbe phase differential in (m, n) individual radiating element battle array, by phase shifter, control and realize beam scanning;

4.3. the phase differential with respect to reference radiation unit by each radiating element, becomes the form of matrix according to the sequential storage of radiating element Position Number, this matrix represents the phase differential of active phase array antenna actinal surface.

5. the fast determination method of a kind of active phase array antenna structure tolerance according to claim 1, is characterized in that: described step 5 is the result of calculation according to step 4, calculates antenna far field Electric Field Distribution:

Draw the electric field value of Antenna Far Field region point (θ, φ); Change (θ, φ) numerical value, double counting process, can draw in certain concrete scope of far-field region electric field value a little, field value is taken the logarithm, and is plotted in coordinate system, can obtain the directional diagram of far field regional extent; Wherein, I _mnfor the radiating element exciting current amplitude after mutual coupling; f _mn(θ, φ) is (m, n) individual radiating element direction.

6. the fast determination method of a kind of active phase array antenna structure tolerance according to claim 1, is characterized in that: described step 6 is the result of calculation according to step 5, calculates the relevant unit for electrical property parameters of antenna and comprises the steps:

6.1. from directional diagram, obtain the electrical properties such as minor level SLL, beam position BS and beam angle HPBW;

6.2. according to Electric Field Distribution E (θ, φ), calculate antenna gain

$G=10{\log }_{10}\frac{4\mathrm{\&pi;}}{{\&Integral;}_0^{2\mathrm{\&pi;}}{\&Integral;}_0^{\mathrm{\&pi;}}{E(\mathrm{\&theta;},\mathrm{\&phi;})]}^2\sin \mathrm{\&theta;d\&theta;d\&phi;}}---\left(8\right);$

6.3. based on Antenna Design index, calculate reducing amount or the deterioration degree of the electrical properties such as antenna gain G, minor level SLL, beam position BS.

7. the fast determination method of a kind of active phase array antenna structure tolerance according to claim 1, is characterized in that: described step 8: revise front flatness, make standard deviation sigma _zmultiple for wavelength X is distinguished value successively in the following order

$[\frac{1}{2},\frac{1}{4},\frac{1}{5},\frac{1}{8},\frac{1}{10},\frac{1}{15},\frac{1}{20},\frac{1}{25},\frac{1}{30},\frac{1}{40},\frac{1}{60},\frac{1}{80},\frac{1}{100}]---\left(9\right).$

8. the fast determination method of a kind of active phase array antenna structure tolerance according to claim 1, is characterized in that: described step 9: revise front installation accuracy, make standard deviation sigma _x, σ _y(σ _x=σ _y) for the multiple of wavelength X, distinguish in the following order successively value

$[\frac{1}{2},\frac{1}{4},\frac{1}{5},\frac{1}{8},\frac{1}{10},\frac{1}{15},\frac{1}{20},\frac{1}{25},\frac{1}{30},\frac{1}{40},\frac{1}{60},\frac{1}{80},\frac{1}{100}]---\left(10\right).$

Images