International Journal of Antennas and Propagation
Volume 2016 (2016), Article ID 7087298, 8 pages
http://dx.doi.org/10.1155/2016/7087298

Research Article

Iterative GA Optimization Scheme for Synthesis of Radiation Pattern of Linear Array Antenna

Sarayoot Todnatee and Chuwong Phongcharoenpanich

Faculty of Engineering, King Mongkut’s Institute of Technology Ladkrabang, Bangkok 10520, Thailand

Received 16 March 2016; Revised 3 June 2016; Accepted 8 June 2016

Academic Editor: Sotirios K. Goudos

Copyright © 2016 Sarayoot Todnatee and Chuwong Phongcharoenpanich. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This research has proposed the iterative genetic algorithm (GA) optimization scheme to synthesize the radiation pattern of an aperiodic (nonuniform) linear array antenna. The aim of the iterative optimization is to achieve a radiation pattern with a side lobe level (SLL) of ≤−20 dB. In the optimization, the proposed scheme iteratively optimizes the array range (spacing) and the number of array elements, whereby the array element with the lowest absolute complex weight coefficient is first removed and then the second lowest and so on. The removal (the element reduction) is terminated once the SLL is greater than −20 dB (>−20 dB) and the elemental increment mechanism is triggered. The results indicate that the proposed iterative GA optimization scheme is applicable to the nonuniform linear array antenna and also is capable of synthesizing the radiation pattern with SLL ≤ −20 dB.

1. Introduction

In addition to the fabrication simplification, an array antenna with minimal number of array elements is capable of generating a desired radiation pattern for unequally spaced and nonuniform excitation strategies. The array antenna with minimal array elements could also achieve high directivity, narrow beamwidth, and low side lobe level (SLL). The reduced number of array elements affords this kind of array antenna a wider variety of applications, in particular where the antenna dimension, weight, and the budget availability play a significant role, for example, in the radar systems, satellite communications, and mobile communications.

Regarding research on the radiation pattern synthesis techniques for uniformly spaced and equally distributed excitation, Dolph-Tschebyscheff and Taylor methods were utilized to derive the antenna array excitation coefficients for the radiation pattern synthesis. The proposed technique however resulted in a large, bulky antenna [1]. On the other hand, in [2, 3], an array antenna with nonuniform range (spacing) between antenna elements was deployed to improve the radiation pattern synthesis by which the reduction in the spatial aperture was achieved. Furthermore, in [4, 5], adjustments were made with the weight coefficients and the array ranges were to successfully improve the radiation pattern synthesis performance of the antenna.

Moreover, the noniterative/iterative technique has been implemented to enhance the radiation pattern synthesis of the nonuniform array antennas. In [6, 7], a noniterative algorithmic scheme based on the matrix pencil method (MPM) was utilized with the nonuniform linear array antenna to reduce the computational time. On the iterative side, both the local and global optimization algorithms have been utilized to optimize the weight coefficients and array range. With the local iterative optimization, the search time is shorter and the resources requirement lowers vis-à-vis the global optimization. The local optimization technique nevertheless requires that a suitable starting point be identified and, in a number of occasions, suffers from the suboptimal solution. Examples of the local optimization techniques are the steepest descent, quasi-Newton, and conjugate gradient techniques [8].

In regard to the global iterative optimization, the evolutionary algorithms (EAs) [9] were utilized to optimize the weight coefficients and array ranges for synthesis of the desired radiation pattern of the linear array antenna. In [10, 11], the simulated annealing (SA) algorithm was applied to the thinned array antenna to achieve the radiation pattern with low SLL. In [12, 13], the particle swarm algorithm (PSO) technique was applied to unequally spaced linear array for the radiation pattern synthesis with suppressed SLL and null control. In [14, 15], the genetic algorithm (GA) was employed to achieve the radiation pattern with low SLL, while in [16] the improved GA, in which the number of array elements was reduced, was proposed and subsequently applied to a linear aperiodic (nonuniform) array antenna to achieve a shortened computational time and low SLL radiation pattern.

Further research has been carried out to achieve the radiation pattern with low SLL and minimal number of array elements required. In [17], a linear sparse array antenna with the suppressed SLL radiation pattern and reduced number of array elements was presented. In determining the optimally minimal number of array elements, the algorithm iteratively optimizes the array range and weight coefficients. Instead of removing the array element in an individual fashion, their proposed method removes in one iteration with a multiple of array elements. In [18, 19], the GA and numerical Moment Method (MM) techniques were applied to successfully synthesize the radiation pattern with low SLL using the reduced number of array elements.

In order to optimize problems concurred in attaining the second objective for sparse nonuniform array synthesis, an iterative algorithmic scheme is proposed in this research, where the number of array elements is minimized in an iterative fashion. The ranges of array element are accordingly designed by using GA. The achievement of minimal array elements is removed by iterative algorithmic scheme being considered to have minimum contribution to weight coefficients. The fitness function (objective function) is used by the only term at minimizing the SLL. This research has proposed an iterative GA optimization scheme to synthesize the suppressed SLL radiation pattern (≤−20 dB) of an aperiodic (nonuniform) linear array antenna. In operation, the scheme iteratively optimizes the array range and the number of array elements, whereby the array element with the lowest absolute complex weight coefficient is first removed and then the second lowest and so on. The removal (the element reduction) is terminated once the SLL is greater than −20 dB (>−20 dB) and the elemental increment mechanism is triggered.

2. Method on Array Synthesis

In the synthesis of radiation pattern of an array antenna, the first step is to determine the array factor (AF). The focus of this research is on the nonuniform and asymmetrical linear array antenna; this section thus assumes a nonuniform and asymmetrical linear array with array elements along the and directions. The AF of the linear sparse array antenna can be expressed as

where

where is the number of array elements along the and directions; and are the positions of each element along the plane, as shown in Figure 1; is the complex weight coefficient of the th element, where and are the amplitude and phase; is the wavelength of the arrival signal incident wave at an azimuth of and elevation of ; and and are the direction cosines given by

The AF from (1) can be expressed in the matrix form as

where and are, respectively, the complex weight coefficient and the position of element matrix variable. To calculate , the matrix form in (4) is rewritten as in

where

The matrix of function in the exponential form can be expressed as

Figure 1: Geometry of nonuniform excitation coefficients and array range of a linear array antenna with elements.

The definition of array synthesis is to design the parameters of the array, that is, the ranges and the number of array elements. The squared response error (SQR) is performed on a pattern close to the desired radiation pattern. The squared response error (SQR) between the array factors of the synthesized radiation pattern () and reference radiation pattern () is then calculated. The SQR is subsequently substituted in the fitness function (fit) of the optimization process, which is later discussed. The squared response error (SQR) can be expressed as

where and , respectively, are the array factors of the synthesized and reference radiation patterns.

For the optimization process to approach the optimal solution, an iterative scheme is applied to the SQR to improve the fit process. The optimal solution is then algorithmically deployed to optimize the element range and the number of array elements. A typical SQR encompasses the main lobe (ML) and side lobe (SL) regions. Specifically, the iteration of SQR in the ML region is carried out through to the th-iteration, which is expressed as

where is the updated value of array structure from the optimization process and

where is the total number of sample points in the ML region.

The iteration of SQR in the SL region through to the th-iteration, in which the difference between and centers around zero, can be expressed as

where is the updated value of array structure from the optimization process and

where is the total number of sample points in the SL region.

As the fitness function (fit) of optimization is to minimize the SLL of the radiation pattern, the fitness function can be expressed as

where

To simplify this optimization problem, only the minimization of SLL is considered in the optimization. The ML width is fixed to be within a given range according to the design specifications.

3. Optimization Problem

In this research, the optimization aim is to reduce the number of array elements by iteratively optimizing the array range for the corresponding minimum weight coefficient () to ultimately achieve the radiation pattern with low SLL. The iteration is carried out using the proposed iterative GA optimization scheme until arriving at the optimally minimal number of array elements at which the synthesized radiation pattern with SLL −20 dB is realized. Thus, the optimization problem can be expressed as

where is the absolute minimum weight coefficient and and are the design specifications and constraints, respectively, of the main lobe (ML) and side lobe (SL).

3.1. Selective Weight Coefficient with Minimized Value

Prior to the reduction of array elements, the weight coefficients () of the array elements are calculated using (18) and, in the optimization, the element with the lowest absolute complex weight coefficient () is first removed and then the second lowest and so on. The removal process is terminated once the SLL is greater than −20 dB (>−20 dB) at which the elemental increment mechanism is triggered. The weight coefficient can be expressed as

where is equivalent to (7) and is the matrix of Generalized Method of Moments (GMM), which corresponds to [19], is the set of direction cosines, and th is the number of elements along the direction. The minimum of the weight coefficient can be calculated by

where

where and are, respectively, the th and th positions of the weight coefficient.

3.2. Array Ranges Reconfiguration

In the reduction of array elements, the proposed GA optimization scheme is implemented in an iterative fashion with a new array structure containing the reduced number of array elements and the reconfigured array range resulting in response to each successive iteration. The number of array elements is iteratively reduced to the optimally minimal number of array elements, subject to (16) and (17). In addition, the reduced array element number contributes to a reconfiguration array range size. The array range results of the reconfiguration fewer elements can be expressed as

where is the th array element, is the th array element of weight coefficient removed, and is the range constant of element positions.

3.3. Increment Array Element

In this research, the increment of array elements is undertaken once the SLL is greater than −20 dB (>−20 dB), by which an additional number of elements are reintroduced into the array. The increment would rectify and restore the target SLL. The array range between the (increment array) and th (current array) elements, which is the middle point of the maximum and the minimum boundary value of array ranges, can be calculated by

where is the th array element of . and are, respectively, the maximum and minimum boundary value of array ranges.

4. Iterative GA Optimization Scheme

This section discusses the iterative GA optimization scheme to identify the optimal array range and number of array elements that generate the radiation pattern with the target suppressed SLL. In the optimization, the proposed algorithmic scheme iteratively locates the optimal array range and minimal array elements which ultimately leads to the low SLL radiation pattern (≤−20 dB). In determination of the minimal number of array elements, the element with the lowest absolute weight coefficient () is first removed and then that of the second lowest and so on. The removal process ceases once the SLL is greater than −20 dB (>−20 dB), which subsequently prompts the increment of array elements.

The major steps of the proposed method are summarized below.

Step 1 (initialization). The parameter value of GA, that is, the population size, crossover, and mutation, must be specified. The initial parameters values of linear array and SLL are the target. In relation to initialization of the array ranges and population size, the array range (array spacing) relative to the population size can be expressed as

where and are, respectively, the maximum and minimum boundary values and is normally distributed random numbers, where is the th population size and is the th array elements.

Step 2 (the algorithmic scheme). Specifically, the iterative optimization starts with the array range of interval [, ] and array elements, both of which are substituted in (23). The array range is varied between and . The optimization goal is to identify a combination of array range and array element number that results in the radiation pattern with SLL ≤ −20 dB. The array element reduction and increment are sequentially implemented subject to the SLL target (≤−20 dB) to arrive at the optimal combination of array range and array element number. The optimization process is removed on array elements (target SLL ≤ −20 dB); the element with the lowest absolute weight coefficient is first removed and then that of the second lowest of one iteration (). Note that, for the first iteration process, the optimization scheme has (previous parameter). Consequently, the array element is reconfigured by (21) of the new ranges array elements. If the solution is unsatisfied on fewer array elements (SLL > −20 dB), the increment would rectify and restore the target optimization. The increment of array element is inserted by and the initial array range is calculated by (22).

Step 3 (the convergence of iteration method). The optimization process is terminated on condition of the following in the order given:

(1)Satisfying solution is found; that is, the optimal array elements and SLL have been obtained.(2)The fitness function during successive generation is improved on unsatisfied solution.

For the first condition, the results obtained in the current parameter of the optimization process are certified as the final solution. For the second condition, the results are obtained in the previous parameter as the final design.

5. Numerical Results

To validate its effectiveness, the proposed iterative GA optimization scheme is simulated using the population size, the crossover and mutation parameters, and the maximum iteration of 100, 0.8, 0.2, and 100, respectively. In addition, the minimum () and maximum () array ranges (spacing) are 0.5 and 1.0. The direction cosines of the and are sets and . The total number of sample points is set to be 1024, by which the number of ML and SL sample points is confined to 960 and 244, respectively.

5.1. Optimal Element Number with Complex Weight Coefficients

In this scenario, the initial parameters, as tabulated in Table 1, are deployed to generate the target radiation pattern with the direction confined to 0° with the ML region fixed and SLL of 10° and −20 dB. The iteration begins with the arbitrarily elemental number of 25 elements and the maximum iteration of 8.

Table 1: Initial parameters for the synthesized target array pattern.

Parameter
Value

Number of elements
25

ML direction ()
0°

ML region ()
5°, 5°

SL region ()
90°, −5°5°, 90°

Desired SLL (dB)
−20

Number of samples NML, NSL
960, 244

Figure 2 depicts the SLL of each iteration run for the maximum iteration of 8. The optimal array number for the target radiation pattern (SLL ≤ −20 dB) is 14 elements, corresponding to the 8th iteration. Interestingly, at the 7th iteration the SLL is greater than −20 dB (>−20 dB) and thereby induces the elemental reduction mechanism to cease and at the same time triggers the incremental mechanism to restore the SLL to the target (≤−20 dB). By comparison, the best SLL is −20.98 dB achieved in the 1st iteration with 25 array elements while the worst SLL is −18.47 dB (the 7th iteration and 13 array elements). Despite a poorer SLL of the 8th iteration (−20.06 dB) relative to the 1st iteration (−20.98 dB), the former is more operationally favorable due to a lower number of array elements required (i.e., 14 versus 25 elements). For the array range (spacing), the GA optimization iteratively optimizes the array spacing given the SLL target of ≤−20 dB. The optimal array range is achieved at 323 generations of the complex weight, as presented in Table 2.

Table 2: Results of iterative process.

Number of iterations
Number of generations
SLL (dB)

Complex weight
Real weight
Complex weight
Real weight

1st
13
23

2nd
14
10

3rd
27
9

4th
19
7

5th
14
22

6th
42
100

7th
100
9

8th
54
—

—

Figure 2: Number of iterations in relation to SLL with the target SLL of ≤−20 dB.

Figure 3 illustrates the synthesized radiation pattern with the optimally minimal array elements of 14 elements. In the figure, the main lobe (ML) is direction confined to 0°() and the synthesized radiation pattern using the proposed iterative GA optimization scheme exhibits the ML region fixed and SLL of 10°() and −20 dB, respectively. The results indicate that the proposed iterative algorithmic scheme is applicable to the linear array antenna for generation of the radiation pattern with SLL ≤ −20 dB.

Figure 3: The synthesized radiation pattern with optimally minimal array elements (14 elements) under complex weight coefficient scenario.

5.2. Optimal Element Number with Real Weight Coefficients

In this setting, the initial parameters to synthesize the radiation pattern with the direction confined to 0°() with the ML region fixed and SLL of 10°() and −20 dB are identical to the previous scenario (Table 1). The initial number of array elements is 25 elements and the maximum iteration is 7. In the optimization, the initial real weight coefficient interval () of the proposed iterative GA scheme is between 0 and 1 () and . In this case, the proposed method simultaneously optimizes the real weight coefficients and array range.

Figure 4 portrays the SLL of individual iteration run for the maximum iteration of 7. The optimal array number for the target radiation pattern synthesis (SLL ≤ −20 dB) is 16 elements, corresponding to the 7th iteration. In the figure, at the 6th iteration the SLL is greater than −20 dB (>−20 dB), prompting the elemental reduction mechanism to terminate and concurrently triggering the incremental mechanism. By comparison, the best SLL is −20.80 dB achieved in the 1st iteration with 25 array elements while the worst SLL is −19.48 dB (the 6th iteration and 15 array elements). Despite a poorer SLL of the 7th iteration (−20.67 dB) relative to the 1st iteration’s (−20.80 dB), the former is more operationally favorable due to a lower number of array elements required (i.e., 16 versus 25 elements). For the array range (spacing), the GA optimization iteratively optimizes the array spacing given the SLL target of ≤20 dB. The optimal array range is achieved at 180 generations of the real weight coefficient (Table 2).

Figure 4: Number of iterations in relation to SLL with the real weight coefficient interval and SLL of and <−20 dB.

Figure 5 depicts the synthesized radiation pattern under the real weight coefficient condition with the optimally minimal array elements of 16 elements. It can be observed that the main lobe (ML) of the radiation pattern is confined to 0°() and that its ML region fixed and SLL are, respectively, 10°() and −20 dB. The findings reveal that the proposed algorithmic scheme under the real weight coefficient condition is also applicable to the linear array antenna for synthesis of the radiation pattern with SLL ≤ −20 dB.

Figure 5: The synthesized radiation pattern with optimally minimal array elements (16 elements) under real weight coefficient scenario.

5.3. Performance Comparisons and Discussions

The performance comparison of the proposed method with the previous work is described in this subsection. The iterative method to reduce the number of elements used has been presented in [2, 17]. In [17], the linear array was synthesized with the number of elements of 25 over the 50 of spatial aperture. The beam pattern is in region of , where = , and the SLL range is spanned between and . The process starts with 25 elements and produces a resultant array with 19 elements with an overall spatial aperture of 21.09. The obtained SLL is −14.49 dB, and the width of the main lobe measured at −3 dB is 0.0214.

The performance is compared with above achievements by method used in [17]. The proposed optimization method is set to find an array element with minimal number (<19 elements) and SLL is set to lower than −14.49 dB for the condition to reduce the array element configuration. In addition, for the fair performance comparison, the real amplitude excitation coefficients are set between and the element spacing of not less than 0.5. These parameters are employed in the simulation process corresponding to [17]. The sample of 1024 points is on region () of beam pattern computation.

For the computational complexity of evolutionary algorithms, the array factor is employed for computation in terms of number of fitness functions. The sampling point of array factor is spanned with the interval of 1024 samples. Equation (13) is used to calculate the fitness function evaluated for 1024 times for each possible design.

The initialization is done with nonuniform linear array, the number of elements is started from 25 elements, and the spacing interval is in the range . The optimization process is checked for the convergence conditions for the best individual with 10 generations. The iterations of the proposed method are run for the reduction process of the 4th iteration (6, 17, 22, and 39) and the 1st iteration run (99) for the compensation. From the comparison results, the proposed iterative scheme has minimum number of elements (18 elements) and is followed by method in [17] (19 elements). In addition, the achievement of the lowest number of array elements is employed with 137 iterations’ run (18 elements), in comparison with 110 iterations’ run (19 elements) for the method in [17]. The number of iterations for evaluated fitness function is shown in Figure 6. The comparison results of the iterative process are summarized in Table 3.

Table 3: Comparison results of iterative process.

Number of iterations
Number of generations
Number of elements
Minimum SLL (dB)

Proposed method
Method [17]
Proposed method
Method [17]
Proposed method
Method [17]

1st
6
70
25
25

2nd
17
80
23
20

3rd
22
110
21
19

4th
39
—
19
—

—

5th
99
—
17
—

—

6th
137
—
18
—

—

Figure 6: The performance comparison of computation complexity in terms of number of fitness functions.

Table 4 compares the spatial aperture, minimum array elements, and minimum SLL belonging to method [2], method [17], and the proposed iterative GA method. By comparison, the proposed iterative scheme has the smallest aperture size (19.15), followed by method [2] (25) and method [17] (21.09). In addition, the proposed algorithmic scheme could achieve the lowest number of array elements (18 elements) and minimum SLL of −14.51 dB, in comparison with 24 and 19 elements, and −14.15 dB and −14.49 dB for method [2] and method [17], respectively.

Table 4: Comparison between method [2], method [17], and proposed iterative GA methods with regard to minimum array elements and minimized SLL.

Specifics
Methodology

Proposed
SA [2]
IGA [17]

Spatial aperture
19.15
25
20.09

Minimum SLL
−14.51 dB
−14.45 dB
−14.49 dB

Min. number of elements
18
24
19

For further comparison of the previous optimization methods, the optimization problem is concurred in attaining the second objective for sparse nonuniform array synthesis. The optimization design problem is to determine the optimum linear array where the SLL is minimized and the beamwidth is minimal. In the comparative consideration with method [20] (NSGA-II), both algorithms are set with a population size of 30 and run for 100 iterations. For NSGA-II [20] the crossover and mutation probabilities are set equal to 0.9 and 0.1, respectively. The obtained results of the proposed method and method in [20] are done for the number of elements of 25. In Figure 7, the obtained beamwidth is 2.7273 and 2.9032 degrees for both algorithms. The comparison of the array factor pattern, minimum beamwidth, and minimum SLL of the proposed method and method in [20] is shown in Figure 8. From the comparison results, the proposed method possesses the obtained beamwidth of 2.7273 degrees and the SLL of −17.45 dB. For the method in [20], the beamwidth is 2.7273, and the SLL is −17.41 dB, respectively.

Figure 7: The beamwidth comparison of linear array with 25 elements.

Figure 8: The comparison of array factor pattern.

6. Conclusions

This research has proposed the iterative GA optimization scheme for nonuniform linear array antennas to synthesize the radiation pattern with SLL ≤ −20 dB. In the optimization, the proposed GA scheme iteratively optimizes the array range and the number of array elements, whereby the element with the lowest absolute complex weight coefficient is first removed and then the second lowest and so on. The removal is terminated once the SLL is greater than −20 dB (>−20 dB) and subsequently the elemental increment mechanism is triggered. The results indicate that the proposed iterative GA optimization scheme is capable of effectively synthesizing the radiation pattern with SLL ≤ −20 dB and is applicable to the nonuniform (aperiodic) linear array antennas.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

References

C. A. Balanis, Antenna Theory: Analysis and Design, John Wiley & Sons, New York, NY, USA, 3rd edition, 2005.

A. Trucco and V. Murino, “Stochastic optimization of linear sparse arrays,” IEEE Journal of Oceanic Engineering, vol. 24, no. 3, pp. 291–298, 1999.

Z. Lin, W. Jia, M. Yao, and L. Hao, “Synthesis of sparse linear arrays using vector mapping and simultaneous perturbation stochastic approximation,” IEEE Antennas and Wireless Propagation Letters, vol. 11, pp. 220–223, 2012.

S. K. Goudos, K. Siakavara, T. Samaras, E. E. Vafiadis, and J. N. Sahalos, “Sparse linear array synthesis with multiple constraints using differential evolution with strategy adaptation,” IEEE Antennas and Wireless Propagation Letters, vol. 10, pp. 670–673, 2011.

R. C. Nongpiur and D. J. Shpak, “Synthesis of linear and planar arrays with minimum element selection,” IEEE Transactions on Signal Processing, vol. 62, no. 20, pp. 5398–5410, 2014.

Y. Liu, Z. Nie, and Q. H. Liu, “Reducing the number of elements in a linear antenna array by the matrix pencil method,” IEEE Transactions on Antennas and Propagation, vol. 56, no. 9, pp. 2955–2962, 2008.

Y. Liu, Q. H. Liu, and Z. Nie, “Reducing the number of elements in the synthesis of shaped-beam patterns by the forward-backward matrix pencil method,” IEEE Transactions on Antennas and Propagation, vol. 58, no. 2, pp. 604–608, 2010.

N. Damavandi and S. Safavi-Naeini, “A hybrid evolutionary programming method for circuit optimization,” IEEE Transactions on Circuits and Systems I: Regular Papers, vol. 52, no. 5, pp. 902–910, 2005.

C. Jacob, Illustrating Evolutionary Computation with Mathematica, Morgan Kaufann, 2001.

V. Murino, A. Trucco, and C. S. Regazzoni, “Synthesis of unequally spaced arrays by simulated annealing,” IEEE Transactions on Signal Processing, vol. 44, no. 1, pp. 119–123, 1996.

C. A. Meijer, “Simulated annealing in the design of thinned arrays having low sidelobe levels,” in Proceedings of the IEEE South African Symposium on Communications and Signal Processing (COMSIG '98), pp. 361–366, IEEE, Rondebosch, South Africa, September 1998.

M. M. Khodier and C. G. Christodoulou, “Linear array geometry synthesis with minimum sidelobe level and null control using particle swarm optimization,” IEEE Transactions on Antennas and Propagation, vol. 53, no. 8, pp. 2674–2679, 2005.

S. K. Goudos, V. Moysiadou, T. Samaras, K. Siakavara, and J. N. Sahalos, “Application of a comprehensive learning particle swarm optimizer to unequally spaced linear array synthesis with sidelobe level suppression and null control,” IEEE Antennas and Wireless Propagation Letters, vol. 9, pp. 125–129, 2010.

K. F. Man, K. S. Tang, and S. Kwong, Genetic Algorithms: Concepts and Designs, Springer, London, UK, 1999.

R. L. Haupt, “Thinned arrays using genetic algorithms,” IEEE Transactions on Antennas and Propagation, vol. 42, no. 7, pp. 993–999, 1994.

L. Cen, Z. L. Yu, W. Ser, and W. Cen, “Linear aperiodic array synthesis using an improved genetic algorithm,” IEEE Transactions on Antennas and Propagation, vol. 60, no. 2, pp. 895–902, 2012.

L. Cen, W. Ser, Z. L. Yu, S. Rahardja, and W. Cen, “Linear sparse array synthesis with minimum number of sensors,” IEEE Transactions on Antennas and Propagation, vol. 58, no. 3, pp. 720–726, 2010.

A. H. Hussein, H. H. Abdullah, A. M. Salem, S. Khamis, and M. Nasr, “Optimum design of linear antenna arrays using a hybrid MoM/GA algorithm,” IEEE Antennas and Wireless Propagation Letters, vol. 10, pp. 1232–1235, 2011.

R. F. Harrington, Field Computation by Moment Methods, IEEE Press, New York, NY, USA, 1993.

K. Deb, A. Pratap, S. Agarwal, and T. Meyarivan, “A fast and elitist multiobjective genetic algorithm: NSGA-II,” IEEE Transactions on Evolutionary Computation, vol. 6, no. 2, pp. 182–197, 2002.