# Shaped Beam Pattern Synthesis of Antenna Arrays Using Composite Differential Evolution with Eigenvector-Based Crossover Operator

Authors: Sotirios K. Goudos
Journal: International Journal of Antennas and Propagation
 (2015)

---

## Abstract

This paper addresses the problem of designing shaped beam patterns with arbitrary arrays subject to constraints. The constraints could include the sidelobe level suppression in specified angular intervals, the mainlobe halfpower beamwidth, and the predefined number of elements. In this paper, we propose a new Differential Evolution algorithm, which combines Composite DE with an eigenvector-based crossover operator (CODE-EIG). This operator utilizes eigenvectors of covariance matrix of individual solutions, which makes the crossover rotationally invariant. We apply this novel design method to shaped beam pattern synthesis for linear and conformal arrays. We compare this algorithm with other popular algorithms and DE variants. The results show CODE-EIG outperforms the other DE algorithms in terms of statistical results and convergence speed.


---

## Body

## 1. Introduction

Array synthesis is a classic and challenging optimization problem, which has been extensively studied using several analytical or stochastic methods [1–4]. Common optimization goals in array synthesis are the sidelobe level suppression and the matching of the main lobe to a desired pattern. Thus, the optimization problem is usually to find a set of element excitations and/or positions that closely match a desired pattern. The desired pattern shape can vary widely depending on the application. The shaped beam pattern synthesis problem has received wide attention over the years, and several methods or techniques have been reported in the literature [5–11]. Among others these methods include Tabu search [5], ant colony optimization [11], and linear programming [7]. Additionally, evolutionary algorithms (EAs) like Genetic Algorithms (GAs), Particle Swarm Optimization (PSO), and Differential Evolution (DE) [12], which is a population-based stochastic global optimization algorithm, have been applied to a variety of array design problems [13–21]. Several DE variants or strategies exist. One of the DE advantages is that very few control parameters have to be adjusted in each algorithm run.Composite DE (CODE) [22] is an adaptive DE variant, which combines three different trial vector generation strategies with three preset control parameter settings. The above combination is performed in a random way in order to generate trial vectors. The main advantage of CODE is that it has a simple structure and thus it is very easy to be implemented in any programming language. In this way significant computational costs spent on searching using a trial-and-error procedure can be avoided. The DE performance depends on control parameters, mutation strategies, and crossover operators. Most of the DE strategies or variants use the binomial crossover operator, which has been found to produce better results than the exponential crossover operator [23]. The authors in [24] propose an alternative crossover operator, namely, the eigenvector-based crossover. This operator utilizes the eigenvector information of the covariance matrix of the population to rotate the coordinate system. Additionally, this crossover operator can be applied to any DE variant. In this paper, we incorporate this scheme to the CODE algorithm and therefore introduce a new algorithm, the Composite DE with Eigenvector-Based Crossover (CODE-EIG).In this paper, we apply CODE-EIG to the shaped beam pattern synthesis problem. More specifically, we utilize the CODE-EIG algorithm for designing different linear arrays cases and a conformal array case. Additionally, in order to evaluate the new algorithm, we compare the algorithm’s performance with other DE strategies, PSO, and GAs.This paper is organized as follows: A brief description of the CODE and the CODE-EIG algorithm is given in Section2. Section 3 describes the problem formulation and the objective functions and presents the numerical results for different array design cases. Finally, the conclusion is given in Section 4.


## 2. The Composite Differential Evolution (CODE) Algorithm

A population in DE consists of NP vectorsx
-
G
,
i, i
=
1,2
,
…
,
N
P, where G is the generation number. The population is initialized randomly from a uniform distribution. Each D-dimensional vector represents a possible solution, which is expressed as(1)
x
-
G
,
i
=
x
G
,
1
i
,
x
G
,
2
i
,
…
,
x
G
,
j
i
,
…
,
x
G
,
D
i
.The population is initialized as follows:(2)
x
0
,
j
i
=
rand
j
0,1
⁡
x
j
,
U
-
x
j
,
L
+
x
j
,
L
j
=
1,2
,
…
,
D
,where x
j
,
L and x
j
,
U are D-dimensional vectors of the lower and upper bounds, respectively, and r
a
n
d
j
[
0,1
) is a uniformly distributed random number within [
0,1
).The initial population evolves in each generation with the use of three operators: mutation, crossover, and selection. Depending on the form of these operators, several DE variants or strategies exist in the literature [12, 25]. The choice of the best DE strategy depends on problem type [23]. In CODE the following three strategies are used for trial vector generation. These include DE/rand/1bin, DE/rand-to-best/2/bin, DE/rand/2/bin, and DE/current-to-rand/1 [26]. In these strategies, a mutant vector v
-
G
+
1
,
i for each target vector x
-
G
,
i is computed byDE/rand/1/bin(3a)
v
-
G
+
1
,
i
=
x
-
G
,
r
1
+
F
x
-
G
,
r
2
-
x
-
G
,
r
3
,
r
1
≠
r
2
≠
r
3
.
DE/rand/2/bin(3b)
v
-
G
+
1
,
i
=
x
-
G
,
r
1
+
F
x
-
G
,
r
2
-
x
-
G
,
r
3
+
F
x
-
G
,
r
4
-
x
-
G
,
r
5
,
r
1
≠
r
2
≠
r
3
≠
r
4
≠
r
5
.
DE/current-to-rand/1(3c)
v
-
G
+
1
,
i
=
x
-
G
,
i
+
rand
0,1
⁡
x
-
G
,
r
1
-
x
-
G
,
i
+
F
x
-
G
,
r
2
-
x
-
G
,
r
3
,
r
1
≠
r
2
≠
r
3
,where r
1, r
2, r
3, r
4, and r
5 are randomly chosen indices from the population, which are different from index i, F is a mutation control parameter, and r
a
n
d
[
0,1
) is a randomly generated number from a uniform distribution within the interval [
0,1
). The basic idea of CODE [22] is to generate a trial vector for each of the three above-mentioned strategies. The control parameters settings are randomly chosen. Therefore, CODE maintains a strategy candidate pool, consisting of the three strategies given in (3a), (3b), and (3c). A parameter candidate pool is also used with the settings given below:(4)
F
=
1.0
,
CR
=
0.1
,
F
=
1.0
,
CR
=
0.9
,
F
=
0.8
,
CR
=
0.2
.These are strategies and settings that have been widely used. Therefore, their properties have been extensively studied. The authors in [22] have modified the DE/rand/2/bin strategy by using a random number from the interval [
0,1
) instead of the first F mutation control parameter. The implementation of the CODE is quite simple. The best trial vector of the three generated is chosen. Then it is compared with the old vector and it replaces it if its fitness function value is smaller.After mutation in each generation for every vector of the population, the binary crossover operator is applied to generate a trial vectoru
-
G
+
1
,
i
=
(
u
G
+
1,1
i
,
u
G
+
1,2
i
,
…
,
u
G
+
1
,
j
i
,
…
,
u
G
+
1
,
D
i
) whose coordinates are given by(5)
u
G
+
1
,
j
i
=
v
G
+
1
,
j
i
,
if
rand
j
0,1
⁡
≤
CR  or
j
=
rn
i
,
x
G
+
1
,
j
i
,
i
f
r
a
n
d
j
0,1
>
CR
,
j
≠
rn
i
,where j
=
1,2
,
…
,
D, r
a
n
d
j
[
0,1
) is a number from a uniform random distribution within the interval [
0,1
), r
n
(
i
) a randomly chosen index from (
1,2
,
…
,
D
), and CR the crossover constant from the interval [
0,1
].In CODE, the crossover operator is not applied for the DE/current-to-rand/1 strategy. CODE uses a greedy selection operator, which for minimization problems is defined by(6)
x
-
G
+
1
,
i
=
u
-
G
+
1
,
i
,
if
f
u
-
G
+
1
,
i
<
f
x
-
G
,
i
,
x
-
G
,
i
,
otherwise
,where f
(
u
-
G
+
1
,
i
), f
(
x
-
G
,
i
) are the fitness values of the trial and the old vector, respectively. Therefore, the newly found trial vector u
-
G
+
1
,
i replaces the old vector x
-
G
,
i only when it produces a lower objective function value than the old one. Otherwise, the old vector remains in the next generation. The stopping criterion for the CODE or DE in general is usually the generation number or the number of objective function evaluations.


### 2.1. CODE with Eigenvector-Based Crossover Operator (CODE-EIG)

The authors in [24] propose an eigenvector-based crossover operator, which utilizes eigenvectors of covariance matrix of individual solutions. Thus, the crossover is rotationally invariant. To avoid losing diversity of population, the offspring can be stochastically born from the parents with either the standard coordinate system or the rotated coordinate system. The authors in [24] also introduce a new parameter to control the probability of selecting one of the coordinate systems. They have shown that this scheme can increase the population diversity and prevent premature convergence. Additionally, another significant advantage of this operator is that it can be applied to any crossover strategy with minimal changes. Therefore, it can enhance any existing DE variant. The main idea is to exchange the information between the target vector and the mutant vector in the eigenvector basis instead of the natural basis. We consider the same population of NP D-dimensional vectors defined in the previous section. The covariance between ith and jth dimension of the population in the Gth generation is given by [24](7)
cov
⁡
i
,
j
=
∑
k
=
1
N
P
x
G
,
i
k
-
μ
G
,
i
x
G
,
j
k
-
μ
G
,
j
N
P
-
1
,where μ
G
,
i, μ
G
,
j are the mean values of the variables in the ith and jth dimension, respectively.To compute the eigenvector basis we need to factorize the covariance matrixC
G
=
(
c
i
j
,
c
i
j
=
c
o
v
(
i
,
j
)
) into a canonical form(8)
C
G
=
Q
G
Λ
G
Q
G
-
1
,where Q
G is the square matrix (D
×
D) whose ith column is the eigenvector q
-
G
,
i of C
G and Λ
G is the diagonal matrix whose diagonal elements are the corresponding eigenvalues. The factorization of a matrix into a canonical form is called eigendecomposition. The authors in [24] use Jacobi’s method [27] for eigendecomposition. When the eigenvector basis is found, the ith target vectors x
-
G
,
i can be expressed by Q
G
·
x
-
G
,
i; the ith mutant vectors v
-
G
+
1
,
i can be expressed by Q
G
·
v
-
G
+
1
,
i. Then, a predefined crossover operator, such as binomial crossover, will exchange some of the elements of the mutant vector with some of the elements of its target vector to form a trial vector. The trial vector is then given by [24](9)
u
-
G
+
1
,
i
=
Q
G
∗
·
xover
Q
G
·
x
-
G
,
i
,
Q
G
·
v
-
G
+
1
,
i
,
if
rand
0,1
⁡
≤
P
,
xover
x
-
G
,
i
,
v
-
G
+
1
,
i
,
otherwise
,where Q
G
∗ is the conjugate transpose of the eigenvector basis Q
G and x
o
v
e
r
x
-
,
y
- is a crossover operator on two vectors x
- and y
-, where P is a new control parameter introduced in [24] which is called eigenvector ratio between 0 and 1 that determines the ratio of the eigenvector-based crossover operator and the other crossover operator. The main feature of this approach is that no matter how the crossover operator exchanges the elements in the eigenvector basis, the crossover behavior will become rotationally invariant in the natural basis.Therefore, to create the CODE-EIG variant we use (3a), (3b), and (3c) to generate the mutant vector and (9) to create the trial vector. Additionally, this new variant requires the setting of the eigenvector ratio P. It must be noted that when P
=
1 then only the eigenvector crossover operator is used, while when P
=
0 then only the binary crossover operator is used.The CODE-EIG algorithm is outlined below:(1)
Initialize the CODE-EIG parameters. That is, set the maximum number of generationsG
m
a
x, the population size number NP, and the eigenvector ratio P.
(2)
Initialize a uniformly distributed random population of NP individuals. SetG
=
1.
(3)
Evaluate objective function value for every vector of the population.
(4)
For each vector of the population, generate three mutant vectorsv
-
G
+
1
,
i
_
1, v
-
G
+
1
,
i
_
2, v
-
G
+
1
,
i
_
3, each with one of the three strategies given in (3a), (3b), and (3c) and with a control parameter setting randomly selected from the parameter candidate pool.
(5)
Apply the eigenvector-based crossover operator according to (9) and create three trial vectors u
-
G
+
1
,
i
_
1, u
-
G
+
1
,
i
_
2, u
-
G
+
1
,
i
_
3.
(6)
Evaluate objective function values for the trial vectors.
(7)
Select the best of the trial vectors.
(8)
Apply the selection operator according to (6).
(9)
Keep the best solution so far.
(10)
If the maximum number of generationsG
m
a
x is reached, then exit; otherwise increment generation number and go to step 4 for the next generation.The CODE-EIG flowchart is depicted in Figure 1.Figure 1
CODE-EIG flowchart.


## 2.1. CODE with Eigenvector-Based Crossover Operator (CODE-EIG)

The authors in [24] propose an eigenvector-based crossover operator, which utilizes eigenvectors of covariance matrix of individual solutions. Thus, the crossover is rotationally invariant. To avoid losing diversity of population, the offspring can be stochastically born from the parents with either the standard coordinate system or the rotated coordinate system. The authors in [24] also introduce a new parameter to control the probability of selecting one of the coordinate systems. They have shown that this scheme can increase the population diversity and prevent premature convergence. Additionally, another significant advantage of this operator is that it can be applied to any crossover strategy with minimal changes. Therefore, it can enhance any existing DE variant. The main idea is to exchange the information between the target vector and the mutant vector in the eigenvector basis instead of the natural basis. We consider the same population of NP D-dimensional vectors defined in the previous section. The covariance between ith and jth dimension of the population in the Gth generation is given by [24](7)
cov
⁡
i
,
j
=
∑
k
=
1
N
P
x
G
,
i
k
-
μ
G
,
i
x
G
,
j
k
-
μ
G
,
j
N
P
-
1
,where μ
G
,
i, μ
G
,
j are the mean values of the variables in the ith and jth dimension, respectively.To compute the eigenvector basis we need to factorize the covariance matrixC
G
=
(
c
i
j
,
c
i
j
=
c
o
v
(
i
,
j
)
) into a canonical form(8)
C
G
=
Q
G
Λ
G
Q
G
-
1
,where Q
G is the square matrix (D
×
D) whose ith column is the eigenvector q
-
G
,
i of C
G and Λ
G is the diagonal matrix whose diagonal elements are the corresponding eigenvalues. The factorization of a matrix into a canonical form is called eigendecomposition. The authors in [24] use Jacobi’s method [27] for eigendecomposition. When the eigenvector basis is found, the ith target vectors x
-
G
,
i can be expressed by Q
G
·
x
-
G
,
i; the ith mutant vectors v
-
G
+
1
,
i can be expressed by Q
G
·
v
-
G
+
1
,
i. Then, a predefined crossover operator, such as binomial crossover, will exchange some of the elements of the mutant vector with some of the elements of its target vector to form a trial vector. The trial vector is then given by [24](9)
u
-
G
+
1
,
i
=
Q
G
∗
·
xover
Q
G
·
x
-
G
,
i
,
Q
G
·
v
-
G
+
1
,
i
,
if
rand
0,1
⁡
≤
P
,
xover
x
-
G
,
i
,
v
-
G
+
1
,
i
,
otherwise
,where Q
G
∗ is the conjugate transpose of the eigenvector basis Q
G and x
o
v
e
r
x
-
,
y
- is a crossover operator on two vectors x
- and y
-, where P is a new control parameter introduced in [24] which is called eigenvector ratio between 0 and 1 that determines the ratio of the eigenvector-based crossover operator and the other crossover operator. The main feature of this approach is that no matter how the crossover operator exchanges the elements in the eigenvector basis, the crossover behavior will become rotationally invariant in the natural basis.Therefore, to create the CODE-EIG variant we use (3a), (3b), and (3c) to generate the mutant vector and (9) to create the trial vector. Additionally, this new variant requires the setting of the eigenvector ratio P. It must be noted that when P
=
1 then only the eigenvector crossover operator is used, while when P
=
0 then only the binary crossover operator is used.The CODE-EIG algorithm is outlined below:(1)
Initialize the CODE-EIG parameters. That is, set the maximum number of generationsG
m
a
x, the population size number NP, and the eigenvector ratio P.
(2)
Initialize a uniformly distributed random population of NP individuals. SetG
=
1.
(3)
Evaluate objective function value for every vector of the population.
(4)
For each vector of the population, generate three mutant vectorsv
-
G
+
1
,
i
_
1, v
-
G
+
1
,
i
_
2, v
-
G
+
1
,
i
_
3, each with one of the three strategies given in (3a), (3b), and (3c) and with a control parameter setting randomly selected from the parameter candidate pool.
(5)
Apply the eigenvector-based crossover operator according to (9) and create three trial vectors u
-
G
+
1
,
i
_
1, u
-
G
+
1
,
i
_
2, u
-
G
+
1
,
i
_
3.
(6)
Evaluate objective function values for the trial vectors.
(7)
Select the best of the trial vectors.
(8)
Apply the selection operator according to (6).
(9)
Keep the best solution so far.
(10)
If the maximum number of generationsG
m
a
x is reached, then exit; otherwise increment generation number and go to step 4 for the next generation.The CODE-EIG flowchart is depicted in Figure 1.Figure 1
CODE-EIG flowchart.


## 3. Numerical Results

In this paper, we compare the proposed CODE-EIG with three other DE algorithms, PSO, and GAs on two different array design cases. These algorithms are the CODE algorithm [26], the popular DE/rand/1/bin strategy (binomial crossover), the DE/rand/1/eig strategy (eigenvector-based mixed with binomial crossover), the Constriction Factor PSO [28], and a real coded GA. In the DE/rand/1/eig and DE/rand/1/bin strategies, F is set to 0.5 and CR to 0.9. In the PSO, the cognitive learning factor c
1 and the social learning factor c
2 are both set to 2.05. For the GA the crossover probability is set to 0.9 and the mutation probability is set to 0.1. The eigenvector ratio is set to 0.5 for both CODE-EIG and DE/rand/1/eig. For the cases where the algorithms are compared, the population size and the number of generations are set equal for all algorithms. The stopping criterion for all algorithms is the generation number. The same initial conditions are used for all algorithms. All algorithms are compiled using the same compiler (Microsoft Visual C++ 2010) in a PC with Intel Core i5-2410M, at 2.30 GHz with 4 GB RAM running Windows 7.


### 3.1. Linear Array Synthesis

We consider anN-element linear array of equally spaced isotropic elements. The array factor is expressed as(10)
AF
u
,
I
-
,
φ
-
=
∑
n
=
1
N
I
n
e
j
2
π
/
λ
n
d
u
+
φ
n
,where λ is the wavelength, I
n and φ
n are the amplitude and the phase of the nth element, respectively, I
- and φ
- are the corresponding vectors, d is the distance between two adjacent elements, u
=
s
i
n
(
ϑ
) is the direction cosine, and ϑ is the steering angle measured from broadside of the array.For a symmetrically excited array (10) becomes(11)
AF
u
,
I
-
,
φ
-
=
2
∑
n
=
1
M
I
n
cos
⁡
2
π
λ
d
n
u
+
φ
n
+
I
0
,where M
=
N
/
2 is the largest integer less than or equal to N
/
2. For an even number of elements I
0
=
0 and for an odd number of elements we set I
0
=
1. We assume that d
=
0.5
λ for all cases.The shaped beam design cases are chosen such as the desired pattern is a cosecant beam and outside the main beam the peak sidelobe level (PSLL) should not exceed predefined values [5, 7, 11]. We consider case 1, where the desired pattern is given by the following expression:(12)
AF
d
e
s
i
r
e
d
u
=
cosecant
u
,
for
0.1
≤
u
≤
0.5
,
0
,
elsewhere,
PSLL
d
u
≤
-
30
dB
,where P
S
L
L
d is the desired peak sidelobe level.One may observe that theP
S
L
L
d is set to −30 dB instead of −25 dB in [5, 11]. The second case uses the same cosecant main beam in the sector of u
∈
[
0.1,0.5
] as in (12), while the P
S
L
L
d is given by [5, 7, 11](13)
PSLL
d
u
=
-
35
dB
,
for
-
1
≤
u
≤
-
0.35
,
-
25
dB
,
for
-
0.35
≤
u
≤
0
,
-
40
dB
,
for
0.6
≤
u
≤
1
.The requirement of matching a desired shaped beam pattern can be expressed as follows:(14)
F
x
-
=
∑
k
=
1
K
w
u
k
AF
d
e
s
i
r
e
d
u
k
-
AF
c
a
l
c
u
k
,
x
-
2
,
w
u
=
500
,
for
0.1
<
u
≤
0.5
,
1000
,
for
u
=
0.1
,
300
,
if
PSLL
c
u
>
PSLL
d
u
,
1
,
otherwise
,where P
S
L
L
c is the calculated peak sidelobe level and x
- is the vector of the unknown variables (the array amplitudes and phases). We have found the above weight values based on a trial-and-error procedure. The population size is set to 100 and the maximum number of generations is set to 1000 for all linear array design cases.The first example case is that of a symmetrically excited array of 20 and 18 elements. Thus, the total number of unknowns is 20 and 18, respectively. The desired pattern is that of case 1. In this case, we compare the algorithms performance forN
=
20 elements. The best value, the worst value, the mean, and the standard deviation of the last generation computed by each algorithm are presented here. Figure 2 shows the best radiation patterns obtained by CODE-EIG for this case. The obtained patterns are in good agreement with the desired one. The 20-element array obtained a PSLL value of −30 dB, while the 18-element array obtained a PLL value of −29.96 dB. Table 1 holds the optimum element amplitudes and phases derived for the first two arrays. The convergence rate plot for this case for N
=
20 is shown in Figure 3. The CODE-EIG strategy converges faster than the other DE algorithms. The DE/rand/1/eig strategy converges faster than the CODE and the DE/rand/1/bin strategy. PSO and GA converge slower than the DE algorithms. Table 2 presents the statistical results for all algorithms. Both CODE-EIG and DE/rand/1/eig strategy produce about the same best result. However, CODE-EIG is better than DE/rand/1/eig in terms of mean and standard deviation values. We also notice that the DE/rand/1/eig results are better than the CODE and DE/rand/1/bin strategy. The PSO and GA obtained results present higher mean and standard deviation values. The DE/rand/1/bin strategy is the fastest among all the algorithms in terms of computational time. The PSO and GA are also faster than the other DE algorithms.Table 1
Case  1. Element amplitudes and phasesN
=
20 and N
=
18 symmetrically excited isotropic elements.
n
I
n
φ
n (deg)
I
n
φ
n (deg)
1
0.909
163.906
0.997
165.307
2
0.762
133.396
0.851
137.129
3
0.535
109.236
0.618
113.545
4
0.348
99.892
0.393
101.486
5
0.288
100.513
0.290
105.145
6
0.279
88.565
0.288
97.232
7
0.226
67.017
0.249
71.910
8
0.144
46.497
0.151
51.300
9
0.073
41.120
0.083
75.256
10
0.132
63.919Table 2
Comparative results for symmetrically excited design Case  1 forN
=
20. The smaller values are in bold.
Algorithm
Best
Worst
Mean
St. dev.
CPU time (sec)
CODE-EIG
34619.53
37476.33
36030.69
1121.19
239.63
CODE
37416.86
39859.79
38629.17
1190.15
243.36
DE/rand/1/eig
34619.57
37497.31
36866.06
1405.14
217.44
DE/rand/1/bin
37406.04
39856.17
38370.45
1168.07
193.36
PSO
37534.00
54611.18
43211.50
5502.29
206.31
GA
36993.32
48823.57
43758.32
3804.47
202.42Figure 2
Case  1. Radiation patterns obtained forN
=
20 and N
=
18 symmetrically excited isotropic elements.Figure 3
Case  1. Convergence rate plot for the 20-element array design case.The second example case is again a symmetrically excited array of 20 and 18 elements. However, the desired pattern is that of case 2. Again, we compare the algorithms performance forN
=
20 elements. Figure 4 presents the patterns of the best designs obtained for both N
=
20 and N
=
18 elements. We notice that the derived patterns closely match the desired one, while the 20-element array is a better approximation than the 18-element one. The maximum SLL of the 20-element case is −25.09 dB while for the 18-element case the SLL value is −24.94 dB. Table 3 presents the optimum amplitudes and phases obtained by the CODE-EIG strategy for both arrays. Figure 5 shows the convergence rate plots for the 20-element case. It is obvious that CODE-EIG converges faster than the other algorithms. CODE is the second faster algorithm for this case. All DE algorithms converge at similar final values. PSO and GA converge slower than the other algorithms and a higher final value. The mean values and the standard deviation values are larger for PSO and GA. Table 4 presents the statistical results of the objective function values for both arrays. CODE-EIG clearly outperforms the other algorithms. It has obtained the lower best, mean, and standard deviation values. The DE/rand/1/eig strategy has produced the best worst value and produces better results than the CODE and the DE/rand/1/bin. Regarding computational time, again the DE/rand/1/bin is faster than the other algorithms. The CODE-EIG is faster than CODE and DE/rand/1/eig in this case.Table 3
Case  2. Element amplitudes and phasesN
=
20 and N
=
18 symmetrically excited isotropic elements.
n
I
n
φ
n (deg)
I
n
φ
n (deg)
1
1.000
162.224
0.999
162.997
2
0.811
129.582
0.765
131.770
3
0.535
107.654
0.473
114.275
4
0.360
105.240
0.358
114.863
5
0.340
99.504
0.351
105.281
6
0.314
78.941
0.286
86.582
7
0.239
59.991
0.182
82.091
8
0.169
59.714
0.137
83.111
9
0.142
59.948
0.109
83.064
10
0.111
55.830Table 4
Comparative results for symmetrical excited design Case  2 forN
=
20. The smaller values are in bold.
Algorithm
Best
Worst
Mean
St. dev.
Mean CPU time (sec)
CODE-EIG
9385.09
14856.80
11553.94
2638.87
212.69
CODE
9502.48
22680.59
15026.70
3014.12
300.09
DE/rand/1/eig
9493.09
14718.99
12051.00
2665.16
226.53
DE/rand/1/bin
14716.18
22681.45
18394.28
3762.84
184.83
PSO
10348.90
28280.86
21872.56
5090.96
188.25
GA
15993.97
28415.97
19632.79
4315.33
190.61Figure 4
Case  2. Radiation patterns obtained forN
=
20 and N
=
18 symmetrically excited isotropic elements.Figure 5
Case  2. Convergence rate plot for the 20-element array design case.Next, the examples that follow design the same two cases but with asymmetrically excited arrays. Therefore, the total number of unknowns is2
N for these cases. The third example is that of an asymmetrically excited array with N
=
16, 17, and 18 elements. The desired pattern is that of case 1. The best obtained patterns by CODE-EIG are depicted in Figure 6. We notice that the three patterns approximate the desired pattern. The peak SLL values are −30.20 dB, −30 dB, and −29.99 dB for N
=
18, 17, and 16, respectively. The element amplitude and phases for the three patterns are presented in Table 5.Table 5
Case  1. Element amplitudes and phases forN
=
18, 17,  and 15 asymmetrically excited isotropic elements.
n
I
n
φ
n (deg)
I
n
φ
n (deg)
I
n
φ
n (deg)
1
0.136
179.895
0.057
0.774
0.139
179.976
2
0.268
180.000
0
177.597
0.120
179.946
3
0.342
149.887
0.108
179.997
0.167
179.975
4
0.370
137.275
0.104
179.914
0.336
179.994
5
0.455
136.661
0.153
179.702
0.572
161.028
6
0.625
134.342
0.277
179.926
0.853
139.058
7
0.868
114.851
0.549
174.629
0.999
107.863
8
0.990
91.554
0.839
144.046
0.972
81.897
9
0.973
66.089
0.985
116.748
0.871
60.612
10
0.833
41.468
1.000
91.827
0.630
42.012
11
0.560
10.548
0.905
67.285
0.430
37.924
12
0.291
0
0.674
47.620
0.411
39.555
13
0.142
2.974
0.472
39.502
0.386
27.793
14
0.108
0.884
0.355
36.326
0.311
0.003
15
0.065
0
0.354
28.654
0.161
0.004
16
0
0.630
0.241
0.021
0.036
2.615
17
0.048
177.224
0.139
0.099
18
0.070
135.981Figure 6
Case  1. Radiation patterns obtained forN
=
18, 17, and 16 asymmetrically excited isotropic elements.The fourth example presents the design case 2 as the desired pattern obtained by an asymmetrically excited linear array of 20 and 18 isotropic elements. The author has not been able to find in the literature an example that uses this design case asymmetrically excited. The best obtained patterns by CODE-EIG are illustrated in Figure7. We notice that the 20-element array matches better the desired pattern than the 18-element one. The peak SLL values are −25.01 dB and −24.99 dB for N
=
20 and N
=
18 elements, respectively. Table 6 holds the corresponding element amplitudes and phases for both arrays.Table 6
Case  2. Element amplitudes and phases forN
=
20  and 18 asymmetrically excited isotropic elements.
n
I
n
φ
n (deg)
I
n
φ
n (deg)
1
0.155
179.993
0.132
179.993
2
0.297
167.574
0.182
179.93
3
0.356
133.836
0.246
177.114
4
0.227
104.702
0.392
174.602
5
0.081
179.853
0.552
155.925
6
0.460
179.979
0.643
134.614
7
0.809
146.989
0.670
117.929
8
1.000
114.094
0.766
107.713
9
0.980
86.945
0.926
90.237
10
0.844
68.696
0.980
63.686
11
0.757
56.667
0.807
31.301
12
0.698
40.979
0.463
0
13
0.595
20.449
0.120
0.454
14
0.413
5.048
0.226
69.332
15
0.282
0.173
0.341
44.786
16
0.244
3.129
0.305
20.131
17
0.196
0
0.194
0.009
18
0.109
0.070
0.108
0.022
19
0.070
24.731
20
0.094
25.046Figure 7
Case  6. Radiation patterns obtained forN
=
20 and 18 asymmetrically excited isotropic elements.


### 3.2. Conformal Array Synthesis

The final design case presents a conformal cylindrical array design. This case has been solved in [4] using simulated annealing. In [7] another design is obtained using linear programming. We assume a cylindrical array consisting of MXN elements in the ϕ and z directions, respectively. For this case M
=
N
=
8. The radius R
0 of the cylindrical surface is 5
λ, while the cylinder height is h
0
=
3.5
λ. The angular region of the array extends from -
ϕ
0 to ϕ
0, where ϕ
0
=
20
°. The array elements are uniformly spaced by 0.5
λ.The far-field radiation pattern of the conformal cylindrical array is given by(15)
FF
ϑ
,
ϕ
,
I
-
,
φ
-
=
∑
m
=
1
M
∑
n
=
1
N
I
m
n
×
e
l
ϑ
,
ϕ
×
e
j
φ
m
n
+
2
π
/
λ
x
m
n
sin
⁡
ϑ
cos
⁡
ϕ
+
y
m
n
sin
⁡
ϑ
sin
⁡
ϕ
+
z
m
n
cos
⁡
ϑ
,where m, n are the element indices in the ϕ and z directions, respectively, and e
l
(
ϑ
,
ϕ
) is the element pattern given by [7](16)
e
l
ϑ
,
ϕ
=
sin
⁡
ϑ
cos
⁡
ϕ
,
for
ϕ
-
ϕ
m
≤
π
2
,
0
,
otherwise
,where ϕ
m is the element azimuthal position.The positions of the array elements are(17)
x
m
n
=
R
0
cos
⁡
ϕ
m
,
y
m
n
=
R
0
sin
⁡
ϕ
m
,
z
m
n
=
-
h
0
2
+
n
×
0.5
λ
.The desired far-field pattern, F
F
d, has a cosecant square dependency in the ϕ
=
0
° plane for an elevation angle ranging between ϑ
s
t
a
r
t
c
o
s and ϑ
e
n
d
c
o
s. The setting of the half-power beamwidth (HPBW) in the ϑ
=
90
° plane to a desired value, H
P
B
W
d
ϑ
=
90
°, within a tolerance percentage is also a requirement. Additionally, the peak sidelobe levels outside the shaped beam region must remain below a predefined desired value, P
S
L
L
d.For this case, we express the objective function as(18)
F
I
-
,
φ
-
=
1
K
1
∑
k
=
1
K
1
F
F
d
ϑ
k
,
0
,
I
-
,
φ
-
-
F
F
c
ϑ
k
,
0
,
I
-
,
φ
-
2
+
Ξ
·
max
⁡
0
,
PSLL
d
-
PSLL
c
+
Ξ
·
max
⁡
0
,
HPBW
c
ϑ
=
90
°
-
HPBW
d
ϑ
=
90
°
HPBW
d
ϑ
=
90
°
-
T
H
P
B
W
,where F
F
c is the calculated far-field pattern in the ϕ
=
0
° plane, K
1 is the number of pattern samples taken for the ϕ
=
0
° plane, P
S
L
L
c is the calculated peak sidelobe level in a set of ϕ
-
c
u
t
s, H
P
B
W
c
ϑ
=
90
° is the calculated half-power beamwidth in the ϑ
=
90
° plane, and T
H
P
B
W is the desired beamwidth tolerance.The design parameters we have selected for this case areϑ
s
t
a
r
t
c
o
s
=
100
°, ϑ
e
n
d
c
o
s
=
140
°, H
P
B
W
d
ϑ
=
90
°
=
16
°, T
H
P
B
W
=
0.01, and P
S
L
L
d
=
-
27 dB (instead of −25 dB in [4, 7]). We set the population size to 100 vectors and the maximum number of iterations is set to 3000. Figures 8(a), 8(b), and 8(c) depict the best obtained radiation pattern with CODE-EIG for the ϕ
=
0
°, ϑ
=
90
°, and a 3D pattern, respectively. One may notice from the ϕ
=
0
° plane pattern that the pattern found has PSLL value of −27.01 dB, which is about 3 dB below the one from [4], while it is about 0.4 dB higher than the one from [7]. Moreover, the ϑ
=
90
° plane pattern has a HPBW about 16
° and the PSLL value in this pattern is below the corresponding values from [4, 7]. Table 7 holds the amplitudes and phases of the best obtained pattern found by CODE-EIG.Table 7
Excitation coefficients for the best conformal design obtained by CODE-EIG. The phase values are in degrees.
m
1,8
2,7
3,6
4,5
n
I
m
n
φ
m
n
I
m
n
φ
m
n
I
m
n
φ
m
n
I
m
n
φ
m
n
1
0.27
179.14
0.15
168.91
0.33
96.05
0.28
104.38
2
0.43
−61.11
0.39
163.30
0.61
147.49
0.74
130.96
3
0.57
−116.48
0.51
−161.67
1.06
−157.97
1.08
155.87
4
0.81
−20.42
1.04
−106.93
1.31
−133.54
1.37
−137.05
5
0.24
49.47
0.69
−53.72
0.87
−50.12
1.05
−108.07
6
0.31
−34.05
0.23
−6.12
0.63
−123.44
0.58
−83.90
7
0.26
50.05
0.53
−76.22
0.57
−33.05
0.64
−92.52
8
0.08
37.82
0.25
30.14
0.36
−23.14
0.47
−46.99Figure 8
Conformal array radiation far-field patterns in (a)ϕ
=
0
∘ plane, (b) ϑ
=
9
0
∘ plane, and (c) 3D pattern.
(a)
(b)
(c)


### 3.3. CODE-EIG Tuning Parameter Sensitivity

The main control parameter of CODE-EIG is the eigenvector ratioP. In order to explore the CODE-EIG sensitivity to control parameter selection we evaluated the first example using different eigenvector ratio settings. The number of elements of the pattern is N
=
20. We set the number of generations to 1000 and the population size to 100. We run the algorithm for 50 independent trials for each case. The best value, the worst value, the mean, and the standard deviation at the last generation are presented here. Table 8 shows the results. CODE-EIG performance is relatively robust to different P values, as long as P is 0.5 or 0.8. Higher or lower eigenvector ratio values result in slightly worse statistical values. This conclusion is in agreement with [24], where the authors show the eigenvector ratio value of 0.5 produces better results.Table 8
Comparative results for symmetrically excited design Case  1 forN
=
20 using different eigenvector ratio settings. The smaller values are in bold.
P
Best
Worst
Mean
St. dev.
0.5
34619.53
37476.33
36030.69
1121.19
0.2
34647.83
37471.51
36779.67
1368.40
0.4
34655.49
37502.64
36573.75
1308.63
0.8
34641.35
37528.71
36510.02
1294.57
0.9
34652.01
37527.25
36806.05
1387.19


## 3.1. Linear Array Synthesis

We consider anN-element linear array of equally spaced isotropic elements. The array factor is expressed as(10)
AF
u
,
I
-
,
φ
-
=
∑
n
=
1
N
I
n
e
j
2
π
/
λ
n
d
u
+
φ
n
,where λ is the wavelength, I
n and φ
n are the amplitude and the phase of the nth element, respectively, I
- and φ
- are the corresponding vectors, d is the distance between two adjacent elements, u
=
s
i
n
(
ϑ
) is the direction cosine, and ϑ is the steering angle measured from broadside of the array.For a symmetrically excited array (10) becomes(11)
AF
u
,
I
-
,
φ
-
=
2
∑
n
=
1
M
I
n
cos
⁡
2
π
λ
d
n
u
+
φ
n
+
I
0
,where M
=
N
/
2 is the largest integer less than or equal to N
/
2. For an even number of elements I
0
=
0 and for an odd number of elements we set I
0
=
1. We assume that d
=
0.5
λ for all cases.The shaped beam design cases are chosen such as the desired pattern is a cosecant beam and outside the main beam the peak sidelobe level (PSLL) should not exceed predefined values [5, 7, 11]. We consider case 1, where the desired pattern is given by the following expression:(12)
AF
d
e
s
i
r
e
d
u
=
cosecant
u
,
for
0.1
≤
u
≤
0.5
,
0
,
elsewhere,
PSLL
d
u
≤
-
30
dB
,where P
S
L
L
d is the desired peak sidelobe level.One may observe that theP
S
L
L
d is set to −30 dB instead of −25 dB in [5, 11]. The second case uses the same cosecant main beam in the sector of u
∈
[
0.1,0.5
] as in (12), while the P
S
L
L
d is given by [5, 7, 11](13)
PSLL
d
u
=
-
35
dB
,
for
-
1
≤
u
≤
-
0.35
,
-
25
dB
,
for
-
0.35
≤
u
≤
0
,
-
40
dB
,
for
0.6
≤
u
≤
1
.The requirement of matching a desired shaped beam pattern can be expressed as follows:(14)
F
x
-
=
∑
k
=
1
K
w
u
k
AF
d
e
s
i
r
e
d
u
k
-
AF
c
a
l
c
u
k
,
x
-
2
,
w
u
=
500
,
for
0.1
<
u
≤
0.5
,
1000
,
for
u
=
0.1
,
300
,
if
PSLL
c
u
>
PSLL
d
u
,
1
,
otherwise
,where P
S
L
L
c is the calculated peak sidelobe level and x
- is the vector of the unknown variables (the array amplitudes and phases). We have found the above weight values based on a trial-and-error procedure. The population size is set to 100 and the maximum number of generations is set to 1000 for all linear array design cases.The first example case is that of a symmetrically excited array of 20 and 18 elements. Thus, the total number of unknowns is 20 and 18, respectively. The desired pattern is that of case 1. In this case, we compare the algorithms performance forN
=
20 elements. The best value, the worst value, the mean, and the standard deviation of the last generation computed by each algorithm are presented here. Figure 2 shows the best radiation patterns obtained by CODE-EIG for this case. The obtained patterns are in good agreement with the desired one. The 20-element array obtained a PSLL value of −30 dB, while the 18-element array obtained a PLL value of −29.96 dB. Table 1 holds the optimum element amplitudes and phases derived for the first two arrays. The convergence rate plot for this case for N
=
20 is shown in Figure 3. The CODE-EIG strategy converges faster than the other DE algorithms. The DE/rand/1/eig strategy converges faster than the CODE and the DE/rand/1/bin strategy. PSO and GA converge slower than the DE algorithms. Table 2 presents the statistical results for all algorithms. Both CODE-EIG and DE/rand/1/eig strategy produce about the same best result. However, CODE-EIG is better than DE/rand/1/eig in terms of mean and standard deviation values. We also notice that the DE/rand/1/eig results are better than the CODE and DE/rand/1/bin strategy. The PSO and GA obtained results present higher mean and standard deviation values. The DE/rand/1/bin strategy is the fastest among all the algorithms in terms of computational time. The PSO and GA are also faster than the other DE algorithms.Table 1
Case  1. Element amplitudes and phasesN
=
20 and N
=
18 symmetrically excited isotropic elements.
n
I
n
φ
n (deg)
I
n
φ
n (deg)
1
0.909
163.906
0.997
165.307
2
0.762
133.396
0.851
137.129
3
0.535
109.236
0.618
113.545
4
0.348
99.892
0.393
101.486
5
0.288
100.513
0.290
105.145
6
0.279
88.565
0.288
97.232
7
0.226
67.017
0.249
71.910
8
0.144
46.497
0.151
51.300
9
0.073
41.120
0.083
75.256
10
0.132
63.919Table 2
Comparative results for symmetrically excited design Case  1 forN
=
20. The smaller values are in bold.
Algorithm
Best
Worst
Mean
St. dev.
CPU time (sec)
CODE-EIG
34619.53
37476.33
36030.69
1121.19
239.63
CODE
37416.86
39859.79
38629.17
1190.15
243.36
DE/rand/1/eig
34619.57
37497.31
36866.06
1405.14
217.44
DE/rand/1/bin
37406.04
39856.17
38370.45
1168.07
193.36
PSO
37534.00
54611.18
43211.50
5502.29
206.31
GA
36993.32
48823.57
43758.32
3804.47
202.42Figure 2
Case  1. Radiation patterns obtained forN
=
20 and N
=
18 symmetrically excited isotropic elements.Figure 3
Case  1. Convergence rate plot for the 20-element array design case.The second example case is again a symmetrically excited array of 20 and 18 elements. However, the desired pattern is that of case 2. Again, we compare the algorithms performance forN
=
20 elements. Figure 4 presents the patterns of the best designs obtained for both N
=
20 and N
=
18 elements. We notice that the derived patterns closely match the desired one, while the 20-element array is a better approximation than the 18-element one. The maximum SLL of the 20-element case is −25.09 dB while for the 18-element case the SLL value is −24.94 dB. Table 3 presents the optimum amplitudes and phases obtained by the CODE-EIG strategy for both arrays. Figure 5 shows the convergence rate plots for the 20-element case. It is obvious that CODE-EIG converges faster than the other algorithms. CODE is the second faster algorithm for this case. All DE algorithms converge at similar final values. PSO and GA converge slower than the other algorithms and a higher final value. The mean values and the standard deviation values are larger for PSO and GA. Table 4 presents the statistical results of the objective function values for both arrays. CODE-EIG clearly outperforms the other algorithms. It has obtained the lower best, mean, and standard deviation values. The DE/rand/1/eig strategy has produced the best worst value and produces better results than the CODE and the DE/rand/1/bin. Regarding computational time, again the DE/rand/1/bin is faster than the other algorithms. The CODE-EIG is faster than CODE and DE/rand/1/eig in this case.Table 3
Case  2. Element amplitudes and phasesN
=
20 and N
=
18 symmetrically excited isotropic elements.
n
I
n
φ
n (deg)
I
n
φ
n (deg)
1
1.000
162.224
0.999
162.997
2
0.811
129.582
0.765
131.770
3
0.535
107.654
0.473
114.275
4
0.360
105.240
0.358
114.863
5
0.340
99.504
0.351
105.281
6
0.314
78.941
0.286
86.582
7
0.239
59.991
0.182
82.091
8
0.169
59.714
0.137
83.111
9
0.142
59.948
0.109
83.064
10
0.111
55.830Table 4
Comparative results for symmetrical excited design Case  2 forN
=
20. The smaller values are in bold.
Algorithm
Best
Worst
Mean
St. dev.
Mean CPU time (sec)
CODE-EIG
9385.09
14856.80
11553.94
2638.87
212.69
CODE
9502.48
22680.59
15026.70
3014.12
300.09
DE/rand/1/eig
9493.09
14718.99
12051.00
2665.16
226.53
DE/rand/1/bin
14716.18
22681.45
18394.28
3762.84
184.83
PSO
10348.90
28280.86
21872.56
5090.96
188.25
GA
15993.97
28415.97
19632.79
4315.33
190.61Figure 4
Case  2. Radiation patterns obtained forN
=
20 and N
=
18 symmetrically excited isotropic elements.Figure 5
Case  2. Convergence rate plot for the 20-element array design case.Next, the examples that follow design the same two cases but with asymmetrically excited arrays. Therefore, the total number of unknowns is2
N for these cases. The third example is that of an asymmetrically excited array with N
=
16, 17, and 18 elements. The desired pattern is that of case 1. The best obtained patterns by CODE-EIG are depicted in Figure 6. We notice that the three patterns approximate the desired pattern. The peak SLL values are −30.20 dB, −30 dB, and −29.99 dB for N
=
18, 17, and 16, respectively. The element amplitude and phases for the three patterns are presented in Table 5.Table 5
Case  1. Element amplitudes and phases forN
=
18, 17,  and 15 asymmetrically excited isotropic elements.
n
I
n
φ
n (deg)
I
n
φ
n (deg)
I
n
φ
n (deg)
1
0.136
179.895
0.057
0.774
0.139
179.976
2
0.268
180.000
0
177.597
0.120
179.946
3
0.342
149.887
0.108
179.997
0.167
179.975
4
0.370
137.275
0.104
179.914
0.336
179.994
5
0.455
136.661
0.153
179.702
0.572
161.028
6
0.625
134.342
0.277
179.926
0.853
139.058
7
0.868
114.851
0.549
174.629
0.999
107.863
8
0.990
91.554
0.839
144.046
0.972
81.897
9
0.973
66.089
0.985
116.748
0.871
60.612
10
0.833
41.468
1.000
91.827
0.630
42.012
11
0.560
10.548
0.905
67.285
0.430
37.924
12
0.291
0
0.674
47.620
0.411
39.555
13
0.142
2.974
0.472
39.502
0.386
27.793
14
0.108
0.884
0.355
36.326
0.311
0.003
15
0.065
0
0.354
28.654
0.161
0.004
16
0
0.630
0.241
0.021
0.036
2.615
17
0.048
177.224
0.139
0.099
18
0.070
135.981Figure 6
Case  1. Radiation patterns obtained forN
=
18, 17, and 16 asymmetrically excited isotropic elements.The fourth example presents the design case 2 as the desired pattern obtained by an asymmetrically excited linear array of 20 and 18 isotropic elements. The author has not been able to find in the literature an example that uses this design case asymmetrically excited. The best obtained patterns by CODE-EIG are illustrated in Figure7. We notice that the 20-element array matches better the desired pattern than the 18-element one. The peak SLL values are −25.01 dB and −24.99 dB for N
=
20 and N
=
18 elements, respectively. Table 6 holds the corresponding element amplitudes and phases for both arrays.Table 6
Case  2. Element amplitudes and phases forN
=
20  and 18 asymmetrically excited isotropic elements.
n
I
n
φ
n (deg)
I
n
φ
n (deg)
1
0.155
179.993
0.132
179.993
2
0.297
167.574
0.182
179.93
3
0.356
133.836
0.246
177.114
4
0.227
104.702
0.392
174.602
5
0.081
179.853
0.552
155.925
6
0.460
179.979
0.643
134.614
7
0.809
146.989
0.670
117.929
8
1.000
114.094
0.766
107.713
9
0.980
86.945
0.926
90.237
10
0.844
68.696
0.980
63.686
11
0.757
56.667
0.807
31.301
12
0.698
40.979
0.463
0
13
0.595
20.449
0.120
0.454
14
0.413
5.048
0.226
69.332
15
0.282
0.173
0.341
44.786
16
0.244
3.129
0.305
20.131
17
0.196
0
0.194
0.009
18
0.109
0.070
0.108
0.022
19
0.070
24.731
20
0.094
25.046Figure 7
Case  6. Radiation patterns obtained forN
=
20 and 18 asymmetrically excited isotropic elements.


## 3.2. Conformal Array Synthesis

The final design case presents a conformal cylindrical array design. This case has been solved in [4] using simulated annealing. In [7] another design is obtained using linear programming. We assume a cylindrical array consisting of MXN elements in the ϕ and z directions, respectively. For this case M
=
N
=
8. The radius R
0 of the cylindrical surface is 5
λ, while the cylinder height is h
0
=
3.5
λ. The angular region of the array extends from -
ϕ
0 to ϕ
0, where ϕ
0
=
20
°. The array elements are uniformly spaced by 0.5
λ.The far-field radiation pattern of the conformal cylindrical array is given by(15)
FF
ϑ
,
ϕ
,
I
-
,
φ
-
=
∑
m
=
1
M
∑
n
=
1
N
I
m
n
×
e
l
ϑ
,
ϕ
×
e
j
φ
m
n
+
2
π
/
λ
x
m
n
sin
⁡
ϑ
cos
⁡
ϕ
+
y
m
n
sin
⁡
ϑ
sin
⁡
ϕ
+
z
m
n
cos
⁡
ϑ
,where m, n are the element indices in the ϕ and z directions, respectively, and e
l
(
ϑ
,
ϕ
) is the element pattern given by [7](16)
e
l
ϑ
,
ϕ
=
sin
⁡
ϑ
cos
⁡
ϕ
,
for
ϕ
-
ϕ
m
≤
π
2
,
0
,
otherwise
,where ϕ
m is the element azimuthal position.The positions of the array elements are(17)
x
m
n
=
R
0
cos
⁡
ϕ
m
,
y
m
n
=
R
0
sin
⁡
ϕ
m
,
z
m
n
=
-
h
0
2
+
n
×
0.5
λ
.The desired far-field pattern, F
F
d, has a cosecant square dependency in the ϕ
=
0
° plane for an elevation angle ranging between ϑ
s
t
a
r
t
c
o
s and ϑ
e
n
d
c
o
s. The setting of the half-power beamwidth (HPBW) in the ϑ
=
90
° plane to a desired value, H
P
B
W
d
ϑ
=
90
°, within a tolerance percentage is also a requirement. Additionally, the peak sidelobe levels outside the shaped beam region must remain below a predefined desired value, P
S
L
L
d.For this case, we express the objective function as(18)
F
I
-
,
φ
-
=
1
K
1
∑
k
=
1
K
1
F
F
d
ϑ
k
,
0
,
I
-
,
φ
-
-
F
F
c
ϑ
k
,
0
,
I
-
,
φ
-
2
+
Ξ
·
max
⁡
0
,
PSLL
d
-
PSLL
c
+
Ξ
·
max
⁡
0
,
HPBW
c
ϑ
=
90
°
-
HPBW
d
ϑ
=
90
°
HPBW
d
ϑ
=
90
°
-
T
H
P
B
W
,where F
F
c is the calculated far-field pattern in the ϕ
=
0
° plane, K
1 is the number of pattern samples taken for the ϕ
=
0
° plane, P
S
L
L
c is the calculated peak sidelobe level in a set of ϕ
-
c
u
t
s, H
P
B
W
c
ϑ
=
90
° is the calculated half-power beamwidth in the ϑ
=
90
° plane, and T
H
P
B
W is the desired beamwidth tolerance.The design parameters we have selected for this case areϑ
s
t
a
r
t
c
o
s
=
100
°, ϑ
e
n
d
c
o
s
=
140
°, H
P
B
W
d
ϑ
=
90
°
=
16
°, T
H
P
B
W
=
0.01, and P
S
L
L
d
=
-
27 dB (instead of −25 dB in [4, 7]). We set the population size to 100 vectors and the maximum number of iterations is set to 3000. Figures 8(a), 8(b), and 8(c) depict the best obtained radiation pattern with CODE-EIG for the ϕ
=
0
°, ϑ
=
90
°, and a 3D pattern, respectively. One may notice from the ϕ
=
0
° plane pattern that the pattern found has PSLL value of −27.01 dB, which is about 3 dB below the one from [4], while it is about 0.4 dB higher than the one from [7]. Moreover, the ϑ
=
90
° plane pattern has a HPBW about 16
° and the PSLL value in this pattern is below the corresponding values from [4, 7]. Table 7 holds the amplitudes and phases of the best obtained pattern found by CODE-EIG.Table 7
Excitation coefficients for the best conformal design obtained by CODE-EIG. The phase values are in degrees.
m
1,8
2,7
3,6
4,5
n
I
m
n
φ
m
n
I
m
n
φ
m
n
I
m
n
φ
m
n
I
m
n
φ
m
n
1
0.27
179.14
0.15
168.91
0.33
96.05
0.28
104.38
2
0.43
−61.11
0.39
163.30
0.61
147.49
0.74
130.96
3
0.57
−116.48
0.51
−161.67
1.06
−157.97
1.08
155.87
4
0.81
−20.42
1.04
−106.93
1.31
−133.54
1.37
−137.05
5
0.24
49.47
0.69
−53.72
0.87
−50.12
1.05
−108.07
6
0.31
−34.05
0.23
−6.12
0.63
−123.44
0.58
−83.90
7
0.26
50.05
0.53
−76.22
0.57
−33.05
0.64
−92.52
8
0.08
37.82
0.25
30.14
0.36
−23.14
0.47
−46.99Figure 8
Conformal array radiation far-field patterns in (a)ϕ
=
0
∘ plane, (b) ϑ
=
9
0
∘ plane, and (c) 3D pattern.
(a)
(b)
(c)


## 3.3. CODE-EIG Tuning Parameter Sensitivity

The main control parameter of CODE-EIG is the eigenvector ratioP. In order to explore the CODE-EIG sensitivity to control parameter selection we evaluated the first example using different eigenvector ratio settings. The number of elements of the pattern is N
=
20. We set the number of generations to 1000 and the population size to 100. We run the algorithm for 50 independent trials for each case. The best value, the worst value, the mean, and the standard deviation at the last generation are presented here. Table 8 shows the results. CODE-EIG performance is relatively robust to different P values, as long as P is 0.5 or 0.8. Higher or lower eigenvector ratio values result in slightly worse statistical values. This conclusion is in agreement with [24], where the authors show the eigenvector ratio value of 0.5 produces better results.Table 8
Comparative results for symmetrically excited design Case  1 forN
=
20 using different eigenvector ratio settings. The smaller values are in bold.
P
Best
Worst
Mean
St. dev.
0.5
34619.53
37476.33
36030.69
1121.19
0.2
34647.83
37471.51
36779.67
1368.40
0.4
34655.49
37502.64
36573.75
1308.63
0.8
34641.35
37528.71
36510.02
1294.57
0.9
34652.01
37527.25
36806.05
1387.19


## 4. Conclusion

We have presented a novel design approach for shaped beam pattern synthesis with multiple constraints based on Differential Evolution. We have introduced a new DE algorithm, namely, the CODE-EIG, which extends the Composite DE with the eigenvector-based crossover operator. The CODE-EIG algorithm is a powerful and robust global optimizer, which selects the appropriate strategy for trial vector generation and control parameters from a predefined candidate pool. An additional control parameter that is included in CODE-EIG is the eigenvector ratio. Thus, CODE-EIG requires only the adjustment of three parameters: the population size, the number of iterations, and the eigenvector ratio.One of the DE advantages is that very few control parameters have to be adjusted in each algorithm run. However, the control parameters involved in DE are highly dependent on the optimization problem. Moreover, the selection of the appropriate strategy for trial vector generation requires additional computational time using a trial-and-error search procedure. Therefore, it is not always an easy task to tune these parameters. The CODE-EIG algorithm self-adapts these parameters. Although the CODE-EIG algorithm requires more computational time than the DE/rand/1/bin strategy, significant time is saved from using the trial-and-error search procedure.We have designed shaped beam patterns of linear and conformal arrays, which have proven the effectiveness of this approach. The numerical results show that compared with the other DE strategies the CODE-EIG converges faster and produces in average better results. However, as with all evolutionary algorithms there is no assurance for locating the global optimum. In our future work, we plan to study possible modifications to CODE-EIG and apply the algorithm to additional antenna design problems.



---

*Source: 295012-2015-09-06.xml*