# A Low Complexity 2D Pattern Synthesis Algorithm for Cylindrical Array

Authors: Chao Liu; Zhizhong Ding; Xiaoping Liu
Journal: International Journal of Antennas and Propagation
 (2013)

---

## Abstract

This paper proposes a 2D pattern synthesis algorithm for cylindrical array. According to the geometric characteristic of cylinder, we can regard a cylindrical array as an equivalent linear array whose elements are identical circular subarrays. Therefore, the beam pattern can be obtained by the product of the array factor of linear array and beam pattern of circular subarray. Then, the 2D beamforming can be realized by two 1D beamforming processes. We can prove that the complex excitation vector of a cylindrical array is the Kronecker product of linear array’s weight vector and circular array’s weight vector. By this algorithm of decomposition and reconstruction, the computational complexity of 2D beamforming could be significantly reduced. Finally, simulation results further illustrate the validity of the proposed method.


---

## Body

## 1. Introduction

Conformal arrays can be commodiously mounted on a curved surface of platforms, such as aerial vehicles and fighter aircraft, and provide aerodynamic shape compatible with the corresponding fuselage. This conformal design leads to an excellent aerodynamic performance and thus has an extremely wide range of applications due to its advantages of low radar cross section (RCS), large surveillance coverage, and volume saving [1]. Although conformal arrays are widely used in many areas, the pattern synthesis of such arrays is still a challenge because of the complexity induced by the unusual configuration. For the arrays, the simple characteristics associated with linear or planar arrays do not hold, and most of the conventional methods would not work.In the recent years, pattern synthesis of conformal arrays has attracted increasing attention and a wide variety of techniques have been developed for the pattern synthesis of conformal arrays. In [2], Bucci et al. propose a method to solve the pattern synthesis problem of arbitrary planar array. They alternately found the projection on the upper and lower bounds of the desired pattern to control the shape of main lobe and the side lobe level (SLL). Afterwards, the method was extended to conformal arrays. Considering the embedded element radiation pattern, Steyskal used it in a conformal wing array [3]. Tseng and Griffiths provide a new iterative method in [4] to achieve the desired beam pattern of arbitrary arrays by controlling the side lobe peaks. But the number of side lobe peaks that it can control is limited. In [5], Vaskelainen presents an iterative least-squares method by assigning different weight values to different directions. He further presents a modified least-squares optimization method with linear constraints in [6] to obtain a prescribed shape of main lobe. In [7], Olen and Compton propose another method based on adaptive array theory. In this method, they first set excessive artificial interference signals in observation area outside the main lobe and then adaptively adjust the intensities of these interference signals to control the side lobe level. In [8], Zhou and Ingram improve the method above. They only allow iterations occurred in the main lobe and among side lobe peaks. The method can provide more convenient main lobe shape control, but finding the location of side lobe peaks might increase computational complexity. In [9], Guo et al. use linearly constrained minimum variance (LCMV) criterion to achieve the desired beam pattern with a flat top main lobe. Dohmen et al. give a synthesis method of conformal array in [10] to design both the copolarized and cross-polarized patterns. In [11], Zou et al. present an adaptive beamforming method with low level of polarization components based on geometric algebra. The synthesis problem of an arbitrary array antenna can be seen as a general optimization problem. Therefore, intelligent optimization algorithms can also be used to solve this problem. Genetic algorithms (GA) are applied to complete this kind of optimization in [12, 13]. Simulated annealing and particle swarm optimization (PSO) also have good performance of global optimization and are used in the area of pattern synthesis [14, 15].Most of these methods are applied to general conformal arrays and need an iterative process. The computational complexity is acceptable in 1D condition. However, the amount of calculation will significantly increase when they are applied to 2D beamforming. Some researchers have realized this problem and began to find fast algorithms [16]. Nevertheless, the structure information is not taken into consideration. Most of the conformal arrays have regular structures which may be helpful to the array beamforming. In [17], the authors indicate that the cylindrical array can be treated as a linear array while beamforming. However, they did not touch upon the 2D beamforming problem.In this paper, we present a new method to solve the 2D beam pattern synthesis for cylindrical arrays. According to the geometry feature of cylinder, we treat the cylindrical array as an equivalent linear array whose elements are identical circular subarrays. Since the linear array and circular array are orthogonal to each other, 2D beam pattern of the array in elevation and azimuth direction will be mainly affected by linear and circular arrays, respectively. Therefore, the 2D beamforming process can be realized by two individual 1D pattern syntheses of the circular and linear arrays successively. Through this process, the new method can greatly decrease the amount of computation of 2D beamforming for cylindrical arrays.


## 2. Problem Formulation

Consider a cylindrical array consisting ofM×N elements; the sketch map can be seen in Figure 1. M and N are the number of elements contained in each line and circle. O(x,y,z) is the global Cartesian coordinate system, while O(x′,y′,z′) is the local system, and ai is the central angle of the ith element. The far-field beam pattern in the generic direction (θ,φ) can be written as
(1)E(θ,φ)=∑i=1MN‍wigi(θi′,φi′)exp{jkViTe(θ,φ)},
in which θ and φ are the azimuth and elevation, wi is the complex excitation voltage of the ith element, k=2π/λ is the phase constant, and gi(θi′,φi′) is the antenna gain of the ith element with its local direction (θi′,φi′) associated with the global direction (θ,φ). Vi=[xi,yi,zi]T is the position vector of the ith element, while e(θ,φ)=[cosθcosφ,sinθcosφ,sinφ]T is the unit vector in the direction (θ,φ).Figure 1
Geometry configuration and array coordinate system of cylindrical array.Since the array is conformal to a curved surface, the antenna elements generally direct their radiation beams toward different directions. Therefore, the transformation between the global coordinate system and the element local coordinate system needs to be carried out to calculate the contribution of each element to the whole conformal array radiation, which is shown as follows:(2)(θ,φ)⟺(x,y,z)⟺(x′,y′,z′)⟺(θ′,φ′).The Euler rotation matrix is a very useful tool for this spatial rotation transformation [18]. The corresponding rotation matrix can be written as
(3)R(D,E,F)=R1R2R3=[cosFsinF0-sinFcosF0001][cosE0-sinE010sinE0cosE]×[cosDsinD0-sinDcosD0001],
where [D,E,F] are three Euler rotation angles of the axes Oz, Oy′′, and O′z′, respectively. And then, the spatial rotation transformation can finally be written as
(4)[cosθ′cosφ′sinθ′cosφ′sinφ′]=R(D,E,F)[cosθcosφsinθcosφsinφ].For a cylindrical array, the maximal radiation of antenna is along the normal direction at most of the time; the coordinate transformation process can be expressed as(5)O(x,y,z)⟹O(x′′,y′′,z)⟹O(x′,y′,z′).Without considering the coordinate shift, the rotation angles in (3) are D=ai,E=π/2,andF=0 in this case, where ai is the corresponding central angle of the ith element.


## 3. 2D Beamforming for Cylindrical Arrays

In a planar array, the array factor and the element pattern are separable, but this condition does not hold for the case of a general conformal array. Consider the structure characteristic of cylindrical array; it consists of a series of identical circular arrays and the whole array can be seen as an equivalent linear array whose elements are these identical circular subarrays. Therefore, the array beam pattern can still be obtained by the principle of pattern multiplication. And then, the 2D beam pattern of cylindrical arrayE(θ,φ) can be expressed by the following formula:
(6)E(θ,φ)=E1(θ,φ)E2(θ,φ),
where E1(θ,φ) is the array factor of the linear array and E2(θ,φ) is the beam pattern of the circular subarray which is used as the element pattern here. Assume that the complex excitation vector of linear array and circular array are w1=[w11,w12,…,w1M]T and w2=[w21,w22,…,w2N]T, respectively. Then, the beam pattern E1(θ,φ) and E2(θ,φ) can be written as
(7)E1(θ,φ)=w1Ta1(θ,φ),(8)E2(θ,φ)=w2T[diag{g(θ,φ)}a2(θ,φ)],
where a1(θ,φ) and a2(θ,φ) are the ideal steering vector of the linear array and the circular array and g(θ,φ)=[g1(θ1′,φ1′),…,gN(θN′,φN′)]T. The symbol diag{g} represents a diagonal matrix with diagonal entries g1,…,gN. Define the equivalent steering vector a2′(θ,φ)=diag{g(θ,φ)}a2(θ,φ), and from (6), (7), and (8), beam pattern in (1) can be rewritten as follows:
(9)E(θ,φ)=(w1Ta1(θ,φ))(w2Ta2′(θ,φ))=(w1T⊗w2T)(a1(θ,φ)⊗a2′(θ,φ))=wTa(θ,φ),
in which a(θ,φ)=a1(θ,φ)⊗a2′(θ,φ) is the steering vector of the whole cylindrical array and its weight vector can be obtained as follows:
(10)w=w1⊗w2,
where “⊗" denotes the Kronecker product.Since the array factor of linear array is identical in azimuth, it can be simplified asE1(φ). Moreover, since the linear array and circular array are orthogonal to each other, the beam pattern in the azimuth direction is mainly determined by the circular subarrays. Therefore, the 2D beam pattern can be decomposed into two 1D beamforming processes, one in the azimuth direction and one in the elevation direction. In order to achieve the desired 2D beam pattern Ed(θ,φ), we first carry out beamforming for circular subarray to achieve the given SLL in azimuth and then obtain the weight vector w2 which can satisfy the following equation:
(11)w2Ta2′(θ,φ=φd)=E2(θ,φ=φd)=Ed(θ,φ=φd).And then,E2(θ=θd,φ) is used as the element pattern in the beamforming for linear array in the direction of elevation to obtain the weight vector w1 which can satisfy the following equation:
(12)w1Ta1(φ)E2(θ=θd,φ)=E1(φ)E2(θ=θd,φ)=Ed(θ=θd,φ).The weight vector for the 1D pattern synthesis can be obtained based on an adaptive array method as follows [7]:
(13)w=μRi+n-1a*(ϕd)=μ(σn2I+∑j=1J‍σj2ajajH)-1a*(ϕd),
where ϕd is the look direction, while σn2 and σj2,j=1,2,…,J, are the power of noise and artificial interference signals, respectively. Through the following iteration, the optimal weight vector can be obtained to approach the desired beam pattern:
(14)σj2(n+1)={0ϕjinmain loberegionmax{0,Γj(n)}ϕjinside loberegion,
where n denotes the nth iteration, Γj(n)=σj2(n)+α[En(ϕj)-Ed(ϕj)], and Ed(ϕ) is the desired beam pattern.Since the radiation of array element is not isotropic, some elements may have little contribution to the look direction. And consequently we ignore these elements in the beamforming process to further decrease the computational complexity. In most cases, elements radiate outwards around the normal direction, so we only use the elements in an area around the look direction, and the subarray becomes a circular arc array as shown in Figure2. The solid ones represent the used elements and the circle ones represent the ignored elements.Figure 2
Elements used in the beamforming for circular array.From the above, the proposed algorithm can be summarized as follows.(1)
Consider a cylindrical array consisting ofM×N identical elements with their identical pattern g(θ,φ). Select the contributive elements as shown in Figure 2, and carry out the beamforming process for the circular arc array in azimuth with the fixed elevation φ=φd. Through the iteration of (13) and (14), the weight vector w2 of circular subarray and its 1D beam pattern E2(θ,φd) can be obtained.(2)
Calculate the beam patternE2(θd,φ) of circular subarray by (8). Use it as the element pattern in the M-element linear array and similarly carry out the beamforming process in elevation with the fixed azimuth θ=θd.(3)
After gettingw2 and w1 for the circular and linear arrays, respectively, the weight vector w of the cylindrical array can be obtained by (10), and then the 2D beam pattern E(θ,φ) can be obtained by (9) as well.From the above, we can see that the practically used array is a part of a complete cylindrical array. Actually, the proposed method can be applied to any conformal array with similar structure that it can be seen as an equivalent linear array whose elements are identical subarrays.


## 4. Results and Discussions

In this section, simulations are provided to illustrate the effectiveness of the proposed method. Consider a uniform cylindrical array withM=8 and N=25 as shown in Figure 1. The space between neighboring elements is half wavelength. Each element is assumed to have the same pattern function g(θi,φi)=sinφi,θi∈[0,2π], and φi∈[0,π/2] in their respective local coordinate system. Since the attenuation of the pattern is more than 10 dB at φi=15° in the local coordinate system and this angle is θ=75° in the global coordinate system, the contributive element we choose to carry out the uniform circular array beamforming is in the area θ∈[θd-75°,θd+75°] and others are ignored with this element pattern. Without loss of generality, we assume that the look direction is (θd,φd)=(0°,0°). Then, the number of actually used elements in the uniform circular array beamforming is 11 in following simulations. The desired beam pattern is defined by the main lobe beamwidth between first nulls (50°) and side lobe level (−40 dB). The interested angle area is [-90°,90°]×[-90°,90°] with its scan step 1°.First, we carry out beamforming for the 11-element arc array. The artificial interferences are set with an interval of 3° and the same below. After 100 iterations, the weight vectorw2 is obtained and its beam pattern E2(θ,φ=0°) is shown in Figure 3. Consequently, E2(θ=0°,φ) can be obtained and is used as the element pattern for the following 8-element uniform linear array beamforming in the elevation direction. The process is similar to the above and iteration number is also set to be 100. Figure 4 shows the final 1D beam pattern in the elevation direction.Figure 3
1D beam pattern in the azimuth direction, withφ=0°.Figure 4
1D beam pattern in the elevation direction, withθ=0°.After getting the weight vectorsw1 and w2, the whole weight vector w of the practically used conformal array can be calculated and so is the 2D beam pattern as shown in Figure 5. Figures 6 and 7 show its projections in azimuth and elevation directions, respectively. Obviously, they coincide with the 1D beam pattern we obtained above, which well supports the analyses above. It can be seen from the figures that the finally obtained beam pattern can successfully achieve the prescribed beam pattern response.Figure 5
2D beam pattern for the cylindrical array.Figure 6
Beam pattern projection in the azimuth direction.Figure 7
Beam pattern projection in the elevation direction.The method in [9] is also performed for the 2D beamforming for cylindrical array in the same condition. With the same interval, there will be 61×61=3721 artificial interferences. Result can be seen in Figure 8 which is similar to the above. However, the amount of computation is more than that of the proposed method. Figure 9 shows the normalized weights obtained by the two methods, where method 1 denotes the proposed method and method 2 denotes the method in [9].Figure 8
2D beam pattern for the cylindrical array.Normalized weights of the cylindrical array.
(a)
(b)In the pattern synthesis, the calculation mainly focuses on the spatial rotation transformation, the update of the power of artificial interferences, and matrix inversion. In our method, the whole process is mainly composed of two 1D beamforming processes. There are 61 interferences for the beamforming of each of azimuth and elevation. Assume that the number of iterations is L. Then, to get the final optimal weight vectorw, we need to calculate rotation transformation 11×61 times, matrix inversion 2×L times and interference power 2×61×L times. Moreover, the size of matrix inversion is small (11×11 and 8×8). Meanwhile, if the method in [9] is directly applied to the 2D pattern synthesis of this cylindrical array, then we need to calculate rotation transformation 88×3721 times, matrix inversion L times and interference power 3721×L times. The size of matrix inversion is 88×88. The number of artificial interferences and size of matrix are the main reason for the large amount of calculation. Our method can significantly reduce both of them and hence reduce the amount of calculation. Moreover, with the increase of array element number, the amount of computation will be rapidly increasing for those direct 2D beamforming algorithms, while our method will not. More scenarios are simulated to verify the proposed algorithm. With the Intel (R) Core (TM) i5 CPU 2.4 GHz and 4 GB RAM, the elapsed time of getting the optimal weight vector for both the proposed method and the method in [9] is shown in Table 1.Table 1
CPU time for different array sizes.
Array size (M×N)
T
1 (s)
T
2 (s)
8
×
25
0.671006
268.723427
12
×
36
0.731100
562.419302
20
×
50
1.116190
2756.073839T
1 and T2 in Table 1 are the elapsed time for the proposed method and the method in [9], respectively. Although the elapsed time may be affected by the computing environment and some parameters of the conformal array, there is no doubt that our method has greatly reduced the amount of calculation without performance degradation.


## 5. Conclusion

In this paper, an efficient 2D beamforming method is proposed for cylindrical arrays. According to the geometry characteristics, a cylindrical array can be seen as an equivalent linear array composed of identical circular subarrays, so that the principle of pattern multiplication can be applied and 2D pattern synthesis process can be realized by two 1D pattern syntheses based on the linear and circular arrays, respectively. This method avoids the complicated calculation of 2D beam pattern and can significantly reduce the amount of computation. Furthermore, the proposed method can be applied to any similar conformal arrays which can be seen as a linear array whose elements are identical subarrays. Effectiveness of the algorithm has been illustrated by the above simulations.



---

*Source: 352843-2013-12-04.xml*