# Conformal Array Pattern Synthesis and Activated Elements Selection Strategy Based on PSOGSA Algorithm

Authors: Bin Sun; Chunheng Liu; Yang Liu; Xiaofang Wu; Yongzhen Li; Xuesong Wang
Journal: International Journal of Antennas and Propagation
 (2015)

---

## Abstract

The pattern synthesis and activated element selection for conformal array is investigated based on hybrid particle swarm optimization-gravitational search algorithm (PSOGSA) in this paper. With the introduction of PSOGSA algorithm which is a novel hybrid optimization technique, the element excitations are optimized to obtain the desired pattern for conformal array in the case of considering uncoupled and coupled element pattern. Numerical simulation and full-wave electromagnetic calculation verify the advantage and efficiency of our method. Then, a novel strategy of activated element selection based on PSOGSA algorithm is proposed for saving the energy consumption in conformal array.


---

## Body

## 1. Introduction

Conformal antenna, due to its wide range of scan angle and little influence on aerodynamic performance of the aerial vehicle, is vastly used in radar, communication, navigation, and remote sensing system [1]. In the past years, antenna pattern synthesis has been an attractive research topic in the antenna community. Many classical synthesis methods of linear and planar array can not be applied to conformal array, since the pattern synthesis of conformal array can not be directly deduced from the product of element factor and array factor. Until now, many iterative numerical computing techniques have been done focusing on the pattern synthesis method for conformal array, such as intersection projector method (APM) [2, 3], least mean square (LMS) [4, 5], and adaptive array method [6].Intellectual evolutionary optimization algorithms are proved to be effective in solving antenna array pattern synthesis problems, because of high robust performance and effective search ability. Many algorithms have been proposed and applied to conformal array synthesis, such as genetic algorithm (GA) [7, 8], simulated annealing algorithm (SA) [9, 10], particle swarm optimization algorithm (PSO) [11, 12], differential evolution algorithm (DE) [13], and Invasive Weed Optimization (IWO) [14]. However, the above independent algorithm often shows some drawbacks of complex operation, locally convergence or slow convergence speed. Recently, several hybrid algorithms have been proposed for conformal array design, which inherits advantages from its constituent algorithm, could overcome the drawbacks, and has more powerful search ability, for example, hybrid IGAIPSO algorithm [15] and hybrid IWOPSO algorithm [16]. However, these algorithms scarcely consider the mutual coupling effect, which is a significant factor in the practical pattern synthesis for antenna array. They also did not dispose element selection problem for low hardware and energy consumption.In this paper, we introduce the PSOGSA algorithm into the pattern synthesis with and without coupled element pattern for conformal array. Moreover, a novel antenna element selection strategy is proposed based on PSOGSA algorithm.Gravitational search algorithm (GSA) is a novel optimization computing technique based on the basic principle of the Newton’s law, first proposed by Rashedi et al. [17]. Then, by combining the ability of exploration in GSA algorithm and the ability of exploration in PSO algorithm, Seydali and Siti developed the hybrid PSOGSA algorithm [18], which shows a better ability to find global optimum with fast convergence and be applied in many fields [19].The remaining sections of this paper are organized as follows. The pattern synthesis model and its cost function for conformal array are present in Section2. The principle of PSOGSA is present in Section 3. Numerical comparison and full-wave calculation for conformal array pattern synthesis are addressed in Section 4. Conformal array element selection strategy is addressed in Section 5. Finally, brief conclusions are discussed in Section 6.


## 2. Conformal Antenna Array Pattern Synthesis

Let us consider a circular antenna array composed ofN elements shown in Figure 1, in which the far field radiation pattern can be expressed as(1)Fθ,φ,θ0,φ0=∑n=0N-1An·fnθ,φ·ej2π/λRsin⁡θcos⁡φ-φn-sin⁡θ0cos⁡φ0-φn,where R is the circular radius, θ0 and φ0 are the direction of the main beam in elevation and azimuth angle, and fn(θ,φ) is the element pattern. An=anejϕn, and an and ϕn are the amplitude and phase of the excitation coefficient, respectively.Figure 1
Far field coordinates for the circle array antenna.The goal of the pattern synthesis for conformal array is to optimize the element excitations for minimizing the difference between the synthesized maximum side lobe level (MSLLsyn) and the desired maximum side lobe level (MSLLdesired) under certain main lobe first null beam width (FNBW). Thus, the cost function defined in this paper is the absolute values of excess between the synthesized MSLLsyn and the desired MSLLdesired under the FNBW constraint as follows: (2)min.C=MSLLsyn-MSLLdesireds.t.FNBW=ΦFNBW.


## 3. Hybrid PSOGSA Algorithm and Its Numerical Simulation

### 3.1. Hybrid PSOGSA Algorithm

In this part, hybrid PSOGSA algorithm is introduced and applied in the pattern synthesis of conformal array. The hybrid algorithm mainly integrates the ability of exploration in PSO algorithm and the ability of exploration in GSA algorithm to find optimum solution more efficiently, which has a better balance between the ability of exploitation and exploration efficiently to find the global optimum. The process of the hybrid algorithm is formulated in detail as follows [18].Step 1 (initiate all the agents randomly).
Suppose a system withN agents; the algorithm starts with randomly placing all agents in search space.Step 2 (evaluate the fitness value for each agent).
All the agents are ranked based on their fitness. Thegbest is the agent with the best fitness value.Step 3 (computing the gravitational force andG(t)).
The gravitational forces from agentj on agent i at a specific time t, is defined as follows:(3)Fijdt=GtMpit×MajtRijt+εxjdt-xidt,where Maj(t) is the active gravitational mass related to agent j, Mpi(t) is the passive gravitational mass related to agent i, ε is a small constant, G(t) is the gravitational constant at time t, and Rij(t) is the Euclidian distance between two agents i and j:(4)Rijt=Xit,Xjt2.
G
(
t
) is calculated as(5)Gt=G0×exp⁡-α×itermax⁡iter,where α and G0 are descending coefficient and initial value respectively, iter is the current iteration, and max⁡iter is the maximum number of iterations.Step 4 (the masses and total force computation).
Inertial masses are calculated according their fitness value, as follows:(6)Mai=Mpi=Mii=Mi,i=1,2,…,N,mit=fiti⁡t-worst⁡tbest⁡t-worst⁡t,Mit=mit∑j=1Nmjt,where fiti⁡t is the fitness of agent i at time t.
For a minimization problem,worst⁡t and best⁡t are defined as follows:(7)best⁡t=minj∈1,2,…,N⁡fitj⁡t,worstt=maxj∈1,2,…,N⁡fitj⁡t.
For a maximization problem,worst(t) and best⁡t are defined as follows:(8)best⁡t=maxj∈1,2,…,N⁡fitj⁡t,worst⁡t=minj∈1,2,…,N⁡fitj⁡t.
If the dimension of the problem isd, the total force that acts on agent i is calculated as the following equation:(9)Fidt=∑j=1,j≠iNrandj⁡Fijdt,where randj is a random number in the interval [0,1].Step 5.
According to the law of motion, the acceleration of an agent is proportional to the result force and inverse to its mass; therefore, the acceleration of all agents is defined as follows:(10)aidt=FidtMit,where Mi(t) is the mass of object i.Step 6.
The velocity and position of agents are calculated as follow:(11)vidt+1=randi×vidt+c1×rand×aidt+c2×rand×gbest-xidt,xidt+1=xidt+vidt+1.Step 7.
Renew the fitness value for all agents. The fitness values of renewed agents are calculated; then, process returns to Step2, repeated until either the maximum number of iteration is reached or the fitness is met.The hybrid algorithm flow chart is summarized as in Figure2.Figure 2
Hybrid PSOGSA algorithm flow chart.


### 3.2. Numerical Results Based on Benchmark Functions

In order to show the efficiency PSOGSA algorithm, some tests are executed as follows. A typical bench mark function is applied to compare the efficiency between hybrid and single algorithm:(12)fx=14000∑i=1nxi2-∏i=1ncos⁡xii+1,xi∈-10,10,where f(x) has an amount of local minima and a global minimum zero. It is a tough work to solve the multimodal optimization problem, and the above function is a good one to test the efficiency of proposed algorithm.The optimized parameters are configured as Table1.Table 1
Parameters value of PSO, GSA, and PSOGSA.
Parameters
PSO
GSA
PSOGSA
i
t
e
r
m
a
x
1000
1000
1000
dim
30
30
30
P
m
a
x
40
—
40
x
i
n
i
[−100, 100]
[−100, 100]
[−100, 100]
c
1
2
—
2
c
2
2
—
2
w
m
a
x
0.9
—
0.9
w
m
i
n
0.4
—
0.4Figure3 presents the simulation of three algorithms, which show that the hybrid algorithm is superior to the single algorithm. The hybrid has a better balance between the ability of globe exploration and local exploitation, which could avoid trapping in local minima and capable of enhanced search ability.Figure 3
Testing results of Rastrigin function.


## 3.1. Hybrid PSOGSA Algorithm

In this part, hybrid PSOGSA algorithm is introduced and applied in the pattern synthesis of conformal array. The hybrid algorithm mainly integrates the ability of exploration in PSO algorithm and the ability of exploration in GSA algorithm to find optimum solution more efficiently, which has a better balance between the ability of exploitation and exploration efficiently to find the global optimum. The process of the hybrid algorithm is formulated in detail as follows [18].Step 1 (initiate all the agents randomly).
Suppose a system withN agents; the algorithm starts with randomly placing all agents in search space.Step 2 (evaluate the fitness value for each agent).
All the agents are ranked based on their fitness. Thegbest is the agent with the best fitness value.Step 3 (computing the gravitational force andG(t)).
The gravitational forces from agentj on agent i at a specific time t, is defined as follows:(3)Fijdt=GtMpit×MajtRijt+εxjdt-xidt,where Maj(t) is the active gravitational mass related to agent j, Mpi(t) is the passive gravitational mass related to agent i, ε is a small constant, G(t) is the gravitational constant at time t, and Rij(t) is the Euclidian distance between two agents i and j:(4)Rijt=Xit,Xjt2.
G
(
t
) is calculated as(5)Gt=G0×exp⁡-α×itermax⁡iter,where α and G0 are descending coefficient and initial value respectively, iter is the current iteration, and max⁡iter is the maximum number of iterations.Step 4 (the masses and total force computation).
Inertial masses are calculated according their fitness value, as follows:(6)Mai=Mpi=Mii=Mi,i=1,2,…,N,mit=fiti⁡t-worst⁡tbest⁡t-worst⁡t,Mit=mit∑j=1Nmjt,where fiti⁡t is the fitness of agent i at time t.
For a minimization problem,worst⁡t and best⁡t are defined as follows:(7)best⁡t=minj∈1,2,…,N⁡fitj⁡t,worstt=maxj∈1,2,…,N⁡fitj⁡t.
For a maximization problem,worst(t) and best⁡t are defined as follows:(8)best⁡t=maxj∈1,2,…,N⁡fitj⁡t,worst⁡t=minj∈1,2,…,N⁡fitj⁡t.
If the dimension of the problem isd, the total force that acts on agent i is calculated as the following equation:(9)Fidt=∑j=1,j≠iNrandj⁡Fijdt,where randj is a random number in the interval [0,1].Step 5.
According to the law of motion, the acceleration of an agent is proportional to the result force and inverse to its mass; therefore, the acceleration of all agents is defined as follows:(10)aidt=FidtMit,where Mi(t) is the mass of object i.Step 6.
The velocity and position of agents are calculated as follow:(11)vidt+1=randi×vidt+c1×rand×aidt+c2×rand×gbest-xidt,xidt+1=xidt+vidt+1.Step 7.
Renew the fitness value for all agents. The fitness values of renewed agents are calculated; then, process returns to Step2, repeated until either the maximum number of iteration is reached or the fitness is met.The hybrid algorithm flow chart is summarized as in Figure2.Figure 2
Hybrid PSOGSA algorithm flow chart.


## 3.2. Numerical Results Based on Benchmark Functions

In order to show the efficiency PSOGSA algorithm, some tests are executed as follows. A typical bench mark function is applied to compare the efficiency between hybrid and single algorithm:(12)fx=14000∑i=1nxi2-∏i=1ncos⁡xii+1,xi∈-10,10,where f(x) has an amount of local minima and a global minimum zero. It is a tough work to solve the multimodal optimization problem, and the above function is a good one to test the efficiency of proposed algorithm.The optimized parameters are configured as Table1.Table 1
Parameters value of PSO, GSA, and PSOGSA.
Parameters
PSO
GSA
PSOGSA
i
t
e
r
m
a
x
1000
1000
1000
dim
30
30
30
P
m
a
x
40
—
40
x
i
n
i
[−100, 100]
[−100, 100]
[−100, 100]
c
1
2
—
2
c
2
2
—
2
w
m
a
x
0.9
—
0.9
w
m
i
n
0.4
—
0.4Figure3 presents the simulation of three algorithms, which show that the hybrid algorithm is superior to the single algorithm. The hybrid has a better balance between the ability of globe exploration and local exploitation, which could avoid trapping in local minima and capable of enhanced search ability.Figure 3
Testing results of Rastrigin function.


## 4. Application of Hybrid PSOGSA Algorithm in Conformal Array

### 4.1. Circular Array Pattern Synthesis without Mutual Coupling Effect

Now consider a half-circular conformal array consisting of 25 elements as shown in Figure4(a), in which all elements located along the circular arc with the radius R=12λ/π. The arc distance between neighbor elements is λ/2.Figure 4
(a) Half-circular array antenna. (b) Microstrip patch antenna.
(a)
(b)The conformal element antenna is microstrip patch antenna, which has a center frequency of 5 GHz and a bandwidth exceeding 10% for return loss less than −10 dB shown in Figure4(b). The pattern is shown in Figure 5(a) calculated via Ansoft HFSS. The unsynthesized pattern of conformal array is presented in Figure 5(b) without considering the mutual coupling effect.Figure 5
(a) Element antenna pattern (dB). (b) conformal array pattern.
(a)
(b)The desired pattern mask is set as follows. The main lobe range is inφM∈[78∘,102∘], and MSLLdesired is set −28 dB. The parameters of optimization algorithms are set in Table 2.Table 2
Parameter setting for each algorithm.
PSO
GSA
IWOPSO
PSOGSA
Param.
Val.
Param.
Val.
Param.
Val.
Param.
Val.
Pop_size
100
Pop_size
100
Pop_size
100
Pop_size
100
C
1
2
G
0
1
Initial_SD
100
C
1
2
C
2
2
—
—
Final_SD
0.005
C
2
2
—
—
—
—
s_max
5
G
0
1
—
—
—
—
s_min
2
—
—
—
—
—
—
C
1
=
C
2
2
—
—The array pattern is synthesized by PSO [11], GSA [17], IWOPSO [16], and PSOGSA algorithm, respectively. After running 1000 times iteration and 10 independent experiments, the average element antenna excitation amplitudes are obtained in Figure 6. Apparently, the amplitude for PSOGSA algorithm is smoother than the other algorithms.Figure 6
Element amplitude distribution for half-circle array.The average convergence curves of cost function value for all algorithms are illustrated in Figure7(a), which indicates that the cost of PSOGSA has the lowest absolute error between the synthesized MSLLsyn and the desired MSLLdesired compared with the other algorithms. Figure 7(b) presents the synthesized array radiation pattern for circular array based on all algorithms. As shown in the figure, the MSLLsyn of PSOGSA algorithm is −27.5 dB, the MSLLsyn of PSO is −26.9 dB, the MSLLsyn of GSA is −27.19 dB, and the MSLLsyn of IWOPSO is −27.18 dB. Thus, PSOGSA algorithm is the most advantageous compared with the other algorithms in the synthesized pattern.Figure 7
(a) Cost function versus iteration times. (b) Radiation pattern for half-circle array.
(a)
(b)


### 4.2. Circular Array Pattern Synthesis with Coupled Element Pattern

In order to obtain the practical pattern of the conformal antenna array, mutual coupling effect has to be considered. Here, 13th element located atφ=90∘ is selected as the representation of the coupled element pattern. The full-wave calculation result via Ansoft HFSS is show in Figure 8(a).Figure 8
(a) Embedded element antenna pattern for 13th element. (b) Simulating each element pattern under coupling effect.
(a)
(b)Compared with Figure5(a), the mutual coupling effect makes the element antenna pattern gain decrease and main lobe with ripples and back lobe level increase as shown in Figure 8(a). And the whole conformal array pattern consisting of the coupled element pattern is present in Figure 8(b).By considering mutual coupling effect, the desiredMSLLdesired is modified as −21 dB; the main lobe is in φM∈[78∘,102∘]. With the same array configuration as in Figure 4, the array pattern was synthesized using coupled element pattern based on PSOGSA algorithm.After the PSOGSA optimization, the excitation amplitude and phase obtained in coupled situation are shown in Table3.Table 3
Each element amplitude and phase derived from MGSA-PSO.
Element sequence
Amplitude
Phase
1
0.0937
360.0000
2
0.0456
180.4226
3
0.0010
3.9177
4
0.1105
193.5056
5
0.3185
32.1019
6
0.1103
242.4684
7
0.2314
107.1652
8
0.4748
348.5075
9
0.5022
248.5256
10
0.4412
168.9301
11
0.7887
111.0829
12
0.5743
75.9739
13
0.9999
64.2038
14
0.7664
75.9739
15
0.7340
111.0829
16
0.7172
168.9301
17
0.5780
248.5256
18
0.4458
348.5075
19
0.3180
107.1652
20
0.2627
242.4684
21
0.2734
32.1019
22
0.0011
193.5056
23
0.0703
3.9177
24
0.0010
180.4226
25
0.0414
360.0000The average convergence curves of the cost function derived from GSA, IWOPSO, and PSOGSA algorithm are illustrated in Figure9(a), and the array pattern synthesized with the coupled element pattern is shown in Figure 9(b).Figure 9
(a) Cost function versus iteration times with coupled pattern. (b) Radiation pattern for the half-circular array with coupled pattern.
(a)
(b)According to Figure9(a), it is indicated that PSOGSA algorithm has the lowest absolute error between the synthesized MSLLsyn and the desired MSLLdesired. According to Figure 9(b), the MSLLsyn of PSOGSA algorithm is −20.0 dB, the MSLLsyn of PSO is −19.6 dB, the MSLLsyn of GSA is −19.7 dB, and the MSLLsyn of IWOPSO is −19.8 dB. Thus, PSOGSA algorithm is the most advantageous compared with the other algorithms in the synthesized pattern. In addition, the mutual coupling effect makes the MSLLsyn increase about 7dB compared with Figure 7(b).


### 4.3. Full-Wave Simulation Experiment via Ansoft HFSS

To further more testify the validation of PSOGSA in conformal array pattern synthesis, full-wave electromagnetic calculation is executed via Ansoft HFSS. The result of electromagnetic calculation is nearly the same as numerical simulation using the excitations from Table3 according to Figure 10(a), which demonstrates the efficiency and accuracy of the PSOGSA algorithm. The little difference between Ansoft HFSS calculation and numerical method is due to the element pattern difference. Clearly, the patterns of 13th element and the 25th element are not the same as shown in Figure 10(b). The 13th element is coupled with two-side 12 elements; however, the 25th element is coupled mostly with the right side 12 elements. We have neglected this when synthesizing the array antenna pattern using PSOGSA algorithm and just consider all element pattern the same as 13th element.Figure 10
(a) Comparison from Ansoft HFSS and numerical simulation. (b) Difference between 13th and 25th pattern.
(a)
(b)


## 4.1. Circular Array Pattern Synthesis without Mutual Coupling Effect

Now consider a half-circular conformal array consisting of 25 elements as shown in Figure4(a), in which all elements located along the circular arc with the radius R=12λ/π. The arc distance between neighbor elements is λ/2.Figure 4
(a) Half-circular array antenna. (b) Microstrip patch antenna.
(a)
(b)The conformal element antenna is microstrip patch antenna, which has a center frequency of 5 GHz and a bandwidth exceeding 10% for return loss less than −10 dB shown in Figure4(b). The pattern is shown in Figure 5(a) calculated via Ansoft HFSS. The unsynthesized pattern of conformal array is presented in Figure 5(b) without considering the mutual coupling effect.Figure 5
(a) Element antenna pattern (dB). (b) conformal array pattern.
(a)
(b)The desired pattern mask is set as follows. The main lobe range is inφM∈[78∘,102∘], and MSLLdesired is set −28 dB. The parameters of optimization algorithms are set in Table 2.Table 2
Parameter setting for each algorithm.
PSO
GSA
IWOPSO
PSOGSA
Param.
Val.
Param.
Val.
Param.
Val.
Param.
Val.
Pop_size
100
Pop_size
100
Pop_size
100
Pop_size
100
C
1
2
G
0
1
Initial_SD
100
C
1
2
C
2
2
—
—
Final_SD
0.005
C
2
2
—
—
—
—
s_max
5
G
0
1
—
—
—
—
s_min
2
—
—
—
—
—
—
C
1
=
C
2
2
—
—The array pattern is synthesized by PSO [11], GSA [17], IWOPSO [16], and PSOGSA algorithm, respectively. After running 1000 times iteration and 10 independent experiments, the average element antenna excitation amplitudes are obtained in Figure 6. Apparently, the amplitude for PSOGSA algorithm is smoother than the other algorithms.Figure 6
Element amplitude distribution for half-circle array.The average convergence curves of cost function value for all algorithms are illustrated in Figure7(a), which indicates that the cost of PSOGSA has the lowest absolute error between the synthesized MSLLsyn and the desired MSLLdesired compared with the other algorithms. Figure 7(b) presents the synthesized array radiation pattern for circular array based on all algorithms. As shown in the figure, the MSLLsyn of PSOGSA algorithm is −27.5 dB, the MSLLsyn of PSO is −26.9 dB, the MSLLsyn of GSA is −27.19 dB, and the MSLLsyn of IWOPSO is −27.18 dB. Thus, PSOGSA algorithm is the most advantageous compared with the other algorithms in the synthesized pattern.Figure 7
(a) Cost function versus iteration times. (b) Radiation pattern for half-circle array.
(a)
(b)


## 4.2. Circular Array Pattern Synthesis with Coupled Element Pattern

In order to obtain the practical pattern of the conformal antenna array, mutual coupling effect has to be considered. Here, 13th element located atφ=90∘ is selected as the representation of the coupled element pattern. The full-wave calculation result via Ansoft HFSS is show in Figure 8(a).Figure 8
(a) Embedded element antenna pattern for 13th element. (b) Simulating each element pattern under coupling effect.
(a)
(b)Compared with Figure5(a), the mutual coupling effect makes the element antenna pattern gain decrease and main lobe with ripples and back lobe level increase as shown in Figure 8(a). And the whole conformal array pattern consisting of the coupled element pattern is present in Figure 8(b).By considering mutual coupling effect, the desiredMSLLdesired is modified as −21 dB; the main lobe is in φM∈[78∘,102∘]. With the same array configuration as in Figure 4, the array pattern was synthesized using coupled element pattern based on PSOGSA algorithm.After the PSOGSA optimization, the excitation amplitude and phase obtained in coupled situation are shown in Table3.Table 3
Each element amplitude and phase derived from MGSA-PSO.
Element sequence
Amplitude
Phase
1
0.0937
360.0000
2
0.0456
180.4226
3
0.0010
3.9177
4
0.1105
193.5056
5
0.3185
32.1019
6
0.1103
242.4684
7
0.2314
107.1652
8
0.4748
348.5075
9
0.5022
248.5256
10
0.4412
168.9301
11
0.7887
111.0829
12
0.5743
75.9739
13
0.9999
64.2038
14
0.7664
75.9739
15
0.7340
111.0829
16
0.7172
168.9301
17
0.5780
248.5256
18
0.4458
348.5075
19
0.3180
107.1652
20
0.2627
242.4684
21
0.2734
32.1019
22
0.0011
193.5056
23
0.0703
3.9177
24
0.0010
180.4226
25
0.0414
360.0000The average convergence curves of the cost function derived from GSA, IWOPSO, and PSOGSA algorithm are illustrated in Figure9(a), and the array pattern synthesized with the coupled element pattern is shown in Figure 9(b).Figure 9
(a) Cost function versus iteration times with coupled pattern. (b) Radiation pattern for the half-circular array with coupled pattern.
(a)
(b)According to Figure9(a), it is indicated that PSOGSA algorithm has the lowest absolute error between the synthesized MSLLsyn and the desired MSLLdesired. According to Figure 9(b), the MSLLsyn of PSOGSA algorithm is −20.0 dB, the MSLLsyn of PSO is −19.6 dB, the MSLLsyn of GSA is −19.7 dB, and the MSLLsyn of IWOPSO is −19.8 dB. Thus, PSOGSA algorithm is the most advantageous compared with the other algorithms in the synthesized pattern. In addition, the mutual coupling effect makes the MSLLsyn increase about 7dB compared with Figure 7(b).


## 4.3. Full-Wave Simulation Experiment via Ansoft HFSS

To further more testify the validation of PSOGSA in conformal array pattern synthesis, full-wave electromagnetic calculation is executed via Ansoft HFSS. The result of electromagnetic calculation is nearly the same as numerical simulation using the excitations from Table3 according to Figure 10(a), which demonstrates the efficiency and accuracy of the PSOGSA algorithm. The little difference between Ansoft HFSS calculation and numerical method is due to the element pattern difference. Clearly, the patterns of 13th element and the 25th element are not the same as shown in Figure 10(b). The 13th element is coupled with two-side 12 elements; however, the 25th element is coupled mostly with the right side 12 elements. We have neglected this when synthesizing the array antenna pattern using PSOGSA algorithm and just consider all element pattern the same as 13th element.Figure 10
(a) Comparison from Ansoft HFSS and numerical simulation. (b) Difference between 13th and 25th pattern.
(a)
(b)


## 5. The Strategy of the Active Elements Selection for Conformal Array Antenna

### 5.1. Active Elements Selection Strategy

For beam scanning angle is specified, only parts of elements make great contribution to the pattern synthesis, and the others contribute so little that can be neglected. In this situation, the strategy of selecting active elements becomes a crucial issue.In this paper, a novel strategy is proposed for active elements selection based on PSOGSA algorithm, which is shown in Figure11.Figure 11
Element selection strategy.Step 1.
Firstly, all the elements are optimized to synthesize the desired array pattern using PSOGSA algorithm. Then, the normalization is performed for each element. Turn off the element excitations whose coefficients are smaller thant(0≤t≤1), and keep the other element excitation activated.Step 2.
For the active elements from Step1, reoptimize the element excitation for the desired pattern and renormalization. Remove the excitation whose coefficient is below t and save the others.Step 3.
Repeat Steps1 and 2 until all the remainder excitation coefficients are bigger than t or reach the specified least element number.Now the conformal array configured as Figure4 is used to illustrate the performance of our active elements selection strategy. For practical phased antenna array, several scan angles should be formed; however, we only take φ=90∘ as design example. Let t=0.1; the beam pattern is synthesized with the strategy proposed, and the entire optimization progress is as Table 4.Table 4
Active elements selection process.
Element sequence
Phi = 90weight ①
Step  1selectionweight ②
Step  2Reoptimization  weight ③
Step  1Selection weight ④
Step  2Reoptimization weight ⑤
Satisfy end condition
1
0.0261
—
—
—
—
—
2
0.0010
—
—
—
—
—
3
0.0299
—
—
—
—
—
4
0.0739
—
—
—
—
—
5
0.1502
0.1502
0.0928
—
—
—
6
0.2184
0.2184
0.1529
0.1527
0.1090
0.1090
7
0.3324
0.3324
0.2562
0.2524
0.1973
0.1973
8
0.4723
0.4723
0.4015
0.3935
0.3418
0.3418
9
0.6027
0.6027
0.5494
0.5411
0.5108
0.5108
10
0.7519
0.7519
0.7165
0.7121
0.7039
0.7039
11
0.8832
0.8832
0.8578
0.8549
0.8549
0.8549
12
1.0000
1.0000
0.9751
0.9773
0.9824
0.9824
13
1.0000
1.0000
0.9883
1.0000
0.9919
0.9919
14
1.0000
1.0000
0.9545
0.9828
0.9524
0.9524
15
0.9170
0.9170
0.8551
0.8816
0.8073
0.8073
16
0.7816
0.7816
0.6937
0.7278
0.6340
0.6340
17
0.6192
0.6192
0.5397
0.5702
0.4537
0.4537
18
0.4948
0.4948
0.3766
0.3999
0.2870
0.2870
19
0.3306
0.3306
0.2421
0.2571
0.1666
0.1666
20
0.2140
0.2140
0.1350
0.1456
0.0780
0.0780
21
0.1289
0.1289
0.0766
—
—
—
22
0.0520
—
—
—
—
—
23
0.0291
—
—
—
—
—
24
0.0010
—
—
—
—
—
25
0.0207
—
—
—
—
—The active elements of each step are shown in Figure12.Figure 12
(a) All elements. (b) Active elements selection in Step1. (c) Active elements selection in Step 2.
(a)
(b)
(c)Figure13(a) presents the pattern synthesized result by weight ①, weight ③, and weight ⑤ based on PSOGSA algorithm. Apparently, the patterns obtained from 25 elements, 17 elements, and 15 elements are nearly the same, which means using only 60% of total elements can obtain nearly the same pattern from the whole array elements activated.Figure 13
(a) Pattern synthesized by weight derived from weight ①, weight ③, and weight ⑤ in Table4. (b) Comparison between Ansoft HFSS and numerical results.
(a)
(b)In order to validate the strategy proposed, full-wave electromagnetic calculation is executed using 15 active elements as shown in Figure13(b). We can see that the array pattern obtained from Ansoft HFSS and numerical method is nearly the same, which demonstrates the accuracy and availability of our selection strategy.


## 5.1. Active Elements Selection Strategy

For beam scanning angle is specified, only parts of elements make great contribution to the pattern synthesis, and the others contribute so little that can be neglected. In this situation, the strategy of selecting active elements becomes a crucial issue.In this paper, a novel strategy is proposed for active elements selection based on PSOGSA algorithm, which is shown in Figure11.Figure 11
Element selection strategy.Step 1.
Firstly, all the elements are optimized to synthesize the desired array pattern using PSOGSA algorithm. Then, the normalization is performed for each element. Turn off the element excitations whose coefficients are smaller thant(0≤t≤1), and keep the other element excitation activated.Step 2.
For the active elements from Step1, reoptimize the element excitation for the desired pattern and renormalization. Remove the excitation whose coefficient is below t and save the others.Step 3.
Repeat Steps1 and 2 until all the remainder excitation coefficients are bigger than t or reach the specified least element number.Now the conformal array configured as Figure4 is used to illustrate the performance of our active elements selection strategy. For practical phased antenna array, several scan angles should be formed; however, we only take φ=90∘ as design example. Let t=0.1; the beam pattern is synthesized with the strategy proposed, and the entire optimization progress is as Table 4.Table 4
Active elements selection process.
Element sequence
Phi = 90weight ①
Step  1selectionweight ②
Step  2Reoptimization  weight ③
Step  1Selection weight ④
Step  2Reoptimization weight ⑤
Satisfy end condition
1
0.0261
—
—
—
—
—
2
0.0010
—
—
—
—
—
3
0.0299
—
—
—
—
—
4
0.0739
—
—
—
—
—
5
0.1502
0.1502
0.0928
—
—
—
6
0.2184
0.2184
0.1529
0.1527
0.1090
0.1090
7
0.3324
0.3324
0.2562
0.2524
0.1973
0.1973
8
0.4723
0.4723
0.4015
0.3935
0.3418
0.3418
9
0.6027
0.6027
0.5494
0.5411
0.5108
0.5108
10
0.7519
0.7519
0.7165
0.7121
0.7039
0.7039
11
0.8832
0.8832
0.8578
0.8549
0.8549
0.8549
12
1.0000
1.0000
0.9751
0.9773
0.9824
0.9824
13
1.0000
1.0000
0.9883
1.0000
0.9919
0.9919
14
1.0000
1.0000
0.9545
0.9828
0.9524
0.9524
15
0.9170
0.9170
0.8551
0.8816
0.8073
0.8073
16
0.7816
0.7816
0.6937
0.7278
0.6340
0.6340
17
0.6192
0.6192
0.5397
0.5702
0.4537
0.4537
18
0.4948
0.4948
0.3766
0.3999
0.2870
0.2870
19
0.3306
0.3306
0.2421
0.2571
0.1666
0.1666
20
0.2140
0.2140
0.1350
0.1456
0.0780
0.0780
21
0.1289
0.1289
0.0766
—
—
—
22
0.0520
—
—
—
—
—
23
0.0291
—
—
—
—
—
24
0.0010
—
—
—
—
—
25
0.0207
—
—
—
—
—The active elements of each step are shown in Figure12.Figure 12
(a) All elements. (b) Active elements selection in Step1. (c) Active elements selection in Step 2.
(a)
(b)
(c)Figure13(a) presents the pattern synthesized result by weight ①, weight ③, and weight ⑤ based on PSOGSA algorithm. Apparently, the patterns obtained from 25 elements, 17 elements, and 15 elements are nearly the same, which means using only 60% of total elements can obtain nearly the same pattern from the whole array elements activated.Figure 13
(a) Pattern synthesized by weight derived from weight ①, weight ③, and weight ⑤ in Table4. (b) Comparison between Ansoft HFSS and numerical results.
(a)
(b)In order to validate the strategy proposed, full-wave electromagnetic calculation is executed using 15 active elements as shown in Figure13(b). We can see that the array pattern obtained from Ansoft HFSS and numerical method is nearly the same, which demonstrates the accuracy and availability of our selection strategy.


## 6. Conclusions

In this paper, PSOGSA algorithm is efficiently applied to pattern synthesis and activated elements election for conformal array. The theory of PSOGSA algorithm is explained and simulated based on Bench function. By considering uncoupled and coupled element pattern, PSOGSA algorithm is applied in conformal pattern synthesis with numerical comparison and full-wave calculation. Then, an element selection strategy is proposed in conformal array with nearly the same radiation pattern. The results of numerical simulation and full-wave calculation demonstrate that PSOGSA algorithm has excellent ability of global optimum search for the conformal pattern synthesis and the active elements selection.



---

*Source: 858357-2015-08-11.xml*