Abstract
To lower the peak sidelobe level (PSLL) of sparse concentric ring arrays, a method with multiple design constraints that embed a function model into modified real genetic algorithm (MGA) and select the grid ring radii as optimization individual to synthesize sparse concentric ring arrays is proposed. The multiple constraints include the array aperture, the minimum element spacing, and the number of elements. The proposed method dynamically calculates the ratio of element on each ring, and it has a faster convergence rate than other algorithms. The MGA uses real number to code the optimization variable, and it reduces the complexity of coding and improves the search efficiency. Finally, the results demonstrate the accuracy and effectiveness of the algorithm.
1. Introduction
Antenna array plays an important role in technology, including mobile communication, radar system, satellite, and medical treatment. For example, in [1], a concentric ring array antenna is designed for reconfigurable isoflux pattern, which uses 61 elements, and the frequency is 2.8 GHz. The antenna arrays can reduce the volume occupation and complex feeding and heat dissipation systems in a satellite antenna system. In [2], a sparse concentric rings array with isoflux coverage to earth surface is designed for low earth orbit (LEO) satellite system, which optimizes the angular separation between the antenna elements and the amplitude excitation across the array to generate the characteristics of the desired isoflux radiation. The proposed antenna arrays can generate a wide isoflux pattern with a good accuracy with respect to the isoflux mask for LEO satellites. In [3], a repeater antenna is designed for wireless body area networks (WBAN) applications, which can obtain the biological information by implanting the device to human-body for treatment, diagnosis, and potential medical applications.
The PSLL is an important criterion to evaluate the performance of the antenna array. Thus, with the multiple design constraints, including the number of elements, the aperture of array, and the minimum spacing of element, it is an important research topic to reduce the PSLL of antenna array in recent years [4]. Study on array synthesis had been carried out for many years and produced many classical methods. These methodologies include Woodward Lawson method and optimization calculation method [5, 6], which can be used to improve pattern performance, and Taylor synthesis method [7], which can be applied to obtain a given pattern. With the development of array antenna, the scale of arrays becomes larger and the structure of arrays becomes complex; the traditional method of array synthesis is no longer applicable. So a lot of intelligence methods were introduced to synthesize array antenna, including genetic algorithm (GA), simulated annealing algorithm (SAA), differential evolution algorithm (DEA), and particle swarm optimization (PSO) [8–11]. There are a few other topics about circular antenna arrays, such as the design of nonuniform circular antenna arrays for low sidelobe [12, 13], the comparison of algorithms about the multiobject design of circular antenna arrays [14], the design of electronically steerable linear arrays [15], and the synthesis of sparse circular antenna arrays for acquiring desired radiation beam pattern [16]. This article selects GA to synthesize the sparse concentric ring arrays for low PSLL. In [17], GA was proposed. In [8], GA was applied to optimize thinned array antenna and acquired efficient result. In [18], a modified real genetic algorithm (MGA) was proposed to synthesize the sparse concentric ring arrays.
Due to the difference of element sparsity degree on each auxiliary ring [19], a function model that presents the relationship between ring radius and the degree of sparsity is proposed in this paper, which is embedded into the process of MGA to synthesize the sparse concentric ring arrays for low PSLL. The proposed method dynamically calculates the ratio of element on each ring, so it has a faster convergence rate than other algorithms. Besides, the MGA uses real number to code the optimization variable, and it reduces the complexity of coding and improves the search efficiency. Comparing the result to that of other algorithms, the proposed method can reduce PSLL efficiently and has a faster convergence.
2. The Optimization Model of Sparse Concentric Ring Array
In a planar array, generate a point randomly; the point is selected as centre point, and then generate several circles with different radius and place some elements on each ring; the resulting planar array is a concentric ring array [20]. The diagram of concentric ring array is shown in Figure 1.
Place an element in the centre of the concentric ring array; the radiation pattern can be described as where = number of elements in ring = number of rings = wavelength = element weights for ring = radius for ring = element location ;
For simplicity, we assume that array antenna meets the ideal conditions, including the element that is isotropic and the element that has uniform excitation amplitude and phase shift. The weight in the same ring is equal, and the elements in each ring are uniform distribution. The direction of main beam of the array points to the normal direction of the array, and the array pattern can be described as where is the sum of element number, is the radius of th element, and is the angle of th element.
If the sparse concentric ring array has multiple constraints, including the number of rings , the maximum array aperture , the minimum element spacing , and the number of elements , the optimization goal is to obtain the minimum PSLL, and the optimal mathematical model can be described as where where is the maximum value of array pattern; , are the regions excluding the main beam; is the element number of th rings. Due to the fact that the elements number of each ring must be an integer, the value of must be rounded up or down to an integer. To meet the requirement of minimum element spacing , the digits to the right of the decimal point are discarded.
3. The Algorithm of Optimization
The flowchart of function model algorithm using MGA is shown in Figure 2.
3.1. Function Model
There is a conclusion in [19]: in the sparse concentric ring array, when the ring is away from the array centre, the more sparsity the element is, the lower the PSLL is. Thus we can propose a function model to present the relation between ring radius and the degree of sparsity. We select the data of [21] as fitting data shown in Table 1.
In Table 1, is the actual element number of th rings, and is the maximum element number th ring can place. is element retention rate of th ring.Different function model is used to fit the data of Table 1, and the result is shown in Table 2.
Plot the fitting data and function model in Figure 3.
In the fitting function model, the is the determination coefficient of fitting degree, and the larger the is, the higher the fitting degree is. In Table 2, the linear model has the maximum , so the linear model is selected as fitting function model. is the radius of circular; the value of is the retention rate of ring with radius . , , are the coefficients of the function model, and they will be different with the variety of the number of elements and the size of the rings. The actual element number of rings with radii can be calculated as follows:
3.2. The Process of MGA
Standard genetic algorithm uses 0,1 binary code for individual, and it has a low efficiency. MGA uses real number code for individual directly, and it has a high freedom of search and efficiency. The main steps of MGA include generating initial population, calculating fitness value, selection, crossover, and mutation.
3.2.1. Set Up the Initial Population
We select ring radii as optimization individual. In order to meet the requirement of minimum element space , an increasing function is used to generate the ring radius as follows:Generate individual repeatedly to make up the population . is the individual number of the population .
3.2.2. Calculate Fitness Value with Function Model
Bring into equation (5), so the maximum element number of each ring can be calculated. Equation (7) can calculate the retention rate of each ring; thus the actual element number in each ring can be acquired. Because the element in each ring is uniform distribution, the angle of th element on the th ring can be calculated as follows:The polar coordinate representation of the individual is as follows:Substitute into equation (4), so the fitness value can be calculated.
3.2.3. Judgement of Termination Condition
If the termination conditions (the number of iterations reaches the limit or the PSLL meets the requirement) are met, output the optimal individual and the minimum PSLL of the current population, end. Otherwise, continue.
3.2.4. Selection Operation
Sort the individual according to fitness value.Select the excellent individual by mean of truncated selection rate .Replicate to make up new population as genetic operator.The number of individuals in is equal to .
3.2.5. Crossover and Mutation
In order to meet the requirement of minimum element spacing , we need to set some constraints in process of crossover and mutation.
Crossover. If the th genes of th and th individuals are selected as crossover gene.
The precrossover genetic operator is The postcrossover genetic operator is where
Mutation. If the th gene of th individuals is selected as mutation gene.
The premutation genetic operator is The postmutation genetic operator is where The individual, after performing the crossover and mutation, forms next generation population .
4. Simulation
In order to verify the efficiency of proposed method, two simulations were performed. We set some common parameters for two simulations, such as selecting equation (4) as objective function and the minimum element spacing ; the radiation pattern in the -region () and -region () is sampled 1024 points. The basic parameters of MGA are set as follows: the number of genetic generations is 200, the individual number of populations is 100, truncate rate is 50%, crossover rate is 20%, mutation rate is 2%, and elitism is also employed.
4.1. Optimizing Sparse Concentric Ring Arrays with the Number of Elements N = 201
We set the same parameter values as [21, 22]; that is, the number of elements = 201, and the aperture = . If ring number = 6, the arrays are almost full and the element spacing is equal to nearly ; there is no room for optimization; thus let the number of rings = 7. Set the coefficient of function model = -0.0072, = 0.0036, and = 0.6973. Five independent experiments were carried out in this simulation, and convergence curve takes the average of 5 independent experiments. The radiation pattern, the element configuration diagram, and the cut surface of = 0, = 0 are from the optimal result of 5 experiments. The element configuration of the best sparse concentric ring array is shown in Figure 4. Figure 5 shows the radiation pattern of the optimal result. The cut surface of = 0 and = 0 is shown in Figure 6. The PSLL of optimal result is -26.21 dB, and it is lower 3.27 dB than [21] and lower 0.76 dB than [22]. Figure 7 shows the convergence characteristics of function model algorithm using MGA, and the proposed method has fast convergence characteristic. Contrast with that of [21, 22] is shown in Table 3.
4.2. Optimizing Sparse Concentric Ring Arrays with the Number of Elements N = 142
According to [21, 22], set the number of elements = 142, the aperture = 9.4, and the number of rings = 6. Set the coefficient of function model = -0.1969, = 0.00228, and = 0.6543. Five independent experiments were carried out in this simulation, and convergence curve takes the average of 5 independent experiments. The radiation pattern, the element configuration diagram, and the cut surface of = 0, = 0 are from the optimal result of 5 experiments. The element configuration of the best sparse concentric ring array is shown in Figure 8. Figure 9 shows the radiation pattern of the optimal result. The cut surface of radiation pattern when = 0 and = 0 is shown in Figure 10. The PSLL of optimal result is -28.19 dB, and it is lower 0.37 dB than [21] and lower 0.12 dB than [22]. Figure 11 shows the convergence characteristics function model algorithm using MGA, and the proposed method has fast convergence characteristic. Contrast with that of [21, 22] is shown in Table 4.
5. Conclusion
In this paper, a function model that presents the relationship between ring radius and the degree of sparsity is proposed, which is embedded into the process of MGA to synthesize the sparse concentric ring arrays for low PSLL. Comparing the optimal result with other literatures proves that the proposed method is an efficient way to reduce the PSLL and has a faster convergence rate. It is reasonable and efficient to select the element number in the rings according to the function model.
Data Availability
The data of this paper can be accessed for the public; the data have been uploaded to github repository and all the data present in the research are available. The chart, docx data used to support the findings of this study, the result of figure, and the program of this research have been deposited to the github repository; these data can be accessed through https://github.com/1367356/Optimization-of-sparse-concentric-ring-arrays-for-low-sidelobe.
Conflicts of Interest
Kesong Chen, Yafei Li, and Jiajia Shi declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The research is funded by joint fund for civil aviation (U1233103).