Open access peer-reviewed chapter

Synthesis of Antenna Systems According to the Desired Amplitude Radiation Characteristics

Written By

Mykhaylo Andriychuk

Submitted: November 18th, 2011Reviewed: April 24th, 2012Published: September 19th, 2012

DOI: 10.5772/46164

Chapter metrics overview

1,785 Chapter Downloads

View Full Metrics

1. Introduction

The problems of antenna synthesis, which use the amplitude RP as input information, are often used in the process of antenna design for many practical applications [4, 16, 27, 31]. In spite of the fact that the respective mathematical problems are ill-posed [32] and they have the branching solutions [29], the antenna synthesis according to the desired amplitude characteristics is very useful and perspective.

As a rule, the branching of solutions depends on the properties of prescribed amplitude RP, geometrical and physical parameters of the considered antenna. The methods of nonlinear functional analysis [23] allowing to localize the branching solutions are applied for investigation of solutions and determination of their number and qualitative characteristics. Such approach too much simplifies determination of the optimal solutions by the numerical methods. The iterative processes for the numerical solving of the corresponding non-linear equations were elaborated in [3, 8, 11].

The Chapter is organized as follows.

In Section 2 we derive the main formulas for RP of antennas and introduce the objective functionals for the synthesis problem. Also in Section 2 we consider the variational statement of problems and derive the fundamental nonlinear equations of the synthesis.

Section 3 contains the application of the proposed approach to several types of antennas. Depending on the restrictions which are imposed on the sought distribution of current or field in the antenna elements and type of antenna, the problems of amplitude-phase, amplitude, and phase synthesis are considered for the specific antennas. The methods of successive approximations are applied for solving the derived non-linear integral equations; the convergence of the elaborated methods is discussed. The direct optimization of the proposed functionals by the gradient methods is performed and successfully applied to solving the amplitude and phase synthesis problems.

In Section 4 conclusions are formulated.


2. The theoretical background of the synthesis problems

In this Section, we present the necessary information about the properties of the electromagnetic (EM) field in far zone, introduce the variational approach to the antenna synthesis problems, as well as discuss the arising nonlinear integral and matrix equations.

2.1. RP of electromagnetic field

The EM field in the non-limited homogeneous medium satisfies the Maxwell equations [25]


where Eis a component of electric field, His a component of magnetic field, Iis the extrinsic current density, ϟis the volume density of electric charge.

The scalar and vector potentials are used for solving the equations (1)-(4). Introducing the vector potential Aeas [34]


we satisfy the equation (4). Substituting (5) into (2), we receive the equation


which testifies that the vector field in the parenthesis of (6) is potential. This yields the equation


where Ϥeis the scalar potential. Substituting (5) and (7) into (1), we receive the following equation for I:


Using Lorentz lemma [15]


we receive the inhomogeneous Helmholtz equation for the vector potential


where k=ϧϓϚis a propagation coefficient. For a free space, the values ϓand Ϛare real and related with the light velocity as c=1/ϓϚ, and coefficient k=ϧ/c=2Ϟ/ϙis the wavenumber, ϙis the length of wave.

The vector potential Aein the arbitrary observation point Q(x,y,z)is determined by formula [26]


where RPQ=(x-x')2+(y-y')2+(z-z')2is the distance between the observation point Q(x,y,z)and integration point P(x',y',z'), I(x',y'z')is the density of current in the domain V.

Using formulas (5), (7), and (11), we receive the solution to system (1)-(4) in the term of electric vector potential


In the process of solving the synthesis problem, the representation of field in far zone is of interest. Using the approximate representation of distance RPQin far zone, we receive the formula for vector potential in this region


where Ϗis the angle between the vectors, directed into observation and integration point, r'and rare the radius-vector of the points Pand Q, respectively.

Substituting (14) into (13) and neglecting by terms of O(1r2)order, we receive the formulas for components of magnetic field Hin term of vector potential Ae. These formulas in the spherical coordinates have form


Consequently, the formulas for Eare the following


Using formulas (15), (16) and relation


we receive the formulas for non-zero Ecomponents of EM field


Function e-ikr4Ϟris a spherical wave and it depends on ronly. The second terms in (18) and (19) are the functions of angular coordinates of the observation point and are determined by the current I(x',y',z')


The functions (20) and (21) characterize the angular distribution of Ecomponents in far zone.

In such a way, the vectors Eand Hin far zone are expressed by formulas [Andriychuk et al., 1993]


The functions fϖ(ϖ,Ϥ)and fϤ(ϖ,Ϥ)are defined as the RPs. Function


is the power RP and it characterizes the angular distribution of power intensity radiation. The functions fϖϖ,Ϥ, fϤϖ,Ϥ, and N(ϖ,Ϥ)are used in the process of formulating and solving the synthesis problems for various antennas.

2.2.Variational. statement of the synthesis problems

Abstracting of the specific type of antenna, we present the functions fϖand fϤin the formulas (22)-(24) by the linear operator A={Aϖ,AϤ}:


acting from some complex Hilbertian space HI, to which the distribution functions of current or field belong, into the complex space Cf2=CχCχof vector-valued continuous functions on the compact χ̣R2(or χ̣R1)[29]. The form and properties of the operators Aϛdepend on a type and geometry of antenna. In many practical applications, one can reduce the synthesis problem to separate consideration of the fϖandfϤcomponents. This allows to reduce the synthesis problem to the scalar one. In that way, we will use more simple formula


for operator expression of the RP f.

The specific form of the operator Adepends on the antenna type. This operator is integral for the continuous antennas. As an example, for a cylindrical antenna with curvilinear generatrix and with current polarized along the cylinder axis, the RP in the transversal plane has form [21]


where fand Iare nonzero components fϖand Izrespectively; Ϥis the angular coordinate of the point in far zone, Ϥ'is the angular coordinate of the point in antenna, r=r(Ϥ')describes the generatrix in polar coordinates, dSϤ'=r2+(dr/dϤ')2is an element of arc. In the case of array, the operator Ais described by a finite sum.

In the previous subsection, we consider the direct external problem of electrodynamics consisting of determination of the asymptotic (RP) of EM field in far zone. The inverse problem, namely determination of such a current Ithat create EM field with the desired RP f, is of specific interest in the antenna design. The characteristic parameters of antenna can be fixed or be subject to determination in process of solving this problem. In such interpretation, the inverse problem is defined as the synthesis problem, namely the problem of determination of the current according to the desired RP.

The RPs are due to satisfy a series of requirements for the main lobe and sidelobes. In one case, the RP with narrow main beam is required, another time, this beam should have the specific wide (for example, cosecant) form; the sidelobes be as low as possible.

The angular distribution of the radiation power is characterized by the amplitude of RP, but not the whole RP. Therefore, only amplitude |f|of function fis interesting in the process of statement and solving the synthesis problem. In this case, the freedom of choice of the phase argfof function fis used for better approximation to the amplitude RP. In the synthesis theory, function |f|is amplitude RP, and function argfis phase RP.

In this way, in the process of synthesis we prescribe not the whole complex function f, but only its amplitude. We denote this function as F, the created (synthesized) by antenna amplitude RP is denoted by |f|. The both functions are real and positive. Of course, these functions can not coincide in any real-world situation. This fact yields to use the variational statement of the synthesis problem. In such statement, one requires not the whole coincidence of the functions |f|and F, but the better approximation of function |f|to Fin a certain sense only. The mean-square deviation of both RPs is used as the criterion of optimization.

Let us introduce the Hilbertian spaces of radiation patterns Hfand currents HI. Let (,)fbe an inner product in Hf, and norm of fis defined as [4]


The functional


is used in most cases as optimization criterion. Dependence of ϡon the sought distribution of the current Iin antenna is determined by formula (26). The additional multiplier sin (29) can be either prescribed (for example, s=1)or determined from the condition ϡ/s=0. In the latter case


Let us introduce the normalization (F,F)f=1, then (29) can be written as



and ||I||is the norm of current, determined by inner product (,)Iin the Hilbertian space HIof currents: ||I||2=(I,I)I. The factor Qcharacterizes the goodness (reactivity) of antenna, and Ϙdetermines its power efficiency, namely the part of power radiated by the RP.

The generalized functional, allowing to diminish the mean-square deviation of RPs and relative norm of current, has the form


Choosing the weight multiplier t, one can regularize the value of mean-square deviation of the RPs and norm of current. Similar functional appears when the method of Lagrange multipliers [16] is used to solving the conditional extremum problems. Simultaneously, one can use some weight function p0in the definition (28) and to improve the approximation to prescribed amplitude RP Fin the appointed angular range.

2.3. The fundamental equations of synthesis

Let us demonstrate how one can receive the respective Lagrange-Euler’s equation [17] in the process of minimization of functional ϡ. Let s=1. It is known [30] that the equality to zero of functional’s gradient is requirement of its extremum. This yields determination of the maximum of the following value


In order to determine the derivative ϡz, one requires to factorize ϡ(I+ϓz)in series relatively to ϓ


and to extract the linear term ϒϡ. Evidently, ϡz=ϒϡ||z||. The increment of amplitude RP |f|should be known in order to calculate ϒϡ. This increment has the form [4]


Substituting this expression into (35), we receive




where A*is an operator adjoint to Ain the following sense:

If Ais an integral operator then A*is the same with the complex conjugate kernel and integration with respect to the second argument.

Using the Cauchy-Bunyakovsky-Schwarz inequality [10] and maximizing ϡz, we receive the expression for gradient of ϡ:


which is used usually for the numerical minimization of ϡ. In order to receive the Lagrange-Euler’s equation one should equate to zero the function z(condition of ϡminimum)


The equation (42) contains fas an unknown function. One can turn out that its solutions represent the unreliazable patterns [20]. In order to avoid such solutions, it is necessary to substitute formula (26) into (42) instead of f. In this way, we receive the integral equation for the optimal current distribution


If operator Ais integral (consequently, A*is also integral) then (43) is the nonlinear Hammerstein equation [38].

One can receive the similar expression for zin the process of Ϙmaximization:


Equating the right part of (44) to zero, we receive the explicit expression for I:


Equations (26) and (45) yield the system of nonlinear integral equations for the optimal current Iand RP fcreated by it.

System of nonlinear equations (26), (45) can be reduced to one nonlinear equation


Since right hand side of (46) is a result of acting by the operator Aon some function, any solution to this equation represents the realizable RP. Once this equation is solved, the optimal distribution of current is determined by formula (45). The equation (46) is more simple than system (26), (45) because determination of its solution does not require additional operation of Aand A*.

Taking into account the above considerations, we receive the following nonlinear Hammerstein equation of the second kind for the RP f


in the case of functional ϡtminimization. As in the case of (46), the solution to this equation (at t0)represents the realizable RP. Having the solution to (47), we determine the optimal current distribution by formula


The integral equations (43), (46), and (47) are the fundamental equations for the synthesis problems according to the prescribed amplitude RP F.

2.4. The numerical solution of the integral equations

The equations (43), (46), and (47) are solved numerically by the method of successive approximations. The new approximation is determined by explicit formula


in the process of Ϙmaximization. The equation


is solved when the synthesis problem is formulated in the term of functional ϡminimization. If functional ϡtis minimized then the respective iterative procedure has form


In the last case, one can use the iterative process


but the convergence domain of this process is limited.

In accordance with the procedure used in Subsection 2.3, the solutions to the nonlinear equations (43), (46), and (47) are the stationary points of the respective functionals. Since the used functionals are nonconvex, the several solutions can appear, what corresponds to existence of several local minima or saddle points. The number of solutions can vary depending on the physical parameters of the problem what requires the special careful analysis of the obtained solutions [11].


3. Application to specific antennas

3.1. Amplitude-phase synthesis of cylindrical antenna

In this Subsection, we consider the application of variational approach to the synthesis problem of cylindrical curvilinear antenna. Functional (34) is used as criterion of optimization, complex function Iis optimizing function. Since both the amplitude and phase of complex function Iare optimizing parameters, the considered problem is the amplitude-phase synthesis problem.

3.1.1. The integral equation approach

Let the generatrix of the antenna has the length 2land be parallel to Ozaxis, and form of cross-section be determined by close curve Swhich be described by formula r=r(Ϥ'), where Ϥ'is the angular coordinate on S.

In many practical applications, the antennas with currents linearly polarized along Ozaxis are of interest. For such case, the RP has only the component fϖ(ϖ,Ϥ), denote it as f(ϖ,Ϥ). On account of formula (18):


Let the current distribution in antenna surface be determined as


then the spatial RP is completely determined by the RP f1(ϖ,Ϥ), created by the distribution of current I1(r(Ϥ'))in S, and RP f2(ϖ)in the longitudinal plane, created by the current distribution in the generatrix of cylinder


Function f2(ϖ)is the RP of a linear antenna with length 2l. Function


is the RP of a plane curvilinear antenna with form S. Consequently, one can reduce the synthesis problem of cylindrical antenna to two independent problems for synthesis of linear and plane curvilinear antennas.

There is a many literature sources on the synthesis problem of linear antennas (see, e.g. [9] and references there). Therefore we consider here the synthesis problem for plane curvilinear antenna.

It is easily seen from formula (53) that the angle ϖdetermines effective electrical scale of antenna. Therefore one can suppose ϖ=Ϟ/2. Omitting the indices "1"’ in the distribution of current and RP, we represent the RP (56) in form


that is, RP is determined due to action of linear bounded operator A.

The amplitude-phase synthesis problem for closed plane curvilinear antenna according to desired amplitude RP Fconsists of determination of such distribution of the current I(Ϥ'), that the amplitude RP |f(Ϥ)|created (synthesized) by it, is the most close to F(Ϥ). The functional (34) is used as the criterion of optimization. The inner products in the spaces of the RPs and currents are defined as


Using (47), (57), and definitions (58) and (59), we receive the nonlinear equation with respect to RP f


where K(Ϥ,Ϥ1)is the kernel of operator AA*


and nonlinear operator Bis determined as


The methods of successive approximations are applied for solving the nonlinear equation (60). The simplest of them


has the limited region of convergence determined by formula


where lis the length of contour S. Once the function fis found, the optimal current Iis determined by formula (48).

The iterative process (51) is more preferable, it yields the converging sequence of functional ϡtwhich satisfies the condition


for arbitrary t.

3.1.2. The gradient methods of optimization

Above we mentioned the methods of successive approximations for solving the arising nonlinear equations. The direct optimization of ϡtfunctional by the gradient methods can be also applied for solving the synthesis problem. The simplest gradient method is defined by the formula


where ϒnis an optimizing multiplier, znis a gradient of functional ϡt(34) on the function In:


where anis a number determined by known values in the n+1-th iteration [4].

The disadvantage of method (66) is that only the information about optimizing function from previous iteration is used, in addition it has the slow convergence at the end of iterative process. The method of conjugate gradients [28]


and proposed in [3] generalized gradient method


do not have such disadvantage. Here hnis a combination of znfrom previous iterations, and rn(1)=In,rn(2)=zn,rn(3)...,rn(M)is a set of some orthogonal functions, ϒnand ϒn(m)are the coefficients subject to determination.

For method (70), the problem of minimization of ϡtis reduced to the solution of nonlinear algebraic system


with unknown ϒn(m). This system is solved effectively by the method of successive approximations substituting in its right hand side the function argfnfrom previous iteration. Such iterative process is converging and similarly to iterative process (51) yields the converging sequence of ϡtnwhich satisfies the condition (65).

3.1.3. The numerical results

The numerical results are shown for the prescribed amplitude RPs F(Ϥ)=sin2(Ϥ/2)and F(Ϥ)=sin128(Ϥ/2)(Fig. 1a and Fig. 1b respectively). The influence of the parameter tin the functional (34) on the quality of synthesis is investigated. One can see that decrease of timproves the proximity of given Fand synthesized |f|RPs. But the norm ||I||of currents grows if tdecreases. In the case of narrow Fit is necessary to diminish tin order to decrease the mean-square deviation of the RPs. The detailed information about the synthesis quality is shown in Table 1.


Table 1.

The values of ϡand Ϙfor two desired RPs F(Ϥ)

Figure 1.

Dependence of synthesis quality on the parametertin the functionalϡt

3.2. The problem of amplitude synthesis for resonant antennas

Resonant antennas are a new type of antennas [24], which allow to form the radiation characteristics satisfying a wide spectrum of practical requirements. Such antennas are formed by several surfaces, one of which is semitransparent. Antennas with one semitransparent and other metal boundary are considered here.

The synthesis problem consists of determination of such parameters of antenna (the geometry of inner boundary and transparency of the outer boundary), which form the amplitude RP or front-to-rear factor (FRF) the most close to the prescribed ones.

3.2.1. Generalization of variational statement

The generalized method of eigen oscillations [1] is the mathematical basis for solving the analysis (direct) problem of resonant antennas. The two-dimensional model of antennas (the case of E- polarization) is considered.

The main constructive parameter of resonant antennas is the cophased field in the outer surface. This field can be considered quite real (i.e., only its amplitude can be considered) since the constant phase shift of field does not change the amplitude RP. In this connection, the synthesis problem for resonant antennas is formulated as the amplitude synthesis problem.

The direct problem consists of determination of the RP f(Ϥ)by the known field v(S)in the outer boundary Sof the given form. The RP created by this field can be presented similarly to (26). The operator Ain the case of circular external boundary has form [37]


where kernel

ϒ0nis the Kronecker delta, function Hn2(ka)is the Hankel function of second kind, ais the circle radius.

In the case of resonant antenna with arbitrary outer boundary, the method of auxiliary sources [2, 12, 18] is used for determination of the RP fby the field v. In this method, the field u(r,Ϥ)outside of antenna is represented approximately by the finite sum


where Rn=r2+rn2-2rrncos(Ϥ-Ϥn)is the distance between an observation point and n-th auxiliary source; r,Ϥand rn,Ϥnare the polar coordinates of a point of observation and n-th source, respectively; anare the unknown coefficients subject to determination in the process of solving the synthesis problem.

The RP is given by


Solving the synthesis problem, we determine the field in outer boundary Sof the antenna and transparency of this boundary. The form of inner boundary S0is determined as a curve of constant phase of the field u(r,Ϥ)[37].

Additionally, the restrictions on a field in some areas of a near zone can be prescribed. The functional


which is generalization of (34), allows to take into account these requirements. Here Ui(Si)are prescribed values of the field’s amplitude in the areas of restriction, ui(Si)are the obtained values of the field. Functions p(Ϥ), pi(Si)are the weight functions, allowing to adjust a degree of proximity of the given and received values of RP and field, tis the parameter limiting norm of the field v.

In the first step of solving the synthesis problem, the field von the outer boundary Sis determined from a condition of minimum of the functional (76). Minimization of functional can be carried out by the gradient methods, or by solving the respective Lagrange-Euler’s equation. In the first case, the generalized gradient method [3]


is used. The gradient zof the functional (76) (by virtue of the requirement of real field) has form


where A*and Bi*are the operators adjoined to Aand Bi, respectively [4]. Each step of iterative process (77) reduces ϡt. Since ϡtis limited from below (ϡt0), the process (77) is converging.

In the second step, the transparency ϟof the outer boundary Sand the form of inner metal boundary S0are determined.

For the antenna with circular outer boundary, the transparency distribution can be presented in the explicit form


where Jnand Nnare the Bessel and Neumann functions, respectively.

In the case of antenna with arbitrary outer boundary, similarly to [37], the distribution of transparency is determined by the formula


3.2.2. The numerical results

The numerical calculations are carried out for the resonant antenna with a given outer elliptic boundary. The prescribed amplitude RP is: F(Ϥ)=sin8(Ϥ/2). In Fig. 2, the results are presented for the antenna with parameters kb=15and different ka: ka=12.75,ka=14.25, where band aare the big and small semiaxes of ellipse. For such antenna the level of side lobes in the synthesized amplitude RP |f|is smaller than –20 dB, and distribution of transparency is smoother in the area of main radiation (the continuous lines in Figs. 2a, 2b correspond to ka=12.75, and the dashed ones correspond to ka=14.25). The outer elliptic boundaries (dashed lines), the found form of inner metallic boundaries (continuous lines), and the inner contour Saof placement of the auxiliary sources (dash-and-dot lines) are shown in Figs. 3a, 3b. The auxiliary sources are distributed uniformly on the Sa.

Figure 2.

Synthesis of resonant antenna with elliptic outer boundary

The distribution of transparency ϟin the area opposite to direction of main radiation has a spasmodic character. Such distribution cannot be realized by the physical reason. Therefore the values of ϟare averaged in this range in order to receive the smooth distribution of ϟ. This leads to some change of field von S, but the numerical calculations show small change of the synthesized amplitude RP |f|. The more smooth distribution of ϟin the area mentioned above can be achieved by increasing the number of auxiliary sources here.

The numerical results for solution of the synthesis problem with restrictions on the field in a near zone are given for the antenna with circular outer boundary. The prescribed amplitude RP is: F(Ϥ)=sin8(Ϥ/2); ka=15. Minimization of a field was carried out in two points Ϥ=Ϟ/2,3Ϟ/2on the additional circle with radius kb=20. These points were allocated in the second summand of the functional (76) using the weight function p1(Ϥ)=ϒ(Ϟ/2,3Ϟ/2); p(Ϥ)1, t=0.01. In Fig. 4a, the prescribed RP F(thick continuous line) and synthesized |f|(thin continuous line) amplitude RP are shown. The amplitude |u1|of obtained field on the circle of restrictions is marked by dashed line.

Figure 3.

Form of antenna boundariesSandS0

Figure 4.

Synthesis of resonant antenna with restriction on the field

It can be seen that the field at restriction points is reduced up to level -37 dB. The synthesized field v(continuous line) and transparency ϟ(dashed line) are presented in Fig. 4b. The form of the inner synthesized metallic boundary is more complicate than in the previous example.

3.3. Waveguide resonant antenna

Synthesis of resonant antenna with waveguide excitation is carried out according to the FRF. The optimizing functional enables to take into account a various requirements to the FRF of antenna in the operating frequency range, as well as outside this range.

3.3.1. The physical description of problem

The geometrical parameters of resonant antenna with waveguide excitation are shown in Fig. 5. In order to create RP enough narrow, the width Lof antenna should be much larger than the wavelength ϙ. The height dis of the order of ϙ/2. Excitation is carried out by a metal single-mode waveguide with semitransparent grid at its end; the width lof waveguide is of the order ϙ/2, and both its length and the length Dof the antenna along the Oxaxis are of the order of L.

Figure 5.

Geometry of the resonant antenna

The RP of antenna has form


The direct (analysis) problem on determination of the electromagnetic field components in the semitransparent aperture is reduced to two separate problems in the planes xOzand yOzrespectively. We consider here the case of E-polarization. The unknown function uis the Eycomponent of electromagnetic field. In the region over the antenna this component satisfies the Helmholtz equation and the boundary conditions


on the all metallic walls;


on the semitransparent upper boundary, the same conditions in the aperture of the exciting waveguide; the condition of radiation on the infinity


and asymptotical condition in the exciting waveguide


The problem of determination of the field u(x,y)in the semitransparent boundary Sis solved in three steps [37]. In the first step, the field in the irregular region of antenna is determined using the cross-section method [19]. In the second step, the field over the exciting waveguide is sought for, and the matching of field in the regular and irregular regions of antenna is fulfilled. The reflection coefficient R1is determined by the fulfillment of the adjoint boundary conditions (85) in the third step.

Under condition of the linear polarization of field in the aperture of exciting waveguide, the RP (81) can be represented as product of two functions, namely the RP of plane antenna with variable height d(y)of the upper wall and transparency ϟ(y)of the lower wall, and RP of the linear antenna in the xOzplane.

In such a way, the RP in the yOzplane can be written down in the form [37]


Numerical calculations can be essentially simplified, if one assumes that the field above the antenna is represented approximately in the form [37]


3.3.2. The objective functionals

In the process of statement of the synthesis problem one requires to provide the best approximation to the prescribed FRF in the operating frequency range [k1,k2]and the minimal values of the FRF outside this range, that is for k[k0,k1)and k(k2,k3].

The variational approach for solving this problem was developed in [37]. Modification of variational statement of the synthesis problem is proposed. The problem consists of determination of functions d(y)and ϟ(y), which maximize FRF ϕ1in the operating range [k1,k2]. In this case, the least value of FRF in this range is specified as a criterion of optimization, and this value is maximized by a choice of functions dand ϟ. That is, the functional is maximized


The additional parameter of optimization 1-|R1(k)|2, where R1(k)is the reflection factor of the main wave in exciting waveguide, with some weight multiplier can be included into functional (88). In this case, the transparency of waveguide aperture can be also used as an additional parameter of optimization.

Minimization of FRF outside the operating frequency range is one of requirements of electromagnetic compatibility for radiating systems [36]. Thus, the value of FRF should remain the largest in the main frequency range. Under these requirements, the following generalization of the variational statement of problem is considered: to find functions dand ϟ, maximizing the functional ϕ1, and, at the same time, minimizing additionally the following functional [6]


that is, minimization of the maximal FRF value outside the [k1,k2]range is required. Of course, it is necessary to take into account the restrictions on the functions d(x)and ϟ(x)owing the physical reason:


For example, d0,dmare the boundary values of height, which provide the single-mode conditions in antenna, and ϟ0,ϟmare the boundary values of transparency, which provide a good quality of antenna in the required ranges. Moreover, the received functions should satisfy the additional conditions of smoothness, which can be formulated as restriction on the second derivatives


3.3.3. The modeling results

The numerical results are shown in Fig. 6. Calculations were carried out for the problem of ϕmaximization in range ±5.0%in neighborhood of the central frequency (kc=6.0). It is necessary to minimize the FRF outside of this range for 5.2<k<5.6and 6.4<k6.8. Parameters of antenna are the following: the length of antenna L=6.0, half-width of excited waveguide l/2=0.3L, number of considered reflected waves in waveguide N1=5, number of waves in a background N2=20.

The optimized values of FRF are marked by solid line in the basic and additional ranges; the dashed line corresponds to not optimized values of ϕin the main and additional ranges.

Figure 6.

Optimized values of FRF in the basic and additional ranges versus the frequency (k)

Process of additional optimization is carried out on the simplified procedure, that is the control of decrease of the FRF in the basic range is omitted [Andriychuk Zamorska, 2004]. Therefore, the values of FRF in the basic range are slightly decreased in comparison with the FRF values for initial problem.

The optimal form d(y)of the lower boundary of antenna and transparency ϟ(y)of the upper boundary are slightly different for the both cases of optimization.

In Fig. 7, the change of the FRF values at three points of main and additional ranges of frequency (two extreme points and middle one) is shown. The width of main range is equal to 8.33%, and width of additional range is equal to 10.0%. In Fig. 7a, curve 1 corresponds to the central value of frequency k=12.0, and curves 2 and 3 correspond to the extreme points k=11.0and k=13.0respectively. In Fig. 7b, curve 1 corresponds to the central frequency k=14.0, and curves 2 and 3 correspond to the frequency values k=13.01and k=15.0. One can see that the main optimization takes place in the first iterations; there is the improvement only in low-order digit in the next steps. Therefore, it is enough to make 3-5 steps for the main range and 9-11 steps for the additional range, in iterative procedure to receive the practically interesting results.

The total number of iterations also depends on the width of the considered frequency ranges.

Figure 7.

The FRF values versus the iteration numberN

3.4. Phase synthesis problem for cylindrical array

The phase distributions of excitation currents in the array’s elements are the optimizing parameters in the problem of phase synthesis. The optimization of considered functionals is reduced to the solution of the corresponding system of nonlinear equations. The gradient methods for direct optimization of functionals are used in practical applications.

3.4.1. The RP of array

A series of simplifications in process of synthesis of the cylindrical array is used in order to reduce the computing time [9, 22]. Separation of variables onto the vertical and horizontal components for the distributions of currents in the array elements, as well as for the RPs, is one of the simplifications.

Thus, the expressions for current Inm(x,y)in radiators and RP f(ϖ,Ϥ)have the form [4]


It is assumed that the radiating elements are flat apertures, for example, end of open waveguide. The RP depends on the angular coordinates ϖand Ϥ. From the practical point of view, the consideration of such models of the arrays is justified by the fact that they allow to receive the values of required radiation characteristics with the accuracy of 2% - 5%, while the time of solution for the respective problems of analysis (determination of the RP of array), considerably decreases. Such approach to the solution of direct electrodynamic problems is effective especially for the arrays with constant coordinate surfaces, e.g., for the plane, cylindrical and conical arrays. Proceeding from the above assumptions, we separate the synthesis problem of such arrays into two synthesis problems for the linear and circular arrays.

The spatial RP of cylindrical array [9] is:


where Nis the number of radiators in circular subarray (identical for all subarrays), Mis the number of circular subarrays, ais radius of cylinder. The currents Inmare complex numbers, by means of which choice the approximation to the given amplitude RP F(ϯ,Ϥ)is carried out.

The functions fnm(ϯ,Ϥ)are the RPs of separate radiators, (zm,a,Ϥn)are the coordinates of radiators. The RP of separate radiators is identical and can be presented in the following form


Following the above assumptions, the RP (94) can be written as




is the RP of linear array, and


is the RP of circular array for each ϖ. Below we consider the synthesis problem for the circular array.

3.4.2. Optimization criteria

The complex currents in the array elements are determined by their amplitudes and phases. We denote these values |In|and Ϧn, respectively. The amplitudes |In|of currents are prescribed together with the amplitude RP Fin the problem of phase synthesis. The phases Ϧnare the optimizing parameters. We use the functionals (29) and (33) for optimization.

The equalities


should be satisfied at the points of functional (33) maximum. This set yields the system of transcendental equations for the phases Ϧof current and phase RP ϥ.

Using normalization of the current Ivalues: ||I||=1, we write down the functional (33) in two equivalent forms


where ϥ=argfis the phase RP. The operator A*is adjoint to Aand determined similarly to [Andriychuk et al., 1993].

At first, we consider the problem of Ϙmaximization. Substituting (99) into (100), we receive the system of nonlinear algebraic equations for optimal phase distribution of currents


In practice, the system (99), (100) is more convenient for the numerical solution. For this reason, we use the following iterative process


to find the phases Ϧ(k+1).

The gradient methods are more convenient to solve the minimization problem for the functional ϡ. The method of conjugated gradients [28] is the most suitable for this purpose. In this method the next approximation of phase vector Ϧ={Ϧn,n=1,...N}is calculated by the formula



z(k)is the gradient of ϡwith respect to the phases of currents in the "point" Ϧ(k). The components of vector z(k)are the following

Parameter ϒ(k)is determined from minimum of ϡbeing a function of this parameter.

In practice, it is necessary to solve the problem of discrete phase synthesis [9, 22], because the arbitrary phase distributions cannot be realized in the array radiators. These distributions are prescribed as a set of discrete values, which are multiple to the given phase discrete value β, that is Ϧn=ϙnβ, the integers ϙnare unknown in this case [5, 7].

The algorithm consisting of two enclosed iterative processes is used for the solution of this problem. The value of phase Ϧnis improved in the n-th step of the internal iterative process, the phases in other radiators remain fixed. At the same time, the mean-square deviation of the synthesized RP and function Fexp(iargf(k)), where argf(k)is the phase RP, which is received in the previous step of external iterative process, is minimized. The phase RP argfis improved in the external iterations.

The internal cycle consists of the successive improvement of phases in the separate radiators changing their number from 1up to N. The value of ϡdecreases in each step of internal cycle.

The new phase RP argf(k+1)is calculated by the found values of {Ϧ}in external cycle. The values of ϡcorresponding to new phase RP also decrease, what provides the convergence of the whole algorithm. In view of the step-type behavior of Ϧnvalues, this convergence exists not only for a sequence of ϡ, but also for the phase distributions. The iterative process is considered completed, if there is no change of Ϧnin the internal cycle.

The problem of discrete phase synthesis is solved in two steps. In the first step, the synthesis (with small accuracy) without the account of phase discrete values is carried out. After that, the found phases are approximated up to the nearest discrete values. In the second step, the described above algorithm of discrete synthesis is used. As a rule, one is enough to make several external cycles in the latter algorithm.

3.4.3. The results of numerical modeling

The numerical results are given for the sector array. The RPs of separate radiators have a cosine form, and mutual coupling of separate radiators is not taken into account [35]. The number of radiators N=32, the radiators are placed in active sector ϐ=90o.

The synthesis results are shown in Fig. 8 for the bi-directional RP


The thick line corresponds to the prescribed amplitude RP F, and thin line corresponds to the synthesized RP |f|. The prescribed Fand synthesized |f|RPs coincide in the main lobe up to level –20 dB, the level of side lobes does not exceed -20 dB (see Fig. 8a).

In the practical applications, the problem of phase scanning [35] is considered for arrays. The array alongside with creating the amplitude RP which is more close to the desired one should provide the moving this RP along the angular coordinate in the scanning process. This moving is carried out by the change of the phase distribution {Ϧn}only; the amplitudes {|In|}of current remain constant. In fact, the problem of phase synthesis is solved for each scanning angle Ϥsseparately.

The difference between the results of continuous and discrete synthesis depends on the value of phase discrete β. The smaller the difference, the smaller this value. In Fig. 8b the values of ϡfor two types of synthesis are shown for the process of scanning, Ϥsis changed in range from 10up to 90. Solid line corresponds to the case of continuous synthesis, and dashed line correspond to the case of discrete synthesis, the value of β=22.50. This value of βgives not big difference for two types of synthesis. So, the difference between Fand |f|in main lobe of RPs does not exceed 1dB, and this difference does not exceed 10dBin the side lobes. Such difference is satisfactory for the engineering practice. Of course, the above mentioned difference grows if the scanning angle Ϥsapproaches to the left or right border of the active sector ϐ.

Figure 8.

Synthesis of bi-directional RPF

3.5. Investigation of branching solution

The problem of the non-uniqueness of solutions for phase synthesis problem is investigated on the example of linear array.

The various modifications of the Newton method have been developed for solving the nonlinear equations in [13], and have been detailed for the synthesis problems in [11]. We consider here the above approach for determination of the number of solutions and investigation of its properties by the example of the nonlinear equation (46), corresponding to functional Ϙ(33).

Let the operator Adescribe the RP (array factor) of linear array [4].


where c=kasinϏ, ais the radius of array, Ϗis the angle in which the RP Fis non-zero, N=2M+1is a total number of array elements, Ϝis the generalized angular coordinate. The operator A*acts in the following way


In this case, the functional (33) can be presented as


3.5.1. The numerical results

The numerical calculations are carried out for the prescribed amplitude radiation pattern F(Ϝ)=cos(ϞϜ/2)and are shown in Figs. 9 and 10. In Fig. 9a, the values of Ϙand ϡare shown for various types of initial approximation of the current’s phase Ϧ(0)(x). For this case, the maximization problem of Ϙis equivalent in some sense to minimization problem of ϡ[Andriychuk et al., 1993].

The solid lines correspond to values of Ϙ, and the dashed lines correspond to values of ϡ. The number of array elements N=11, parameter cchanges from c=0up to c=2. For the values of Ncwhich do not exceed Nc=5all types of solutions give the same values of Ϙand ϡ. At Nc2Ϟthe branching of solutions appears, and optimal value for Ϙand ϡfunctionals gives the solution with even phase Ϧ. The prescribed amplitude radiation pattern Fand synthesized |f|are shown in Fig. 9b. The amplitudes |f|in the considerable extent differ from the amplitude Fbecause of small value of cparameter (c=1.6).

Figure 9.

The values of optimizing functionals and synthesized RPs for various type of solutions

The optimal values of sought phase distributions Ϧare shown in Fig. 10a, the given current amplitude distribution is |I|1. The optimal phase distributions Ϧkeep the parity properties of corresponding initial approximations Ϧ0. The optimal values of Ϙand ϡprovide the solution with phase distribution Ϧ0(x)=cos(x)(see Fig. 9).

The quality of approximation to prescribed amplitude pattern Ftoo much depends on the parameter c(see Fig. 10b). At c=3.14the Ϙand ϡvalues are noticeably smaller than for c=1.6, although the value of Nis larger in case of Fig. 9.

Figure 10.

The optimal phase distributions atN=11(a), and synthesis results atc=3.14(b)

3.6. Synthesis of waveguide array

The mutual coupling of the separate elements of array is taken into account in the process of solution of a direct problem (analysis problem) [7, 14].

3.6.1. Statement of synthesis problem

The objective functional is formulated as [4]


where N=2M+1is a number of exciting waveguides, F(ϖ,Ϥ)is the prescribed amplitude RP, |f(ϖ,Ϥ)|is the amplitude of synthesized RP, In(xn,yn)are the currents in the waveguide apertures. Geometry of waveguide array is shown in Fig. 11.

Figure 11.

Geometry of plane periodical waveguide array

The determination of the currents In(xn,yn)in the waveguide apertures (the solution of analysis problem) results in solution of the integral equation system [14].

The RP (array factor) [9] of array is:


where xn0and yn0are the coordinates of central points of apertures, anis the complex excitation coefficient for n-th waveguide.

Introducing the generalized angular coordinates Ϝ1and Ϝ2, we receive the expression for the RP




and finally


where a={a1,a2,...aN}, fn(Ϝ1,Ϝ2)are the RPs of separate waveguides.

The expression (116) indicates that the calculation of array factor f(Ϝ1,Ϝ2)using the excitation coefficients anis realized by the linear operator A. The coefficients anwill be the optimization parameters in the synthesis problem. Solving the synthesis problem, we minimize the functional ϡt(34).

3.6.2. The modeling results

The results of numerical calculations are presented for the waveguide arrays consisting of 15 and 31 radiators; kl=1.2, kL=18.75, kis wavenumber, land Lare width and length of waveguide aperture respectively. The prescribed amplitude RP is: F(ϖ,Ϥ)=F1(ϖ)F2(Ϥ), where


In Fig. 12, the dependence of the synthesis results on the value Nof waveguides is shown. It is easy to see, that the synthesized amplitude RP has narrower main lobe if Nincreases. The level of the first sidelobe is -30.26dB and -30.71dB respectively. The low level of sidelobes and velocity of its decrease is very important characteristic of the synthesized amplitude RPs. As rule, one requires the level of first sidelobe not greater than -20dB and not very slow decreasing the next sidelobes. The above mentioned characteristic for the synthesized amplitude RP in the plane Ϝ2Ozare shown in Table 2. The amplitude RP |f|at N=31has lower first sidelobe -30.71dB and faster decrease of the far sidelobes.


Table 2.

The level of sidelobes (in dB) corresponding to array with various N

Figure 12.

The synthesized amplitude RPs for variousN


4. Conclusions

In the process of solving the amplitude-phase synthesis problem, the influence of weight multiplier ton the synthesis results was investigated. It turned out that the mean-square deviation of the prescribed Fand synthesized |f|amplitude RPs diminishes if tdecreases. At the same time the norm ||I||of current grows. This testify that one should take into account the above fact in order to receive the solution with smaller mean-square deviation of the RPs or with small norm of current. The elaborated iterative procedures guarantee convergence of the successive approximation methods for numerical solving the arising non-linear equations.

The proposed approach for solving the synthesis problems of resonant antennas is universal, and it provides the possibility to synthesize antennas with the arbitrary form of external boundary. The calculation time of the RP of antenna is small enough what it is very important in the process of solution of the synthesis problem. The used variational statement of the synthesis problem also allows to take into account restrictions on the field at the given points (areas) of a near zone.

The synthesis of the resonant antenna with waveguide excitation gives the possibility to take into account the various requirements to the FRF in the operating frequency range. The developed algorithms enable to optimize the values of FRF in the single range, as well as in the several frequency ranges. The values of the objective parameters dand ϟ, which are received in the process of numerical calculations, are constructive characteristics of resonant antenna and they can be directly used in the antenna design.

The variational approach for solving the phase synthesis problem can be applied effectively for the plane, cylindrical and conical arrays. It allows to decrease the computational time, at the same time the accuracy of determination of the array characteristics is sufficient for practice. The branching solutions are investigated for the case of linear array. It is shown that one can receive the solutions with various properties starting the iterative process with different initial approximations.

The considered optimization problems of waveguide array give the possibility to take into account the requirements to the amplitude RP and amplitude-phase distribution of field in the aperture of exciting waveguides. The developed algorithms enable to achieve the minimal mean-square deviation ϡof the prescribed and synthesized amplitude RPs, and to optimize simultaneously the restrictions on the phase or amplitude characteristics of the excited fields.


  1. 1.AgranovichM. S.KatsenelenbaumB. Z.SivovA. N.VoitovichN. N.1999Synthesis of Antenna Systems According to the Desired Amplitude Radiation CharacteristicsBerlin, Wiley-VCH.
  2. 2.AlexidzeM. A.1978Solution of Boundary Problems by Decomposition on Non-orthogonal Functions, (in Russian). Moscow, Nauka.
  3. 3.AndriychukM. I.VoitovichN. N.1985Synthesis of a Closed Planar Antenna with Given Amplitude Pattern. Radio Eng. and Electron. Phys.,3053540
  4. 4.AndriychukM. I.VoitovichN. N.SavenkoP. A.TkachukV. P.1993The Antenna Synthesis According to Amplitude Radiation Pattern: Numerical Methods and Algorithms, (inRussian). Kiev, Nauk. Dumka.
  5. 5.AndriychukM. I.1998The Analytical-numerical Solution Method of the Nonlinear Problems of the Antenna Phase Synthesis. Proc. of IIIrd Intern. Seminar/Workshop onDirect and Inverse Problems of Electromagnetic and AcousticWave Theory (DIPED-98), Tbilisi,Georgia,8689
  6. 6.AndriychukM. I.ZamorskaO. F.2004The Antenna Systems Synthesis According to the Amplitude Characteristics under Condition of Electromagnetic Compatibility. 2004 Second Intern. Workshop on Ultrawideband and Ultrashort Impulse Signals, Sevastopol,Ukraine. Sept. 2004,135137
  7. 7.AndriychukM. I.KlakovychL. M.SavenkoP. O.TkachM. D.2005Numerical Solution of Nonlinear Synthesis Problems of Antenna Arrays with Regard for Coupling of Radiators. Proc. of Vth International Conference on Antenna Theory and Techniques, 24-27May 2005, Kyiv, Ukraine,213216
  8. 8.AndriychukM.ZamorskaO.2006Waveguide Antenna Synthesis According to the Amplitude Radiation Characteristics in the Frequency Band. Proc. of IstEuropeanConference on Antennas and Propagation.610November 2006. Nice, France, 6 p.
  9. 9.BalanisC. A.2005Synthesis of Antenna Systems According to the Desired Amplitude Radiation Characteristicsrd Edition. Hoboken, NJ, Wiley.
  10. 10.BityutskovV. I.2001Bunyakovskii Inequality, in M. Hazewinkel, Encyclopedia ofMathematics. Springer.
  11. 11.BulatsykO. O.KatsenelenbaumB. Z.TopolyukY. P.VoitovichN. N.2010Synthesis of Antenna Systems According to the Desired Amplitude Radiation CharacteristicsWeinheim, Wiley-VCH.
  12. 12.CupradzeV. D.1967About Approximate Solution of Mathematical Physics Problem, (in Russian). Uspekhi Mat. Nauk,22259107
  13. 13.DeuflhardP.2004Synthesis of Antenna Systems According to the Desired Amplitude Radiation CharacteristicsSeries Computational Mathematics, 35, Springer.
  14. 14.DmitrievV. I.BerezinaN. I.1986Numerical Methods of Solution of the Synthesis Problems for the Radiating Systems, (In Russian). Moscow, MGU.
  15. 15.FeldYa. N.1992On the Quadratic Lemma in Electrodynamics. Sov. Phys.- Dokl., 37,235236
  16. 16.FourikisN.2000Synthesis of Antenna Systems According to the Desired Amplitude Radiation CharacteristicsPhased Array Systems. Ascot park, Academic Press.
  17. 17.GelfandI. M.FominS. V.2000Calculus of Variations. Dover Publ.
  18. 18.KaklamaniD. I.AnastassiuH. T.2002Synthesis of Antenna Systems According to the Desired Amplitude Radiation CharacteristicsIEEE Antennas & Propagation Magazine,4464864
  19. 19.KatsenelenbaumB.MarcaderdelRio. L.PereyaslavetsM.atal. (1998). Theory of Nonuniform Waveguides, London, IEE Series.
  20. 20.KatsenelenbaumB. Z.2003Electromagnetic Fields- Restriction and Approximation. Weinheim, Wiley-VCH.
  21. 21.KatsenelenbaumB. Z.2006Synthesis of Antenna Systems According to the Desired Amplitude Radiation CharacteristicsWeinheim, Wiley-VCH.
  22. 22.KhzmalyanA. D.KondratyevA. S.2003Synthesis of Antenna Systems According to the Desired Amplitude Radiation CharacteristicsIEEE Trans. Antennas and Propag.,AP-512264272
  23. 23.KrasnoselskyM. A.ZabreikoP. P.1975Geometric Methods of Nonlinear Analysis, (in Russian). Moscow, Nauka.
  24. 24.KumarA.HristovH. D.1989Synthesis of Antenna Systems According to the Desired Amplitude Radiation CharacteristicsBoston- London, Artech House.
  25. 25.LandauL.LifshitzE.1984Synthesis of Antenna Systems According to the Desired Amplitude Radiation CharacteristicsLondon, Pergamon Press.
  26. 26.MarkovG. T.PetrovB. M.GrudinskayaG. P.1979Electrodynamics and Propagation of Waves, (in Russian). Moscow, Sov. Radio.
  27. 27.MilliganT. A.2005Synthesis of Antenna Systems According to the Desired Amplitude Radiation CharacteristicsHoboken, NJ, John Wiley & Sons.
  28. 28.PolakE.1971Computational Methods in Application. A Unified Approach. New York, Academic Press.
  29. 29.SavenkoP. O.2002Nonlinear Synthesis Problems of Radiating Systems. Lviv, Ukraine, IAPMM.
  30. 30.SmirnovV. I.1958Higher Mathematics, (in Russian).4Moscow, GIFML.
  31. 31.StyeskalH. l.HerdJ. S.1990Mutual Coupling Compensation in Small Array Antennas. IEEE Trans. on AP,381219711975
  32. 32.TikhonovA. N.ArseninV. Y.1977Synthesis of Antenna Systems According to the Desired Amplitude Radiation CharacteristicsNew York, Wiley.
  33. 33.VapnyarskiiI. B.2001Lagrange Multipliers, in M. Hazewinkel, Encyclopedia of Mathematics, Springer.
  34. 34.VeinsteinL. A.1998The Electromagnetic Waves, (in Russian). Moscow, Radio iSvyaz.
  35. 35.VendikO. G.ParnesM. D.2002Antennas with Electrical Scanning, (in Russian). Saint Petersburg, Science Press.
  36. 36.Vladimirov et al.1985Electromagnetic Compatibility of Radioelectronic Devices and Systems, (in Russian). Moscow, Radio iSvyaz.
  37. 37.VoitovichN. N.KatsenelenbaumB. Z.KorshunovaE. N.atal. (1989Electrodynamics of Antennas with Semitransparent Surfaces, (in Russian). Moscow, Nauka.
  38. 38.WeinbergM. M.TrenoginV. A.1969Theory of the Branching of Solutions of Nonlinear Equations, (in Russian). Moscow, Nauka.

Written By

Mykhaylo Andriychuk

Submitted: November 18th, 2011Reviewed: April 24th, 2012Published: September 19th, 2012