Abstract

In this paper, an effective synthesis method for sparse uniformly excited concentric ring array (CRA) element positions featuring a minimum sidelobe level (SLL) is presented. This method is based on the chaos sparrow search algorithm (CSSA), which can search for the optimal solution under multiple constraints. By improving the constraint on the number of array elements on each ring, the solution range of the optimal solution is further reduced, and the global search ability of the algorithm is effectively improved by introducing tent mapping into the algorithm to initialize the population. Numerical examples are presented to assess the effectiveness and reliability of the proposed method, showing that it can achieve better results than existing methods in the design of sparse concentric ring array arrangements.

1. Introduction

Concentric ring array (CRA), which has the symmetrical structure of a circular array, ensures that its beam and antenna gain performance are basically maintained in a certain range, and it can relatively achieve mutual coupling balance. These excellent properties have made it the subject of intensive investigations in recent years [14]. In particular, array synthesis patterns with arbitrary element positions have attracted more attention [5, 6]. Compared with a uniform array antenna with the same aperture, the sparse array antenna can achieve the same resolution with fewer array elements; then, the cost of the antenna system is reduced [7]. At the same time, sparse array antennas can achieve a low sidelobe level (SLL) without amplitude weighting. In the synthesis of sparse uniformly excited concentric ring arrays, the location of the element is an important factor in the optimal performance because of its ability to obtain the minimum peak sidelobe level (PSLL), but it also leads to complex nonlinear optimization problems [8, 9].

To solve the nonlinear problem of sparse uniformly excited concentric ring array design, several synthesis techniques have been proposed. A method based on space tapering has been suggested by Willey [10], which can be used by controlling the element density taper to match the amplitude taper of the minimum SLL. For the fast synthesis of CRA, a deterministic approach that takes inspiration from density-tapered techniques to optimize the location of array elements has been presented in [11, 12]. The normalized Taylor amplitude taper of the conventional filled array was considered a distribution function to determine the location of the array elements. However, the results obtained by analytical methods are often not globally optimal, and their applications in array synthesis are mostly limited to the suppression of sidelobes. With the advent of intelligent optimization methods, such as genetic algorithms (GAs) [1316], particle swarm optimization (PSO) [17], and differential evolution (DE) [18], they have been successfully applied to the synthesis of sparse CRA due to their global search capability. A hybrid approach (HA) for CRA synthesizing to suppress the PSLL has been proposed in [19], which is convenient for synthesizing uniform amplitude concentric ring arrays with certain constraints. Although intelligent optimization algorithms have the advantage of solving sparse array nonlinear synthesis and achieving better optimization results, their solving efficiency is low, and it is easy to fall into a local optimum.

The sparrow search algorithm (SSA) imitates the foraging and antipredation behavior of sparrows [20], which has the advantages of high search accuracy, strong versatility, and easy implementation compared with other intelligent algorithms [21, 22]. To improve the global search ability and prevent falling into the local optimal solution, a chaos sparrow search algorithm (CSSA) is proposed by introducing a tent-mapping chaos operator into the SSA, which effectively guarantees the uniformity of the initial population and enhances population diversity. This method achieves the suppression of PSLL by optimizing the position distribution of uniformly excited CRA elements. The simulation results demonstrate the effectiveness and superiority of the proposed method.

The paper is organized as follows: Section 2 formulates the problem, and Section 3 describes the procedure of synthesis. Numerical examples and results are presented in Section 4 to fully demonstrate the effectiveness and superiority of this proposed method. Finally, Section 5 concludes this paper.

2. Synthesis of CRA

In Figure 1, a sparse concentric array with a single element at the center and rings from the center to the edge is shown, and the starting element of each ring is arranged on the X-axis. The radius of each ring is , and the number of elements on the corresponding ring is . Generally, the array factor can be characterized as follows [19]:where is the wave number, , is the operating wavelength, is the angular position of the mth element in the nth ring, , and and are the azimuth and elevation angles, respectively.

Assuming that the diameter of the outer ring is D, the minimum spacing of array elements is . To improve the search efficiency, all elements of the array need to satisfy the minimum spacing constraint between adjacent elements and the maximum aperture constraint:

The array radius vector can be represented as follows [19]:where is the constant part and is the variable part.

After determining the ring radius, the maximum number of array elements on a ring can be expressed by the ring radius and minimum array element spacing as follows [13]:

Theoretically, can be any integer between 1 and , but the scarcity of array elements will result in array directivity and gain losses. Therefore, constraints should be set to the minimum number of array elements on the ring. To make the element distribution of the sparse array satisfy the feature of gradual sparsity from the center to the edge, on the sparse array is expressed aswhere is the sparse coefficient, which is used to control the sparse rate of elements on each ring, is determined by the amplitude distribution on the corresponding full-array ring, i.e., sparse coefficient , and is the minimum number of array elements on the ring. To increase the optimization space of , is set to a random number between and .

The fitness function is defined as follows [19]:where represents the array factor, is the maximum value of the main beam, and is position of the first null point.

3. CSSA Analysis

To improve the efficiency of the convergence performance of the sparse concentric ring array optimization problem while reducing the PSLL, the CSSA is used to solve the synthesis of sparse concentric ring arrays. Since the position updating of the CSSA is jumpy and discontinuous, the local optimum can be effectively avoided. The specific optimization steps are as follows.

Step 1. Population initialization.
Assuming that the SSA involves a group of sparrows with a spatial dimension , the position distribution of sparrows can be expressed asTo improve the global search ability of the algorithm and avoid the decline of population diversity in later iterations, tent mapping is chosen as the initialized population for the chaotic sequence of the optimization algorithm. The chaotic self-mapping of the tent is represented aswhich can be abbreviated asRandom variables are introduced to the tent mapping to prevent the destruction of the randomness, ergodicity, and regularity of chaotic variables:where is the maximum number of iterations.
After the chaotic variable is generated by (10), it is introduced to the solution space of variables, and the initial population is generated aswhere and are the minimum and maximum values of , respectively.

Step 2. Population update based on the SSA.
In the SSA, a discoverer-joiner sparrow population model is established, and some sparrows are randomly selected as guards. The discoverer position update formula is as follows [22]:where is the position of the sparrow of the generation in the dimension, is a uniform random number in (0, 1], is the maximum number of iterations, indicates the warning value, indicates the alert threshold, is a random number satisfying the standard normal distribution, and is a matrix of whose elements are all 1. When , the sparrow population is in a safe state, and the discoverer continues to forage. When , danger is found, and all sparrows must fly to safety immediately.
In addition to discoverers, the remaining sparrows in the population are joiners, and their position update formula is as follows [21]:where is the best position of the discoverers and indicates the global worst position in the current iteration. When , it means that the ith joiner needs to change its search area to continue searching for food. Conversely, it means that the ith joiner converges to the global optimal position and performs random foraging around.
The formula for updating the position of the guard position is as follows [21]:where is the current global best position, is a step control parameter satisfies a distribution of N (0, 1), is a random number between [−1, 1], is a very small constant to avoid zero denominators, is the fitness value of the individual sparrow at present, and and are the best and worst fitness values at present, respectively. The position of these sparrows is random at first. When , the sparrows are vulnerable to predators at the edge of their foraging area. When , it suggests that sparrows at the center of the group are conscious of the danger and demand to approach others quickly to adjust their foraging strategies.

4. Numerical Examples and Results

To validate the superiority and stability of the proposed sparse array synthesis method based on the CSSA, we show three examples of unequally spaced concentric ring arrays with different array apertures, where 100 independent runs for each case are conducted. In all experiments, the antenna elements are rotationally symmetric on the rings, the minimum element spacing is set to 0.5λ, and the CRA’s amplitude is equal. In the algorithm, the proportion factor of the number of discoverers to the population , the guard scale factor , and the safety value in (12).

4.1. Example A:

In the first example, the array including a total of 8 concentric rings with an array radius of 4.98λ is considered, and the distributed antenna elements are equally spaced on each ring to form an annular grid array. The population size is 100, and the number of generations is 1000. An eight-ring uniform array is used as the original reference array, and a Taylor distribution of −30 dB is chosen as the array excitation function; then, the sparse coefficient kn= {0.98, 0.91, 0.81, 0.70, 0.57, 0.46, 0.38, and 0.34}.

Figures 2(a) and 2(b) show the adaptation curves of the SSA and CSSA for 10 independent experiments, respectively. It can be seen that the adaptation curve of the CSSA converges faster and has better convergence than the SSA. The value of the optimal solution obtained from 10 independent experiments using the SSA ranged from −28.21 dB to −29.83 dB, with an average value of −29.08 dB. The value of the optimal solution obtained from 10 independent experiments using the CSSA ranged from −29.08 dB to −30.07 dB, with an average value of −29.77 dB, which is about 0.69 dB lower than that of the SSA. The results show that the CSSA has better stability than the SSA. Although Figure 2(b) shows that nearly 500 iterations are necessary to reach convergence, numerical results suggested that after 500 iterations, performance remains good. Therefore, the maximum number of iterations T = 500 is chosen in the later numerical experiments.

Figure 3(a) shows the optimal element distribution of an acquired array with a single element at the center in 100 independent trials. It has 8 rings with 192 elements. The radius of each ring is 0.50, 1.00, 1.50, 2.00, 2.61, 3.34, 4.17, and 4.98. The corresponding number of elements is 6, 12, 18, 21, 31, 37, 31, and 35. Figure 3(b) shows the radiation pattern in the plane, and the PSLL is −30.08 dB, which is better than when using GA (−22.94 dB) [13], MGA (−23.74 dB) [14], and HA (−29.03 dB) [19]. The 3 dB main beamwidth obtained by the CSSA is 0.1129, which is slightly wider than HA (0.1126) [19]. The radiation pattern in three dimensions is illustrated in Figure 3(c). Finally, the total computation time for 500 iterations to complete a single trial is approximately 512 s on a laptop with 8 GB RAM running MATLAB R2019b.

4.2. Example B:

Another synthesis experiment with an array radius of 4.7λ has been conducted with the proposed CSSA. A seven-ring uniform array is used as the original reference array, and a Taylor distribution of −30 dB is chosen as the array excitation function; then, the sparse coefficient kn = {0.97, 0.89, 0.77, 0.63, 0.50, 0.39, and 0.34}. We compared the CSSA with the GA and HA, respectively. GA [13] obtained the best PSLL of −27.82 dB by only optimizing the ring radius and strictly constraining the element number on the corresponding ring. HA [19] obtained the best PSLL of −27.92 dB by both optimizing the ring radius and the element number on the ring. The best PSLL obtained by the CSSA is −29.16 dB when the number of rings and the total amount of elements  = 161, which is 1.34 dB lower than when using GA and 1.24 dB lower than when using HA. Figure 4(a) shows the best element distribution of an acquired array. The radius of each ring is 0.59, 1.13, 1.66, 2.26, 3.01, 3.80, and 4.70. The corresponding number of elements is 7, 13, 19, 25, 32, 29, and 35. Figure 4(b) depicts the radiation pattern in the plane, while Figure 4(c) depicts the radiation pattern in three dimensions. Finally, the total computation time for this example to complete a single trial is approximately 381 s.

4.3. Example C:

In the third example, the proposed CSSA method has been used to carry out a synthesis experiment with an array radius of 4.3λ. The best array is built with 147 elements distributed in 7 rings, and the best element distribution of the obtained array is shown in Figure 5(a). Figure 5(b) shows the comparison of array radiation patterns between HA and CSSA in the plane. It can be concluded that the CSSA has obtained a 1.13 dB lower PSLL than HA and that the 3 dB main beamwidth obtained by the CSSA is slightly higher than that obtained by HA. Figure 5(c) depicts the corresponding 3-D view of the radiation pattern. Finally, the total computation time for this example to complete a single trial is approximately 377 s.

Finally, Table 1 presents the performance comparison results between the proposed CSSA and some state-of-the-art algorithms when R = 4.98λ, 4.7λ, and 4.3λ, which include the number of rings , the total number of elements , the best PSLL, the direction coefficient, and the 3 dB beamwidth. Table 2 illustrates the radius and the number of elements on the corresponding ring acquired by the CSSA. Numerical examples show that the CSSA is able to obtain lower PSLL values under the same fitness conditions.

Based on the above analysis results, it can be concluded that optimizing the array element positions by the CSSA can obtain a better PSLL suppression effect while maintaining a higher directivity coefficient and a narrower 3 dB beamwidth. In the three simulation examples, the proposed method has preferable efficiency and good convergence. In addition, the comparison of the proposed method with state-of-the-art algorithms demonstrates its superiority in array sparse optimization.

5. Conclusion

We proposed an effective method for the synthesis of sparse CRA with multiple constraints, which effectively achieved SLL suppression of equally uniformly excited CRA. The numerical results demonstrate that the proposed CSSA method achieved a better SLL suppression effect than other methods when the apertures were fixed. Moreover, the proposed method has the potential to solve other array models with multiple constraints, such as rectangular planar arrays and conformal unequally spaced arrays, which have a certain reference value for the optimal design and engineering applications of array antennas in the future.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China under Grant no. 62271011 and the Chongqing Postgraduate Research and Innovation Project under Grant no. CYS21321.