A Phase Angle-Modulated Bat Algorithm with Application to Antenna Topology Optimization

: This paper proposes a phase angle-modulated bat algorithm (P-AMBA) for high-dimensional binary optimization. The idea was to reduce the optimization time by introducing angle modulation technology to reduce the optimization dimensions. Different from the original angle-modulated bat algorithm (AMBA), the control of the trigonometric generating function cosine wave is by introducing new parameters, thereby improving the perturbation ability of the function curve near the x -axis. P-AMBA can explore more 0/1 solutions, and it has advantages in optimizing convergence speed and global search capabilities. The numerical results of the 0–1 knapsack problem tests show that P-AMBA is superior to the contrast algorithms on optimization ability and optimization time. Finally, the experimental result of a compact dual-band planar monopole antenna design showed the effectiveness of P-AMBA in engineering applications.


Introduction
With the rapid development of modern society, big data have become a research hotspot and a large number of data sets need to be optimized. High-dimensional optimization problems exist in many application scenarios, such as complex neural network weight optimization [1], text mining and analysis [2], sensor networking [3], vehicle routing [4], and antenna topology optimization [5]. To solve such optimization problems, meta-heuristic algorithms have been proposed as a reliable solution [6]. As an improvement of heuristic algorithms, meta-heuristic algorithms are the combination of random algorithms and local search algorithms. Meta-heuristic algorithms can find the optimal solution to the problem under limited time and space conditions. A large number of meta-heuristic algorithms have been proposed including the genetic algorithm (GA) [7], particle swarm optimization (PSO) [8], differential evolution (DE) [9], artificial bee colony (ABC) [10], the dragonfly algorithm (DA) [11], the fireworks algorithm (FA) [12], brain storm optimization (BSO) [13], harmony search (HS) [14], moth search (MS) [15], and the bat algorithm (BA) [16].
The BA, proposed by Yang in 2010, is a meta-heuristic algorithm developed according to the principle of bat echolocation [16]. Compared with other meta-heuristic algorithms, BA has two characteristics: (i) frequency tuning can facilitate the control of the optimization direction of the algorithm, and (ii) dynamic conversion of global search and local search can effectively avoid premature convergence of the algorithm. The advantages of BA are its simplicity, fewer parameters, robustness, and ease of implementation. Due to its superiority, BA has been applied in various fields, such as wireless sensor network deployment in the 3D environment [17], data clustering [18], the low-carbon job shop scheduling problem [19], optimization of reactive power dispatch [20], and the blood vessel segmentation problem [21]. BA was first proposed to address the optimization problem of the continuous

The Bat Algorithm
The BA is a heuristic search algorithm that was proposed by Yang in 2010 based on swarm intelligence [16]. It is an effective method to search for the global optimal solution. It simulates bats' behavior of using ultrasound to detect and locate obstacles or prey, as shown in Figure 1. The bionic principle of BA is to map the population of individual bats to the n-dimensional problem space and then to simulate the optimization process. The fitness function value of the solution is used to measure the advantages and disadvantages of the bat's position. The survival process of the most suitable individual is similar to the iterative process of replacing the poorer feasible solutions with better feasible solutions in the optimization and search process. Compared with other algorithms, BA is superior in terms of accuracy and effectiveness, and it does not have many parameters that need to be adjusted.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 3 of 20 effectiveness of the proposed P-AMBA. Finally, Section 5 concludes the paper and provides suggestions for future studies.

The Bat Algorithm
The BA is a heuristic search algorithm that was proposed by Yang in 2010 based on swarm intelligence [16]. It is an effective method to search for the global optimal solution. It simulates bats' behavior of using ultrasound to detect and locate obstacles or prey, as shown in Figure 1. The bionic principle of BA is to map the population of individual bats to the n-dimensional problem space and then to simulate the optimization process. The fitness function value of the solution is used to measure the advantages and disadvantages of the bat's position. The survival process of the most suitable individual is similar to the iterative process of replacing the poorer feasible solutions with better feasible solutions in the optimization and search process. Compared with other algorithms, BA is superior in terms of accuracy and effectiveness, and it does not have many parameters that need to be adjusted.
where X ( ) k i t and V ( ) k i t indicate the k-th dimension position and velocity of the i-th bat at iteration t, respectively. X gbest is the best solution thus far, and F i indicates the frequency of the i-th bat that is updated in each course of iterations as follows: where β is the random number uniformly distributed in [0,1]. According to Equations (1) and (2), different frequencies can improve the diversification of artificial bats and reach the optimal solution. To improve the local search capability, each bat has a local search strategy. When In the BA, the optimization parameter is n-dimensional, and the number of bat populations is n. Each artificial bat has a frequency vector X k i (i = 1, 2, 3, . . . , n, k = 1, 2, 3, . . . , N), a velocity vector V k i (i = 1, 2, 3, . . . , n, k = 1, 2, 3, . . . , N), and a position vector F i (i = 1, 2, 3, . . . , n), which are updated at iteration t as Equations (1)-(3): where X k i (t) and V k i (t) indicate the k-th dimension position and velocity of the i-th bat at iteration t, respectively. X gbest is the best solution thus far, and F i indicates the frequency of the i-th bat that is updated in each course of iterations as follows: where β is the random number uniformly distributed in [0,1]. According to Equations (1) and (2), different frequencies can improve the diversification of artificial bats and reach the optimal solution. To improve the local search capability, each bat has a local search strategy. When rand > r, the new position X k i (t) new changes with a random walk: where ε ∈ [−1, 1], rand ∈ [0, 1], r is pulse emission rate, and A i (t) is the loudness value of all particles at iteration t. As a result of the introduction of pulse emission rate (r) and loudness (A), the bat algorithm can be seen as a balance between particle swarm search and intensive local search. With the update of the two parameters, the optimization enters the exploration stage from the exploitation state. For each iteration of the algorithm, the pulse emission rate (r) and loudness (A) are updated as follows: where α and γ are constants, with α being similar to the cooling factor in the simulated annealing (SA) [36]. The loudness and pulse emission rate are updated in each iteration until the artificial bat finds the optimal solution or reaches the maximum number of iterations. For any α > 0;γ < 1: The pseudo-code of the bat algorithm is shown in Algorithm 1: , r is pulse emission rate, and ( ) i A t is the loudness value of all particles at iteration t. As a result of the introduction of pulse emission rate (r) and loudness (A), the bat algorithm can be seen as a balance between particle swarm search and intensive local search. With the update of the two parameters, the optimization enters the exploration stage from the exploitation state. For each iteration of the algorithm, the pulse emission rate (r) and loudness (A) are updated as follows: where α and γ are constants, with α being similar to the cooling factor in the simulated annealing (SA) [36]. The loudness and pulse emission rate are updated in each iteration until the artificial bat finds the optimal solution or reaches the maximum number of iterations. For any 0 α > ; γ < 1 : The pseudo-code of the bat algorithm is shown in Algorithm 1:

Binary Bat Algorithm
The emergence of the bat algorithm was as a means to solve the continuous numerical optimization problem and to find the necessary 0/1 string in the feature selection problem. Additionally, the BBA introduces the transfer function to change the position update equation [22]. The velocity and frequency of the BBA are defined in continuous space. Compared with BA, BBA was additionally modified as follows: • A transformation function was added to map a continuous position vector to a discrete bit vector.

•
The random walk Equation (4) that improves the local search ability was deleted, as it is not suitable for the binary optimization algorithm.

Binary Bat Algorithm
The emergence of the bat algorithm was as a means to solve the continuous numerical optimization problem and to find the necessary 0/1 string in the feature selection problem. Additionally, the BBA introduces the transfer function to change the position update equation [22]. The velocity and frequency of the BBA are defined in continuous space. Compared with BA, BBA was additionally modified as follows:

•
A transformation function was added to map a continuous position vector to a discrete bit vector.

•
The random walk Equation (4) that improves the local search ability was deleted, as it is not suitable for the binary optimization algorithm.
In the BBA, using Equation (8) as the transfer function, the probability of the individual velocity of the BBA population is mapped to its position vector element flip. The greater the absolute value of the independent variable speed, the greater the dependent variable result, and the greater the probability of element flipping. The location update rules are as follows: where X k i (t) and V k i (t) indicate the k-th dimension position and velocity of the i-th bat at iteration t, and (X k i (t)) −1 indicates the complement of X k i (t).

AMBA
Even though traditional BBA can be optimized in the binary search space, the optimization time is longer, and the optimization ability is limited [24]. Thus, the idea of angle modulation was introduced. The AMBA is an improvement of the BBA [22], and its core idea is to introduce a trigonometric function as a bit vector generator to map the real-valued solutions into binary ones. This method was first used as the angle-modulated (AM) technique for communication signal processing [37][38][39].
This angle modulation idea uses a compound sin/cos generating function: where x is a sequence of equal intervals, the number of which needs to be determined by the number of optimized dimensions, the coefficient a represents the horizontal shift of the function; b represents the maximum frequency of the sine wave; c represents the frequency of the cosine wave (which changes the rate at which the frequency of the sine wave changes); and d represents the vertical shift of the function. The default value of the four parameters in the generating function is a = 0, b = 1, c = 0.5, and d = 0, and a function diagram is shown in Figure 2.

AMBA
Even though traditional BBA can be optimized in the binary search space, optimization time is longer, and the optimization ability is limited [24]. Thus, the ide angle modulation was introduced. The AMBA is an improvement of the BBA [22], and core idea is to introduce a trigonometric function as a bit vector generator to map the r valued solutions into binary ones. This method was first used as the angle-modula (AM) technique for communication signal processing [37][38][39].
This angle modulation idea uses a compound sin/cos generating function: where x is a sequence of equal intervals, the number of which needs to be determined the number of optimized dimensions, the coefficient a represents the horizontal shif the function; b represents the maximum frequency of the sine wave; c represents frequency of the cosine wave (which changes the rate at which the frequency of the s wave changes); and d represents the vertical shift of the function. The default value of four parameters in the generating function is a = 0, b = 1, c = 0.5, and d = 0, and a funct diagram is shown in Figure 2. After introducing the angle modulation strategy, a high-dimensional binary sea space can be obtained by optimizing the four-dimensional continuous search space. updated range of the four variables a, b, c, and d is as follows: After introducing the angle modulation strategy, a high-dimensional binary search space can be obtained by optimizing the four-dimensional continuous search space. The updated range of the four variables a, b, c, and d is as follows: The change from a high-dimensional search space to four-dimensional search space greatly reduces the optimization time under the premise of ensuring optimization accuracy.

P-AMBA
In AMBA, after each iteration of the four-dimensional parameters, the richness of the discrete search space mapped by the trigonometric generating function is not enough to meet the optimization requirements. To maximize the impact of the parameter changes on the generating function, the P-AMBA proposed in this paper adds a phase adjustment to the cosine wave and an amplitude adjustment to the sine wave. P-AMBA can increase the global search ability of the optimization process to a certain extent. The generating function is described as follows: When x takes the number of equal intervals from 1 to n dimensions, the function g(x) will generate the corresponding value. By converting g(x)(x = 0, 1, 2, · · ·, n) into a 0/1 bit string through the conversion of Equation (12), we can achieve a mapping relationship. The rule of the conversion function is to take 0 when g(x) ≤ 0 and to take 1 when g(x) > 0.
The original generating function lacks the control of the amplitude of the sine wave, so that the update of the vertical displacement d has less influence on the function curve. To make up for the lack of amplitude, the parameter e that controls the amplitude of the sine wave is added. Similarly, the cosine wave lacks phase adjustment. Adding a parameter g after the cosine wave can increase the fluctuation ability of the function near the zero points. The improved function adds two parameters, so the individual bat in the bat algorithm is six-dimensional: A comparison of the figures for the generating function between P-AMBA and AMBA is represented in Figure 3. It can be seen from the figure that when a parameter e is added, the parameter d that controls the phase of the entire function has a greater impact on the up and down positions of the function. This can better control the position of the function and obtain the desired 0/1 solution set. At the same time, by adding the parameter g to control the phase of the cosine wave, the function is more complicated. Each iterative update of the value of the six-dimensional parameter will result in a greater change in the overall mapping relationship of the function, and it will further enrich the exploration of possible solutions. P-AMBA has more parameters than BBA and AMBA, which may increase the complexity of the algorithm in the optimization process. However, the operation of adding new parameters can significantly improve the optimization ability in the high-dimensional binary search space. After introducing the idea of angle modulation, the optimization process is essentially a continuous optimization process applied under binary demand. The key step is to substitute six parameters into the trigonometric generating function to generate an n-dimensional bit string and then to calculate the fitness function. The pseudo-code of the P-AMBA is shown in Algorithm 2. To better demonstrate the mapping process, the After introducing the idea of angle modulation, the optimization process is essentially a continuous optimization process applied under binary demand. The key step is to substitute six parameters into the trigonometric generating function to generate an ndimensional bit string and then to calculate the fitness function. The pseudo-code of the P-AMBA is shown in Algorithm 2. To better demonstrate the mapping process, the principle through which P-AMBA obtains discrete n-dimensional search space results by optimizing the continuous six-dimensional search space is shown in Figure 4.

Experimental Results and Discussions
The processor used for the simulation experiment environment was an Intel(R) Core

Experimental Results and Discussions
The processor used for the simulation experiment environment was an Intel(R) Core (TM) i5-8265U 1.80 GHz (with 8.0 GB RAM and a 64-bit Windows10 operating system), and the simulation software was Matlab 2016a and HFSS 15.

Zero-One Knapsack Problems
The 0-1 knapsack problem is a common and typical combinatorial optimization nondeterministic Polynomial complete (NP-C) problem, which is often used for numerical testing of binary optimization algorithms. The knapsack problem was first proposed by Dantzig in 1957 [40,41]. Classical 0-1 knapsack problems can be described as follows: where p i and w i represent the profit and weight of the individual items i, N is the number of individuals, and x represents the decision variable of the i-th item. If the item i is packed, x takes 1; otherwise, x takes 0. The backpack problem can be described as follows: On the premise of ensuring that the total weight f (w) of the objects in the backpack does not exceed the maximum capacity C that the backpack can bear, we find the distribution strategy with the highest total value f (p) of the items. The profit and weight of the items are positive.
In this paper, 12 different dimensional test examples of the knapsack problem were selected to compare P-AMBA with other optimization algorithms. The data of the classic knapsack problem test cases k1-k5 are shown in Table 1, where 'D' represents the dimension of the backpack, that is, how many items are there. 'Parameter (w, p, C)' represents the weight and value of the item and the maximum load capacity of the backpack, respectively, while 'Opt' is the optimal value for this instance. Meanwhile, cases k6-k12 are high-scale knapsack problems. The item data are generated by a random generator to compare the application of several optimization algorithms in a high-dimensional environment. The weight w i (i = 1, 2, · · · , N) ∈ [5,20], and the profit p i (i = 1, 2, · · · , N) ∈ [50, 100] [42]. The relevant data of these randomly generated backpack items are shown in Table 2. The 'Total values' denotes the total values of all items. In the k6-k12 cases, the formula 3 4 ∑ N i=1 w i was used to calculate the total capacity of items in the backpack.
For the example test of the knapsack problem, the classic BPSO [23], the traditional BBA [22], the AMBA [27], and the A-AMBA [28] (the latter two also used the idea of angle modulation) were selected for comparison with the proposed P-AMBA. The relevant parameter settings of the optimization algorithms in the verification tests are shown in Table 3. The parameter range is set according to [28], which is suitable for the knapsack problem. In order to objectively compare the optimization capabilities of different optimization algorithms, the parameter settings of this paper are all initial settings. The population size was set to 30, and the number of iterations was set to 500. To avoid accidental deviations, the algorithm was run 30 times independently for each function. Table 4 shows the test results of the five optimization algorithms after the optimization, including the maximum, minimum, average, standard deviation, and optimization time. In the table, 'No.' and 'Alg.' represent the test examples and the optimization algorithm types of the knapsack problem, respectively. 'Best' represents the maximum value of the 30 population optimal values, which is the optimal value that the optimization algorithm can find in the optimization process. 'Worst' represents the worst value of the optimal value of the 30 populations. 'Mean' is the average optimal value of the 30 optimization results. 'Std' is the standard deviation, reflecting the degree of dispersion of the 30 optimization results. 'Time' is the average time consumed for each optimization, reflecting the optimization speed of the optimization algorithm. Moreover, the average value of the iterative data of the examples k7-k12 in 30 optimization tests was selected, and the iterative convergence of the five algorithms was compared, as shown in Figure 5.    The optimal, worst, and mean values in Table 4 show that the higher the dimensionality of the optimization problem, the more prominent the optimization ability of P-AMBA. By comparing the optimization time of the five algorithms, although the BPSO algorithm took less time, it did not perform as well as the other algorithms on high-dimensional problems. With the introduction of angle modulation technology, the time consumption of the bat algorithms was greatly reduced. In a high-dimensional environment, the time consumption of the P-AMBA algorithm was only half that of the traditional binary bat algorithm.
The comparison of the convergence curves in Figure 5 obviously shows that P-AMBA has a strong ability to optimize high-dimensional problems. After 500 iterations, a better solution set of optimal values was found by P-AMBA in comparison to the other optimization algorithms. Different from AMBA and A-AMBA, in P-AMBA, in addition to the effect of the amplitude amplifying the vertical displacement coefficient, new parameters were also introduced to adjust the phase of the cosine wave, thereby further affecting the disturbance of the function and adding more 0/1 combination possibilities.

Compact Dual-Band Planar Monopole Antenna Design
Antennas are widely used in mobile communications, radar, and satellite communications as they play an important role in the transmitting and receiving of signals. Since antenna performance determine the quality of a wireless communication system and different antenna configurations have different antenna performances, the optimization of the antenna configurations has attracted a lot of attention. Traditional antenna optimization starts from a given initial geometric structure and only adjusts the size of the relevant geometric parameters. However, this kind of optimization method is not conducive to exploring the possibility of novel antenna structures. The topology optimization used in this paper, also called pixel optimization, is a technique that can automatically generate an antenna topology suitable for the design requirements in a given design area. It divides the design space into several pixels represented by a matrix with "1" (conductor) and "0" (air). Through some stochastic optimization techniques, we can find the optimal pixel distribution mode (i.e., the optimal antenna topology) that satisfies the design requirements.  The optimal, worst, and mean values in Table 4 show that the higher the dimensionality of the optimization problem, the more prominent the optimization ability of P-AMBA. By comparing the optimization time of the five algorithms, although the BPSO algorithm took less time, it did not perform as well as the other algorithms on high-

Antenna Topology Optimization Problem Formulation
In antenna topology optimization, the design area is divided into pixel blocks of the same size, and air and conductor are filled in each pixel block according to the binary mapping relationship. Therefore, the problem of antenna topology designs can be expressed as a binary optimization problem: where x = [x 1 ,x 2 , . . . ,x n ] is an n-dimensional binary design variable, representing a binary code mapped from the specific antenna geometric structure; x * is the obtained optimal antenna design result; f (x) is the objective function of the antenna design, such as reflection coefficients, gain, or efficiency. Figure 6 shows the mapping process of the antenna structure. If the pixel is filled with the conductor, it is mapped as 1; otherwise, it is mapped as 0.
automatically generate an antenna topology suitable for the design requirements in a given design area. It divides the design space into several pixels represented by a matrix with "1" (conductor) and "0" (air). Through some stochastic optimization techniques, we can find the optimal pixel distribution mode (i.e., the optimal antenna topology) that satisfies the design requirements.

Antenna Topology Optimization Problem Formulation
In antenna topology optimization, the design area is divided into pixel blocks of the same size, and air and conductor are filled in each pixel block according to the binary mapping relationship. Therefore, the problem of antenna topology designs can be expressed as a binary optimization problem: (16) where x = [x1,x2,…,xn] is an n-dimensional binary design variable, representing a binary code mapped from the specific antenna geometric structure; x * is the obtained optimal antenna design result; f(x) is the objective function of the antenna design, such as reflection coefficients, gain, or efficiency. Figure 6 shows the mapping process of the antenna structure. If the pixel is filled with the conductor, it is mapped as 1; otherwise, it is mapped as 0.

Antenna Topology Optimization Design
Multi-band antennas meet the needs of different communication standards with only one device. The advantage of a multi-band device is that it has adaptability and scalability, and the number of frequency bands can differ. Therefore, multi-band antennas are very popular in the field of antenna optimization designs.
To verify the effectiveness of the P-AMBA, herein, we present a design example of a dual-band planar monopole antenna and compare it with the BPSO [23], BBA [22], AMBA [27], and A-AMBA [28] optimization algorithms. The initial antenna geometry is given in Figure 7, where the left and right pictures are the top views and side views, respectively. It consists of a 14.6 × 12.8 mm rectangular patch and a 16.8 × 4 mm rectangular ground plane. The antenna is fed by a 50Ω microstrip line and printed on a 16.8 × 24.6 mm FR4 substrate with a 1.6 mm thickness, a 4.4 permittivity, and a 0.02 loss tangent. The antenna was designed to achieve the reflection coefficients lower than −10 dB (i.e., S11 < −10 dB) in the two frequency bands of 2.

Antenna Topology Optimization Design
Multi-band antennas meet the needs of different communication standards with only one device. The advantage of a multi-band device is that it has adaptability and scalability, and the number of frequency bands can differ. Therefore, multi-band antennas are very popular in the field of antenna optimization designs.
To verify the effectiveness of the P-AMBA, herein, we present a design example of a dual-band planar monopole antenna and compare it with the BPSO [23], BBA [22], AMBA [27], and A-AMBA [28] optimization algorithms. The initial antenna geometry is given in Figure 7, where the left and right pictures are the top views and side views, respectively. It consists of a 14.6 × 12.8 mm rectangular patch and a 16.8 × 4 mm rectangular ground plane. The antenna is fed by a 50Ω microstrip line and printed on a 16.8 × 24.6 mm FR4 substrate with a 1.6 mm thickness, a 4.4 permittivity, and a 0.02 loss tangent. The antenna was designed to achieve the reflection coefficients lower than −10 dB (i.e., S 11 < −10 dB) in the two frequency bands of 2.3−2.5 GHz and 5.1−5.9 GHz, covering the entire 2.4/5.2/5.8 GHz WLAN bands. The fitness function can be expressed as: where f i denotes the i th sampling frequencies within the given operating bands; S 11 (f i ) is the reflection coefficient of sample f i ; m is the total number of samples; and F is the fitness value.
where fi denotes the i th sampling frequencies within the given operating the reflection coefficient of sample fi; m is the total number of samples; and value. The initial antenna setting discretized the rectangular patch into 8 × 56 pixels. The size of each pixel was 2 × 2 mm. To ensure that the two s were only in contact at the opposite corners could pass the current, overlapped at the edge with an overlapping width of 0.2 mm [43]. In parameter settings, the population size was set as 10, and the maximum it was set as 30. To ensure that the antenna has enough metal sheets to ac performance, the curve of the angle modulation function was controlled a as much as possible, so the initial range of the parameters d and g was set The setting range of the other parameters of the algorithm is the same as the above knapsack problem. After testing, these parameter ranges can be antenna designs.
A schematic diagram of the application of phase angle-modulated antenna optimization is shown in Figure 8. By optimizing the six-dimensi binary data of the 8 × 7 matrix corresponding to the antenna model cou Then, the HFSS was called for electromagnetic simulation, and the obta values were calculated as the fitness values. Finally, an iterative optimization algorithm was performed. The initial antenna setting discretized the rectangular patch into 8 × 7 pixels, that is, 56 pixels. The size of each pixel was 2 × 2 mm. To ensure that the two sub patches that were only in contact at the opposite corners could pass the current, each pixel was overlapped at the edge with an overlapping width of 0.2 mm [43]. In the algorithm parameter settings, the population size was set as 10, and the maximum iteration number was set as 30. To ensure that the antenna has enough metal sheets to achieve radiation performance, the curve of the angle modulation function was controlled above the y-axis as much as possible, so the initial range of the parameters d and g was set within [−0.5, 1]. The setting range of the other parameters of the algorithm is the same as the test data of the above knapsack problem. After testing, these parameter ranges can be well applied to antenna designs.
A schematic diagram of the application of phase angle-modulated technology in antenna optimization is shown in Figure 8. By optimizing the six-dimensional vector, the binary data of the 8 × 7 matrix corresponding to the antenna model could be mapped. Then, the HFSS was called for electromagnetic simulation, and the obtained return loss values were calculated as the fitness values. Finally, an iterative update of the optimization algorithm was performed.
A comparison of the convergence performance of the five optimization algorithms is shown in Figure 9, where the abscissa and the ordinate represent the number of iterations and the fitness value, respectively. It can be seen that the P-AMBA proposed in this paper achieved a good balance between exploitation and exploration, showing a strong optimization ability in the early stage of iteration and not falling into the local optimum prematurely. The final antenna topology and S 11 results optimized by the P-AMBA are shown in Figure 10. The obtained antenna structure met the design objective requirements of the dual-band operation.
X a X b X c X d X e X g

Electromagnetic simulation with HFSS
Get the return loss value (S11)

Equation (17)
Fitness value calculation A comparison of the convergence performance of the five optimization algorithms is shown in Figure 9, where the abscissa and the ordinate represent the number of iterations and the fitness value, respectively. It can be seen that the P-AMBA proposed in this paper achieved a good balance between exploitation and exploration, showing a strong optimization ability in the early stage of iteration and not falling into the local optimum prematurely. The final antenna topology and S11 results optimized by the P-AMBA are shown in Figure 10. The obtained antenna structure met the design objective requirements of the dual-band operation.  A comparison of the convergence performance of the five optimization algorithms is shown in Figure 9, where the abscissa and the ordinate represent the number of iterations and the fitness value, respectively. It can be seen that the P-AMBA proposed in this paper achieved a good balance between exploitation and exploration, showing a strong optimization ability in the early stage of iteration and not falling into the local optimum prematurely. The final antenna topology and S11 results optimized by the P-AMBA are shown in Figure 10. The obtained antenna structure met the design objective requirements of the dual-band operation.  The antenna topology design is different from the traditional antenna structure design. It starts from a simple initial geometry and exhibits superior electrical performance with a compact size. This section used a design example of a dual-band planar monopole antenna to confirm the advantages of P-AMBA in global optimization capabilities and the convergence rate. The results showed that the antenna structure designed by this method is novel and can be used for antenna designs in limited surrounding environment or unexplored areas. The antenna topology design is different from the traditional antenna structure design. It starts from a simple initial geometry and exhibits superior electrical performance with a compact size. This section used a design example of a dual-band planar monopole antenna to confirm the advantages of P-AMBA in global optimization capabilities and the convergence rate. The results showed that the antenna structure designed by this method is novel and can be used for antenna designs in limited surrounding environment or unexplored areas.

Conclusions
This paper proposed an improved binary bat algorithm for antenna topology designs. The idea was to reduce the optimization time by introducing angle modulation technology to reduce the optimization dimensions. Compared with other optimization algorithms, it also has advantages in global search capability and convergence speed. The knapsack problem was used to verify the optimization ability of P-AMBA in highdimensional problems. The results showed that the higher the dimension, the more obvious the advantage of this algorithm. Finally, through a topology design example of a compact dual-band planar monopole antenna, the superiority of the proposed P-AMBA in the optimization of high-dimensional discrete search space was verified.
In the antenna optimization example, since much time was spent on the antenna electromagnetic simulation process, the advantages of the algorithm in terms of optimization time were not obvious. Subsequent research will aim at neural networksbased antenna designs by introducing surrogate models instead of electromagnetic simulation tools to achieve fast antenna topology optimization. Moreover, a multiobjective antenna design will also be the task of follow-up research.

Conclusions
This paper proposed an improved binary bat algorithm for antenna topology designs. The idea was to reduce the optimization time by introducing angle modulation technology to reduce the optimization dimensions. Compared with other optimization algorithms, it also has advantages in global search capability and convergence speed. The knapsack problem was used to verify the optimization ability of P-AMBA in high-dimensional problems. The results showed that the higher the dimension, the more obvious the advantage of this algorithm. Finally, through a topology design example of a compact dual-band planar monopole antenna, the superiority of the proposed P-AMBA in the optimization of high-dimensional discrete search space was verified.
In the antenna optimization example, since much time was spent on the antenna electromagnetic simulation process, the advantages of the algorithm in terms of optimization time were not obvious. Subsequent research will aim at neural networks-based antenna designs by introducing surrogate models instead of electromagnetic simulation tools to achieve fast antenna topology optimization. Moreover, a multi-objective antenna design will also be the task of follow-up research.