Optimization of Planar Antenna Arrays Using the Firefly Algorithm

Brazilian Microwave and Optoelectronics Society-SBMO received 30 July 2018; for review 2 Aug 2018; accepted 12 Dec 2018 Brazilian Society of Electromagnetism-SBMag © 2019 SBMO/SBMag ISSN 2179-1074 Abstract— This paper describes the application of the Firefly Algorithm for the optimization of linear and planar antenna arrays. In order to demonstrate the potential of this technique, three applications are reported. The first one is the synthesis of a radiation pattern with isoflux distribution with a non-uniformly spaced linear array composed of isotropic antennas. The second application is the optimization of a non-uniformly spaced planar array composed of isotropic radiators installed on an 8-U nanosatellite. Finally, the optimization of a 3×3 planar antenna array for beam steering with simultaneous side lobe level control is described, whereby the beamforming coefficients were allowed to assume only discrete values. In the three cases considered, good agreement between the desired (mask) and the optimized patterns has been obtained.


I. INTRODUCTION
The beamforming is an important topic in wireless communications, since realistic antennas provide spatial filtering, i.e. an antenna receives/radiates different power levels depending on the direction of signal incidence/radiation [1].In order to achieve beamforming, antenna arrays are applied normally in linear, planar, circular or spherical configurations.Thus, according to [1], in addition to the geometric arrangement, there are four additional design parameters: the distance between the array elements, the amplitude and the phase of the excitation coefficients, and the elemental radiation pattern.Due to the fact that the design parameters allow infinite combinations, sophisticated (or complex) patterns can only be synthesized by using a proper optimization approach.
In order to deal with this issue, many optimization methods can be used to find the best set of design parameters that generates a desired radiation pattern, also called mask.In most cases, heuristic methods are appropriate, such as the Genetic Algorithm (GA) [2]- [7], the Particle Swarm Optimization (PSO) [8]- [12], the Ant Colony Optimization (ACO) [13]- [17], and the Cat Swarm Optimization (CSO) [18] to cite only a few.There are also optimization methods that are based on the Calculation Theory, such as the Sequential Quadratic Programming (SQP) and the Gradient Method (GM) [19], [20].When optimization methods are applied in beamforming, they try to find optimal values for the design parameters in such a way as to minimize the deviation of the synthesized Optimization of Planar Antenna Arrays Using the Firefly Algorithm Eduardo Yoshimoto and Marcos V. T. Heckler Laboratório de Eletromagnetismo, Micro-Ondas e Antenas -LEMA Universidade Federal do Pampa (UNIPAMPA) -Campus Alegrete, Alegrete -RS, Brazil yoshi.edu@gmail.com;marcos.heckler@unipampa.edu.brradiation pattern from the desired pattern (mask).
In [21], two annular ring antenna arrays of isotropic elements were modeled in order to produce an isoflux radiation pattern for satellites operating in Geostationary Orbit (GEO) and Medium-Earth Orbit (MEO).Differential evolution was used for beamshaping and, at the same time, to minimize the sidelobe level (SLL) [22].The optimization yielded both arrays with diameter equal to 4 λ, which allowed low deviation from the desired main beam shape and SLL of -26.97 dB and -23.8 dB for the satellites operating in GEO and MEO, respectively.
In [23], an annular ring antenna array was also designed to produce an isoflux radiation pattern.The antenna array operates in X-band with circular polarization obtained by using sequential rotation of the elements and has been integrated onto one of the faces of a 1-U nanosatellite.The magnitudes and phases of the excitation coefficients and the radii of the two concentric rings were optimized using GA.The synthesized pattern showed deviation from the mask that has been attributed to the absence of mutual coupling in the model used for the synthesis.
Another heuristic applicable to pattern synthesis is the Firefly Algorithm (FA) optimization method, which was introduced by [24] in 2009.This method is based on the behavior of fireflies in searching a partner for mating, attracting preys or as a mechanism of protection against predators.The firefly that presents the best response to a particular goal generates bioluminescence with larger intensity, which attracts the rest of the colony.The attraction is weighted according to the distance between the fireflies and depends on the decaying nature of light propagation.In addition, in order to allow the fireflies to explore fully the search space, randomness must be included in the formulation.The efficiency of the method is demonstrated by [24], whereby optimizations of multi-dimensional mathematical functions and multiple optimal values have been demonstrated.The FA exhibited better performance in comparison to GA and PSO in 90 % of the analyzed cases.Other applications of the FA have been reported in [25]- [29].
Due to these advantages, this paper describes the application of FA for beamforming scenarios typical for space applications.The next section presents the mathematical background for the FA implementation.With this formulation, three different applications of the proposed technique are presented: the synthesis of a radiation pattern with isoflux distribution with a non-uniformly spaced linear array of isotropic antennas is described in section III, the optimization of a non-uniformly spaced planar array of isotropic radiators installed on an 8-U nanosatellite is presented in section IV, and the modelling of a 3×3 planar antenna array for beam steering with simultaneous sidelobe level control is discussed in section V.In this last case, the beamforming coefficients were allowed to assume only discrete values.
II. MATHEMATICAL BACKGROUND FOR THE FA IMPLEMENTATION Given the good performance of the FA in solving multidimensional problems, the equation presented by [24] is briefly described below, where the update of the position of the i-th non-optimal firefly towards the best located j-th firefly, i.e. the firefly who is supposed to be the closest to the optimum solution, is given by where t represents the iteration, β 0 is the attractiveness between the fireflies, which can be in the range [0, 1], γ is the light absorption coefficient, which can be in the range [0.01, 100], r ij is the line-of-sight distance between the fireflies, is the decaying rate due to the light propagation, α is the randomness applied to the i-th firefly in the range [0, 1] and ε i is a uniformly distributed random scalar value between -1 and 1.The parameters β 0 and γ can be used to control the fireflies searching capability.For instance, if large values (in the respective allowed ranges) are set to β 0 and γ, than the fireflies tend to focus on local search.On the other hand, the fireflies can be driven to perform global search by attributing a large value (in the allowed range) to the parameter α.
The term r ij is calculated by whered is the amount of variables to be optimized.The dynamics described by ( 1) can be summarized in Fig. 1.The function to be optimized (minimized) is the error function, which, in this work, has been defined as where mask E  is the electric field matrix referring to the desired mask and E  corresponds to the with where N x and N y correspond to the number of elements along the x and y directions, respectively, a mn and φ mn are the amplitudes and phases of the excitation coefficients that should weight the field components E θmn and E ϕmn radiated by the mn-th array element, k 0 is the propagation constant in free space, and x m and y n are the coordinates of the mn-th array element.According to [26], for minimization problems, the brightness of each firefly is inversely proportional to the value of the error function, whereas for maximization problems, the brightness is directly proportional.For the present case, the FA has been implemented in a computer code in MATLAB following the flowchart shown in Fig. 3.

III. FA APPLICATION FOR NON-UNIFORMLY SPACED LINEAR ANTENNA ARRAY COMPOSED OF ISOTROPIC ELEMENTS
The first proposed application for the FA considers a linear antenna array with 12 isotropic elements.The elements are arranged along the x axis according to Fig. 4. The positions of the elements have been adjusted to the dimensions of an 8-U nanosatellite, which is shown in Fig. 5 flying over the Earth.The cubic structure, as well as the open flaps have edge sides of 20 cm.The radiation pattern to be synthesized should follow an isoflux distribution, which is plotted as a black dotted curve in Fig. 6, so as to allow illuminating the Earth surface with constant power distribution.
The main pattern constraints are shown in Fig. 5, with SLL ≤ -15 dB below the maximum of the radiation pattern.To start the FA optimization, 15 fireflies were randomly allocated in the search space.The parameters β 0 , α and γ were varied during the optimization process as given in Table I.
The resulting excitation coefficients are listed in Table II and the optimized radiation pattern is shown in Fig. 6.It is notorious that the synthesized radiation pattern satisfied the mask with irrelevant deviations around θ = ±45°.The regions without mask tracing are not considered in the calculation of (3).The value of the error function at the beginning of the optimization was 7.32 and, after 1000 iterations, the optimization was terminated with error = 0.18, yielding a reduction to 2.46 % of the initial value.

IV. FA APPLICATION FOR A PLANAR ANTENNA ARRAY OF ISOTROPIC ELEMENTS INSTALLED ONTO AN 8-U NANOSATELLITE
For the scenario depicted in Fig. 5, the antenna should present a nearly rotation symmetric pattern following the contour defined in Fig. 6 for any ϕ.This results in the two-dimensional mask for the upper hemisphere of the array (0°≤ϕ≤ 360° and 0°≤θ≤ 90°) as shown in Fig. 7.In order to fulfill such a mask, a planar antenna array is needed.In a first scenario, the antennas could be placed onto the open flaps, and this is the case considered in this section.A second scenario would be to consider the antennas to be integrated onto the bottom surface of the 8-U cubic structure, which is the case to be analyzed in the next section.The antenna array topology considered in this case is shown in Fig. 8.The nanosatellite structure itself and the design frequency (2.26 GHz) yield an edge-to-edge distance between subarrays equal to In Fig. 9, the angular aperture between the grating lobes is close to the angular aperture required for the isoflux distribution (Fig. 7).Thus, the FA was applied to the synthesis of an isoflux radiation pattern with SLL ≤ -15 dB.For this case, the FA has been started with 10 fireflies with random values for excitation coefficients (magnitude and phase), whilst attractiveness, randomness and light absorption coefficients are varied according to Table III during the FA run.The synthesized isoflux radiation pattern can be seen in Figs. 10 and 11 as 3D polar plot with color scale and 2D rectangular plot for ϕ = 90°, respectively.The shape of the synthesized isoflux radiation pattern was obtained successfully with deviations from the mask only in small regions of the main lobe.The value of the error function at the beginning of the optimization was 619.9 and could be reduced to 47.5 (7.66 % of its initial value).These results confirm the good convergence behavior of FA in the optimization of multidimensional problems (128 variables in the present case).

COEFFICIENTS
In this case, a 3×3 planar array composed of realistic patch antennas with right-handed circular polarization (RHCP) and dimensions compatible with the bottom face of an 8-U CubeSat (200 mm of edge size) was considered.The geometry of the array element is based on the antenna designed in [31] and is detailed in Fig. 12.The geometry of the antenna array is shown in Fig. 13.
This array should be able to steer the main beam 45° off from boresight and with SLL ≤ -20 dB.For such an array, the worst case is normally achieved when beamsteering is done to ϕ max = θ max = 45°, whereby large sidelobes tend to appear.Since the array geometry is implemented now in microstrip technology, mutual coupling must be considered.This has been achieved by modelling the array in ANSYS HFSS and by extracting the embedded pattern of each element.The FA has been configured to work with 50 fireflies that have been randomly placed in the search space.The other parameters are listed in Table IV.In addition, the algorithm was configured to perform the optimization only for a discrete set of the excitation coefficients with 4-bit representation for both magnitude and phase.Thus, the magnitude can assume values between 0-15 dB, with steps of 1 dB, and the phase can range from 0° to 337.5°, with steps of 22.5°.The gains in the desired directions and the achieved SLL are summarized in Table IX.The failure to comply with the direction of maximum gain and the constraint SLL ≤ -20 dB can be explained due to the discretization of the excitation coefficients.Better performance seems to be only possible if the number of array elements or the number of bits for representation of the excitation coefficients is increased.However, this is not feasible for the proposed application due to size, weight and system complexity restrictions of the nanosatellite under consideration.
To assess how the beam steering affects the reflection coefficients of the array elements, Fig. 18 shows the variation of the reflection coefficients Γ (m,n) for the four patterns shown above.The elements (1,1), (2,1) and (2,2)stand for antennas located on the corner, edge and on the middle of the array, respectively.Based on the plotted curves, the beamsteering seems to have little impact on the reflection coefficients.

VI. CONCLUSION
In this paper, the use of an optimization code based on the Firefly Algorithm (FA) was demonstrated for the synthesis of a non-uniformly spaced linear antenna array, a non-uniformly spaced planar array and for a uniformly spaced planar microstrip array.Initially, the mathematical background needed to understand and implement the FA was presented.As main results, the synthesized radiation patterns exhibited good agreement with the desired masks for the three cases considered.The main contribution of this paper is to demonstrate the efficiency of the FA for the optimization of complex radiation patterns, which have been synthesized for a variety of array geometries.

Fig. 3 .
Fig.3.Flowchart of the computer code based on the FA.

Fig. 4 .
Fig.4.Positions of the elements on the linear antenna array along the x axis (design frequency equal to 2.26 GHz, d x = 0.3λ 0 and d a = 0.6λ 0 ).

Fig. 6 .
Fig.6.Desired and synthesized isoflux radiation pattern for the scenario shown in Fig. 5.

Fig. 7 .
Fig.7.Isoflux radiation pattern for the characteristics shown in Fig. 5.The angular scale corresponds to the ϕ coordinate of the spherical coordinate system and the radial scale to the θ coordinate.The color scale indicates the desired electric field intensity.

Fig. 9 .
Fig.9.Radiation pattern produced by the antenna array of Fig. 8 for broadside radiation without SLL control.

Fig. 10 .Fig. 11 .
Fig.10.Isoflux radiation pattern produced by the antenna array of Fig. 8 with the excitation coefficients found by the FA optimization method.

Fig. 13 .
Fig.13.Right-handed circularly polarized 3×3 planar microstrip array (L array = 199.5 mm, the distance of the center-to-center of elements is equal to 66.5 mm).The figures between parenthesis stand for the identification of the mn-th array element.
technique in a realistic scenario, the FA has been applied to determine the beamforming coefficients to steer the main beam to the following directions (SLL ≤ -20 dB): θ max = 0° and ϕ max = 0°, θ max = 15° and ϕ max = 45°, θ max = 30° and ϕ max = 45°, and θ max = 45° and ϕ max = 45°.The discrete excitation coefficients found by the FA are listed in Tables V-VIII.These coefficients produce the radiation patterns shown in Figs.14-17, which have been calculated with ANSYS HFSS.

TABLE II .
EXCITATION COEFFICIENTS OBTAINED WITH FA FOR THE LINEAR ANTENNA ARRAY.

TABLE III .
FA PARAMETERS FOR THE OPTIMIZATION OF THE PLANAR ANTENNA ARRAY.

TABLE IV .
FA PARAMETERS FOR THE OPTIMIZATION OF THE 3×3 PLANAR MICROSTRIP ARRAY.

TABLE V .
EXCITATION COEFFICIENTS OPTIMIZED TO POINT THE MAIN BEAM TO θ max = 0° AND ϕ max = 0° WITH SLL ≤ -20 dB.MAGNITUDES ARE GIVEN IN dB.

TABLE VI .
EXCITATION COEFFICIENTS OPTIMIZED TO POINT THE MAIN BEAM TO θ max = 15° AND ϕ max = 45° WITH SLL ≤ -20 dB.MAGNITUDES ARE GIVEN IN dB.

TABLE VII .
EXCITATION COEFFICIENTS OPTIMIZED TO POINT THE MAIN BEAM TO θ max = 30° AND ϕ max = 45° WITH SLL ≤ -20 dB.MAGNITUDES ARE GIVEN IN dB.

TABLE VIII .
EXCITATION COEFFICIENTS OPTIMIZED TO POINT THE MAIN BEAM TO θ max = 45° AND ϕ max = 45° WITH SLL ≤ -20 dB.MAGNITUDES ARE GIVEN IN dB.

TABLE IX .
COMPARISON BETWEEN THE ACHIEVED GAINS AND SIDELOBE LEVELS.