# Fast 3D Pattern Synthesis with Polarization and Dynamic Range Ratio Control for Conformal Antenna Arrays

Authors: Massimiliano Comisso; Roberto Vescovo
Journal: International Journal of Antennas and Propagation
 (2014)

---

## Abstract

This paper proposes an iterative algorithm for the 3D synthesis of the electric far-field pattern of a conformal antenna array in the presence of requirements on both the polarization and the dynamic range ratio (DRR) of the excitations. Thanks to the use of selectable weights, the algorithm allows a versatile control of the DRR and of the polarization in a given angular region and requires a low CPU time to provide the array excitations. Furthermore, a modified version of the algorithm is developed to enable the optimization of the polarization state by phase-only control. Numerical results are presented to verify the usefulness of the proposed approach for the joint pattern and polarization synthesis of conformal arrays with reduced or even unitary DRR.


---

## Body

## 1. Introduction

The polarization properties of an antenna array can be exploited to acquire considerable benefits, such as mitigating fading effects, improving the performance of a radar system, and compensating the transmitter/receiver polarization mismatch. To obtain a desired polarization, an approach widely adopted by the antenna developers relies on the design of the physical antenna system, thus operating on the array geometry [1], on the technology of the single array element [2, 3], and on the position of the feeds [4] or, more in general, using combinations of these techniques [5–8]. When the physical antenna system is given, including the elements, the geometry, and the feeding network, the designer can act on the array excitations by properly selecting their amplitudes and phases in order to satisfy the polarization requirements. Accordingly, different solutions have been presented in the literature for minimizing the cross-polar pattern, when a linear polarization is desired [9–11], or for enabling the generation of more specific circular or even elliptical polarizations [12, 13]. Furthermore, in addition to the polarization requirements, the array must be often able to generate a desired pattern in the presence of constraints on the dynamic range ratio (DRR) of the excitations. For the case in which no constraints on the DRR are imposed, the authors of [9] describe an algorithm based on the alternating projection method and the use of an Euler rotation matrix to provide low cross-polarization patterns for conformal arrays. In [10], the array synthesis in the presence of requirements on the cross-polarization level and the side lobe level (SLL) is performed by combining particle swarm optimization and alternate projection methods. A technique for simultaneously synthesizing the copolar and the cross-polar array amplitude patterns in the presence of DRR requirements is proposed in [11]. A least-square optimization approach is adopted in [12] for synthesizing the pattern of a conformal array imposing a given polarization and a specified value for each excitation amplitude. The joint synthesis of an upper bounded pattern and a desired polarization in a given angular region for a conformal array in absence of DRR constraints is addressed in [13] by adopting a convex formulation of the synthesis problem, which is solved by a second-order cone program.This overview reveals that the problem of controlling the DRR together with the polarization of the electric far-field vector of the array is still open. To deal with such problem, this study proposes a versatile method for the 3D pattern synthesis of a conformal array, when DRR reduction and polarization control in a given angular region are simultaneously required. The method extends the approach used in [11], whose application is limited to linear polarizations, to enable the generation of any desired polarization while maintaining a low synthesis CPU time. Moreover, with respect to [11], a modified version of the method is developed to solve the problem of joint pattern and polarization synthesis by phase-only control.The paper is organized as follows. Section2 introduces the problem. Section 3 presents the algorithm and its modification for phase-only control. Section 4 discusses the numerical results. Section 5 summarizes the main conclusions.


## 2. Problem Formulation

In the 3D space, referred to as a spherical coordinate systemO
(
r
,
θ
,
ϕ
), a conformal array of N elements and the column vector of its complex excitations a
=
[
a
1
,
…
,
a
N
]
T is considered, with (
·
)
T denoting the transpose operator. Such excitations generate the electric far-field patterns E
θ
(
a
;
r
^
) and E
ϕ
(
a
;
r
^
), where r
^ is the unit vector of the direction of observation, while E
θ
(
a
;
r
^
) and E
ϕ
(
a
;
r
^
) are the components, along the standard unit vectors θ
^ and ϕ
^, respectively, of the electric far-field vector E
(
a
;
r
^
); that is, E
(
a
;
r
^
)
=
E
θ
(
a
;
r
^
)
θ
^
+
E
ϕ
(
a
;
r
^
)
ϕ
^. Thus, we can write(1a)
E
θ
(
a
;
r
^
)
=
∑
n
=
1
N
‍
a
n
E
~
n
(
r
^
)
e
j
(
2
π
/
λ
)
r
-
n
·
r
^
=
∑
n
=
1
N
‍
a
n
f
~
n
(
r
^
)
=
a
T
f
~
(
r
^
)
,
(1b)
E
ϕ
(
a
;
r
^
)
=
∑
n
=
1
N
‍
a
n
E
^
n
(
r
^
)
e
j
(
2
π
/
λ
)
r
-
n
·
r
^
=
∑
n
=
1
N
‍
a
n
f
^
n
(
r
^
)
=
a
T
f
^
(
r
^
)
, where E
~
n
(
r
^
) and E
^
n
(
r
^
) are the components along θ
^ and ϕ
^, respectively, of the electric far-field vector of the nth array element, r
-
n is the position of the nth element, λ is the wavelength, f
~
(
r
^
)
=
[
f
~
1
(
r
^
)
,
…
,
f
~
N
(
r
^
)
]
T, and f
^
(
r
^
)
=
[
f
^
1
(
r
^
)
,
…
,
f
^
N
(
r
^
)
]
T. Consider now a positive function E
(
r
^
) normalized to unity, and a pair of coefficients (
δ
,
ξ
), with δ
≥
0 and -
π
≤
ξ
≤
π, which specifies a desired polarization for E
(
a
;
r
^
). The objective is to find an excitation vector a such that(2a)
|
E
θ
(
a
;
r
^
)
|
approximates
E
(
r
^
)
;
(2b)
E
(
a
;
r
^
)
is
polarized
according
to
(
δ
,
ξ
)
;
that
is
,
E
θ
(
a
;
r
^
)
=
E
ϕ
(
a
;
r
^
)
δ
e
j
ξ
in
a
given
angular
region
,
(2c)
DRR
(
a
)
=
max
⁡
n
{
|
a
n
|
}
min
⁡
n
{
|
a
n
|
}
is
reduced
.As remarked in [11, 13–15], one can observe that the expressions in (1a) and (1b) include the mutual coupling effects between the array elements, provided that the single-element patterns E
~
n
(
r
^
) and E
^
n
(
r
^
) are numerically calculated by using a suitable software. Precisely, following [15], the patterns E
~
n
(
r
^
) and E
^
n
(
r
^
) can be evaluated by considering the nth element as fed by unity voltage and zero phase, and the other N
-
1 array elements without feeds. Thus, E
~
n
(
r
^
) and E
^
n
(
r
^
) do not identify the isolated element patterns but identify the equivalent element patterns accounting for the coupling with the remaining array elements. The possibility of using this approach is enabled by the linear relationships in ((1a) and (1b)) between the electric far-field patterns of the single-element and those generated by the array [11, 13–15].


## 3. Synthesis Algorithm

To solve the above formulated problem, we introduce the following cost function:(3)
F
(
a
)
=
∫
Θ
‍
|
|
E
θ
(
a
;
r
^
)
|
-
E
(
r
^
)
w
(
r
^
)
|
2
w
(
r
^
)
d
θ
d
ϕ
+
∫
Θ
ML
‍
|
E
θ
(
a
;
r
^
)
-
E
ϕ
(
a
;
r
^
)
δ
e
j
ξ
|
2
v
(
r
^
)
d
θ
d
ϕ
+
b
∑
m
,
n
=
1
N
‍
(
|
a
n
|
-
|
a
m
|
)
2
N
,
where Θ and Θ
ML denote the solid angular regions that identify the entire domain of interest and the main lobe, respectively, w
(
r
^
) and v
(
r
^
) are positive weight functions, and b is a nonnegative weight. Minimizing the first term in (3) allows the amplitude of E
θ
(
a
;
r
^
) to approximate E
(
r
^
). Minimizing the second term allows one to obtain a pattern E
ϕ
(
a
;
r
^
) approximately scaled by 1
/
δ and with a phase shift -
ξ with respect to E
θ
(
a
;
r
^
), thus providing the desired polarization. Minimizing the third term reduces the DRR of the excitations.Now, starting from the desired positive functionE
(
r
^
), let us consider the complex function E
(
r
^
)
e
j
Φ
(
r
^
), where Φ
(
r
^
) is an arbitrary real function. Then, minimizing the first integral in (3) is equivalent to minimizing the new integral ∫
Θ
‍
|
E
θ
(
a
;
r
^
)
-
E
(
r
^
)
e
j
Φ
(
r
^
)
w
-
1
(
r
^
)
|
2
w
(
r
^
)
d
θ
d
ϕ. Note that the function Φ
(
r
^
) minimizing the latter integral coincides with the phase function of E
θ
(
a
;
r
^
). Thus, our problem can be solved by minimizing the functional
(4)
d
2
[
a
,
Φ
(
r
^
)
]
=
∫
Θ
‍
|
E
θ
(
a
;
r
^
)
-
E
(
r
^
)
e
j
Φ
(
r
^
)
w
(
r
^
)
|
2
w
(
r
^
)
d
θ
d
ϕ
+
∫
Θ
ML
‍
|
E
θ
(
a
;
r
^
)
-
E
ϕ
(
a
;
r
^
)
δ
e
j
ξ
|
2
v
(
r
^
)
d
θ
d
ϕ
+
b
∑
m
,
n
=
1
N
‍
(
|
a
n
|
-
|
a
m
|
)
2
N
.
Thus, the power synthesis problem at hand is reformulated as a field synthesis problem consisting in the minimization of (4). The unknowns a and Φ
(
r
^
) minimizing d
2 can be found recalling the procedure used in [11]. As a first step, (4) can be rewritten as
(5)
d
2
[
a
,
Φ
(
r
^
)
]
=
∫
Θ
‍
|
E
θ
(
a
;
r
^
)
|
2
w
(
r
^
)
d
θ
d
ϕ
-
∫
Θ
‍
E
(
r
^
)
[
E
θ
(
a
;
r
^
)
e
-
j
Φ
(
r
^
)
kkkkkkkkkkk
+
E
θ
*
(
a
;
r
^
)
e
j
Φ
(
r
^
)
]
d
θ
d
ϕ
+
∫
Θ
ML
‍
|
E
θ
(
a
;
r
^
)
-
E
ϕ
(
a
;
r
^
)
δ
e
j
ξ
|
2
v
(
r
^
)
d
θ
d
ϕ
+
2
b
N
(
N
∑
n
=
1
N
‍
|
a
n
|
2
-
∑
m
,
n
=
1
N
‍
|
a
n
a
m
|
)
+
Υ
2
,
where (
·
)
* is the complex conjugate and
(6)
Υ
2
=
∫
Θ
‍
E
2
(
r
^
)
w
(
r
^
)
d
θ
d
ϕ
.
Assume now the following notations: (
·
)
H is the Hermitian operator, R
(
·
) is the real part, ∥
·
∥
F is the Frobenius norm, and ∥
·
∥
1 is the sum of the moduli of the elements of a matrix or vector, thus representing the entrywise absolute norm for a matrix and the Manhattan norm for a vector. Substituting (1a) and (1b) into (5) and manipulating, one obtains
(7)
d
2
[
a
,
Φ
(
r
^
)
]
=
a
H
A
a
-
2
R
(
g
H
a
)
-
2
b
N
∥
a
a
H
∥
1
+
2
b
∥
a
∥
F
2
+
Υ
2
,
where A
=
[
A
m
n
], g
=
[
g
1
,
…
,
g
N
]
T, and(8a)
A
m
n
=
∫
Θ
‍
f
~
n
*
(
r
^
)
f
~
m
(
r
^
)
w
(
r
^
)
d
θ
d
ϕ
+
∫
Θ
ML
‍
{
-
δ
[
f
^
n
*
(
r
^
)
f
~
m
(
r
^
)
e
-
j
ξ
+
f
~
n
*
(
r
^
)
f
^
m
(
r
^
)
e
j
ξ
]
kkkkkkkkk
+
f
~
n
*
(
r
^
)
f
~
m
(
r
^
)
+
δ
2
f
^
n
*
(
r
^
)
f
^
m
(
r
^
)
[
f
^
n
*
(
r
^
)
f
~
m
(
r
^
)
e
-
j
ξ
+
f
~
n
*
(
r
^
)
f
^
m
(
r
^
)
e
j
ξ
]
}
v
(
r
^
)
d
θ
d
ϕ
(8b)
g
n
=
∫
Θ
‍
f
~
n
*
(
r
^
)
E
(
r
^
)
e
j
Φ
(
r
^
)
d
θ
d
ϕ which can be minimized iteratively by subdividing each iteration in two parts. In the first part, given Φ
(
r
^
), d
2 is minimized with respect to a by modifying the excitations one at a time. In particular, observing that A
m
n
=
A
n
m
*, (7) can be manipulated to put into evidence its dependence on the generic pth component a
p. Thus, defining ζ
p
=
|
a
p
| and ψ
p
=
arg
(
a
p
), one obtains
(9)
d
2
[
a
,
Φ
(
r
^
)
]
=
2
{
|
T
p
|
cos
⁡
[
ψ
p
-
arg
(
T
p
)
]
-
2
b
N
∥
a
p
∥
1
}
ζ
p
+
[
A
p
p
+
2
b
(
1
-
1
N
)
]
ζ
p
2
+
χ
p
,
where T
p
=
A
p
a
p
-
g
p, with A
p denoting the pth row of the matrix A and a
p
=
[
a
1
,
…
,
a
p
-
1
,
0
,
a
p
+
1
,
…
,
a
N
]
T, and where χ
p is a term independent of a
p. Now, (9) can be minimized with respect to the argument ψ
p and to the amplitude ζ
p. Precisely, imposing that the multiplicative coefficient of ζ
p be minimum, one obtains ψ
p
=
arg
(
T
p
)
+
(
2
κ
+
1
)
π, where κ is an arbitrary integer, while imposing that the derivative of d
2 with respect to ζ
p be zero, one obtains ζ
p
=
(
N
|
T
p
|
+
2
b
∥
a
p
∥
1
)
/
[
N
A
p
p
+
2
b
(
N
-
1
)
]. Combining these two expressions, the value of a
p
=
ζ
p
e
j
ψ
p minimizing (9) can be written in compact form as
(10)
a
p
=
-
2
b
e
j
arg
(
T
p
)
∥
a
p
∥
1
+
N
T
p
N
A
p
p
+
2
b
(
N
-
1
)
.In the second part of each iteration, givena, d
2 is minimized with respect to Φ
(
r
^
). Precisely, by substituting (8b) into (7), one obtains
(11)
d
2
[
a
,
Φ
(
r
^
)
]
=
η
-
2
R
[
∫
Θ
‍
a
T
f
~
(
r
^
)
E
(
r
^
)
e
-
j
Φ
(
r
^
)
d
θ
d
ϕ
]
which, being the term η independent of Φ
(
r
^
), is minimized for
(12)
Φ
(
r
^
)
=
arg
[
a
T
f
~
(
r
^
)
]
.Summarizing, the proposed algorithm for polarization and DRR control develops as follows. Moving from a starting pointa
0, the function Φ
0
(
r
^
) is evaluated by (12). At each successive iteration q
=
1,2
,
…
,
g
q and a
q are updated by (8b) and (10), respectively, for each p
=
1
,
…
,
N, and then Φ
q
(
r
^
) is updated by (12). Since, at each step, both (10) and (12) minimize d
2
(
a
,
Φ
), it results in d
2
[
a
0
,
Φ
0
(
r
^
)
]
≥
d
2
[
a
1
,
Φ
1
(
r
^
)
]
≥
⋯
≥
d
2
[
a
q
,
Φ
q
(
r
^
)
]
≥
⋯. Therefore, the sequence {
d
q
2
}
=
{
d
2
[
a
q
,
Φ
q
(
r
^
)
]
} is nonincreasing and thus convergent. This sequence is terminated when (
d
q
2
-
d
q
+
1
2
)
/
d
q
+
1
2
≤
ε, where ε is a suitable threshold. Note that the proposed method has the considerable advantage of solving a problem characterized by multiple requirements (pattern shaping, polarization selection, and DRR reduction) with a low computational cost, thanks to the availability of the closed forms (10) and (12). However, it is a local search algorithm; thus the final solution may be dependent on the starting point a
0. This latter aspect will be discussed in Section 4 by numerically analyzing the impact of the initial conditions on the final solution.


### 3.1. Phase-Only Control

The above presented algorithm allows one to control the DRR of the array excitations by properly selecting the nonnegative weightb in (4). However, in order to simplify the feeding network, some applications require that the DRR be exactly equal to one, so that the joint pattern and polarization synthesis is realized by phase-only control. To this aim, this subsection proposes a modified version of the algorithm in which all the excitations have the same amplitude.The problem can be reformulated by considering (7) with b
=
0 and a
n
=
ζ
e
j
ψ
n for n
=
1
,
…
,
N, where ζ
=
|
a
1
|
=
⋯
=
|
a
N
|. Thus, the objective becomes that of finding the common excitation amplitude ζ, the column vector of the N excitation phases ψ
=
[
ψ
1
,
…
,
ψ
N
]
T, and the function Φ
(
r
^
), in such a way as to minimize the functional
(13)
d
′
2
[
ζ
,
ψ
,
Φ
(
r
^
)
]
=
ζ
2
Ψ
H
A
Ψ
-
2
ζ
R
(
g
H
Ψ
)
+
Υ
2
,
where Ψ
=
[
e
j
ψ
1
,
…
,
e
j
ψ
N
]
T. An iterative approach can still be adopted to minimize (13), but now subdividing each iteration into three parts. In the first part, given ψ and Φ
(
r
^
), d
′
2 is minimized with respect to ζ. To this aim, one can observe that d
′
2 is a parabolic function of ζ; thus, d
′
2 is minimum for
(14)
ζ
=
R
(
g
H
Ψ
)
Ψ
H
A
Ψ
.
To this purpose, it is important to observe that both quantities Ψ
H
A
Ψ and R
(
g
H
Ψ
) are always positive because a
H
A
a and R
(
g
H
a
) in (7) are positive for each a
∈
C
N. In fact, a
H
A
a is the sum of the first and the third integral in (5), which are positive; furthermore, 2
R
(
g
H
a
) is equal to the second integral in (5), which is positive when recalling that the optimal function Φ
(
r
^
) coincides with the phase of E
θ
(
a
;
r
^
). This ensures that (14) always yields a positive value that minimizes (13).In the second part of each iteration, givenζ and Φ
(
r
^
), d
′
2 is minimized with respect to ψ by modifying the excitation phases one at a time. Precisely, (13) can be manipulated to put into evidence its dependence on the generic pth phase ψ
p, obtaining
(15)
d
′
2
[
ζ
,
ψ
,
Φ
(
r
^
)
]
=
ρ
p
+
2
ζ
|
ν
p
|
cos
⁡
[
ψ
p
-
arg
(
ν
p
)
]
,
where ρ
p is a term independent of ψ
p and ν
p
=
A
p
Ψ
p
-
g
p, with Ψ
p
=
[
e
j
ψ
1
,
…
,
e
j
ψ
p
-
1
,
0
,
e
j
ψ
p
+
1
,
…
,
e
j
ψ
N
]
T. Hence, the functional in (15) is minimum for
(16)
ψ
p
=
arg
(
ν
p
)
+
(
2
κ
+
1
)
π
,
where κ is an arbitrary integer.In the third part of each iteration, givenζ and ψ, d
′
2 is minimized with respect to Φ
(
r
^
). In particular, by substituting (8b) into (13) and identifying as η
′ a term independent of Φ
(
r
^
), d
′
2 can be expressed as
(17)
d
′
2
[
ζ
,
ψ
,
Φ
(
r
^
)
]
=
η
′
-
2
ζ
R
[
∫
Θ
‍
Ψ
T
f
~
(
r
^
)
E
(
r
^
)
e
-
j
Φ
(
r
^
)
d
θ
d
ϕ
]
which is minimum for
(18)
Φ
(
r
^
)
=
arg
[
Ψ
T
f
~
(
r
^
)
]
.The modified algorithm develops as follows. Moving from the starting common amplitudeζ
0 and the starting vector of the excitation phases ψ
0, the function Φ
0
(
r
^
) is evaluated by (18). At each successive iteration q
=
1,2
,
…, once g
q is updated by (8b) for p
=
1
,
…
,
N, the quantities ζ
q, ψ
q, and Φ
q
(
r
^
) are updated by (14), (16) for p
=
1
,
…
,
N, and (18), respectively. The sequence {
d
q
′
2
}
=
{
d
′
2
[
ζ
q
,
ψ
q
,
Φ
q
(
r
^
)
]
}, which is still nonincreasing and hence convergent, is terminated when (
d
q
′
2
-
d
q
+
1
′
2
)
/
d
q
+
1
′
2
≤
ε
′, where ε
′ is a suitable threshold.


## 3.1. Phase-Only Control

The above presented algorithm allows one to control the DRR of the array excitations by properly selecting the nonnegative weightb in (4). However, in order to simplify the feeding network, some applications require that the DRR be exactly equal to one, so that the joint pattern and polarization synthesis is realized by phase-only control. To this aim, this subsection proposes a modified version of the algorithm in which all the excitations have the same amplitude.The problem can be reformulated by considering (7) with b
=
0 and a
n
=
ζ
e
j
ψ
n for n
=
1
,
…
,
N, where ζ
=
|
a
1
|
=
⋯
=
|
a
N
|. Thus, the objective becomes that of finding the common excitation amplitude ζ, the column vector of the N excitation phases ψ
=
[
ψ
1
,
…
,
ψ
N
]
T, and the function Φ
(
r
^
), in such a way as to minimize the functional
(13)
d
′
2
[
ζ
,
ψ
,
Φ
(
r
^
)
]
=
ζ
2
Ψ
H
A
Ψ
-
2
ζ
R
(
g
H
Ψ
)
+
Υ
2
,
where Ψ
=
[
e
j
ψ
1
,
…
,
e
j
ψ
N
]
T. An iterative approach can still be adopted to minimize (13), but now subdividing each iteration into three parts. In the first part, given ψ and Φ
(
r
^
), d
′
2 is minimized with respect to ζ. To this aim, one can observe that d
′
2 is a parabolic function of ζ; thus, d
′
2 is minimum for
(14)
ζ
=
R
(
g
H
Ψ
)
Ψ
H
A
Ψ
.
To this purpose, it is important to observe that both quantities Ψ
H
A
Ψ and R
(
g
H
Ψ
) are always positive because a
H
A
a and R
(
g
H
a
) in (7) are positive for each a
∈
C
N. In fact, a
H
A
a is the sum of the first and the third integral in (5), which are positive; furthermore, 2
R
(
g
H
a
) is equal to the second integral in (5), which is positive when recalling that the optimal function Φ
(
r
^
) coincides with the phase of E
θ
(
a
;
r
^
). This ensures that (14) always yields a positive value that minimizes (13).In the second part of each iteration, givenζ and Φ
(
r
^
), d
′
2 is minimized with respect to ψ by modifying the excitation phases one at a time. Precisely, (13) can be manipulated to put into evidence its dependence on the generic pth phase ψ
p, obtaining
(15)
d
′
2
[
ζ
,
ψ
,
Φ
(
r
^
)
]
=
ρ
p
+
2
ζ
|
ν
p
|
cos
⁡
[
ψ
p
-
arg
(
ν
p
)
]
,
where ρ
p is a term independent of ψ
p and ν
p
=
A
p
Ψ
p
-
g
p, with Ψ
p
=
[
e
j
ψ
1
,
…
,
e
j
ψ
p
-
1
,
0
,
e
j
ψ
p
+
1
,
…
,
e
j
ψ
N
]
T. Hence, the functional in (15) is minimum for
(16)
ψ
p
=
arg
(
ν
p
)
+
(
2
κ
+
1
)
π
,
where κ is an arbitrary integer.In the third part of each iteration, givenζ and ψ, d
′
2 is minimized with respect to Φ
(
r
^
). In particular, by substituting (8b) into (13) and identifying as η
′ a term independent of Φ
(
r
^
), d
′
2 can be expressed as
(17)
d
′
2
[
ζ
,
ψ
,
Φ
(
r
^
)
]
=
η
′
-
2
ζ
R
[
∫
Θ
‍
Ψ
T
f
~
(
r
^
)
E
(
r
^
)
e
-
j
Φ
(
r
^
)
d
θ
d
ϕ
]
which is minimum for
(18)
Φ
(
r
^
)
=
arg
[
Ψ
T
f
~
(
r
^
)
]
.The modified algorithm develops as follows. Moving from the starting common amplitudeζ
0 and the starting vector of the excitation phases ψ
0, the function Φ
0
(
r
^
) is evaluated by (18). At each successive iteration q
=
1,2
,
…, once g
q is updated by (8b) for p
=
1
,
…
,
N, the quantities ζ
q, ψ
q, and Φ
q
(
r
^
) are updated by (14), (16) for p
=
1
,
…
,
N, and (18), respectively. The sequence {
d
q
′
2
}
=
{
d
′
2
[
ζ
q
,
ψ
q
,
Φ
q
(
r
^
)
]
}, which is still nonincreasing and hence convergent, is terminated when (
d
q
′
2
-
d
q
+
1
′
2
)
/
d
q
+
1
′
2
≤
ε
′, where ε
′ is a suitable threshold.


## 4. Numerical Results

The performance of the proposed algorithm is evaluated by considering a conformal array consisting ofN
=
15
×
7 elements located on a cylindrical arc of radius R
=
4.46
λ and aperture α
=
9
0
° (Figure 1). Each element is a quadrifilar helix antenna of radius 0.1
λ and turn spacing 0.2
λ with two turns and axis parallel to the y-axis. The generic helix lying in the kth arc (k
=
1
,
…
,
7) and hth column (h
=
1
,
…
,
15) of the array is placed at the point of cylindrical coordinates [
R
sin
[
(
22
-
h
)
π
/
28
]
,
R
cos
⁡
[
(
22
-
h
)
π
/
28
]
,
(
2
-
k
/
2
)
λ
]. This configuration is implemented in the SuperNEC 2.7 electromagnetic simulator to derive the single-element patterns E
~
n
(
r
^
) and E
^
n
(
r
^
), in order to include the mutual coupling effects by adopting the approach described in [15].Figure 1
Cylindrical arc array ofN
=
15
×
7 quadrifilar helixes.The algorithm is implemented in MATLAB on an Intel Core2 Quad Q9300 @2.50 GHz Sun Ultra 24 workstation using just one core of the processor. A discretization step of 4 degrees in both the elevation and the azimuth domains is adopted; hence the desired directions and beamwidths are selected as multiples of 4 degrees. All results are derived by settingε
=
1
0
-
5 in absence of phase-only control, ε
′
=
1
0
-
7 in presence of phase-only control, w
(
r
^
)
=
[
E
(
r
^
)
]
-
μ with μ
=
6 [11], and v
(
r
^
)
=
c
·
1
0
β with β
=
7 and c taking different values.As a first example, we consider a desired patternE
(
r
^
) having a maximum side-lobe level SLL
=
-
10 dB, and a maximum at [
θ
0
,
ϕ
0
]
=
[
8
8
°
,
8
0
°
], thus involving a desired direction close to the azimuth plane. The desired radiated far-field vector must be right-hand circularly polarized (δ
=
1, ξ
=
-
π
/
2) within the main lobe Θ
ML
=
[
7
6
°
,
10
0
°
]
×
[
7
2
°
,
8
8
°
], and the starting points are a
0
=
[
0
,
…
,
0
]
T in absence of phase-only control, and ζ
0
=
0, ψ
0
=
[
0
,
…
,
0
]
T in presence of phase-only control. Figure 2 presents the 2D cuts, corresponding to θ
0 and ϕ
0, of the power patterns |
E
θ
|
2 (Figures 2(a) and 2(b)) and |
E
ϕ
|
2 (Figures 2(c) and 2(d)), and of the axial ratios with a magnification of the main lobe region (Figures 2(e) and 2(f)), while Table 1 reports the performance of the proposed algorithm for different values of the parameters. We remark that a 3D synthesis has been performed; thus the values in the table refer to the entire 3D domain, while the curves in Figure 2 refer to 2D cuts to maintain readability. The results show that in all the considered cases the pattern synthesized approximates the desired one. It is important to note that, even if the quadrifilar helix is an antenna designed to provide a circular polarization, when several helixes are used to form a conformal array, satisfying the requirements on the shape of the radiated far-field pattern might lead to a deterioration of the desired polarization state (case b
=
c
=
0 in Figures 2(e) and 2(f)). This effect can be considerably reduced using the proposed algorithm. In fact, in the presented application, the polarization within Θ
ML can be controlled by properly selecting the parameter c. Precisely, when c
=
1, the axial ratio inside the main lobe is equal to 0.1 dB and hence a right-hand circular polarization is substantially obtained, while, additionally, when b
=
1, a DRR reduction in the order of 87% is achieved. The results reveal also that the possibility of setting the values of the parameters b and c to match a specific requirement provides a considerable versatility to the proposed method, which furthermore requires a low CPU time, considering that just few seconds have been sufficient to perform the 3D synthesis of an array with N
=
105 elements. The last column of Table 1 shows that the modified algorithm is able to satisfy the requirement on the SLL maintaining an axial ratio of 0.2 dB when the synthesis must be performed by phase only control. In this case, the computational time increases to few minutes because a lower threshold ε
′ is required. This behavior is confirmed by Table 2, which shows the impact of the starting point a
0 on the performance of the presented algorithm for b
=
c
=
1 and of the starting point (
ζ
0
,
ψ
0
) on the performance of the modified algorithm for c
=
0.3. The values in the table, which are written in the form x
-
±
σ, where x
- denotes the mean value and σ the standard deviation, are obtained by performing one hundred realizations. The amplitudes ζ
0 and ζ
p
0, for p
=
1
,
…
,
N, are generated according to a uniform distribution inside the interval [
0,1
], while the phases ψ
p
0, for p
=
1
,
…
,
N, are generated according to a uniform distribution inside the interval [
-
π
,
π
]. The results in the table show, in particular for the case b
=
c
=
1, that the standard deviations referred to as the SLL and the axial ratio are sufficiently low, thus revealing that, for different starting points, even if different acceptable solutions are obtained, they substantially converge to a similar performance.Table 1
Performance of the proposed algorithm.
b
=
0
b
=
0
b
=
1
Phase-only control
c
=
0
c
=
1
c
=
1
withc
=
0.3
SLL [dB]
−15.4
−10.2
−10.0
−10.0
Maximum axial ratio in the main lobe [dB]
2.7
0.1
0.1
0.2
DRR
99.8
60.8
7.9
1.0
CPU time [s]
4.6
9.3
6.1
312.1Table 2
Impact of the starting point on the results.
b
=
1
Phase-only control
c
=
1
withc
=
0.3
SLL [dB]
−10.97 ± 0.81
−8.94 ± 1.26
Maximum axial ratio in the main lobe [dB]
0.17 ± 0.07
0.48 ± 0.30
DRR
11.89 ± 4.36
1.00 ± 0.00
CPU time [s]
6.99 ± 1.86
206.82 ± 52.17Patterns and axial ratios obtained using the proposed algorithm: (a)|
E
θ
|
2 for ϕ
=
ϕ
0
=
8
0
°, (b) |
E
θ
|
2 for θ
=
θ
0
=
8
8
°, (c) |
E
ϕ
|
2 for ϕ
=
ϕ
0
=
8
0
°, (d) |
E
ϕ
|
2 for θ
=
θ
0
=
8
8
°, (e) axial ratio for ϕ
=
ϕ
0
=
8
0
°, and (f) axial ratio for θ
=
θ
0
=
8
8
°.
(a)
(b)
(c)
(d)
(e)
(f)As a second example, we consider a two-dimensional scanning application, which requires the generation of a family of6
×
6
=
36 patterns able to scan the angular region [
8
0
°
,
10
0
°
]
×
[
8
0
°
,
10
0
°
] by steps of 4 degrees, according to the adopted discretization. The generic desired pattern E
i
,
l
(
r
^
) of this family has a maximum at [
θ
0
i
,
ϕ
0
l
]
=
[
8
0
°
+
i
·
4
°
,
8
0
o
+
l
·
4
°
] (i
,
l
=
0
,
…
,
5), and SLL
=
-
10 dB. The desired radiated far-field vector must still be right-hand circularly polarized within the main lobe Θ
ML
i
,
l
=
[
θ
0
i
-
1
2
°
,
θ
0
i
+
1
2
°
]
×
[
ϕ
0
l
-
8
°
,
ϕ
0
l
+
8
°
], and the starting points are a
0
=
[
0
,
…
,
0
]
T in absence of phase-only control and ζ
0
=
0, ψ
0
=
[
0
,
…
,
0
]
T in presence of phase-only control. Table 3 shows the obtained results, still referred to as the 3D domain and reported in the form x
-
±
σ, where the mean value x
- and the standard deviation σ are evaluated by averaging over the 36 patterns necessary to perform the scanning. The values in the table are derived considering four cases: b
=
c
=
0 and b
=
c
=
0.2 in absence of phase-only control; c
=
0 and c
=
0.1 in presence of phase-only control. Besides, to provide a view of the synthesized patterns for the most challenging case, in which the scanning is performed by both polarization and phase-only control, Figure 3 presents the 2D cuts, corresponding to θ
=
8
8
° and ϕ
=
8
8
°, of the power patterns |
E
θ
|
2 (Figures 3(a) and 3(b)) and |
E
ϕ
|
2 (Figures 3(c) and 3(d)).Table 3
Results obtained for the two-dimensional scanning of the region [80°, 100°] × [80°, 100°].
b
=
c
=
0
b
=
c
=
0.2
SLL [dB]
−15.00 ± 0.70
−10.75 ± 2.59
Maximum axial ratio in the main lobe [dB]
1.94 ± 0.59
0.42 ± 0.25
DRR
65.56 ± 49.27
32.44 ± 15.73
CPU time [s]
4.68 ± 1.24
4.48 ± 0.95
Phase-only control
c
=
0
c
=
0.1
SLL [dB]
−14.03 ± 1.12
−9.89 ± 2.44
Maximum axial ratio in the main lobe [dB]
1.71 ± 0.95
0.80 ± 0.29
DRR
1.00 ± 0.00
1.00 ± 0.00
CPU time [s]
131.50 ± 24.01
146.12 ± 32.92Patterns obtained for a two-dimensional scanning of the region[
8
0
°
,
10
0
°
]
×
[
8
°
,
10
0
°
] using the proposed algorithm by phase-only control with c
=
0.1: (a) |
E
θ
|
2 for ϕ
=
8
8
°, (b) |
E
θ
|
2 for θ
=
8
8
°, (c) |
E
ϕ
|
2 for ϕ
=
8
8
°, and (d) |
E
ϕ
|
2 for θ
=
8
8
°.
(a)
(b)
(c)
(d)These results reveal the satisfactory scanning capabilities of the proposed algorithm, which can be maintained also when improvements on the purity of the polarization and reductions of the DRR are required. In particular, one may observe from Figure3 that the main lobe remains properly shaped even when a unitary DRR is imposed. In this regard, it is worth noticing that, during the scanning process, as the desired direction becomes farther from the direction identified by the y-axis, which is parallel to the main axis of each quadrifilar helix that composes the array, more difficulties arise in achieving a low axial ratio. In fact, in absence of polarization requirements, the scanned region can be considerably extended both in the elevation and the azimuth domain, obtaining satisfactory patterns both in absence and in presence of phase-only control. This possibility is confirmed by Table 4, which presents the results derived maintaining, with respect to the previous case, the same starting points and pattern requirements in terms of beamwidth, shape of the main lobe, and SLL, but setting the parameter c equal to zero and considering a two-dimensional scanning of the region [
6
4
°
,
11
6
°
]
×
[
6
4
°
,
11
6
°
] by steps of 4 degrees. The SLL values in the table and the 2D cuts of the power patterns |
E
θ
|
2 and |
E
ϕ
|
2, reported in Figures 4(a) and 4(b) and Figures 4(c) and 4(d), respectively, show that the main lobe remains properly shaped even with a phase-only control approach and a larger scanning region, which almost completely covers the azimuthal aperture of the cylindrical arc array.Table 4
Results obtained for the two-dimensional scanning of the region [64°, 116°] × [64°, 116°].
b
=
0.2
Phase-only control
c
=
0
withc
=
0
SLL [dB]
−13.91 ± 1.37
−12.04 ± 1.82
Maximum axial ratio in the main lobe [dB]
2.72 ± 1.39
2.38 ± 1.80
DRR
33.04 ± 21.10
1.00 ± 0.00
CPU time [s]
4.62 ± 1.18
140.19 ± 34.06Patterns obtained for a two-dimensional scanning of the region[
6
4
°
,
11
6
°
]
×
[
6
4
°
,
11
6
°
] using the proposed algorithm by phase-only control with c
=
0: (a) |
E
θ
|
2 for ϕ
=
8
8
°, (b) |
E
θ
|
2 for θ
=
8
8
°, (c) |
E
ϕ
|
2 for ϕ
=
8
8
°, and (d) |
E
ϕ
|
2 for θ
=
8
8
°.
(a)
(b)
(c)
(d)


## 5. Conclusions

A fast and versatile method for the joint pattern and polarization synthesis of conformal arrays when DRR or phase-only control is required has been presented. The method can be adapted to specific problems by properly selecting the available parameters, whose setting is not critical, as it requires just few attempts. The proposed algorithm can be applied both to arrays designed to support multiple polarizations, in order to select the desired one, and to arrays providing a unique polarization, in order to improve its purity by a proper choice of the excitations. Satisfactory results have also been obtained for two-dimensional scanning applications that involve the simultaneous control of the patterns' shape, of the polarization, and of the DRR.



---

*Source: 106741-2014-07-17.xml*